Stability analysis of an inverse coefficients problem in a system of partial differential equations
Houcine Meftahi, Chayma Nssibi
TL;DR
The paper tackles the inverse coefficients problem in linear elasticity by aiming to recover $\lambda(x)$, $\mu(x)$ and $\rho(x)$ from the Neumann-to-Dirichlet map $\Lambda_{\lambda,\mu,\rho}$ on $d\ge 2$ dimensions. It first treats the density reconstruction with known Lamé parameters, proving Lipschitz stability via a monotonicity relation and the Runge approximation, and provides a constructive Lipschitz estimate for piecewise-constant densities. It then extends to simultaneous recovery of all coefficients under a priori bounds and finite-dimensional subspace constraints, employing monotonicity and localized potentials to obtain Lipschitz stability, including a constructive scheme. The results enable robust numerical reconstruction and provide theoretical guarantees for finite-dimensional parameter spaces, with potential impact on non-destructive tissue assessment and related elastic imaging problems. Overall, the work advances the understanding of stability in multiparameter inverse elasticity problems beyond scalar conductivity analogies by leveraging monotonicity and localized potentials in arbitrary dimensions.
Abstract
In this study, we address the inverse problem of recovering the Lamé parameters ($λ, μ$) and the density $ρ$ of a medium from the Neumann-to-Dirichlet map for any dimension $d\geq 2$. This inverse problem finds its motivation in the reconstruction of mechanical properties of tissues in medical diagnostics. We first assume that the Lamé parameters ($λ, μ$) are know and we look for the inverse problem of recovering the density $ρ$. In this context, we derive a constrcutive Lipschitz stability estimate in terms of the Neumann to Dirichlet map in the case of piecewise constant parameters. Then, we look for the inverse problem of recovering $λ$, $μ$ and $ρ$ simultameousely. We establish Lipschitz stability estimate, provided that the parameters $λ$, $μ$ and $ρ$ have upper and lower bounds and belong to a known finite-dimensional subspace. The proofs hinge on monotonicity relations between the parameters and the Neumann-to-Dirichlet operator, coupled with the techniques of localized potentials.
