Fast, high-order solution of surface scattering problems

We present a new algorithm for the numerical solution of problems of acoustic scattering by surfaces in three-dimensional space. This algorithm evaluates the scattered field, through fast, high-order solution of the boundary integral equation. The high-order of the solver is achieved through use of partition of unity together with analytical resolution of kernel singularities. The acceleration in turn, results from a novel approach which, based on high-order "two-face" equivalent source approximations, reduces the evaluation of far interactions to evaluation of 3-D FFTs. 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.