A mixed finite element method for the forward problem of electrical impedance tomography with the shunt model

We propose a mixed finite element method for the forward problem of electrical impedance tomography with the shunt model. The method is based on a constrained minimization formulation in terms of the current density, where parts of the Lagrange multiplier are interpreted as the voltages on the electrodes. We justify the well-posedness of this continuous formulation, estimate the regularity of its solution, and use Raviart–Thomas elements for its discretization. The resulting mixed method is then proved to be stable, and converge especially with respect to the voltages on the electrodes. Simulation on a 2D model also gives consistent results.