Subject-specific computational simulation of left ventricular flow based on magnetic resonance imaging

Impact of left atrial wall motion assumptions in fluid simulations on proposed predictors of thrombus formation

Atrial fibrillation (AF) poses a significant risk of stroke due to thrombus formation, which primarily occurs in the left atrial appendage (LAA). Medical image-based computational fluid dynamics (CFD) simulations can provide valuable insight into patient-specific hemodynamics and could potentially enhance personalized assessment of thrombus risk. However, the importance of accurately representing the left atrial (LA) wall dynamics has not been fully resolved. In this study, we compared four modeling scenarios; rigid walls, a generic wall motion based on a reference motion, a semi-generic wall motion based on patient-specific motion, and patient-specific wall motion based on medical images. We considered a LA geometry acquired from 4D computed tomography during AF, systematically performed convergence tests to assess the numerical accuracy of our solution strategy, and quantified the differences between the four approaches. The results revealed that wall motion had no discernible impact on LA cavity hemodynamics, nor on the markers that indicate thrombus formation. However, the flow patterns within the LAA deviated significantly in the rigid model, indicating that the assumption of rigid walls may lead to errors in the estimated risk factors. In contrast, the generic, semi-generic, and patient-specific cases were qualitatively similar. The results highlight the crucial role of wall motion on hemodynamics and predictors of thrombus formation, and also demonstrate the potential of using a generic motion model as a surrogate for the more complex patient-specific motion. While the present study considered a single case, the employed CFD framework is entirely open-source and designed for adaptability, allowing for integration of additional models and generic motions.

A CFD study on the interplay of torsion and vortex guidance by the mitral valve on the left ventricular wash-out making use of overset meshes (Chimera technique)

Cardiovascular disease often occurs with silent and gradual alterations of cardiac blood flow that can lead to the onset of chronic pathological conditions. Image-based patient-specific Computational Fluid Dynamics (CFD) models allow for an extensive quantification of the flow field beyond the direct capabilities of medical imaging techniques that could support the clinicians in the early diagnosis, follow-up, and treatment planning of patients. Nonetheless, the large and impulsive kinematics of the left ventricle (LV) and the mitral valve (MV) pose relevant modeling challenges. Arbitrary Lagrangian-Eulerian (ALE) based computational fluid dynamics (CFD) methods struggle with the complex 3D mesh handling of rapidly moving valve leaflets within the left ventricle (LV). We, therefore, developed a Chimera-based (overset meshing) method to build a patient-specific 3D CFD model of the beating LV which includes a patient-inspired kinematic model of the mitral valve (LVMV). Simulations were performed with and without torsion. In addition, to evaluate how the intracardiac LV flow is impacted by the MV leaflet kinematics, a third version of the model without the MV was generated (LV with torsion). For all model versions, six cardiac cycles were simulated. All simulations demonstrated cycle-to-cycle variations that persisted after six cycles but were albeit marginal in terms of the magnitude of standard deviation of velocity and vorticity which may be related to the dissipative nature of the numerical scheme used. The MV was found to have a crucial role in the development of the intraventricular flow by enhancing the direct flow, the apical washout, and the propagation of the inlet jet towards the apical region. Consequently, the MV is an essential feature in the patient-specific CFD modeling of the LV. The impact of torsion was marginal on velocity, vorticity, wall shear stress, and energy loss, whereas it resulted to be significant in the evaluation of particle residence times. Therefore, including torsion could be considered in patient-specific CFD models of the LV, particularly when aiming to study stasis and residence time. We conclude that, despite some technical limitations encountered, the Chimera technique is a promising alternative for ALE methods for 3D CFD models of the heart that include the motion of valve leaflets.

Non-Newtonian blood rheology impacts left atrial stasis in patient-specific simulations

The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects are subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered when calculating patient-specific blood stasis indices by computational fluid dynamics.

Numerical simulations of flow patterns in the human left ventricle model with a novel dynamic mesh morphing approach based on radial basis function

We present a new numerical simulation framework for prediction of flow patterns in the human left ventricle model. In this study, a radial basis function (RBF) mesh morphing method is developed and applied within the finite-volume computational fluid dynamics (CFD) approach. The numerical simulations are designed to closely mimic details of recent tomographic particle image velocimetry (TomoPIV) experiments. The numerically simulated dynamic motions of the left ventricle and tri-leaflet biological mitral valve are emulated through the RBF morphing method. The arbitrary Lagrangian-Eulerian (ALE) based CFD is performed with the RBF-defined deforming wall boundaries. The results obtained show a good agreement with experiments, confirming the reliability and accuracy of the developed simulation framework.

