English

Convexification numerical algorithm for a 2D inverse scattering problem with backscatter data

Numerical Analysis 2020-02-25 v2 Numerical Analysis

Abstract

This paper is concerned with the inverse scattering problem which aims to determine the spatially distributed dielectric constant coefficient of the 2D Helmholtz equation from multifrequency backscatter data associated with a single direction of the incident plane wave. We propose a globally convergent convexification numerical algorithm to solve this nonlinear and ill-posed inverse problem. The key advantage of our method over conventional optimization approaches is that it does not require a good first guess about the solution. First, we eliminate the coefficient from the Helmholtz equation using a change of variables. Next, using a truncated expansion with respect to a special Fourier basis, we approximately reformulate the inverse problem as a system of quasilinear elliptic PDEs, which can be numerically solved by a weighted quasi-reversibility approach. The cost functional for the weighted quasi-reversibility method is constructed as a Tikhonov-like functional that involves a Carleman Weight Function. Our numerical study shows that, using a version of the gradient descent method, one can find the minimizer of this Tikhonov-like functional without any advanced \emph{a priori} knowledge about it.

Keywords

Cite

@article{arxiv.2002.08427,
  title  = {Convexification numerical algorithm for a 2D inverse scattering problem with backscatter data},
  author = {Trung Truong and Dinh-Liem Nguyen and Michael Klibanov},
  journal= {arXiv preprint arXiv:2002.08427},
  year   = {2020}
}

Comments

25 pages