Comparative numerical study on left ventricular fluid dynamics after dilated cardiomyopathy

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.

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.

Surrogate models provide new insights on metrics based on blood flow for the assessment of left ventricular function

Recent developments on the grading of cardiac pathologies suggest flow-related metrics for a deeper evaluation of cardiac function. Blood flow evaluation employs space-time resolved cardiovascular imaging tools, possibly integrated with direct numerical simulation (DNS) of intraventricular fluid dynamics in individual patients. If a patient-specific analysis is a promising method to reproduce flow details or to assist virtual therapeutic solutions, it becomes impracticable in nearly-real-time during a routine clinical activity. At the same time, the need to determine the existence of relationships between advanced flow-related quantities of interest (QoIs) and the diagnostic metrics used in the standard clinical practice requires the adoption of techniques able to generalize evidences emerging from a finite number of single cases. In this study, we focus on the left ventricular function and use a class of reduced-order models, relying on the Polynomial Chaos Expansion (PCE) technique to learn the dynamics of selected QoIs based on a set of synthetic cases analyzed with a high-fidelity model (DNS). The selected QoIs describe the left ventricle blood transit and the kinetic energy and vorticity at the peak of diastolic filling. The PCE-based surrogate models provide straightforward approximations of these QoIs in the space of widely used diagnostic metrics embedding relevant information on left ventricle geometry and function. These surrogates are directly employable in the clinical analysis as we demonstrate by assessing their robustness against independent patient-specific cases ranging from healthy to diseased conditions. The surrogate models are used to perform global sensitivity analysis at a negligible computational cost and provide insights on the impact of each diagnostic metric on the QoIs. Results also suggest how common flow transit parameters are principally dictated by ejection fraction.

Noninvasive Evaluation of Intraventricular Flow Dynamics by the HyperDoppler Technique: First Application to Normal Subjects, Athletes, and Patients with Heart Failure

Background: HyperDoppler is a new echocardiographic color Doppler-based technique that can assess intracardiac flow dynamics. The aim of this study was to verify the feasibility and reproducibility of this technique in unselected patients and its capability to differentiate measures of vortex flow within the left ventricle (LV) in normal sedentary subjects, athletes, and patients with heart failure. Methods: Two hundred unselected, consecutive patients presenting at the echocardiographic laboratory, 50 normal subjects, 30 athletes, and 50 patients with chronic heart failure and LV ejection fraction <50% were enrolled. Images were acquired using a MyLab X8 echo-scanner. Area, intensity, depth, length, and kinetic energy dissipation (KED) of vortex flow were measured. Results: The HyperDoppler technique feasibility was 94.5%. According to the intraclass correlation coefficient evaluations, repeatability and reproducibility of vortex flow measures were good for vortex area (0.82, 0.85), length (0.83, 0.82), and depth (0.87, 0.84) and excellent for intensity (0.92, 0.90) and KED (0.98, 0.98). Combining different vortex flow measures, the LV flow profile of healthy sedentary individuals, athletes, and heart failure patients could be differentiated. Conclusions: HyperDoppler is a feasible, reliable, and practical technique for the assessment of LV flow dynamics and may distinguish normal subjects and patients with heart failure.

Intraventricular Flow Simulations in Singular Right Ventricles Reveal Deteriorated Washout and Low Vortex Formation

Purpose

Patients with a functionally univentricular heart represent one of the most common severe cardiac lesions with a prevalence of 3 per 10,000 live births. Hemodynamics of the singular ventricle is a major research topic in cardiology and there exists a relationship between fluid dynamical features and cardiac behavior in health and disease. The aim of the present work was to compare intraventricular flow in single right ventricle (SRV) patients and subjects with healthy left hearts (LV) through patient-specific CFD simulations.

Methods

Three-dimensional real-time echocardiographic images were obtained for five SRV patients and two healthy subjects and CFD simulations with a moving mesh methodology were performed. Intraventricular vortex formation and vortex formation time (VFT) as well as the turbulent kinetic energy (TKE) and ventricular washout were evaluated.

Results

The results show significantly lower values for the VFT and the TKE in SRV patients compared with healthy LV subjects. Furthermore, vortex formation does not progress to the apex in SRV patients. These findings were confirmed by a significantly lower washout in SRV patients.

Conclusions

The study pinpoints the intriguing role of intraventricular flows to characterize performance of SRVs that goes beyond standard clinical metrics such as ejection fraction.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00598-9.

Comparative Analysis of Right Ventricle Fluid Dynamics

The right and left sides of the human heart operate with a common timing and pump the same amount of blood. Therefore, the right ventricle (RV) presents a function that is comparable to the left ventricle (LV) in terms of flow generation; nevertheless, the RV operates against a much lower arterial pressure (afterload) and requires a lower muscular strength. This study compares the fluid dynamics of the normal right and left ventricles to better understand the role of the RV streamlined geometry and provide some physics-based ground for the construction of clinical indicators for the right side. The analysis is performed by image-based direct numerical simulation, using the immersed boundary technique including the simplified models of tricuspid and mitral valves. Results demonstrated that the vortex formation process during early diastole is similar in the two ventricles, then the RV vorticity rapidly dissipates in the subvalvular region while the LV sustains a weak circulatory pattern at the center of the chamber. Afterwards, during the systolic contraction, the RV geometry allows an efficient transfer of mechanical work to the propelled blood; differently from the LV, this work is non-negligible in the global energetic balance. The varying behavior of the RV, from reservoir to conduct, during the different phases of the heartbeat is briefly discussed in conjunction to the development of possible dysfunctions.

