A fluid dynamics study in a 50 cc pulsatile ventricular assist device: influence of heart rate variability

Although left ventricular assist devices (LVADs) have had success in supporting severe heart failure patients, thrombus formation within these devices still limits their long term use. Research has shown that thrombosis in the Penn State pulsatile LVAD, on a polyurethane blood sac, is largely a function of the underlying fluid mechanics and may be correlated to wall shear rates below 500 s(-1). Given the large range of heart rate and systolic durations employed, in vivo it is useful to study the fluid mechanics of pulsatile LVADs under these conditions. Particle image velocimetry (PIV) was used to capture planar flow in the pump body of a Penn State 50 cubic centimeters (cc) LVAD for heart rates of 75-150 bpm and respective systolic durations of 38-50%. Shear rates were calculated along the lower device wall with attention given to the uncertainty of the shear rate measurement as a function of pixel magnification. Spatial and temporal shear rate changes associated with data collection frequency were also investigated. The accuracy of the shear rate calculation improved by approximately 40% as the resolution increased from 35 to 12 μm/pixel. In addition, data collection in 10 ms, rather than 50 ms, intervals was found to be preferable. Increasing heart rate and systolic duration showed little change in wall shear rate patterns, with wall shear rate magnitude scaling by approximately the kinematic viscosity divided by the square of the average inlet velocity, which is essentially half the friction coefficient. Changes in in vivo operating conditions strongly influence wall shear rates within our device, and likely play a significant role in thrombus deposition. Refinement of PIV techniques at higher magnifications can be useful in moving towards better prediction of thrombosis in LVADs.     

Computational simulations of flow dynamics and blood damage through a bileaflet mechanical heart valve scaled to pediatric size and flow

Despite pressing needs, there are currently no FDA approved prosthetic valves available for use in the pediatric population. This study is performed for predictive assessment of blood damage in bileaflet mechanical heart valves (BMHVs) with pediatric sizing and flow conditions. A model of an adult-sized 23 mm St. Jude Medical (SJM) Regent^™ valve is selected for use in simulations, which is scaled in size for a 5-year old child and 6-month old infant. A previously validated lattice-Boltzmann method (LBM) is used to simulate pulsatile flow with thousands of suspended platelets for cases of adult, child, and infant BMHV flows. Adult BMHV flows demonstrate more disorganized small-scale flow features, but pediatric flows are associated with higher fluid shear stresses. Platelet damage in the pediatric cases is higher than in adult flow, highlighting thrombus complication dangers of pediatric BMHV flows. This does not necessarily suggest clinically important differences in thromboembolic potential. Highly damaged platelets in pediatric flows are primarily found far downstream of the valve, as there is less flow recirculation in pediatric flows. In addition, damage levels are well below expected thresholds for platelet activation. The extent of differences here documented between the pediatric and adult cases is of concern, demanding particular attention when pediatric valves are designed and manufactured. However, the differences between the pediatric and adult cases are not such that development of pediatric sized valves is untenable. This study may push for eventual approval of prosthetic valves resized for the pediatric population. Further studies will be necessary to determine the validity and potential thrombotic and clinical implications of these findings.     

Turbulence characteristics downstream of bileaflet aortic valve prostheses

This study was focused on a series of in vitro tests on the turbulent flow characteristics of three bileaflet aortic valves: St. Jude Medical (SJM), CarboMedics (CM), and Edwards Tekna (modified Duromedics, DM). The flow fields of the valves were measured in a pulsatile flow model with a laser-Doppler anemometer (LDA) at the aortic sinus area downstream of the valves. The heart rate was set at 70 beats per minute, the cardiac output was maintained at 5 liters per minute, and the aortic pressure wave forms were kept within the physiological range. Cycle-resolved analysis was applied to obtain turbulence data, including mean velocity, Reynolds stresses, autocorrelation coefficients, energy spectral density functions, and turbulence scales. The Reynolds shear stresses of all three valves induced only minor damage to red blood cells, but directly damaged the platelets, increasing the possibility of thrombosis. The smallest turbulence length scale, which offers a more reliable estimate of the effects of turbulence on blood cell damage, was three times the size of red blood cells and five times the size of platelets. This suggests that there is more direct interaction with the blood cells, thus causing more damage.     

In vitro flow dynamics of four prosthetic aortic valves: a comparative analysis

