Abstract:Composite basis functions for pyramidal elements on the spaces $H^1(\Omega)$, $H(\mathrm{curl},\Omega)$, $H(\mathrm{div},\Omega)$ and $L^2(\Omega)$ are presented. In particular, we construct the lowest-order composite pyramidal elements and show that they respect the de Rham diagram, i.e. we have an exact sequence and satisfy the commuting property. Moreover, the finite elements are fully compatible with the standard finite elements for the lowest-order Raviart-Thomas-Nédélec sequence on tetrahedral and hexahedral elements. That is to say, the new elements have the same degrees of freedom on the shared interface with the neighbouring hexahedral or tetrahedra elements, and the basis functions are conforming in the sense that they maintain the required level of continuity (full, tangential component, normal component, ...) across the interface. Furthermore, we study the approximation properties of the spaces as an initial partition consisting of tetrahedra, hexahedra and pyramid elements is successively subdivided and show that the spaces result in the same (optimal) order of approximation in terms of the mesh size $h$ as one would obtain using purely hexahedral or purely tetrahedral partitions.
What problem does this paper attempt to address?
This paper attempts to solve the meshing and basis function selection problems encountered in finite element analysis in complex three - dimensional applications, especially the construction of the lowest - order composite finite element exact sequence for pyramidal elements. Specifically:
1. **Meshing problem**: When dealing with complex three - dimensional domains, hexahedral elements are efficient but difficult to generate; while tetrahedral elements are easy to generate but not as efficient as hexahedral elements. To combine the advantages of both, a tetrahedron - hexahedron - pyramid (THP) hybrid mesh partition is used. Among them, the pyramid elements are used to connect hexahedral and tetrahedral elements.
2. **Basis function selection problem**: For pyramid elements, choosing appropriate basis functions to ensure compatibility with adjacent hexahedral and tetrahedral elements is a challenge. The article points out that in order to maintain the required continuity at the interface (such as full continuity, tangential component continuity, normal component continuity, etc.), basis functions that meet these requirements must be constructed.
3. **Composite finite element method**: The article proposes a method based on composite or macro - element basis functions to construct the lowest - order finite element spaces \(H^1(\Omega)\), \(H(\text{curl},\Omega)\), \(H(\text{div},\Omega)\) and \(L^2(\Omega)\) on pyramid elements. This method avoids using rational basis functions by subdividing the pyramid elements into two tetrahedra and ensures that these new elements are fully compatible with the standard tetrahedral and hexahedral lowest - order Raviart - Thomas - Nédélec sequences.
4. **Guarantee of mathematical properties**: The constructed lowest - order composite pyramid elements not only respect the de Rham diagram, that is, they have an exact sequence and satisfy the commutativity property, but also can maintain the optimal approximation order after the initial THP partition is gradually subdivided.
In summary, this paper aims to solve the difficult problem of basis function selection for pyramid elements in THP hybrid meshes by introducing the lowest - order composite finite element exact sequence, thereby improving the effectiveness and accuracy of finite element analysis for complex three - dimensional problems.