BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael V. Klibanov

Michael V. Klibanov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.NAmath.APNumerical Analysismath.OCphysics.class-phq-fin.MF

Catalog footprint

What is connected

46works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

46 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A globally convergent method for a 3-D inverse medium problem for the generalized Helmholtz equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A globally convergent numerical method for a 1-d inverse medium problem with experimental data

In this paper, a reconstruction method for the spatially distributed dielectric constant of a medium from the back scattering wave field in the frequency domain is considered. Our approach is to propose a globally convergent algorithm, which does not require any knowledge of a small neighborhood of the solution of the inverse problem in advance. The Quasi-Reversibility Method (QRM) is used in the algorithm. The convergence of the QRM is proved via a Carleman estimate. The method is tested on both computationally simulated and experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A globally convergent numerical method for a 3D coefficient inverse problem with a single measurement of multi-frequency data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Carleman weight functions for a globally convergent numerical method for ill-posed Cauchy problems for some quasilinear PDEs

In a series of publications of the second author, including some with coauthors, globally strictly convex Tikhonov-like functionals were constructed for some nonlinear ill-posed problems. The main element of such a functional is the presence of the Carleman Weight Function. Compared with previous publications, the main novelty of this paper is that the existence of the regularized solution (i.e. the minimizer) is proved rather than assumed. The method works for both ill-posed Cauchy problems for some quasilinear PDEs of the second order and for some Coefficient Inverse Problems. However, to simplify the presentation, we focus here only on ill-posed Cauchy problems. Along with the theory, numerical results are presented for the case of a 1-D quasilinear parabolic PDE with the lateral Cauchy data given on one edge of the interval (0,1).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical solution of an ill-posed Cauchy problem for a quasilinear parabolic equation using a Carleman weight function

This is the first publication in which an ill-posed Cauchy problem for a quasi- linear PDE is solved numerically by a rigorous method. More precisely, we solve the side Cauchy problem for a 1-d quasilinear parabolc equation. The key idea is to minimize a strictly convex cost functional with the Carleman Weight Function in it. Previous publications about numerical solutions of ill-posed Cauchy problems were considering only linear equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Single measurement experimental data for an inverse medium problem inverted by a multi-frequency globally convergent numerical method

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nanostructures imaging via numerical solution of a 3-d inverse scattering problem without the phase information

Inverse scattering problems without the phase information arise in imaging of nanostructures whose sizes are hundreds of nanometers as well as in imaging of biological cells. The governing equation is the 3-d generalized Helmholtz equation with the unknown coefficient, which represents the spatially distributed dielectric constant. It is assumed that only the modulus of the complex valued wave field is measured on a frequency range. The phase is not measured. Two rigorous numerical methods are tested and their performances are compared for realistic ranges of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical studies of an adaptive finite element method applied to the reconstruction of shapes of buried objects from experimental data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On a problem of S.L. Sobolev

In 1930 Sergey L. Sobolev [7,8] has proposed a construction of the solution of the Cauchy problem for the hyperbolic equation of the second order with variable coefficients in 3-d. Although Sobolev did not construct the fundamental solution, his construction was modified later by Romanov [4,5] to obtain the fundamental solution. However, these works impose a restrictive assumption of the regularity of geodesic lines in a large domain. In addition, it is unclear how to realize those methods numerically. In this paper a simple construction of a function, which is associated in a clear way with the fundamental solution of the acoustic equation with the variable speed in 3-d, is proposed. Conditions on geodesic lines are not imposed. An important feature of this construction is that it lends itself to effective computations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Profitable forecast of prices of stock options on real market data via the solution of an ill-posed problem for the Black-Scholes equation

A new mathematical model for the Black-Scholes equation is proposed to forecast option prices. This model includes new interval for the price of the underlying stock as well as new initial and boundary conditions. Conventional notions of maturity time and strike prices are not used. The Black-Scholes equation is solved as a parabolic equation with the reversed time, which is an ill-posed problem. Thus, a regularization method is used to solve it. This idea is verified on real market data for twenty liquid options. A trading strategy is proposed. This strategy indicates that our method is profitable on at least those twenty options. We conjecture that our method might lead to significant profits of those financial institutions which trade large amounts of options. We caution, however, that detailed further studies are necessary to verify this conjecture.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reconstruction procedures for two inverse scattering problems without the phase information

