Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

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.     

Flow-pressure drop measurement and calculation in a tapered femoral artery of a dog

This study describes the in vivo measurement of pressure drop and flow during the cardiac cycle in the femoral artery of a dog, and the computer simulation of the experiment based on the use of the measured flow, vessel dimensions and blood viscosity. In view of the experimental uncertainty in obtaining the accurate velocity profile at the wall region, the velocity pulse at the center was measured and numerical calculations were performed for the center Une instantaneous velocity and within the two limits of spatial distribution of inlet flow conditions: uniform and parabolic. Temporal and spatial variations of flow parameters, i.e., velocity profile, shear rate, non-Newtonian viscosity, wall shear stress, and pressure drop were calculated. There existed both positive and negative shear rates during a pulse cycle, i.e., the arterial wall experiences zero shear three times during a cardiac cycle. For the parabolic inlet condition, the taper of the artery not only increased the magnitude of the positive and negative shear rates, but caused a steep gradient in shear rate, a phenomenon which in turn affects wall shear stress and pressure. In contrast, for the uniform inlet condition, the flow through the tapered artery was predominantly the developing type, which resulted in reduction in magnitude of wall shear rate along the axial direction.     

The effect of vortex formation on left ventricular filling and mitral valve efficiency

A new mechanism for quantifying the filling energetics in the left ventricle (LV) and past mechanical heart valves (MHV) is identified and presented. This mechanism is attributed to vortex formation dynamics past MHV leaflets. Recent studies support the conjecture that the natural healthy left ventricle (LV) performs in an optimum, energy-preserving manner by redirecting the flow with high efficiency. Yet to date, no quantitative proof has been presented. The present work provides quantitative results and validation of a theory based on the dynamics of vortex ring formation, which is governed by a critical formation number (FN) that corresponds to the dimensionless time at which the vortex ring has reached its maximum circulation content, in support of this hypothesis. Herein, several parameters (vortex ring circulation, vortex ring energy, critical FN, hydrodynamic efficiencies, vortex ring propagation speed) have been quantified and presented as a means of bridging the physics of vortex formation in the LV. In fact, the diastolic hydrodynamic efficiencies were found to be 60, 41, and 29%, respectively, for the porcine, anti-anatomical, and anatomical valve configurations. This assessment provides quantitative proof of vortex formation, which is dependent of valve design and orientation, being an important flow characteristic and associated to LV energetics. Time resolved digital particle image velocimetry with kilohertz sampling rate was used to study the ejection of fluid into the LV and resolve the spatiotemporal evolution of the flow. The clinical significance of this study is quantifying vortex formation and the critical FN that can potentially serve as a parameter to quantify the LV filling process and the performance of heart valves.     

PIV-measured versus CFD-predicted flow dynamics in anatomically realistic cerebral aneurysm models

Computational fluid dynamics (CFD) modeling of nominally patient-specific cerebral aneurysms is increasingly being used as a research tool to further understand the development, prognosis, and treatment of brain aneurysms. We have previously developed virtual angiography to indirectly validate CFD-predicted gross flow dynamics against the routinely acquired digital subtraction angiograms. Toward a more direct validation, here we compare detailed, CFD-predicted velocity fields against those measured using particle imaging velocimetry (PIV). Two anatomically realistic flow-through phantoms, one a giant internal carotid artery (ICA) aneurysm and the other a basilar artery (BA) tip aneurysm, were constructed of a clear silicone elastomer. The phantoms were placed within a computer-controlled flow loop, programed with representative flow rate waveforms. PIV images were collected on several anterior-posterior (AP) and lateral (LAT) planes. CFD simulations were then carried out using a well-validated, in-house solver, based on micro-CT reconstructions of the geometries of the flow-through phantoms and inlet/outlet boundary conditions derived from flow rates measured during the PIV experiments. PIV and CFD results from the central AP plane of the ICA aneurysm showed a large stable vortex throughout the cardiac cycle. Complex vortex dynamics, captured by PIV and CFD, persisted throughout the cardiac cycle on the central LAT plane. Velocity vector fields showed good overall agreement. For the BA, aneurysm agreement was more compelling, with both PIV and CFD similarly resolving the dynamics of counter-rotating vortices on both AP and LAT planes. Despite the imposition of periodic flow boundary conditions for the CFD simulations, cycle-to-cycle fluctuations were evident in the BA aneurysm simulations, which agreed well, in terms of both amplitudes and spatial distributions, with cycle-to-cycle fluctuations measured by PIV in the same geometry. The overall good agreement between PIV and CFD suggests that CFD can reliably predict the details of the intra-aneurysmal flow dynamics observed in anatomically realistic in vitro models. Nevertheless, given the various modeling assumptions, this does not prove that they are mimicking the actual in vivo hemodynamics, and so validations against in vivo data are encouraged whenever possible.     