Left Ventricular Deformation and Vortex Analysis in Heart Failure: From Ultrasound Technique to Current Clinical Application

Heart failure (HF) is a leading cause of cardiovascular morbidity and mortality. However, its symptoms and signs are not specific or can be absent. In this context, transthoracic echocardiography plays a key role in diagnosing the various forms of HF, guiding therapeutic decision making and monitoring response to therapy. Over the last few decades, new ultrasound modalities have been introduced in the field of echocardiography, aiming at better understanding the morpho-functional abnormalities occurring in cardiovascular diseases. However, they are still struggling to enter daily and routine use. In our review article, we turn the spotlight on some of the newest ultrasound technologies; in particular, analysis of myocardial deformation by speckle tracking echocardiography, and intracardiac flow dynamics by color Doppler flow mapping, highlighting their promising applications to HF diagnosis and management. We also focus on the importance of these imaging modalities in the selection of responses to cardiac resynchronization therapy.

Time-Resolved Echo-Particle Image/Tracking Velocimetry Measurement of Interactions Between Native Cardiac Output and Veno-Arterial ECMO Flows

Determination of optimal hemodynamic and pressure-volume loading conditions for patients undergoing veno-arterial extracorporeal membrane oxygenation (VA-ECMO) would benefit from understanding the impact of ECMO flow rates (QE) on the native cardiac output in the admixing zone, i.e., aortic root. This study characterizes the flow in the aortic root of a pig with severe myocardial ischemia using contrast-enhanced ultrasound particle image/tracking velocimetry (echo-PIV/PTV). New methods for data preprocessing are introduced, including autocontouring to remove surrounding tissues, followed by blind deconvolution to identify the centers of elongated bubble traces in images with low signal to noise ratio. Calibrations based on synthetic images show that this procedure increases the number of detected bubbles and reduces the error in their locations by 50%. Then, an optimized echo-PIV/PTV procedure, which integrates image enhancement with velocity measurements, is used for characterizing the time-resolved two-dimensional (2D) velocity distributions. Phase-averaged and instantaneous flow fields show that the ECMO flow rate influences the velocity and acceleration of the cardiac output during systole, and secondary flows during diastole. When QE is 3.0 L/min or higher, the cardiac ejection velocity, phase interval with open aortic valve, velocity-time integral (VTI), and mean arterial pressure (MAP) increase with decreasing QE, all indicating sufficient support. For lower QE, the MAP and VTI decrease as QE is reduced, and the deceleration during transition to diastole becomes milder. Hence, for this specific case, the optimal ECMO flow rate is 3.0 L/min.

Impact of synchronous atrioventricular delay optimization on left ventricle flow force angle evaluated by echocardiographic particle image velocimetry

PURPOSE

To evaluate the improvement in electrical synchrony and left ventricle (LV) hemodynamics provided by combining the dynamic atrioventricular delay (AVD) of SyncAV^TM CRT and the multiple LV pacing sites of MultiPoint pacing (MPP).

METHODS

Patients with LBBB and QRS duration (QRSd) > 140 ms implanted with a CRT-D or CRT-P device and quadripolar LV lead were enrolled in this prospective study. During a post-implant follow-up visit, QRSd was measured from 12-lead surface electrograms by experts blinded to pacing configurations. QRSd reduction relative to intrinsic rhythm was evaluated during biventricular pacing (BiV) and MPP for two AVDs: nominal (140/110 ms paced/sensed) and SyncAV (patient-optimized SyncAV offset [10-60 ms] minimizing QRSd). Echocardiography particle imaging velocimetry (Echo-PIV) analysis was performed for each configuration. The resulting hemodynamic force LV flow angle (φ) was analyzed, which ranges from 0^o (predominantly base-apex forces) to 90^o (predominantly transverse forces). Higher angles indicate more energy dissipation at lateral walls due to transverse flow; lower angles indicate healthier flow aligned with the longitudinal base-apex path of the pressure gradient.

RESULTS

Twelve patients (58% male, 17% ischemic, 32±7% ejection fraction, 165 ± 18 ms intrinsic QRSd) completed QRSd and Echo-PIV assessment. Relative to intrinsic rhythm, BiV and MPP with nominal AVD reduced QRSd by 10 ± 9% and 12 ± 9%, respectively. BiV+SyncAV and MPP+SyncAV further reduced QRSd by 19 ± 8%, (p < 0.05 vs. BiV with nominal AVD) and 23 ± 9% (p < 0.05 vs BiV+SyncAV), respectively. Echo-PIV showed similar sequential hemodynamic improvements. LV flow angular orientation during intrinsic activation (46 ± 3^o) reduced with BiV+SyncAV (37 ± 4^o, p < 0.05 vs intrinsic) and further with MPP+SyncAV (34 ± 4^o, p < 0.05 vs BiV+SyncAV).

