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

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

Year:    2023

Author:    Qiwei Feng, Bin Han, Michelle Michelle

Communications in Computational Physics, Vol. 34 (2023), Iss. 3 : pp. 672–712

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 is sixth-order consistent for a constant wavenumber, and fifth-order consistent 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 $kh$ is near 1 with $k$ being the wavenumber and $h$ the mesh size.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.OA-2023-0062

Communications in Computational Physics, Vol. 34 (2023), Iss. 3 : pp. 672–712

Published online:    2023-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    41

Keywords:    Helmholtz equation finite difference pollution effect interface pollution minimization mixed boundary conditions corner treatment.

Author Details

Qiwei Feng

Bin Han

Michelle Michelle

  1. Fast High Order Algorithm for Three-Dimensional Helmholtz Equation Involving Impedance Boundary Condition with Large Wave Numbers

    Tong, Chengjie | Fang, Xianqi | Zhao, Meiling

    American Journal of Computational Mathematics, Vol. 13 (2023), Iss. 02 P.211

    https://doi.org/10.4236/ajcm.2023.132011 [Citations: 1]
  2. Sixth-order hybrid finite difference methods for elliptic interface problems with mixed boundary conditions

    Feng, Qiwei | Han, Bin | Minev, Peter

    Journal of Computational Physics, Vol. 497 (2024), Iss. P.112635

    https://doi.org/10.1016/j.jcp.2023.112635 [Citations: 0]