Acoustic scattering : high frequency boundary element methods and unified transform methods

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. De-Santo and co-authors, for problems of acoustic scattering by diffraction gratings.

Testing a linear propagation module on some acoustic scattering problems

Abstract not available

Species identification of Chinese sturgeon using acoustic descriptors and ascertaining their spatial distribution in the spawning ground of Gezhouba Dam

This study has developed an improved subjective approach of classification in conjunction with Step wise DFA analysis to discriminate Chinese sturgeon signals from other targets. The results showed that all together 25 Chinese sturgeon echo-signals were detected in the spawning ground of Gezhouba Dam during the last 3 years, and the identification accuracy reached 90.9%. In Stepwise DFA, 24 out of 67 variables were applied in discrimination and identification. PCA combined with DFA was then used to ensure the significance of the 24 variables and detailed the identification pattern. The results indicated that we can discriminate Chinese sturgeon from other fish species and noise using certain descriptors such as the behaviour variables, echo characteristics and acoustic cross-section characteristics. However, identification of Chinese sturgeon from sediments is more difficult and needs a total of 24 variables. This is due to the limited knowledge about the acoustic-scattering properties of the substrate regions. Based on identified Chinese sturgeon individuals, 18 individuals were distributed in the region between the site of Gezhouba Dam and Miaozui reach, with a surface area of about 3.4 km(2). Seven individuals were distributed in the region between Miaozui and Yanshouba reach, with a surface area of about 13 km(2).

Scattering poles near the real axis for two strictly convex obstacles

Iantchenko, A., (2007) 'Scattering poles near the real axis for two strictly convex obstacles', Annales of the Institute Henri Poincar? 8 pp.513-568 RAE2008

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.

Boundary integral methods in high frequency scattering

In this article we review recent progress on the design, analysis and implementation of numerical-asymptotic boundary integral methods for the computation of frequency-domain acoustic scattering in a homogeneous unbounded medium by a bounded obstacle. The main aim of the methods is to allow computation of scattering at arbitrarily high frequency with finite computational resources.

The domain derivative in rough-surface scattering and rigorous estimates for first-order perturbation theory

We investigate Fréchet differentiability of the scattered field with respect to variation in the boundary in the case of time–harmonic acoustic scattering by an unbounded, sound–soft, one–dimensional rough surface. We rigorously prove the differentiability of the scattered field and derive a characterization of the Fréchet derivative as the solution to a Dirichlet boundary value problem. As an application of these results we give rigorous error estimates for first–order perturbation theory, justifying small perturbation methods that have a long history in the engineering literature. As an application of our rigorous estimates we show that a plane acoustic wave incident on a sound–soft rough surface can produce an unbounded scattered field.

On the stability and convergence of the finite section method for integral equation formulations of rough surface scattering

We consider the Dirichlet and Robin boundary value problems for the Helmholtz equation in a non-locally perturbed half-plane, modelling time harmonic acoustic scattering of an incident field by, respectively, sound-soft and impedance infinite rough surfaces.Recently proposed novel boundary integral equation formulations of these problems are discussed. It is usual in practical computations to truncate the infinite rough surface, solving a boundary integral equation on a finite section of the boundary, of length 2A, say. In the case of surfaces of small amplitude and slope we prove the stability and convergence as A→∞ of this approximation procedure. For surfaces of arbitrarily large amplitude and/or surface slope we prove stability and convergence of a modified finite section procedure in which the truncated boundary is ‘flattened’ in finite neighbourhoods of its two endpoints. Copyright © 2001 John Wiley & Sons, Ltd.

A Nyström method for a class of integral equations on the real line with applications to scattering by diffraction gratings and rough surfaces

We propose a Nystr¨om/product integration method for a class of second kind integral equations on the real line which arise in problems of two-dimensional scalar and elastic wave scattering by unbounded surfaces. Stability and convergence of the method is established with convergence rates dependent on the smoothness of components of the kernel. The method is applied to the problem of acoustic scattering by a sound soft one-dimensional surface which is the graph of a function f, and superalgebraic convergence is established in the case when f is infinitely smooth. Numerical results are presented illustrating this behavior for the case when f is periodic (the diffraction grating case). The Nystr¨om method for this problem is stable and convergent uniformly with respect to the period of the grating, in contrast to standard integral equation methods for diffraction gratings which fail at a countable set of grating periods.

Scattering by infinite one-dimensional rough surfaces

