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.     

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.     

PIV measurements of flow in a centrifugal blood pump: steady flow

Magnetically suspended left ventricular assist devices have only one moving part, the impeller. The impeller has absolutely no contact with any of the fixed parts, thus greatly reducing the regions of stagnant or high shear stress that surround a mechanical or fluid bearing. Measurements of the mean flow patterns as well as viscous and turbulent stresses were made in a shaft-driven prototype of a magnetically suspended centrifugal blood pump at several constant flow rates (3-9 L/min) using particle image velocimetry (PIV). The chosen range of flow rates is representative of the range over which the pump may operate while implanted. Measurements on a three-dimensional measurement grid within several regions of the pump, including the inlet, blade passage, exit volute, and diffuser are reported. The measurements are used to identify regions of potential blood damage due to high shear stress and/or stagnation of the blood, both of which have been associated with blood damage within artificial heart valves and diaphragm-type pumps. Levels of turbulence intensity and Reynolds stresses that are comparable to those in artificial heart valves are reported. At the design flow rate (6 L/min), the flow is generally well behaved (no recirculation or stagnant flow) and stress levels are below levels that would be expected to contribute to hemolysis or thrombosis. The flow at both high (9 L/min) and low (3 L/min) flow rates introduces anomalies into the flow, such as recirculation, stagnation, and high stress regions. Levels of viscous and Reynolds shear stresses everywhere within the pump are below reported threshold values for damage to red cells over the entire range of flow rates investigated; however, at both high and low flow rate conditions, the flow field may promote activation of the clotting cascade due to regions of elevated shear stress adjacent to separated or stagnant flow.     

A reevaluation and discussion on the threshold limit for hemolysis in a turbulent shear flow

One major factor of red blood cell damage induced by artificial heart valves is the magnitude of turbulent shear stresses in the flow field. An often-cited threshold for hemolysis is 400 N/m(2) (Sallam and Hwang, Biorheology 21 (1984) 783). This value, however, was measured with a one-component laser Doppler anemometer and was not calculated from the major principal Reynolds shear stress, therefore underestimating the threshold. This study performed flow field measurements under similar conditions, using a two-component laser Doppler anemometer and determining major principal Reynolds shear stress, to reevaluate the hemolytic threshold as 800 N/m(2) with an exposure time of 1 ms. The Kolmogorov length scales, approximately 9 microm, are capable of inflicting direct damage to the red blood cells. These results will serve as a more accurate reference in future heart valve design and testing.     

Turbulent stresses in the region of a Hancock porcine bioprosthetic aortic valve

The purpose of this study was to measure stresses associated with turbulence (Reynolds stresses), in the region of a 29-mm-dia porcine bioprosthetic valve (Hancock, Model 242). Studies were performed in an in vitro pulse duplicating system with the valve mounted in the aortic position. The Reynolds stresses were calculated from velocities obtained with a two channel laser Doppler anemometer. The largest Reynolds shear stress and normal stress occurred at the highest stroke volume used (80 mL). Averaged over ejection they were 38 dynes/cm2 and 380 dynes/cm2, respectively. The maximal instantaneous Reynolds shear stress was 2500 dynes/cm2 and the maximal instantaneous Reynolds normal stress was 6800 dynes/cm2. Stresses of these magnitudes are in the range reported to damage platelets.     

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.     

Turbulence downstream of subcoronary stentless and stented aortic valves

Regions of turbulence downstream of bioprosthetic heart valves may cause damage to blood components, vessel wall as well as to aortic valve leaflets. Stentless aortic heart valves are known to posses several hemodynamic benefits such as larger effective orifice areas, lower aortic transvalvular pressure difference and faster left ventricular mass regression compared with their stented counterpart. Whether this is reflected by diminished turbulence formation, remains to be shown. We implanted either stented pericardial valve prostheses (Mitroflow), stentless valve prostheses (Solo or Toronto SPV) in pigs or they preserved their native valves. Following surgery, blood velocity was measured in the cross sectional area downstream of the valves using 10MHz ultrasonic probes connected to a dedicated pulsed Doppler equipment. As a measure of turbulence, Reynolds normal stress (RNS) was calculated at two different blood pressures (baseline and 50% increase). We found no difference in maximum RNS measurements between any of the investigated valve groups. The native valve had significantly lower mean RNS values than the Mitroflow (p=0.004), Toronto SPV (p=0.008) and Solo valve (p=0.02). There were no statistically significant differences between the artificial valve groups (p=0.3). The mean RNS was significantly larger when increasing blood pressure (p=0.0006). We, thus, found no advantages for the stentless aortic valves compared with stented prosthesis in terms of lower maximum or mean RNS values. Native valves have a significantly lower mean RNS value than all investigated bioprostheses.     

