Computational analysis of the effect of valvular regurgitation on ventricular mechanics using a 3D electromechanics model

In silico screening method for non-responders to cardiac resynchronization therapy in patients with heart failure: a pilot study

Background

Cardiac resynchronization therapy (CRT) is an effective treatment option for patients with heart failure (HF) and left ventricular (LV) dyssynchrony. However, the problem of some patients not responding to CRT remains unresolved. This study aimed to propose a novel in silico method for CRT simulation.

Methods

Three-dimensional heart geometry was constructed from computed tomography images. The finite element method was used to elucidate the electric wave propagation in the heart. The electric excitation and mechanical contraction were coupled with vascular hemodynamics by the lumped parameter model. The model parameters for three-dimensional (3D) heart and vascular mechanics were estimated by matching computed variables with measured physiological parameters. CRT effects were simulated in a patient with HF and left bundle branch block (LBBB). LV end-diastolic (LVEDV) and end-systolic volumes (LVESV), LV ejection fraction (LVEF), and CRT responsiveness measured from the in silico simulation model were compared with those from clinical observation. A CRT responder was defined as absolute increase in LVEF ≥ 5% or relative increase in LVEF ≥ 15%.

Results

A 68-year-old female with nonischemic HF and LBBB was retrospectively included. The in silico CRT simulation modeling revealed that changes in LVEDV, LVESV, and LVEF by CRT were from 174 to 173 mL, 116 to 104 mL, and 33 to 40%, respectively. Absolute and relative ΔLVEF were 7% and 18%, respectively, signifying a CRT responder. In clinical observation, echocardiography showed that changes in LVEDV, LVESV, and LVEF by CRT were from 162 to 119 mL, 114 to 69 mL, and 29 to 42%, respectively. Absolute and relative ΔLVESV were 13% and 31%, respectively, also signifying a CRT responder. CRT responsiveness from the in silico CRT simulation model was concordant with that in the clinical observation.

Conclusion

This in silico CRT simulation method is a feasible technique to screen for CRT non-responders in patients with HF and LBBB.

Transmural Distribution of Coronary Perfusion and Myocardial Work Density Due to Alterations in Ventricular Loading, Geometry and Contractility

Myocardial supply changes to accommodate the variation of myocardial demand across the heart wall to maintain normal cardiac function. A computational framework that couples the systemic circulation of a left ventricular (LV) finite element model and coronary perfusion in a closed loop is developed to investigate the transmural distribution of the myocardial demand (work density) and supply (perfusion) ratio. Calibrated and validated against measurements of LV mechanics and coronary perfusion, the model is applied to investigate changes in the transmural distribution of passive coronary perfusion, myocardial work density, and their ratio in response to changes in LV contractility, preload, afterload, wall thickness, and cavity volume. The model predicts the following: (1) Total passive coronary flow varies from a minimum value at the endocardium to a maximum value at the epicardium transmurally that is consistent with the transmural distribution of IMP; (2) Total passive coronary flow at different transmural locations is increased with an increase in either contractility, afterload, or preload of the LV, whereas is reduced with an increase in wall thickness or cavity volume; (3) Myocardial work density at different transmural locations is increased transmurally with an increase in either contractility, afterload, preload or cavity volume of the LV, but is reduced with an increase in wall thickness; (4) Myocardial work density-perfusion mismatch ratio at different transmural locations is increased with an increase in contractility, preload, wall thickness or cavity volume of the LV, and the ratio is higher at the endocardium than the epicardium. These results suggest that an increase in either contractility, preload, wall thickness, or cavity volume of the LV can increase the vulnerability of the subendocardial region to ischemia.

Computational analysis of the effect of KCNH2 L532P mutation on ventricular electromechanical behaviors

The KCNH2 L532P mutation is an alteration in the I_Kr channel that is associated with short QT syndrome and atrial fibrillation in zebrafish. In preliminary studies, the electrophysiological effects of the hERG L532P mutation were investigated using a mathematical model in a single-cell and 2D sheet medium. The objective of this study was to quantify the effects of the KCNH2 L532P mutation on the 3D ventricular electrophysiological behavior and the mechanical pumping responses. We used a realistic three-dimensional ventricular electrophysiological-mechanical model, which was adjusted into two conditions: the wild-type (WT) condition, i.e., the original case of the Tusscher et al. model, and the L532P mutation condition, with modification of the original I_Kr equation. The action potential duration (APD) in the mutant ventricle was reduced by 73% owing to the significant increase of the I_Kr current density. In the 3D simulation, the L532P mutation maintained the sustainability of reentrant waves; however, the reentry was terminated in the WT condition. The contractility of the ventricle with L532P mutation was significantly reduced compared with that in WT which results in sustain shivering heart during reentry condition. The reduction of the contractility was associated with the shortening APD which simultaneously shortened the duration of the Ca²⁺ channel opening. In conclusion, the ventricle with KCNH2 L532P mutation is prone to reentry generation with a sustained chaotic condition, and the mutation significantly reduced the pumping performance of the ventricles.

