Image-Based Simulations Show Important Flow Fluctuations in a Normal Left Ventricle: What Could be the Implications?

Image-Based Computational Fluid Dynamics to Compare Two Repair Techniques for Mitral Valve Prolapse

OBJECTIVE : The treatment of mitral valve prolapse involves two distinct repair techniques: chordal replacement (Neochordae technique) and leaflet resection (Resection technique). However, there is still a debate in the literature about which is the optimal one. In this context, we performed an image-based computational fluid dynamic study to evaluate blood dynamics in the two surgical techniques. METHODS : We considered a healthy subject (H) and two patients (N and R) who underwent surgery for prolapse of the posterior leaflet and were operated with the Neochordae and Resection technique, respectively. Computational Fluid Dynamics (CFD) was employed with prescribed motion of the entire left heart coming from cine-MRI images, with a Large Eddy Simulation model to describe the transition to turbulence and a resistive method for managing valve dynamics. We created three different virtual scenarios where the operated mitral valves were inserted in the same left heart geometry of the healthy subject to study the differences attributed only to the two techniques. RESULTS : We compared the three scenarios by quantitatively analyzing ventricular velocity patterns and pressures, transition to turbulence, and the ventricle ability to prevent thrombi formation. From these results, we found that the operative techniques affected the ventricular blood dynamics in different ways, with variations attributed to the reduced mobility of the Resection posterior leaflet. Specifically, the Resection technique resulted in turbulent forces, related with the risk of hemolysis formation, up to 640 Pa, while the other two scenarios exhibited a maximum of 240 Pa. Moreover, in correspondence of the ventricular apex, the Resection technique reduced the areas with low velocity to 15%, whereas the healthy case and the Neochordae case maintained these areas at 30 and 48%, respectively. Our findings suggest that the Neochordae technique developed a more physiological flow with respect to the Resection technique. CONCLUSION: Resection technique gives rise to a different direction of the mitral jet during diastole increasing the ability to washout the ventricular apex preventing from thrombi formation, but at the same time it promotes turbulence formation that is associated with ventricular effort and risk of hemolysis.

Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.

Assessing turbulent effects in ascending aorta in presence of bicuspid aortic valve

Aortic valves with bicuspids have two rather than three leaflets, which is a congenital heart condition. About 0.5-2% of people have a bicuspid aortic valve. Blood flow through the aorta is commonly believed to be laminar, although aortic valve disorders can cause turbulent transitions. Understanding the impact of turbulence is crucial for foreseeing how the disease will progress. The study's objective was use large eddy simulation to provide a thorough analysis of the turbulence in bicuspid aortic valve dysfunction. Using a large eddy simulation, the blood flow patterns of the bicuspid and tricuspid aortic valves were compared, and significant discrepancies were found. The velocity field in flow in bicuspid configurations was asymmetrically distributed toward the ascending aorta. In tricuspid aortic valve (TAV) the flow, on the other hand, was symmetrical within the same aortic segment. Moreover, we looked into standard deviation, Q-criterion, viscosity ratio and wall shear stresses for each cases to understand transition to turbulence. Our findings indicate that in the bicuspid aortic valve (BAV) case, the fluid-dynamic abnormalities increase. The global turbulent kinetic energy and time-averaged wall shear stress for the TAV and BAV scenarios were also examined. We discovered that the global turbulent kinetic energy was higher in the BAV case compared to TAV, in addition to the increased wall shear stress induced by the BAV in the ascending aorta.

Generic framework for quantifying the influence of the mitral valve on ventricular blood flow

Blood flow within the left ventricle provides important information regarding cardiac function in health and disease. The mitral valve strongly influences the formation of flow structures and there exist various approaches for the representation of the valve in numerical models of left ventricular blood flow. However, a systematic comparison of the various mitral valve models is missing, making a priori decisions considering the overall model's context of use impossible. Within this study, a benchmark setup to compare the influence of mitral valve modeling strategies on intraventricular flow features was developed. Then, five mitral valve models of increasing complexity: no modeling, static wall, 2D and 3D porous medium with time-dependent porosity, and one-way fluid-structure interaction (FSI) were compared with each other. The flow features velocity, kinetic energy, transmitral pressure drop, vortex formation, flow asymmetry as well as computational cost and ease-of-implementation were evaluated. The one-way FSI approach provides the highest level of flow detail, which is accompanied by the highest numerical costs and challenges with the implementation. As an alternative, the porous medium approach with the expansion including time-dependent porosity provides good results with up to 10% deviations in the flow features (except the transmitral pressure drop) in comparison to the FSI model and only a fraction (11%) of numerical costs. However, jet propagation speed is highly underestimated by all alternative approaches to the FSI model. Taken together, our benchmark setup allows a quantitative comparison of various mitral valve modeling approaches and is provided to the scientific community for further testing and expansion.

Assessing a heterogeneous model for accounting for endovascular devices in hemodynamic simulations of cerebral aneurysms

