Michael V. Klibanov

Michael V. Klibanov, Aleksandr E. Kolesov

A version of the so-called "convexification" numerical method for a coefficient inverse scattering problem for the 3D Hemholtz equation is developed analytically and tested numerically. Backscattering data are used, which result from a single direction of the propagation of the incident plane wave on an interval of frequencies. The method converges globally. The idea is to construct a weighted Tikhonov-like functional. The key element of this functional is the presence of the so-called Carleman Weight Function (CWF). This is the function which is involved in the Carleman estimate for the Laplace operator. This functional is strictly convex on any appropriate ball in a Hilbert space for an appropriate choice of the parameters of the CWF. Thus, both the absence of local minima and convergence of minimizers to the exact solution are guaranteed. Numerical tests demonstrate a good performance of the resulting algorithm. Unlikeprevious the so-called tail functions globally convergent method, we neither do not impose the smallness assumption of the interval of wavenumbers, nor we do not iterate with respect to the so-called tail functions.

Michael V. Klibanov, Aleksandr E. Kolesov, Anders Sullivan et al.

A new version of the convexification method is developed analytically and tested numerically for a 1-D coefficient inverse problem in the frequency domain. Unlike the previous version, this one does not use the so-called "tail function", which is a complement of a certain truncated integral with respect to the wave number. Globally strictly convex cost functional is constructed with the Carleman Weight Function. Global convergence of the gradient projection method to the correct solution is proved. Numerical tests are conducted for both computationally simulated and experimental data.

Michael V. Klibanov, Aleksandr E. Kolesov, Lam Nguyen et al.

A new numerical method is proposed for a 1-D inverse medium scattering problem with multi-frequency data. This method is based on the construction of a weighted cost functional. The weight is a Carleman Weight Function (CWF). In other words, this is the function, which is present in the Carleman estimate for the undelying differential operator. The presence of the CWF makes this functional strictly convex on any a priori chosen ball with the center at $\left\{ 0\right\} $ in an appropriate Hilbert space. Convergence of the gradient minimization method to the exact solution starting from any point of that ball is proven. Computational results for both computationally simulated and experimental data show a good accuracy of this method.

Larisa Beilina, Nguyen Trung Thành, Michael V. Klibanov et al.

We consider a two-stage numerical procedure for imaging of objects buried in dry sand using time-dependent backscattering experimental radar measurements. These measurements are generated by a single point source of electric pulses and are collected using a microwave scattering facility which was built at the University of North Carolina at Charlotte. Our imaging problem is formulated as the inverse problem of the reconstruction of the spatially distributed dielectric permittivity $\varepsilon_\mathrm{r}\left(\mathbf{x}\right), \ \mathbf{x}\in \mathbb{R}^{3}$, which is an unknown coefficient in Maxwell's equations. On the first stage an approximately globally convergent method is applied to get a good first approximation for the exact solution. On the second stage a local adaptive finite element method is applied to refine the solution obtained on the first stage. The two-stage numerical procedure results in accurate imaging of all three components of interest of targets: shapes, locations and refractive indices. In this paper we briefly describe methods and present new reconstruction results for both stages.

Dinh-Liem Nguyen, Michael V. Klibanov, Loc H. Nguyen et al.

This paper is concerned with the numerical solution to a 3D coefficient inverse problem for buried objects with multi-frequency experimental data. The measured data, which are associated with a single direction of an incident plane wave, are backscatter data for targets buried in a sandbox. These raw scattering data were collected using a microwave scattering facility at the University of North Carolina at Charlotte. We develop a data preprocessing procedure and exploit a newly developed globally convergent inversion method for solving the inverse problem with these preprocessed data. It is shown that dielectric constants of the buried targets as well as their locations are reconstructed with a very good accuracy. We also prove a new analytical result which rigorously justifies an important step of the so-called "data propagation" procedure.

Dinh-Liem Nguyen, Michael V. Klibanov, Loc H. Nguyen et al.

