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.     

Comparative study of the function of the Abiomed polyurethane heart valve for use in left ventricular assist devices

Hydrodynamic testing of the Abiomed polyurethane trileaflet valve has been carried out to establish performance data of valve function. A Medtronic Hall tilting disk, a Carbomedics bileaflet, a Hancock II bioprosthesis and an Abiomed polyurethane trileaflet valve, all size 27 mm, underwent both pulsatile and steady-flow hydrodynamic testing. Results of the variation of pressure difference with RMS pulsatile flow and steady flow, and effective orifice area, showed that the Abiomed valve had significantly poorer opening characteristics than the tissue valve and the two mechanicalvalves. The Abiomed valve's performance was seen to be related to its construction and manufacture. This study highlights some of the problems associated with the design and development of synthetic trileaflet heart valve prostheses.     

Velocity and turbulence measurements past mitral valve prostheses in a model left ventricle

Thrombogenesis and hemolysis have both been linked to the flow dynamics past heart valve prostheses. To learn more about the particular flow dynamics past mitral valve prostheses in the left ventricle under controlled experimental conditions, an in vitro study was performed. The experimental methods included velocity and turbulent shear stress measurements past caged-ball, tilting disc, bileaflet, and polyurethane trileaflet mitral valves in an acrylic rigid model of the left ventricle using laser Doppler anemometry. The results indicate that all four prosthetic heart valves studied create at least mildly disturbed flow fields. The effect of the left ventricular geometry on the flow development is to produce a stabilizing vortex which engulfs the entire left ventricular cavity, depending on the orientation of the valve. The measured turbulent shear stress magnitudes for all four valves did not exceed the reported value for hemolytic damage. However, the measured turbulent shear stresses were near or exceeded the critical shear stress reported in the literature for platelet lysis, a known precursor to thrombus formation.     

Numerical simulation of steady flow in a two-dimensional total artificial heart model.

