Search

A system of two coupled nonlinear parabolic partial differential equations with two opposite directions of time is considered. In fact, this is the so-called "Mean Field Games System" (MFGS), which is derived in the mean field games (MFG) theory. This theory has numerous applications in social sciences. The topic of Coefficient Inverse Problems (CIPs) in the MFG theory is in its infant age, both in theory and computations. A numerical method for this CIP is developed. Convergence analysis ensures the global convergence of this method. Numerical experiments are presented., Comment: arXiv admin note: text overlap with arXiv:2310.08878

The time dependent experimental data are always collected at discrete grids with respect to the time t. The step size h of such a grid is always separated from zero by a certain positive number. The same is true for all computations, which are always done on discrete grids with their grid step sizes being not too small. These applied considerations prompt us to introduce a new type of Ill-Posed Problems and Coefficient Inverse Problems (CIP)for parabolic equations. In these problems the t-derivatives of corresponding parabolic operators are written in finite differences with the grid step size being separated from zero. We call this the "t-finite difference framework" (TFD). We address three long standing open questions within the TFD framework. Finally, a numerical method is developed for the CIP of monitoring epidemics. The global convergence of this method is proven.

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.

H\"older stability estimate and uniqueness are proven for a retrospective problem of Mean Field Games with a non-quadratic Hamiltonian. The previous result was only for the quadratic Hamiltonian. The main tool is the apparatus of Carleman estimates. The Mean Field Games theory has a broad range of applications in socio-economic sciences., Comment: The paper needs to be modified

The globally convergent convexification numerical method is constructed for a Coefficient Inverse Problem for the Mean Field Games System. A coefficient characterizing the global interaction term is recovered from the single measurement data. In particular, a new Carleman estimate for the Volterra integral operator is proven, and it stronger than the previously known one. Numerical results demonstrate accurate reconstructions from noisy data.

The convexification numerical method with the rigorously established global convergence property is constructed for a problem for the Mean Field Games System of the second order. This is the problem of the retrospective analysis of a game of infinitely many rational players. In addition to traditional initial and terminal conditions, one extra terminal condition is assumed to be known. Carleman estimates and a Carleman Weight Function play the key role. Numerical experiments demonstrate a good performance for complicated functions. Various versions of the convexification have been actively used by this research team for a number of years to numerically solve coefficient inverse problems.

A Coefficient Inverse Problem (CIP) of the determination of a coefficient of the Mean Field Games System (MFGS) of the second order is considered. The input data are generated by a single measurement event. Lateral Cauchy data, i.e. Dirichlet and Neumann boundary data are given for solutions of both equations forming the MFGS. In addition, solutions of both these equations are given at a fixed moment of time. H\"older stability estimates are obtained for both complete and incomplete lateral Cauchy data. These estimates imply uniqueness of our CIP. The apparatus of Carleman estimates is the main mathematical tool of this paper.

The mean field games (MFG) theory has broad application in mathematical modeling of social phenomena. The Mean Field Games System (MFGS) is the key to the MFG theory. This is a system of two nonlinear parabolic partial differential equations with two opposite directions of time 0

Our objective is to calculate the derivatives of data corrupted by noise. This is a challenging task as even small amounts of noise can result in significant errors in the computation. This is mainly due to the randomness of the noise, which can result in high-frequency fluctuations. To overcome this challenge, we suggest an approach that involves approximating the data by eliminating high-frequency terms from the Fourier expansion of the given data with respect to the polynomial-exponential basis. This truncation method helps to regularize the issue, while the use of the polynomial-exponential basis ensures accuracy in the computation. We demonstrate the effectiveness of our approach through numerical examples in one and two dimensions.