The heterogeneous model developed by Berod et al [Int J Numer Method Biomed Eng 38, 2021] for representing the hemodynamic effects of endovascular prostheses is applied to a series of 10 patient specific cerebral aneurysms, 6 being treated by flow diverters, 4 being equipped with WEBs. Two markers correlated with the medical outcome of the treatment are used to assess the potential of the model, namely the saccular mean velocity and the inflow rate at the neck of the aneurysm. The comparison with the corresponding wire-resolved simulations is very favorable in both cases, and the model-based simulations also retrieve the jetting-type flows generated downstream of the struts. Noteworthy, the very same model was used for representing the flow diverters and the WEBs, showing the versatility and robustness of the heterogeneous modeling of the devices.

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.

Accelerated sequences of 4D flow MRI using GRAPPA and compressed sensing: A comparison against conventional MRI and computational fluid dynamics

PURPOSE

To evaluate hemodynamic markers obtained by accelerated GRAPPA (R = 2, 3, 4) and compressed sensing (R = 7.6) 4D flow MRI sequences under complex flow conditions.

METHODS

The accelerated 4D flow MRI scans were performed on a pulsatile flow phantom, along with a nonaccelerated fully sampled k-space acquisition. Computational fluid dynamics simulations based on the experimentally measured flow fields were conducted for additional comparison. Voxel-wise comparisons (Bland-Altman analysis, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>L</mml:mi> <mml:mn>2</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$$ {L}_2 $$</mml:annotation></mml:semantics> </mml:math> -norm metric), as well as nonderived quantities (velocity profiles, flow rates, and peak velocities), were used to compare the velocity fields obtained from the different modalities.

RESULTS

4D flow acquisitions and computational fluid dynamics depicted similar hemodynamic patterns. Voxel-wise comparisons between the MRI scans highlighted larger discrepancies at the voxels located near the phantom's boundary walls. A trend for all MR scans to overestimate velocity profiles and peak velocities as compared to computational fluid dynamics was noticed in regions associated with high velocity or acceleration. However, good agreement for the flow rates was observed, and eddy-current correction appeared essential for consistency of the flow rates measurements with respect to the principle of mass conservation.

CONCLUSION

GRAPPA (R = 2, 3) and highly accelerated compressed sensing showed good agreement with the fully sampled acquisition. Yet, all 4D flow MRI scans were hampered by artifacts inherent to the phase-contrast acquisition procedure. Computational fluid dynamics simulations are an interesting tool to assess these differences but are sensitive to modeling parameters.

Full-volume three-component intraventricular vector flow mapping by triplane color Doppler

Objective. Intraventricular vector flow mapping (iVFM) is a velocimetric technique for retrieving two-dimensional velocity vector fields of blood flow in the left ventricular cavity. This method is based on conventional color Doppler imaging, which makes iVFM compatible with the clinical setting. We have generalized the iVFM for a three-dimensional reconstruction (3D-iVFM). Approach. 3D-iVFM is able to recover three-component velocity vector fields in a full intraventricular volume by using a clinical echocardiographic triplane mode. The 3D-iVFM problem was written in the spherical (radial, polar, azimuthal) coordinate system associated to the six half-planes produced by the triplane mode. As with the 2D version, the method is based on the mass conservation, and free-slip boundary conditions on the endocardial wall. These mechanical constraints were imposed in a least-squares minimization problem that was solved through the method of Lagrange multipliers. We validated 3D-iVFM in silico in a patient-specific CFD (computational fluid dynamics) model of cardiac flow and tested its clinical feasibility in vivo in patients and in one volunteer. Main results. The radial and polar components of the velocity were recovered satisfactorily in the CFD setup (correlation coefficients, r  = 0.99 and 0.78). The azimuthal components were estimated with larger errors ( r  = 0.57) as only six samples were available in this direction. In both in silico and in vivo investigations, the dynamics of the intraventricular vortex that forms during diastole was deciphered by 3D-iVFM. In particular, the CFD results showed that the mean vorticity can be estimated accurately by 3D-iVFM. Significance. Our results tend to indicate that 3D-iVFM could provide full-volume echocardiographic information on left intraventricular hemodynamics from the clinical modality of triplane color Doppler.

A Pipeline for the Generation of Synthetic Cardiac Color Doppler

Color Doppler imaging (CDI) is the modality of choice for simultaneous visualization of myocardium and intracavitary flow over a wide scan area. This visualization modality is subject to several sources of error, the main ones being aliasing and clutter. Mitigation of these artifacts is a major concern for better analysis of intracardiac flow. One option to address these issues is through simulations. In this article, we present a numerical framework for generating clinical-like CDI. Synthetic blood vector fields were obtained from a patient-specific computational fluid dynamics CFD model. Realistic texture and clutter artifacts were simulated from real clinical ultrasound cineloops. We simulated several scenarios highlighting the effects of 1) flow acceleration; 2) wall clutter; and 3) transmit wavefronts, on Doppler velocities. As a comparison, an "ideal" color Doppler was also simulated, without these harmful effects. This synthetic dataset is made publicly available and can be used to evaluate the quality of Doppler estimation techniques. Besides, this approach can be seen as a first step toward the generation of comprehensive datasets for training neural networks to improve the quality of Doppler imaging.

