Dirac Assisted Tree Method for 1D Heterogeneous Helmholtz Equations with Arbitrary Variable Wave Numbers

DOI10.1016/J.CAMWA.2021.06.015MaRDI QIDQ6352853arXivWikidataFDO

Authors Bin Han, Michelle Michelle, Yau shu Wong

Publication date 2 November 2020

Published in arXiv

Full work available at URL https://doi.org/10.1016/J.CAMWA.2021.06.015

Mathematics Subject Classification ID

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Finite difference methods for boundary value problems involving PDEs (65N06) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Error bounds for boundary value problems involving PDEs (65N15) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)

arXiv classification

cs.NAmath.NA

Abstract: In this paper we introduce a new method called the Dirac Assisted Tree (DAT) method, which can handle 1D heterogeneous Helmholtz equations with arbitrarily large variable wave numbers. DAT breaks an original global problem into many parallel tree-structured small local problems, which are linked together to form a global solution by solving small linking problems. To solve the local problems in DAT, we propose a compact finite difference method (FDM) with arbitrarily high accuracy order and low numerical dispersion for piecewise smooth coefficients and variable wave numbers. This compact FDM is particularly appealing for DAT, because the local problems and their fluxes in DAT can be computed with high accuracy. DAT with such compact FDMs can solve heterogeneous Helmholtz equations with arbitrarily large variable wave numbers accurately by solving small linear systems - $4imes4$ matrices in the extreme case - with tridiagonal coefficient matrices in a parallel fashion. Several numerical examples are provided to illustrate the effectiveness of DAT using the $M$th order compact FDMs with $M=6,8$ for numerically solving heterogeneous Helmholtz equations with variable wave numbers. We shall also discuss how to solve some special 2D Helmholtz equations using DAT.

This page was built for publication: Dirac Assisted Tree Method for 1D Heterogeneous Helmholtz Equations with Arbitrary Variable Wave Numbers

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6352853)