CONCLUSION

These results suggest that SyncAV may improve electrical synchrony and influence LV flow patterns in patients suffering from heart failure compared to conventional CRT with a fixed AVD, with further improvement observed by combining with MPP.

Mixed impact of torsion on LV hemodynamics: A CFD study based on the Chimera technique

Image-based patient-specific Computational Fluid Dynamics (CFD) models of the Left Ventricle (LV) can be used to quantify hemodynamics-based biomarkers that can support the clinicians in the early diagnosis, follow-up and treatment planning of patients, beyond the capabilities of the current imaging modalities. We propose a workflow to build patient-specific CFD models of the LV with moving boundaries based on the Chimera technique to overcome the convergence issues previously encountered by means of the Arbitrarian Lagrangian Eulerian approach. The workflow was tested while investigating whether the torsional motion has an impact on LV fluid dynamics. Starting from 3D cine MRI scans of a healthy volunteer, six cardiac cycles were simulated in three CFD LV models: with no, physiological, and exaggerated torsion. The Chimera technique was robust in handling the impulsive motion of the LV endocardium, allowing to notice cycle-to-cycle variations in every simulated case. Torsion affected slightly velocity, vorticity, WSS. It did not affect energy loss and induced a double-sided effect in terms of residence time: the particles ejected in one beat decreased, whereas the motility of the particles remaining in the LV was affected only in the exaggerated torsion case, indicating that implementation of torsion can be discarded in case of physiological levels. Nonetheless, caution is warranted when interpreting these results given the absence of the mitral valve, the papillary muscles, and the trabeculae. The effects of the mitral valve will be evaluated within an Fluid Structure Interaction simulation framework as further development of the current model.

Analysis of mitral valve regurgitation by computational fluid dynamics

The clinical syndrome of mitral insufficiency is a common consequence of mitral valve (MV) prolapse, when the MV leaflets do not seal the closed orifice and blood regurgitates back to the atrium during ventricular contraction. There are different types of MV prolapse that may influence the degree of regurgitation also in relation to the left ventricle (LV) geometry. This study aims to provide some insight into the fluid dynamics of MV insufficiency in view of improving the different measurements available in the clinical setting. The analysis is performed by a computational fluid dynamics model coupled with an asymptotic model of the MV motion. The computational fluid dynamics solution uses the immersed boundary method that is efficiently integrated with clinical imaging technologies. Healthy and dilated LVs obtained by multislice cardiac MRI are combined with simplified models of healthy and pathological MVs deduced from computed tomography and 4D-transesophageal echocardiography. The results demonstrated the properties of false regurgitation, blood that did not cross the open MV orifice and returns into the atrium during the backward motion of the MV leaflets, whose entity should be accounted when evaluating small regurgitation. The regurgitating volume is found to be proportional to the effective orifice area, with the limited dependence of the LV geometry and type of prolapse. These affect the percentage of old blood returning to the atrium which may be associated with thrombogenic risk.

Prediction of Cardiac Resynchronization Therapy Response in Dilated Cardiomyopathy Using Vortex Flow Mapping on Cine Magnetic Resonance Imaging

Background: We investigated the association between left ventricle ejection fraction (LVEF) and vortex flow (VF), and whether cardiac resynchronization therapy (CRT) response can be predicted using VF mapping (VFM) in patients with dilated cardiomyopathy (DCM).

Methods and Results: Cardiac magnetic resonance imaging data for 20 patients with heart failure (HF) with LVEF ≥40% and 25 patients with DCM with LVEF <40%, scheduled for CRT, were retrospectively analyzed. The maximum VF (MVF) on short-axis, long-axis and 4-chamber LV cine imaging were calculated using VFM. Summed MVF was used as a representative value for each case and was significantly greater for patients with DCM than for patients with HF with LVEF ≥40% (25.2±19.2% vs. 12.1±15.4%, P<0.005). Summed MVF was significantly greater for CRT responders (n=12, 35.8±22.7%) than for non-responders (n=13, 15.8±8.7%, P=0.04) during the mean follow-up period of 38.4 months after CRT. Patients with summed MVF ≥31.3% had a significantly higher major adverse cardiac event-free rate than those with MVF <31.3% (log-rank=4.51, P<0.05).

Conclusions: On VFM analysis, LV VF interrupted efficient ejection in HF. Summed MVF can predict CRT response in DCM.

Impact of myocardial infarction on intraventricular vortex and flow energetics assessed using computational simulations

