A hybrid approach to solve the high-frequency Helmholtz equation with source singularity in smooth heterogeneous media

From MaRDI portal
Publication:2425268

DOI10.1016/J.JCP.2018.03.011zbMATH Open1415.65255arXiv1710.02307OpenAlexW2963063923MaRDI QIDQ2425268FDOQ2425268

Jun Fang, Hongkai Zhao, Leonardo Zepeda-Núñez, Jianliang Qian

Publication date: 26 June 2019

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

Abstract: We propose a hybrid approach to solve the high-frequency Helmholtz equation with point source terms in smooth heterogeneous media. The method is based on the ray-based finite element method (ray-FEM), whose original version can not handle the singularity close to point sources accurately. This pitfall is addressed by combining the ray-FEM, which is used to compute the smooth far-field of the solution accurately, with a high-order asymptotic expansion close to the point source, which is used to properly capture the singularity of the solution in the near-field. The method requires a fixed number of grid points per wavelength to accurately represent the wave field with an asymptotic convergence rate of mathcalO(omega1/2), where omega is the frequency parameter in the Helmholtz equation. In addition, a fast sweeping-type preconditioner is used to solve the resulting linear system. We present numerical examples in 2D to show both accuracy and efficiency of our method as the frequency increases. In particular, we provide numerical evidence of the convergence rate, and we show empirically that the overall complexity is mathcalO(omega2) up to a poly-logarithmic factor.


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




Recommendations




Cites Work


Cited In (9)

Uses Software





This page was built for publication: A hybrid approach to solve the high-frequency Helmholtz equation with source singularity in smooth heterogeneous media

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