Physical conductivities are positive and bounded:
Why the bounds matter.
- Lower bound: it makes the forward problem uniformly elliptic and coercive (Lax-Milgram Theorem). Perfect insulators () would disconnect the domain.
- Upper bound: it keeps the problem uniformly elliptic, gives continuity of the bilinear form, and excludes perfect conductors.
- Regularisation: bounds rule out oscillating minimising sequences and restore compactness in combination with a regulariser (see Kohn-Vogelius Functional and Implicit Regularization).
- Prior knowledge: typical biological tissue lies roughly between and S/m (bone and fat low, blood and cerebrospinal fluid high). The best metallic conductor, silver, is about S/m. After normalising to a reference, simulations often use with a small .
Enforcing them.
- Bound-constrained optimisers: L-BFGS-B, or projected gradient steps with = clipping.
- As the regulariser in ADMM, whose prox is clipping.
- Reparametrisation: , or a sigmoid, turns the problem into an unconstrained one. This changes the geometry of the problem and the gradient ( chain factor).
In ModularEIT.jl: minimize.
References
- L. Borcea (2002). Electrical impedance tomography. Inverse Problems 18(6), R99–R136. doi:10.1088/0266-5611/18/6/201
- C. Gabriel, S. Gabriel, E. Corthout (1996). The dielectric properties of biological tissues: I. Literature survey. Phys. Med. Biol. 41(11), 2231–2249. doi:10.1088/0031-9155/41/11/001
- R. H. Byrd, P. Lu, J. Nocedal, C. Zhu (1995). A Limited Memory Algorithm for Bound Constrained Optimization. SIAM J. Sci. Comput. 16(5), 1190–1208. doi:10.1137/0916069