Flow energetics have been proposed as early indicators of progressive left ventricular (LV) functional impairment in patients with myocardial infarction (MI), but its correlation with individual MI parameters has not been fully explored. Using electro-fluid-structure interaction LV models, this study investigated the correlation between four MI parameters: infarct size, infarct multiplicity, regional enhancement of contractility at the viable myocardium area (RECVM), and LV mechanical dyssynchrony (LVMD) with intraventricular vortex and flow energetics. In LV with small infarcts, our results showed that infarct appearance amplified the energy dissipation index (DI), where substantial viscous energy loss was observed in areas with high flow velocity and near the infarct-vortex interface. The LV with small multiple infarcts and RECVM showed remarkable DI increment during systole and diastole. In correlation analysis, the systolic kinetic energy fluctuation index (E') was positively related to ejection fraction (EF) (R²  = 0.982) but negatively correlated with diastolic E' (R²  = 0.970). Diastolic E' was inversely correlated with vortex kinetic energy (R²  = 0.960) and vortex depth (R²  = 0.876). We showed an excessive systolic DI could differentiate infarcted LV with normal EF from healthy LV. Strong flow acceleration, LVMD, and vortex-infarct interactions were predominant factors that induced excessive DI in infarcted LVs. Instead of causing undesired flow turbulence, high systolic E' suggested the existence of energetic flow acceleration, while high diastolic E' implied an inefficient diastolic filling. Thus, systolic E' is not a suitable early indicator for progressive LV dysfunction in MI patients, while diastolic E' may be a useful index to indicate diastolic impairment in these patients.

Usefulness of Left Ventricular Vortex Flow Analysis for Predicting Clinical Outcomes in Patients with Chronic Heart Failure: A Quantitative Vorticity Imaging Study Using Contrast Echocardiography

The goal of the study described here was to evaluate whether left ventricular vortex flow parameters, as assessed by contrast echocardiography, enhance prediction of major adverse cardiac events (MACE) in patients with chronic heart failure and systolic dysfunction. A total of 75 patients with contrast echocardiography and systolic dysfunction (ejection fraction ≤45%) were prospectively enrolled and underwent vortex flow analysis with particle image velocimetry using contrast echocardiography. Vortex flow parameters, including kinetic energy fluctuation (KEF), were evaluated. Patients were followed up for a primary endpoint of MACE that comprised hospital admission for cardiovascular causes and cardiac deaths. Across a median 277-d follow-up, 29 patients (38.7%) experienced MACE. Among these, the incidence of diabetes and the E/e' ratio were significantly higher in patients with MACE than in those without, whereas the hemoglobin level and ejection fraction were significantly lower. KEF was significantly lower in patients with MACE. In the multivariate analysis, higher KEF was associated with a lower risk of MACE (hazard ratio = 0.18, 95% confidence interval: 0.04-0.97, p = 0.046). The addition of KEF to a model with conventional parameters (e.g., age, diabetes, ejection fraction and the E/e' ratio) significantly improved the model's discrimination. Elevations in the quantitative left ventricular vortex flow parameter, KEF, as determined by contrast echocardiography, are associated with a lower risk of MACE and improved functional status among patients with chronic heart failure.

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.

Synergy between Diastolic Mitral Valve Function and Left Ventricular Flow Aids in Valve Closure and Blood Transport during Systole

Highly resolved three-dimensional (3D) fluid structure interaction (FSI) simulation using patient-specific echocardiographic data can be a powerful tool for accurately and thoroughly elucidating the biomechanics of mitral valve (MV) function and left ventricular (LV) fluid dynamics. We developed and validated a strongly coupled FSI algorithm to fully characterize the LV flow field during diastolic MV opening under physiologic conditions. Our model revealed that distinct MV deformation and LV flow patterns developed during different diastolic stages. A vortex ring that strongly depended on MV deformation formed during early diastole. At peak E wave, the MV fully opened, with a local Reynolds number of ~5500, indicating that the flow was in the laminar-turbulent transitional regime. Our results  showed that during diastasis, the vortex structures caused the MV leaflets to converge, thus increasing mitral jet’s velocity. The vortex ring became asymmetrical, with the vortex structures on the anterior side being larger than on the posterior side. During the late diastolic stages, the flow structures advected toward the LV outflow tract, enhancing fluid transport to the aorta. This 3D-FSI study demonstrated the importance of leaflet dynamics, their effect on the vortex ring, and their influence on MV function and fluid transport within the LV during diastole.

Cardiac fluid dynamics meets deformation imaging

Cardiac function is about creating and sustaining blood in motion. This is achieved through a proper sequence of myocardial deformation whose final goal is that of creating flow. Deformation imaging provided valuable contributions to understanding cardiac mechanics; more recently, several studies evidenced the existence of an intimate relationship between cardiac function and intra-ventricular fluid dynamics. This paper summarizes the recent advances in cardiac flow evaluations, highlighting its relationship with heart wall mechanics assessed through the newest techniques of deformation imaging and finally providing an opinion of the most promising clinical perspectives of this emerging field. It will be shown how fluid dynamics can integrate volumetric and deformation assessments to provide a further level of knowledge of cardiac mechanics.

Optimized Time-Resolved Echo Particle Image Velocimetry– Particle Tracking Velocimetry Measurements Elucidate Blood Flow in Patients With Left Ventricular Thrombus

