An efficient high-order Nystr\"om scheme for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface

This text proposes a fast, rapidly convergent Nystr\"{o}m method for the solution of the Lippmann-Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by inhomogeneous obstacles, while allowing the material properties to jump across the interface. The method works with overlapping coordinate charts as a description of the given scatterer. In particular, it employs "partitions of unity" to simplify the implementation of high-order quadratures along with suitable changes of parametric variables to analytically resolve the singularities present in the integral operator to achieve desired accuracies in approximations. To deal with the discontinuous material interface in a high-order manner, a specialized quadrature is used in the boundary region. The approach further utilizes an FFT based strategy that uses equivalent source approximations to accelerate the evaluation of large number of interactions that arise in the approximation of the volumetric integral operator and thus achieves a reduced computational complexity of $O(N \log N)$ for an $N$-point discretization. A detailed discussion on the solution methodology along with a variety of numerical experiments to exemplify its performance in terms of both speed and accuracy are presented in this paper

Improved convergence of fast integral equation solvers for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface

In recent years, several fast solvers for the solution of the Lippmann-Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by penetrable inhomogeneous obstacles, have been proposed. While many of these fast methodologies exhibit rapid convergence for smoothly varying scattering configurations, the rate for most of them reduce to either linear or quadratic when material properties are allowed to jump across the interface. A notable exception to this is a recently introduced Nystr\"{o}m scheme [J. Comput. Phys., 311 (2016), 258--274] that utilizes a specialized quadrature in the boundary region for a high-order treatment of the material interface. In this text, we present a solution framework that relies on the specialized boundary integrator to enhance the convergence rate of other fast, low order methodologies without adding to their computational complexity of $O(N \log N)$ for an $N$-point discretization. In particular, to demonstrate the efficacy of the proposed framework, we explain its implementation to enhance the order to convergence of two schemes, one introduced by Duan and Rokhlin [J. Comput. Phys., 228(6) (2009), 2152--2174] that is based on a pre-corrected trapezoidal rule while the other by Bruno and Hyde [J. Comput. Phys., 200(2) (2004), 670--694] which relies on a suitable decomposition of the Green's function via Addition theorem. In addition to a detailed description of these methodologies, we also present a comparative performance study of the improved versions of these two and the Nystr\"{o}m solver in [J. Comput. Phys., 311 (2016), 258--274] through a wide range of numerical experiments

An efficient high-order Nyström scheme for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface

This text proposes a fast, rapidly convergent Nyström method for the solution of the Lippmann–Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by inhomogeneous obstacles, while allowing the material properties to jump across the interface. The method works with overlapping coordinate charts as a description of the given scatterer. In particular, it employs “partitions of unity” to simplify the implementation of high-order quadratures along with suitable changes of parametric variables to analytically resolve the singularities present in the integral operator to achieve desired accuracies in approximations. To deal with the discontinuous material interface in a high-order manner, a specialized quadrature is used in the boundary region. The approach further utilizes an FFT based strategy that uses equivalent source approximations to accelerate the evaluation of large number of interactions that arise in the approximation of the volumetric integral operator and thus achieves a reduced computational complexity of O(N log N) for an N-point discretization. A detailed discussion on the solution methodology along with a variety of numerical experiments to exemplify its performance are presented in this paper