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

by   Hadrien Montanelli, et al.

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.


Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme

We present a fast, high-order accurate and adaptive boundary integral sc...

Hybrid high-order method for singularly perturbed fourth-order problems on curved domains

We propose a novel hybrid high-order method (HHO) to approximate singula...

Efficient high-order singular quadrature schemes in magnetic fusion

Several problems in magnetically confined fusion, such as the computatio...

Efficient high-order accurate Fresnel diffraction via areal quadrature and the nonuniform FFT

We present a fast algorithm for computing the diffracted field from arbi...

Modified BDF2 schemes for subdiffusion models with a singular source term

The aim of this paper is to study the time stepping scheme for approxima...

Computing singular elements modulo squares

The group of singular elements was first introduced by Helmut Hasse and ...

Accurate Spectral Collocation Computation of High Order Eigenvalues for Singular Schrödinger Equations

We are concerned with the study of some classical spectral collocation m...