Contrast ultrasound is a widely used clinical tool to obtain real-time qualitative blood flow assessments in the heart, liver, etc. Echocardiographic particle image velocimetry (echo-PIV) is a technique for obtaining quantitative velocity maps from contrast ultrasound images. However, unlike optical particle image velocimetry (PIV), routine echo images are prone to nonuniform spatiotemporal variations in tracer distribution, making analysis difficult for standard PIV algorithms. This study introduces optimized procedures that integrate image enhancement, PIV, and particle tracking velocimetry (PTV) to obtain reliable time-resolved two-dimensional (2D) velocity distributions. During initial PIV analysis, multiple results are obtained by varying processing parameters. Optimization involving outlier removal and smoothing is used to select the correct vector. These results are used in a multiparameter PTV procedure. To demonstrate their clinical value, the procedures are implemented to obtain velocity and vorticity distributions over multiple cardiac cycles using images acquired from four left ventricular thrombus (LVT) patients. Phase-averaged data elucidate flow structure evolution over the cycle and are used to calculate penetration depth and strength of left ventricular (LV) vortices, as well as apical velocity induced by them. The present data are consistent with previous time-averaged results for the minimum vortex penetration depth associated with LVT occurrence. However, due to decay and fragmentation of LV vortices, as they migrate away from the mitral annulus, in two cases with high penetration, there is still poor washing near the resolved clot throughout the cycle. Hence, direct examination of entire flow evolution may be useful for assessing risk of LVT relapse before prescribing anticoagulants.

Evaluation of the anatomic and hemodynamic abnormalities in tricuspid atresia before and after surgery using computational fluid dynamics

Analysis of hemodynamics inside tricuspid atresia (TA) chamber is essential to the understanding of TA for optimal treatment. In this study, we introduced a combined computational fluid dynamics (CFD) to simulate blood flow in the left ventricle (LV) to study the diastolic flow changes in TA.Real-time 3-dimentional echocardiography loops (ECHO) were acquired in normal control group, in TA patients before surgery (pre-op group) and after surgery (post-op group). ECHO loops were reconstructed and simulated by CFD, the geometric, volumetric changes, and vortices in the LV were studies and compare among 3 groups.Compared with the control group, pre-op TA patients demonstrated significant LV remodeling, manifesting with smaller LV length, larger diameter, width and spherical index, as well as lager volumes; post-op TA group showed revisions in values of both geometric and volumetric measurements. CDF also demonstrated the abnormality of vortices in the pre-op TA patients and the alteration of existence and measurements of vortex in postoperation group.Echo-based CFD modeling can show the abnormality of TA in both LV geometric, volumetric measurements and intracardiac vortices; and CFD is capable to demonstrate the alterations of LV after Fontan and Glenn surgical procedure.

High-speed, high-frequency ultrasound, in utero vector-flow imaging of mouse embryos

Real-time imaging of the embryonic murine cardiovascular system is challenging due to the small size of the mouse embryo and rapid heart rate. High-frequency, linear-array ultrasound systems designed for small-animal imaging provide high-frame-rate and Doppler modes but are limited in regards to the field of view that can be imaged at fine-temporal and -spatial resolution. Here, a plane-wave imaging method was used to obtain high-speed image data from in utero mouse embryos and multi-angle, vector-flow algorithms were applied to the data to provide information on blood flow patterns in major organs. An 18-MHz linear array was used to acquire plane-wave data at absolute frame rates ≥10 kHz using a set of fixed transmission angles. After beamforming, vector-flow processing and image compounding, effective frame rates were on the order of 2 kHz. Data were acquired from the embryonic liver, heart and umbilical cord. Vector-flow results clearly revealed the complex nature of blood-flow patterns in the embryo with fine-temporal and -spatial resolution.

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.

Complex blood flow patterns in an idealized left ventricle: A numerical study

In this paper, we study the blood flow dynamics in a three-dimensional (3D) idealized left ventricle of the human heart whose deformation is driven by muscle contraction and relaxation in coordination with the action of the mitral and aortic valves. We propose a simplified but realistic mathematical treatment of the valves function based on mixed time-varying boundary conditions (BCs) for the Navier-Stokes equations modeling the flow. These switchings in time BCs, from natural to essential and vice versa, model either the open or the closed configurations of the valves. At the numerical level, these BCs are enforced by means of the extended Nitsche's method (Tagliabue et al., Int. J. Numer. Methods Fluids, 2017). Numerical results for the 3D idealized left ventricle obtained by means of Isogeometric Analysis are presented, discussed in terms of both instantaneous and phase-averaged quantities of interest and validated against those available in the literature, both experimental and computational. The complex blood flow patterns are analysed to describe the characteristic fluid properties, to show the transitional nature of the flow, and to highlight its main features inside the left ventricle. The sensitivity of the intraventricular flow patterns to the mitral valve properties is also investigated.

Left ventricular hemodynamic forces as a marker of mechanical dyssynchrony in heart failure patients with left bundle branch block