The velocity fields downstream of four prosthetic heart valves were mapped in vitro over the entire cross-section of a model aortic root using laser Doppler anemometry. THe Bj?rk-Shiley 60 degrees convexo-concave tilting disc valve, the Smeloff-Cutter caged ball valve, the St. Jude Medical bileaflet valve, and the Ionescu-Shiley standard bioprosthesis were examined under both steady and pulsatile flows. Velocity profiles under steady flow conditions were a good approximation for pulsatile profiles only during midsystole. The pulsatile flow characteristics of the four valves showed variation in large scale flow structures. Comparison of the valves according to pressure drop, shear stress and maximum velocities are also provided.     

Flush mounted hot film anemometer measurement of wall shear stress distal to a tri-leaflet valve for Newtonian and non-Newtonian blood analog fluids

Wall shear stress has been measured by flush-mounted hot film anemometry distal to an Ionescu-Shiley tri-leaflet valve under pulsatile flow conditions. Both Newtonian (aqueous glycerol) and non-Newtonian (aqueous polyacrylamide) blood analog fluids were investigated. Significant differences in the axial distribution of wall shear stress between the two fluids are apparent in flows having nearly identical Reynolds numbers. The Newtonian fluid exhibits a (peak) wall shear rate which is maximized near the valve seat (30 mm) and then decays to a fully developed flow value (by 106 mm). In contrast, the shear rate of the non-Newtonian fluid at 30 mm is less than half that of the Newtonian fluid and at 106 mm is more than twice that of the Newtonian fluid. It is suggested that non-Newtonian rheology influences valve flow patterns either through alterations in valve opening associated with low shear separation zones behind valve leaflets, or because of variations in the rate of jet spreading. More detailed studies are required to clarify the mechanisms. The Newtonian wall shear stresses for this valve are low. The highest value observed anywhere in the aortic chamber was 2.85 N/m2 at a peak Reynolds number of 3694.     

Superior performance of continuous over pulsatile flow ventricular assist devices in the single ventricle circulation: A computational study

This study compares the physiological responses of systemic-to-pulmonary shunted single ventricle patients to pulsatile and continuous flow ventricular assist devices (VADs). Performance differences between pulsatile and continuous flow VADs have been clinically observed, but the underlying mechanism remains poorly understood. Six systemic-to-pulmonary shunted single ventricle patients (mean BSA=0.30 m²) were computationally simulated using a lumped-parameter network tuned to match patient specific clinical data. A first set of simulations compared current clinical implementation of VADs in single ventricle patients. A second set modified pulsatile flow VAD settings with the goal to optimize cardiac output (CO). For all patients, the best-case continuous flow VAD CO was at least 0.99 L/min greater than the optimized pulsatile flow VAD CO (p=0.001). The 25 and 50 mL pulsatile flow VADs exhibited incomplete filling at higher heart rates that reduced CO as much as 9.7% and 37.3% below expectations respectively. Optimization of pulsatile flow VAD settings did not achieve statistically significant (p<0.05) improvement to CO. Results corroborate clinical experience that continuous flow VADs produce higher CO and superior ventricular unloading in single ventricle patients. Impaired filling leads to performance degradation of pulsatile flow VADs in the single ventricle circulation.     

Two-component laser velocimeter measurements downstream of heart valve prostheses in pulsatile flow

Elevated turbulent shear stresses resulting from disturbed blood flow through prosthetic heart valves can cause damage to red blood cells and platelets. The purpose of this study was to measure the turbulent shear stresses occurring downstream of aortic prosthetic valves during in-vitro pulsatile flow. By matching the indices of refraction of the blood analog fluid and model aorta, correlated, simultaneous two-component laser velocimeter measurements of the axial and radial velocity components were made immediately downstream of two aortic prosthetic valves. Velocity data were ensemble averaged over 200 or more cycles for a 15-ms window opened at peak systolic flow. The systolic duration for cardiac flows of 8.4 L/min was 200 ms. Ensemble-averaged total shear stress levels of 2820 dynes/cm2 and 2070 dynes/cm2 were found downstream of a trileaflet valve and a tilting disk valve, respectively. These shear stress levels decreased with axial distance downstream much faster for the tilting disk valve than for the trileaflet valve.     

Blood flow dynamics in saccular aneurysm models of the basilar artery

