Three-dimensional CFD/MRI modeling reveals that ventricular surgical restoration improves ventricular function by modifying intraventricular blood flow

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

CT-Based Analysis of Left Ventricular Hemodynamics Using Statistical Shape Modeling and Computational Fluid Dynamics

Background

Cardiac computed tomography (CCT) based computational fluid dynamics (CFD) allows to assess intracardiac flow features, which are hypothesized as an early predictor for heart diseases and may support treatment decisions. However, the understanding of intracardiac flow is challenging due to high variability in heart shapes and contractility. Using statistical shape modeling (SSM) in combination with CFD facilitates an intracardiac flow analysis. The aim of this study is to prove the usability of a new approach to describe various cohorts.

Materials and Methods

CCT data of 125 patients (mean age: 60.6 ± 10.0 years, 16.8% woman) were used to generate SSMs representing aneurysmatic and non-aneurysmatic left ventricles (LVs). Using SSMs, seven group-averaged LV shapes and contraction fields were generated: four representing patients with and without aneurysms and with mild or severe mitral regurgitation (MR), and three distinguishing aneurysmatic patients with true, intermediate aneurysms, and globally hypokinetic LVs. End-diastolic LV volumes of the groups varied between 258 and 347 ml, whereas ejection fractions varied between 21 and 26%. MR degrees varied from 1.0 to 2.5. Prescribed motion CFD was used to simulate intracardiac flow, which was analyzed regarding large-scale flow features, kinetic energy, washout, and pressure gradients.

Results

SSMs of aneurysmatic and non-aneurysmatic LVs were generated. Differences in shapes and contractility were found in the first three shape modes. Ninety percent of the cumulative shape variance is described with approximately 30 modes. A comparison of hemodynamics between all groups found shape-, contractility- and MR-dependent differences. Disturbed blood washout in the apex region was found in the aneurysmatic cases. With increasing MR, the diastolic jet becomes less coherent, whereas energy dissipation increases by decreasing kinetic energy. The poorest blood washout was found for the globally hypokinetic group, whereas the weakest blood washout in the apex region was found for the true aneurysm group.

Conclusion

The proposed CCT-based analysis of hemodynamics combining CFD with SSM seems promising to facilitate the analysis of intracardiac flow, thus increasing the value of CCT for diagnostic and treatment decisions. With further enhancement of the computational approach, the methodology has the potential to be embedded in clinical routine workflows and support clinicians.

CT-Based Simulation of Left Ventricular Hemodynamics: A Pilot Study in Mitral Regurgitation and Left Ventricle Aneurysm Patients

Background

Cardiac CT (CCT) is well suited for a detailed analysis of heart structures due to its high spatial resolution, but in contrast to MRI and echocardiography, CCT does not allow an assessment of intracardiac flow. Computational fluid dynamics (CFD) can complement this shortcoming. It enables the computation of hemodynamics at a high spatio-temporal resolution based on medical images. The aim of this proposed study is to establish a CCT-based CFD methodology for the analysis of left ventricle (LV) hemodynamics and to assess the usability of the computational framework for clinical practice.

Materials and Methods

The methodology is demonstrated by means of four cases selected from a cohort of 125 multiphase CCT examinations of heart failure patients. These cases represent subcohorts of patients with and without LV aneurysm and with severe and no mitral regurgitation (MR). All selected LVs are dilated and characterized by a reduced ejection fraction (EF). End-diastolic and end-systolic image data was used to reconstruct LV geometries with 2D valves as well as the ventricular movement. The intraventricular hemodynamics were computed with a prescribed-motion CFD approach and evaluated in terms of large-scale flow patterns, energetic behavior, and intraventricular washout.

Results

In the MR patients, a disrupted E-wave jet, a fragmentary diastolic vortex formation and an increased specific energy dissipation in systole are observed. In all cases, regions with an impaired washout are visible. The results furthermore indicate that considering several cycles might provide a more detailed view of the washout process. The pre-processing times and computational expenses are in reach of clinical feasibility.

Conclusion

