scholarly journals Convexification-based globally convergent numerical method for a 1D coefficient inverse problem with experimental data

2021 ◽  
Vol 0 (0) ◽  
pp. 0
Author(s):  
Michael V. Klibanov ◽  
Thuy T. Le ◽  
Loc H. Nguyen ◽  
Anders Sullivan ◽  
Lam Nguyen

<p style='text-indent:20px;'>To compute the spatially distributed dielectric constant from the backscattering computationally simulated ane experimentally collected data, we study a coefficient inverse problem for a 1D hyperbolic equation. To solve this inverse problem, we establish a new version of the Carleman estimate and then employ this estimate to construct a cost functional, which is strictly convex on a convex bounded set of an arbitrary diameter in a Hilbert space. The strict convexity property is rigorously proved. This result is called the convexification theorem and it is the central analytical result of this paper. Minimizing this cost 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 starting from an arbitrary point of that bounded set. 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, which are based on the optimization, our convexification method converges globally in the sense that it delivers a good approximation of the exact solution without requiring a good initial guess. Results of numerical studies of both computationally simulated and experimentally collected data are presented.</p>


Computation ◽  
2020 ◽  
Vol 8 (3) ◽  
pp. 73
Author(s):  
Dmitriy Klyuchinskiy ◽  
Nikita Novikov ◽  
Maxim Shishlenin

We investigate the mathematical model of the 2D acoustic waves propagation in a heterogeneous domain. The hyperbolic first order system of partial differential equations is considered and solved by the Godunov method of the first order of approximation. This is a direct problem with appropriate initial and boundary conditions. We solve the coefficient inverse problem (IP) of recovering density. IP is reduced to an optimization problem, which is solved by the gradient descent method. The quality of the IP solution highly depends on the quantity of IP data and positions of receivers. We introduce a new approach for computing a gradient in the descent method in order to use as much IP data as possible on each iteration of descent.









Sign in / Sign up

Export Citation Format

Share Document