A transversely isotropic thermo-hyperelastic constitutive model of myocardial tissue with a three-state cell death dynamics for cardiac radiofrequency ablation

• Published in: Journal of the mechanics and physics of solids 2022-04, Vol.161, p.104810, Article 104810
• Main Authors: Molinari, Leonardo, Gerardo-Giorda, Luca, Gizzi, Alessio
• Format: Article
• Language: English
• Subjects: Ablation; Apoptosis; Cell death; Cellular death dynamics; Constitutive model; Constitutive models; Constraint modelling; Couplings; Design optimization; Dynamic models; Finite element analysis; Multiphysics couplings; Myocardium; Myocardium anisotropy; Power control; Radio frequency; Radiofrequency catheter ablation; Thermoelasticity
• ISSN: 0022-5096; 1873-4782
• DOI: 10.1016/j.jmps.2022.104810

Constitutive models of radiofrequency catheter ablation (RFCA) lack a reliable three-dimensional microstructural representation of the myocardium. Besides, complex multiphysics couplings arise during tissue heating. In light of procedure design and optimization, we propose a generalized, thermodynamically consistent, transverse isotropic thermo-hyperelastic constitutive model of the myocardium accounting for local anisotropies and multiscale dynamics. We advance, in particular, a biophysical rationale formulating a continuum damage approach based on a three-state hyperthermic cell death dynamical model. We further introduce a fully coupled thermo-mechanical model by means of the multiplicative decomposition of the deformation gradient and by defining phenomenological temperature dependencies of the material parameters consistent with the general theory of thermoelasticity. The overall multiphysics and multiscale model is numerically solved within an idealized tissue domain, adopting an accurate finite element scheme under the additional constraint of constant power control. Our modeling approach fills the gap of constitutive modeling of cardiac RFCA and results in a better matching of ablating volumes with respect to the current models in the literature. We show for the first time the elliptical shape of the lesion occurring because of anisotropic thermo-electric conduction and the occurrence of residual strains due to tissue phase changes after ablation. We conclude by discussing limitations, potential future applications and perspectives of our modeling approach towards generalized constitutive remodeling theories.