Computing weakly singular and near-singular integrals in high-order boundary elements

11/25/2021

∙

by   Hadrien Montanelli, et al.

∙

0

∙

share

We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with high-order boundary elements. These are based on the computation of the preimage of the singularity on the reference element using Newton's method, singularity subtraction with high-order Taylor-like asymptotic expansions, the continuation approach, and transplanted Gauss quadrature. We demonstrate the accuracy of our method with several numerical experiments, including the scattering by two nearby half-spheres.