Computing weakly singular and near-singular integrals over curved boundary elements
Numerical Analysis
2022-06-28 v3 Numerical Analysis
Abstract
We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with curved boundary elements. These are based on the computation of the preimage of the singularity in the reference element's space using Newton's method, singularity subtraction, the continuation approach, and transplanted Gauss quadrature. We demonstrate the accuracy of our method for quadratic basis functions and quadratic triangles with several numerical experiments, including the scattering by two half-spheres.
Cite
@article{arxiv.2111.13151,
title = {Computing weakly singular and near-singular integrals over curved boundary elements},
author = {Hadrien Montanelli and Matthieu Aussal and Houssem Haddar},
journal= {arXiv preprint arXiv:2111.13151},
year = {2022}
}