Abstract:We propose a new geometrically unfitted finite element method based on discontinuous Trefftz ansatz spaces. Trefftz methods allow for a reduction in the number of degrees of freedom in discontinuous Galerkin methods, thereby, the costs for solving arising linear systems significantly. This work shows that they are also an excellent way to reduce the number of degrees of freedom in an unfitted setting. We present a unified analysis of a class of geometrically unfitted discontinuous Galerkin methods with different stabilisation mechanisms to deal with small cuts between the geometry and the mesh. We cover stability and derive a-priori error bounds, including errors arising from geometry approximation for the class of discretisations for a model Poisson problem in a unified manner. The analysis covers Trefftz and full polynomial ansatz spaces, alike. Numerical examples validate the theoretical findings and demonstrate the potential of the approach.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to propose and analyze a new unfitted finite element method based on the discontinuous Trefftz basis function space. Specifically, the authors hope to solve the following key problems:
1. **Reduce the degrees of freedom**: The traditional discontinuous Galerkin (DG) method introduces a large number of degrees of freedom due to its discretization method, which significantly increases the computational cost of solving the linear system. By using Trefftz basis functions, the number of degrees of freedom can be reduced, thereby reducing the computational cost.
2. **Handle complex geometries**: When dealing with problems with complex geometries, traditional body - fitted methods require frequent mesh generation, adaptation, and re - partitioning, which increases the computational burden. The geometric unfitted method can avoid these problems by embedding the computational domain into the background mesh.
3. **Stability and error estimation**: For geometric unfitted methods, when the boundary cuts through the mesh cells, poorly - shaped cells may be generated, which will affect the numerical stability. To this end, the author introduced different stabilization mechanisms (such as ghost penalty and element aggregation) to ensure the stability of the numerical method and derived a priori error estimates, including geometric approximation errors.
4. **Unified analysis framework**: The author provides a unified analysis framework that covers different types of geometric unfitted discontinuous Galerkin methods, including methods using full - polynomial basis functions and Trefftz basis functions, and combines different stabilization mechanisms.
### Main contributions of the paper
- **Describe geometric unfitted DG and Trefftz DG methods**: Combine three different stabilization mechanisms (ghost penalty, patch - wise ghost penalty, and element aggregation) and analyze them with and without geometric approximation errors.
- **Unified a priori error estimates**: Provide a priori error estimates including geometric errors for these geometric unfitted (Trefftz) DG methods.
- **Numerical experiments and implementation details discussion**: Present the numerical experiment results of these methods and discuss the implementation details.
### Model problem
As a model boundary value problem, the author considers the Poisson equation defined on the open bounded region \(\Omega\subseteq\mathbb{R}^d\) (\(d = 2,3\)):
\[
\begin{cases}
-\Delta u=f&\text{in }\Omega,\\
u = g&\text{on }\Gamma,
\end{cases}
\]
where \(f\in L^2(\Omega)\) and \(g\in H^{1/2}(\Gamma)\) are given functions, and \(\Gamma:=\partial\Omega\) is the boundary of the region \(\Omega\).
### Summary
This paper mainly focuses on developing and analyzing geometric unfitted Trefftz discontinuous Galerkin methods, aiming to reduce the degrees of freedom, improve computational efficiency, handle complex geometries, and ensure the stability and accuracy of the numerical method.