Patient-specific CFD models for intraventricular flow analysis from 3D ultrasound imaging: Comparison of three clinical cases

BackgroundAs the intracardiac flow field is affected by changes in shape and motility of the heart, intraventricular flow features can provide diagnostic indications. Ventricular flow patterns differ depending on the cardiac condition and the exploration of different clinical cases can provide insights into how flow fields alter in different pathologies.MethodsIn this study, we applied a patient-specific computational fluid dynamics model of the left ventricle and mitral valve, with prescribed moving boundaries based on transesophageal ultrasound images for three cardiac pathologies, to verify the abnormal flow patterns in impaired hearts. One case (P1) had normal ejection fraction but low stroke volume and cardiac output, P2 showed low stroke volume and reduced ejection fraction, P3 had a dilated ventricle and reduced ejection fraction.ResultsThe shape of the ventricle and mitral valve, together with the pathology influence the flow field in the left ventricle, leading to distinct flow features. Of particular interest is the pattern of the vortex formation and evolution, influenced by the valvular orifice and the ventricular shape. The base-to-apex pressure difference of maximum 2 mmHg is consistent with reported data.ConclusionWe used a CFD model with prescribed boundary motion to describe the intraventricular flow field in three patients with impaired diastolic function. The calculated intraventricular flow dynamics are consistent with the diagnostic patient records and highlight the differences between the different cases. The integration of clinical images and computational techniques, therefore, allows for a deeper investigation intraventricular hemodynamics in patho-physiology.     

An approach to the simulation of fluid-structure interaction in the aortic valve

A pair of finite element models has been employed to study the interaction of blood flow with the operation of the aortic valve. A three-dimensional model of the left ventricle with applied wall displacements has been used to generate data for the spatially and time-varying blood velocity profile across the aortic aperture. These data have been used as the inlet loading conditions in a three-dimensional model of the aortic valve and its surrounding structures. Both models involve fluid-structure interaction and simulate the cardiac cycle as a dynamic event. Confidence in the models was obtained by comparison with data obtained in a pulse duplicator. The results show a circulatory flow being generated in the ventricle which produces a substantially axial flow through the aortic aperture. The aortic valve behaves in an essentially symmetric way under the action of this flow, so that the pressure difference across the leaflets is approximately uniform. This work supports the use of spatially uniform but temporally variable pressure distributions across the leaflets in dry or structural models of aortic valves. The study is a major advance through its use of truly three-dimensional geometry, spatially non-uniform loading conditions for the valve leaflets and the successful modelling of progressive contact of the leaflets in a fluid environment.     

Fluid dynamics of the left ventricular filling in dilated cardiomyopathy

Modifications in diastolic function occur in a broad range of cardiovascular diseases and there is an increasing evidence that abnormalities in left ventricular function may contribute significantly to the symptomatology. The flow inside the left ventricle during the diastole is here investigated by numerical solution of the Navier-Stokes equations under the axisymmetric assumption. The equation are written in a body-fitted, moving prolate spheroid, system of coordinates and solved using a fractional step method. The system is forced by a given volume time-law derived from clinical data, and varying the two-degrees-of-freedom ventricle geometry on the basis of a simple model. The solution under healthy conditions is analysed in terms of vorticity dynamics, showing that the flow field is characterised by the presence of a vortex wake; it is attached to the mitral valve during the accelerating phase of the E-wave, and it detaches and translate towards the ventricle apex afterwards. The flow evolution is discussed, results are also reported as an M-mode representation of colour-coded Doppler velocity maps. In the presence of ventricle dilatation the mitral jet extends farther inside the ventricle, propagation velocity decreases, and the fluid stagnates longer at the apex.     

Computational analysis of blood flow in an integrated model of the left ventricle and the aorta