We analyze in this paper the performance of a newly developed globally convergent numerical method for a coefficient inverse problem for the case of multi-frequency experimental backscatter data associated to a single incident wave. These data were collected using a microwave scattering facility at the University of North Carolina at Charlotte. The challenges for the inverse problem under the consideration are not only from its high nonlinearity and severe ill-posedness but also from the facts that the amount of the measured data is minimal and that these raw data are contaminated by a significant amount of noise, due to a non-ideal experimental setup. This setup is motivated by our target application in detecting and identifying explosives. We show in this paper how the raw data can be preprocessed and successfully inverted using our inversion method. More precisely, we are able to reconstruct the dielectric constants and the locations of the scattering objects with a good accuracy, without using any advanced \emph{a priori} knowledge of their physical and geometrical properties.

Michael V. Klibanov, Jingzhi Li, Wenlong Zhang

We propose a new numerical method to reconstruct the isotropic electrical conductivity from measured restricted Dirichlet-to-Neumann map data in electrical impedance tomography (EIT) model. "Restricted Dirichlet-to-Neumann (DtN) map data" means that the Dirichlet and Neumann boundary data for EIT are generated by a point source running either along an interval of a straight line or along a curve located outside of the domain of interest. We "convexify" the problem via constructing a globally strictly convex Tikhonov-like functional using a Carleman Weight Function. In particular, two new Carleman estimates are established. Global convergenceto the correct solution of the gradient projection method for this functional is proven. Numerical examples demonstrate a good performance of this numerical procedure.

Michael V. Klibanov, Dinh-Liem Nguyen, Loc H. Nguyen

This paper is concerned with a numerical method for a 3D coefficient inverse problem with phaseless scattering data. These are multi-frequency data generated by a single direction of the incident plane wave. Our numerical procedure consists of two stages. The first stage aims to reconstruct the (approximate) scattered field at the plane of measurements from its intensity. We present an algorithm for the reconstruction process and prove a uniqueness result of this reconstruction. After obtaining the approximate scattered field, we exploit a newly developed globally convergent numerical method to solve the coefficient inverse problem with the phased scattering data. The latter is the second stage of our algorithm. Numerical examples are presented to demonstrate the performance of our method. Finally, we present a numerical study which aims to show that, under a certain assumption, the solution of the scattering problem for the 3D scalar Helmholtz equation can be used to approximate the component of the electric field which was originally incident upon the medium.

Michael V. Klibanov, Anatoly G. Yagola

The key tool of this paper is a new Carleman estimate for an arbitrary parabolic operator of the second order for the case of reversed time data. This estimate works on an arbitrary time interval. On the other hand, the previously known Carleman estimate for the reversed time case works only on a sufficiently small time interval. First, a stability estimate is proven. Next, the quasi-reversibility numerical method is proposed for an arbitrary time interval for the linear case. This is unlike a sufficiently small time interval in the previous work. The convergence rate for the quasi-reversibility method is established. Finally, the quasilinear parabolic equation with reversed time is considered. A weighted globally strictly convex Tikhonov-like functional is constructed. The weight is the Carleman Weight Function which is involved in that Carleman estimate. The global convergence of the gradient projection method to the exact solution is proved for this functional.

Michael V. Klibanov, Dinh-Liem Nguyen, Loc H. Nguyen et al.

The goal of this paper is to reconstruct spatially distributed dielectric constants from complex-valued scattered wave field by solving a 3D coefficient inverse problem for the Helmholtz equation at multi-frequencies. The data are generated by only a single direction of the incident plane wave. To solve this inverse problem, a globally convergent algorithm is analytically developed. We prove that this algorithm provides a good approximation for the exact coefficient without any \textit{a priori} knowledge of any point in a small neighborhood of that coefficient. This is the main advantage of our method, compared with classical approaches using optimization schemes. Numerical results are presented for both computationally simulated data and experimental data. Potential applications of this problem are in detection and identification of explosive-like targets.