We consider the Dirichlet boundary-value problem for the Helmholtz equation in a non-locally perturbed half-plane. This problem models time-harmonic electromagnetic scattering by a one-dimensional, infinite, rough, perfectly conducting surface; the same problem arises in acoustic scattering by a sound-soft surface. ChandlerWilde & Zhang have suggested a radiation condition for this problem, a generalization of the Rayleigh expansion condition for diffraction gratings, and uniqueness of solution has been established. Recently, an integral equation formulation of the problem has also been proposed and, in the special case when the whole boundary is both Lyapunov and a small perturbation of a flat boundary, the unique solvability of this integral equation has been shown by Chandler-Wilde & Ross by operator perturbation arguments. In this paper we study the general case, with no limit on surface amplitudes or slopes, and show that the same integral equation has exactly one solution in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including the incident plane wave, the Dirichlet boundary-value problem for the scattered field has a unique solution.

A high frequency $hp$ boundary element method for scattering by convex polygons

In this paper we propose and analyze a hybrid $hp$ boundary element method for the solution of problems of high frequency acoustic scattering by sound-soft convex polygons, in which the approximation space is enriched with oscillatory basis functions which efficiently capture the high frequency asymptotics of the solution. We demonstrate, both theoretically and via numerical examples, exponential convergence with respect to the order of the polynomials, moreover providing rigorous error estimates for our approximations to the solution and to the far field pattern, in which the dependence on the frequency of all constants is explicit. Importantly, these estimates prove that, to achieve any desired accuracy in the computation of these quantities, it is sufficient to increase the number of degrees of freedom in proportion to the logarithm of the frequency as the frequency increases, in contrast to the at least linear growth required by conventional methods.

A high frequency boundary element method for scattering by a class of nonconvex obstacles

In this paper we propose and analyse a hybrid numerical-asymptotic boundary element method for the solution of problems of high frequency acoustic scattering by a class of sound-soft nonconvex polygons. The approximation space is enriched with carefully chosen oscillatory basis functions; these are selected via a study of the high frequency asymptotic behaviour of the solution. We demonstrate via a rigorous error analysis, supported by numerical examples, that to achieve any desired accuracy it is sufficient for the number of degrees of freedom to grow only in proportion to the logarithm of the frequency as the frequency increases, in contrast to the at least linear growth required by conventional methods. This appears to be the first such numerical analysis result for any problem of scattering by a nonconvex obstacle. Our analysis is based on new frequency-explicit bounds on the normal derivative of the solution on the boundary and on its analytic continuation into the complex plane.

The submarine volcano eruption off El Hierro Island: Effects on the scattering migrant biota and the evolution of the pelagic communities

[EN] The submarine volcano eruption off El Hierro Island (Canary Islands) on 10 October 2011 promoted dramatic perturbation of the water column leading to changes in the distribution of pelagic fauna. To study the response of the scattering biota, we combined acoustic data with hydrographic profiles and concurrent sea surface turbidity indexes from satellite imagery. We also monitored changes in the plankton and nekton communities through the eruptive and post-eruptive phases. Decrease of oxygen, acidification, rising temperature and deposition of chemicals in shallow waters resulted in a reduction of epipelagic stocks and a disruption of diel vertical migration (nocturnal ascent) of mesopelagic organisms. Furthermore, decreased light levels at depth caused by extinction in the volcanic plume resulted in a significant shallowing of the deep acoustic scattering layer. Once the eruption ceased, the distribution and abundances of the pelagic biota returned to baseline levels. There was no evidence of a volcano-induced bloom in the plankton community.

First-principles calculation of carrier-phonon scattering in n-type Si1−xGex alloys

First-principles electronic structure methods are used to find the rates of inelastic intravalley and intervalley n-type carrier scattering in Si1-xGex alloys. Scattering parameters for all relevant Delta and L intra- and intervalley scattering are calculated. The short-wavelength acoustic and the optical phonon modes in the alloy are computed using the random mass approximation, with interatomic forces calculated in the virtual crystal approximation using density functional perturbation theory. Optical phonon and intervalley scattering matrix elements are calculated from these modes of the disordered alloy. It is found that alloy disorder has only a small effect on the overall inelastic intervalley scattering rate at room temperature. Intravalley acoustic scattering rates are calculated within the deformation potential approximation. The acoustic deformation potentials are found directly and the range of validity of the deformation potential approximation verified in long-wavelength frozen phonon calculations. Details of the calculation of elastic alloy scattering rates presented in an earlier paper are also given. Elastic alloy disorder scattering is found to dominate over inelastic scattering, except for almost pure silicon (x approximate to 0) or almost pure germanium (x approximate to 1), where acoustic phonon scattering is predominant. The n-type carrier mobility, calculated from the total (elastic plus inelastic) scattering rate, using the Boltzmann transport equation in the relaxation time approximation, is in excellent agreement with experiments on bulk, unstrained alloys..