A fast, higher-order solver for scattering by penetrable bodies in three dimensions
详细信息 查看全文
- 作者：Hyde ; E. McKay ; Bruno ; Oscar P.
- 关键词：35J05 ; 65D30 ; 65N12 ; 65R20 ; 65T50 ; 78A45 ; 78M25 ; Helmholtz equation ; Lippmann??????Schwinger equation ; Fourier series ; Fast Fourier transform ; High-order quadrature
- 刊名：Journal of Computational Physics
- 出版年：2005
- 出版时间：January 1, 2005
- 年：2005
- 卷：202
- 期：1
- 页码：236-261
- 全文大小：696 K
文摘
In this paper, we introduce a new fast, higher-order solver for scattering by inhomogeneous media in three dimensions. As in previously existing methods, the low complexity of our integral equation method, Formula Not Shown operations for an N point discretization, is obtained through extensive use of the fast Fourier transform (FFT) for the evaluation of convolutions. However, the present approach obtains significantly higher-order accuracy than these previous approaches, yielding, at worst, third-order far field accuracy (or substantially better for smooth scatterers), even for discontinuous and complex refractive index distributions (possibly containing severe geometric singularities such as corners and cusps). The increased order of convergence of our method results from (i) a partition of unity decomposition of the Greenaaas function into a smooth part with unbounded support and a singular part with compact support, and (ii) replacement of the (possibly discontinuous) scatterer by an appropriate aaaFourier smoothedaaa scatterer; the resulting convolutions can then be computed with higher-order accuracy by means of Formula Not Shown FFTs. We present a parallel implementation of our approach, and demonstrate the methodaaas efficiency and accuracy through a variety of computational examples. For a very large scatterer considered earlier in the literature (with a volume of 3648aa3, where aa is the wavelength), using the same number of points per wavelength and in computing times comparable to those required by the previous approach, the present algorithm produces far-field values whose errors are two orders of magnitude smaller than those reported previously.