Tangential-Normal Decompositions of Finite Element Differential Forms

Long Chen,Xuehai Huang
2024-10-27
Abstract:The paper introduces a novel tangential-normal ($t$-$n$) decomposition for finite element differential forms. Its main contribution is the development of a $t$-$n$ basis where the degrees of freedom and shape functions are explicitly dual to each other. This duality simplifies the assembly of stiffness matrices and enhances the efficiency of interpolation and numerical integration in finite element methods. Additionally, the well-documented Lagrange element basis can be used to expedite implementation. This paper focuses on the full polynomial spaces, excluding trimmed polynomial differential forms, and provides a new perspective on constructing finite element differential forms.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the finite element method, how to deal with the discretization problem of differential forms more efficiently. Specifically, the author introduced a novel tangential - normal (t - n) decomposition method to handle finite - element differential forms. The main contribution of this method lies in the development of a set of t - n bases, under which the degrees of freedom and the shape functions are explicitly dual, thus simplifying the assembly of the stiffness matrix and improving the efficiency of interpolation and numerical integration. ### Detailed Interpretation: 1. **Research Background**: - **Finite Element Exterior Calculus (FEEC)** is a framework that combines the finite - element method with differential geometry to solve partial differential equations (PDEs), especially equations involving differential forms. FEEC preserves key geometric and topological properties, ensuring the stability and accuracy of numerical solutions. 2. **Main Problem**: - In the finite - element method, dealing with the discretization of differential forms is a complex and computationally intensive process. Traditional basis selections often lead to an unclear relationship between the degrees of freedom and the shape functions, increasing the complexity of stiffness matrix assembly and affecting the efficiency of interpolation and numerical integration. 3. **Solution**: - **Tangential - Normal (t - n) Decomposition**: The author proposed a new t - n decomposition method, which decomposes the differential form into tangential and normal parts. Through this decomposition, a set of bases can be constructed such that there is an explicit dual relationship between the degrees of freedom and the shape functions. - **Dual Bases**: Two sets of bases are defined in the paper, one for defining the degrees of freedom (DoFs) and the other for the shape functions. These two sets of bases are dual, which means they correspond to each other mathematically, simplifying the calculation process. - **Application of Lagrange Element Bases**: The author utilized the existing Lagrange element bases to accelerate the implementation of this new method, further improving the computational efficiency. 4. **Theoretical Contribution**: - A new perspective for constructing finite - element differential forms was proposed, especially in the construction in the full polynomial space without involving trimmed polynomial differential forms. - Through t - n decomposition, the handling of differential forms in high - dimensional spaces was simplified, especially in the application on the symplectic complex. 5. **Practical Application**: - This new t - n decomposition method can significantly improve the efficiency and accuracy of the finite - element method in dealing with complex geometric structures and high - order differential equations, and is especially suitable for physical simulations and engineering calculations that require precise preservation of geometric and topological properties. ### Summary: The core problem of this paper is to improve the discretization process of finite - element differential forms by introducing the tangential - normal (t - n) decomposition method, thereby improving the computational efficiency and the accuracy of numerical solutions. This method not only simplifies the assembly of the stiffness matrix but also improves the efficiency of interpolation and numerical integration, providing a more effective numerical solution tool for complex physical and engineering problems.