Pollution-free and fast hybridizable discontinuous Galerkin solvers for the high-frequency Helmholtz equation

• Main Authors: Núñez, Leonardo Zepeda (Author), Hewett, Russell (Author), Taus, Matthias F (Contributor), Demanet, Laurent (Contributor)
• Other Authors: Massachusetts Institute of Technology. Department of Mathematics (Contributor)
• Format: Article
• Language: English
• Published: Society of Exploration Geophysicists, 2018-06-12T14:44:56Z.
• Subjects: Article
• Online Access: Get fulltext

In this work we propose a hybridizable discontinuous Galerkin (hdG) discretization of the high-frequency Helmholtz equation in the presence of point sources and highly heterogeneous and discontinuous wave speed models. We show that it delivers solutions that are provably second-order accurate and do not suffer from the pollution error, as long as a slightly higher order hdG method is used where the polynomial degree is chosen such that p = O(logw). These results hold even if the discontinuities in the wave speed are not resolved by the hdG mesh, as long as the integration procedure used in the assembly of the stiffness matrix respects the discontinuities. Further, we show that the associated linear systems can be solved using a modification of the method of polarized traces resulting in a method with linear complexity up to a poly-logarithmic factor, or sub-linear complexity in a parallel environment. To our knowledge and surprise, this note contains the first instance of a numerical method that is at the same time fast (O(N) runtime) and accurate (second-order, pollution-free) in the context of models of geophysical interest. Keywords: finite element, frequency-domain, numerical, wave equation, acoustic
TOTAL (Firm)
United States. Air Force. Office of Scientific Research (Grant FA9550-12-1-0328)
United States. Air Force. Office of Scientific Research (Grant FA9550-15-1-0078)
United States. Office of Naval Research (Grant N00014-16-1-2122)
National Science Foundation (U.S.) (Grant DMS-1255203)