Prediction of Cardiac Mechanical Performance From Electrical Features During Ventricular Tachyarrhythmia Simulation Using Machine Learning Algorithms

In ventricular tachyarrhythmia, electrical instability features including action potential duration, dominant frequency, phase singularity, and filaments are associated with mechanical contractility. However, there are insufficient studies on estimated mechanical contractility based on electrical features during ventricular tachyarrhythmia using a stochastic model. In this study, we predicted cardiac mechanical performance from features of electrical instability during ventricular tachyarrhythmia simulation using machine learning algorithms, including support vector regression (SVR) and artificial neural network (ANN) models. We performed an electromechanical tachyarrhythmia simulation and extracted 12 electrical instability features and two mechanical properties, including stroke volume and the amplitude of myocardial tension (ampTens). We compared predictive performance according to kernel types of the SVR model and the number of hidden layers of the ANN model. In the SVR model, the prediction accuracies of stroke volume and ampTens were the highest when using the polynomial kernel and linear kernel, respectively. The predictive performance of the ANN model was better than that of the SVR model. The prediction accuracies were the highest when the ANN model consisted of three hidden layers. Accordingly, we propose the ANN model with three hidden layers as an optimal model for predicting cardiac mechanical contractility in ventricular tachyarrhythmia. The results of this study are expected to be used to indirectly estimate the hemodynamic response from the electrical cardiac map measured by the optical mapping system during cardiac surgery, as well as cardiac contractility under normal sinus rhythm conditions.

Effect of myocardial heterogeneity on ventricular electro-mechanical responses: a computational study

BACKGROUND

The heart wall exhibits three layers of different thicknesses: the outer epicardium, mid-myocardium, and inner endocardium. Among these layers, the mid-myocardium is typically the thickest. As indicated by preliminary studies, heart-wall layers exhibit various characteristics with regard to electrophysiology, pharmacology, and pathology. Construction of an accurate three-dimensional (3D) model of the heart is important for predicting physiological behaviors. However, the wide variability of myocardial shapes and the unclear edges between the epicardium and soft tissues are major challenges in the 3D model segmentation approach for identifying the boundaries of the epicardium, mid-myocardium, and endocardium. Therefore, this results in possible variations in the heterogeneity ratios between the epicardium, mid-myocardium, and endocardium. The objective of this study was to observe the effects of different thickness ratios of the epicardium, mid-myocardium, and endocardium on cardiac arrhythmogenesis, reentry instability, and mechanical responses during arrhythmia.

METHODS

We used a computational method and simulated three heterogeneous ventricular models: Model 1 had the thickest M cell layer and thinnest epicardium and endocardium. Model 2 had intermediate layer thicknesses. Model 3 exhibited the thinnest mid-myocardium and thickest epicardium and endocardium. Electrical and mechanical simulations of the three heterogeneous models were performed under normal sinus rhythm and reentry conditions.

RESULTS

Model 1 exhibited the highest probability of terminating reentrant waves, and Model 3 exhibited to experience greater cardiac arrhythmia. In the reentry simulation, at 8 s, Model 3 generated the largest number of rotors (eight), while Models 1 and 2 produced five and seven rotors, respectively. There was no significant difference in the cardiac output obtained during the sinus rhythm. Under the reentry condition, the highest cardiac output was generated by Model 1 (19 mL/s), followed by Model 2 (9 mL/s) and Model 3 (7 mL/s).

CONCLUSIONS

A thicker mid-myocardium led to improvements in the pumping efficacy and contractility and reduced the probability of cardiac arrhythmia. Conversely, thinner M cell layers generated more unstable reentrant spiral waves and hindered the ventricular pumping.

Intra-aortic balloon pump postcardiac surgery: A literature review

