Abstract:The recently developed quadrature by expansion (QBX) technique accurately evaluates the layer potentials with singular, weakly or nearly singular, or even hyper singular kernels in the integral equation reformulations of partial differential equations. The idea is to form a local complex polynomial or partial wave expansion centered at a point away from the boundary to avoid the singularity in the integrand, and then extrapolate the expansion at points near or even exactly on the boundary. In this paper, in addition to the local complex Taylor polynomial expansion, we derive new representations of the Laplace layer potentials using both the local complex polynomial and plane wave expansions. Unlike in the QBX, the local complex polynomial expansion in the new quadrature by two expansions (QB2X) method only collects the far-field contributions and its number of expansion terms can be analyzed using tools from the classical fast multipole method. The plane wave type expansion in the QB2X method better captures the layer potential features near the boundary. It is derived by applying the Fourier extension technique to the density and boundary geometry functions and then analytically utilizing the Residue Theorem for complex contour integrals. The internal connections of the layer potential with its density function and curvature on the boundary are explicitly revealed in the plane wave expansion and its error is bounded by the Fourier extension errors. We present preliminary numerical results to demonstrate the accuracy of the QB2X representations and to validate our analysis.

Testing of Quadrature Formulas for the Direct Value of the Normal Derivative of a Single-Layer Potential at the Boundary of a Thin Body

Improved Quadrature Formulas for the Direct Value of the Normal Derivative of a Single-Layer Potential

The Numerical Steepest Descent Path Method for Calculating Physical Optics Integrals on Smooth Conducting Quadratic Surfaces

Quadrature by Two Expansions: Evaluating Laplace Layer Potentials using Complex Polynomial and Plane Wave Expansions

Recursive reduction quadrature for the evaluation of Laplace layer potentials in three dimensions

On the accurate evaluation of unsteady Stokes layer potentials in moving two-dimensional geometries

Quadrature by Two Expansions for Evaluating Helmholtz Layer Potentials

Spectrally-accurate quadratures for evaluation of layer potentials close to the boundary for the 2D Stokes and Laplace equations

The adjoint double layer potential on smooth surfaces in and the Neumann problem

A local target specific quadrature by expansion method for evaluation of layer potentials in 3D

The adjoint double layer potential on smooth surfaces in $\mathbb{R}^3$ and the Neumann problem

A Fast Method for the Off-Boundary Evaluation of Laplace Layer Potentials by Convolution Sums

A High-Precision Calculation Method for Interface Normal and Curvature on an Unstructured Grid.

Solution of the Dirichlet problem by a finite difference analog of the boundary integral equation

A Weak Form Quadrature Element Formulation for Geometrically Exact Thin Shell Analysis

Formulas for Numerical Differentiation on a Uniform Mesh in the Presence of a Boundary Layer

Singularity swap quadrature for nearly singular line integrals on closed curves in two dimensions

Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation

The Helmholtz Dirichlet and Neumann problems on piecewise smooth open curves