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

arXiv2007.06293MaRDI QIDQ6344958FDOQ6344958

Authors: Ambuj Pandey, Akash Anand 

Publication date: 13 July 2020

Abstract: 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"{o}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(NlogN)$. 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.

This page was built for publication: Fourier smoothed pre-corrected trapezoidal rule for solution of Lippmann-Schwinger integral equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6344958)