A new family of stable mixed finite elements for the 3D Stokes equations 论文

2004Mathematics of Computation引用 274
Advanced Numerical Methods in Computational MathematicsAdvanced Mathematical Modeling in EngineeringElectromagnetic Simulation and Numerical Methods

摘要

A natural mixed-element approach for the Stokes equations in the velocity-pressure formulation would approximate the velocity by continuous piecewise-polynomials and would approximate the pressure by discontinuous piecewise-polynomials of one degree lower. However, many such elements are unstable in 2D and 3D. This paper is devoted to proving that the mixed finite elements of this <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="bold upper P Subscript k"> <mml:semantics> <mml:msub> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:mi mathvariant="bold">P</mml:mi> </mml:mrow> <mml:mi>k</mml:mi> </mml:msub> <mml:annotation encoding="application/x-tex">\mathbf {P}_k</mml:annotation> </mml:semantics> </mml:math> </inline-formula> - <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper P Subscript k minus 1"> <mml:semantics> <mml:msub> <mml:mi>P</mml:mi> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:mi>k</mml:mi> <mml:mo> − </mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msub> <mml:annotation encoding="application/x-tex">P_{k-1}</mml:annotation> </mml:semantics> </mml:math> </inline-formula> type when <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="k greater-than-or-equal-to 3"> <mml:semantics> <mml:mrow> <mml:mi>k</mml:mi> <mml:mo> ≥ </mml:mo> <mml:mn>3</mml:mn> </mml:mrow> <mml:annotation encoding="application/x-tex">k \ge 3</mml:annotation> </mml:semantics> </mml:math> </inline-formula> satisfy the stability condition—the Babuška-Brezzi inequality on macro-tetrahedra meshes where each big tetrahedron is subdivided into four subtetrahedra. This type of mesh simplifies the implementation since it has no restrictions on the initial mesh. The new element also suits the multigrid method.