Superconvergence analysis of DG-FEM based on the polynomial preserving recovery for Helmholtz equation with high wave number
For computational scientists solving high-frequency Helmholtz problems, this provides rigorous superconvergence theory enabling more efficient error estimation and mesh adaptation.
The paper proves superconvergence of the linear discontinuous Galerkin finite element method with polynomial preserving recovery for the 2D Helmholtz equation, deriving error estimates with explicit dependence on wave number, penalty parameter, and mesh condition. Numerical examples confirm the theoretical results.
We study superconvergence property of the linear discontinuous Galerkin finite element method with the polynomial preserving recovery (PPR) and Richardson extrapolation for the two dimensional Helmholtz equation. The error estimate with explicit dependence on the wave number $k$, the penalty parameter $μ$ and the mesh condition parameter $α$ is derived. First, we prove that under the assumption $k(kh)^2\leq C_0$ ($h$ is the mesh size) and certain mesh condition, the estimate between the finite element solution and the linear interpolation of the exact solution is superconvergent under the $\norme{\cdot}$-seminorm. Second, we prove a superconvergence result for the recovered gradient by PPR. Furthermore, we estimate the error between the finite element gradient and recovered gradient, which motivate us to define the a posteriori error estimator. Finally, Some numerical examples are provided to confirm the theoretical results of superconvergence analysis. All theoretical findings are verified by numerical tests.