Nystrom methods for high-order CQ solutions of the wave equation in two dimensions

An investigation of high order Convolution Quadratures (CQ) methods for the solution of the wave equation in unbounded domains in two dimensions is presented. These rely on Nystrom discretizations for the solution of the ensemble of associated Laplace domain modified Helmholtz problems. Two classes of CQ discretizations are considered: one based on linear multistep methods and the other based on Runge-Kutta methods. Both are used in conjunction with Nystrom discretizations based on Alpert and QBX quadratures of Boundary Integral Equation (BIE) formulations of the Laplace domain Helmholtz problems with complex wavenumbers. CQ in conjunction with BIE is an excellent candidate to eventually explore numerical homogenization to replace a chaff cloud by a dispersive lossy dielectric that produces the same scattering. To this end, a variety of accuracy tests are presented that showcase the high-order in time convergence (up to and including fifth order) that the Nystrom CQ discretizations are capable of delivering for a variety of two dimensional single and multiple scatterers. Particular emphasis is given to Lipschitz boundaries and open arcs with both Dirichlet and Neumann boundary conditions

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

An efficient high-order algorithm for acoustic scattering from penetrable thin structures in three dimensions

This paper presents a high-order accelerated algorithm for the solution of the integral-equation formulation of volumetric scattering problems. The scheme is particularly well suited to the analysis of “thin” structures as they arise in certain applications (e.g., material coatings); in addition, it is also designed to be used in conjunction with existing low-order FFT-based codes to upgrade their order of accuracy through a suitable treatment of material interfaces. The high-order convergence of the new procedure is attained through a combination of changes of parametric variables (to resolve the singularities of the Green function) and “partitions of unity” (to allow for a simple implementation of spectrally accurate quadratures away from singular points). Accelerated evaluations of the interaction between degrees of freedom, on the other hand, are accomplished by incorporating (two-face) equivalent source approximations on Cartesian grids. A detailed account of the main algorithmic components of the scheme are presented, together with a brief review of the corresponding error and performance analyses which are exemplified with a variety of numerical results

Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D

We present an effective harmonic density interpolation method for the numerical evaluation of singular and nearly singular Laplace boundary integral operators and layer potentials in two and three spatial dimensions. The method relies on the use of Green's third identity and local Taylor-like interpolations of density functions in terms of harmonic polynomials. The proposed technique effectively regularizes the singularities present in boundary integral operators and layer potentials, and recasts the latter in terms of integrands that are bounded or even more regular, depending on the order of the density interpolation. The resulting boundary integrals can then be easily, accurately, and inexpensively evaluated by means of standard quadrature rules. A variety of numerical examples demonstrate the effectiveness of the technique when used in conjunction with the classical trapezoidal rule (to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type quadrature rule (to integrate over surfaces given as unions of non-overlapping quadrilateral patches) in three-dimensions