A heterogeneous model of endovascular devices for the treatment of intracranial aneurysms

Numerical computations of hemodynamics inside intracranial aneurysms treated by endovascular braided devices such as flow-diverters contribute to understanding and improving such treatment procedures. Nevertheless, these simulations yield high computational and meshing costs due to the heterogeneity of length scales between the dense weave of the fine struts of the device and the arterial volume. Homogeneous strategies developed over the last decade to circumvent this issue substitute local dissipations due to the wires with a global effect in the form of a pressure-drop across the device surface. However, these methods cannot accurately reproduce the flow-patterns encountered near the struts, the latter strongly dictating the intra-saccular flow environment. In this work, a versatile theoretical framework which aims at correctly reproducing the local flow heterogeneities due to the wires while keeping memory consumption, meshing and computational times as low as possible is introduced. This model reproduces the drag forces exerted by the device struts onto the fluid, thus producing local and heterogeneous effects on the flow. Extensive validation for various flow and geometric configurations using an idealized device is performed. To further illustrate the method capabilities, a real patient-specific aneurysm endovascularly treated with a flow-diverter is used, enabling quantitative comparisons with classical approaches for both intra-saccular velocities and computational costs reduction. The proposed heterogeneous model endeavors to bridge the gap between computational fluid dynamics and clinical applications and ushers in a new era of numerical treatment planning with minimally costing computational tools.

Image-Based Computational Hemodynamics Analysis of Systolic Obstruction in Hypertrophic Cardiomyopathy

Hypertrophic Cardiomyopathy (HCM) is a pathological condition characterized by an abnormal thickening of the myocardium. When affecting the medio-basal portion of the septum, it is named Hypertrophic Obstructive Cardiomyopathy (HOCM) because it induces a flow obstruction in the left ventricular outflow tract. In any type of HCM, the myocardial function can become compromised, possibly resulting in cardiac death. In this study, we investigated with computational analysis the hemodynamics of patients with different types of HCM. The aim was quantifying the effects of this pathology on the intraventricular blood flow and pressure gradients, and providing information potentially useful to guide the indication and the modality of the surgical treatment (septal myectomy). We employed an image-based computational approach, integrating fluid dynamics simulations with geometric and functional data, reconstructed from standard cardiac cine-MRI acquisitions. We showed that with our approach we can better understand the patho-physiological behavior of intraventricular blood flow dynamics due to the abnormal morphological and functional aspect of the left ventricle. The main results of our investigation are: (a) a detailed patient-specific analysis of the blood velocity, pressure and stress distribution associated to HCM; (b) a computation-based classification of patients affected by HCM that can complement the current clinical guidelines for the diagnosis and treatment of HOCM.

Physics-constrained intraventricular vector flow mapping by color Doppler

Color Doppler by transthoracic echocardiography creates two-dimensional fan-shaped maps of blood velocities in the cardiac cavities. It is a one-component velocimetric technique since it only returns the velocity components parallel to the ultrasound beams. Intraventricular vector flow mapping (iVFM) is a method to recover the blood velocity vectors from the Doppler scalar fields in an echocardiographic three-chamber view. We improved ouriVFM numerical scheme by imposing physical constraints. TheiVFM consisted in minimizing regularized Doppler residuals subject to the condition that two fluid-dynamics constraints were satisfied, namely planar mass conservation, and free-slip boundary conditions. The optimization problem was solved by using the Lagrange multiplier method. A finite-difference discretization of the optimization problem, written in the polar coordinate system centered on the cardiac ultrasound probe, led to a sparse linear system. The single regularization parameter was determined automatically for non-supervision considerations. The physics-constrained method was validated using realistic intracardiac flow data from a patient-specific computational fluid dynamics (CFD) model. The numerical evaluations showed that theiVFM-derived velocity vectors were in very good agreement with the CFD-based original velocities, with relative errors ranged between 0.3% and 12%. We calculated two macroscopic measures of flow in the cardiac region of interest, the mean vorticity and mean stream function, and observed an excellent concordance between physics-constrainediVFM and CFD. The capability of physics-constrainediVFM was finally tested within vivocolor Doppler data acquired in patients routinely examined in the echocardiographic laboratory. The vortex that forms during the rapid filling was deciphered. The physics-constrainediVFM algorithm is ready for pilot clinical studies and is expected to have a significant clinical impact on the assessment of diastolic function.

A deep-learning approach for direct whole-heart mesh reconstruction

