Abstract:We study the numerical reconstruction problem in acousto-electric tomography (AET) of recovering the conductivity distribution in a bounded domain from multiple interior power density data. The Two-Point-Gradient-$\Theta$ (TPG-$\Theta$) in Kaczmarz type is proposed, with a general convex penalty term $\Theta$, the algorithm can be utilized in AET problem for recovering sparse and discontinuous conductivity distributions. We establish the convergence of such iterative regularized method. Extensive numerical experiments are presented to illustrate the feasibility and effectiveness of the proposed approach.
What problem does this paper attempt to address?
This paper aims to solve the numerical reconstruction problem in acousto - electric tomography (AET), specifically to recover the conductivity distribution within a bounded domain from multiple internal power density data. AET is a hybrid imaging technique that combines electrical impedance tomography (EIT) and ultrasonic modulation, thereby providing additional internal information and significantly improving the quality of conductivity reconstruction.
### Specific Problem Description
1. **Background and Requirements**
- Conductivity values vary greatly among different soft tissue types, and accurate imaging of it can provide valuable information about the physiological and pathological states of tissues.
- As an emerging technology, EIT aims to reconstruct the conductivity distribution in the body from voltage and current measurements on the body surface. However, EIT is a severely ill - posed problem, so convergent and stable algorithms need to be proposed to reconstruct conductivity.
2. **Limitations of Existing Methods**
- Traditional EIT methods are not effective in dealing with complex conductivity distributions, especially in cases where the conductivity distribution is sparse or discontinuous.
- EIT methods combined with other physical phenomena (such as magnetic resonance, ultrasonic waves) (such as magnetic resonance EIT, acousto - electric tomography, etc.) have improved, but still face challenges.
3. **Core Contributions of the Paper**
- A new Kaczmarz - type iterative regularization method - the Two - Point - Gradient - Θ (TPG - Θ) method is proposed to stably approximate the conductivity distribution in AET.
- By introducing a general convex penalty term Θ, this method can effectively recover sparse and discontinuous conductivity distributions in the AET problem.
- The convergence of this iterative regularization method is established, and its feasibility and effectiveness are verified through extensive numerical experiments.
### Mathematical Model and Method
The mathematical formulation of the AET problem is as follows:
Given a bounded domain \(\Omega\subset\mathbb{R}^d\) and conductivity \(\sigma > 0\), recover the conductivity distribution from \(I\) internal power density measurements \(H_i(\sigma)=\sigma|\nabla u_i(\sigma)|^2\). Where \(u_i(\sigma)\) is the weak solution of the following elliptic equation:
\[
\begin{cases}
-\nabla\cdot(\sigma\nabla u_i) = 0, & \text{in }\Omega,\\
\sigma\frac{\partial u_i}{\partial\nu}=f_i, & \text{on }\partial\Omega,
\end{cases}
\]
where \(\nu\) represents the unit outer normal vector on \(\partial\Omega\).
### Proposed Method
The TPG - Θ method proposed in the paper is based on the Kaczmarz - type iterative regularization framework, combined with the two - and - a - half - gradient acceleration technique. The specific steps are as follows:
1. **Initialization**
- Set the initial guesses \(\xi_0\) and \(\sigma_0\).
2. **Iterative Update**
- For each \(i = 0,1,\ldots,I - 1\), calculate the combination parameter \(\lambda_{n,i}\), and define
\[
\zeta_{n,i}=\xi_{n,i}+\lambda_{n,i}(\xi_{n,i}-\xi_{n,i - 1}),
\]
\[
z_{n,i}=\arg\min_{\sigma\in L^2(\Omega)}\left\{\Theta(\sigma)-\langl
\]