Preconditioning the 2D Helmholtz equation with polarized traces
Name
polarized_traces.pdf
Size
1.27 MB
Format
Adobe PDF
Checksum (MD5)
a4a6eded57d3e7f8ad61b96133a17bfd
Author(s)
Zepeda Nunez, Leonardo Andres
Hewett, Russell
Demanet, Laurent
Date Issued
2014
Journal
SEG Technical Program Expanded Abstracts 2014
Publisher
Society of Exploration Geophysicists
Citation
Zepeda-Núñez*, Leonardo, Russell J. Hewett, and Laurent Demanet. “Preconditioning the 2D Helmholtz Equation with Polarized Traces.” SEG Technical Program Expanded Abstracts 2014 (August 5, 2014).
Version
Author's final manuscript
Abstract
We present a domain decomposition solver for the 2D Helmholtz equation, with a special choice of integral transmission condition that involves polarizing the waves into oneway components. This refinement of the transmission condition is the key to combining local direct solves into an efficient iterative scheme, which can then be deployed in a highperformance computing environment. The method involves an expensive, but embarrassingly parallel precomputation of local Green's functions, and a fast online computation of layer potentials in partitioned low-rank form. The online part has sequential complexity that scales sublinearly with respect to the number of volume unknowns, even in the high-frequency regime. The favorable complexity scaling continues to hold in the context of low-order finite difference schemes for standard community models such as BP and Marmousi2, where convergence occurs in 5 to 10 GMRES iterations.
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/SEGAM2014-1275.1