Automated construction of surface geometries of cardiac structures from volumetric medical images is important for a number of clinical applications. While deep-learning-based approaches have demonstrated promising reconstruction precision, these approaches have mostly focused on voxel-wise segmentation followed by surface reconstruction and post-processing techniques. However, such approaches suffer from a number of limitations including disconnected regions or incorrect surface topology due to erroneous segmentation and stair-case artifacts due to limited segmentation resolution. We propose a novel deep-learning-based approach that directly predicts whole heart surface meshes from volumetric CT and MR image data. Our approach leverages a graph convolutional neural network to predict deformation on mesh vertices from a pre-defined mesh template to reconstruct multiple anatomical structures in a 3D image volume. Our method demonstrated promising performance of generating whole heart reconstructions with as good or better accuracy than prior deep-learning-based methods on both CT and MR data. Furthermore, by deforming a template mesh, our method can generate whole heart geometries with better anatomical consistency and produce high-resolution geometries from lower resolution input image data. Our method was also able to produce temporally-consistent surface mesh predictions for heart motion from CT or MR cine sequences, and therefore can potentially be applied for efficiently constructing 4D whole heart dynamics. Our code and pre-trained networks are available at https://github.com/fkong7/MeshDeformNet.

On delayed transition to turbulence in an eccentric stenosis model for clean vs. noisy high-fidelity CFD

Recent comparisons between experiments and computational fluid dynamics (CFD) simulations of flow in the Food and Drug Administration (FDA) standardized nozzle geometry have highlighted the potential sensitivity of axisymmetric CFD models to small perturbations induced by mesh and inlet velocity, particularly for Reynolds numbers (Re) in the transitional regime. This evokes the classic experiment of Reynolds on transition to turbulence in a straight pipe, which can be delayed, apparently indefinitely, if special care is taken to control for external influences. Such idealized experiments are, however, extremely difficult to perform and, in the context of cardiovascular modeling, belie the "noise" inherent in typical experimental and physiological systems. Previous high-fidelity CFD of a canonical eccentric (i.e., non-axisymmetric) stenosis model showed transition occurring for steady flow at Re ~ 700-800, with modest delay caused by the introduction of shear-thinning rheology. On the other hand, recent experimental measurements of steady flowing blood and blood-mimicking fluids in this same stenosis model report transition for Re ~ 400-500. Taking a cue from the FDA nozzle controversy, the present study demonstrates that the addition of small-magnitude random noise at the inlet brings the eccentric-stenosis CFD results more in-line with experiments, and reveals a more gradual transition towards turbulence. This highlights that, even in non-axisymmetric idealized geometries, unnaturally "clean" high-fidelity CFD may impede not only good agreement with experiments, but also understanding of the onset and character of blood flow instabilities as they may exist, naturally, in the vasculature.

Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls

Congenital abnormalities in girls and women with Turner syndrome (TS), alongside an underlying predisposition to obesity and hypertension, contribute to an increased risk of cardiovascular disease and ultimately reduced life expectancy. We observe that children with TS present a greater variance in aortic arch morphology than their healthy counterparts, and hypothesize that their hemodynamics is also different. In this study, computational fluid dynamic (CFD) simulations were performed for four TS girls, and three age-matched healthy girls, using patient-specific inlet boundary conditions, obtained from phase-contrast MRI data. The visualization of multidirectional blood flow revealed an increase in vortical flow in the arch, supra-aortic vessels, and descending aorta, and a correlation between the presence of aortic abnormalities and disturbed flow. Compared to the relatively homogeneous pattern of time-averaged wall shear stress (TAWSS) on the healthy aortae, a highly heterogeneous distribution with elevated TAWSS values was observed in the TS geometries. Visualization of further shear stress parameters, such as oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS), revealed dissimilar heterogeneity in the oscillatory and multidirectional nature of the aortic flow. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% were obese or overweight), which is believed to accelerate the initiation and progression of endothelial dysfunction, these findings may be an indication of atherosclerotic disease manifesting earlier in life in TS patients. Age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk in TS children.

Numerical simulation of time-resolved 3D phase-contrast magnetic resonance imaging

A numerical approach is presented to efficiently simulate time-resolved 3D phase-contrast Magnetic resonance Imaging (or 4D Flow MRI) acquisitions under realistic flow conditions. The Navier-Stokes and Bloch equations are simultaneously solved with an Eulerian-Lagrangian formalism. A semi-analytic solution for the Bloch equations as well as a periodic particle seeding strategy are developed to reduce the computational cost. The velocity reconstruction pipeline is first validated by considering a Poiseuille flow configuration. The 4D Flow MRI simulation procedure is then applied to the flow within an in vitro flow phantom typical of the cardiovascular system. The simulated MR velocity images compare favorably to both the flow computed by solving the Navier-Stokes equations and experimental 4D Flow MRI measurements. A practical application is finally presented in which the MRI simulation framework is used to identify the origins of the MRI measurement errors.

Turbulence and turbulent flow structures in a ventricular assist device-A numerical study using the large-eddy simulation

