AbstractWe study boundary element methods for time-harmonic scattering in $${\mathbb {R}}^n$$
R
n
($$n=2,3$$
n
=
2
,
3
) by a fractal planar screen, assumed to be a non-empty bounded subset $$\Gamma $$
Γ
of the hyperplane $$\Gamma _\infty ={\mathbb {R}}^{n-1}\times \{0\}$$
Γ
∞
=
R
n
-
1
×
{
0
}
. We consider two distinct cases: (i) $$\Gamma $$
Γ
is a relatively open subset of $$\Gamma _\infty $$
Γ
∞
with fractal boundary (e.g. the interior of the Koch snowflake in the case $$n=3$$
n
=
3
); (ii) $$\Gamma $$
Γ
is a compact fractal subset of $$\Gamma _\infty $$
Γ
∞
with empty interior (e.g. the Sierpinski triangle in the case $$n=3$$
n
=
3
). In both cases our numerical simulation strategy involves approximating the fractal screen $$\Gamma $$
Γ
by a sequence of smoother “prefractal” screens, for which we compute the scattered field using boundary element methods that discretise the associated first kind boundary integral equations. We prove sufficient conditions on the mesh sizes guaranteeing convergence to the limiting fractal solution, using the framework of Mosco convergence. We also provide numerical examples illustrating our theoretical results.
A fast multipole accelerated BEM for 3-D elastic wave computation
