A Caltech Library Service

Discrete Geometric Structures in Homogenization and Inverse Homogenization with Application to EIT

Desbrun, Mathieu and Donaldson, Roger D. and Owhadi, Houman (2009) Discrete Geometric Structures in Homogenization and Inverse Homogenization with Application to EIT. ACM Technical Reports, 2009-02. California Institute of Technology , Pasadena, CA. (Unpublished)

PDF - Submitted Version
See Usage Policy.


Use this Persistent URL to link to this item:


We introduce a new geometric approach for the homogenization and inverse homogenization of the divergence form elliptic operator with rough conductivity coefficients σ(x) in dimension two. We show that conductivity coefficients are in one-to-one correspondence with divergence-free matrices and convex functions s(x) over the domain Ω. Although homogenization is a non-linear and non-injective operator when applied directly to conductivity coefficients, homogenization becomes a linear interpolation operator over triangulations of Ω when re-expressed using convex functions, and is a volume averaging operator when re-expressed with divergence-free matrices. We explicitly give the transformations which map conductivity coefficients into divergence-free matrices and convex functions, as well as their respective inverses. Using optimal weighted Delaunay triangulations for linearly interpolating convex functions, we apply this geometric framework to obtain an optimally robust homogenization algorithm for arbitrary rough coefficients, extending the global optimality of Delaunay triangulations with respect to a discrete Dirichlet energy to weighted Delaunay triangulations. Next, we consider inverse homogenization, that is, the recovery of the microstructure from macroscopic information, a problem which is known to be both non-linear and severly ill-posed. We show how to decompose this reconstruction into a linear ill-posed problem and a well-posed non-linear problem. We apply this new geometric approach to Electrical Impedance Tomography (EIT) in dimension two. It is known that the EIT problem admits at most one isotropic solution. If an isotropic solution exists, we show how to compute it from any conductivity having the same boundary Dirichlet-to-Neumann map. This is of practical importance since the EIT problem always admits a unique solution in the space of divergence-free matrices and is stable with respect to G-convergence in that space (this property fails for isotropic matrices). As such, we suggest that the space of convex functions is the natural space to use to parameterize solutions of the EIT problem.

Item Type:Report or Paper (Technical Report)
Related URLs:
URLURL TypeDescription Paper
Desbrun, Mathieu0000-0003-3424-6079
Owhadi, Houman0000-0002-5677-1600
Group:Applied & Computational Mathematics
Series Name:ACM Technical Reports
Issue or Number:2009-02
Record Number:CaltechAUTHORS:20111011-163848887
Persistent URL:
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:27172
Deposited On:19 Oct 2011 18:16
Last Modified:03 Oct 2019 03:21

Repository Staff Only: item control page