Fourier smoothed pre-corrected trapezoidal rule for solution of Lippmann-Schwinger integral equation

07/13/2020
by   Ambuj Pandey, et al.
0

For the numerical solution of the Lippmann-Schwinger equation, while the pre-corrected trapezoidal rule converges with high-order for smooth compactly supported densities, it exhibits only the linear convergence in the case of discontinuity in material properties across the interface. In this short article, we propose a Nyström solver based on "Fourier smoothed pre-corrected trapezoidal rule" that converges with second order for such scattering problems while maintaining the computational complexity of O(N log N). Moreover, the method is not only very simple to implement, it is also applicable to problems with geometrically complex inhomogeneities including those with corners and cusps. We present a variety of numerical experiments including comparative studies with competing approaches reported in [J. Comput. Phys., 200(2) (2004), 670–694] by Bruno and Hyde, and in [J. Fourier Anal. Appl., 11(4) (2005), 471-487 ] by Andersson and Holst to exemplify its performance in terms of speed and accuracy. This Fourier smoothed numerical integration scheme can also be adapted to other problems of interest where the convolution integral with discontinuous density is required to be computed.

READ FULL TEXT

page 7

page 8

research
07/08/2021

Fast accurate approximation of convolutions with weakly singular kernel and its applications

In this article, we present an O(N log N) rapidly convergent algorithm f...
research
01/16/2020

A high-order integral equation-based solver for the time-dependent Schrodinger equation

We introduce a numerical method for the solution of the time-dependent S...
research
10/07/2022

Corrected Trapezoidal Rule-IBIM for linearized Poisson-Boltzmann equation

In this paper, we solve the linearized Poisson-Boltzmann equation, used ...
research
04/16/2021

A short note on the accuracy of the discontinuous Galerkin method with reentrant faces

We study the convergence of the discontinuous Galerkin (DG) method appli...
research
12/18/2019

IsoGeometric Approximations for Electromagnetic Problems in Axisymmetric Domains

We propose a numerical method for the solution of electromagnetic proble...
research
07/27/2020

Zeta Correction: A New Approach to Constructing Corrected Trapezoidal Quadrature Rules for Singular Integral Operators

We introduce a new quadrature method for the discretization of boundary ...
research
10/12/2020

Efficient high-order accurate Fresnel diffraction via areal quadrature and the nonuniform FFT

We present a fast algorithm for computing the diffracted field from arbi...

Please sign up or login with your details

Forgot password? Click here to reset