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

C. D. Riyanti, A. Kononov, Y. A. Erlangga, C. Vuik, C. W. Oosterlee, R. E. Plessix, W. A. Mulder

Research output: Contribution to journalArticlepeer-review

61 Citations (Scopus)

Abstract

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.

Original languageEnglish
Pages (from-to)431-448
Number of pages18
JournalJournal of Computational Physics
Volume224
Issue number1
DOIs
Publication statusPublished - May 20 2007

Keywords

  • Helmholtz equation
  • Krylov subspace method
  • Multigrid method
  • Preconditioner

ASJC Scopus subject areas

  • Numerical Analysis
  • Modelling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation'. Together they form a unique fingerprint.

Cite this