Spectral Optimization Problems Controlling Wave Phenomena

Title:

Spectral Optimization Problems Controlling Wave Phenomena

Author(s):

Osting, Braxton

Thesis Advisor(s):

Weinstein, Michael I.
Keyes, David E.

Date:

2011

Type:

Dissertations

Department:

Applied Physics and Applied Mathematics

Notes:

Ph.D., Columbia University.

Abstract:

Design problems seek a material arrangement or shape which fully harnesses the physical properties of the material(s) to create an environment in which a particular phenomena is most (or least) pronounced. Mathematically, design problems are formulated as PDEconstrained optimization problems to find the material arrangement that maximizes an objective function which expresses the desired behavior. The PDE constraint describes the relationship between the material and the phenomena of interest. The focus of this thesis is four design problems where the PDE constraint is a timeindependent wave equation and the objective function governs some aspect of wave motion. We consider the shape optimization of functions of DirichletLaplacian eigenvalues associated with the set of starshaped, symmetric, bounded planar regions with smooth boundary. The boundary of such a region is represented using a Fouriercosine series and the optimization problem is solved numerically using a quasiNewton method. The method is applied to maximizing two particular nonsmooth functions of the eigenvalues: (a) the ratio of the nth to first eigenvalues and (b) the ratio of the nth eigenvalue gap to first eigenvalue. Both are generalizations of the PaynePólyaWeinberger ratio. The optimal values of these ratios and regions for which they are attained, for n ≤ 13, are presented and interpreted as a study of the range of the DirichletLaplacian eigenvalues. For both spectral functions and each n, the optimal region has multiplicity two nth eigenvalue. We consider a system governed by the wave equation with index of refraction n(x), taken to be variable within a bounded region of ddimensional space and constant outside. The solution of the timedependent wave equation with spatiallylocalized initial data spreads and decays with advancing time. The rate of spatially localized energy decay can be measured in terms of the eigenvalues of the scattering resonance problem, a nonselfadjoint eigenvalue problem consisting of the timeharmonic wave (Helmholtz) equation with outgoing radiation condition at infinity. Specifically, the rate of energy escape 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 n* within an admissible class which has a scattering frequency with minimal imaginary part. The admissible class is defined in terms of the compact support of n(x)1 and pointwise upper and lower (material) bounds on n(x): 0 < n ≤ n(x) ≤ n+ < ∞. We formulate this problem as a constrained optimization problem and prove that an optimal structure, n* exists. Furthermore, n*(x) is piecewise constant and achieves the material bounds, i.e., n*(x) is n or n+ almost everywhere. In one dimension, we establish a connection between n*(x) and the wellknown class of Bragg structures, where n(x) is constant on intervals whose length is onequarter of the effective wavelength. Consider a system governed by the timedependent Schroedinger equation in its ground state. When subjected to weak parametric forcing by an "ionizing field" (timevarying), the state decays with advancing time due to coupling of the bound state to radiation modes. The decayrate of this metastable state is governed by Fermi's Golden Rule (FGR), which depends on the potential V and the details of the forcing. We pose the potential design problem: find V* which minimizes FGR (maximizes the lifetime of the state) over an admissible class of potentials with fixed spatial support. We formulate this problem as a constrained optimization problem and prove that an admissible optimal solution exists. Then, using quasiNewton methods, we compute locally optimal potentials. These have the structure of a truncated periodic potential with a localized defect. In contrast to optimal structures for other spectral optimization problems, the optimizing potentials appear to be interior points of the constraint set and to be smooth. The multiscale structures that emerge incorporate the physical mechanisms of energy confinement via material contrast and interference effects. An analysis of locally optimal potentials reveals local optimality is attained via two mechanisms: (i) decreasing the density of states near a resonant frequency in the continuum and (ii) tuning the oscillations of extended states to make FGR, an oscillatory integral, small. Finally, we explore the performance of optimal potentials via simulations of the timeevolution. We consider a general class of twodimensional passive propagation media, represented as a planar graph where nodes are capacitors connected to a common ground and edges are inductors. Capacitances and inductances are fixed in time but vary in space. Kirchhoff's laws give the time dynamics of voltage and current in the system. By harmonically forcing input nodes and collecting the resulting steadystate signal at output nodes, we obtain a linear, analog device that transforms the inputs to outputs. We pose the lattice synthesis problem: given a linear transformation, find the inductances and capacitances for an inductorcapacitor circuit that can perform this transformation. Formulating this as an optimization problem, we numerically demonstrate its solvability using gradientbased methods. By solving the lattice synthesis problem for various desired transformations, we design several devices that can be used for signal processing and filtering. In addition to these spectral optimization problems, we study several problems on wave propagation, diffraction, and scattering. The focus is on the behavior of timeharmonic solutions to continuous and discrete wave equations.

Subject(s):

Applied mathematics
Physics
Materials science
 Item views:

411
 Metadata:

View