We are concerned with the mathematical study of the Mean Field Games system (MFGS). In the conventional setup, the MFGS is a system of two coupled nonlinear parabolic PDEs of the second order in a backward-forward manner, namely one terminal and one initial conditions are prescribed respectively for the value function and the population density. In this paper, we show that uniqueness of solutions to the MFGS can be guaranteed if, among all four possible terminal and initial conditions, eitheir only two terminal or only two initial conditions are given. In both cases H\"older stability estimates are proven. This means that he accuracies of the solutions are estimated in terms of the given data. Moreover, these estimates readily imply uniqueness of corresponding problems for the MFGS. The main mathematical apparatus to establish those results is two new Carleman estimates, which may find application in other contexts associated with coupled parabolic PDEs.

The second order Mean Field Games system (MFGS) in a bounded domain with the lateral Cauchy data is considered. This means that both Dirichlet and Neumann boundary data for the solution the MFGS are given. Two H\"older stability estimates for two slightly diffeent cases are derived. These estimates indicate how stable the solution of the MFGS is with respect to the possible noise in the lateral Cauchy data. Our stability estimates imply uniqueness. The key mathematical apparatus is the apparatus of two new Carleman estimates.

An overdetermination is introduced in an initial condition for the second order mean field games system (MFGS). This makes the resulting problem close to the classical ill-posed Cauchy problems for PDEs. Indeed, in such a problem and overdetermination in boundary conditions usually takes place. A Lipschitz stability estimate is obtained. This estimate implies uniqueness. A new Carleman estimate is derived. The second estimate the author calls "quasi-Carleman estimate", since it contains two test functions rather than a single one in conventional Carleman estimates. These two estimates play the key role. Carleman estimates were not applied to the MFGS prior to the recent work of Klibanov and Aveboukh in arXiv: 2302.10709, 2023. Applications are discussed.

A Coefficient Inverse Problem for the radiative transport equation is considered. The globally convergent numerical method, the so-called convexification, is developed. For the first time, the viscosity solution is considered for a boundary value problem for the resulting system of two coupled partial differential equations. A Lipschitz stability estimate is proved for this boundary value problem using a Carleman estimate for the Laplace operator. Next, the global convergence analysis is provided via that Carleman estimate. Results of numerical experiments demonstrate a high computational efficiency of this approach.

A retrospective analysis process for the mean field games system (MFGS) is considered. For the first time, Carleman estimates are applied to the analysis of the MFGS. Two new Carleman estimates are derived. They allow to obtain the Lipschitz stability estimate with respect to the possible error in the input initial and terminal data for a retrospective problem for MFGS. This stability estimate, in turn implies uniqueness theorem for the problem under the consideration. The idea of using Carleman estimates to obtain stability and uniqueness results came from the field of Ill-Posed and Inverse Problems.

This paper is about Holder and Lipschitz stability estimates and uniqueness theorems for some coefficient inverse problems and associated inverse source problems for a general linear parabolic equation of the second order with variable coefficients. The data for the inverse problem are given at the final moment of time {t=T}. In addition, both Dirichlet and Neumann boundary conditions are given either on a part or on the entire lateral boundary. Thus, if these boundary conditions are given only at a part of the boundary, then even if the target coefficient is known, still the forward problem is not a classical initial boundary value problem.

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.

The first globally convergent numerical method for a Coefficient Inverse Problem (CIP) for the Riemannian Radiative Transfer Equation (RRTE) is constructed. This is a version of the so-called \textquotedblleft convexification" method, which has been pursued by this research group for a number of years for some other CIPs for PDEs. Those PDEs are significantly different from the RRTE. The presence of the Carleman Weight Function (CWF) in the numerical scheme is the key element which insures the global convergence. Convergence analysis is presented along with the results of numerical experiments, which confirm the theory. RRTE governs the propagation of photons in the diffuse medium in the case when they propagate along geodesic lines between their collisions. Geodesic lines are generated by the spatially variable dielectric constant of the medium.

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.

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\"older 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.

