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

From MaRDI portal
Publication:2375151

DOI10.1016/J.JCP.2016.03.023zbMATH Open1349.65571arXiv1504.01609OpenAlexW2284793936MaRDI QIDQ2375151FDOQ2375151


Authors: Christiaan C. Stolk Edit this on Wikidata


Publication date: 5 December 2016

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

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 hk, 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 108 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.


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




Recommendations




Cites Work


Cited In (21)

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)