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

11/25/2021
by   Hadrien Montanelli, et al.
0

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.

READ FULL TEXT

Please sign up or login with your details

Forgot password? Click here to reset