Abstract:We construct finite element Stokes complexes on tetrahedral meshes in three-dimensional space. In the lowest order case, the finite elements in the complex have 4, 18, 16, and 1 degrees of freedom, respectively. As a consequence, we obtain gradcurl-conforming finite elements and inf-sup stable Stokes pairs on tetrahedra which fit into complexes. We show that the new elements lead to convergent algorithms for solving a gradcurl model problem as well as solving the Stokes system with precise divergence-free condition. We demonstrate the validity of the algorithms by numerical experiments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to construct a finite - element Stokes complex on three - dimensional tetrahedral meshes, in order to maintain the exact divergence - free condition (i.e., incompressibility) in numerical discretization, and to provide an effective method for solving high - order curl problems and Stokes systems.
Specifically, the authors construct a simple discrete sub - complex:
\[ 0 \rightarrow \mathbb{R} \rightarrow \Sigma_h \rightarrow V_h \rightarrow \Sigma^+_h \rightarrow W_h \rightarrow 0, \]
where each space has 4, 18, 16 and 1 degrees of freedom respectively. By introducing the modified Bernardi - Raugel bubble function and Poincaré operator, they successfully reduce the polynomial order and avoid the extensive use of macro - element structures.
In addition, these new finite elements can be applied to solve high - order curl problems and Stokes equations that satisfy the exact divergence - free condition. In particular, for high - order curl problems, these elements ensure the convergence of the algorithm and simplify the previous complex constructions.
In summary, the main contribution of this paper is to provide a new finite - element construction method, which enables more effective numerical simulation and maintains the accuracy of physical properties when dealing with certain problems in fluid mechanics and electromagnetics.