Abstract:This paper presents an efficient weak Galerkin (WG) finite element method with reduced stabilizers for solving the time-harmonic Maxwell equations on both convex and non-convex polyhedral meshes. By employing bubble functions as a critical analytical tool, the proposed method enhances efficiency by partially eliminating the stabilizers traditionally used in WG methods. This streamlined WG method demonstrates stability and effectiveness on convex and non-convex polyhedral meshes, representing a significant improvement over existing stabilizer-free WG methods, which are typically limited to convex elements within finite element partitions. The method achieves an optimal error estimate for the exact solution in a discrete $H^1$ norm, and additionally, an optimal $L^2$ error estimate is established for the WG solution. Several numerical experiments are conducted to validate the method's efficiency and accuracy.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: to develop an efficient weak Galerkin finite element method (WG method) applicable to non - convex polyhedral meshes for solving Maxwell equations without convexity constraints. Specifically, this method aims to reduce the dependence on stabilization terms and, by using bubble functions as a key analytical tool, enable the method to be applied to convex and non - convex polyhedral meshes.
### Core problems of the paper
1. **Reduce stabilization terms**:
- In the traditional weak Galerkin method, stabilization terms are necessary, but they increase computational complexity and programming difficulty.
- This paper proposes a new method. By using high - degree polynomials to calculate the discrete weak curl, it partially eliminates the need for stabilization terms while retaining the stabilization terms of the discrete weak gradient.
2. **Applicable to non - convex polyhedral meshes**:
- Existing weak Galerkin methods without stabilization terms are usually limited to convex elements.
- This method bypasses the convexity limitation by introducing bubble functions as a key tool, making the method applicable to a wider range of geometric shapes, including non - convex polyhedral meshes.
3. **Theoretical and numerical verification**:
- It is theoretically proven that this method has the optimal error estimate in the discrete H1 norm and the L2 error estimate is established.
- The effectiveness and accuracy of this method are verified through a series of numerical experiments.
### Background of Maxwell equations
Maxwell equations describe the behavior of electromagnetic fields. Their time - harmonic forms are:
\[
\nabla \times E=-i\omega \mu H,
\]
\[
\nabla \times H = i\omega \epsilon E+\sigma E,
\]
\[
\nabla \cdot (\epsilon E)=\rho,
\]
\[
\nabla \cdot (\mu H) = 0.
\]
where:
- \(E\) represents the electric field strength,
- \(H\) represents the magnetic field strength,
- \(D = \epsilon E\) represents the displacement current density,
- \(B=\mu H\) represents the magnetic flux density,
- \(j=\sigma E\) represents the current density,
- \(\mu\) is the magnetic permeability,
- \(\sigma\) is the electrical conductivity,
- \(\rho\) is the charge density,
- \(\epsilon\) is the dielectric constant.
### Main contributions of the method
1. **Reduce stabilization terms**:
- By using high - degree polynomials to calculate the discrete weak curl, one stabilization term is eliminated while the stabilization term of the discrete weak gradient is retained, thus reducing the size of the stiffness matrix and global sparsity and decreasing the programming complexity.
2. **Applicable to non - convex polyhedral meshes**:
- Using bubble functions as a key tool enables the method to handle non - convex polyhedral meshes, expanding the application range.
3. **Theoretical analysis and numerical verification**:
- The optimal error estimate in the discrete H1 norm and the L2 error estimate are theoretically proven.
- The effectiveness and accuracy of the method are verified through numerical experiments.
In conclusion, this paper proposes an efficient weak Galerkin finite element method, which not only reduces the dependence on stabilization terms but also expands the applicable mesh types, significantly improving the efficiency and flexibility of solving Maxwell equations.