The proposed CCT-based CFD method allows to compute patient-specific intraventricular hemodynamics and thus complements the informative value of CCT. The method can be applied to any CCT data of common quality and represents a fair balance between model accuracy and overall expenses. With further model enhancements, the computational framework has the potential to be embedded in clinical routine workflows, to support clinical decision making and treatment planning.

Clinical and hemodynamic outcomes of the Dor procedure in adults with ischemic cardiomyopathy

Ischemic cardiomyopathy continues to be a major contributor to congestive heart failure, which places a significant burden on our healthcare system. Improving medications and different coronary revascularization strategies are the mainstays in the management of ischemic cardiomyopathy. Although medications and mechanical circulatory support are playing an ever-increasing role, cardiac transplantation remains the gold standard for treating advanced heart failure. Given the small number of available and suitable donor hearts, transplantation is limited for the majority of patients. Surgical ventricular restoration has repeatedly been suggested as a viable alternative in managing heart failure in select patients, as it is believed that surgically returning the ventricle to its original dimensions is possible and associated with favorable outcomes. The purpose of this manuscript is to comprehensively review the current literature on various surgical strategies for ventricular restoration. We also contextualize the published data with respect to ventricular function, volume, structure, arrhythmias, mitral regurgitation, and clinical outcomes.

Effects of Choice of Medical Imaging Modalities on a Non-invasive Diagnostic and Monitoring Computational Framework for Patients With Complex Valvular, Vascular, and Ventricular Diseases Who Undergo Transcatheter Aortic Valve Replacement

Due to the high individual differences in the anatomy and pathophysiology of patients, planning individualized treatment requires patient-specific diagnosis. Indeed, hemodynamic quantification can be immensely valuable for accurate diagnosis, however, we still lack precise diagnostic methods for numerous cardiovascular diseases including complex (and mixed) valvular, vascular, and ventricular interactions (C3VI) which is a complicated situation made even more challenging in the face of other cardiovascular pathologies. Transcatheter aortic valve replacement (TAVR) is a new less invasive intervention and is a growing alternative for patients with aortic stenosis. In a recent paper, we developed a non-invasive and Doppler-based diagnostic and monitoring computational mechanics framework for C3VI, called C3VI-DE that uses input parameters measured reliably using Doppler echocardiography. In the present work, we have developed another computational-mechanics framework for C3VI (called C3VI-CT). C3VI-CT uses the same lumped-parameter model core as C3VI-DE but its input parameters are measured using computed tomography and a sphygmomanometer. Both frameworks can quantify: (1) global hemodynamics (metrics of cardiac function); (2) local hemodynamics (metrics of circulatory function). We compared accuracy of the results obtained using C3VI-DE and C3VI-CT against catheterization data (gold standard) using a C3VI dataset (N = 49) for patients with C3VI who undergo TAVR in both pre and post-TAVR with a high variability. Because of the dataset variability and the broad range of diseases that it covers, it enables determining which framework can yield the most accurate results. In contrast with C3VI-CT, C3VI-DE tracks both the cardiac and vascular status and is in great agreement with cardiac catheter data.

Blood flow patterns estimation in the left ventricle with low-rate 2D and 3D dynamic contrast-enhanced ultrasound

BACKGROUND AND OBJECTIVE

Left ventricle (LV) dysfunction always occurs at early heart-failure stages, producing variations in the LV flow patterns. Cardiac diagnostics may therefore benefit from flow-pattern analysis. Several visualization tools have been proposed that require ultrafast ultrasound acquisitions. However, ultrafast ultrasound is not standard in clinical scanners. Meanwhile techniques that can handle low frame rates are still lacking. As a result, the clinical translation of these techniques remains limited, especially for 3D acquisitions where the volume rates are intrinsically low.

METHODS

To overcome these limitations, we propose a novel technique for the estimation of LV blood velocity and relative-pressure fields from dynamic contrast-enhanced ultrasound (DCE-US) at low frame rates. Different from other methods, our method is based on the time-delays between time-intensity curves measured at neighbor pixels in the DCE-US loops. Using Navier-Stokes equation, we regularize the obtained velocity fields and derive relative-pressure estimates. Blood flow patterns were characterized with regard to their vorticity, relative-pressure changes (dp/dt) in the LV outflow tract, and viscous energy loss, as these reflect the ejection efficiency.