Impact of Aortoseptal Angle Abnormalities and Discrete Subaortic Stenosis on Left-Ventricular Outflow Tract Hemodynamics: Preliminary Computational Assessment

Discrete subaortic stenosis (DSS) is an obstruction of the left ventricular outflow tract (LVOT) due to the formation of a fibromuscular membrane upstream of the aortic valve. DSS is a major risk factor for aortic regurgitation (AR), which often persists after surgical resection of the membrane. While the etiology of DSS and secondary AR is largely unknown, the frequent association between DSS and aortoseptal angle (AoSA) abnormalities has supported the emergence of a mechanobiological pathway by which hemodynamic stress alterations on the septal wall could trigger a biological cascade leading to fibrosis and membrane formation. The resulting LVOT flow disturbances could activate the valve endothelium and contribute to AR. In an effort to assess this hypothetical mechano-etiology, this study aimed at isolating computationally the effects of AoSA abnormalities on septal wall shear stress (WSS), and the impact of DSS on LVOT hemodynamics. Two-dimensional computational fluid dynamics models featuring a normal AoSA (N-LV), a steep AoSA (S-LV), and a steep AoSA with a DSS lesion (DSS-LV) were designed to compute the flow in patient-specific left ventricles (LVs). Boundary conditions consisted of transient velocity profiles at the mitral inlet and LVOT outlet, and patient-specific LV wall motion. The deformation of the DSS lesion was computed using a two-way fluid-structure interaction modeling strategy. Turbulence was accounted for via implementation of the k-ω turbulence model. While the N-LV and S-LV models generated similar LVOT flow characteristics, the DSS-LV model resulted in an asymmetric LVOT jet-like structure, subaortic stenotic conditions (up to 2.4-fold increase in peak velocity, 45% reduction in effective jet diameter vs. N-LV/S-LV), increased vorticity (2.8-fold increase) and turbulence (5- and 3-order-of-magnitude increase in turbulent kinetic energy and Reynolds shear stress, respectively). The steep AoSA subjected the septal wall to a 23% and 69% overload in temporal shear magnitude and gradient, respectively, without any substantial change in oscillatory shear index. This study reveals the existence of WSS overloads on septal wall regions prone to DSS lesion formation in steep LVOTs, and the development of highly turbulent, stenotic and asymmetric flow in DSS LVOTs, which support a possible mechano etiology for DSS and secondary AR.

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.

Patient-Specific Left Ventricular Flow Simulations From Transthoracic Echocardiography: Robustness Evaluation and Validation Against Ultrasound Doppler and Magnetic Resonance Imaging

The combination of medical imaging with computational fluid dynamics (CFD) has enabled the study of 3-D blood flow on a patient-specific level. However, with models based on gated high-resolution data, the study of transient flows, and any model implementation into routine cardiac care, is challenging. This paper presents a novel pathway for patient-specific CFD modelling of the left ventricle (LV), using 4-D transthoracic echocardiography (TTE) as input modality. To evaluate the clinical usability, two sub-studies were performed. First, a robustness evaluation was performed, where repeated models with alternating input variables were generated for six subjects and changes in simulated output quantified. Second, a validation study was carried out, where the pathway accuracy was evaluated against pulsed-wave Doppler (100 subjects), and 2-D through-plane phase-contrast magnetic resonance imaging measurements over seven intraventricular planes (6 subjects). The robustness evaluation indicated a model deviation of <12%, with highest regional and temporal deviations at apical segments and at peak systole, respectively. The validation study showed an error of <11% (velocities <10 cm/s) for all subjects, with no significant regional or temporal differences observed. With the patient-specific pathway shown to provide robust output with high accuracy, and with the pathway dependent only on 4-D TTE, the method has a high potential to be used within future clinical studies on 3-D intraventricular flow patterns. To this, future model developments in the form of e.g., anatomically accurate LV valves may further enhance the clinical value of the simulations.

Left ventricular flow propagation velocity measurement: Is it cast in stone?