Numerical flow simulations that analyze the turbulent flow characteristics within a turbopump are important for optimizing the efficiency of such machines. In the case of ventricular assist devices (VADs), turbulent flow characteristics must be also examined in order to improve hemocompatibility. Turbulence increases the shear stresses in the VAD flow, which can lead to an increased damage to the transported blood components. Therefore, an understanding of the turbulent flow patterns and their significance for the numerical blood damage prediction is particularly important for flow optimizations in VADs in order to identify and thus minimize flow regions where blood could be damaged due to high turbulent stresses. Nevertheless, the turbulence occurring in VADs and the local turbulent structures that lead to increased turbulent stresses have not yet been analyzed in detail in these machines. Therefore, this study aims to investigate the turbulence in an axial VAD in a comprehensive and double tracked way. First, the flow in an axial VAD was computed using the large-eddy simulation method, and it was verified that the majority of the turbulence was directly resolved by the simulation. Then, the turbulent flow state of the VAD was quantified globally. For this purpose, a self-designed evaluation method, the power loss analysis, was used. Subsequently, local flow regions and flow structures were identified where significant turbulent stresses prevail. It will be shown that the identified regions are universal and will also occur in other axial blood pumps as well, for example, in the HeartMate II.

Automating Model Generation for Image-Based Cardiac Flow Simulation

Computational fluid dynamics (CFD) modeling of left ventricle (LV) flow combined with patient medical imaging data has shown great potential in obtaining patient-specific hemodynamics information for functional assessment of the heart. A typical model construction pipeline usually starts with segmentation of the LV by manual delineation followed by mesh generation and registration techniques using separate software tools. However, such approaches usually require significant time and human efforts in the model generation process, limiting large-scale analysis. In this study, we propose an approach toward fully automating the model generation process for CFD simulation of LV flow to significantly reduce LV CFD model generation time. Our modeling framework leverages a novel combination of techniques including deep-learning based segmentation, geometry processing, and image registration to reliably reconstruct CFD-suitable LV models with little-to-no user intervention.1 We utilized an ensemble of two-dimensional (2D) convolutional neural networks (CNNs) for automatic segmentation of cardiac structures from three-dimensional (3D) patient images and our segmentation approach outperformed recent state-of-the-art segmentation techniques when evaluated on benchmark data containing both magnetic resonance (MR) and computed tomography(CT) cardiac scans. We demonstrate that through a combination of segmentation and geometry processing, we were able to robustly create CFD-suitable LV meshes from segmentations for 78 out of 80 test cases. Although the focus on this study is on image-to-mesh generation, we demonstrate the feasibility of this framework in supporting LV hemodynamics modeling by performing CFD simulations from two representative time-resolved patient-specific image datasets.

Contrast-Enhanced High-Frame-Rate Ultrasound Imaging of Flow Patterns in Cardiac Chambers and Deep Vessels

Cardiac function and vascular function are closely related to the flow of blood within. The flow velocities in these larger cavities easily reach 1 m/s, and generally complex spatiotemporal flow patterns are involved, especially in a non-physiologic state. Visualization of such flow patterns using ultrasound can be greatly enhanced by administration of contrast agents. Tracking the high-velocity complex flows is challenging with current clinical echographic tools, mostly because of limitations in signal-to-noise ratio; estimation of lateral velocities; and/or frame rate of the contrast-enhanced imaging mode. This review addresses the state of the art in 2-D high-frame-rate contrast-enhanced echography of ventricular and deep-vessel flow, from both technological and clinical perspectives. It concludes that current advanced ultrasound equipment is technologically ready for use in human contrast-enhanced studies, thus potentially leading to identification of the most clinically relevant flow parameters for quantifying cardiac and vascular function.

An image-based computational hemodynamics study of the Systolic Anterior Motion of the mitral valve

Systolic Anterior Motion (SAM) of the mitral valve - often associated with Hypertrophic Obstructive Cardiomyopathy (HOCM) - is a cardiac pathology in which a functional subaortic stenosis is induced during systole by the mitral leaflets partially obstructing the outflow tract of the left ventricle. Its assessment by diagnostic tests is often difficult, possibly underestimating its severity and thus increasing the risk of heart failure. In this paper, we propose a new computational pipeline, based on cardiac cine Magnetic Resonance Imaging (cine-MRI) data, for the assessment of SAM. The pipeline encompasses image processing of the left ventricle and the mitral valve, and numerical investigation of cardiac hemodynamics by means of Computational Fluid Dynamics (CFD) in a moving domain with image-based prescribed displacement. Patient-specific geometry and motion of the left ventricle are considered in view of an Arbitrary Lagrangian-Eulerian approach for CFD, while the reconstructed mitral valve is immersed in the computational domain by means of a resistive method. We assess clinically relevant flow and pressure indicators in a parametric study for different degrees of SAM severity, in order to provide a better quantitative evaluation of the pathological condition. Moreover, we provide specific indications for its possible surgical treatment, i.e. septal myectomy.

Augmented resistive immersed surfaces valve model for the simulation of cardiac hemodynamics with isovolumetric phases

