In this paper we investigate staggered discontinuous Galerkin method for the Helmholtz equation with large wave number on general polygonal meshes. The method is highly flexible by allowing rough grids such as the trapezoidal grids and highly distorted grids, and at the same time, is numerical flux free. Furthermore, it allows hanging nodes, which can be simply treated as additional vertices. By exploiting a modified duality argument, the stability and convergence can be proved under the condition that [Formula presented], where m is the polynomial order, κ is the wave number, h is the mesh size and C0 is a positive constant independent of κ,h. Error estimates for both the scalar and vector variables in L2 norm are established. Several numerical experiments are tested to verify our theoretical results and to present the capability of our method for capturing singular solutions.
Bibliographical noteFunding Information:
The research of Eric T. Chung is partially supported by the Hong Kong RGC General Research Fund (Project numbers 14304217 and 14302018 ) and CUHK Faculty of Science Direct Grant 2019–20. The research of Eun-Jae Park was supported by the National Research Foundation of Korea (NRF) grant funded by the Ministry of Science and ICT ( NRF-2015R1A5A1009350 and NRF-2019R1H1A2080222 ).
© 2020 Elsevier Ltd
All Science Journal Classification (ASJC) codes
- Modelling and Simulation
- Computational Theory and Mathematics
- Computational Mathematics