Quadrature error estimates for layer potentials evaluated near curved surfaces in three dimensions

by   Ludvig af Klinteberg, et al.

The quadrature error associated with a regular quadrature rule for evaluation of a layer potential increases rapidly when the evaluation point approaches the surface and the integral becomes nearly singular. Error estimates are needed to determine when the accuracy is insufficient and a more costly special quadrature method should be utilized. The final result of this paper are such quadrature error estimates for the composite Gauss-Legendre rule and the global trapezoidal rule, when applied to evaluate layer potentials defined over smooth curved surfaces in R^3. The estimates have no unknown coefficients and can be efficiently evaluated given the discretization of the surface, invoking a local one-dimensional root-finding procedure. They are derived starting with integrals over curves, using complex analysis involving contour integrals, residue calculus and branch cuts. By complexifying the parameter plane, the theory can be used to derive estimates also for curves in in R^3. These results are then used in the derivation of the estimates for integrals over surfaces. In this procedure, we also obtain error estimates for layer potentials evaluated over curves in R^2. Such estimates combined with a local root-finding procedure for their evaluation were earlier derived for the composite Gauss-Legendre rule for layer potentials written on complex form [4]. This is here extended to provide quadrature error estimates for both complex and real formulations of layer potentials, both for the Gauss-Legendre and the trapezoidal rule. Numerical examples are given to illustrate the performance of the quadrature error estimates. The estimates for integration over curves are in many cases remarkably precise, and the estimates for curved surfaces in R^3 are also sufficiently precise, with sufficiently low computational cost, to be practically useful.



There are no comments yet.


page 10

page 19

page 25

page 27

page 31

page 33

page 34


Neglecting discretization corrections in regularized singular or nearly singular integrals

A method for computing singular or nearly singular integrals on closed s...

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 Two Expansions: Evaluating Laplace Layer Potentials using Complex Polynomial and Plane Wave Expansions

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

On the Approximation of Local Expansions of Laplace Potentials by the Fast Multipole Method

In this paper, we present a generalization of the classical error bounds...

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

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

Layer potential identities and subtraction techniques

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

Delineating Bone Surfaces in B-Mode Images Constrained by Physics of Ultrasound Propagation

Bone surface delineation in ultrasound is of interest due to its potenti...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.