RESULTS

We evaluated the proposed method on 18 patients (9 responders and 9 non-responders) who underwent cardiac resynchronization therapy (CRT). After CRT, the responder group evidenced a significant (p<0.05) increase in vorticity and peak dp/dt, and a non-significant decrease in viscous energy loss. No significant difference was found in the non-responder group. Relative feature variation before and after CRT evidenced a significant difference (p<0.05) between responders and non-responders for vorticity and peak dp/dt. Finally, the method feasibility is also shown with 3D DCE-US.

CONCLUSIONS

Using the proposed method, adequate visualization and quantification of blood flow patterns are successfully enabled based on low-rate DCE-US of the LV, facilitating the clinical adoption of the method using standard ultrasound scanners. The clinical value of the method in the context of CRT is also shown.

The effect of chordae tendineae on systolic flow

When using Computational Fluid Dynamics to simulate ventricular blood flow in the heart, it has been common practice to neglect the effect of the sub-valvular apparatus and the trabeculae on the flow conditions. In this study, we analyze the effect of neglecting the chordae tendineae on the fluid flow and pressure drop. To test the assumption we use a previously developed dynamic 3D model of the left ventricle, aorta and valves that is based on 3D echocardiographic recordings. To this model we add the chordae tendineae as a sub-grid model. The previously developed 3D model for the left ventricle during systole is based on real-time three-dimensional echocardiography (RT3DE) recordings of a 30 years old female volunteer. The segmented ventricular wall does not include details of the aorta and the mitral valve, so these were reconstructed. The subgrid model for the flow across the chordae tendineae is based on the Actuator Line Method, which means that they are represented by drag coefficients. The analysis shows that the effect of the chordae tendineae on the pressure drop and work efficiency of the normal heart during systole is minor, and it seems that for simulating ventricular fluid flow and pressure drop during systole, one can follow the current practice and ignore the chordae. However, there can be local effects such as small vortices behind the chordae. Whether such effects are important for a particular application must be evaluated for the given case.

Application of Patient-Specific Computational Fluid Dynamics in Coronary and Intra-Cardiac Flow Simulations: Challenges and Opportunities

The emergence of new cardiac diagnostics and therapeutics of the heart has given rise to the challenging field of virtual design and testing of technologies in a patient-specific environment. Given the recent advances in medical imaging, computational power and mathematical algorithms, patient-specific cardiac models can be produced from cardiac images faster, and more efficiently than ever before. The emergence of patient-specific computational fluid dynamics (CFD) has paved the way for the new field of computer-aided diagnostics. This article provides a review of CFD methods, challenges and opportunities in coronary and intra-cardiac flow simulations. It includes a review of market products and clinical trials. Key components of patient-specific CFD are covered briefly which include image segmentation, geometry reconstruction, mesh generation, fluid-structure interaction, and solver techniques.

Assessment of human left ventricle flow using statistical shape modelling and computational fluid dynamics

Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.

Modeling Left Ventricular Blood Flow Using Smoothed Particle Hydrodynamics

This study aims to investigate the capability of smoothed particle hydrodynamics (SPH), a fully Lagrangian mesh-free method, to simulate the bulk blood flow dynamics in two realistic left ventricular (LV) models. Three dimensional geometries and motion of the LV, proximal left atrium and aortic root are extracted from cardiac magnetic resonance imaging and multi-slice computed tomography imaging data. SPH simulation results are analyzed and compared with those obtained using a traditional finite volume-based numerical method, and to in vivo phase contrast magnetic resonance imaging and echocardiography data, in terms of the large-scale blood flow phenomena usually clinically measured. A quantitative comparison of the velocity fields and global flow parameters between the in silico models and the in vivo data shows a reasonable agreement, given the inherent uncertainties and limitations in the modeling and imaging techniques. The results indicate the capability of SPH as a promising tool for predicting clinically relevant large-scale LV flow information.

Heart blood flow simulation: a perspective review