In this paper, a numerical simulation of steady laminar and turbulent flow in a two-dimensional model for the total artificial heart is presented. A trileaflet polyurethane valve was simulated at the outflow orifice while the inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-epsilon-E model. The SIMPLER algorithm was used to solve the problem in primitive variables. The numerical solutions of the simulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller with the flap valve than with the trileaflet valve. These results also suggest a correlation between high turbulent stresses and the presence of thrombus in the vicinity of the valves in the total artificial hearts. The computed velocity vectors and turbulent stresses were comparable with previously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solutions suggests that geometries similar to the flap valve (or a tilting disk valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.     

Numerical simulation of steady turbulent flow through trileaflet aortic heart valves—II. Results on five models

Turbulent flow simulations are run for five aortic trileaflet valve geometries, ranging from a valve leaflet orifice area of 1.1 cm² (Model A1—very stenotic) to 5.0 cm² (Model A5—natural valve). The simulated data compares well with experimental measurements made downstream of various aortic trileaflet valves by Woo (PhD Thesis, 1984). The location and approximate width and length of recirculation regions are correctly predicted. The less stenotic valve models reattach at the end of the aortic sinus region, 1.1 diameters downstream of the valve. The central jet exiting the less stenotic valve models is not significantly different from fully developed flow, and therefore recovers very quickly downstream of the reattachment point. The more stenotic valves disturb the flow to a greater degree, generating recirculation regions large enough to escape the sinuses and reattach further downstream. Peak turbulent shear stress values downstream of the aortic valve models which approximated prosthetic valves are 125 and 300 N m⁻², very near experimental observations of 150 to 350 N m⁻². The predicted Reynolds stress profiles also present the correct shape, a double peak profile, with the location of the peak occuring at the location of maximum velocity gradient, which occurs near the recirculation region. The pressure drop across model A2 (leaflet orifice area 1.6 cm²) is 20 mmHg at 1.6 diameters downstream. This compares well with values ranging from 19.5 to 26.2 mmHg for valves of similar orifice areas. The pressure drop decreases with decreasing valve stenosis, to a negligible value across the least stenotic valve model. Based on the good agreement between experimental measurements of velocity, shear stress and pressure drop, compared to the simulated data, the model has the potential to be a valuable tool in the analysis of heart valve designs.     

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.     

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.     

Comparison of the hemodynamic and thrombogenic performance of two bileaflet mechanical heart valves using a CFD/FSI model

The hemodynamic and the thrombogenic performance of two commercially available bileaflet mechanical heart valves (MHVs)--the ATS Open Pivot Valve (ATS) and the St. Jude Regent Valve (SJM), was compared using a state of the art computational fluid dynamics-fluid structure interaction (CFD-FSI) methodology. A transient simulation of the ATS and SJM valves was conducted in a three-dimensional model geometry of a straight conduit with sudden expansion distal the valves, including the valve housing and detailed hinge geometry. An aortic flow waveform (60 beats/min, cardiac output 4 l/min) was applied at the inlet. The FSI formulation utilized a fully implicit coupling procedure using a separate solver for the fluid problem (FLUENT) and for the structural problem. Valve leaflet excursion and pressure differences were calculated, as well as shear stress on the leaflets and accumulated shear stress on particles released during both forward and backward flow phases through the open and closed valve, respectively. In contrast to the SJM, the ATS valve opened to less than maximal opening angle. Nevertheless, maximal and mean pressure gradients and velocity patterns through the valve orifices were comparable. Platelet stress accumulation during forward flow indicated that no platelets experienced a stress accumulation higher than 35 dyne x s/cm2, the threshold for platelet activation (Hellums criterion). However, during the regurgitation flow phase, 0.81% of the platelets in the SJM valve experienced a stress accumulation higher than 35 dyne x s/cm2, compared with 0.63% for the ATS valve. The numerical results indicate that the designs of the ATS and SJM valves, which differ mostly in their hinge mechanism, lead to different potential for platelet activation, especially during the regurgitation phase. This numerical methodology can be used to assess the effects of design parameters on the flow induced thrombogenic potential of blood recirculating devices.     

Techniques of aortic valve repair

Similar to mitral repair, newer methods of aortic valve reconstruction are achieving excellent outcomes with an 85% to 90% freedom from valve-related complications at 10 years. The goal of this review is to illustrate these newer and more stable techniques of aortic valve repair. Most patients with aortic insufficiency from either trileaflet or bicuspid aortic valves are candidates for repair, in addition to selected patients with mixed aortic stenosis/insufficiency and aortic root aneurysms. Initially, aggressive commissural annuloplasty is performed to reduce measured valve diameter to 19 to 21 mm. Leaflet prolapse is corrected with plication stitches placed in the free edge of each leaflet adjacent to the Nodulus Arantius. In this regard, the leaflet free edge functions as the chorda tendinea of the aortic valve, and shortening with plication stitches raises the leaflet to a proper "effective height." Leaflet defects are augmented with gluteraldehyde-fixed autologous pericardium, and mild-to-moderate strategically placed spicules of calcium are removed with the cavitron ultrasonic surgical aspirator. Using these methods, most insufficient aortic valves, and many with mixed lesions, can be satisfactorily repaired. Six cases are illustrated in this review, spanning the spectrum of pathologies from annular dilatation without leaflet defects, to standard congenital bicuspid valve with prolapse, to trileaflet prolapse, to unusual bicuspid pathology with calcification, to a moderately calcified trileaflet valve with mixed lesions, and to aortic root aneurysms with severe aortic insufficiency. All valves were repaired using the techniques described above with trivial residual leak and minimal gradients. All repairs have been followed with yearly echocardiography, and valve reconstruction with these methods is now quite stable with excellent late outcomes. Most insufficient aortic valves now can undergo stable repair with minimal late valve-related complications. Greater application of aortic valve repair seems indicated.     

Three-dimensional visualization of stationary flow behind artificial heart valves]

