Hybrid weakly over-penalised symmetric interior penalty method on anisotropic meshes

Hiroki Ishizaka

2024-06-12

Abstract:In this study, we investigate a hybrid-type anisotropic weakly over-penalised symmetric interior penalty method for the Poisson equation on convex domains. Compared with the well-known hybrid discontinuous Galerkin methods, our approach is simple and easy to implement. Our primary contributions are the proposal of a new scheme and the demonstration of a proof for the consistency term, which allows us to estimate the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart--Thomas finite element space and a discontinuous space. In numerical experiments, we compare the calculation results for standard and anisotropic mesh partitions.

Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the problem of how to effectively apply the Hybrid Weakly Over - Penalised Symmetric Interior Penalty (HWOPSIP) method when solving the Poisson equation on anisotropic meshes. Specifically, the main research problems include:

1. **Propose a new scheme**: Introduce a new HWOPSIP method, which combines the relationship between the Raviart - Thomas (RT) finite - element space and the discontinuous space, thereby simplifying the implementation and improving the computational efficiency.

2. **Prove the consistency term**: Provide a proof of the consistency term for the new method, which is crucial for estimating the anisotropic consistency error. The key lies in using the relationship between RT finite - element interpolation and the discontinuous space to derive the optimal consistency error estimate.

3. **Numerical experiment comparison**: Verify the effectiveness and superiority of the new method by comparing the computational results under standard meshing and anisotropic meshing through numerical experiments.

4. **Stability and error estimate**: Analyze the stability and error estimate of the proposed HWOPSIP method in the energy norm and \(L^2\) norm, and pay special attention to the performance under the semi - regular condition (i.e., the maximum - angle condition).

5. **Handle the challenges of anisotropic meshes**: In view of the characteristics of anisotropic meshes, solve the problem that the optimal estimate cannot be obtained by using the trace inequality, and ensure the computational accuracy and stability on these complex meshes.

In summary, this research aims to develop an efficient and stable HWOPSIP method suitable for anisotropic meshes and verify its performance through theoretical analysis and numerical experiments. This helps to provide more accurate and reliable numerical solutions when dealing with physical problems with anisotropic behavior.