A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation

Journal Article · Sun May 20 00:00:00 EDT 2007 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2007.03.033· OSTI ID:20991585

Riyanti, C D ^[1]; Kononov, A ^[2]; Erlangga, Y A ^[3]; Vuik, C ^[1]; Oosterlee, C W ^[1]; Plessix, R -E ^[4]; Mulder, W A ^[4]

---

1. Delft Institute of Applied Mathematics, Delft University of Technology, Mekelweg 4, 2628 CD Delft (Netherlands)
2. Computational Physics Group, Delft University of Technology (Netherlands)
3. Scientific Computing, Technical University, Berlin (Germany)
4. Shell International E and P, P.O. Box 60, 2280 AB Rijswijk (Netherlands)

We investigate the parallel performance of an iterative solver for 3D heterogeneous Helmholtz problems related to applications in seismic wave propagation. For large 3D problems, the computation is no longer feasible on a single processor, and the memory requirements increase rapidly. Therefore, parallelization of the solver is needed. We employ a complex shifted-Laplace preconditioner combined with the Bi-CGSTAB iterative method and use a multigrid method to approximate the inverse of the resulting preconditioning operator. A 3D multigrid method with 2D semi-coarsening is employed. We show numerical results for large problems arising in geophysical applications.

• OSTI ID: 20991585
• Journal Information: Journal of Computational Physics, Vol. 224, Issue 1; Other Information: DOI: 10.1016/j.jcp.2007.03.033; PII: S0021-9991(07)00147-7; Copyright (c) 2007 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
EQUATIONS
HELMHOLTZ THEOREM
ITERATIVE METHODS
PERFORMANCE
SEISMIC WAVES
THREE-DIMENSIONAL CALCULATIONS
TWO-DIMENSIONAL CALCULATIONS
WAVE PROPAGATION