Intra-aortic balloon pump (IABP) has been the most commonly used mechanical assist circulatory device in many postcardiotomy low output disorders for decades. Mechanism of IABP is based on its inflation in time of the diastolic pressure in the aortic root resulting increase in the blood and oxygen amount of the coronary artery and its deflation in left ventricular afterload during the systolic period. Prophylactic and postoperative application of IABP has been suggested by researchers, which has been commonly used in high-risk patients undertaking coronary artery bypass grafting surgery or percutaneous coronary intervention. Other researchers put forward the idea of the percutaneous IABP insertion throughout the left axillary artery as a reliable and relatively well-tolerated approach and also as a recovery tool to bridge patients with end-stage heart failure to heart transplantation. The current review was aimed to give further insight into routine IABP application by presenting the basic principles and trends in the incidence, management, role of IABP recovery, and long-lasting mortality outcomes in patients with cardiovascular disorders and discussing previous and current evidence.

Computational analysis of the electromechanical performance of mitral valve cerclage annuloplasty using a patient-specific ventricular model

We aimed to propose a novel computational approach to predict the electromechanical performance of pre- and post-mitral valve cerclage annuloplasty (MVCA). Furthermore, we tested a virtual estimation method to optimize the left ventricular basement tightening scheme using a pre-MVCA computer model. The present model combines the three-dimensional (3D) electromechanics of the ventricles with the vascular hemodynamics implemented in a lumped parameter model. 3D models of pre- and post-MVCA were reconstructed from the computed tomography (CT) images of two patients and simulated by solving the electromechanical-governing equations with the finite element method. Computed results indicate that reduction of the dilated heart chambers volume (reverse remodeling) appears to be dependent on ventricular stress distribution. Reduced ventricular stresses in the basement after MVCA treatment were observed in the patients who showed reverse remodeling of heart during follow up over 6 months. In the case who failed to show reverse remodeling after MVCA, more virtual tightening of the ventricular basement diameter than the actual model can induce stress unloading, aiding in heart recovery. The simulation result that virtual tightening of the ventricular basement resulted in a marked increase of myocardial stress unloading provides in silico evidence for a functional impact of MVCA treatment on cardiac mechanics and post-operative heart recovery. This technique contributes to establishing a pre-operative virtual rehearsal procedure before MVCA treatment by using patient-specific cardiac electromechanical modeling of pre-MVCA.

Influence of the KCNQ1 S140G Mutation on Human Ventricular Arrhythmogenesis and Pumping Performance: Simulation Study

The KCNQ1 S140G mutation, which is involved in I_Ks current, affects atrial fibrillation. However, little is known about its effect on the mechanical behavior of the heart. Therefore, we assessed the influence of the KCNQ1 S140G mutation on ventricular electrophysiological stability and mechanical pumping performance using a multi-scale model of cardiac electromechanics. An image-based electromechanical model was used to assess the effect on electrical propagation and arrhythmogenesis of the KCNQ1 S140G mutation. In addition, it was used to compare the mechanical response under the wild-type (WT) and S140G mutation conditions. The intracellular calcium transient obtained from the electrophysiological model was applied as an input parameter to a mechanical model to implement excitation–contraction coupling. The I_Ks current equation was modified to account for expression of the KCNQ1 S140G mutation, and it included a scaling factor (ϕ) for mutant expressivity. The WT and S140G mutation conditions were compared at the single-cell and three-dimensional (3D) tissue levels. The action potential duration (APD) was reduced by 60% by the augmented I_Ks current under the S140G mutation condition, which resulted in shorter QT interval. This reduced the 3D sinus rhythm wavelength by 60% and the sustained re-entry by 56%. However, pumping efficiency of mutant ventricles was superior in sinus rhythm condition. In addition, the shortened wavelength in cardiac tissue allowed a re-entrant circuit to form and increased the probability of sustaining ventricular tachycardia and ventricular fibrillation. In contrast, under the WT condition, a normal wavelength (20.8 cm) was unlikely to initiate and sustain re-entry in the cardiac tissue. Subsequently, the S140G mutant ventricles developed a higher dominant frequency distribution range (2.0–5.3 Hz) than the WT condition (2.8–3.7 Hz). In addition, stroke volume of mutant ventricles was reduced by 65% in sustained re-entry compared to the WT condition. In conclusion, signs of the S140G mutation might be difficult to identify in sinus rhythm even though the mutant ventricles show shortened QT interval. This suggests that the KCNQ1 S140G mutation increases the risk of death by sudden cardiac arrest. In addition, the KCNQ1 S140G mutation can induce ventricular arrhythmia and lessen ventricular contractility under re-entrant conditions.

