We consider a system governed by the wave equation with index of refraction @math , taken to be variable within a bounded region @math and constant in @math . The solution of the time-dependent wave equation with initial data, which is localized in @math , spreads and decays with advancing time. This rate of decay can be measured (for @math , and more generally, @math odd) in terms of the eigenvalues of the scattering resonance problem, a non--self-adjoint eigenvalue problem governing the time-harmonic solutions of the wave (Helmholtz) equation which are outgoing at @math . Specifically, the rate of energy escape from @math is governed by the complex scattering eigenfrequency, which is closest to the real axis. We study the structural design problem: Find a refractive index profile @math within an admissible class which has a scattering frequency with minimal imaginary part. The admissible class is defined in terms of the compact suppo...
Results on the existence of optimal scattering resonances and general bounds on the imaginary parts of scattering resonances for Schr "odinger operators can be found in @cite_21 @cite_5 @cite_37 . Very recently, optimal designs have been considered in @cite_26 . Our results for the Helmholtz equation make use of some of the arguments introduced in these papers.
@cite_5 @cite_37 @cite_21 @cite_26
2164119610 2069009833 1517670975 2038870672
The Helmholtz equation governing wave propagation and scattering phenomena is difficult to solve numerically. Its discretization with piecewise linear finite elements results in typically large linear systems of equations. The inherently parallel domain decomposition methods constitute hence a promising class of preconditioners. An essential element of these methods is a good coarse space. Here, the Helmholtz equation presents a particular challenge, as even slight deviations from the optimal choice can be devastating. In this paper, we present a coarse space that is based on local eigenproblems involving the Dirichlet-to-Neumann operator. Our construction is completely automatic, ensuring good convergence rates without the need for parameter tuning. Moreover, it naturally respects local variations in the wave number and is hence suited also for heterogeneous Helmholtz problems. The resulting method is parallel by design and its efficiency is demonstrated on 2D homogeneous and heterogeneous numerical examples. We prove an asymptotic formula for the number of scattering resonances in a strip near the real axis when the trapped set is r-normally hyperbolic with r large and a pinching condition on the normal expansion rates holds. Our dynamical assumptions are stable under smooth perturbations and motivated by the setting of black holes. The key tool is a Fourier integral operator which microlocally projects onto the resonant states in the strip. In addition to Weyl law, this operator provides new information about microlocal concentration of resonant states. Preface to the Classics Edition Preface Symbols 1. The Riesz-Fredholm theory for compact operators 2. Regularity properties of surface potentials 3. Boundary-value problems for the scalar Helmholtz equation 4. Boundary-value problems for the time-harmonic Maxwell equations and the vector Helmholtz equation 5. Low frequency behavior of solutions to boundary-value problems in scattering theory 6. The inverse scattering problem: exact data 7. Improperly posed problems and compact families 8. The determination of the shape of an obstacle from inexact far-field data 9. Optimal control problems in radiation and scattering theory References Index. In this article we describe recent progress on the design, analysis and implementation of hybrid numerical-asymptotic boundary integral methods for boundary value problems for the Helmholtz equation that model time harmonic acoustic wave scattering in domains exterior to impenetrable obstacles. These hybrid methods combine conventional piecewise polynomial approximations with high-frequency asymptotics to build basis functions suitable for representing the oscillatory solutions. They have the potential to solve scattering problems accurately in a computation time that is (almost) independent of frequency and this has been realized for many model problems. The design and analysis of this class of methods requires new results on the analysis and numerical analysis of highly oscillatory boundary integral operators and on the high-frequency asymptotics of scattering problems. The implementation requires the development of appropriate quadrature rules for highly oscillatory integrals. This article contains a historical account of the development of this currently very active field, a detailed account of recent progress and, in addition, a number of original research results on the design, analysis and implementation of these methods.