In order to reduce the complexity of heart hemodynamics simulations, uncoupling approaches are often considered for the modeling of the immersed valves as an alternative to complex fluid-structure interaction (FSI) models. A possible shortcoming of these simplified approaches is the difficulty to correctly capture the pressure dynamics during the isovolumetric phases. In this work, we propose an enhanced resistive immersed surfaces (RIS) model of cardiac valves, which overcomes this issue. The benefits of the model are investigated and tested in blood flow simulations of the left heart where the physiological behavior of the intracavity pressure during the isovolumetric phases is recovered without using fully coupled fluid-structure models and without important alteration of the associated velocity field.

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.

Quantitative analysis of intraventricular flow-energetics and vortex in ischaemic hearts

OBJECTIVE

This study investigated the intraventricular flow dynamics in ischaemic heart disease patients.

PATIENTS AND METHODS

Fourteen patients with normal ejection fraction and 16 patients with reduced ejection fraction were compared with 20 healthy individuals. Phase-contrast MRI was used to assess intraventricular flow variables and speckle-tracking echocardiography to assess myocardial strain and left ventricular (LV) dyssynchrony. Infarct size was acquired using delayed-enhancement MRI.

RESULTS

The results obtained showed no significant differences in intraventricular flow variables between the healthy group and the patients with normal ejection fraction group, whereas considerable reductions in kinetic energy (KE) fluctuation index, E' (P<0.001) and vortex KE (P=0.003) were found in the patients with reduced ejection fraction group. In multivariate analysis, only vortex KE and infarct size were significantly related to LV ejection fraction (P<0.001); furthermore, vortex KE was correlated negatively with energy dissipation, energy dissipation index (r=-0.44, P=0.021).

CONCLUSION

This study highlights that flow energetic indices have limited applicability as early predictors of LV progressive dysfunction, whereas vortex KE could be an alternative to LV performance.

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.

Reconciling PC-MRI and CFD: An in-vitro study

Several well-resolved 4D Flow MRI acquisitions of an idealized rigid flow phantom featuring an aneurysm, a curved channel as well as a bifurcation were performed under pulsatile regime. The resulting hemodynamics were processed to remove MRI artifacts. Subsequently, they were compared with CFD predictions computed on the same flow domain, using an in-house high-order low dissipative flow solver. Results show that reaching a good agreement is not straightforward but requires proper treatments of both techniques. Several sources of discrepancies are highlighted and their impact on the final correlation evaluated. While a very poor correlation (r²  = 0.63) is found in the entire domain between raw MRI and CFD data, correlation as high as r²  = 0.97 is found when artifacts are removed by post-processing the MR data and down sampling the CFD results to match the MRI spatial and temporal resolutions. This work demonstrates that, in a well-controlled environment, both PC-MRI and CFD might bring reliable and correlated flow quantities when a proper methodology to reduce the errors is followed.

High Frame Rate Ultrasound Particle Image Velocimetry for Estimating High Velocity Flow Patterns in the Left Ventricle

Echocardiographic determination of multicomponent blood flow dynamics in the left ventricle remains a challenge. In this paper, we compare contrast enhanced, high frame rate (HFR) (1000 frames/s) echo-particle image velocimetry (ePIV) against optical particle image velocimetry (oPIV, gold standard), in a realistic left ventricular (LV) phantom. We find that ePIV compares well to oPIV, even for the high velocity inflow jet (normalized RMSE = 9% ± 1%). In addition, we perform the method of proper orthogonal decomposition, to better qualify and quantify the differences between the two modalities. We show that ePIV and oPIV resolve very similar flow structures, especially for the lowest order mode with a cosine similarity index of 86%. The coarser resolution of ePIV does result in increased variance and blurring of smaller flow structures when compared to oPIV. However, both modalities are in good agreement with each other for the modes that constitute the bulk of the kinetic energy. We conclude that HFR ePIV can accurately estimate the high velocity diastolic inflow jet and the high energy flow structures in an LV setting.

Impact of Pulmonary Venous Inflow on Cardiac Flow Simulations: Comparison with In Vivo 4D Flow MRI

Blood flow simulations are making their way into the clinic, and much attention is given to estimation of fractional flow reserve in coronary arteries. Intracardiac blood flow simulations also show promising results, and here the flow field is expected to depend on the pulmonary venous (PV) flow rates. In the absence of in vivo measurements, the distribution of the flow from the individual PVs is often unknown and typically assumed. Here, we performed intracardiac blood flow simulations based on time-resolved computed tomography on three patients, and investigated the effect of the distribution of PV flow rate on the flow field in the left atrium and ventricle. A design-of-experiment approach was used, where PV flow rates were varied in a systematic manner. In total 20 different simulations were performed per patient, and compared to in vivo 4D flow MRI measurements. Results were quantified by kinetic energy, mitral valve velocity profiles and root-mean-square errors of velocity. While large differences in atrial flow were found for varying PV inflow distributions, the effect on ventricular flow was negligible, due to a regularizing effect by mitral valve. Equal flow rate through all PVs most closely resembled in vivo measurements and is recommended in the absence of a priori knowledge.

The Influence of Rotary Blood Pump Speed Modulation on the Risk of Intraventricular Thrombosis

