Researcher profile

Michael V. Klibanov

Michael V. Klibanov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2024arXiv

Spatiotemporal Monitoring of Epidemics via Solution of a Coefficient Inverse Problem

Let S,I and R be susceptible, infected and recovered populations in a city affected by an epidemic. The SIR model of Lee, Liu, Tembine, Li and Osher, \emph{SIAM J. Appl. Math.},~81, 190--207, 2021 of the spatiotemoral spread of epidemics is considered. This model consists of a system of three nonlinear coupled parabolic Partial Differential Equations with respect to the space and time dependent functions S,I and R. For the first time, a Coefficient Inverse Problem (CIP) for this system is posed. The so-called \textquotedblleft convexification" numerical method for this inverse problem is constructed. The presence of the Carleman Weight Function (CWF) in the resulting regularization functional ensures the global convergence of the gradient descent method of the minimization of this functional to the true solution of the CIP, as long as the noise level tends to zero. The CWF is the function, which is used as the weight in the Carleman estimate for the corresponding Partial Differential Operator. Numerical studies demonstrate an accurate reconstruction of unknown coefficients as well as S,I,R functions inside of that city. As a by-product, uniqueness theorem for this CIP is proven. Since the minimal measured input data are required, then the proposed methodology has a potential of a significant decrease of the cost of monitoring of epidemics.

preprint2023arXiv

A comparative study of two globally convergent numerical methods for acoustic tomography

The comparative study of two globally convergent numerical methods for acoustic tomography is carried out in two dimensions. These are the boundary control method and the quasi-reversibility method. The novelty is that in the latter a nonlinear inverse problem is reduced to a family of the linear integral equations of the first kind via the Lavrentiev approach and this reduction is used within the quasi-reversibility method. The analysis of its stability is carried out via Carleman estimates. The computational effectiveness of these methods is tested in the numerical experiments with the smooth and discontinuous coefficients to be recovered from the tomographic data.

preprint2022arXiv

A Hölder Stability Estimate for a 3D Coefficient Inverse Problem for a Hyperbolic Equation With a Plane Wave

A 3D coefficient inverse problem for a hyperbolic equation with non-overdetermined data is considered. The forward problem is the Cauchy problems with the initial condition the delta function concentrated at a single plane (i.e. the plane wave). A certain associated operator is written in finite differences with respect to two out of three spatial variables, i.e. "partial finite differences". The grid step size is bounded from the below by a fixed number. A Carleman estimate is applied to obtain, for the first time, a Hölder stability estimates for this problem. Another new result is an estimate from the below of the amplitude of the first term of the expansion of the solution of the forward problem near the characteristic wedge.

preprint2022arXiv

Carleman estimates and the contraction principle for an inverse source problem for nonlinear hyperbolic equation

The main aim of this paper is to solve an inverse source problem for a general nonlinear hyperbolic equation. Combining the quasi-reversibility method and a suitable Carleman weight function, we define a map of which fixed point is the solution to the inverse problem. To find this fixed point, we define a recursive sequence with an arbitrary initial term by the same manner as in the classical proof of the contraction principle. Applying a Carleman estimate, we show that the sequence above converges to the desired solution with the exponential rate. Therefore, our new method can be considered as an analog of the contraction principle. We rigorously study the stability of our method with respect to noise. Numerical examples are presented.

preprint2022arXiv

