Score: 1

Finite element approximation for quantitative photoacoustic tomography in a diffusive regime

Published: May 8, 2025 | arXiv ID: 2505.05361v1

By: Giovanni S. Alberti, Siyu Cen, Zhi Zhou

Potential Business Impact:

Find hidden things inside bodies using light and sound.

Business Areas:
Laser Hardware, Science and Engineering

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.

Country of Origin
🇭🇰 🇮🇹 Hong Kong, Italy

Page Count
24 pages

Category
Mathematics:
Numerical Analysis (Math)