A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
arXiv:1806.11277 · doi:10.1016/j.jcp.2023.112622
Abstract
In this work we extend the shifted Laplacian approach to the elastic Helmholtz equation. The shifted Laplacian multigrid method is a common preconditioning approach for the discretized acoustic Helmholtz equation. In some cases, like geophysical seismic imaging, one needs to consider the elastic Helmholtz equation, which is harder to solve: it is three times larger and contains a nullity-rich grad-div term. These properties make the solution of the equation more difficult for multigrid solvers. The key idea in this work is combining the shifted Laplacian with approaches for linear elasticity. We provide local Fourier analysis and numerical evidence that the convergence rate of our method is independent of the Poisson's ratio. Moreover, to better handle the problem size, we complement our multigrid method with the domain decomposition approach, which works in synergy with the local nature of the shifted Laplacian, so we enjoy the advantages of both methods without sacrificing performance. We demonstrate the efficiency of our solver on 2D and 3D problems in heterogeneous media.
Published in the Journal of Computational Physics
References in corpus (4)
- A convergent Born series for solving the inhomogeneous Helmholtz equation in arbitrarily large media
- L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation
- Fast Multipole Method For 3-D Helmholtz Equation In Layered Media
- Multidirectionnal sweeping preconditioners with non-overlapping checkerboard domain decomposition for Helmholtz problems