An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched 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...

Full description

Saved in:
Bibliographic Details
Main Authors: Dongsheng Cheng, Jianjun Chen, Guangqing Long
Format: Article
Language:English
Published: Wiley 2021-01-01
Series:Discrete Dynamics in Nature and Society
Online Access:http://dx.doi.org/10.1155/2021/2539272
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1832551645789028352
author Dongsheng Cheng
Jianjun Chen
Guangqing Long
author_facet Dongsheng Cheng
Jianjun Chen
Guangqing Long
author_sort Dongsheng Cheng
collection DOAJ
description 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.
format Article
id doaj-art-0085c4d276bf40e2b475d05a9c1bf6b4
institution Kabale University
issn 1607-887X
language English
publishDate 2021-01-01
publisher Wiley
record_format Article
series Discrete Dynamics in Nature and Society
spelling doaj-art-0085c4d276bf40e2b475d05a9c1bf6b42025-02-03T06:00:49ZengWileyDiscrete Dynamics in Nature and Society1607-887X2021-01-01202110.1155/2021/2539272An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface BoundaryDongsheng Cheng0Jianjun Chen1Guangqing Long2School of Software EngineeringSchool of Computer Science and EngineeringGuangxi Key Lab of Human Machine Interaction and Intelligent DecisionIn 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.http://dx.doi.org/10.1155/2021/2539272
spellingShingle Dongsheng Cheng
Jianjun Chen
Guangqing Long
An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
Discrete Dynamics in Nature and Society
title An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
title_full An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
title_fullStr An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
title_full_unstemmed An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
title_short An Optimal Fourth-Order Finite Difference Scheme for the Helmholtz Equation Based on the Technique of Matched Interface Boundary
title_sort optimal fourth order finite difference scheme for the helmholtz equation based on the technique of matched interface boundary
url http://dx.doi.org/10.1155/2021/2539272
work_keys_str_mv AT dongshengcheng anoptimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary
AT jianjunchen anoptimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary
AT guangqinglong anoptimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary
AT dongshengcheng optimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary
AT jianjunchen optimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary
AT guangqinglong optimalfourthorderfinitedifferenceschemeforthehelmholtzequationbasedonthetechniqueofmatchedinterfaceboundary