A Globally Convergent Algorithm for a PDE-Constrained Optimization Problem Arising in Electrical Impedance Tomography

Abstract

We study the convergence properties of an algorithm for the inverse problem of electrical impedance tomography, which can be reduced to a partial differential equation (PDE) constrained optimization problem. The direct problem consists of the potential equation div( u) = 0 in a circle, with Neumann condition describing the behavior of the electrostatic potential in a medium with conductivity given by the function (x, y). We suppose that at each time a current i is applied to the boundary of the circle (Neumann's data), and that it is possible to measure the corresponding potential i (Dirichlet data). The inverse problem is to find (x, y) given a finite number of Cauchy pairs measurements ( i , i ), i = 1, , N . The problem is formulated as a least squares problem, and the developed algorithm solves the continuous problem using descent iterations in its corresponding finite element approximations. Wolfe's conditions are used to ensure the global convergence of the optimization algorithm for the continuous problem. Although exact data are assumed, measurement errors in data and regularization methods shall be considered in a future work.

Nonlinear Electrical Impedance Tomography Method Using a Complete Electrode Model for the Characterization of Heterogeneous Domains

“…The proposed inverse EIT problem deals with the Laplace equation, Robin, and Neumann boundary conditions in the PDE-constrained optimization. There are similar EIT approaches developed so far, but they mostly deal with Neumann or Dirichlet boundary conditions to solve the forward problem [17][18][19]. Some recent studies have included the Robin boundary condition to reflect the contact impedance between the domain and electrode boundaries, but their inverse formulation and optimization schemes are different from the present study [20,21].…”