High-Order Compact Finite Difference Methods for Solving the High-Dimensional Helmholtz Equations

In this paper, the sixth-order compact finite difference schemes for solving two-dimensional (2D) and three-dimensional (3D) Helmholtz equations are proposed. Firstly, the sixth-order compact difference operators for the second-order derivatives are applied to approximate the Laplace operator. Meanwhile, with the original differential equation, the sixth-order compact difference schemes are proposed. However, the truncation errors of the proposed scheme obviously depend on the unknowns, source function and wavenumber. Thus, we correct the truncation error of the sixth-order compact scheme to obtain an improved sixth-order compact scheme that is more accurate. Theoretically, the convergence and stability of the present improved method are proved. Finally, numerical tests verify that the improved schemes are more accurate.