This is a continuation of two recent publications of the authors about reconstruction procedures for 3-d phaseless inverse scattering problems. The main novelty of this paper is that the Born approximation for the case of the wave-like equation is not considered. It is shown here that the phaseless inverse scattering problem for the 3-d wave-like equation in the frequency domain leads to the well known Inverse Kinematic Problem. Uniqueness theorem follows. Still, since the Inverse Kinematic Problem is very hard to solve, a linearization is applied. More precisely, geodesic lines are replaced with straight lines. As a result, an approximate explicit reconstruction formula is obtained via the inverse Radon transform. The second reconstruction method is via solving a problem of the integral geometry using integral equations of the Abel type.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two reconstruction procedures for a 3-d phaseless inverse scattering problem for the generalized Helmholtz equation

The 3-d inverse scattering problem of the reconstruction of the unknown dielectric permittivity in the generalized Helmholtz equation is considered. The main difference with the conventional inverse scattering problems is that only the modulus of the scattering wave field is measured. The phase is not measured. The initializing wave field is the incident plane wave. On the other hand, in the previous recent works of the authors about the "phaseless topic" the case of the point source was considered [20,21,22]. Two reconstruction procedures are developed for a linearized case. However, the linearization is not the Born approximation. This means that, unlike the Born approximation, our linearization does not break down when the frequency tends to the infinity. Applications are in imaging of nanostructures and biological cells.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Carleman estimates for the regularization of ill-posed Cauchy problems

This is a survey, which is a continuation of the previous survey of the author about applications of Carleman estimates to Inverse Problems, J. Inverse and Ill-Posed Problems, 21, 477-560, 2013. It is shown here that Tikhonov functionals for some ill-posed Cauchy problems for linear PDEs can be generated by unbounded linear operators of those PDEs. These are those operators for which Carleman estimates are valid, e.g. elliptic, parabolic and hyperbolic operators of the second order. Convergence rates of minimizers are established using Carleman estimates. Generalizations to nonlinear inverse problems, such as problems of reconstructions of obstacles and coefficient inverse problems are discussed as well.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Globally convergent and adaptive finite element methods in imaging of buried objects from experimental backscattering radar measurements

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Imaging of buried objects from experimental backscattering time dependent measurements using a globally convergent inverse algorithm

We consider the problem of imaging of objects buried under the ground using backscattering experimental time dependent measurements generated by a single point source or one incident plane wave. In particular, we estimate dielectric constants of those objects using the globally convergent inverse algorithm of Beilina and Klibanov. Our algorithm is tested on experimental data collected using a microwave scattering facility at the University of North Carolina at Charlotte. There are two main challenges working with this type of experimental data: (i) there is a huge misfit between these data and computationally simulated data, and (ii) the signals scattered from the targets may overlap with and be dominated by the reflection from the ground's surface. To overcome these two challenges, we propose new data preprocessing steps to make the experimental data to be approximately the same as the simulated ones, as well as to remove the reflection from the ground's surface. Results of total 25 data sets of both non blind and blind targets indicate a good accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the first solution of a long standing problem: uniqueness of the phaseless quantum inverse scattering problem in 3-d

After publishing his recent paper in SIAM J. Appl. Math, 74, 392-410, 2014 the author has realized that actually he has addressed in that paper, for the first time, a long standing open question being unaware about this. This question is about the uniqueness of a 3-d inverse scattering problem without the phase information. Thus, it makes sense in the current paper to explicitly make the latter statement and to formulate corresponding uniqueness theorems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reconstruction formula for a 3-d phaseless inverse scattering problem for the Schrodinger equation