To study the effects of intraventricular flow dynamics on the aortic flow, we created an integrated model of the left ventricle and aorta and conducted a computer simulation of diastolic and systolic blood flow within this model. The results demonstrated that the velocity profile at the aortic annulus changed dynamically, and was influenced by the intraventricular flow dynamics. The profile was almost flat in early systole but became nonuniform as systole progressed, and was skewed toward the posterior side in midsystole and toward the anterior side in later systole. At a distance from the aortic annulus, a different velocity profile was induced by the twisting and torsion of the aorta. In the ascending aorta, the fastest flow was initially located in the posteromedial sector, and it moved to the posterior section along the circumference as systole progressed. The nonuniformity of the aortic inflow gave rise to a complex wall shear stress (WSS) distribution in the aorta. A comparison of the WSS distribution obtained in this integrated analysis with that obtained in flow calculations using an isolated aorta model with Poiseuille and flat inlet conditions showed that intraventricular flow affected the WSS distribution in the ascending aorta. These results address the importance of an integrated analysis of flow in the left ventricle and aorta.     

Model and influence of mitral valve opening during the left ventricular filling

The flow inside a model left ventricle during filling (diastole) is simulated by the numerical solution of the equations of motion under the axisymmetric approximation. The left ventricle is taken with a truncated ellipsoid geometry, and a simple conceptual model is introduced to simulate the presence of the moving mitral valve. A relevant role during the left ventricle diastolic flow, as already discussed by other authors, is played by the travelling vortex wake that is formed from the transmitral jet during the early filling acceleration phase. The presence of a moving valve is found to produce a non-simultaneous spatial development of the entering bulk flow and a slightly more complex vortex wake structure; the results are discussed in comparison with fixed valve ones. They are analysed also in terms of M-mode representation suggesting a physical interpretation of the pattern detected in the clinical measurements that extends the one given previously on the basis of fixed valve models.     

Jet collisions and vortex reversal in the human left ventricle

Unnatural dynamics of the notorious vortex in the left ventricle is often associated with cardiac disease. Understanding how different cardiac diseases alter the flow physics in the left ventricle may therefore provide a powerful tool for disease detection. In this work, the fluid dynamics in the left ventricle subject to different severities of aortic regurgitation is experimentally investigated by performing time-resolved particle image velocimetry in a left heart duplicator. Diastolic vortex reversal was observed in the left ventricle accompanied by an increase in viscous energy dissipation. Vortex dynamics and energy dissipation may provide useful insights on sub-optimal flow patterns in the left ventricle.     

Intraventricular flow patterns and stasis in the LVAD-assisted heart

Left ventricular assist device (LVAD) support disrupts the natural blood flow path through the heart, introducing flow patterns associated with thrombosis, especially in the presence of medical devices. The aim of this study was to quantitatively evaluate the flow patterns in the left ventricle (LV) of the LVAD-assisted heart, with a focus on alterations in vortex development and stasis. Particle image velocimetry of a LVAD-supported LV model was performed in a mock circulatory loop. In the Pre-LVAD flow condition, a vortex ring initiating from the LV base migrated toward the apex during diastole and remained in the LV by the end of ejection. During LVAD support, vortex formation was relatively unchanged although vortex circulation and kinetic energy increased with LVAD speed, particularly in systole. However, as pulsatility decreased and aortic valve opening ceased, a region of fluid stasis formed near the left ventricular outflow tract. These findings suggest that LVAD support does not substantially alter vortex dynamics unless cardiac function is minimal. The altered blood flow introduced by the LVAD results in stasis adjacent to the LV outflow tract, which increases the risk of thrombus formation in the heart.     

Influence of dilated cardiomyopathy and a left ventricular assist device on vortex dynamics in the left ventricle

Together with new developments in mechanical cardiac support, the analysis of vortex dynamics in the left ventricle has become an increasingly important topic in literature. The aim of this study was to develop a method to investigate the influence of a left ventricular assist device (LVAD) on vortex dynamics in a failing ventricle. An axisymmetric fluid dynamics model of the left ventricle was developed and coupled to a lumped parameter model of the complete circulation. Simulations were performed for healthy conditions and dilated cardiomyopathy (DCM). Vortex structures in these simulations were analysed by means of automated detection. Results show that the strength of the leading vortex ring is lower in a DCM ventricle than in a healthy ventricle. The LVAD further influences the maximum strength of the vortex and also causes the vortex to disappear earlier in time with increasing LVAD flows. Understanding these phenomena by means of the method proposed in this study will contribute to enhanced diagnostics and monitoring during cardiac support.     

Influence of dilated cardiomyopathy and a left ventricular assist device on vortex dynamics in the left ventricle