Cardiovascular disease (CVD), the leading cause of death today, incorporates a wide range of cardiovascular system malfunctions that affect heart functionality. It is believed that the hemodynamic loads exerted on the cardiovascular system, the left ventricle (LV) in particular, are the leading cause of CVD initiation and propagation. Moreover, it is believed that the diagnosis and prognosis of CVD at an early stage could reduce its high mortality and morbidity rate. Therefore, a set of robust clinical cardiovascular assessment tools has been introduced to compute the cardiovascular hemodynamics in order to provide useful insights to physicians to recognize indicators leading to CVD and also to aid the diagnosis of CVD. Recently, a combination of computational fluid dynamics (CFD) and different medical imaging tools, image-based CFD (IB-CFD), has been widely employed for cardiovascular functional assessment by providing reliable hemodynamic parameters. Even though the capability of CFD to provide reliable flow dynamics in general fluid mechanics problems has been widely demonstrated for many years, up to now, the clinical implications of the IB-CFD patient-specific LVs have not been applicable due to its limitations and complications. In this paper, we review investigations conducted to numerically simulate patient-specific human LV over the past 15 years using IB-CFD methods. Firstly, we divide different studies according to the different LV types (physiological and different pathological conditions) that have been chosen to reconstruct the geometry, and then discuss their contributions, methodologies, limitations, and findings. In this regard, we have studied CFD simulations of intraventricular flows and related cardiology insights, for (i) Physiological patient-specific LV models, (ii) Pathological heart patient-specific models, including myocardial infarction, dilated cardiomyopathy, hypertrophic cardiomyopathy and hypoplastic left heart syndrome. Finally, we discuss the current stage of the IB-CFD LV simulations in order to mimic realistic hemodynamics of patient-specific LVs. We can conclude that heart flow simulation is on the right track for developing into a useful clinical tool for heart function assessment, by (i) incorporating most of heart structures' (such as heart valves) operations, and (ii) providing useful diagnostic indices based hemodynamic parameters, for routine adoption in clinical usage.

Correcting motion in multiplanar cardiac magnetic resonance images

BACKGROUND

Misalignment in cardiac magnetic resonance (CMR) images can adversely affect three-dimensional left ventricle modelling and downstream quantitative analysis. Currently, there are two types of approaches for dealing with realignment and motion distortion problems, one image based and the other geometry based. Image-based approaches are limited by the inherent non-homogeneity and anisotropy of CMR images. Geometry-based approaches rely on idealized models and over-simplified assumptions. This study was motivated by the need for a robust and effective approach for correcting motion related distortions due to misalignment in CMR images.

METHODS

A cine cardiac magnetic resonance image sequence was acquired using our routine clinical imaging protocol. The left ventricular endocardium was delineated manually with software assistance on all long and short-axis images. Long and short-axis contours were projected onto a patient-based coordinate system and then realigned using iterative registration. The realigned contour points were used to reconstruct the shape of the left ventricle for quantitative validation.

RESULTS

The method was tested on five myocardial infarction patients whose images showed substantial misalignment. Realignment time was about 16 seconds per case, using a 2.5 GHz CPU desktop with obvious elimination of the distortion in the reconstructed model. Using the long-axis contour as a reference in evaluating the reconstructed models, it was apparent that the models with realigned contours had better accuracy than the non-realigned ones.

CONCLUSION

This study presents a novel, geometry-based method for correcting motion distortions in CMR images. The method incorporates (1) manual delineation, (2) registration based on a generalized, iterative closest point algorithm, and (3) reconstruction of the shape of the left ventricle for quantitative validation. The effectiveness of our approach is corroborated both visually and by quantitative assessment. We envision the use of our method in current clinical practice as a means of improving accuracy in the evaluation of cardiac function.

The numerical analysis of non-Newtonian blood flow in human patient-specific left ventricle