Blood flow dynamics under physiologically realistic pulsatile conditions plays an important role in the growth, rupture, and surgical treatment of intracranial aneurysms. The temporal and spatial variations of wall pressure and wall shear stress in the aneurysm are hypothesized to be correlated with its continuous expansion and eventual rupture. In addition, the assessment of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils. This paper describes the flow dynamics in two representative models of a terminal aneurysm of the basilar artery under Newtonian and non-Newtonian fluid assumptions, and compares their hemodynamics with that of a healthy basilar artery. Virtual aneurysm models are investigated numerically, with geometric features defined by beta = 0 deg and beta = 23.2 deg, where beta is the tilt angle of the aneurysm dome with respect to the basilar artery. The intra-aneurysmal pulsatile flow shows complex ring vortex structures for beta = 0 deg and single recirculation regions for beta = 23.2 deg during both systole and diastole. The pressure and shear stress on the aneurysm wall exhibit large temporal and spatial variations for both models. When compared to a non-Newtonian fluid, the symmetric aneurysm model (beta = 0 deg) exhibits a more unstable Newtonian flow dynamics, although with a lower peak wall shear stress than the asymmetric model (beta = 23.2 deg). The non-Newtonian fluid assumption yields more stable flows than a Newtonian fluid, for the same inlet flow rate. Both fluid modeling assumptions, however, lead to asymmetric oscillatory flows inside the aneurysm dome.     

Flush-mounted hot film anemometer accuracy in pulsatile flow

The accuracy of a flush-mounted hot film anemometer probe for wall shear stress measurements in physiological pulsatile flows was evaluated in fully developed pulsatile flow in a rigid straight tube. Measured wall shear stress waveform based on steady flow anemometer probe calibrations were compared to theoretical wall shear stress waveforms based on well-established theory and measured flow rate waveforms. The measured and theoretical waveforms were in close agreement during systole (average deviation of 14 percent at peak systole). As expected, agreement was poor during diastole because of flow reversal and diminished frequency response at low shear rate.     

Influence of cardiac flow rate on turbulent shear stress from a prosthetic heart valve

Elevated turbulent shear stresses associated with sufficient exposure times are potentially damaging to blood constituents. Since these conditions can be induced by mechanical heart valves, the objectives of this study were to locate the maximum turbulent shear stress in both space and time and to determine how the maximum turbulent shear stress depends on the cardiac flow rate in a pulsatile flow downstream of a tilting disk valve. Two-component, simultaneous, correlated laser velocimeter measurements were recorded at four different axial locations and three different flow rates in a straight tube model of the aorta. All velocity data were ensemble averaged within a 15 ms time window located at approximately peak systolic flow over more than 300 cycles. Shear stresses as high as 992 dynes/cm2 were found 0.92 tube diameters downstream of the monostrut, disk valve. The maximum turbulent shear stress was found to scale with flow rate to the 0.72 power. A repeatable starting vortex was shed from the disk at the beginning of each cycle.     

Pulsatile flow past aortic valve bioprostheses in a model human aorta

Pulsatile flow development past tissue valve prostheses in a model human aorta has been studied using qualitative flow visualization and quantitative laser-Doppler techniques. Experiments were conducted both in steady and physiological pulsatile flow situations and the measurements included the pressure drop across the valve, the instantaneous flow rate as well as the velocity profiles and turbulent stresses downstream to the valves. Our study shows that the velocity profiles with pericardial valves are closer to those measured past natural aortic valves. The porcine valves with a smaller valve opening area produce a narrower and stronger jet downstream from the valve with relatively larger turbulent axial stresses in the boundary of the jet. Our study suggests that the pericardial valves with turbulent stresses comparable to those of caged ball and tilting disc valves are preferable from a hemodynamic point of view.     

Laser anemometry measurements of pulsatile flow past aortic valve prostheses

Experimental results are presented on physiological pulsatile flow past caged ball and tilting disc aortic valve prostheses mounted in an axisymmetric chamber incorporated in a mock circulatory system. The measurements of velocity profiles and turbulent normal stresses during several times in a cardiac cycle were obtained using laser-Doppler anemometry. Our results show that with increased angle of opening for the tilting disc valves, a large but locally confined vortex is observed along the wall in the minor flow region throughout most of the cardiac cycle. The turbulent normal stresses measured downstream to the tilting disc in the minor flow region parallel to the tilt axis were found to be larger than those measured downstream to the caged ball valves. Comparison of measurements with steady flow at flow rates comparable to peak pulsatile flow rate show that the turbulent normal stresses are larger by a factor of two in pulsatile flow with a frequency of 1.2 Hz.     

The congenital bicuspid aortic valve can experience high-frequency unsteady shear stresses on its leaflet surface