Convexification for a CIP for the RTE]{Convexification Numerical Method for a Coefficient Inverse Problem for the Radiative Transport Equation

An $\left( n+1\right) -$D coefficient inverse problem for the radiative stationary transport equation is considered for the first time. A globally convergent so-called convexification numerical \ method is developed and its convergence analysis is provided. The analysis is based on a Carleman estimate. In particular, convergence analysis implies a certain uniqueness theorem. Extensive numerical studies in the 2-D case are presented.

preprint2021arXiv

Can a single PDE govern well the propagation of the electric wave field in a heterogeneous medium in 3D?

It is demonstrated in this paper that the propagation of the electric wave field in a heterogeneous medium in 3D can sometimes be governed well by a single PDE, which is derived from the Maxwell's equations. The corresponding component of the electric field dominates two other components. This justifies some past results of the second author with coauthors about numerical solutions of coefficient inverse problems with experimental electromagnetic data. In addition, since it is simpler to work in applications with a single PDE rather than with the complete Maxwell's system, then the result of this paper might be useful to researchers working on applied issues of the propagation of electromagnetic waves in inhomogeneous media.

preprint2021arXiv

Carleman contraction mapping for a 1D inverse scattering problem with experimental time-dependent data

It is shown that the contraction mapping principle with the involvement of a Carleman Weight Function works for a Coefficient Inverse Problem for a 1D hyperbolic equation. Using a Carleman estimate, the global convergence of the corresponding numerical method is established. Numerical studies for both computationally simulated and experimentally collected data are presented. The experimental part is concerned with the problem of computing dielectric constants of explosive-like targets in the standoff mode using severely underdetermined data.

preprint2020arXiv

An inverse problem of a simultaneous reconstruction of the dielectric constant and conductivity from experimental backscattering data

This report extends our recent progress in tackling a challenging 3D inverse scattering problem governed by the Helmholtz equation. Our target application is to reconstruct dielectric constants, electric conductivities and shapes of front surfaces of objects buried very closely under the ground. These objects mimic explosives, like, e.g., antipersonnel land mines and improvised explosive devices. We solve a coefficient inverse problem with the backscattering data generated by a moving source at a fixed frequency. This scenario has been studied so far by our newly developed convexification method that consists in a new derivation of a boundary value problem for a coupled quasilinear elliptic system. However, in our previous work only the unknown dielectric constants of objects and shapes of their front surfaces were calculated. Unlike this, in the current work performance of our numerical convexification algorithm is verified for the case when the dielectric constants, the electric conductivities and those shapes of objects are unknown. By running several tests with experimentally collected backscattering data, we find that we can accurately image both the dielectric constants and shapes of targets of interests including a challenging case of targets with voids. The computed electrical conductivity serves for reliably distinguishing conductive and non-conductive objects. The global convergence of our numerical procedure is shortly revisited.

preprint2020arXiv

Convexification and experimental data for a 3D inverse scattering problem with the moving point source

Inverse scattering problems of the reconstructions of physical properties of a medium from boundary measurements are substantially challenging ones. This work aims to verify the performance on experimental data of a newly developed convexification method for a 3D coefficient inverse problem for the case of objects buried in a sandbox a fixed frequency and the point source moving along an interval of a straight line. Using a special Fourier basis, the method of this work strongly relies on a new derivation of a boundary value problem for a system of coupled quasilinear elliptic equations. This problem, in turn, is solved via the minimization of a Tikhonov-like functional weighted by a Carleman Weight Function. The global convergence of the numerical procedure is established analytically. The numerical verification is performed using experimental data, which are raw backscatter data of the electric field. These data were collected using a microwave scattering facility at The University of North Carolina at Charlotte.

preprint2020arXiv

Convexification for a 1D Hyperbolic Coefficient Inverse Problem with Single Measurement Data

A version of the convexification numerical method for a Coefficient Inverse Problem for a 1D hyperbolic PDE is presented. The data for this problem are generated by a single measurement event. This method converges globally. The most important element of the construction is the presence of the Carleman Weight Function in a weighted Tikhonov-like functional. This functional is strictly convex on a certain bounded set in a Hilbert space, and the diameter of this set is an arbitrary positive number. The global convergence of the gradient projection method is established. Computational results demonstrate a good performance of the numerical method for noisy data.

preprint2020arXiv

Convexification for an Inverse Parabolic Problem

A convexification-based numerical method for a Coefficient Inverse Problem for a parabolic PDE is presented. The key element of this method is the presence of the so-called Carleman Weight Function in the numerical scheme. Convergence analysis ensures the global convergence of this method, as opposed to the local convergence of the conventional least squares minimization techniques. Numerical results demonstrate a good performance.

preprint2019arXiv

Convergent numerical methods for parabolic equations with reversed time via a new Carleman estimate

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.

preprint2013arXiv

A globally convergent algorithm for the frequency sounding and Slichter-Langer-Tikhonov problem of electrical prospecting

The paper presents a globally convergent algorithm for solving coefficient inverse problems. Being rooted in the globally convergent numerical method (SIAM J. Sci. Comput., 31, No.1 (2008), pp. 478-509) for solving multidimensional coefficient inverse problems, it has two distinctive features: the new iterative and refinement procedures. These novelties enhance, sometimes significantly, both the spatial and contrast resolutions. The computational effectiveness of the proposed technique is demonstrated in numerical experiments with two applied coefficient inverse problems: electromagnetic or acoustic frequency sounding and electrical prospecting of layered media. The Slichter-Langer-Tikhonov formulation is exploited as a mathematical model of the latter.

preprint2013arXiv

Reconstruction of the refractive index from experimental backscattering data using a globally convergent inverse method

The problem to be studied in this work is within the context of coefficient identification problems for the wave equation. More precisely, we consider the problem of reconstruction of the refractive index (or equivalently, the dielectric constant) of an inhomogeneous medium using one backscattering boundary measurement. The goal of this paper is to analyze the performance of a globally convergent algorithm of Beilina and Klibanov on experimental data acquired in the Microwave Laboratory at University of North Carolina at Charlotte. The main challenge working with experimental data is the the huge misfit between these data and computationally simulated data. We present data pre-processing steps to make the former somehow look similar to the latter. Results of both non-blind and blind targets are shown indicating good reconstructions even for high contrasts between the targets and the background medium.

preprint2013arXiv

The Gel'fand-Levitan-Krein method and the globally convergent method for experimental data

Comparison of numerical performances of two methods for coefficient inverse problems is described. The first one is the classical Gel'fand-Levitan-Krein equation method, and the second one is the recently developed approximately globally convergent numerical method. This comparison is performed for both computationally simulated and experimental data.

preprint2012arXiv

A new approximate mathematical model for global convergence for a coefficient inverse problem with backscattering data

An approximately globally convergent numerical method for a 3d Coefficient Inverse Problem for a hyperbolic equation with backscattering data is presented. A new approximate mathematical model is presented. An approximation is used only on the first iteration and amounts to the truncation of a certain asymptotic series. A significantly new element of the convergence analysis is that the so-called "tail functions" are estimated. Numerical results in 2d and 3d cases are presented, including the one for a quite heterogeneous medium.

preprint2012arXiv

The adaptivity refines approximate solutions of ill-posed problems due to the relaxation property

Adaptive Finite Element Method (adaptivity) is known to be an effective numerical tool for some ill-posed problems. The key advantage of the adaptivity is the image improvement with local mesh refinements. A rigorous proof of this property is the central part of this paper. In terms of Coefficient Inverse Problems with single measurement data, the authors consider the adaptivity as the second stage of a two-stage numerical procedure. The first stage delivers a good approximation of the exact coefficient without an advanced knowledge of a small neighborhood of that coefficient. This is a necessary element for the adaptivity to start iterations from. Numerical results for the two-stage procedure are presented for both computationally simulated and experimental data.