Lipschitz stability estimate and reconstruction of Lam\'e parameters in linear elasticity

TLDR: In this article, the authors consider the inverse problem of recovering an isotropic elastic tensor from the Neumann-to-Dirichlet map and prove a Lipschitz stability estimate for Lame parameters with certain regularity assumptions.

Abstract: In this paper, we consider the inverse problem of recovering an isotropic elastic tensor from the Neumann-to-Dirichlet map. To this end, we prove a Lipschitz stability estimate for Lame parameters with certain regularity assumptions. In addition, we assume that the Lame parameters belong to a known finite subspace with a priori known bounds and that they fulfill a monotonicity property. The proof relies on a monotonicity result combined with the techniques of localized potentials. To numerically solve the inverse problem, we propose a Kohn-Vogelius-type cost functional over a class of admissible parameters subject to two boundary value problems. The reformulation of the minimization problem via the Neumann-to-Dirichlet operator allows us to obtain the optimality conditions by using the Frechet differentiability of this operator and its inverse. The reconstruction is then performed by means of an iterative algorithm based on a quasi-Newton method. Finally, we give and discuss several numerical examples.