Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method

被引：55

作者：

Capolino, F. ^{[1
]}

Wilton, D. R.

Johnson, W. A.

机构：

[1] Univ Siena, Dept Informat Engn, I-53100 Siena, Italy

[2] Univ Houston, Dept Elect & Comp Engn, Houston, TX 77004 USA

[3] Sandia Natl Labs, Electromagnet & Plasma Phys Anal Dept, Albuquerque, NM 87185 USA

来源：

JOURNAL OF COMPUTATIONAL PHYSICS | 2007年 / 223卷 / 01期

关键词：

arrays; series acceleration; fast methods; green function; gratings; numerical methods; periodic structures;

D O I：

10.1016/j.jcp.2006.09.013

中图分类号：

TP39 [计算机的应用];

学科分类号：

081203 ; 0835 ;

摘要：

The Ewald method is applied to accelerate the evaluation of the Green's function (GF) of an infinite equispaced linear array of point sources with linear phasing. Only a few terms are needed to evaluate Ewald sums, which are cast in terms of error functions and exponential integrals, to high accuracy. It is shown analytically that the choice of the standard "optimal" Ewald splitting parameter E-0 causes overflow errors at high frequencies (period large compared to the wavelength), and convergence rates are analyzed. A recipe for selecting the Ewald splitting parameter is provided. (c) 2006 Elsevier Inc. All rights reserved.

引用

页码：250 / 261

页数：12

共 23 条

[1]

Abramowitz M., 1970, HDB MATH FUNCTIONS

[2]

Bressan M, 2000, IEEE MTT-S, P1097, DOI 10.1109/MWSYM.2000.863548

[3] Efficient computation of the 2-D Green's function for 1-D periodic structures using the Ewald method [J].