On the monodimensional approach to the estimation of the highest reynolds shear stress in a turbulent flow

The measurement of the Reynolds stress tensor, or at least of some of its components, is a necessary step to assess if the turbulence associated with the flow near prosthetic devices can damage blood constituents. Because of the intrinsic three dimensionality of turbulence, in general, a three-component anemometer should be used to measure directly the components of the Reynolds stress tensor. However, this can be practically unfeasible, especially in vivo; therefore, it is interesting to investigate the possibility of characterizing the turbulent flows that may occur in the circulatory system with the monodimensional data that a less complete equipment (e.g., a pulsed ultrasound Doppler) can yield. From the general expression of the Reynolds stress tensor, the highest shear stress can be deduced, as well as the Reynolds normal stress in the main flow direction. The relation between these two quantities, which is an issue already addressed in previous works, can thus be rigorously formulated in terms of some characteristic parameters of the Reynolds stress tensor, the principal normal stresses and the angles that the directions that define them form with the main flow direction. An experimental verification of the ratio of the two above-mentioned quantitites for the flow across bileaflet valves, investigated by means of two-dimensional laser Doppler anemometry, will illustrate the limitations of the monodimensional approach estimating the maximum load on blood constituents.     

Turbulent shear stress measurements in the vicinity of aortic heart valve prostheses

A two dimensional laser Doppler anemometer system has been used to measure the turbulent shear fields in the immediate downstream vicinity of a variety of mechanical and bioprosthetic aortic heart valves. The measurements revealed that all the mechanical valves studied, created regions of elevated levels of turbulent shear stress during the major portion of systole. The tissue bioprostheses also created elevated levels of turbulence, but they were confined to narrow regions in the bulk of the flow field. The newer generation of bioprostheses create turbulent shear stresses which are considerably lower than those created by the older generation tissue valve designs. All the aortic valves studied (mechanical and tissue) create turbulent shear stress levels which are capable of causing sub-lethal and/or lethal damage to blood elements.     

Estimation of turbulent shear stresses in pulsatile flow immediately downstream of two artificial aortic valves in vitro

Measuring turbulent shear stresses is of major importance in artificial heart valve evaluation. Bi- and unidirectional fluid velocity measurements enable calculation of Reynolds shear stress ( ) and Reynolds normal stress ( ). τ is important due to the relation to hemolysis and thrombus formation, but σ is the only obtainable parameter in vivo. Therefore, determination of a correlation factor between τ and σ is pertinent.

In a pulsatile flow model, laser Doppler (LDA) and hot-film (HFA) anemometry were used for simultaneous bi- and unidirectional fluid velocity measurements downstream of a Hall Kaster and a Hancock Porcine aortic valve. Velocities were registered in two flow field locations and at four cardiac outputs. The velocity signals were subjected to analog signal processing prior to digital turbulence analysis, as a basis for calculation of τ and σ.

A correlation factor of 0.5 with a correlation coefficient of 0.97 was found between the maximum Reynolds shear stress and Reynolds normal stress, implying . In vitro estimation of turbulent shear stresses downstream of artificial aortic valves, based on the axial velocity component alone, seems possible.     

Improved flow cytometric method to enumerate residual cells: minimal linear detection limits for platelets,erythrocytes, and leukocytes