The bicuspid aortic valve (BAV) is a common congenital malformation of the aortic valve (AV) affecting 1% to 2% of the population. The BAV is predisposed to early degenerative calcification of valve leaflets, and BAV patients constitute 50% of AV stenosis patients. Although evidence shows that genetic defects can play a role in calcification of the BAV leaflets, we hypothesize that drastic changes in the mechanical environment of the BAV elicit pathological responses from the valve and might be concurrently responsible for early calcification. An in vitro model of the BAV was constructed by surgically manipulating a native trileaflet porcine AV. The BAV valve model and a trileaflet AV (TAV) model were tested in an in vitro pulsatile flow loop mimicking physiological hemodynamics. Laser Doppler velocimetry was used to make measurements of fluid shear stresses on the leaflet of the valve models using previously established methodologies. Furthermore, particle image velocimetry was used to visualize the flow fields downstream of the valves and in the sinuses. In the BAV model, flow near the leaflets and fluid shear stresses on the leaflets were much more unsteady than for the TAV model, most likely due to the moderate stenosis in the BAV and the skewed forward flow jet that collided with the aorta wall. This additional unsteadiness occurred during mid- to late-systole and was composed of cycle-to-cycle magnitude variability as well as high-frequency fluctuations about the mean shear stress. It has been demonstrated that the BAV geometry can lead to unsteady shear stresses under physiological flow and pressure conditions. Such altered shear stresses could play a role in accelerated calcification in BAVs.     

Pressure drop and flow rate measurements in a human aortic bifurcation cast for steady and pulsatile flow

Pressure drop and flow rate measurements in a rigid cast of a human aortic bifurcation under both steady and physiological pulsatile flow conditions are reported. Integral momentum and mechanical energy balances are used to calculate impedance, spatially averaged wall shear stress and viscous dissipation rate from the data. In the daughter branches, steady flow impedance is within 30% of the Poiseuille flow prediction, while pulsatile flow impedance is within a factor of 2 of fully developed, oscillatory, straight tube flow theory (Womersley theory). Estimates of wall shear stress are in accord with measurements obtained from velocity profiles. Mean pressure drop and viscous dissipation rate are elevated in pulsatile flow relative to steady flow at the mean flow rate, and the exponents of their Reynolds number dependence are in accord with available theory.     

Pulsatile extracorporeal circulation during on-pump cardiac surgery enhances aortic wall shear stress

Controversy on superiority of pulsatile versus non-pulsatile extracorporeal circulation in cardiac surgery still continues. Stroke as one of the major adverse events during cardiopulmonary bypass is, in the majority of cases, caused by mobilization of aortic arteriosclerotic plaques that is inducible by pathologically elevated wall shear stress values. The present study employs computational fluid dynamics to evaluate the aortic blood flow and wall shear stress profiles under the influence of antegrade or retrograde perfusion with pulsatile versus non-pulsatile extracorporeal circulation. While, compared to physiological flow, a non-pulsatile perfusion resulted in generally decreased blood velocities and only moderately increased shear forces (48 Pa versus 20 Pa antegradely and 127 Pa versus 30 Pa retrogradely), a pulsatile perfusion extensively enhanced the occurrence of turbulences, maximum blood flow speed and maximum wall shear stress (1020 Pa versus 20 Pa antegradely and 1178 Pa versus 30 Pa retrogradely). Under these circumstances arteriosclerotic embolism has to be considered. Further simulations and experimental work are necessary to elucidate the impact of our findings on the scientific discourse of pulsatile versus non-pulsatile extracorporeal circulation.     

A numerical and experimental investigation of transitional pulsatile flow in a stenosed channel

In the present paper, a closely coupled numerical and experimental investigation of pulsatile flow in a prototypical stenotic site is presented. Detailed laser Doppler velocimetry measurements upstream of the stenosis are used to guide the specification of velocity boundary conditions at the inflow plane in a series of direct numerical simulations (DNSs). Comparisons of the velocity statistics between the experiments and DNS in the post-stenotic area demonstrate the great importance of accurate inflow conditions, and the sensitivity of the post-stenotic flow to the disturbance environment upstream. In general, the results highlight a borderline turbulent flow that sequentially undergoes transition to turbulence and relaminarization. Before the peak mass flow rate, the strong confined jet that forms just downstream of the stenosis becomes unstable, forcing a role-up and subsequent breakdown of the shear layer. In addition, the large-scale structures originating from the shear layer are observed to perturb the near wall flow, creating packets of near wall hairpin vortices.     

Flow-induced wall shear stress in abdominal aortic aneurysms: Part II--pulsatile flow hemodynamics

