Related papers: Adapting free-space fast multipole method for laye…
The fast multipole method (FMM) has had great success in reducing the computational complexity of solving the boundary integral form of the Helmholtz equation. We present a formulation of the Helmholtz FMM that uses Fourier basis functions…
This paper considers the approximation of spatial convolution with a given radial integral kernel. Previous studies have demonstrated that approximating spatial convolution using a system of partial differential equations (PDEs) can…
In a homogeneous medium, the far-field generated by a localized source can be expanded in terms of multipoles; the coefficients are determined by the moments of the localized charge distribution. We show that this structure survives to some…
A matrix basis formulation is introduced to represent the 3 x 3 dyadic Green's functions in the frequency domain for the Maxwell's equations and the elastic wave equation in layered media. The formulation can be used to decompose the…
Directional wave speeds variations in anisotropic elastic solids enables material characterisation capabilities, such as determination of elastic constants and volumetric measurement of crystallographic texture. However, achieving such…
Frequency domain Mie solutions to scattering from spheres have been used for a long time. However, deriving their transient analogue is a challenge as it involves an inverse Fourier transform of the spherical Hankel functions (and their…
In this paper we present a new highly efficient calculation method for the far field amplitude pattern that arises from scattering problems governed by the d-dimensional Helmholtz equation and, by extension, Schr\"odinger's equation. The…
The Marchenko method is a powerful tool for reconstructing full-wavefield Green's functions using surface-recorded seismic data. These Green's functions can then be utilized to produce subsurface images that are not affected by artifacts…
This paper introduces a fast algorithm, applicable throughout the electromagnetic spectrum, for the numerical solution of problems of scattering by periodic surfaces in two-dimensional space. The proposed algorithm remains highly accurate…
In this paper, we propose a fast multipole method (FMM) for 3-D linearized Poisson-Boltzmann (PB) equation in layered media. The main framework of the algorithm is analogous to the FMM for Helmholtz and Laplace equation in layered media…
A linear algebraic method named the shifted conjugate-orthogonal-conjugate-gradient method is introduced for large-scale electronic structure calculation. The method gives an iterative solver algorithm of the Green's function and the…
The desire for wide-field of view, large fractional bandwidth, high sensitivity, high spectral and temporal resolution has driven radio interferometry to the point of big data revolution where the data is represented in at least three…
The spectral decomposition of the elastic wave operator in a layered isotropic half-space is derived by means of standard functional analytic methods. Particular attention is paid to the coupled $P$-$SV$ waves. The problem is formulated…
This article presents a new high-order accurate algorithm for finding a particular solution to a linear, constant-coefficient partial differential equation (PDE) by means of a convolution of the volumetric source function with the Green's…
We describe an efficient algorithm for computing the matrix vector products that appear in the numerical resolution of boundary integral equations in 2 space dimension. This work is an extension of the so-called Sparse Cardinal Sine…
We present an accurate and efficient boundary integral (BI) method for simulating the deformation of drops and bubbles in Stokes flow with soluble surfactant. Soluble surfactant advects and diffuses in bulk fluids while adsorbing and…
We present a solver for plane wave scattering from a periodic dielectric grating with a large number $M$ of inclusions lying in each period of its middle layer.Such composite material geometries have a growing role in modern photonic…
The present study seeks to investigate mathematical structures of a multi-frequency subspace migration weighted by the natural logarithmic function for imaging of thin electromagnetic inhomogeneities from measured far-field pattern. To this…
We present a new boundary integral formulation for time-harmonic wave diffraction from two-dimensional structures with many layers of arbitrary periodic shape, such as multilayer dielectric gratings in TM polarization. Our scheme is robust…
In this paper, we consider band-structure calculations governed by the Helmholtz or Maxwell equations in piecewise homogeneous periodic materials. Methods based on boundary integral equations are natural in this context, since they…