Stephen Langdon

Timo Betcke, Simon N. Chandler-Wilde, Ivan G. Graham et al.

We consider the classical coupled, combined-field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle. In recent work, we have proved lower and upper bounds on the $L^2$ condition numbers for these formulations, and also on the norms of the classical acoustic single- and double-layer potential operators. These bounds to some extent make explicit the dependence of condition numbers on the wave number $k$, the geometry of the scatterer, and the coupling parameter. For example, with the usual choice of coupling parameter they show that, while the condition number grows like $k^{1/3}$ as $k\to\infty$, when the scatterer is a circle or sphere, it can grow as fast as $k^{7/5}$ for a class of `trapping' obstacles. In this paper we prove further bounds, sharpening and extending our previous results. In particular we show that there exist trapping obstacles for which the condition numbers grow as fast as $\exp(γk)$, for some $γ>0$, as $k\to\infty$ through some sequence. This result depends on exponential localisation bounds on Laplace eigenfunctions in an ellipse that we prove in the appendix. We also clarify the correct choice of coupling parameter in 2D for low $k$. In the second part of the paper we focus on the boundary element discretisation of these operators. We discuss the extent to which the bounds on the continuous operators are also satisfied by their discrete counterparts and, via numerical experiments, we provide supporting evidence for some of the theoretical results, both quantitative and asymptotic, indicating further which of the upper and lower bounds may be sharper.

David P. Hewett, Stephen Langdon, Simon N. Chandler-Wilde

We propose and analyse a hybrid numerical-asymptotic $hp$ boundary element method for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom $N$ increases, and that to achieve any desired accuracy it is sufficient to increase $N$ in proportion to the square of the logarithm of the frequency as the frequency increases (standard boundary element methods require $N$ to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary "breakwater" problem of propagation through an aperture in an infinite sound-hard screen.

Samuel P. Groth, David P. Hewett, Stephen Langdon

We present a novel hybrid numerical-asymptotic boundary element method for high frequency acoustic and electromagnetic scattering by penetrable (dielectric) convex polygons. Our method is based on a standard reformulation of the associated transmission boundary value problem as a direct boundary integral equation for the unknown Cauchy data, but with a nonstandard numerical discretization which efficiently captures the high frequency oscillatory behaviour. The Cauchy data is represented as a sum of the classical geometrical optics approximation, computed by a beam tracing algorithm, plus a contribution due to diffraction, computed by a Galerkin boundary element method using oscillatory basis functions chosen according to the principles of the Geometrical Theory of Diffraction. We demonstrate with a range of numerical experiments that our boundary element method can achieve a fixed accuracy of approximation using only a relatively small, frequency-independent number of degrees of freedom. Moreover, for the scattering scenarios we consider, the inclusion of the diffraction term provides an order of magnitude improvement in accuracy over the geometrical optics approximation alone.

Andrew Gibbs, Stephen Langdon, Andrea Moiola

For problems of time-harmonic scattering by polygonal obstacles, embedding formulae provide a useful means of computing the far-field coefficient induced by any incident plane wave, given the far-field coefficient of a relatively small set of canonical problems. The number of such problems to be solved depends only on the geometry of the scatterer. Whilst the formulae themselves are exact in theory, any implementation will inherit numerical error from the method used to solve the canonical problems. This error can lead to numerical instabilities. Here, we present an effective approach to identify and regulate these instabilities. This approach is subsequently extended to the case where the incident wave is a Herglotz wave function, and we suggest how this could potentially remove frequency dependence of a T-matrix method.

Simon N. Chandler-Wilde, Stephen Langdon

We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A.S. Fokas and co-authors, analysing particular instances of this method, proposed by J.A. DeSanto and co-authors, for problems of acoustic scattering by diffraction gratings.