Abstract:Modeling the behavior of biological tissues and organs often necessitates the knowledge of their shape in the absence of external loads. However, when their geometry is acquired in-vivo through imaging techniques, bodies are typically subject to mechanical deformation due to the presence of external forces, and the load-free configuration needs to be reconstructed. This paper addresses this crucial and frequently overlooked topic, known as the inverse elasticity problem (IEP), by delving into both theoretical and numerical aspects, with a particular focus on cardiac mechanics. In this work, we extend Shield's seminal work to determine the structure of the IEP with arbitrary material inhomogeneities and in the presence of both body and active forces. These aspects are fundamental in computational cardiology, and we show that they may break the variational structure of the inverse problem. In addition, we show that the inverse problem might have no solution even in the presence of constant Neumann boundary conditions and a polyconvex strain energy functional. We then present the results of extensive numerical tests to validate our theoretical framework, and to characterize the computational challenges associated with a direct numerical approximation of the IEP. Specifically, we show that this framework outperforms existing approaches both in terms of robustness and optimality, such as Sellier's iterative procedure, even when the latter is improved with acceleration techniques. A notable discovery is that multigrid preconditioners are, in contrast to standard elasticity, not efficient, where a one-level additive Schwarz and generalized Dryja–Smith–Widlund provide a much more reliable alternative. Finally, we successfully address the IEP for a full-heart geometry, demonstrating that the IEP formulation can compute the stress-free configuration in real-life scenarios where Sellier's algorithm proves inadequate.

On the reconstruction of cavities in a nonlinear model arising from cardiac electrophysiology

On a nonlinear model in domains with cavities arising from cardiac electrophysiology

On a semilinear elliptic boundary value problem arising in cardiac electrophysiology

On an inverse problem with applications in cardiac electrophysiology

On the determination of ischemic regions in the monodomain model of cardiac electrophysiology from boundary measurements

On the inverse problem of detecting cardiac ischemias: theoretical analysis and numerical reconstruction

Reconstructing relaxed configurations in elastic bodies: Mathematical formulations and numerical methods for cardiac modeling

A non-conforming-in-space numerical framework for realistic cardiac electrophysiological outputs

Nonlinear Fourier analysis for discontinuous conductivities: computational results

Reconstruction of Cracks in Calderón's Inverse Conductivity Problem Using Energy Comparisons

A staggered-in-time and non-conforming-in-space numerical framework for realistic cardiac electrophysiology outputs

Reconstructing unknown inclusions for the biharmonic equation

Cardiac Motion and Elasticity Characterization with Iterative Sequential Hinfinity Criteria

Boundary integral formulation of the cell-by-cell model of cardiac electrophysiology

Reconstruction for the Calderón problem with Lipschitz conductivities

Reconstruction of unknown nonlinear operators in semilinear elliptic models using optimal inputs

Reconstruction of Excitation Waves from Mechanical Deformation using Physics-Informed Neural Networks