4 research outputs found
Recommended from our members
Diffraction on the Two-Dimensional Square Lattice
We solve the thin-slit diffraction problem for two-dimensional lattice waves. More precisely, for the discrete Helmholtz equation on the semi-infinite square lattice with data prescribed on the left boundary (the aperture), we use lattice Green's functions and a discrete Sommerfeld outgoing radiation condition to derive the exact solution everywhere in the lattice. The solution is a discrete convolution that can be evaluated in closed form for the wave number . For other wave numbers, we give a recursive algorithm for computing the convolution kernel
Recommended from our members
Spectral Optimization Problems Controlling Wave Phenomena
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 PDE-constrained 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 time-independent wave equation and the objective function governs some aspect of wave motion. We consider the shape optimization of functions of Dirichlet-Laplacian eigenvalues associated with the set of star-shaped, symmetric, bounded planar regions with smooth boundary. The boundary of such a region is represented using a Fourier-cosine series and the optimization problem is solved numerically using a quasi-Newton method. The method is applied to maximizing two particular nonsmooth functions of the eigenvalues: (a) the ratio of the n-th to first eigenvalues and (b) the ratio of the n-th eigenvalue gap to first eigenvalue. Both are generalizations of the Payne-Pólya-Weinberger 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 Dirichlet-Laplacian eigenvalues. For both spectral functions and each n, the optimal region has multiplicity two n-th 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 d-dimensional space and constant outside. The solution of the time-dependent wave equation with spatially-localized 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 non-selfadjoint eigenvalue problem consisting of the time-harmonic 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 well-known class of Bragg structures, where n(x) is constant on intervals whose length is one-quarter of the effective wavelength. Consider a system governed by the time-dependent Schroedinger equation in its ground state. When subjected to weak parametric forcing by an "ionizing field" (time-varying), the state decays with advancing time due to coupling of the bound state to radiation modes. The decay-rate 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 quasi-Newton 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 multi-scale 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 time-evolution. We consider a general class of two-dimensional 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 steady-state 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 inductor-capacitor circuit that can perform this transformation. Formulating this as an optimization problem, we numerically demonstrate its solvability using gradient-based 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 time-harmonic solutions to continuous and discrete wave equations
Hypergeometric expression for the resolvent of the discrete Laplacian in low dimensions
We present an explicit formula for the resolvent of the discrete Laplacian on
the square lattice, and compute its asymptotic expansions around thresholds in
low dimensions. As a by-product we obtain a closed formula for the fundamental
solution to the discrete Laplacian. For the proofs we express the resolvent in
a general dimension in terms of the Appell--Lauricella hypergeometric function
of type outside a disk encircling the spectrum. In low dimensions it
reduces to a generalized hypergeometric function, for which certain
transformation formulas are available for the desired expansions
A boundary algebraic formulation for plane strain elastodynamic scattering
Solving of elastodynamic problems arises in many scientific fields such as wave propagation in the ground, non-destructive testing, vibration design of buildings, or vibroacoustics in general. An integral formulation based on boundary algebraic equations is presented here. This formulation leads to a numerical method with a discretised boundary. An important advantage of the method over the standard boundary element method (BEM) is that no contour (2D) or surface (3D) integral needs to be computed. This feature is helpful in order to obtain a discrete version of the combined field integral equations (designed to damp numerically the fictitious eigenfrequencies) without difficulties caused by the evaluation of hypersingular integrals. The key aspects are: (i) the approach deals with discrete equations from the very beginning; (ii) discrete (instead of continuous) tensor Green's functions are considered (the methodology to evaluate them is demonstrated); (iii) the boundary must be described by means of a regular square grid. In order to overcome the drawback of this third condition the boundary integral is coupled, if needed, with a thin layer of finite elements. This improves the description of curved geometries and reduces numerical errors.
The properties of the method are demonstrated by means of numerical examples: the scattering of waves by objects and holes in an unbounded elastic medium, and an interior elastic problem.Peer ReviewedPostprint (author's final draft