A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory

DOI10.1016/J.JCP.2016.03.023MaRDI QIDQ2375151zbMATHOpenAlexFDO

Publication date 5 December 2016

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1504.01609

domain decomposition finite difference method Helmholtz equation multigrid numerical dispersion

Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Finite difference methods for boundary value problems involving PDEs (65N06)

Abstract: We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a discrete operator to be applied to the source and the wavefields are constructed. Their coefficients are piecewise polynomial functions of

h k

, chosen such that phase and amplitude errors are minimal. The phase errors of the scheme are very small, approximately as small as those of the 2-D quasi-stabilized FEM method and substantially smaller than those of alternatives in 3-D, assuming the same number of gridpoints per wavelength is used. In numerical experiments, accurate solutions are obtained in constant and smoothly varying media using meshes with only five to six points per wavelength and wave propagation over hundreds of wavelengths. When used as a coarse level discretization in a multigrid method the scheme can even be used with downto three points per wavelength. Tests on 3-D examples with up to

1 0^{8}

degrees of freedom show that with a recently developed hybrid solver, the use of coarser meshes can lead to corresponding savings in computation time, resulting in good simulation times compared to the literature.

Recommendations

Cites work

Cited in

(24)

Describes a project that uses

Uses Software

This page was built for publication: A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2375151)