The visualization and quantitative analysis of flow offers a possibility for the hydrodynamic characterization of artificial heart valves. Different types of valves can be compared if velocity profile and the turbulent shear stress caused by the prosthesis are known. The tracer technique was selected, since it permits visualization also of turbulent flow through the valve. With the aid of a simple optical device the three-dimensional flow pattern behind the valve is determinable. The main features of the method are: The regions of interest can easily be identified. Velocity profiles can be determined and shear stress and turbulence intensities estimated. The experimental setup is simple, calibration is not necessary, and it can be used for turbulent flows. The method can be used only with transparent fluids and vessels; measurements in blood are not possible. Because of the large number of measuring points required the method is very time-consuming. The use of an automatic picture analyzing system would make it possible to increase the number of pictures processed, and thus increase resolution. The velocity profile of a three-finger-valve, the TAD 29, was established at a distance of 20 mm from the ring, and compared with known profiles from the literature. The valve has an opening angle of 70 degrees. All typical regions for the flow of an artificial heart valve, such as jet, stagnation gone, backflow and turbulence were demonstrated.     

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.     

Hydrodynamic in vitro comparison of a new artificial heart valve with the Bj?rk-Shiley valve (standard)

Measurements performed to compare a newly developed tilting disc valve with the Bj?rk-Shiley valve included velocity profiles downstream of the heart valves, valve-induced flow turbulence and pressure drop across the opened valves. The velocity profiles measured with pulsed Doppler ultrasound are similar, although they do not permit a quantitative comparison of the valves. The interpretation of the 90 degrees-component of Doppler signals as a measure of the turbulence permits a quantitative comparison without the need for extensive measurements. However, only large vortices are recorded, so that our turbulent shear stresses are lower than these reported in the literature. The pressure drop across the opened valve is a measure of the energy loss, and important parameters for the valve can be derived from it. The pressure drop is dependent on the test conditions, and is therefore not a characteristic constant of the valve. The transformation of the power law Q = C delta P beta into a relation between Re- and Eu-number gives a nondimensional similarity number that is characteristic for tilting disc valves. Its verification requires more investigations, involving variation of valve size and the viscosity of the test fluid.     

Probing the effect of morphology on lymphatic valve dynamic function

The lymphatic system is vital to the circulatory and immune systems, performing a range of important functions such as transport of interstitial fluid, fatty acid, and immune cells. Lymphatic vessels are composed of contractile walls and lymphatic valves, allowing them to pump lymph against adverse pressure gradients and to prevent backflow. Despite the importance of the lymphatic system, the contribution of mechanical and geometric changes of lymphatic valves and vessels in pathologies of lymphatic dysfunction, such as lymphedema, is not well understood. We develop a fully coupled fluid–solid, three-dimensional computational model to interrogate the various parameters thought to influence valve behavior and the consequences of these changes to overall lymphatic function. A lattice Boltzmann model is used to simulate the lymph, while a lattice spring model is used to model the mechanics of lymphatic valves. Lymphatic valve functions such as enabling lymph flow and preventing backflow under varied lymphatic valve geometries and mechanical properties are investigated to provide an understanding of the function of lymphatic vessels and valves. The simulations indicate that lymphatic valve function is optimized when valves are of low aspect ratio and bending stiffness, so long as these parameters are maintained at high enough values to allow for proper valve closing. This suggests that valve stiffening could have a profound effect on overall lymphatic pumping performance. Furthermore, dynamic valve simulations showed that this model captures the delayed response of lymphatic valves to dynamic flow conditions, which is an essential feature of valve operation. Thus, our model enhances our understanding of how lymphatic pathologies, specifically those exhibiting abnormal valve morphologies such as has been suggested to occur in cases of primary lymphedema, can lead to lymphatic dysfunctions.     

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.     

Effect of systolic flow rate on the prediction of effective prosthetic valve orifice area

