A Primal Staggered Discontinuous Galerkin Method on Polytopal Meshes

L. Chen,X. Huang,E. Park,R. Wang

2024-10-31

Abstract:This paper introduces a novel staggered discontinuous Galerkin (SDG) method tailored for solving elliptic equations on polytopal meshes. Our approach utilizes a primal-dual grid framework to ensure local conservation of fluxes, significantly improving stability and accuracy. The method is hybridizable and reduces the degrees of freedom compared to existing approaches. It also bridges connections to other numerical methods on polytopal meshes. Numerical experiments validate the method's optimal convergence rates and computational efficiency.

Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the problem of efficiently solving elliptic equations on polyhedral meshes. Specifically, the author proposes a new Staggered Discontinuous Galerkin (SDG) method, which utilizes the primal - dual mesh framework to ensure local flux conservation, thereby significantly improving stability and accuracy. The following are the main contributions of the paper:

1. **New discretization scheme**:
- This method uses the primal - dual mesh framework on polyhedral meshes. The primal scalar variable is continuous on the polyhedral mesh, while the flux vector variable is continuous on the face - based staggered dual mesh.
- Specifically, the primal variable is in the form of a single polynomial on each polyhedron, which significantly reduces the number of degrees of freedom.

2. **Hybridization**:
- By introducing Lagrange multipliers to enforce the continuity of the normal flux on the primal face, the flux can be eliminated element - by - element, thus simplifying the saddle - point system into a symmetric positive - definite system.
- This method can be regarded as the discretization of the primal form of the weak gradient.

3. **Connection with other methods**:
- For the case of the lowest order \(k = 0\) (i.e., the \(P_0 - P_1\) type), after barycentric refinement on the triangular mesh, the resulting matrix is the same as the Crouzeix - Raviart non - conforming finite - element matrix, except for the right - hand side.
- The hybridized \(P_0 - P_1\) SDG can be regarded as the generalization of the linear non - conforming finite element on general polyhedra.
- The primal SDG method in the hybridized version can be regarded as the unstabilized weak Galerkin method or the unstabilized non - conforming virtual element method for the primal form of the Poisson equation.

4. **Local conservation**:
- The standard continuous Galerkin method lacks a locally conservative flux approximation, which is crucial for applications such as streamline construction and coupled flow - transport simulations.
- This method retains the local conservation of the SDG framework, ensuring mass conservation on each primal polyhedral unit while improving stability and accuracy.

5. **Numerical experiments**:
- The author conducted numerical experiments to verify the optimal convergence rate and computational efficiency of this method, demonstrating its advantages in terms of convergence speed and computational efficiency.

In conclusion, this paper proposes a novel Staggered Discontinuous Galerkin method, aiming to solve the problem of solving elliptic equations on polyhedral meshes. It ensures local conservation through the primal - dual mesh framework and improves computational efficiency by reducing the number of degrees of freedom through the hybridization technique.