This study aims to investigate the measurement of left ventricular flow propagation velocity, V _p, using phase contrast magnetic resonance imaging and to assess the discrepancies resulting from inflow jet direction and individual left ventricular size. Three V _p measuring techniques, namely non-adaptive (NA), adaptive positions (AP) and adaptive vectors (AV) method, were suggested and compared. We performed the comparison on nine healthy volunteers and nine post-infarct patients at four measurement positions, respectively, at one-third, one-half, two-thirds and the conventional 4 cm distances from the mitral valve leaflet into the left ventricle. We found that the V _p measurement was affected by both the inflow jet direction and measurement positions. Both NA and AP methods overestimated V _p, especially in dilated left ventricles, while the AV method showed the strongest correlation with the isovolumic relaxation myocardial strain rate (r = 0.53, p < 0.05). Using the AV method, notable difference in mean V _p was also observed between healthy volunteers and post-infarct patients at positions of: one-half (81 ± 31 vs. 58 ± 25 cm/s), two-thirds (89 ± 32 vs. 45 ± 15 cm/s) and 4 cm (98 ± 23 vs. 47 ± 13 cm/s) distances. The use of AV method and measurement position at one-half distance was found to be the most suitable method for assessing diastolic dysfunction given varying left ventricular sizes and inflow jet directions.

Patient-Specific Simulation of Cardiac Blood Flow From High-Resolution Computed Tomography

Cardiac hemodynamics can be computed from medical imaging data, and results could potentially aid in cardiac diagnosis and treatment optimization. However, simulations are often based on simplified geometries, ignoring features such as papillary muscles and trabeculae due to their complex shape, limitations in image acquisitions, and challenges in computational modeling. This severely hampers the use of computational fluid dynamics in clinical practice. The overall aim of this study was to develop a novel numerical framework that incorporated these geometrical features. The model included the left atrium, ventricle, ascending aorta, and heart valves. The framework used image registration to obtain patient-specific wall motion, automatic remeshing to handle topological changes due to the complex trabeculae motion, and a fast interpolation routine to obtain intermediate meshes during the simulations. Velocity fields and residence time were evaluated, and they indicated that papillary muscles and trabeculae strongly interacted with the blood, which could not be observed in a simplified model. The framework resulted in a model with outstanding geometrical detail, demonstrating the feasibility as well as the importance of a framework that is capable of simulating blood flow in physiologically realistic hearts.

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.

Patient-specific CFD simulation of intraventricular haemodynamics based on 3D ultrasound imaging

BACKGROUND

The goal of this paper is to present a computational fluid dynamic (CFD) model with moving boundaries to study the intraventricular flows in a patient-specific framework. Starting from the segmentation of real-time transesophageal echocardiographic images, a CFD model including the complete left ventricle and the moving 3D mitral valve was realized. Their motion, known as a function of time from the segmented ultrasound images, was imposed as a boundary condition in an Arbitrary Lagrangian-Eulerian framework.

RESULTS

The model allowed for a realistic description of the displacement of the structures of interest and for an effective analysis of the intraventricular flows throughout the cardiac cycle. The model provides detailed intraventricular flow features, and highlights the importance of the 3D valve apparatus for the vortex dynamics and apical flow.

CONCLUSIONS

The proposed method could describe the haemodynamics of the left ventricle during the cardiac cycle. The methodology might therefore be of particular importance in patient treatment planning to assess the impact of mitral valve treatment on intraventricular flow dynamics.

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.

A novel methodology for personalized simulations of ventricular hemodynamics from noninvasive imaging data

Current state-of-the-art imaging techniques can provide quantitative information to characterize ventricular function within the limits of the spatiotemporal resolution achievable in a realistic acquisition time. These imaging data can be used to personalize computer models, which in turn can help treatment planning by quantifying biomarkers that cannot be directly imaged, such as flow energy, shear stress and pressure gradients. To date, computer models have typically relied on invasive pressure measurements to be made patient-specific. When these data are not available, the scope and validity of the models are limited. To address this problem, we propose a new methodology for modeling patient-specific hemodynamics based exclusively on noninvasive velocity and anatomical data from 3D+t echocardiography or Magnetic Resonance Imaging (MRI). Numerical simulations of the cardiac cycle are driven by the image-derived velocities prescribed at the model boundaries using a penalty method that recovers a physical solution by minimizing the energy imparted to the system. This numerical approach circumvents the mathematical challenges due to the poor conditioning that arises from the imposition of boundary conditions on velocity only. We demonstrate that through this technique we are able to reconstruct given flow fields using Dirichlet only conditions. We also perform a sensitivity analysis to investigate the accuracy of this approach for different images with varying spatiotemporal resolution. Finally, we examine the influence of noise on the computed result, showing robustness to unbiased noise with an average error in the simulated velocity approximately 7% for a typical voxel size of 2mm(3) and temporal resolution of 30ms. The methodology is eventually applied to a patient case to highlight the potential for a direct clinical translation.

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.

