Finite element approximation for quantitative photoacoustic tomography in a diffusive regime
By: Giovanni S. Alberti, Siyu Cen, Zhi Zhou
Potential Business Impact:
Find hidden things inside bodies using light and sound.
In this paper, we focus on the numerical analysis of quantitative photoacoustic tomography. Our goal is to reconstruct the optical coefficients, i.e., the diffusion and absorption coefficients, using multiple internal observational data. The foundation of our numerical algorithm lies in solving an inverse diffusivity problem and a direct problem associated with elliptic equations. The stability of the inverse problem depends critically on a non-zero condition in the internal observations, a condition that can be met using randomly chosen boundary excitation data. Utilizing these randomly generated boundary data, we implement an output least squares formulation combined with finite element discretization to solve the inverse problem. In this scenario, we provide a rigorous error estimate in $L^2(\Omega)$ norm for the numerical reconstruction using a weighted energy estimate, inspired by the analysis of a newly proposed conditional stability result. The resulting error estimate serves as a valuable guide for selecting appropriate regularization parameters and discretization mesh sizes according to the noise levels present in the data. Several numerical experiments are presented to support our theoretical results and illustrate the effectiveness of our numerical scheme.
Similar Papers
A Fourier finite volume approach for the optical inverse problem of quantitative photoacoustic tomography
Computational Physics
Helps doctors see inside bodies better.
A shape-optimization approach for inverse diffusion problems using a single boundary measurement
Numerical Analysis
Find hidden shapes from just one measurement.
Numerical Analysis of Simultaneous Reconstruction of Initial Condition and Potential in Subdiffusion
Numerical Analysis
Find hidden information in slow-moving processes.