Show simple item record

dc.contributor.advisor Borcea, Liliana
dc.creatorGuevara Vasquez, Fernando
dc.date.accessioned 2018-12-03T18:32:44Z
dc.date.available 2018-12-03T18:32:44Z
dc.date.issued 2006
dc.identifier.urihttps://hdl.handle.net/1911/103693
dc.description.abstract Electric impedance tomography (EIT) consists in finding the conductivity inside a body from electrical measurements taken at its surface. This is a severely ill-posed problem: any numerical inversion scheme requires some form of regularization. We present inversion schemes that address the instability of the problem by proper sparse parametrization of the unknown conductivity. To guide us, we consider a consistent finite difference approach to an inverse Sturm-Liouville problem arising in EIT for layered media. The method first solves a model reduction problem for the differential equation where the reduced model parameters are essentially averages of the conductivity over the cells of a grid depending on the conductivity. Fortunately the dependence is weak. Thus one can efficiently estimate conductivity averages by using the grid for a reference conductivity. This simple inversion method converges to the true solution as the number of measurements increases. We analyze the sensitivity of the reduced model parameters to small changes in the conductivity, and introduce a Newton-type iteration to improve the reconstructions of the simple inversion method. As an added bonus, our method can benefit from a priori information if available. We generalize both methods to the 2D EIT problem by considering finite volumes discretizations of size determined by the measurement precision, but where the node locations are to be determined adaptively. This discretization can be viewed as a resistor network, where the resistors are essentially averages of the conductivity over grid cells. We show that the model reduction problem of finding the smallest resistor network (of fixed topology) that can predict meaningful measurements of the Dirichlet-to-Neumann map is uniquely solvable for a broad class of measurements. We propose a simple inversion method that, as in the simple method for the inverse Sturm-Liouville problem, is based on an interpretation of the resistors as conductivity averages over grid cells, and an iterative method that improves such reconstructions by using sensitivity information on the changes in the resistors due to small changes in the conductivity. A priori information can also be incorporated to the latter method.
dc.format.extent 143 pp
dc.language.iso eng
dc.subjectMathematics
Pure sciences
Electric impedance tomography
Elliptic partial differential equations
dc.title On the parameterization of ill-posed inverse problems arising from elliptic partial differential equations
dc.identifier.digital 305267471
dc.type.genre Thesis
dc.type.material Text
thesis.degree.department Mathematics
thesis.degree.discipline Natural Sciences
thesis.degree.grantor Rice University
thesis.degree.level Doctoral
thesis.degree.name Doctor of Philosophy
dc.identifier.callno THESIS MATH. SCI. 2007 GUEVARA VASQUEZ
dc.identifier.citation Guevara Vasquez, Fernando. "On the parameterization of ill-posed inverse problems arising from elliptic partial differential equations." (2006) Diss., Rice University. https://hdl.handle.net/1911/103693.


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record