Increasing demand for quality control of blood products requires more sensitive methods to enumerate residual cells. Presently, the reported threshold (in cells per microliter) is 400 for red blood cells, 30-500 for platelets, and 1 for leukocytes. To examine precision and linearity in enumerating residual platelets and red blood cells, EDTA-anticoagulated blood from healthy donors was serially diluted with serum, stained in TruCount tubes using a no-lyse/no-wash procedure and a monoclonal antibody cocktail against the CD42a (FL1) and glycophorin-A (FL2) epitopes, and analyzed by flow cytometry. Leukocyte counts were determined in separate tubes. Cell preparation and analysis were performed once for 20 blood samples each and 20 times using the same specimen. Acquisition from the same tube was performed separately for platelets (threshold on FL1) and red blood cells (threshold on FL2). Multiparameter analysis was used for data evaluation. Linear results were obtained for platelets per microliter between 3,410 and 5 and for red blood cells per microliter between 54,000 and 3. For the lower cell concentrations, the coefficient of variation was 16.7% for platelets and 10.9% for red blood cells. The presented method allows the distinction between physiologically intact and ghost red blood cells. The method represents a reliable, sensitive, and accurate approach to quantify platelets and red blood cells in diluted blood. It can be applied to enumerate residual cells in plasma products and meets the increasing demand for quality control in blood components.     

Turbulent stress measurements downstream of six mechanical aortic valves in a pulsatile flow model

In a pulsatile flow model aortic Bj?rk-Shiley Standard, Convex-Concave and Monostrut valves were investigated together with the Hall-Kaster (Medtronic-Hall), St Jude Medical and Starr-Edwards Silastic Ball valve using hot-film anemometry. Three-dimensional visualization of average systolic Reynolds normal stresses (RNS) reflected the design of the valves. Mean average RNS were used for comparison of the fluid dynamic performance along with Velocity Energy Ratio (VER100) and Turbulence Energy Ratio (TER) as a relative turbulence intensity for pulsatile flow. Mean average RNS ranged from 13.2 to 37.6 Nm-2 for all the valves with the highest levels for the Bj?rk-Shiley Standard and Starr-Edwards Ball valve and lowest values for the St Jude Medical valve and with the Hall-Kaster (Medtronic-Hall), Bj?rk-Shiley Convex-Concave and Monostrut valves in between.     

A comparative computational study of blood flow through prosthetic heart valves using the finite element method

The steady flow of blood through three common types of prosthetic heart valves was simulated numerically using the finite element method. The velocity, pressure and stress fields were obtained for the disk-type, tilting-disk and ball-type prosthetic heart valves in aortic position, for increasing Reynolds numbers up to 900, 1500 and 2000 respectively. Computer graphics of calculated velocities are presented, showing in detail the accelerated flow, recirculation and stagnation areas developed in the prosthesis. Maximum wall shear stresses were found at 0.5, 1.4, and 1.2 diameters from the sewing ring downstream for the disk, tilting-disk and ball valves being the values of 55, 18 and 33 dyn cm-2 respectively. In the vicinity of the occluder, maximum shear stresses of 38, 30 and 47 dyn cm-2 respectively were computed. The flow characteristics and performance for each valve are compared, the results are presented in terms of energy loss and maximum shear stress. The velocity and stress fields are compared with in vitro evaluations found in the literature.     

Estimation of turbulent shear stresses in pulsatile flow immediately downstream of two artificial aortic valves in vitro

Measuring turbulent shear stresses is of major importance in artificial heart valve evaluation. Bi- and unidirectional fluid velocity measurements enable calculation of Reynolds shear stress ( ) and Reynolds normal stress ( ). τ is important due to the relation to hemolysis and thrombus formation, but σ is the only obtainable parameter in vivo. Therefore, determination of a correlation factor between τ and σ is pertinent.In a pulsatile flow model, laser Doppler (LDA) and hot-film (HFA) anemometry were used for simultaneous bi- and unidirectional fluid velocity measurements downstream of a Hall Kaster and a Hancock Porcine aortic valve. Velocities were registered in two flow field locations and at four cardiac outputs. The velocity signals were subjected to analog signal processing prior to digital turbulence analysis, as a basis for calculation of τ and σ.A correlation factor of 0.5 with a correlation coefficient of 0.97 was found between the maximum Reynolds shear stress and Reynolds normal stress, implying . In vitro estimation of turbulent shear stresses downstream of artificial aortic valves, based on the axial velocity component alone, seems possible.     

Augmentation of sickling process due to turbulent blood flow.

