Wave-number-explicit bounds in time-harmonic scattering

In this paper we consider the problem of scattering of time-harmonic acoustic waves by a bounded sound soft obstacle in two and three dimensions, studying dependence on the wave number in two classical formulations of this problem. The first is the standard variational/weak formulation in the part of the exterior domain contained in a large sphere, with an exact Dirichletto-Neumann map applied on the boundary. The second formulation is as a second kind boundary integral equation in which the solution is sought as a combined single- and double-layer potential. For the variational formulation we obtain, in the case when the obstacle is starlike, explicit upper and lower bounds which show that the inf-sup constant decreases like k −1 as the wave number k increases. We also give an example where the obstacle is not starlike and the inf-sup constant decreases at least as fast as k −2. For the boundary integral equation formulation, if the boundary is also Lipschitz and piecewise smooth, we show that the norm of the inverse boundary integral operator is bounded independently of k if the coupling parameter is chosen correctly. The methods we use also lead to explicit bounds on the solution of the scattering problem in the energy norm when the obstacle is starlike in which the dependence of the norm of the solution on the wave number and on the geometry are made explicit

On stability of discretizations of the Helmholtz equation (extended version)

We review the stability properties of several discretizations of the Helmholtz equation at large wavenumbers. For a model problem in a polygon, a complete $k$-explicit stability (including $k$-explicit stability of the continuous problem) and convergence theory for high order finite element methods is developed. In particular, quasi-optimality is shown for a fixed number of degrees of freedom per wavelength if the mesh size $h$ and the approximation order $p$ are selected such that $kh/p$ is sufficiently small and $p = O(\log k)$, and, additionally, appropriate mesh refinement is used near the vertices. We also review the stability properties of two classes of numerical schemes that use piecewise solutions of the homogeneous Helmholtz equation, namely, Least Squares methods and Discontinuous Galerkin (DG) methods. The latter includes the Ultra Weak Variational Formulation

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

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

On Open Scattering Channels for Manifolds with Ends

In the framework of time-dependent geometric scattering theory, we study the existence and completeness of the wave operators for perturbations of the Riemannian metric for the Laplacian on a complete manifold of dimension $n$. The smallness condition for the perturbation is expressed (intrinsically and coordinate free) in purely geometric terms using the harmonic radius; therefore, the size of the perturbation can be controlled in terms of local bounds on the injectivity radius and the Ricci-curvature. As an application of these ideas we obtain a stability result for the scattering matrix with respect to perturbations of the Riemannian metric. This stability result implies that a scattering channel which interacts with other channels preserves this property under small perturbations.Comment: updated version, now 43 page