Rotary left ventricular assist devices (LVADs) are commonly operated at a constant speed, attenuating blood flow pulsatility. Speed modulation of rotary LVADs has been demonstrated to improve vascular pulsatility and pump washout. The effect of LVAD speed modulation on intraventricular flow dynamics is not well understood, which may have an influence on thromboembolic events. This study aimed to numerically evaluate intraventricular flow characteristics with a speed modulated LVAD. A severely dilated anatomical left ventricle was supported by a HeartWare HVAD in a three-dimensional multiscale computational fluid dynamics model. Three LVAD operating scenarios were evaluated: constant speed and sinusoidal co- and counter-pulsation. In all operating scenarios, the mean pump speed was set to restore the cardiac output to 5.0 L/min. Co- and counter-pulsation was speed modulated with an amplitude of 750 rpm. The risk of thrombosis was evaluated based on blood residence time, ventricular washout, kinetic energy densities, and a pulsatility index map. Blood residence time for co-pulsation was on average 1.8 and 3.7% lower than constant speed and counter-pulsation mode, respectively. After introducing fresh blood to displace preexisting blood for 10 cardiac cycles, co-pulsation had 1.5% less old blood in comparison to counter-pulsation. Apical energy densities were 84 and 27% higher for co-pulsation in comparison to counter-pulsation and constant speed mode, respectively. Co-pulsation had an increased pulsatility index around the left ventricular outflow tract and mid-ventricle. Improved flow dynamics with co-pulsation was caused by increased E-wave velocities which minimized blood stasis. In the studied scenario and from the perspective of intraventricular flow dynamics, co-pulsation of rotary LVADs could minimize the risk of intraventricular thrombosis.

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.

Left Ventricular Trabeculations Decrease the Wall Shear Stress and Increase the Intra-Ventricular Pressure Drop in CFD Simulations

The aim of the present study is to characterize the hemodynamics of left ventricular (LV) geometries to examine the impact of trabeculae and papillary muscles (PMs) on blood flow using high performance computing (HPC). Five pairs of detailed and smoothed LV endocardium models were reconstructed from high-resolution magnetic resonance images (MRI) of ex-vivo human hearts. The detailed model of one LV pair is characterized only by the PMs and few big trabeculae, to represent state of art level of endocardial detail. The other four detailed models obtained include instead endocardial structures measuring ≥1 mm² in cross-sectional area. The geometrical characterizations were done using computational fluid dynamics (CFD) simulations with rigid walls and both constant and transient flow inputs on the detailed and smoothed models for comparison. These simulations do not represent a clinical or physiological scenario, but a characterization of the interaction of endocardial structures with blood flow. Steady flow simulations were employed to quantify the pressure drop between the inlet and the outlet of the LVs and the wall shear stress (WSS). Coherent structures were analyzed using the Q-criterion for both constant and transient flow inputs. Our results show that trabeculae and PMs increase the intra-ventricular pressure drop, reduce the WSS and disrupt the dominant single vortex, usually present in the smoothed-endocardium models, generating secondary small vortices. Given that obtaining high resolution anatomical detail is challenging in-vivo, we propose that the effect of trabeculations can be incorporated into smoothed ventricular geometries by adding a porous layer along the LV endocardial wall. Results show that a porous layer of a thickness of 1.2·10⁻² m with a porosity of 20 kg/m² on the smoothed-endocardium ventricle models approximates the pressure drops, vorticities and WSS observed in the detailed models.

Fluid-structure interaction of a pulsatile flow with an aortic valve model: A combined experimental and numerical study

The complex fluid-structure interaction problem associated with the flow of blood through a heart valve with flexible leaflets is investigated both experimentally and numerically. In the experimental test rig, a pulse duplicator generates a pulsatile flow through a biomimetic rigid aortic root where a model of aortic valve with polymer flexible leaflets is implanted. High-speed recordings of the leaflets motion and particle image velocimetry measurements were performed together to investigate the valve kinematics and the dynamics of the flow. Large eddy simulations of the same configuration, based on a variant of the immersed boundary method, are also presented. A massively parallel unstructured finite-volume flow solver is coupled with a finite-element solid mechanics solver to predict the fluid-structure interaction between the unsteady flow and the valve. Detailed analysis of the dynamics of opening and closure of the valve are conducted, showing a good quantitative agreement between the experiment and the simulation regarding the global behavior, in spite of some differences regarding the individual dynamics of the valve leaflets. A multicycle analysis (over more than 20 cycles) enables to characterize the generation of turbulence downstream of the valve, showing similar flow features between the experiment and the simulation. The flow transitions to turbulence after peak systole, when the flow starts to decelerate. Fluctuations are observed in the wake of the valve, with maximum amplitude observed at the commissure side of the aorta. Overall, a very promising experiment-vs-simulation comparison is shown, demonstrating the potential of the numerical method.

Ventricular flow dynamics with varying LVAD inflow cannula lengths: In-silico evaluation in a multiscale model