Aleksandr E. Kolesov, Michael V. Klibanov, Loc H. Nguyen et al.

The recently developed globally convergent numerical method for an inverse medium problem for the Helmholtz equation is tested on experimental data. The data were originally collected in the time domain, whereas the method works in the frequency domain with the multi-frequency data. Due to a huge discrepancy between the collected and computationally simulated data, the straightforward Fourier transform of the experimental data does not work. Hence, it is necessary to develop a heuristic data preprocessing procedure. This procedure is described. The preprocessed data are used as the input for the inversion algorithm. Numerical results demonstrate good accuracy in the reconstruction of both refracive indices and locations of targets. Furthermore, the reconstruction errors for refractive indices of dielectric targets are significantly less than errors of a posteriori direct measurements.

Michael V. Klibanov

Numerical issues for the 3D travel time tomography problem with non-overdetemined data are considered. Truncated Fourier series with respect to a special orthonormal basis of functions depending on the source position is used. In addition, truncated trigonometric Fourier series with respect to two out of three spatial variables is used. First, the Lipschitz stability estimate is obtained. Next, a globally convergent numerical method is constructed using a Carleman estimate for an integral operator.

Shi Chen, Michael V. Klibanov, Kevin McGoff et al.

We apply a convexification-based numerical method to forecast public sentiment dynamics using Mean Field Games (MFGs). The theoretical foundation for the convexification approach, established in our prior work, guarantees global convergence to the unique solution to the MFG system. The present work demonstrates the practical potential of this framework using real-world sentiment data extracted from social media public discussion during the COVID-19 pandemic. The results show that the MFG model with appropriate parameters and convexification yields sentiment density predictions that align closely with observed data and satisfy the governing equations. While current parameter selection relies on manual calibration, our findings establish the first proof-of-concept evidence that MFG models can capture complex temporal patterns in public sentiment, laying the groundwork for future work on systematic parameter identification methods, i.e. solutions of coefficient inverse problems for the MFG system.

Michael V. Klibanov, Aleksandr E. Kolesov, Dinh-Liem Nguyen

We present in this paper a novel numerical reconstruction method for solving a 3D coefficient inverse problem with scattering data generated by a single direction of the incident plane wave. This inverse problem is well-known to be a highly nonlinear and ill-posed problem. Therefore, optimization-based reconstruction methods for solving this problem would typically suffer from the local-minima trapping and require strong a priori information of the solution. To avoid these problems, in our numerical method, we aim to construct a cost functional with a globally strictly convex property, whose minimizer can provide a good approximation for the exact solution of the inverse problem. The key ingredients for the construction of such functional are an integro-differential formulation of the inverse problem and a Carleman weight function. Under a (partial) finite difference approximation, the global strict convexity is proven using the tool of Carleman estimates. The global convergence of the gradient projection method to the exact solution is proven as well. We demonstrate the efficiency of our reconstruction method via a numerical study of experimental backscatter data for buried objects.

Michael V. Klibanov, Hui Liu, Loc H. Nguyen

A 3-D inverse medium problem in the frequency domain is considered. Another name for this problem is Coefficient Inverse Problem. The goal is to reconstruct spatially distributed dielectric constants from scattering data. Potential applications are in detection and identification of explosive-like targets. A single incident plane wave and multiple frequencies are used. A new numerical method is proposed. A theorem is proved, which claims that a small neigborhood of the exact solution of that problem is reached by this method without any advanced knowledge of that neighborhood. We call this property of that numerical method "global convergence". Results of numerical experiments for the case of the backscattering data are presented.

Michael V. Klibanov, Jingzhi Li, Tian Niu et al.

