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

Abstract

Accurate evaluation of layer potentials is crucial when boundary integral equation methods are used to solve partial differential equations. Quadrature by expansion (QBX) is a recently introduced method that can offer high accuracy for singular and nearly singular integrals, using truncated expansions to locally represent the potential. The QBX method is typically based on a spherical harmonics expansion which when truncated at order $p$ has $O(p^2)$ terms. This expansion can equivalently be written with $p$ terms, however paying the price that the expansion coefficients will depend on the evaluation/target point. Based on this observation, we develop a target specific QBX method, and apply it to Laplace's equation on multiply connected domains. The method is local in that the QBX expansions only involve information from a neighborhood of the target point. An analysis of the truncation error in the QBX expansions is presented, practical parameter choices are discussed and the method is validated and tested on various problems.

References in corpus (4)

• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• Fast Ewald summation for free-space Stokes potentials
• Adaptive quadrature by expansion for layer potential evaluation in two dimensions
• Estimation of quadrature errors in layer potential evaluation using quadrature by expansion

Cited by in corpus (15)

• A scalable computational platform for particulate Stokes suspensions
• A highly accurate boundary integral equation method for surfactant-laden drops in 3D
• A Fast Algorithm for Quadrature by Expansion in Three Dimensions
• Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D
• A Fast Algorithm with Error Bounds for Quadrature by Expansion
• Regularized Single and Double Layer Integrals in 3D Stokes Flow
• Efficient high-order singular quadrature schemes in magnetic fusion
• Optimization of Fast Algorithms for Global Quadrature by Expansion Using Target-Specific Expansions
• A robust solver for elliptic PDEs in 3D complex geometries
• Fast multipole methods for evaluation of layer potentials with locally-corrected quadratures
• Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries
• Efficient Convergent Boundary Integral Methods for Slender Bodies
• A new mixed potential representation for the equations of unsteady, incompressible flow
• Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation
• On the Approximation of Local Expansions of Laplace Potentials by the Fast Multipole Method