Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D

Abstract

We present an effective harmonic density interpolation method for the numerical evaluation of singular and nearly singular Laplace boundary integral operators and layer potentials in two and three spatial dimensions. The method relies on the use of Green's third identity and local Taylor-like interpolations of density functions in terms of harmonic polynomials. The proposed technique effectively regularizes the singularities present in boundary integral operators and layer potentials, and recasts the latter in terms of integrands that are bounded or even more regular, depending on the order of the density interpolation. The resulting boundary integrals can then be easily, accurately, and inexpensively evaluated by means of standard quadrature rules. A variety of numerical examples demonstrate the effectiveness of the technique when used in conjunction with the classical trapezoidal rule (to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type quadrature rule (to integrate over surfaces given as unions of non-overlapping quadrilateral patches) in three-dimensions.

References in corpus (6)

• Non-singular boundary integral methods for fluid mechanics applications
• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• A robust and non-singular formulation of the boundary integral method for the potential problem
• Adaptive quadrature by expansion for layer potential evaluation in two dimensions
• A local target specific quadrature by expansion method for evaluation of layer potentials in 3D
• A Chebyshev-based rectangular-polar integral solver for scattering by general geometries described by non-overlapping patches

Cited by in corpus (8)

• A scalable computational platform for particulate Stokes suspensions
• Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme
• General-purpose kernel regularization of boundary integral equations via density interpolation
• Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries
• A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
• Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation
• A non-singular boundary element method for interactions between acoustical field sources and structures
• A complex-scaled boundary integral equation for time-harmonic water waves