On the implementation of perfectly matched layers in a 3D fourth‐order velocity‐stress finite‐difference scheme

2002 ◽  
Author(s):  
Carey Marcinkovich ◽  
Kim B. Olsen
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Fourth Order ◽  
Finite Difference Scheme ◽  
Perfectly Matched Layers ◽  
Matched Layers

scholarly journals An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary

2021 ◽  
Vol 2021 ◽  
pp. 1-16
Author(s):  
Dongsheng Cheng ◽  
Jianjun Chen ◽  
Guangqing Long
Keyword(s):  
Finite Difference ◽  
Boundary Conditions ◽  
Difference Scheme ◽  
Helmholtz Equation ◽  
Fourth Order ◽  
Finite Difference Scheme ◽  
Weighted Average ◽  
Numerical Dispersion ◽  
Taylor Formula ◽  
Interface Boundary

In this paper, a new optimal fourth-order 21-point finite difference scheme is proposed to solve the 2D Helmholtz equation numerically, with the technique of matched interface boundary (MIB) utilized to treat boundary problems. For the approximation of Laplacian, two sets of fourth-order difference schemes are derived firstly based on the Taylor formula, with a total of 21 grid points involved. Then, a weighted combination of the two schemes is employed in order to reduce the numerical dispersion, and the weights are determined by minimizing the dispersion. Similarly, for the discretization of the zeroth-order derivative term, a weighted average of all the 21 points is implemented to obtain the fourth-order accuracy. The new scheme is noncompact; hence, it encounters great difficulties in dealing with the boundary conditions, which is crucial to the order of convergence. To tackle this issue, the matched interface boundary (MIB) method is employed and developed, which is originally used to accommodate free edges in the discrete singular convolution analysis. Convergence analysis and dispersion analysis are performed. Numerical examples are given for various boundary conditions, which show that new scheme delivers a fourth order of accuracy and is efficient in reducing the numerical dispersion as well.

Sign in / Sign up

Export Citation Format

Share Document