A fast, high-order algorithm for the solution of surface scattering problems: Basic implementation, tests, and applications

Oscar P. Bruno

doi:10.1006/jcph.2001.6714

2001-05-01

231

PRE

Oscar P. Bruno

Kunyansky, LA

DOI：10.1006/jcph.2001.6714

原文链接

原文求助

摘要

En 中文

We present a new algorithm for the numerical solution of problems of acoustic scattering by surfaces in three-dimensional space. This algorithm evaluates scattered fields through fast, high-order solution of the corresponding boundary integral equation. The high-order accuracy of our solver is achieved through use of partitions of unity together with analytical resolution of kernel singularities. The acceleration, in turn, results from use of a novel approach which, based on high-order two-face equivalent source approximations, reduces the evaluation of far interactions to evaluation of 3-D fast Fourier transforms (FFTs), This approach is faster and substantially more accurate, and it runs on dramatically smaller memories than other FFT and k-space methods. The present algorithm computes one matrix-vector multiplication in O(N-6/5 log N) to O(N-4/3 log N) operations, where N is the number of surface discretization points. The latter estimate applies to smooth surfaces, for which our high-order algorithm provides accurate solutions with small values of N; the former, more favorable count is valid for highly complex surfaces requiring significant amounts of subwavelength sampling. Further, our approach exhibits super-algebraic convergence. it can be applied to smooth and nonsmooth scatterers, and it does not suffer from accuracy breakdowns of any kind. In this paper we introduce the main algorithmic components in our approach, and we demonstrate its performance with a variety of numerical results. In particular, we show that the present algorithm can evaluate accurately in a personal computer scattering from bodies of acoustical sizes Of several hundreds. (C) 2001 Academic Press.

Keyword:

FAST MULTIPOLE METHOD

ELECTROMAGNETIC SCATTERING

HELMHOLTZ-EQUATION

INTEGRAL-EQUATIONS

RAPID SOLUTION

LAPLACE