Evaluation of Inner Products of Implicitly-defined Finite Element Functions on Multiply Connected Planar Mesh Cells
Jeffrey S. Ovall,Samuel E. Reynolds
2023-03-14
Abstract:Recent advancements in finite element methods allows for the implementation of mesh cells with curved edges. In the present work, we develop the tools necessary to employ multiply connected mesh cells, i.e. cells with holes, in planar domains. Our focus is efficient evaluation the $H^1$ semi-inner product and $L^2$ inner product of implicitly-defined finite element functions of the type arising in boundary element based finite element methods (BEM-FEM) and virtual element methods (VEM). These functions may be defined by specifying a polynomial Laplacian and a continuous Dirichlet trace. We demonstrate that these volumetric integrals can be reduced to integrals along the boundaries of mesh cells, thereby avoiding the need to perform any computations in cell interiors. The dominating cost of this reduction is solving a relatively small Nystrom system to obtain a Dirichlet-to-Neumann map, as well as the solution of two more Nystrom systems to obtain an ``anti-Laplacian'' of a harmonic function, which is used for computing the $L^2$ inner product. We demonstrate that high-order accuracy can be achieved with several numerical examples.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to effectively calculate the \(H^1\) semi - inner product and the \(L^2\) inner product of implicitly defined finite - element functions in multiply - connected planar grid cells (i.e., grid cells with holes). Specifically:
1. **Background and Motivation**:
- Recent progress in the finite - element method has made it possible to implement grid cells with curved boundaries.
- This paper aims to develop the necessary tools to use multiply - connected grid cells (i.e., cells with holes) in planar domains, especially for implicitly defined finite - element functions that occur in the boundary element method (BEM - FEM) and the virtual element method (VEM).
2. **Core of the Problem**:
- Calculate the \(H^1\) semi - inner product \(\int_K \nabla v\cdot\nabla w \, dx\) and the \(L^2\) inner product \(\int_K vw \, dx\) of implicitly defined finite - element functions in multiply - connected grid cells.
- These functions are defined by specifying a polynomial Laplace operator and continuous Dirichlet traces.
3. **Key Points of the Solution**:
- Simplify these volume integrals into integrals along the boundaries of the grid cells, thus avoiding any calculations inside the cells.
- The main computational cost lies in solving the relatively small Nyström system to obtain the Dirichlet - to - Neumann mapping and solving two other Nyström systems to obtain the "inverse Laplace" of harmonic functions.
- High - order accuracy is demonstrated through numerical examples.
4. **Specific Steps**:
- **\(H^1\) Semi - Inner Product**: Transform the volume integral into a boundary integral and use the Dirichlet - to - Neumann mapping to calculate the normal derivative.
- **\(L^2\) Inner Product**: Similarly, transform the volume integral into a boundary integral and use the "inverse Laplace" of harmonic functions for calculation.
- **Handling Multiply - Connectedness**: For grid cells with holes, adjust the calculation methods of the Dirichlet - to - Neumann mapping and the "inverse Laplace".
5. **Conclusion**:
- Through the above methods, the \(H^1\) semi - inner product and the \(L^2\) inner product in multiply - connected grid cells can be effectively calculated without the need for complex two - dimensional integrations inside the cells.
In summary, this paper solves the problem of efficiently calculating the inner products of finite - element functions in complex geometries (such as grid cells with holes), providing theoretical and technical support for a wider range of applications.