A Hausdorff-measure boundary element method for acoustic scattering by fractal screens

A. M. Caetano,S. N. Chandler-Wilde,A. Gibbs,D. P. Hewett,A. Moiola

DOI: https://doi.org/10.1007/s00211-024-01399-7

IF: 2.5

2024-02-27

Numerische Mathematik

Abstract:Sound-soft fractal screens can scatter acoustic waves even when they have zero surface measure. To solve such scattering problems we make what appears to be the first application of the boundary element method (BEM) where each BEM basis function is supported in a fractal set, and the integration involved in the formation of the BEM matrix is with respect to a non-integer order Hausdorff measure rather than the usual (Lebesgue) surface measure. Using recent results on function spaces on fractals, we prove convergence of the Galerkin formulation of this "Hausdorff BEM" for acoustic scattering in ( ) when the scatterer, assumed to be a compact subset of , is a d -set for some , so that, in particular, the scatterer has Hausdorff dimension d . For a class of fractals that are attractors of iterated function systems, we prove convergence rates for the Hausdorff BEM and superconvergence for smooth antilinear functionals, under certain natural regularity assumptions on the solution of the underlying boundary integral equation. We also propose numerical quadrature routines for the implementation of our Hausdorff BEM, along with a fully discrete convergence analysis, via numerical (Hausdorff measure) integration estimates and inverse estimates on fractals, estimating the discrete condition numbers. Finally, we show numerical experiments that support the sharpness of our theoretical results, and our solution regularity assumptions, including results for scattering in by Cantor sets, and in by Cantor dusts.

mathematics, applied

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to effectively solve the scattering problem in acoustic wave scattering when the scatterer is a fractal screen with zero surface measure. Specifically, the paper focuses on the situation when the scatterer \(\Gamma\) itself is a fractal set, especially when \(\Gamma\) is the attractor of a self - similar iterative function system. In this case, the scatterer may not have a boundary in the traditional sense, so the traditional boundary element method (BEM) is no longer applicable.

To deal with such problems, the author proposes a new boundary element method - Hausdorff measure boundary element method (Hausdorff BEM). The core idea of this method is to use Hausdorff measure instead of the traditional Lebesgue surface measure for integration, and each BEM basis function is supported on a fractal set. In this way, the author is able to deal with fractal scatterers with non - integer dimensions and proves the convergence of this method in \(R^{n + 1}\) (\(n=1,2\)).

### Specific problem description

1. **Background**:
- Traditional acoustic wave scattering problems usually assume that the scatterer has a smooth boundary.
- When the scatterer is a fractal screen, the traditional BEM method is no longer applicable because the fractal screen may not have a boundary in the traditional sense.

2. **Objective**:
- Develop a new boundary element method (Hausdorff BEM) that can effectively solve the acoustic wave scattering problem in the case of a fractal screen.
- Prove the convergence of this method and provide numerical experiments to verify the theoretical results.

3. **Method**:
- Use Hausdorff measure instead of Lebesgue surface measure for integration.
- Each BEM basis function is supported on a fractal set.
- Implement discretization through the Galerkin method and prove its convergence.

4. **Main contributions**:
- Propose a new Hausdorff BEM method, which is suitable for the acoustic wave scattering problem of fractal screens.
- Prove the convergence of this method under specific conditions and give the convergence rate.
- Provide numerical experiments to verify the correctness of the theoretical results.

### Mathematical formulas

- **Hausdorff measure**:
\[
H^\alpha(E) := \lim_{\delta \to 0} \left( \inf \sum_{i = 1}^\infty (\text{diam}(U_i))^\alpha \right)
\]
where \(\text{diam}(U_i)\) represents the diameter of the set \(U_i\).

- **Hausdorff dimension**:
\[
\text{dim}_H(E) := \sup \{ \alpha \in \mathbb{R}^+ : H^\alpha(E) = \infty \} = \inf \{ \alpha \in \mathbb{R}^+ : H^\alpha(E) = 0 \}
\]

- **Definition of fractal set**:
\[
c_1 r^d \leq H^d(\Gamma \cap B_r(x)) \leq c_2 r^d, \quad x \in \Gamma, \quad 0 < r \leq 1
\]
where \(B_r(x)\) represents a closed ball centered at \(x\) with radius \(r\).

- **Integral inequality**:
\[
C_1 \int_0^{\text{diam}(\Gamma)} r^{d - 1} f(r) \, dr \leq \int_\Gamma f(|x - y|) \, dH^d(y) \leq C_2 \int_0^{\text{diam}(\Gamma)} r^{d - 1} f(r) \, dr
\]

- **Trace operator**:
\[