Taylor expansion based fast multipole method for 3-D Helmholtz equations in layered media

JOURNAL OF COMPUTATIONAL PHYSICS

卷:401

Article

Wang, Bo^1,2  Chen, Duan³  Zhang, Bo⁴  Zhang, Wenzhong²  Cho, Min Hyung⁵  Cai, Wei² 

[1] Hunan Normal Univ, Sch Math & Stat, MOE LCSM, Changsha 410081, Hunan, Peoples R China

[2] Southern Methodist Univ, Dept Math, Dallas, TX 75275 USA

[3] Univ North Carolina Charlotte, Dept Math & Stat, Charlotte, NC 28223 USA

[4] Indiana Univ, Dept Comp Sci, Bloomington, IN 47408 USA

[5] Univ Massachusetts, Dept Math Sci, Lowell, MA 01854 USA

关键词: Fast multipole method;    Layered media;    Helmholtz equation;    Taylor expansion;   

DOI  :  10.1016/j.jcp.2019.109008

来源: Elsevier

【 摘 要 】

In this paper, we develop a Taylor expansion (TE) based fast multipole method (FMM) for low frequency 3D Helmholtz Green's function in layered media. Two forms of Taylor expansions, with either non-symmetric or symmetric derivatives of layered media Green's functions, are used for the implementations of the proposed TE-FMM. In the implementation with non-symmetric derivatives, an algorithm based on discrete complex image approximations and recurrence formulas is shown to be very efficient and accurate in computing the high order derivatives. Meanwhile, the implementation based on symmetric derivatives is more robust and pre-computed tables for the high order derivatives in translation operators are used. Numerical tests in layered media have validated the accuracy and O(N) complexity of the proposed algorithms. (C) 2019 Elsevier Inc. All rights reserved.

【 授权许可】

Free