High-order close evaluation of Laplace layer potentials: A differential geometric approach

by   Hai Zhu, et al.

This paper presents a new approach for solving the close evaluation problem in three dimensions, commonly encountered while solving linear elliptic partial differential equations via potential theory. The goal is to evaluate layer potentials close to the boundary over which they are defined. The approach introduced here converts these nearly-singular integrals on a patch of the boundary to a set of non-singular line integrals on the patch boundary using the Stokes theorem on manifolds. A function approximation scheme based on harmonic polynomials is designed to express the integrand in a form that is suitable for applying the Stokes theorem. As long as the data – the boundary and the density function – is given in a high-order format, the double-layer potential and its derivatives can be evaluated with high-order accuracy using this scheme both on and off the boundary. In particular, we present numerical results demonstrating seventh-order convergence on a smooth, warped torus example achieving 10-digit accuracy in evaluating double layer potential at targets that are arbitrarily close to the boundary.


page 2

page 13

page 14

page 15


Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme

We present a fast, high-order accurate and adaptive boundary integral sc...

Layer potential identities and subtraction techniques

When using boundary integral equation methods, we represent solutions of...

Accurate quadrature of nearly singular line integrals in two and three dimensions by singularity swapping

The method of Helsing and co-workers evaluates Laplace and related layer...

Quadrature by fundamental solutions: kernel-independent layer potential evaluation for large collections of simple objects

Well-conditioned boundary integral methods for the solution of elliptic ...

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

The recently developed quadrature by expansion (QBX) technique accuratel...

An adaptive kernel-split quadrature method for parameter-dependent layer potentials

Panel-based, kernel-split quadrature is currently one of the most effici...

Fast multipole methods for evaluation of layer potentials with locally-corrected quadratures

While fast multipole methods (FMMs) are in widespread use for the rapid ...