A parallel fast algorithm for computing the Helmholtz integral operator in 3-D layered media

Authors:
Min Hyung Cho;Wei Cai
Affiliations:
Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223-0001, USA;Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223-0001, USA and Institute of Natural Sciences, Shanghai Jiaotong University, Shanghai 200240, ...
Venue:
Journal of Computational Physics
Year:
2012

Citing 8
Cited 0

Rapid solution of integral equations of scattering theory in two dimensions

Journal of Computational Physics
BI-CGSTAB: a fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems

SIAM Journal on Scientific and Statistical Computing
Adaptive multiresolution collocation methods for initial boundary value problems of nonlinear PDEs

SIAM Journal on Numerical Analysis
An adaptive spline wavelet ADI (SW-ADI) method for two-dimensional reaction-diffusion equations

Journal of Computational Physics
Fast calculations of dyadic Green's functions for electromagnetic scattering in a multilayered medium

Journal of Computational Physics
A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications

Journal of Computational Physics
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,

Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,
A wideband fast multipole method for the Helmholtz equation in three dimensions

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

In this paper, we propose a parallel fast algorithm for computing the product of the discretized Helmholtz integral operator in layered media and a vector in O(N"qN"z^2N"xN"ylogN"xN"y) operations. Here N"xN"yN"z is the number of sources and N"q is the number of quadrature points used in the evaluation of the Sommerfeld integral in the definition of layered media Green's function (for problems in thin-layer media, N"z=O(1)). Such a product forms the key step of many iterative solvers (such as the Krylov subspace based GMRES and BiCGSTAB) for linear systems arising from the integral equation methods for the Helmholtz equations. The fast solver is based on two important techniques which reduce the cost of quadrature summations in the Sommerfeld contour integral for Green's functions in 3-D layered media. The first technique is the removal of surface pole effects along the real axis integration contour by identifying the pole locations with a discrete wavelet transform; In the second technique, we apply a window-based high frequency filter to shorten the contour length. As a result, the integral operator for the 3-D layered media can be efficiently written as a sum of 2-D Hankel cylindrical integral operators, and the latter can be calculated by either a tree-code or a 2-D wideband fast multipole method in a fast manner. Numerical results show the efficiency and parallelism of the proposed fast algorithm.