A primal finite element scheme of the Hodge Laplace problem

Shuo Zhang
DOI: https://doi.org/10.48550/arXiv.2208.00575
2022-08-01
Abstract:In this paper, a unified family, for any $n\geqslant 2$ and $1\leqslant k\leqslant n-1$, of nonconforming finite element schemes are presented for the primal weak formulation of the $n$-dimensional Hodge-Laplace equation on $H\Lambda^k\cap H^*_0\Lambda^k$ and on the simplicial subdivisions of the domain. The finite element scheme possesses an $\mathcal{O}(h)$-order convergence rate for sufficiently regular data, and an $\mathcal{O}(h^s)$-order rate on any $s$-regular domain, $0<s\leqslant 1$, no matter what topology the domain has.
Numerical Analysis
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to solve the numerical solution problem of the original weak form of the n - dimensional Hodge - Laplace equation on domains with arbitrary topological structures. Specifically, the author proposes a class of non - conforming finite element schemes to deal with this problem, and these schemes are applicable to differential forms of any dimension \( n \geq 2 \) and any order \( 1 \leq k \leq n - 1 \). #### Background and motivation The Hodge - Laplace equation is a core problem in many applied sciences, including fields such as electromagnetics and fluid - structure interactions. Traditional conforming finite element methods may lead to singular solutions or inaccurate results, especially when the solution is not smooth enough. To solve this problem, the author proposes a new non - conforming finite element method to ensure that correct approximation results can be obtained even when the solution is not smooth enough. #### Main challenges 1. **Singular solution problem**: For non - smooth solutions, traditional conforming finite element methods may converge to the wrong limit. 2. **Topological complexity**: It is necessary to consider the properties of solutions under different topological structures, especially when the solution contains non - trivial harmonic forms. 3. **Numerical stability**: Ensure the stability and convergence of the numerical scheme, especially when dealing with complex geometric shapes and boundary conditions. #### Proposed method The author proposes a new non - conforming finite element space, with the following main features: - **Local basis functions**: All finite element functions are defined by local shape function spaces and satisfy specific continuity conditions. - **No penalty or stabilization terms**: There is no need to introduce penalty or stabilization terms, thus simplifying the calculation process. - **Discrete Poincaré inequality**: The discrete Poincaré inequality is proved through the partial adjoint operator theory, which is crucial for dealing with non - trivial topological structures. #### Results and contributions This method has the following advantages: - For sufficiently smooth data, the convergence rate of the scheme is \( O(h) \). - For any \( s \)-regular domain (\( 0 < s \leq 1 \)), the convergence rate of the scheme is \( O(h^s) \), regardless of the topological structure. - Through the indirect error estimation method, the effectiveness and accuracy of this method are verified. In conclusion, this paper provides a new non - conforming finite element method that can effectively deal with the numerical solution problem of the Hodge - Laplace equation, especially suitable for cases with complex topological structures and non - smooth solutions.