Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4967409 | Journal of Computational Physics | 2017 | 15 Pages |
Abstract
We develop a frequency-domain iterative solver for numerical simulation of acoustic waves in 3D heterogeneous media. It is based on the application of a unique preconditioner to the Helmholtz equation that ensures convergence for Krylov subspace iteration methods. Effective inversion of the preconditioner involves the Fast Fourier Transform (FFT) and numerical solution of a series of boundary value problems for ordinary differential equations. Matrix-by-vector multiplication for iterative inversion of the preconditioned matrix involves inversion of the preconditioner and pointwise multiplication of grid functions. Our solver has been verified by benchmarking against exact solutions and a time-domain solver.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Mikhail Belonosov, Maxim Dmitriev, Victor Kostin, Dmitry Neklyudov, Vladimir Tcheverda,