In the previous paper (Inverse Problems, 32, 015010, 2016), a new heuristic mathematical model was proposed for accurate forecasting of prices of stock options for 1-2 trading days ahead of the present one. This new technique uses the Black-Scholes equation supplied by new intervals for the underlying stock and new initial and boundary conditions for option prices. The Black-Scholes equation was solved in the positive direction of the time variable, This ill-posed initial boundary value problem was solved by the so-called Quasi-Reversibility Method (QRM). This approach with an added trading strategy was tested on the market data for 368 stock options and good forecasting results were demonstrated. In the current paper, we use the geometric Brownian motion to provide an explanation of that effectivity using computationally simulated data for European call options. We also provide a convergence analysis for QRM. The key tool of that analysis is a Carleman estimate.

A version of the convexification globally convergent numerical method is constructed for a coefficient inverse problem for a wave-like partial differential equation. The presence of the Carleman Weight Function in the corresponding Tikhonov-like cost functional ensures the global strict convexity of this functional. Numerical results are presented to illustrate the effectiveness and efficiency of the proposed method.

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.

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.

To compute the spatially distributed dielectric constant from the backscattering data, we study a coefficient inverse problem for a 1D hyperbolic equation. To solve the inverse problem, we establish a new version of Carleman estimate and then employ this estimate to construct a cost functional which is strictly convex on a convex bounded set with an arbitrary diameter in a Hilbert space. The strict convexity property is rigorously proved. This result is called the convexification theorem and is considered as the central analytical result of this paper. Minimizing this convex functional by the gradient descent method, we obtain the desired numerical solution to the coefficient inverse problems. We prove that the gradient descent method generates a sequence converging to the minimizer and we also establish a theorem confirming that the minimizer converges to the true solution as the noise in the measured data and the regularization parameter tend to zero. Unlike the methods that are based on optimization, our convexification method converges globally in the sense that it delivers a good approximation of the exact solution without requiring any initial guess. Results of numerical studies of both computationally simulated and experimental data are presented.

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., Comment: 20 pages

For the first time, we develop a convergent numerical method for the llinear integral equation derived by M.M. Lavrent'ev in 1964 with the goal to solve a coefficient inverse problem for a wave-like equation in 3D. The data are non overdetermined. Convergence analysis is presented along with the numerical results. An intriguing feature of the Lavrent'ev equation is that, without any linearization, it reduces a highly nonlinear coefficient inverse problem to a linear integral equation of the first kind. Nevertheless, numerical results for that equation, which use the data generated for that coefficient inverse problem, show a good reconstruction accuracy. This is similar with the classical Gelfand-Levitan equation derived in 1951, which is valid in the 1D case.

An inverse scattering problem for SAR data in application to through-the-wall imaging is addressed. In contrast with the conventional algorithms of SAR imaging, that work with the linearized mathematical model based on the Born approximation, the fully nonlinear case is considered here. To avoid the local minima problem, the so-called "convexification" globally convergent inversion scheme is applied to approximate the distribution of the slant range (SR) dielectric constant in the 3-D domain. The benchmark scene of this paper comprises a homogeneous dielectric wall and different dielectric targets hidden behind it. The results comprise two dimensional images of the SR dielectric constant of the scene of interest. Numerical results are obtained by the proposed inversion method for both the computationally simulated and experimental data. Our results show that the values, cross-range sizes and locations of SR dielectric constants for targets hidden behind the wall are close to those of real targets. Numerical comparison with the solution of the linearized inverse scattering problem provided by the Born approximation, commonly used in conventional SAR imaging, shows a significantly better accuracy of our results.

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., Comment: Coefficient inverse problem, multiple point sources, experimental data, Carleman weight, global convergence, Fourier series. arXiv admin note: text overlap with arXiv:2003.11513