The effect of myocardial action potential duration on cardiac pumping efficacy: a computational study

BACKGROUND AND AIMS

Although studies on the relation between arrhythmias and the action potential duration (APD) have been carried out, most of them are based only on electrophysiological factors of the heart and lack experiments that consider cardiac mechanical and electromechanical characteristics. Therefore, we conducted this study to clarify the relevance of the shortening of APD of a cell in relation to the mechanical contraction activity of the heart and the associated risk of arrhythmia.

METHODS

The human ventricular model used in this study has two dynamic characteristics: electrophysiological conduction and mechanical contraction. The model simulating electrophysiological characteristics was consisted of lumped parameter circuit that can mimic the phenomenon of ion exchange through the cell membrane of myocyte and consisted of 214,319 tetrahedral finite elements. In contrast, the model simulating mechanical contraction characteristics was constructed to mimic cardiac contraction by means of the crossbridge of a myofilament and consisted of 14,720 hermite-based finite elements to represent a natural 3D curve of the cardiac surface. First, we performed a single cell simulation and the electrophysiological simulation according to the change of the APD by changing the electrical conductivity of the I Ks channel. Thus, we confirmed the correlation between APD and intracellular Ca2+ concentration. Then, we compared mechanical response through mechanical simulation using Ca2+ data from electrical simulation.

RESULTS

The APD and the sum of the intracellular Ca2+ concentrations showed a positive correlation. The shortened APD reduced the conduction wavelength of ventricular cells by shortening the plateau and early repolarization in myocardial cells. The decrease in APD reduced ventricular pumping efficiency by more than 60% as compared with the normal group (normal conditions). This change is caused by the decline of ventricular output owing to reduced ATP consumption during the crossbridge of myofilaments and decreased tension.

CONCLUSION

The shortening of APD owing to increased electrical conductivity of a protein channel on myocardial cells likely decreases the wavelength and the pumping efficiency of the ventricles. Additionally, it may increase tissue sensitivity to ventricular fibrillation, including reentry, and cause symptoms such as dyspnea and dizziness.

Effect of KCNQ1 G229D mutation on cardiac pumping efficacy and reentrant dynamics in ventricles: Computational study

There is growing interest in genetic arrhythmia since mutations in gene which encodes the ion channel underlie numerous arrhythmias. Hasegawa et al reported that G229D mutation in KCNQ1 underlies atrial fibrillation due to significant shortening of action potential duration (APD) in atrial cells. Here, we predicted whether KCNQ1 G229D mutation affects ventricular fibrillation generation, although it shortens APD slightly compared with the atrial cell. We analyzed the effects of G229D mutation on electrical and mechanical ventricle behaviors (not considered in previous studies). We compared action potential shapes under wild-type and mutant conditions. Electrical wave propagations through ventricles were analyzed during sinus rhythm and reentrant conditions. I_Ks enhancement due to G229D mutation shortened the APD in the ventricular cells (6%, 0.3%, and 8% for endo, M, and epi-cells, respectively). The shortened APD contributed to 7% shortening of QT intervals, 29% shortening of wavelengths, 20% decrease in intraventricular pressure, and increase in end-systolic volume 17%, end-diastolic volume 7%, and end-diastolic pressure 11%, which further resulted in reduction in stroke volume as well as cardiac output (28%), ejection fraction 33% stroke work 44%, and ATP consumption 28%. In short, using computational model of the ventricle, we predicted that G229D mutation decreased cardiac pumping efficacy and increased the vulnerability of ventricular fibrillation.

Computational analysis of the effect of mitral and aortic regurgitation on the function of ventricular assist devices using 3D cardiac electromechanical model