The Gorlin equation for the hemodynamic assessment of valve area is commonly used in cardiac catheterization laboratories. A study was performed to test the prediction capabilities of the Gorlin formula as well as the Aaslid and Gabbay formula for the effective orifice area of prosthetic heart valves. Pressure gradient, flow, and valve opening area measurements were performed on four 27 mm valve prostheses (two mechanical bileaflet designs, St. Jude and Edwards-Duromedics, an Edwards pericardial tissue valve, and a trileaflet polyurethane valve) each mounted in the aortic position of an in vitro pulse duplicator. With the known valve orifice area, a different discharge coefficient was computed for each of the four valves and three orifice area formulas. After some theoretical considerations, it was proposed that the discharge coefficient would be a function of the flow rate through the valve. All discharge coefficients were observed to increase with increasing systolic flow rate. An empirical relationship of discharge coefficient as a linear function of systolic flow rate was determined through a regression analysis, with a different relationship for each valve and each orifice area formula. Using this relationship in the orifice area formulas improved the accuracy of the prediction of the effective orifice area with all three formulas performing equally well.     

Numerical simulation of steady turbulent flow through trileaflet aortic heart valves—I. Computational scheme and methodology

A numerical simulation model and technique are described to simulate steady turbulent blood flow through trileaflet tissue valves of varying degrees of stenosis. The aortic trileaflet tissue valve design was chosen as the subject of this study, since it is the only popular valve in current clinical use which is approximately axisymmetric. An axisymmetric geometry is computationally more convenient since it involves only two dimensional equations. The geometry and dimensions of the aorta were designed from angiographic studies and measurements made from cadavers. The valve dimensions were obtained from tissue leaflet photography studies conducted on tissue bioprostheses of varying degrees of stenosis. The nonrectangular nature of this valve necessitated the use of a body conforming grid. Thompson's method coupled with a Chimera grid system was chosen for this purpose. The Chimera grid was used to avoid a grid with highly skewed cells. Turbulence was simulated by using the k- model with the wall function method. This decision was made after comparing the k- model's performance with that of lower order models, and after considering the increased computer time requirements and decreased stability of more complex models, such as the Reynolds stress model. The results of the study which are very encouraging and compare favorably with in vitro experimental data, are described in Part II of the paper.     

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.     

An in vitro comparative study of St. Jude Medical and Edwards-Duromedics bileaflet valves using laser anemometry

An in vitro comparative study of St. Jude (SJ) and Edwards-Duromedics (DM) Bileaflet valves was performed under steady and physiological pulsatile flow conditions in an axisymmetric chamber using Laser Doppler Anemometry (LDA). LDA measurements were conducted in two different orientations; in the first orientation, the LDA traverse was perpendicular and, in the second orientation, parallel to the tilt axis of the leaflets. The axial velocities were measured in both orientations at two different locations distal to the valves. The velocity profiles at peak systole show the presence of stronger vortex in the sinus region for flow past SJ valve in the first orientation compared to the DM valve. Velocity profile distal to the SJ valve in second orientation was relatively flat where as for the DM valve, a jet-like flow was present. The differences found in the velocity profiles between the two valves can be attributed to the differences in geometry with thicker leaflets, smaller angle of leaflets opening and the presence of the leaflet curvature for the DM valve. The results obtained in this study do not show any fluid dynamic advantages due to the curved leaflet geometry of the DM valve.     

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.     

Geometric modeling of functional trileaflet aortic valves: development and clinical applications

The dimensions of the aortic valve components condition its ability to prevent blood from flowing back into the heart. While the theoretical parameters for best trileaflet valve performance have already been established, an effective approach to describe other less optimal, but functional models has been lacking. Our goal was to establish a method to determine by how much the dimensions of the aortic valve components can vary while still maintaining proper function. Measurements were made on silicone rubber casts of human aortic valves to document the range of dimensional variability encountered in normal adult valves. Analytical equations were written to describe a fully three-dimensional geometric model of a trileaflet valve in both the open and closed positions. A complete set of analytical, numerical and graphical tools was developed to explore a range of component dimensions within functional aortic valves. A list of geometric guidelines was established to ensure safe operation of the valve during the cardiac cycle, with practical safety margins. The geometry-based model presented here allows determining quickly if a certain set of valve component dimensions results in a functional valve. This is of great interest to designers of new prosthetic heart valve models, as well as to surgeons involved in valve-sparing surgery.