In continuing the investigation of AAA hemodynamics, unsteady flow-induced stresses are presented for pulsatile blood flow through the double-aneurysm model described in Part I. Physiologically realistic aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50< or =Re(m) < or =300. Hemodynamic disturbance is evaluated for a modified set of indicator functions which include wall pressure (p(w)), wall shear stress (tau(w)), Wall Shear Stress Gradient (WSSG), time-average wall shear stress (tau(w)*), and time-average Wall Shear Stress Gradient WSSG*. At peak flow, the highest shear stress and WSSG levels are obtained at the distal end of both aneurysms, in a pattern similar to that of steady flow. The maximum values of wall shear stresses and wall shear stress gradients are evaluated as a function of the time-average Reynolds number resulting in a fourth order polynomial correlation. A comparison between numerical predictions for steady and pulsatile flow is presented, illustrating the importance of considering time-dependent flow for the evaluation of hemodynamic indicators.     

Flow-induced Wall Shear Stress in Abdominal Aortic Aneurysms: Part II - Pulsatile Flow Hemodynamics

In continuing the investigation of AAA hemodynamics, unsteady flow-induced stresses are presented for pulsatile blood flow through the double-aneurysm model described in Part I. Physiologically realistic aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 h Re m h 300. Hemodynamic disturbance is evaluated for a modified set of indicator functions which include wall pressure ( p w ), wall shear stress ( w ), Wall Shear Stress Gradient (WSSG), time-average wall shear stress ( w *), and time-average Wall Shear Stress Gradient WSSG *. At peak flow, the highest shear stress and WSSG levels are obtained at the distal end of both aneurysms, in a pattern similar to that of steady flow. The maximum values of wall shear stresses and wall shear stress gradients are evaluated as a function of the time-average Reynolds number resulting in a fourth order polynomial correlation. A comparison between numerical predictions for steady and pulsatile flow is presented, illustrating the importance of considering time-dependent flow for the evaluation of hemodynamic indicators.     

In vitro pulsatile flow velocity and shear stress measurements in the vicinity of mechanical mitral heart valve prostheses

A three beam laser Doppler anemometer system was used to study the flow fields created by various types of mitral heart valve prostheses under physiological pulsatile flow conditions. The prosthetic valves studied were: Beall caged disc valve, Bjork-Shiley tilting disc valve, Medtronic-Hall tilting disc valve and St. Jude bileaflet valve. The results indicate that all four prosthetic valve designs studied create very disturbed flow fields with elevated turbulent shear stresses and regions of flow separation and/or stagnation. The observed elevated turbulent shear stresses could cause sublethal and/or lethal damage to red cells and platelets. The regions of flow separation and/or stagnation, could lead to thrombus formation and/or tissue overgrowth on the valve structure, as observed on clinically recovered prosthetic valves.     

Two-dimensional color-mapping of turbulent shear stress distribution downstream of two aortic bioprosthetic valves in vitro.

Since artificial heart valve related complications such as thrombus formation, hemolysis and calcification are considered related to flow disturbances caused by the inserted valve, a thorough hemodynamic characterization of heart valve prostheses is essential. In a pulsatile flow model, fluid velocities were measured one diameter downstream of a Hancock Porcine (HAPO) and a Ionescu-Shiley Pericardial Standard (ISPS) aortic valve. Hot-film anemometry (HFA) was used for velocity measurements at 41 points in the cross-sectional area of the ascending aorta. Three-dimensional visualization of the velocity profiles, at 100 different instants during one mean pump cycle, was performed. Turbulence analysis was performed as a function of time by calculating the axial turbulence energy within 50 ms overlapping time windows during the systole. The turbulent shear stresses were estimated by using the correlation equation between Reynolds normal stress and turbulent (Reynolds) shear stress. The turbulent shear stress distribution was visualized by two-dimensional color-mapping at different instants during one mean pump cycle. Based on the velocity profiles and the turbulent shear stress distribution, a relative blood damage index (RBDI) was calculated. It has the feature of combining the magnitude and exposure time of the estimated shear stresses in one index, covering the entire cross-sectional area. The HAPO valve showed a skewed jet-type velocity profile with the highest velocities towards the left posterior aortic wall. The ISPS valve revealed a more parabolic-shaped velocity profile during systole. The turbulent shear stresses were highest in areas of high or rapidly changing velocity gradients. For the HAPO valve the maximum estimated turbulent shear stress was 194 N m-2 and for the ISPS valve 154 Nm-2. The RBDI was the same for the two valves. The turbulent shear stresses had magnitudes and exposure times that might cause endothelial damage and sublethal or lethal damage to blood corpuscules. The RBDI makes comparison between different heart valves easier and may prove important when making correlation with clinical observations.