A viscoactive constitutive modeling framework with variational updates for the myocardium.

We present a constitutive modeling framework for contractile cardiac mechanics by formulating a single variational principle from which incremental stress-strain relations and kinetic rate equations for active contraction and relaxation can all be derived. The variational framework seamlessly incorporates the hyperelastic behavior of the relaxed and contracted tissue along with the rate - and length - dependent generation of contractile force. We describe a three-element, Hill-type model that unifies the active tension and active deformation approaches. As in the latter approach, we multiplicatively decompose the total deformation gradient into active and elastic parts, with the active deformation parametrizing the contractile Hill element. We adopt as internal variables the fiber, cross-fiber, and sheet normal stretch ratios. The kinetics of these internal variables are modeled via definition of a kinetic potential function derived from experimental force-velocity relations. Additionally, we account for dissipation during tissue deformation by adding a Newtonian viscous potential. To model the force activation, the kinetic equations are coupled with the calcium transient obtained from a cardiomyocyte electrophysiology model. We first analyze our model at the material point level using stress and strain versus time curves for different viscosity values. Subsequently, we couple our constitutive framework with the finite element method (FEM) and study the deformation of three-dimensional tissue slabs with varying cardiac myocyte orientation. Finally, we simulate the contraction and relaxation of an ellipsoidal left ventricular model and record common kinematic measures, such as ejection fraction, and myocardial tissue volume changes.

Effect of intra-myocardial Algisyl-LVR™ injectates on fibre structure in porcine heart failure.

Recent preclinical trials have shown that alginate injections are a promising treatment for ischemic heart disease. Although improvements in heart function and global structure have been reported following alginate implants, the underlying structure is poorly understood. Using high resolution ex vivo MRI and DT-MRI of the hearts of normal control swine (n = 8), swine with induced heart failure (n = 5), and swine with heart failure and alginate injection treatment (n = 6), we visualized and quantified the fibre distribution and implant material geometry. Our findings show that the alginate injectates form solid ellipsoids with a retention rate of 68.7% ±â€¯21.3% (mean ±â€¯SD) and a sphericity index of 0.37 ±â€¯0.03. These ellipsoidal shapes solidified predominantly at the mid-wall position with an inclination of -4.9°â€¯±â€¯31.4° relative to the local circumferential direction. Overall, the change to left ventricular wall thickness and myofiber orientation was minor and was associated with heart failure and not the presence of injectates. These results show that alginate injectates conform to the pre-existing tissue structure, likely expanding along directions of least resistance as mass is added to the injection sites. The alginate displaces the myocardial tissue predominantly in the longitudinal direction, causing minimal disruption to the surrounding myofiber orientations.

Regional segmentation of ventricular models to achieve repolarization dispersion in cardiac electrophysiology modeling.

The electrocardiogram (ECG) is one of the most significant outputs of a computational model of cardiac electrophysiology because it relates the numerical results to clinical data and is a universal tool for diagnosing heart diseases. One key features of the ECG is the T-wave, which is caused by longitudinal and transmural heterogeneity of the action potential duration (APD). Thus, in order to model a correct wave of repolarization, different cell properties resulting in different APDs must be assigned across the ventricular wall and longitudinally from apex to base. To achieve this requirement, a regional parametrization of the heart is necessary. We propose a robust approach to obtain the transmural and longitudinal segmentation in a general heart geometry without relying on ad hoc procedures. Our approach is based on auxiliary harmonic lifting analyses, already used in the literature to generate myocardial fiber orientations. Specifically, the solution of a sequence of Laplace boundary value problems allows parametrically controlled segmentation of both heart ventricles. The flexibility and simplicity of the proposed method is demonstrated through several representative examples, varying the locations and extents of the epicardial, midwall, and endocardial layers. Effects of the control parameters on the T-wave morphology are illustrated via computed ECGs.

Passive ventricular mechanics modelling using MRI of structure and function.

Patients suffering from dilated cardiomyopathy or myocardial infarction can develop left ventricular (LV) diastolic impairment. The LV remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions and this remodeling process can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element model was developed to incorporate physiological and mechanical information derived from in vivo magnetic resonance imaging (MRI) tissue tagging, in vivo LV cavity pressure recording and ex vivo diffusion tensor MRI (DTMRI) of a canine heart. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion (the primary eigenvector) in each voxel of a DTMRI directly correlates with the myocardial fibre orientation. This model was customized to the geometry of the canine LV during diastasis by fitting the segmented epicardial and endocardial surface data from tagged MRI using nonlinear finite element fitting techniques. Myofibre orientations, extracted from DTMRI of the same heart, were incorporated into this geometric model using a free form deformation methodology. Pressure recordings, temporally synchronized to the tissue tagging MRI data, were used to simulate the LV deformation during diastole. Simulation of the diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. This integrated physiological model will allow more insight into the regional passive diastolic mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction in pathological conditions.