Armin Lechleiter

Armin Lechleiter, Peter Monk

We consider time domain acoustic scattering from a penetrable medium with a variable sound speed. This problem can be reduced to solving a time domain volume Lippmann-Schwinger integral equation. Using convolution quadrature in time and trigonometric collocation in space we can compute an approximate solution. We prove that the time domain Lippmann-Schwinger equation has a unique solution and prove conditional convergence and error estimates for the fully discrete solution for smooth sound speeds. Preliminary numerical results show that the method behaves well even for discontinuous sound speeds.

Armin Lechleiter, Dinh Liem Nguyen

Transverse magnetic (TM) scattering of an electromagnetic wave from a periodic dielectric diffraction grating can mathematically be described by a volume integral equation. This volume integral equation, however, in general fails to feature a weakly singular integral operator. Nevertheless, after a suitable periodization, the involved integral operator can be efficiently evaluated on trigonometric polynomials using the fast Fourier transform (FFT) and iterative methods can be used to solve the integral equation. Using Fredholm theory, we prove that a trigonometric Galerkin discretization applied to the periodized integral equation converges with optimal order to the solution of the scattering problem. The main advantage of this FFT-based discretization scheme is that the resulting numerical method is particularly easy to implement, avoiding for instance the need to evaluate quasiperiodic Green's functions.

Armin Lechleiter, Ruming Zhang

Scattering problems for periodic structures have been studied a lot in the past few years. A main idea for numerical solution methods is to reduce such problems to one periodicity cell. In contrast to periodic settings, scattering from locally perturbed periodic surfaces is way more challenging. In this paper, we introduce and analyze a new numerical method to simulate scattering from locally perturbed periodic structures based on the Bloch transform. As this transform is applied only in periodic domains, we firstly rewrite the scattering problem artificially in a periodic domain. With the help of the Bloch transform, we secondly transform this problem into a coupled family of quasiperiodic problems posed in the periodicity cell. A numerical scheme then approximates the family of quasiperiodic solutions (we rely on the finite element method) and backtransformation provides the solution to the original scattering problem. In this paper, we give convergence analysis and error bounds for a Galerkin discretization in the spatial and the quasiperiodicity's unit cells. We also provide a simple and efficient way for implementation that does not require numerical integration in the quasiperiodicity, together with numerical examples for scattering from locally perturbed periodic surfaces computed by this scheme.

Armin Lechleiter, Ruming Zhang

Periodic surface structures are nowadays standard building blocks of optical devices. If such structures are illuminated by aperiodic time-harmonic incident waves as, e.g., Gaussian beams, the resulting surface scattering problem must be formulated in an unbounded layer including the periodic surface structure. An obvious recipe to avoid the need to discretize this problem in an unbounded domain is to set up an equivalent system of quasiperiodic scattering problems in a single (bounded) periodicity cell via the Floquet-Bloch transform. The solution to the original surface scattering problem then equals the inverse Floquet-Bloch transform applied to the family of solutions to the quasiperiodic problems, which simply requires to integrate these solutions in the quasiperiodicity parameter. A numerical scheme derived from this representation hence completely avoids the need to tackle differential equations on unbounded domains. In this paper, we provide rigorous convergence analysis and error bounds for such a scheme when applied to a two-dimensional model problem, relying upon a quadrature-based approximation to the inverse Floquet-Bloch transform and finite element approximations to quasiperiodic scattering problems. Our analysis essentially relies upon regularity results for the family of solutions to the quasiperiodic scattering problems in suitable mixed Sobolev spaces. We illustrate our error bounds as well as efficiency of the numerical scheme via several numerical examples.

Armin Lechleiter, Ruming Zhang

This paper concerns the inverse scattering problem to reconstruct a local perturbation in a periodic structure. Unlike the periodic problems, the periodicity for the scattered field no longer holds, thus classical methods, which reduce quasi-periodic fields in one periodic cell, are no longer available. Based on the Floquet-Bloch transform, a numerical method has been developed to solve the direct problem, that leads to a possibility to design an algorithm for the inverse problem. The numerical method introduced in this paper contains two steps. The first step is initialization, that is to locate the support of the perturbation by a simple method. This step reduces the inverse problem in an infinite domain into one periodic cell. The second step is to apply Newton-CG method to solve the associated optimization problem. The perturbation is then approximated by a finite spline basis. Numerical examples are given at the end of this paper, shows the efficiency of the numerical method.

Thies Gerken, Armin Lechleiter

We add a time-dependent potential to the inhomogeneous wave equation and consider the task of reconstructing this potential from measurements of the wave field. This dynamic inverse problem becomes more involved compared to static parameters, as, e.g. the dimensions of the parameter space do considerably increase. We give a specifically tailored existence and uniqueness result for the wave equation and compute the Fréchet derivative of the solution operator, for which also show the tangential cone condition. These results motivate the numerical reconstruction of the potential via successive linearization and regularized Newton-like methods. We present several numerical examples showing feasibility, reconstruction quality, and time efficiency of the resulting algorithm.