Sharp Penalty Term and Time Step Bounds for the Interior Penalty Discontinuous Galerkin Method for Linear Hyperbolic Problems

S. Geevers,J.J.W. van der Vegt

DOI: https://doi.org/10.1137/16M1091290

2018-07-06

Abstract:We present sharp and sufficient bounds for the interior penalty term and time step size to ensure stability of the Symmetric Interior Penalty Discontinuous Galerkin (SIPDG) method combined with an explicit time-stepping scheme. These conditions hold for generic meshes, including unstructured non-conforming heterogeneous meshes of mixed element types, and apply to a large class of linear hyperbolic problems, including the acoustic wave equation, the (an)isotropic elastic wave equations and Maxwell's equations. The penalty term bounds are computed element-wise, while bounds for the time step size are computed at weighted submeshes requiring only a small number of elements and faces. Numerical results illustrate the sharpness of these bounds.

Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve the problem of how to ensure numerical stability when the Symmetric Interior - Penalty Discontinuous Galerkin (SIPDG) method is used to solve linear hyperbolic problems. Specifically, the paper attempts to find strict and sufficient bounds for the interior penalty term and the time step to ensure the stability of the SIPDG method when using an explicit time - stepping scheme. These conditions are applicable to general grids, including unstructured, non - conforming heterogeneous grids with mixed cell types, and are applicable to a large class of linear hyperbolic problems, such as the acoustic wave equation, (isotropic or anisotropic) elastic wave equations, and Maxwell's equations.

### Core problems of the paper

1. **Boundary of the interior penalty term**: To ensure the stability of the numerical method, the interior penalty term needs to be large enough. If the penalty term is set too small, the SIPDG scheme will become unstable. Therefore, the paper derives the lower bound of the interior penalty term applicable to various grid types.

2. **Limit of the time step**: Similarly, to maintain numerical stability, the time step also needs to be small enough. If the time step is too large, it will also lead to numerical instability. Therefore, the paper also derives the upper bound of the time step.

3. **Applicable to multiple linear hyperbolic problems**: The method proposed in the paper is not only applicable to specific hyperbolic problems, but can also be applied to a wide range of linear hyperbolic problems, including but not limited to the acoustic wave equation, elastic wave equations, and Maxwell's equations.

4. **Applicable to complex grids**: The conditions proposed in the paper are applicable not only to regular grids, but also to unstructured, non - conforming heterogeneous grids, which makes the method more flexible and applicable when dealing with complex geometries and material interfaces.

### Summary of mathematical formulas

- **Lower bound of the penalty term**:
\[
\eta_e \geq c_\kappa \kappa^*_e
\]
where $\kappa^*_e$ is a constant obtained through matrix decomposition and eigenvalue calculation, and $c_\kappa \geq 1$ is an arbitrary constant.

- **Upper bound of the time step**:
\[
\Delta t \leq \frac{C_{\text{CFL}}}{\sqrt{\lambda_{\max}(A_h)}}
\]
where $C_{\text{CFL}}$ is the constant in the CFL condition, and $\lambda_{\max}(A_h)$ is the maximum eigenvalue of the discrete elliptic operator.

Through these strict mathematical derivations, the paper provides a reliable numerical stability guarantee for the SIPDG method when solving linear hyperbolic problems.