Pollution-free and fast hybridizable discontinuous Galerkin solvers for the high-frequency Helmholtz equation
Name
hdg-polarized-traces2.pdf
Size
420.56 KB
Format
Adobe PDF
Checksum (MD5)
dc8488dfe470ce3de48cb8d98a698d79
Author(s) • • •
Núñez, Leonardo Zepeda
Hewett, Russell
Taus, Matthias F
Demanet, Laurent
Date Issued
September 2017
Journal
SEG Technical Program Expanded Abstracts 2017
Publisher
Society of Exploration Geophysicists
Citation
Taus, Matthias, et al. "Pollution-Free and Fast Hybridizable Discontinuous Galerkin Solvers for the High-Frequency Helmholtz Equation." SEG Technical Program Expanded Abstracts 2017, pp. 4068–73.
Version
Author's final manuscript
Abstract
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
MIT Department
Massachusetts Institute of Technology. Department of Mathematics
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1190/SEGAM2017-17728116.1