Integral equation methods for acoustic scattering by fractals
Abstract
We study sound-soft time-harmonic acoustic scattering by general scatterers, including fractal scatterers, in 2D and 3D space. For an arbitrary compact scatterer $\Gamma$ we reformulate the Dirichlet boundary value problem for the Helmholtz equation as a first kind integral equation (IE) on $\Gamma$ involving the Newton potential. The IE is well-posed, except possibly at a countable set of frequencies, and reduces to existing single-layer boundary IEs when $\Gamma$ is the boundary of a bounded Lipschitz open set, a screen, or a multi-screen. When $\Gamma$ is uniformly of $d$-dimensional Hausdorff dimension in a sense we make precise (a $d$-set), the operator in our equation is an integral operator on $\Gamma$ with respect to $d$-dimensional Hausdorff measure, with kernel the Helmholtz fundamental solution, and we propose a piecewise-constant Galerkin discretization of the IE, which converges in the limit of vanishing mesh width. When $\Gamma$ is the fractal attractor of an iterated function system of contracting similarities we prove convergence rates under assumptions on $\Gamma$ and the IE solution, and describe a fully discrete implementation using recently proposed quadrature rules for singular integrals on fractals. We present numerical results for a range of examples and make our software available as a Julia code.
- Publication:
-
arXiv e-prints
- Pub Date:
- September 2023
- DOI:
- 10.48550/arXiv.2309.02184
- arXiv:
- arXiv:2309.02184
- Bibcode:
- 2023arXiv230902184C
- Keywords:
-
- Mathematics - Numerical Analysis;
- 28A80;
- 65R20;
- 78A45