Kinematic, Dynamic, and Energy Characteristics of Diastolic Flow in the Left Ventricle

Blood flow characteristics in the normal left ventricle are studied by using the magnetic resonance imaging, the Navier-Stokes equations, and the work-energy equation. Vortices produced during the mitral valve opening and closing are modeled in a two-dimensional analysis and correlated with temporal variations of the Reynolds number and pressure drop. Low shear stress and net pressures on the mitral valve are obtained for flow acceleration and deceleration. Bernoulli energy flux delivered to blood from ventricular dilation is practically balanced by the energy influx and the rate change of kinetic energy in the ventricle. The rates of work done by shear and energy dissipation are small. The dynamic and energy characteristics of the 2D results are comparable to those of a 3D model.

FLUID DYNAMIC CHARACTERISTICS OF SYSTOLIC BLOOD FLOW OF THE LEFT VENTRICLE

Ejection of blood from the left ventricle to the aorta is studied using two-dimensional Navier–Stokes equations, the work-energy equation and the magnetic resonance imaging of a normal ventricular motion. Vortex shedding in the sinuses of Valsalva is dominated by the aortic jet, flow acceleration and valve motion. Momentums produced by ventricular contraction are in concert with vortices in the ventricle for blood ejection. Shear stresses and net pressures on the aortic valve are calculated during valve opening and closing. The rate of work done by shear and the energy dissipation in the ventricle are small. The Bernoulli energy flux delivered to blood from ventricular contraction is practically balanced by energy flux at the aortic root and the rate change of kinetic energy in the ventricle.

A semi-automated method for patient-specific computational flow modelling of left ventricles

Patient-specific computational fluid dynamics (CFD) modelling of the left ventricle (LV) is a promising technique for the visualisation of ventricular flow patterns throughout a cardiac cycle. While significant progress has been made in improving the physiological quality of such simulations, the methodologies involved for several key steps remain significantly operator-dependent to this day. This dependency limits both the efficiency of the process as well as the consistency of CFD results due to the labour-intensive nature of current methods as well as operator introduced uncertainties in the modelling process. In order to mitigate this dependency, we propose a semi-automated method for patient-specific computational flow modelling of the LV. Using magnetic resonance imaging derived coarse geometry data of a patient's LV endocardium shape throughout a cardiac cycle, we then proceed to refine the geometry to eliminate rough edges before reconstructing meshes for all time frames and finally numerically solving for the intra-ventricular flow. Using a sample of patient-specific volunteer data, we demonstrate that our semi-automated, minimal operator involvement approach is capable of yielding CFD results of the LV that are comparable to other clinically validated LV flow models in the literature.

Fluid-structure interaction of an aortic heart valve prosthesis driven by an animated anatomic left ventricle

We develop a novel large-scale kinematic model for animating the left ventricle (LV) wall and use this model to drive the fluid-structure interaction (FSI) between the ensuing blood flow and a mechanical heart valve prosthesis implanted in the aortic position of an anatomic LV/aorta configuration. The kinematic model is of lumped type and employs a cell-based, FitzHugh-Nagumo framework to simulate the motion of the LV wall in response to an excitation wavefront propagating along the heart wall. The emerging large-scale LV wall motion exhibits complex contractile mechanisms that include contraction (twist) and expansion (untwist). The kinematic model is shown to yield global LV motion parameters that are well within the physiologic range throughout the cardiac cycle. The FSI between the leaflets of the mechanical heart valve and the blood flow driven by the dynamic LV wall motion and mitral inflow is simulated using the curvilinear immersed boundary (CURVIB) method [1, 2] implemented in conjunction with a domain decomposition approach. The computed results show that the simulated flow patterns are in good qualitative agreement with in vivo observations. The simulations also reveal complex kinematics of the valve leaflets, thus, underscoring the need for patient-specific simulations of heart valve prosthesis and other cardiac devices.

