An official website of the United States government
This paper addresses the challenging and interesting inverse problem of reconstructing the spatially varying dielectric constant of a medium from phaseless backscattering measurements generated by single-point illumination. The underlying mathematical model is governed by the three-dimensional Helmholtz equation, and the available data consist solely of the magnitude of the scattered wave field. To address the nonlinearity and servere ill-posedness of this phaseless inverse scattering problem, we introduce a robust, globally convergent numerical framework combining several key regularization strategies. Our method first employs a phase retrieval step based on the Wentzel--Kramers--Brillouin (WKB) ansatz, where the lost phase information is reconstructed by solving a nonlinear optimization problem. Subsequently, we implement a Fourier-based dimension reduction technique, transforming the original problem into a more stable system of elliptic equations with Cauchy boundary conditions. To solve this resulting system reliably, we apply the Carleman convexification approach, constructing a strictly convex weighted cost functional whose global minimizer provides an accurate approximation of the true solution. Numerical simulations using synthetic data with high noise levels demonstrate the effectiveness and robustness of the proposed method, confirming its capability to accurately recover both the geometric location and contrast of hidden scatterers.
more »
« less
Solving the inverse scattering problem via Carleman-based contraction mapping
This paper addresses the inverse scattering problem in a domain $$\Omega$$. The input data, measured outside $$\Omega$$, involve the waves generated by the interaction of plane waves with various directions and unknown scatterers that are fully occluded inside $$\Omega$$. The output of this problem is the spatial dielectric constant of these scatterers. Our approach to solving this problem consists of two primary stages. Initially, we eliminate the unknown dielectric constant from the governing equation, resulting in a system of partial differential equations. Subsequently, we develop the Carleman contraction mapping method to effectively tackle this system. It is noteworthy to highlight the robustness of this method. It does not require a precise initial guess of the true solution, and its computational cost is relatively inexpensive. Some numerical examples are presented.
more »
« less
- Award ID(s):
- 2208159
- PAR ID:
- 10612298
- Publisher / Repository:
- preprint arXiv:2404.04145
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Motivated by the physics literature on “photonic doping” of scatterers made from “epsilon‐near‐zero” (ENZ) materials, we consider how the scattering of time‐harmonic TM electromagnetic waves by a cylindrical ENZ region is affected by the presence of a “dopant” in which the dielectric permittivity is not near zero. Mathematically, this reduces to analysis of a 2D Helmholtz equation with a piecewise‐constant, complex valued coefficientathat is nearly infinite (say with ) in . We show (under suitable hypotheses) that the solutionudepends analytically on δ near 0, and we give a simple PDE characterization of the terms in its Taylor expansion. For the application to photonic doping, it is the leading‐order corrections in δ that are most interesting: they explain why photonic doping is only mildly affected by the presence of losses, and why it is seen even at frequencies where the dielectric permittivity is merely small. Equally important: our results include a PDE characterization of the leading‐order electric field in the ENZ region as , whereas the existing literature on photonic doping provides only the leading‐order magnetic field.more » « less
-
Controlling the propagation of optical fields in three dimensions using arrays of discrete dielectric scatterers is an active area of research. These arrays can create optical elements with functionalities unrealizable in conventional optics. Here, we present an inverse design method based on the inverse Mie scattering problem for producing three-dimensional optical field patterns. Using this method, we demonstrate a device that focuses 1.55-μm light into a depth-variant discrete helical pattern. The reported device is fabricated using two-photon lithography and has a footprint of 144 μm by 144 μm, the largest of any inverse-designed photonic structure to date. This inverse design method constitutes an important step toward designer free-space optics, where unique optical elements are produced for user-specified functionalities.more » « less
-
Inertia-gravity waves in the atmosphere and ocean are transported and refracted by geostrophic turbulent currents. Provided that the wave group velocity is much greater than the speed of geostrophic turbulent currents, kinetic theory can be used to obtain a comprehensive statistical description of the resulting interaction (Savva et al. , J. Fluid Mech. , vol. 916, 2021, A6). The leading-order process is scattering of wave energy along a surface of constant frequency, $$\omega$$ , in wavenumber space. The constant- $$\omega$$ surface corresponding to the linear dispersion relation of inertia-gravity waves is a cone extending to arbitrarily high wavenumbers. Thus, wave scattering by geostrophic turbulence results in a cascade of wave energy to high wavenumbers on the surface of the constant- $$\omega$$ cone. Solution of the kinetic equations shows establishment of a wave kinetic energy spectrum $$\sim k_h^{-2}$$ , where $$k_h$$ is the horizontal wavenumber.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Posted Content:
- The DOI is not currently available.