Valvular insufficiency affects cardiac responses and the pumping efficacy of left ventricular assist devices (LVADs) when patients undergo LVAD therapy. Knowledge of the effect of valvular regurgitation on the function of LVADs is important when treating heart failure patients. The goal of this study was to examine the effect of valvular regurgitation on the ventricular mechanics of a heart under LVAD treatment and the pumping efficacy of an LVAD using a computational model of the cardiovascular system. For this purpose, a 3D electromechanical model of failing ventricles in a human heart was coupled with a lumped-parameter model of valvular regurgitation and an LVAD-implanted vascular system. We used the computational model to predict cardiac responses with respect to the severity of valvular regurgitation in the presence of LVAD treatment. An LVAD could reduce left ventricle (LV) pressure (up to 34%) and end-diastolic ventricular volume (up to 80%) and maintain cardiac output at the estimated flow rate from the LVAD under the condition of mitral regurgitation (MR); however, the opposite would occur under the condition of aortic regurgitation (AR). Considering these physiological responses, we conclude that AR, and not MR, diminishes the pumping function of LVADs.

Computational prediction of the effects of the intra-aortic balloon pump on heart failure with valvular regurgitation using a 3D cardiac electromechanical model

Intra-aortic balloon pump (IABP) is normally contraindicated in significant aortic regurgitation (AR). It causes and aggravates pre-existing AR while performing well in the event of mitral regurgitation (MR). Indirect parameters, such as the mean systolic pressure, product of heart rate and peak systolic pressure, and pressure–volume are used to quantify the effect of IABP on ventricular workload. However, to date, no studies have directly quantified the reduction in workload with IABP. The goal of this study is to examine the effect of IABP therapy on ventricular mechanics under valvular insufficiency by using a computational model of the heart. For this purpose, the 3D electromechanical model of the failing ventricles used in previous studies was coupled with a lumped parameter model of valvular regurgitation and the IABP-treated vascular system. The IABP therapy was disturbed in terms of reducing the myocardial tension generation and contractile ATP consumption by valvular regurgitation, particularly in the AR condition. The IABP worsened the problem of ventricular expansion induced as a result of the regurgitated blood volume during the diastole under the AR condition. The IABP reduced the LV stroke work in the AR, MR, and no regurgitation conditions. Therefore, the IABP helped the ventricle to pump blood and reduced the ventricular workload. In conclusion, the IABP partially performed its role in the MR condition. However, it was disturbed by the AR and worsened the cardiovascular responses that followed the AR. Therefore, this study computationally proved the reason for the clinical contraindication of IABP in AR patients.

The effect of electrical conductivity of myocardium on cardiac pumping efficacy: a computational study

Background and aims

The existence of non-excitable cells in the myocardium leads to the increasing conduction non-uniformity and decreasing myocardial electrical conductivity. Slowed myocardial conduction velocity (MCV) believed to enhance the probability of cardiac arryhthmia and alter the cardiac mechanical pumping efficacy, even in sinus rhythm. Though several studies on the correlation between MCV and cardiac electrical instabilities exist, there has been no study concerning correlation or causality between MCV and cardiac mechanical pumping efficacy, due to the limitation in clinical methods to document and evaluate cardiac mechanical responses directly. The goal of this study was to examine quantitatively the cardiac pumping efficacy under various MCV conditions using three-dimensional (3D) electromechanical model of canine’s failing ventricle.

Methods

The electromechanical model used in this study composed of the electrical model coupled with the mechanical contraction model along with a lumped model of the circulatory system. The electrical model consisted of 241,725 nodes and 1,298,751 elements of tetrahedral mesh, whereas the mechanical model consisted of 356 nodes and 172 elements of hexahedral mesh with Hermite basis. First, we performed the electrical simulation for five different MCV conditions, from 30 to 70 cm/s with 10 cm/s interval during sinus pacing. Then, we compared the cardiac electrical and mechanical responses of each MCV condition, such as the electrical activation time (EAT), pressure, volume, and energy consumption of the myocardium. The energy consumption of the myocardium was calculated by integrating ATP consumption rate of each node in myofilament model.

Results

The result showed that under higher MCV conditions, the EAT, energy consumption, end diastolic and systolic volume are gradually decreased. Meanwhile, the systolic pressure, stroke volume, stroke work, and stroke work to ATP are increased as the MCV values increased. The cardiac functions and performances are more efficient under higher MCV conditions by consuming smaller energy (ATP) while carrying more works.

Conclusion

In conclusion, this study reveals that MCV has strong correlation with the cardiac pumping efficacy. The obtained results provide useful information to estimate the effect of MCV on the electro-physiology and hemodynamic responses of the ventricle and can be used for further study about arrhythmogeneis and heart failure.

Mechanisms of astrocytic K(+) clearance and swelling under high extracellular K(+) concentrations