Review on CFD simulation in heart with dilated cardiomyopathy and myocardial infarction

The heart is a sophisticated functional organ that plays a crucial role in the blood circulatory system. Hemodynamics within the heart chamber can be indicative of exert cardiac health. Due to the limitations of current cardiac imaging modalities, computational fluid dynamics (CFD) have been widely used for the purposes of cardiac function assessment and heart disease diagnosis, as they provide detailed insights into the cardiac flow field. An understanding of ventricular hemodynamics and pathological severities can be gained through studies that employ the CFD method. In this research the hemodynamics of two common myocardial diseases, dilated cardiomyopathy (DCM) and myocardial infarction (MI) were investigated, during both the filling phase and the whole cardiac cycle, through a prescribed geometry and fluid structure interaction (FSI) approach. The results of the research indicated that early stage disease identification and the improvement of cardiac assisting devices and therapeutic procedures can be facilitated through the use of the CFD method.

INTRA-LEFT VENTRICULAR FLOW DISTRIBUTIONS IN DIASTOLIC AND SYSTOLIC PHASES, BASED ON ECHO VELOCITY FLOW MAPPING OF NORMAL SUBJECTS AND HEART FAILURE PATIENTS, TO CHARACTERIZE LEFT VENTRICULAR PERFORMANCE OUTCOMES OF HEART FAILURE

Heart failure (HF), one of the most common diseases in the world, causes left ventricular dysfunction (LV) and high mortality. HF patients are stratified into two groups based on their LV ejection fraction (EF) — HF with normal EF (HFNEF) and with reduced EF (HFREF). EF is a commonly used measure of LV contractile performance. Despite preserved EF, a complex mixture of systolic and diastolic dysfunction and variable degrees of LV remodelling underlying HFNEF poses challenges to diagnose and provide pharmacological treatment for HFNEF. In recent years, the velocity flow mapping (VFM) technique has been developed to generate flow velocity vector fields by post-processing color Doppler echocardiographic (echo) images. We aim to obtain the intra-LV blood flow patterns for patients with HFNEF, HFREF, and normal subjects, in order to characterize the LV performance outcomes of normal subjects and HF patients.

Two subjects from each group of HFNEF, HFREF, and normal underwent echo scans. Velocity vector distributions throughout the cardiac cycle were then analysed using the VFM technique. In each subject, the outflow rate during systole, inflow rate during diastole, as well as wall stress-based pressure-normalized contractility index, dσ*/dt max , were computed and compared among the groups.

This study demonstrated the use of VFM to visualize LV blood flow patterns in HF patients and normal subjects. Different patterns of flow distributions were observed in these subjects. In HFREF patients, dσ*/dt max , the peak outflow rate and peak inflow rate during early filling were markedly reduced. In HFNEF patients, peak outflow rates were increased compared to those of normal subjects.

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian-Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.

An image-based model of the whole human heart with detailed anatomical structure and fiber orientation

Many heart anatomy models have been developed to study the electrophysiological properties of the human heart. However, none of them includes the geometry of the whole human heart. In this study, an anatomically detailed mathematical model of the human heart was firstly reconstructed from the computed tomography images. In the reconstructed model, the atria consisted of atrial muscles, sinoatrial node, crista terminalis, pectinate muscles, Bachmann's bundle, intercaval bundles, and limbus of the fossa ovalis. The atrioventricular junction included the atrioventricular node and atrioventricular ring, and the ventricles had ventricular muscles, His bundle, bundle branches, and Purkinje network. The epicardial and endocardial myofiber orientations of the ventricles and one layer of atrial myofiber orientation were then measured. They were calculated using linear interpolation technique and minimum distance algorithm, respectively. To the best of our knowledge, this is the first anatomically-detailed human heart model with corresponding experimentally measured fibers orientation. In addition, the whole heart excitation propagation was simulated using a monodomain model. The simulated normal activation sequence agreed well with the published experimental findings.

Assessment of the arterial input function for estimation of coronary flow reserve by single photon emission computed tomography: comparison of two different approaches

PURPOSE

Attempts to estimate coronary flow reserve (CFR) with single photon emission computed tomography (SPECT) tracers have been recently made. We compared two different methods for the estimation of CFR by SPECT imaging.

METHODS

