Abstract:We develop a cut Discontinuous Galerkin method (cutDGM) for a diffusion-reaction equation in a bulk domain which is coupled to a corresponding equation on the boundary of the bulk domain. The bulk domain is embedded into a structured, unfitted background mesh. By adding certain stabilization terms to the discrete variational formulation of the coupled bulk-surface problem, the resulting cutDGM is provably stable and exhibits optimal convergence properties as demon- strated by numerical experiments. We also show both theoretically and numerically that the system matrix is well-conditioned, irrespective of the relative position of the bulk domain in the background mesh.
What problem does this paper attempt to address?
This paper attempts to solve the numerical solution problem of coupled body - surface partial differential equations (PDEs), especially when there is an interaction between the body domain and the surface domain on its boundary. Specifically, the author proposes a new Cut Discontinuous Galerkin Method (cutDGM) for discretizing the coupled body - surface problem on a given bounded domain Ω. The main challenges of this method include:
1. **Coupled PDE systems of different topological dimensions**: The coupling between the body domain and the surface domain requires that the numerical method can handle domains of different dimensions.
2. **Fast and robust mesh generation for complex geometries**: Extremely complex surface geometries occur in many practical applications, such as the complex fracture networks in porous media models.
3. **Geometries evolving over time**: Even if the initial geometry is relatively simple, it may change significantly over time and even undergo topological changes.
To address these challenges, the author introduces the idea of the Cut Finite Element Method (CutFEM), embeds the geometry into a fixed background mesh, and adds certain stabilization terms to ensure the stability of the method regardless of the position of the geometry in the background mesh. These stabilization terms provide control over the local variation of the discrete function, thus ensuring the well - conditionedness and optimal convergence of the system.
### Specific problem description
The paper aims to solve the following problems:
- **Develop a new numerical method**: Propose a Cut Discontinuous Galerkin Method (cutDGM) for discretizing the coupled body - surface PDE system.
- **Ensure the stability of the numerical method**: By adding appropriate stabilization terms, ensure that the method remains stable under different geometric configurations.
- **Prove optimal convergence**: Through theoretical analysis and numerical experiments, prove that the proposed method has optimal convergence properties.
- **Ensure the well - conditionedness of the system matrix**: Regardless of the position of the geometry in the background mesh, the system matrix should remain well - conditioned.
### Mathematical formula representation
The mathematical formulas involved in the paper are represented in Markdown format, for example:
- Laplace - Beltrami operator:
\[
\Delta_\Gamma=\nabla_\Gamma\cdot\nabla_\Gamma
\]
- Bilinear form in weak form:
\[
a(u, v)=c_\Omega a_\Omega(u_\Omega, v_\Omega)+c_\Gamma a_\Gamma(u_\Gamma, v_\Gamma)+a_{\Omega\Gamma}(u, v)
\]
Through these methods, the author has successfully solved several key problems in the numerical solution of the coupled body - surface PDE system and provided a solid foundation for future related research.