The inverse scattering problem of the reconstruction of the unknown potential with compact support in the 3-d Schrödinger equation is considered. Only the modulus of the scattering complex valued wave field is known, whereas the phase is unknown. It is shown that the unknown potential can be reconstructed via the inverse Radon transform. Therefore, a long standing problem posed in 1977 by K. Chadan and P.C. Sabatier in their book "Inverse Problems in Quantum Scattering Theory" is solved.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reconstruction of shapes and refractive indices from backscattering experimental data using the adaptivity

We consider the inverse problem of the reconstruction of the spatially distributed dielectric constant $\varepsilon_{r}\left(\mathbf{x}\right), \ \mathbf{x}\in \mathbb{R}^{3}$, which is an unknown coefficient in the Maxwell's equations, from time-dependent backscattering experimental radar data associated with a single source of electric pulses. The refractive index is $n\left(\mathbf{x}\right) =\sqrt{\varepsilon_{r}\left(\mathbf{x}\right)}.$ The coefficient $\varepsilon_{r}\left(\mathbf{x}\right) $ is reconstructed using a two-stage reconstruction procedure. In the first stage an approximately globally convergent method proposed is applied to get a good first approximation of the exact solution. In the second stage a locally convergent adaptive finite element method is applied, taking the solution of the first stage as the starting point of the minimization of the Tikhonov functional. This functional is minimized on a sequence of locally refined meshes. It is shown here that all three components of interest of targets can be simultaneously accurately imaged: refractive indices, shapes and locations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recovering of dielectric constants of explosives via a globally strictly convex cost functional

The inverse problem of estimating dielectric constants of explosives using boundary measurements of one component of the scattered electric field is addressed. It is formulated as a coefficient inverse problem for a hyperbolic differential equation. After applying the Laplace transform, a new cost functional is constructed and a variational problem is formulated. The key feature of this functional is the presence of the Carleman Weight Function for the Laplacian. The strict convexity of this functional on a bounded set in a Hilbert space of an arbitrary size is proven. This allows for establishing the global convergence of the gradient descent method. Some results of numerical experiments are presented.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Globally strongly convex cost functional for a coefficient inverse problem

A Carleman Weight Function (CWF) is used to construct a new cost functional for a Coefficient Inverse Problems for a hyperbolic PDE. Given a bounded set of an arbitrary size in a certain Sobolev space, one can choose the parameter of the CWF in such a way that the constructed cost functional will be strongly convex on that set. Next, convergence of the gradient method, which starts from an arbitrary point of that set, is established. Since restrictions on the size of that set are not imposed, then this is the global convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reconstruction from blind experimental data for an inverse problem for a hyperbolic equation

We consider the problem of reconstruction of dielectrics from blind backscattered experimental data. Experimental data were collected by a device, which was built at University of North Carolina at Charlotte. This device sends electrical pulses into the medium and collects the time resolved backscattered data on a part of a plane. The spatially distributed dielectric constant $\varepsilon_{r}(\mathbf{x}),\mathbf{x}\in \mathbb{R}^{3}$ is the unknown coefficient of a wave-like PDE. This coefficient is reconstructed from those data in blind cases. To do this, a globally convergent numerical method is used.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optical imaging of phantoms from real data by an approximately globally convergent inverse algorithm

A numerical method for an inverse problem for an elliptic equation with the running source at multiple positions is presented. This algorithm does not rely on a good first guess for the solution. The so-called "approximate global convergence" property of this method is shown here. The performance of the algorithm is verified on real data for Diffusion Optical Tomography. Direct applications are in near-infrared laser imaging technology for stroke detection in brains of small animals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermoacoustic tomography with an arbitrary elliptic operator

Thermoacoustic tomography is a term for the inverse problem of determining of one of initial conditions of a hyperbolic equation from boundary measurements. In the past publications both stability estimates and convergent numerical methods for this problem were obtained only under some restrictive conditions imposed on the principal part of the elliptic operator. In this paper logarithmic stability estimates are obatined for an arbitrary variable principal part of that operator. Convergence of the Quasi-Reversibility Method to the exact solution is also established for this case. Both complete and incomplete data collection cases are considered.

0 citations0 reviewsNo code linkedOpen access