Estimation of quadrature errors for layer potentials evaluated near surfaces with spherical topology

Numerical simulations with rigid particles, drops, or vesicles constitute some examples that involve 3D objects with spherical topology. When the numerical method is based on boundary integral equations, the error in using a regular quadrature rule to approximate the layer potentials that appear in...

Full description

Saved in:
Bibliographic Details
Published in:Advances in computational mathematics 2023-12, Vol.49 (6), Article 87
Main Authors: Sorgentone, Chiara, Tornberg, Anna-Karin
Format: Article
Language:English
Subjects:
Boundary integral method
Close evaluation
Computational mathematics
Computational Mathematics and Numerical Analysis
Computational Science and Engineering
Discretization
Error estimate
Estimates
Gaussian grid
Integral equations
Layer potentials
Mathematical and Computational Biology
Mathematical Modeling and Industrial Mathematics
Mathematics
Mathematics and Statistics
Nearly singular
Numerical methods
Quadrature
Quadratures
Spherical topology
Topology
Visualization
Citations: Items that this one cites
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Numerical simulations with rigid particles, drops, or vesicles constitute some examples that involve 3D objects with spherical topology. When the numerical method is based on boundary integral equations, the error in using a regular quadrature rule to approximate the layer potentials that appear in the formulation will increase rapidly as the evaluation point approaches the surface and the integrand becomes sharply peaked. To determine when the accuracy becomes insufficient, and a more costly special quadrature method should be used, error estimates are needed. In this paper, we present quadrature error estimates for layer potentials evaluated near surfaces of genus 0, parametrized using a polar and an azimuthal angle, discretized by a combination of the Gauss-Legendre and the trapezoidal quadrature rules. The error estimates involve no unknown coefficients, but complex-valued roots of a specified distance function. The evaluation of the error estimates in general requires a one-dimensional local root-finding procedure, but for specific geometries, we obtain analytical results. Based on these explicit solutions, we derive simplified error estimates for layer potentials evaluated near spheres; these simple formulas depend only on the distance from the surface, the radius of the sphere, and the number of discretization points. The usefulness of these error estimates is illustrated with numerical examples.
ISSN:1019-7168
1572-9044
1572-9044
DOI:10.1007/s10444-023-10083-7