skip to main content


Title: Fast numerical solutions to direct and inverse scattering for bi-anisotropic periodic Maxwell's equations
This paper is concerned with the numerical solution to the direct and inverse electromagnetic scattering problem for bi-anisotropic periodic structures. The direct problem can be reformulated as an integro-di erential equation. We study the existence and uniqueness of solution to the latter equation and analyze a spectral Galerkin method to solve it. This spectral method is based on a periodization technique which allows us to avoid the evaluation of the quasiperiodic Green's tensor and to use the fast Fourier transform in the numerical implementation of the method. For the inverse problem, we study the orthogonality sampling method to reconstruct the periodic structures from scattering data generated by only two incident fields. The sampling method is fast, simple to implement, regularization free, and very robust against noise in the data. Numerical examples for both direct and inverse problems are presented to examine the efficiency of the numerical solvers.  more » « less
Award ID(s):
1812693
NSF-PAR ID:
10381995
Author(s) / Creator(s):
;
Editor(s):
Nguyen, Dinh-Liem; Nguyen, Loc; Nguyen, Thi-Phong
Date Published:
Journal Name:
AMS Contemporary Mathematics
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract This paper is concerned with the inverse scattering problem for the three-dimensional Maxwell equations in bi-anisotropic periodic structures.The inverse scattering problem aims to determine the shape of bi-anisotropic periodic scatterers from electromagnetic near-field data at a fixed frequency.The factorization method is studied as an analytical and numerical tool for solving the inverse problem.We provide a rigorous justification of the factorization method which results in the unique determination and a fast imaging algorithm for the periodic scatterer.Numerical examples for imaging three-dimensional periodic structures are presented to examine the efficiency of the method. 
    more » « less
  2. Abstract This paper is concerned with the inverse problem of determining the shape of penetrable periodic scatterers from scattered field data. We propose a sampling method with a novel indicator function for solving this inverse problem. This indicator function is very simple to implement and robust against noise in the data. The resolution and stability analysis of the indicator function is analyzed. Our numerical study shows that the proposed sampling method is more stable than the factorization method and more efficient than the direct or orthogonality sampling method in reconstructing periodic scatterers. 
    more » « less
  3. We consider the inverse problem of determining the geometry of penetrable objects from scattering data generated by one incident wave at a fixed frequency. We first study an orthogonality sampling type method which is fast, simple to implement, and robust against noise in the data. This sampling method has a new imaging functional that is applicable to data measured in near field or far field regions. The resolution analysis of the imaging functional is analyzed where the explicit decay rate of the functional is established. A connection with the orthogonality sampling method by Potthast is also studied. The sampling method is then combined with a deep neural network to solve the inverse scattering problem. This combined method can be understood as a network using the image computed by the sampling method for the first layer and followed by the U-net architecture for the rest of the layers. The fast computation and the knowledge from the results of the sampling method help speed up the training of the network. The combination leads to a significant improvement in the reconstruction results initially obtained by the sampling method. The combined method is also able to invert some limited aperture experimental data without any additional transfer training. 
    more » « less
  4. null (Ed.)
    This paper is concerned with the inverse scattering problem which aims to determine the spatially distributed dielectric constant coefficient of the 2D Helmholtz equation from multifrequency backscatter data associated with a single direction of the incident plane wave. We propose a globally convergent convexification numerical algorithm to solve this nonlinear and ill-posed inverse problem. The key advantage of our method over conventional optimization approaches is that it does not require a good first guess about the solution. First, we eliminate the coefficient from the Helmholtz equation using a change of variables. Next, using a truncated expansion with respect to a special Fourier basis, we approximately reformulate the inverse problem as a system of quasilinear elliptic PDEs, which can be numerically solved by a weighted quasi-reversibility approach. The cost functional for the weighted quasi-reversibility method is constructed as a Tikhonov-like functional that involves a Carleman Weight Function. Our numerical study shows that, using a version of the gradient descent method, one can find the minimizer of this Tikhonov-like functional without any advanced a priori knowledge about it. 
    more » « less
  5. Abstract

    We formulate the inverse spectral theory of infinite gap Hill’s operators with bounded periodic potentials as a Riemann–Hilbert problem on a typically infinite collection of spectral bands and gaps. We establish a uniqueness theorem for this Riemann–Hilbert problem, which provides a new route to establishing unique determination of periodic potentials from spectral data. As the potentials evolve according to the Korteweg–de Vries Equation (KdV) equation, we use integrability to derive an associated Riemann–Hilbert problem with explicit time dependence. Basic principles from the theory of Riemann–Hilbert problems yield a new characterization of spectra for periodic potentials in terms of the existence of a solution to a scalar Riemann–Hilbert problem, and we derive a similar condition on the spectrum for the temporal periodicity for an evolution under the KdV equation.

     
    more » « less