The first globally convergent numerical method is developed for a coefficient inverse problem (CIP) for the $n-$d, $n\geq 2$ wave equation with the unknown potential in the most challenging case when the $δ-$ function is present in the initial condition with a single location of the point source. In fact, an approximate mathematical model for that CIP is derived. That globally convergent numerical method is developed for this model. This is a new version of the so-called convexification numerical method. Uniqueness theorem is proven as well within the framework of that approximate mathematical model. The question about uniqueness of this CIP was first posed by a famous mathematician I. M. Gelfand in 1954 as an $n-$d ($n=2,3$) extension of the fundamental theorem of V.A. Marchenko in the 1-d case (1950). Based on a Carleman estimate, convergence analysis is carried out. This analysis ensures the global convergence of the proposed numerical method, i.e. it is not necessary to have a good first guess for the solution. Exhaustive computational experiments with noisy data demonstrate a high reconstruction accuracy of complicated structures. In particular, this accuracy points towards a high adequacy of that approximate mathematical model.

Larisa Beilina, Nguyen Trung Thành, Michael V. Klibanov et al.

We perform extended studies of an adaptive finite element method applied to the reconstruction of shapes of buried objects from experimental backscattering data. We use experimental data which are collected by a microwave scattering facility located at the University of North Carolina at Charlotte, USA. Our numerical tests show accurate imaging of three components of interest of targets: shapes, locations and refractive indices.

Michael V. Klibanov

For the first time, a globally convergent numerical method is developed and Lipschitz stability estimate is obtained for the challenging problem of travel time tomography in 3D for formally determined incomplete data. The semidiscrete case is considered meaning that finite differences are involved with respect to two out of three variables. First, Lipschitz stability estimate is derived, which implies uniqueness. Next, a weighted globally strictly convex Tikhonov-like functional is constructed using a Carleman-like weight function for a Volterra integral operator. The gradient projection method is constructed to minimize this functional. It is proven that this method converges globally to the exact solution if the noise in the data tends to zero.

Alexey V. Smirnov, Michael V. Klibanov, Loc H. Nguyen

A new numerical method to solve an inverse source problem for the radiative transfer equation involving the absorption and scattering terms, with incomplete data, is proposed. No restrictive assumption on those absorption and scattering coefficients is imposed. The original inverse source problem is reduced to boundary value problem for a system of coupled partial differential equations of the first order. The unknown source function is not a part of this system. Next, we write this system in the fully discrete form of finite differences. That discrete problem is solved via the quasi-reversibility method. We prove the existence and uniqueness of the regularized solution. Especially, we prove the convergence of regularized solutions to the exact one as the noise level in the data tends to zero via a new discrete Carleman estimate. Numerical simulations demonstrate good performance of this method even when the data is highly noisy.

Michael V. Klibanov, Loc H. Nguyen

A new numerical method for X-ray tomography for a specific case of incomplete Radon data is proposed. Potential applications are in checking out bulky luggage in airports. This method is based on the analysis of the transport PDE governing the X-ray tomography rather than on the conventional integral formulation. The quasi-reversibility method is applied. Convergence analysis is performed using a new Carleman estimate. Numerical results are presented and compared with the inversion of the Radon transform using the well-known filtered back projection algorithm. In addition, it is shown how to use our method to study the inversion of the attenuated X-ray transform for the same case of incomplete data.

Michael V. Klibanov

By our definition, "restricted Dirichlet-to-Neumann map" (DN) means that the Dirichlet and Neumann boundary data for a Coefficient Inverse Problem (CIP) are generated by a point source running along an interval of a straight line. On the other hand, the conventional DN data can be generated, at least sometimes, by a point source running along a hypersurface. CIPs with the restricted DN data are non-overdetermined in the $n-$D case with $n \geq 2$. We develop, in a unified way, a general and a radically new numerical concept for CIPs with restricted DN data for a broad class of PDEs of the second order, such as, e.g. elliptic, parabolic and hyperbolic ones. Namely, using Carleman Weight Functions, we construct globally convergent numerical methods. Hölder stability and uniqueness are also proved. The price we pay for these features is a well acceptable one in the Numerical Analysis: we truncate a certain Fourier-like series with respect to some functions depending only on the position of that point source. At least three applications are: imaging of land mines, crosswell imaging and electrical impedance tomography.