In response to the elevation of extracellular K(+) concentration ([K(+)]out), astrocytes clear excessive K(+) to maintain conditions necessary for neural activity. K(+) clearance in astrocytes occurs via two processes: K(+) uptake and K(+) spatial buffering. High [K(+)]out also induces swelling in astrocytes, leading to edema and cell death in the brain. Despite the importance of astrocytic K(+) clearance and swelling, the underlying mechanisms remain unclear. Here, we report results from a simulation analysis of astrocytic K(+) clearance and swelling. Astrocyte models were constructed by incorporating various mechanisms such as intra/extracellular ion concentrations of Na(+), K(+), and Cl(-), cell volume, and models of Na,K-ATPase, Na-K-Cl cotransporter (NKCC), K-Cl cotransporter, inwardly-rectifying K(+) (KIR) channel, passive Cl(-) current, and aquaporin channel. The simulated response of astrocyte models under the uniform distribution of high [K(+)]out revealed significant contributions of NKCC and Na,K-ATPase to increases of intracellular K(+) and Cl(-) concentrations, and swelling. Moreover, we found that, under the non-uniform distribution of high [K(+)]out, KIR channels localized at synaptic clefts absorbed excess K(+) by depolarizing the equivalent potential of K(+) (E K) above membrane potential, while K(+) released through perivascular KIR channels was enhanced by hyperpolarizing E K and depolarizing membrane potential. Further analysis of simulated drug effects revealed that astrocyte swelling was modulated by blocking each of the ion channels and transporters. Our simulation analysis revealed controversial mechanisms of astrocytic K(+) clearance and swelling resulting from complex interactions among ion channels and transporters.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

Mathematical analysis of the effects of valvular regurgitation on the pumping efficacy of continuous and pulsatile left ventricular assist devices

•

We numerically investigated the physiological relationship between the severity of regurgitation and the effect of a left ventricular assist device (LVAD) on cardiovascular system responses.

•

Under conditions of mitral regurgitation, the effects of both pulsatile and continuous LVAD treatment on ventricular unloading were significant.

•

Under conditions of aortic regurgitation (AR), the effects of the LVADs on ventricular unloading were not significant. The effects of LVAD treatment decreased according to the severity of AR.

Background

A left ventricular assist device (LVAD) is normally contraindicated in significant aortic regurgitation (AR) and requires intraoperative valve repair or exclusion. Nevertheless, AR can coexist with an LVAD, so a valid question when asked might still be of clinical significance. The purpose of this study is to analyze the effects of valve regurgitation on the pumping efficacy of continuous and pulsatile LVADs with a computational method.

Methods

A cardiovascular model was developed based on the Windkessel model, which reflects the hemodynamic flow resistance and the blood wall elasticity. Using the Windkessel model, important cardiovascular components, such as the right atrium, right ventricle, pulmonary artery, pulmonary vein, left atrium (LA), left ventricle (LV), aorta, and branching blood vessels, were expressed.

Results

In the case of AR, continuous and pulsatile LVADs improved cardiac output and reduced mechanical load slightly. In the case of mitral regurgitation, the LVADs improved cardiac output (cardiac outputs were about 5 L/min regardless of the severity of regurgitation) and reduced afterload significantly.

Conclusion

AR reduced both continuous and pulsatile LVAD function significantly while mitral regurgitation did not affect their pumping efficacy.

A New MRI-Based Model of Heart Function with Coupled Hemodynamics and Application to Normal and Diseased Canine Left Ventricles

A methodology for the simulation of heart function that combines an MRI-based model of cardiac electromechanics (CE) with a Navier–Stokes-based hemodynamics model is presented. The CE model consists of two coupled components that simulate the electrical and the mechanical functions of the heart. Accurate representations of ventricular geometry and fiber orientations are constructed from the structural magnetic resonance and the diffusion tensor MR images, respectively. The deformation of the ventricle obtained from the electromechanical model serves as input to the hemodynamics model in this one-way coupled approach via imposed kinematic wall velocity boundary conditions and at the same time, governs the blood flow into and out of the ventricular volume. The time-dependent endocardial surfaces are registered using a diffeomorphic mapping algorithm, while the intraventricular blood flow patterns are simulated using a sharp-interface immersed boundary method-based flow solver. The utility of the combined heart-function model is demonstrated by comparing the hemodynamic characteristics of a normal canine heart beating in sinus rhythm against that of the dyssynchronously beating failing heart. We also discuss the potential of coupled CE and hemodynamics models for various clinical applications.