Recently, various non-invasive tools such as the magnetic resonance image (MRI), ultrasound imaging (USI), computed tomography (CT), and the computational fluid dynamics (CFD) have been widely utilized to enhance our current understanding of the physiological parameters that affect the initiation and the progression of the cardiovascular diseases (CVDs) associated with heart failure (HF). In particular, the hemodynamics of left ventricle (LV) has attracted the attention of the researchers due to its significant role in the heart functionality. In this study, CFD owing its capability of predicting detailed flow field was adopted to model the blood flow in images-based patient-specific LV over cardiac cycle. In most published studies, the blood is modeled as Newtonian that is not entirely accurate as the blood viscosity varies with the shear rate in non-linear manner. In this paper, we studied the effect of Newtonian assumption on the degree of accuracy of intraventricular hemodynamics. In doing so, various non-Newtonian models and Newtonian model are used in the analysis of the intraventricular flow and the viscosity of the blood. Initially, we used the cardiac MRI images to reconstruct the time-resolved geometry of the patient-specific LV. After the unstructured mesh generation, the simulations were conducted in the CFD commercial solver FLUENT to analyze the intraventricular hemodynamic parameters. The findings indicate that the Newtonian assumption cannot adequately simulate the flow dynamic within the LV over the cardiac cycle, which can be attributed to the pulsatile and recirculation nature of the flow and the low blood shear rate.

Cardiac MRI based numerical modeling of left ventricular fluid dynamics with mitral valve incorporated

Recent numerical studies were focused on the modeling of flow in patient-specific left ventricle (LV); however, the mitral valve (MV) was usually excluded. In this study, both patient-specific LV and MV were modeled to achieve a more realistic intraventricular flow. Cardiac MRI images were acquired from a pulmonary arterial hypertension (PAH) patient and a healthy volunteer, and manual segmentation was conducted to reconstruct three-dimensional (3D) LV and MV geometries at each frame. Based on these 3D geometries, vortex formation time (VFT) was derived, and the values were 4.0 and 6.5 for the normal subject and the PAH patient, respectively. Based on studies in the literature, VTF in the healthy subject fell within the normal range, while that in the PAH patient exceeded the threshold for normality. The vortex structures in the LV clearly showed that the vortex ring was initiated from the tips of the MV instead of the mitral annulus. The excessive VFT during the rapid filling phase in the PAH patient resulted in a trailing flow structure behind the primary vortex ring, which was not observed in the normal subject. It can be deduced from this study that incorporating the MV into a patient-specific model is necessary to produce more reasonable VFT and intraventricular flow.

Variational Reconstruction of Left Cardiac Structure from CMR Images

Cardiovascular Disease (CVD), accounting for 17% of overall deaths in the USA, is the leading cause of death over the world. Advances in medical imaging techniques make the quantitative assessment of both the anatomy and function of heart possible. The cardiac modeling is an invariable prerequisite for quantitative analysis. In this study, a novel method is proposed to reconstruct the left cardiac structure from multi-planed cardiac magnetic resonance (CMR) images and contours. Routine CMR examination was performed to acquire both long axis and short axis images. Trained technologists delineated the endocardial contours. Multiple sets of two dimensional contours were projected into the three dimensional patient-based coordinate system and registered to each other. The union of the registered point sets was applied a variational surface reconstruction algorithm based on Delaunay triangulation and graph-cuts. The resulting triangulated surfaces were further post-processed. Quantitative evaluation on our method was performed via computing the overlapping ratio between the reconstructed model and the manually delineated long axis contours, which validates our method. We envisage that this method could be used by radiographers and cardiologists to diagnose and assess cardiac function in patients with diverse heart diseases.

Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

This work presents a numerical simulation of intraventricular flow after the implantation of a bileaflet mechanical heart valve at the mitral position. The left ventricle was simplified conceptually as a truncated prolate spheroid and its motion was prescribed based on that of a healthy subject. The rigid leaflet rotation was driven by the transmitral flow and hence the leaflet dynamics were solved using fluid-structure interaction approach. The simulation results showed that the bileaflet mechanical heart valve at the mitral position behaved similarly to that at the aortic position. Sudden area expansion near the aortic root initiated a clockwise anterior vortex, and the continuous injection of flow through the orifice resulted in further growth of the anterior vortex during diastole, which dominated the intraventricular flow. This flow feature is beneficial to preserving the flow momentum and redirecting the blood flow towards the aortic valve. To the best of our knowledge, this is the first attempt to numerically model intraventricular flow with the mechanical heart valve incorporated at the mitral position using a fluid-structure interaction approach. This study facilitates future patient-specific studies.