Left bundle branch block (LBBB) causes left ventricular (LV) dyssynchrony which is often associated with heart failure. A significant proportion of heart failure patients do not demonstrate clinical improvement despite cardiac resynchronization therapy (CRT). How LBBB-related effects on LV diastolic function may contribute to those therapeutic failures has not been clarified. We hypothesized that LV hemodynamic forces calculated from 4D flow MRI could serve as a marker of diastolic mechanical dyssynchrony in LBBB hearts. MRI data were acquired in heart failure patients with LBBB or matched patients without LBBB. LV pressure gradients were calculated from the Navier-Stokes equations. Integration of the pressure gradients over the LV volume rendered the hemodynamic forces. The findings demonstrate that the LV filling forces are more orthogonal to the main LV flow direction in heart failure patients with LBBB compared to those without LBBB during early but not late diastole. The greater the conduction abnormality the greater the discordance of LV filling force with the predominant LV flow direction (r² = 0.49). Such unique flow-specific measures of mechanical dyssynchrony may serve as an additional tool for considering the risks imposed by conduction abnormalities in heart failure patients and prove to be useful in predicting response to CRT.

A Newton-Krylov method with an approximate analytical Jacobian for implicit solution of Navier-Stokes equations on staggered overset-curvilinear grids with immersed boundaries

The explicit and semi-implicit schemes in flow simulations involving complex geometries and moving boundaries suffer from time-step size restriction and low convergence rates. Implicit schemes can be used to overcome these restrictions, but implementing them to solve the Navier-Stokes equations is not straightforward due to their non-linearity. Among the implicit schemes for nonlinear equations, Newton-based techniques are preferred over fixed-point techniques because of their high convergence rate but each Newton iteration is more expensive than a fixed-point iteration. Krylov subspace methods are one of the most advanced iterative methods that can be combined with Newton methods, i.e., Newton-Krylov Methods (NKMs) to solve non-linear systems of equations. The success of NKMs vastly depends on the scheme for forming the Jacobian, e.g., automatic differentiation is very expensive, and matrix-free methods without a preconditioner slow down as the mesh is refined. A novel, computationally inexpensive analytical Jacobian for NKM is developed to solve unsteady incompressible Navier-Stokes momentum equations on staggered overset-curvilinear grids with immersed boundaries. Moreover, the analytical Jacobian is used to form preconditioner for matrix-free method in order to improve its performance. The NKM with the analytical Jacobian was validated and verified against Taylor-Green vortex, inline oscillations of a cylinder in a fluid initially at rest, and pulsatile flow in a 90 degree bend. The capability of the method in handling complex geometries with multiple overset grids and immersed boundaries is shown by simulating an intracranial aneurysm. It was shown that the NKM with an analytical Jacobian is 1.17 to 14.77 times faster than the fixed-point Runge-Kutta method, and 1.74 to 152.3 times (excluding an intensively stretched grid) faster than automatic differentiation depending on the grid (size) and the flow problem. In addition, it was shown that using only the diagonal of the Jacobian further improves the performance by 42 - 74% compared to the full Jacobian. The NKM with an analytical Jacobian showed better performance than the fixed point Runge-Kutta because it converged with higher time steps and in approximately 30% less iterations even when the grid was stretched and the Reynold number was increased. In fact, stretching the grid decreased the performance of all methods, but the fixed-point Runge-Kutta performance decreased 4.57 and 2.26 times more than NKM with a diagonal Jacobian when the stretching factor was increased, respectively. The NKM with a diagonal analytical Jacobian and matrix-free method with an analytical preconditioner are the fastest methods and the superiority of one to another depends on the flow problem. Furthermore, the implemented methods are fully parallelized with parallel efficiency of 80-90% on the problems tested. The NKM with the analytical Jacobian can guide building preconditioners for other techniques to improve their performance in the future.

Numerical prediction of thrombus risk in an anatomically dilated left ventricle: the effect of inflow cannula designs

BACKGROUND

Implantation of a rotary blood pump (RBP) can cause non-physiological flow fields in the left ventricle (LV) which may trigger thrombosis. Different inflow cannula geometry can affect LV flow fields. The aim of this study was to determine the effect of inflow cannula geometry on intraventricular flow under full LV support in a patient specific model.

METHODS

Computed tomography angiography imaging of the LV was performed on a RBP candidate to develop a patient-specific model. Five inflow cannulae were evaluated, which were modelled on those used clinically or under development. The inflow cannulae are described as a crown like tip, thin walled tubular tip, large filleted tip, trumpet like tip and an inferiorly flared cannula. Placement of the inflow cannula was at the LV apex with the central axis intersecting the centre of the mitral valve. Full support was simulated by prescribing 5 l/min across the mitral valve. Thrombus risk was evaluated by identifying regions of stagnation. Rate of LV washout was assessed using a volume of fluid model. Relative haemolysis index and blood residence time was calculated using an Eulerian approach.

RESULTS

The inferiorly flared inflow cannula had the lowest thrombus risk due to low stagnation volumes. All cannulae had similar rates of LV washout and blood residence time. The crown like tip and thin walled tubular tip resulted in relatively higher blood damage indices within the LV.

