Numerical Quadrature for Singular Integrals on Fractals

We present and analyse numerical quadrature rules for evaluating smooth and singular integrals on self-similar fractal sets. The integration domain $\Gamma\subset\mathbb{R}^n$ is assumed to be the compact attractor of an iterated function system of contracting similarities satisfying the open set condition. Integration is with respect to the Hausdorff measure $\mathcal{H}^d$, where $d$ is the Hausdorff dimension of $\Gamma$, and both single and double integrals are considered. Our focus is on composite quadrature rules in which integrals over $\Gamma$ are decomposed into sums of integrals over suitable partitions of $\Gamma$ into self-similar sub-components. For certain singular integrands of logarithmic or algebraic type we show how in the context of such a partitioning the self-similarity of $\Gamma$ can be exploited to express the singular integral exactly in terms of regular integrals. For the evaluation of these regular integrals we adopt a composite barycentre rule, which for sufficiently regular integrands exhibits second-order convergence with respect to the maximum diameter of the sub-components. As an application we show how this approach, combined with a singularity-subtraction technique, can be used to accurately evaluate the singular double integrals that arise in Hausdorff-measure Galerkin boundary element methods for acoustic wave scattering by fractal screens.