Together with new developments in mechanical cardiac support, the analysis of vortex dynamics in the left ventricle has become an increasingly important topic in literature. The aim of this study was to develop a method to investigate the influence of a left ventricular assist device (LVAD) on vortex dynamics in a failing ventricle. An axisymmetric fluid dynamics model of the left ventricle was developed and coupled to a lumped parameter model of the complete circulation. Simulations were performed for healthy conditions and dilated cardiomyopathy (DCM). Vortex structures in these simulations were analysed by means of automated detection. Results show that the strength of the leading vortex ring is lower in a DCM ventricle than in a healthy ventricle. The LVAD further influences the maximum strength of the vortex and also causes the vortex to disappear earlier in time with increasing LVAD flows. Understanding these phenomena by means of the method proposed in this study will contribute to enhanced diagnostics and monitoring during cardiac support.     

Transthoracic Doppler echocardiography – noninvasive diagnostic window for coronary flow reserve assessment

Background

Strain Rate Imaging shows the filling phases of the left ventricle to consist of a wave of myocardial stretching, propagating from base to apex. The propagation velocity of the strain rate wave is reduced in delayed relaxation. This study examined the relation between the propagation velocity of strain rate in the myocardium and the propagation velocity of flow during early filling.

Methods

12 normal subjects and 13 patients with treated hypertension and normal systolic function were studied. Patients and controls differed significantly in diastolic early mitral flow measurements, peak early diastolic tissue velocity and peak early diastolic strain rate, showing delayed relaxation in the patient group. There were no significant differences in EF or diastolic diameter.

Results

Strain rate propagation velocity was reduced in the patient group while flow propagation velocity was increased. There was a negative correlation (R = -0.57) between strain rate propagation and deceleration time of the mitral flow E-wave (R = -0.51) and between strain rate propagation and flow propagation velocity and there was a positive correlation (R = 0.67) between the ratio between peak mitral flow velocity/strain rate propagation velocity and flow propagation velocity.

Conclusion

The present study shows strain rate propagation to be a measure of filling time, but flow propagation to be a function of both flow velocity and strain rate propagation. Thus flow propagation is not a simple index of diastolic function in delayed relaxation.     

Targeted particle tracking in computational models of human carotid bifurcations

A significant and largely unsolved problem of computational fluid dynamics (CFD) simulation of flow in anatomically relevant geometries is that very few calculated pathlines pass through regions of complex flow. This in turn limits the ability of CFD-based simulations of imaging techniques (such as MRI) to correctly predict in vivo performance. In this work, I present two methods designed to overcome this filling problem, firstly, by releasing additional particles from areas of the flow inlet that lead directly to the complex flow region ("preferential seeding") and, secondly, by tracking particles both "downstream" and "upstream" from seed points within the complex flow region itself. I use the human carotid bifurcation as an example of complex blood flow that is of great clinical interest. Both idealized and healthy volunteer geometries are investigated. With uniform seeding in the inlet plane (in the common carotid artery (CCA)) of an idealized bifurcation geometry, approximately half the particles passed through the internal carotid artery (ICA) and half through the external carotid artery. However, of those particles entering the ICA, only 16% passed directly through the carotid bulb region. Preferential seeding from selected regions of the CCA was able to increase this figure to 47%. In the second method, seeding of particles within the carotid bulb region itself led to a very high proportion (97%) of pathlines running from CCA to ICA. Seeding of particles in the bulb plane of three healthy volunteer carotid bifurcation geometries led to much better filling of the bulb regions than by particles seeded at the inlet alone. In all cases, visualization of the origin and behavior of recirculating particles led to useful insights into the complex flow patterns. Both seeding methods produced significant improvements in filling the carotid bulb region with particle tracks compared with uniform seeding at the inlet and led to an improved understanding of the complex flow patterns. The methods described may be combined and are generally applicable to CFD studies of fluid and gas flow and are, therefore, of relevance in hemodynamics, respiratory mechanics, and medical imaging science.     

Sensitivity Analysis of Left Ventricle with Dilated Cardiomyopathy in Fluid Structure Simulation

Dilated cardiomyopathy (DCM) is the most common myocardial disease. It not only leads to systolic dysfunction but also diastolic deficiency. We sought to investigate the effect of idiopathic and ischemic DCM on the intraventricular fluid dynamics and myocardial wall mechanics using a 2D axisymmetrical fluid structure interaction model. In addition, we also studied the individual effect of parameters related to DCM, i.e. peak E-wave velocity, end systolic volume, wall compliance and sphericity index on several important fluid dynamics and myocardial wall mechanics variables during ventricular filling. Intraventricular fluid dynamics and myocardial wall deformation are significantly impaired under DCM conditions, being demonstrated by low vortex intensity, low flow propagation velocity, low intraventricular pressure difference (IVPD) and strain rates, and high-end diastolic pressure and wall stress. Our sensitivity analysis results showed that flow propagation velocity substantially decreases with an increase in wall stiffness, and is relatively independent of preload at low-peak E-wave velocity. Early IVPD is mainly affected by the rate of change of the early filling velocity and end systolic volume which changes the ventriculo:annular ratio. Regional strain rate, on the other hand, is significantly correlated with regional stiffness, and therefore forms a useful indicator for myocardial regional ischemia. The sensitivity analysis results enhance our understanding of the mechanisms leading to clinically observable changes in patients with DCM.     

