Absolutely stable local discontinuous Galerkin methods for the Helmholtz equation with large wave number

Abstract:

This paper develops and analyzes two local discontinuous Galerkin (LDG) methods using piecewise linear polynomials for the Helmholtz equation with the first order absorbing boundary condition in the high frequency regime. It is shown that the proposed LDG methods are stable for all positive wave number $k$ and all positive mesh size $h$. Energy norm and $L^2$-norm error estimates are derived for both LDG methods in all mesh parameter regimes including pre-asymptotic regime (i.e., $k^2 h\gtrsim 1$). To analyze the proposed LDG methods, they are recast and treated as (nonconforming) mixed finite element methods. The crux of the analysis is to show that the sesquilinear form associated with each LDG method satisfies a coercivity property in all mesh parameter regimes. These coercivity properties then easily infer the desired discrete stability estimates for the solutions of the proposed LDG methods. In return, the discrete stabilities not only guarantee the well-posedness of the LDG methods but also play a crucial role in the error analysis. Numerical experiments are also presented in the paper to validate the theoretical results and to compare the performance of the proposed two LDG methods.