A mixed quasi-trace surface finite element method for the Laplace-Beltrami problem

Alan Demlow
2023-09-22
Abstract:Trace finite element methods have become a popular option for solving surface partial differential equations, especially in problems where surface and bulk effects are coupled. In such methods a surface mesh is formed by approximately intersecting the continuous surface on which the PDE is posed with a three-dimensional (bulk) tetrahedral mesh. In classical $H^1$-conforming trace methods, the surface finite element space is obtained by restricting a bulk finite element space to the surface mesh. It is not clear how to carry out a similar procedure in order to obtain other important types of finite element spaces such as $H({\rm div})$-conforming spaces. Following previous work of Olshanskii, Reusken, and Xu on $H^1$-conforming methods, we develop a ``quasi-trace'' mixed method for the Laplace-Beltrami problem. The finite element mesh is taken to be the intersection of the surface with a regular tetrahedral bulk mesh as previously described, resulting in a surface triangulation that is highly unstructured and anisotropic but satisfies a classical maximum angle condition. The mixed method is then employed on this mesh. Optimal error estimates with respect to the bulk mesh size are proved along with superconvergent estimates for the projection of the scalar error and a postprocessed scalar approximation.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a mixed quasi - trace surface finite element method for the Laplace - Beltrami problem. Specifically, the author aims to deal with the following issues: 1. **Limitations of existing methods**: Traditional trace finite element methods perform well in solving surface partial differential equations, but these methods mainly focus on \(H^1\) - conforming spaces. For more complex function spaces such as \(H(\text{div})\) and \(H(\text{curl})\), how to construct appropriate finite element spaces remains a challenge. 2. **Need for mixed finite element methods**: For certain physical phenomena, such as two - phase flow problems, it is necessary to consider the partial differential equation systems on both the three - dimensional domain and the surface simultaneously. In this case, mixed finite element methods can provide better stability and accuracy. However, most of the existing mixed methods are based on standard finite element spaces and cannot be directly applied to surface problems. 3. **Impact of geometric errors**: When using a polyhedron to approximate the surface \(\gamma\), geometric errors will affect the convergence order of the solution. Therefore, how to effectively handle geometric errors in the mixed quasi - trace method is a key issue. 4. **Handling of unstructured grids**: Since the surface grids formed by the intersection of the surface and the volume grids are usually highly unstructured, this brings additional complexity to the solution process. How to construct an effective finite element method on such grids and ensure its stability is an important research direction. ### Main contributions of the paper - **Developing the mixed quasi - trace method**: The author proposes a new mixed quasi - trace surface finite element method, which is based on Raviart - Thomas (RT) and Brezzi - Douglas - Marini (BDM) elements and can work effectively on highly unstructured grids. - **Optimal error estimates**: Prove the optimal error estimates with respect to the volume grid size, and give super - convergence estimates for the projected scalar errors. - **Handling of geometric errors**: Analyze in detail the impact of geometric errors on the solution, and propose corresponding error estimates to ensure reliable numerical results even on irregularly - shaped grids. - **Numerical verification**: Verify the correctness of the theoretical analysis through numerical experiments, demonstrating the effectiveness and superiority of the new method. In summary, this paper overcomes the limitations of existing methods in dealing with complex surface partial differential equations by introducing the mixed quasi - trace surface finite element method, and provides new tools and ideas for research in related fields.