1887

Abstract

Summary

We propose a preconditioned iterative method for solving the Helmholtz equation in heterogeneous media. Our method is based on Krylov type linear solvers, similarly to several other iterative solver approaches. The distinctive feature of our method is the use of multilevel preconditioner. Firstly, as a preconditioner we suggest to use complex damped Helmholtz operator for vertically-inhomogeneous media. We represent the initial problem as a perturbation of a preconditioner. As a result, a matrix-by-vector multiplication of the preconditioned system may be evaluated via 2D FFT in x and y directions followed by solution of a number of small banded systems of linear equations in z directions. A second level preconditioner based on Neumann series decomposition provides fast convergence in complex situation. After decomposition of the computational domain along one lateral direction we parallelize our algorithm to be used on high performance computer systems with hybrid architecture (MPI+OpenMP).

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20140266
2014-04-07
2024-04-26
Loading full text...

Full text loading...

References

  1. Colella, P., Bell, J., Keen, N., Ligocki, T., Lijewski, M. and van Straalen, B.
    [2007] Performance and scaling of locally-structured grid methods for partial differential equations. Journal of Physics: Conference Series, 78 012013.
    [Google Scholar]
  2. Engquist, B. and Ying, L.
    [2011] Sweeping preconditioner for the Helmholtz equation: moving perfectily matched layers, Multiscale Model. Simul., 9(2), 686–710.
    [Google Scholar]
  3. Erlangga, Y.A. [2008] Advances in iterative methods and preconditioners for the Helmholtz equation. Arch. Comput. Methods Eng., 15, 37–66.
    [Google Scholar]
  4. Neklyudov, D. and Tcheverda, V.
    [2010a] A Helmholtz iterative solver with semi-analytical preconditioner for seismic imaging and inversion. International Conference & Exhibition, Saint Petersburg-2010, extended abstract, P121.
    [Google Scholar]
  5. [2010b] A Helmholtz Iterative Solver without of Finite– difference Approximations. 72nd EAGE Conference and Exhibition, Barcelona, Spain, 2010, extended abstracts, G006.
    [Google Scholar]
  6. Plessix, R.E.
    [2007] A Helmholtz iterative solver for 3D seismic-imaging problems. Geophysics, 72(5), SM185–SM194.
    [Google Scholar]
  7. Trefethen, L. and Bau, D.
    [1997] Numerical Linear Algebra. SIAM.
    [Google Scholar]
  8. Vainberg, B.R.
    [1966] Principles of radiation, vanishing attenuation and limit amplitude in the general theory of partial differential equations. Usp. Mat. Nauk, 21 3 (129), 115–194. (in Russian).
    [Google Scholar]
  9. Virieux, J. and Operto, S.
    [2009] An overview of full-waveform inversion in exploration geophysics. Geophysics, 74(6), WCC1–WCC26.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20140266
Loading
/content/papers/10.3997/2214-4609.20140266
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error