Left ventricular assist devices are associated with thromboembolic events, which are potentially caused by altered intraventricular flow. Due to patient variability, differences in apical wall thickness affects cannula insertion lengths, potentially promoting unfavourable intraventricular flow patterns which are thought to be correlated to the risk of thrombosis. This study aimed to present a 3D multiscale computational fluid dynamic model of the left ventricle (LV) developed using a commercial software, Ansys, and evaluate the risk of thrombosis with varying inflow cannula insertion lengths in a severely dilated LV. Based on a HeartWare HVAD inflow cannula, insertion lengths of 5, 19, 24 and 50 mm represented cases of apical hypertrophy, typical ranges of apical thicknesses and an experimental length, respectively. The risk of thrombosis was evaluated based on blood washout, residence time, instantaneous blood stagnation and a pulsatility index. By introducing fresh blood to displace pre-existing blood in the LV, after 5 cardiac cycles, 46.7%, 45.7%, 45.1% and 41.8% of pre-existing blood remained for insertion lengths of 5, 19, 24 and 50 mm, respectively. Compared to the 50 mm insertion, blood residence time was at least 9%, 7% and 6% higher with the 5, 19 and 24 mm insertion lengths, respectively. No instantaneous stagnation at the apex was observed directly after the E-wave. Pulsatility indices adjacent to the cannula increased with shorter insertion lengths. For the specific scenario studied, a longer insertion length, relative to LV size, may be advantageous to minimise thrombosis by increasing LV washout and reducing blood residence time.

Impact of LVAD Implantation Site on Ventricular Blood Stagnation

Treatment of end-stage heart failure includes cardiac transplantation or ventricular assist device (VAD) therapy. Although increasingly prevalent, current VAD therapy has inherent complications, including thrombosis. Studies have demonstrated that VAD implantation alters intracardiac blood flow, creating areas of stagnation that predispose to thrombus formation. Two potential surgical configurations exist for VAD implantation: through the apical or diaphragmatic surfaces of the heart. We hypothesized that diaphragmatic implantation causes more stagnation than apical implantation. We also hypothesized that intermittent aortic valve (AV) opening reduces stagnation of blood inside the left ventricle (LV) when compared with a closed AV. To test these hypotheses, a human LV geometry was recreated in silico and a VAD inflow cannula was virtually implanted in each configuration. A computational indicator-dilution study was conducted where "virtually dyed blood" was washed out of the LV by injecting blood with no dye. Simulations demonstrated a substantial reduction in stagnation with intermittent AV opening. In addition, virtual dye was cleared slightly faster in the apical configuration. Simulations from our study demonstrate the clinical importance of VAD management to allow intermittent opening of the AV to prevent subvalvular stagnation, and also suggests that apical configuration might be more hemodynamically favorable.

Introducing the pro-coagulant contact system in the numerical assessment of device-related thrombosis

Thrombosis is a major concern in blood-coated medical devices. Contact activation, which is the initial part of the coagulation cascade in device-related thrombosis, is not considered in current thrombus formation models. In the present study, pro-coagulant reactions including the contact activation system are coupled with a fluid solver in order to evaluate the potential of the contact system to initiate thrombin production. The biochemical/fluid model is applied to a backward-facing step configuration, a flow configuration that frequently appears in medical devices. In contrast to the in vivo thrombosis models in which a specific thrombotic zone (injury region) is set a priori by the user to initiate the coagulation reaction, a reactive surface boundary condition is applied to the whole device wall. Simulation results show large thrombin concentration in regions related to recirculation zones without the need of an a priori knowledge of the thrombus location. The numerical results align well with the regions prone to thrombosis observed in experimental results reported in the literature. This approach could complement thrombus formation models that take into account platelet activity and thrombus growth to optimize a wide range of medical devices.

Fully-coupled fluid-structure interaction simulation of the aortic and mitral valves in a realistic 3D left ventricle model

In this study, we present a fully-coupled fluid-structure interaction (FSI) framework that combines smoothed particle hydrodynamics (SPH) and nonlinear finite element (FE) method to investigate the coupled aortic and mitral valves structural response and the bulk intraventricular hemodynamics in a realistic left ventricle (LV) model during the entire cardiac cycle. The FSI model incorporates valve structures that consider native asymmetric leaflet geometries, anisotropic hyperelastic material models and human material properties. Comparison of FSI results with subject-specific echocardiography data demonstrates that the SPH-FE approach is able to quantitatively predict the opening and closing times of the valves, the mitral leaflet opening and closing angles, and the large-scale intraventricular flow phenomena with a reasonable agreement. Moreover, comparison of FSI results with a LV model without valves reveals substantial differences in the flow field. Peak systolic velocities obtained from the FSI model and the LV model without valves are 2.56 m/s and 1.16 m/s, respectively, compared to the Doppler echo data of 2.17 m/s. The proposed SPH-FE FSI framework represents a further step towards modeling patient-specific coupled LV-valve dynamics, and has the potential to improve our understanding of cardiovascular physiology and to support professionals in clinical decision-making.

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.