This paper is concerned with the convergence of a series associated with a certain version of the convexification method. That version has been recently developed by the research group of the first author for solving coefficient inverse problems. The convexification method aims to construct a globally convex Tikhonov-like functional with a Carleman Weight Function in it. In the previous works the construction of the strictly convex weighted Tikhonov-like functional assumes a truncated Fourier series (i.e. a finite series instead of an infinite one) for a function generated by the total wave field. In this paper we prove a convergence property for this truncated Fourier series approximation. More precisely, we show that the residual of the approximate PDE obtained by using the truncated Fourier series tends to zero in $L^{2}$ as the truncation index in the truncated Fourier series tends to infinity. The proof relies on a convergence result in the $H^{1}$-norm for a sequence of $L^{2}$-orthogonal projections on finite-dimensional subspaces spanned by elements of a special Fourier basis. However, due to the ill-posed nature of coefficient inverse problems, we cannot prove that the solution of that approximate PDE, which results from the minimization of that Tikhonov-like functional, converges to the correct solution., Comment: The proof of Lemma 3 contains an uncorrectable mistake

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., Comment: 35 pages, 22 figures, 2 tables

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., Comment: 22 pages, 4 figures, accepted for publication by Inverse Problems and Imaging

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.

We propose a new numerical method to solve the linearized problem of travel time tomography with incomplete data. Our method is based on the technique of the truncation of the Fourier series with respect to a special basis of L2. This way we derive a boundary value problem for a system of coupled partial diffeerential equations (PDEs) of the first order. This system is solved by the quasi-reversibility method. Hence, the spatially dependent Fourier coefficients of the solution to the linearized Eikonal equation are obtained. The convergence of this method is established. Numerical results for highly noisy data are presented.

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.

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.

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., Comment: 22 pages, 4 figures, to be published in SIAM Journal on Applied Mathematics

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.

A new numerical method to solve an inverse source problem for the Helmholtz equation in inhomogenous media is proposed. This method reduces the original inverse problem to a boundary value problem for a coupled system of elliptic PDEs, in which the unknown source function is not involved. The Dirichlet boundary condition is given on the entire boundary of the domain of interest and the Neumann boundary condition is given on a part of this boundary. To solve this problem, the quasi-reversibility method is applied. Uniqueness and existence of the minimizer are proven. A new Carleman estimate is established. Next, the convergence of those minimizers to the exact solution is proven using that Carleman estimate. Results of numerical tests are presented.

We propose a new approach to constructing globally strictly convex objective functional in a 1-D inverse medium scattering problem using multi-frequency backscattering data. The global convexity of the proposed objective functional is proved using a Carleman estimate. Due to its convexity, no good first guess is required in minimizing this objective functional. We also prove the global convergence of the gradient projection algorithm and derive an error estimate for the reconstructed coefficient. Numerical results show reasonable reconstruction accuracy for simulated data., Comment: 30 pages

An inverse scattering problem for the 3D acoustic equation in time domain is considered. The unknown spatially distributed speed of sound is the subject of the solution of this problem. A single location of the point source is used. Using a Carleman Weight Function, a globally strictly convex cost functional is constructed. A new Carleman estimate is proven. Global convergence of the gradient projection method is proven. Numerical experiments are conducted.

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.

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.

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., Comment: 30 pages, 5 figures

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., Comment: 3 Pages, 7 Figures

We propose in this paper a globally numerical method to solve a phaseless coefficient inverse problem: how to reconstruct the spatially distributed refractive index of scatterers from the intensity (modulus square) of the full complex valued wave field at an array of light detectors located on a measurement board. The propagation of the wave field is governed by the 3D Helmholtz equation. Our method consists of two stages. On the first stage, we use asymptotic analysis to obtain an upper estimate for the modulus of the scattered wave field. This estimate allows us to approximately reconstruct the wave field at the measurement board using an inversion formula. This reduces the phaseless inverse scattering problem to the phased one. At the second stage, we apply a recently developed globally convergent numerical method to reconstruct the desired refractive index from the total wave obtained at the first stage. Unlike the optimization approach, the two-stage method described above is global in the sense that it does not require a good initial guess of the true solution. We test our numerical method on both computationally simulated and experimental data. Although experimental data are noisy, our method produces quite accurate numerical results.

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., Comment: 27 pages, 6 figures

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.