CONCLUSION

Changes in intraventricular flow due to variances in cannula geometry resulted in different stagnation volumes. Cannula geometry does not appreciably affect LV washout rates and blood residence time. The patient specific, full support computational fluid dynamic model provided a repeatable platform to investigate the effects of inflow cannula geometry on intraventricular flow.

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.

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.

Automated Three-Dimensional Reconstruction of the Left Ventricle From Multiple-Axis Echocardiography

Two-dimensional echocardiography (echo) is the method of choice for noninvasive evaluation of the left ventricle (LV) function owing to its low cost, fast acquisition time, and high temporal resolution. However, it only provides the LV boundaries in discrete 2D planes, and the 3D LV geometry needs to be reconstructed from those planes to quantify LV wall motion, acceleration, and strain, or to carry out flow simulations. An automated method is developed for the reconstruction of the 3D LV endocardial surface using echo from a few standard cross sections, in contrast with the previous work that has used a series of 2D scans in a linear or rotational manner for 3D reconstruction. The concept is based on a generalized approach so that the number or type (long-axis (LA) or short-axis (SA)) of sectional data is not constrained. The location of the cross sections is optimized to minimize the difference between the reconstructed and measured cross sections, and the reconstructed LV surface is meshed in a standard format. Temporal smoothing is implemented to smooth the motion of the LV and the flow rate. This software tool can be used with existing clinical 2D echo systems to reconstruct the 3D LV geometry and motion to quantify the regional akinesis/dyskinesis, 3D strain, acceleration, and velocities, or to be used in ventricular flow simulations.

Hemodynamics in the Left Atrium and Its Effect on Ventricular Flow Patterns

In the present study, we investigate the hemodynamics inside left atrium (LA) and understand its impact on the development of ventricular flow patterns. We construct the heart model using dynamic-computed tomographic images and perform simulations using an immersed boundary method based flow solver. We show that the atrial hemodynamics is characterized by a circulatory flow generated by the left pulmonary veins (LPVs) and a direct stream from the right pulmonary veins (RPVs). The complex interaction of the vortex rings formed from each of the PVs leads to vortex breakup and annihilation, thereby producing a regularized flow at the mitral annulus. A comparison of the ventricular flow velocities between the physiological and a simplified pipe-based atrium model shows that the overall differences are limited to about 10% of the peak mitral flow velocity. The implications of this finding on the functional morphology of the left heart as well the computational and experimental modeling of ventricular hemodynamics are discussed.

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.

Assessment of viscous energy loss and the association with three-dimensional vortex ring formation in left ventricular inflow: In vivo evaluation using four-dimensional flow MRI

PURPOSE

To evaluate viscous energy loss and the association with three-dimensional (3D) vortex ring formation in left ventricular (LV) blood flow during diastolic filling.

THEORY AND METHODS

Thirty healthy volunteers were compared with 32 patients with corrected atrioventricular septal defect as unnatural mitral valve morphology and inflow are common in these patients. 4DFlow MRI was acquired from which 3D vortex ring formation was identified in LV blood flow at peak early (E)-filling and late (A)-filling and characterized by its presence/absence, orientation, and position from the lateral wall. Viscous energy loss was computed over E-filling, A-filling, and complete diastole using the Navier-Stokes energy equations.

RESULTS

Compared with healthy volunteers, viscous energy loss was significantly elevated in patients with disturbed vortex ring formation as characterized by a significantly inclined orientation and/or position closer to the lateral wall. Highest viscous energy loss was found in patients without a ring-shaped vortex during E-filling (on average more than double compared with patients with ring-shape vortex, P < 0.003). Altered A-filling vortex ring formation was associated with significant increase in total viscous energy loss over diastole even in the presence of normal E-filling vortex ring.

CONCLUSION

Altered vortex ring formation during LV filling is associated with increased viscous energy loss. Magn Reson Med 77:794-805, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Changes in electrical activation modify the orientation of left ventricular flow momentum: novel observations using echocardiographic particle image velocimetry

AIMS

Changes in electrical activation sequence are known to affect the timing of cardiac mechanical events. We aim to demonstrate that these also modify global properties of the intraventricular blood flow pattern. We also explore whether such global changes present a relationship with clinical outcome.

METHODS AND RESULTS

We investigated 30 heart failure patients followed up after cardiac resynchronization therapy (CRT). All subjects underwent echocardiography before implant and at follow-up after 6+ months. Left ventricular mechanics was investigated at follow-up during active CRT and was repeated after a temporary interruption <5 min later. Strain analysis, performed by speckle tracking, was used to assess the entity of contraction (global longitudinal strain) and its synchronicity (standard deviation of time to peak of radial strain). Intraventricular fluid dynamics, by echographic particle image velocimetry, was used to evaluate the directional distribution of global momentum associated with blood motion. The discontinuation of CRT pacing reflects into a reduction of deformation synchrony and into the deviation of blood flow momentum from the base-apex orientation with the development of transversal flow-mediated haemodynamic forces. The deviation of flow momentum presents a significant correlation with the degree of volumetric reduction after CRT.