Fourteen patients with coronary artery disease underwent dipyridamole 99mTc-sestamibi SPECT and intracoronary Doppler within 5 days. Myocardial blood flow (MBF) was estimated by measurement of first transit counts in the right pulmonary artery (PA) and left ventricular (LV) chamber, and myocardial counts from SPECT images. Estimated CFR was expressed as the ratio of stress MBF to rest MBF.

RESULTS

Rest and stress MBF obtained using first transit counts from PA were higher compared to that from LV chamber (rest: 1.05 ± 0.38 vs 0.87 ± 0.34 counts/pixel per s, respectively, p < 0.01 and stress: 1.34 ± 0.45 vs 0.91 ± 0.20 counts/pixel per s, respectively, p < 0.05). In the study vessels, CFR by Doppler was 1.39 ± 0.42, and SPECT CFR obtained using first transit counts from PA and LV chamber were 1.36 ± 0.43 and 1.16 ± 0.39, respectively (p across categories NS). A significant relationship between SPECT CFR obtained using first transit counts from PA and CFR by Doppler was found (r = 0.85, p < 0.001). No relationship between SPECT CFR obtained using first transit counts from LV chamber and CFR by intracoronary Doppler was OBSERVED (R = 0.43, P = NS).

CONCLUSION

SPECT-estimated CFR obtained using first transit counts from right PA is more accurate and correlates better with the results of intracoronary Doppler than estimated CFR obtained using arterial input function from LV chamber.

Multi-Physics MRI-Based Two-Layer Fluid-Structure Interaction Anisotropic Models of Human Right and Left Ventricles with Different Patch Materials: Cardiac Function Assessment and Mechanical Stress Analysis

Multi-physics right and left ventricle (RV/LV) fluid-structure interaction (FSI) models were introduced to perform mechanical stress analysis and evaluate the effect of patch materials on RV function. The FSI models included three different patch materials (Dacron scaffold, treated pericardium, and contracting myocardium), two-layer construction, fiber orientation, and active anisotropic material properties. The models were constructed based on cardiac magnetic resonance (CMR) images acquired from a patient with severe RV dilatation and solved by ADINA. Our results indicate that the patch model with contracting myocardium leads to decreased stress level in the patch area, improved RV function and patch area contractility.

CFD simulation of flow through heart: a perspective review

The heart is an organ which pumps blood around the body by contraction of muscular wall. There is a coupled system in the heart containing the motion of wall and the motion of blood fluid; both motions must be computed simultaneously, which make biological computational fluid dynamics (CFD) difficult. The wall of the heart is not rigid and hence proper boundary conditions are essential for CFD modelling. Fluid-wall interaction is very important for real CFD modelling. There are many assumptions for CFD simulation of the heart that make it far from a real model. A realistic fluid-structure interaction modelling the structure by the finite element method and the fluid flow by CFD use more realistic coupling algorithms. This type of method is very powerful to solve the complex properties of the cardiac structure and the sensitive interaction of fluid and structure. The final goal of heart modelling is to simulate the total heart function by integrating cardiac anatomy, electrical activation, mechanics, metabolism and fluid mechanics together, as in the computational framework.

Image-Based Patient-Specific Ventricle Models with Fluid-Structure Interaction for Cardiac Function Assessment and Surgical Design Optimization

Recent advances in medical imaging technology and computational modeling techniques are making it possible that patient-specific computational ventricle models be constructed and used to test surgical hypotheses and replace empirical and often risky clinical experimentation to examine the efficiency and suitability of various reconstructive procedures in diseased hearts. In this paper, we provide a brief review on recent development in ventricle modeling and its potential application in surgical planning and management of tetralogy of Fallot (ToF) patients. Aspects of data acquisition, model selection and construction, tissue material properties, ventricle layer structure and tissue fiber orientations, pressure condition, model validation and virtual surgery procedures (changing patient-specific ventricle data and perform computer simulation) were reviewed. Results from a case study using patient-specific cardiac magnetic resonance (CMR) imaging and right/left ventricle and patch (RV/LV/Patch) combination model with fluid-structure interactions (FSI) were reported. The models were used to evaluate and optimize human pulmonary valve replacement/insertion (PVR) surgical procedure and patch design and test a surgical hypothesis that PVR with small patch and aggressive scar tissue trimming in PVR surgery may lead to improved recovery of RV function and reduced stress/strain conditions in the patch area.