Weak Galerkin Method for Electrical Impedance Tomography

Ying Liang,Jun Zou

2023-10-02

Abstract:In this work, we propose and analyse a weak Galerkin method for the electrical impedance tomography based on a bounded variation regularization. We use the complete electrode model as the forward system that is approximated by a weak Galerkin method with lowest order. The error estimates are studied for the forward problem, which are used to establish the convergence of this weak Galerkin algorithm for the inverse problem. Numerical examples are presented to verify the effectiveness and efficiency of the weak Galerkin algorithm for the electrical impedance tomography.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the inverse problem in Electrical Impedance Tomography (EIT). EIT is a technique for estimating the physical conductivity inside an object by measuring the voltages on the object's surface, and has attracted wide - spread interest in fields such as non - destructive testing, breast cancer detection, stroke classification, and geophysical exploration. However, EIT is a highly ill - posed and nonlinear problem, especially the reconstruction process is very unstable with respect to boundary measurements. Therefore, developing accurate and flexible forward solvers and using appropriate regularization techniques in the inverse process are crucial for dealing with noisy measurements and complex geometric boundaries associated with the complete electrode model.

To address these challenges, this paper proposes and analyzes a weak Galerkin method based on Bounded Variation (BV) regularization. Specifically, the main contributions of this paper include:

1. **Weak Galerkin method for the forward problem**: Use the lowest - order weak Galerkin formulation to approximate the system and demonstrate the advantages of this method in recovering piecewise - constant conductivity. Derive error estimates for the direct problem and use these estimates to establish the convergence of the inverse - problem algorithm.

2. **BV regularization for the inverse problem**: Adopt a variational approach in the inverse process and seek the reconstruction of physical conductivity by minimizing an appropriate functional that includes BV regularization. BV regularization can preserve the sharp features of the parameters to be determined and is suitable for the recovery of conductivity that is piecewise - regular and has sharp gradients as expected in many applications.

3. **Numerical experiments**: Provide numerical experiment results to verify the effectiveness and efficiency of the weak Galerkin method.

In summary, this paper aims to provide an effective and stable method for solving the EIT inverse problem by combining the weak Galerkin method and BV regularization.