Recovery of region boundaries of piecewise constant coefficients of an elliptic PDE from boundary data

In this study we consider the recovery of smooth region boundaries of piece wise constant coefficients of an elliptic PDE, -del.a del Phi + b Phi = f, from data on the exterior boundary partial derivative Omega. The assumption made is that the values of the coefficients (a, b) are known cc priori but the information about the geometry of the smooth region boundaries where a and b are discontinuous is missing. For the full characterization of (a, b ) it is then sufficient to find the region boundaries separating different values of the coefficients. This leads to a nonlinear ill-posed inverse pro blem. In this study we propose a numerical algorithm that is based on the f inite-element method and subdivision of the discretization elements. We for mulate the forward problem as a mapping from a set of coefficients represen ting boundary shapes to data on partial derivative Omega, and derive the Ja cobian of this forward mapping. Then an iterative algorithm which seeks a b oundary configuration minimizing the residual norm between measured and pre dicted data is implemented. The method is illustrated first for a general e lliptic PDE and then applied to optical tomography where the goal is to fin d the diffusion and absorption coefficients of the object by transilluminat ing the object with visible or near-infrared light. Numerical test results for this specific application are given with synthetic data.