Spatial velocity profile in mouse embryonic aorta and Doppler-derived volumetric flow: a preliminary model

Characterizing embryonic circulatory physiology requires accurate cardiac output and flow data. Despite recent applications of high-frequency ultrasound Doppler to the study of embryonic circulation, current Doppler analysis of volumetric flow is relatively crude. To improve Doppler derivation of volumetric flow, we sought a preliminary model of the spatial velocity profile in the mouse embryonic dorsal aorta using ultrasound biomicroscopy (UBM)-Doppler data. Embryonic hematocrit is 0.05-0.10 so rheologic properties must be insignificant. Low Reynolds numbers (<500) and Womersley parameters (<0.76) suggest laminar flow. UBM demonstrated a circular dorsal aortic cross section with no significant tapering. Low Dean numbers (<100) suggest the presence of minimal skewing of the spatial velocity profile. The inlet length allows for fully developed flow. There is no apparent aortic wall pulsatility. Extrapolation of prior studies to these vessel diameters (300-350 microm) and flow velocities (~50-200 mm/s) suggests parabolic spatial velocity profiles. Therefore, mouse embryonic dorsal aortic blood flow may correspond to Poiseuille flow in a straight rigid tube with parabolic spatial velocity profiles. As a first approximation, these results are an important step toward precise in utero ultrasound characterization of blood flow within the developing mammalian circulation.     

Computer simulation of intraventricular flow and pressure gradients during diastole

A two-dimensional axisymmetric computer model is developed for the simulation of the filling flow in the left ventricle (LV). The computed results show that vortices are formed during the acceleration phases of the filling waves. During the deceleration phases these are amplified and convected into the ventricle. The ratio of the maximal blood velocity at the mitral valve (peak E velocity) to the flow wave propagation velocity (WPV) of the filling wave is larger than 1. This hemodynamic behavior is also observed in experiments in vitro (Steen and Steen, 1994, Cardiovasc. Res., 28, pp. 1821-1827) and in measurements in vivo with color M-mode Doppler echocardiography (Stugaard et al., 1994, J. Am. Coll. Cardiol., 24, 663-670). Computed intraventricular pressure profiles are similar to observed profiles in a dog heart (Courtois et al., 1988, Circulation, 78, pp. 661-671). The long-term goal of the computer model is to study the predictive value of noninvasive parameters (e.g., velocities measured with Doppler echocardiography) on invasive parameters (e.g., pressures, stiffness of cardiac wall, time constant of relaxation). Here, we show that higher LV stiffness results in a smaller WPV for a given peak E velocity. This result may indicate an inverse relationship between WPV and LV stiffness, suggesting that WPV may be an important noninvasive index to assess LV diastolic stiffness, LV diastolic pressure and thus atrial pressure (preload).     

Time-resolved DPIV analysis of vortex dynamics in a left ventricular model through bileaflet mechanical and porcine heart valve prostheses

The performance of the heart after a mitral valve replacement operation greatly depends on the flow character downstream of the valve. The design and implanting orientation of valves may considerably affect the flow development. A study of the hemodynamics of two orientations, anatomical and anti-anatomical, of the St. Jude Medical (SJM) bileaflet valve are presented and compared with those of the SJM Biocor porcine valve, which served also to represent the natural valve. We document the velocity field in a flexible, transparent (LV) using time-resolved digital particle image velocimetry (TRDPIV). Vortex formation and vortex interaction are two important physical phenomena that dominate the filling and emptying of the ventricle. For the three configurations, the following effects were examined: mitral valve inlet jet asymmetry, survival of vortical structures upstream of the aortic valve, vortex-induced velocities and redirection of theflow in abidance of the Biot-Savart law, domain segmentation, resonant times of vortical structures, and regions of stagnantflow. The presence of three distinct flow patterns, for the three configurations, was identified by the location of vortical structures and level of coherence corresponding to a significant variation in the turbulence level distribution inside the LV. The adverse effect of these observations could potentially compromise the efficiency of the LV and result in flow patterns that deviate from those in the natural heart.