Sixth-Order Compact Finite Difference Method for 2D Helmholtz Equations with Singular Sources and Reduced Pollution Effect

From MaRDI portal
Publication:6142995

DOI10.4208/CICP.OA-2023-0062arXiv2112.07154OpenAlexW4387393792MaRDI QIDQ6142995FDOQ6142995

Author name not available (Why is that?)

Publication date: 4 January 2024

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

Abstract: Due to its highly oscillating solution, the Helmholtz equation is numerically challenging to solve. To obtain a reasonable solution, a mesh size that is much smaller than the reciprocal of the wavenumber is typically required (known as the pollution effect). High order schemes are desirable, because they are better in mitigating the pollution effect. In this paper, we present a high order compact finite difference method for 2D Helmholtz equations with singular sources, which can also handle any possible combinations of boundary conditions (Dirichlet, Neumann, and impedance) on a rectangular domain. Our method achieves a sixth order consistency for a constant wavenumber, and a fifth order consistency for a piecewise constant wavenumber. To reduce the pollution effect, we propose a new pollution minimization strategy that is based on the average truncation error of plane waves. Our numerical experiments demonstrate the superiority of our proposed finite difference scheme with reduced pollution effect to several state-of-the-art finite difference schemes, particularly in the critical pre-asymptotic region where extsfkh is near 1 with extsfk being the wavenumber and h the mesh size.


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




Recommendations




Cites Work


Cited In (5)





This page was built for publication: Sixth-Order Compact Finite Difference Method for 2D Helmholtz Equations with Singular Sources and Reduced Pollution Effect

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