2,987 research outputs found

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

    Get PDF
    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(NlogN)O(N \log N) for an NN-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

    Fast and Accurate Computation of Time-Domain Acoustic Scattering Problems with Exact Nonreflecting Boundary Conditions

    Full text link
    This paper is concerned with fast and accurate computation of exterior wave equations truncated via exact circular or spherical nonreflecting boundary conditions (NRBCs, which are known to be nonlocal in both time and space). We first derive analytic expressions for the underlying convolution kernels, which allow for a rapid and accurate evaluation of the convolution with O(Nt)O(N_t) operations over NtN_t successive time steps. To handle the onlocality in space, we introduce the notion of boundary perturbation, which enables us to handle general bounded scatters by solving a sequence of wave equations in a regular domain. We propose an efficient spectral-Galerkin solver with Newmark's time integration for the truncated wave equation in the regular domain. We also provide ample numerical results to show high-order accuracy of NRBCs and efficiency of the proposed scheme.Comment: 22 pages with 9 figure
    corecore