High Order Compact Scheme and the Moving Mesh Method for Helmholtz Equation with Variable Wave Numbers

Abstract

High order compact scheme has extensively been used for solving the Helmholtz equation on uniform mesh. Although uniform spacing is computationally advantageous for isotropic problems, the non-uniform mesh will be superior to the problems with the anisotropic properties or boundary layers. In this paper, we proposed a compact scheme on nonuniform mesh which has at least third order accuracy for the 2D Helmholtz equation with variable wave number. A novel adaptive mesh algorithm is also presented to generate the orthogonal nonuniform mesh. The 2D Helmholtz equation is solved by combining the HOC scheme on nonuniform mesh with the moving mesh method. The multigrid method is used to solve the discrete algebraic system. Numerical experiments are executed to verify the accuracy and efficiency of the presented method.