CONCLUSION

Changes in electrical activation alter the orientation of blood flow momentum. The long-term CRT outcome correlates with the degree of re-alignment of haemodynamic forces. These preliminary results suggest that flow orientation could be used for optimizing the biventricular pacing setting. However, larger prospective studies are needed to confirm this hypothesis.

The Role of Shape and Heart Rate on the Performance of the Left Ventricle

The left ventricle function is to pump the oxygenated blood through the circulatory system. Ejection fraction is the main noninvasive parameter for detecting heart disease (healthy >55%), and it is thought to be the main parameter affecting efficiency. However, the effects of other parameters on efficiency have yet to be investigated. We investigate the effect of heart rate and left ventricle shape by carrying out 3D numerical simulations of a left ventricle at different heart rates and perturbed geometries under constant, normal ejection fraction. The simulation using the immersed boundary method provide the 3D flow and pressure fields, which enable direct calculation of a new hemodynamic efficiency (H-efficiency) parameter, which does not depend on any reference pressure. The H-efficiency is defined as the ratio of flux of kinetic energy (useful power) to the total cardiac power into the left ventricle control volume. Our simulations show that H-efficiency is not that sensitive to heart rate but is maximized at around normal heart rate (72 bpm). Nevertheless, it is more sensitive to the shape of the left ventricle, which affects the H-efficiency by as much as 15% under constant ejection fraction.

Asymptotic Model of Fluid-Tissue Interaction for Mitral Valve Dynamics

The vortex formation process inside the left ventricle is intrinsically connected to the dynamics of the mitral leaflets while they interact with the flow crossing the valve during diastole. The description of the dynamics of a natural mitral valve still represents a challenging issue, especially because its material properties are not measurable in vivo. Medical imaging can provide some indications about the geometry of the valve, but not about its mechanical properties. In this work, we introduce a parametric model of the mitral valve geometry, whose motion is described in the asymptotic limit under the assumption that it moves with the flow, without any additional resistance other than that given by its shape, and without the need to specify its material properties. The mitral valve model is coupled with a simple description of the left ventricle geometry, and their dynamics is solved numerically together with the equations ruling the blood flow. The intra-ventricular flow is analyzed in its relationship with the valvular motion. It is found that the initial valve opening anticipates the peak velocity of the Early filling wave with little influence of the specific geometry; while subsequent closure and re-opening are more dependent on the intraventricular vortex dynamics and thus on the leaflets' geometry itself. The limitations and potential applications of the proposed model are discussed.

Computational Study of Hemodynamic Effects of Abnormal E/A Ratio on Left Ventricular Filling

Three-dimensional numerical simulations are employed to investigate the hemodynamic effects of abnormal E/A ratios on left ventricular filling. The simulations are performed in a simplified geometric model of the left ventricle (LV) in conjunction with a specified endocardial motion. The model has been carefully designed to match the important geometric and flow parameters under the physiological conditions. A wide range of E/A ratios from 0 to infinity is employed with the aim to cover all the possible stages of left ventricle diastolic dysfunction (DD). The effects of abnormal E/A ratios on vortex dynamics, flow propagation velocity, energy consumption as well as flow transport and mixing are extensively discussed. Our results are able to confirm some common findings reported by the previous studies, and also uncover some interesting new features.

Intraventricular vortex properties in nonischemic dilated cardiomyopathy

Vortices may have a role in optimizing the mechanical efficiency and blood mixing of the left ventricle (LV). We aimed to characterize the size, position, circulation, and kinetic energy (KE) of LV main vortex cores in patients with nonischemic dilated cardiomyopathy (NIDCM) and analyze their physiological correlates. We used digital processing of color-Doppler images to study flow evolution in 61 patients with NIDCM and 61 age-matched control subjects. Vortex features showed a characteristic biphasic temporal course during diastole. Because late filling contributed significantly to flow entrainment, vortex KE reached its maximum at the time of the peak A wave, storing 26 ± 20% of total KE delivered by inflow (range: 1-74%). Patients with NIDCM showed larger and stronger vortices than control subjects (circulation: 0.008 ± 0.007 vs. 0.006 ± 0.005 m(2)/s, respectively, P = 0.02; KE: 7 ± 8 vs. 5 ± 5 mJ/m, P = 0.04), even when corrected for LV size. This helped confining the filling jet in the dilated ventricle. The vortex Reynolds number was also higher in the NIDCM group. By multivariate analysis, vortex KE was related to the KE generated by inflow and to chamber short-axis diameter. In 21 patients studied head to head, Doppler measurements of circulation and KE closely correlated with phase-contract magnetic resonance values (intraclass correlation coefficient = 0.82 and 0.76, respectively). Thus, the biphasic nature of filling determines normal vortex physiology. Vortex formation is exaggerated in patients with NIDCM due to chamber remodeling, and enlarged vortices are helpful for ameliorating convective pressure losses and facilitating transport. These findings can be accurately studied using ultrasound.