The purpose of this study was to determine if the fluid mechanical stresses associated with turbulent blood flow can contribute to the sickling process. Blood from seven patients with sickle cell disease was subjected to intermediate and high levels of turbulent flow in vitro. Turbulence was quantitated by hot film anemometry. Control samples showed 20 +/- 3% sickled cells. Cells subjected to intermediate levels of turbulent flow showed 26 +/- 4% sickling (P less than 0.01); and blood subjected to high intensities of turbulence showed 31 +/- 4% sickling (P less than 0.01). A quantitative count by electronmicroscopy, performed in one patient, showed polymerization of the hemoglobin indicative of sickling in more cells subjected to turbulence than in the control sample. A turbulence-reducing agent, polyethylene oxide, diminished the augmentation of the sickling process as it reduced turbulence at comparable Reynolds numbers. These results support the hypothesis that a deleterious effect upon hemoglobin SS erythrocytes may occur due to the mechanical stresses of turbulent flow. The agitation associated with turbulent flow presumably modifies the stabilizing factors of the intracellular colloidal solution of hemoglobin, thereby contributing to sol-gel transformation. Such hydrodynamic stresses may supplement the previously described factors which contribute to sickle cell crises.     

Spectrum analysis of turbulence in the canine ascending aorta measured with a hot-film anemometer

We measured turbulence velocity in the canine ascending aorta using a hot-film anemometer. Blood flow velocity was measured at various points across the ascending aorta approximately 1.5-2 times the diameter downstream from the aortic valve. The turbulence spectrum was calculated and its characteristics were examined in connection with the mean Reynolds number and/or measuring positions. In the higher wave number range the values of the turbulence spectra were higher at larger mean Reynolds number. In the higher wave number range, the values of the turbulence spectra were higher at points closer to the centerline of the aorta, when the mean Reynolds number was relatively large. The patterns of the turbulence spectra at various points outside the boundary layer on the aortic wall were similar.     

Numerical and experimental investigations of pulsatile blood flow pattern through a dysfunctional mechanical heart valve

Around 250,000 heart valve replacements are performed every year around the world. Due their higher durability, approximately 2/3 of these replacements use mechanical prosthetic heart valves (mainly bileaflet valves). Although very efficient, these valves can be subject to valve leaflet malfunctions. These malfunctions are usually the consequence of pannus ingrowth and/or thrombus formation and represent serious and potentially fatal complications. Hence, it is important to investigate the flow field downstream of a dysfunctional mechanical heart valve to better understand its impact on blood components (red blood cells, platelets and coagulation factors) and to improve the current diagnosis techniques. Therefore, the objective of this study will be to numerically and experimentally investigate the pulsatile turbulent flow downstream of a dysfunctional bileaflet mechanical heart valve in terms of velocity field, vortex formation and potential negative effect on blood components. The results show that the flow downstream of a dysfunctional valve was characterized by abnormally elevated velocities and shear stresses as well as large scale vortices. These characteristics can predispose to blood components damage. Furthermore, valve malfunction led to an underestimation of maximal transvalvular pressure gradient, using Doppler echocardiography, when compared to numerical results. This could be explained by the shifting of the maximal velocity towards the normally functioning leaflet. As a consequence, clinicians should try, when possible, to check the maximal velocity position not only at the central orifice but also through the lateral orifices. Finding the maximal velocity in the lateral orifice could be an indication of valve dysfunction.     

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.     

Velocity measurements in steady flow through axisymmetric stenoses at moderate Reynolds numbers

The velocity field in the neighborhood of axisymmetric constrictions in rigid tubes was investigated using laser Doppler anemometry and flow visualization. Upstream flow conditions were steady; and Reynolds numbers were in the range 500-2000, values which are representative of the larger arteries in humans. Stenoses of 25, 50 and 75% area reduction were studied. Velocity profiles are presented in sufficient detail to allow comparison with computational biofluid dynamics models. Wall shear stresses were estimated from the near wall velocity gradient, and the nature of observed poststenotic flow disturbances is discussed. Results indicate that flow disturbances of discrete oscillation frequency may be more valuable than turbulence as an indicator of early stages of stenosis development. Additionally, despite the fact that poststenotic turbulence exists for the higher degrees of stenosis and Reynolds numbers, the resulting wall shear stresses are only three to four times greater than the Poiseuille value and are considerably less than the wall shear stress within the stenosis itself.