Experimental technique of measuring dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact cause and mechanism of the progression of AV calcification is unknown, although mechanical forces have been known to play a role. It is thus important to characterize the mechanical environment of the AV. In the current study, we establish a methodology of measuring shear stresses experienced by the aortic surface of the AV leaflets using an in vitro valve model and adapting the laser Doppler velocimetry (LDV) technique. The valve model was constructed from a fresh porcine aortic valve, which was trimmed and sutured onto a plastic stented ring, and inserted into an idealized three-lobed sinus acrylic chamber. Valve leaflet location was measured by obtaining the location of highest back-scattered LDV laser light intensity. The technique of performing LDV measurements near to biological surfaces as well as the leaflet locating technique was first validated in two phantom flow systems: (1) steady flow within a straight tube with AV leaflet adhered to the wall, and (2) steady flow within the actual valve model. Dynamic shear stresses were then obtained by applying the techniques on the valve model in a physiologic pulsatile flow loop. Results show that aortic surface shear stresses are low during early systole (<5 dyn/cm2) but elevated to its peak during mid to late systole at about 18-20 dyn/cm2. Low magnitude shear stress (<5 dyn/cm2) was observed during early diastole and dissipated to zero over the diastolic duration. Systolic shear stress was observed to elevate only with the formation of sinus vortex flow. The presented technique can also be used on other in vitro valve models such as congenitally geometrically malformed valves, or to investigate effects of hemodynamics on valve shear stress. Shear stress data can be used for further experiments investigating effects of fluid shear stress on valve biology, for conditioning tissue engineered AV, and to validate numerical simulations.     

Three-dimensional visualization of velocity fields downstream of six mechanical aortic valves in a pulsatile flow model

Velocity fields downstream of 27 mm Bj?rk-Shiley Standard, Bj?rk-Shiley Convex-Concave, Bj?rk-Shiley Monostrut, Hall-Kaster (Medtronic-Hall), St. Jude Medical and Starr-Edwards Silastic Ball aortic valves were studied in a pulsatile mock circulation. Stroke volume was 70 cm3 and frequency 71 min-1 and 88 min-1. Fluid velocity was measured by a catheter mounted hot-film anemometer probe in a glycerol water mixture one and two diameters downstream of the aortic valve. Velocity fields were dynamically visualized by a three-dimensional technique and revealed qualitative independence of frequency. All profiles were flat in the acceleration phase of systole. From peak systole and throughout the systolic deceleration phase profiles characteristic of the individual valves appeared. The pivoting and tilting disc valves caused a skewed velocity profile with highest velocities downstream of the major orifice and lowest velocities downstream of the minor orifice. The differences between the three investigated Bj?rk-Shiley valves were remarkable. The St. Jude Medical valve generated velocity peaks downstream of the two major orifices and the central slit, and lower velocities in the hinge areas. A rather flat profile with central hollowing was seen downstream of the Starr-Edwards Ball valve. All velocity profiles were more or less dampened two diameters downstream.     

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.     

The impedance of curved artery models

The impedance (pressure drop/flow rate) of four curved artery models has been determined experimentally for steady and periodic flows simulating conditions in the aortic arch. Steady flow results indicate that very short entry lengths are required for flow development in curved artery models, and impedance is elevated above straight tube values by a factor of 3-4 for mean flow conditions in the aortic arch. Results for periodic flow with a nonzero mean show a significant elevation of mean flow impedance relative to values for steady flow at the mean flow rate--a factor of 2-3 for aortic arch flow conditions. The impedance of the first harmonic of periodic flows follows straight tube theory at high values of the unsteadiness parameter in agreement with available theory for curved tubes. The implications of the impedance measurements for wall shear stress in the aortic arch are discussed.     

Platelet activation of mechanical versus bioprosthetic heart valves during systole

Thrombus formation is a major concern for recipients of mechanical heart valves (MHVs), which requires them to take anticoagulant drugs for the rest of their lives. Bioprosthetic heart valves (BHVs) do not require life-long anticoagulant therapy but deteriorate after 10–15 years. The thrombus formation is initiated by the platelet activation which is thought to be mainly generated in MHVs by the flow through the hinge and the leakage flow during the diastole. However, our results show that the activation in the bulk flow during the systole phase might play an essential role as well. This is based on our results obtained by comparing the thrombogenic performance of a MHV and a BHV (as control) in terms of shear induced platelet activation under exactly the same conditions. Three different mathematical activation models including linear level of activation, damage accumulation, and Soares model are tested to quantify the platelet activation during systole using the previous simulations of the flow through MHV and BHV in a straight aorta under the same physiologic flow conditions. Results indicate that the platelet activation in the MHV at the beginning of the systole phase is slightly less than the BHV. However, at the end of the systole phase the platelet activation by the bulk flow for the MHV is several folds (1.41, 5.12, and 2.81 for linear level of activation, damage accumulation, and Soares model, respectively) higher than the BHV for all tested platelet activation models.     

Computational fluid dynamics modeling and analysis of the effect of 3-D distortion of the human aortic arch

An idealized CFD model and a realistic one were used to investigate the effect of the 3-D distortion of the aortic arch on the blood flow and its pathophysiological significance with respect to the pathogenesis of the aortic aneurysm. From the results of the flow simulations, the distortion of the centerline of the pipe was shown to affect significantly the flow structure. A right-handed vortex at the descending arch, and a left-handed one at the end of the arch tended to develop in the realistic model. But the secondary flow did not become a single helix. The top of the arch was the region where complex spatial and temporal WSS distributed. It was also observed that the direction of WSS had a significant circumferential component at the top of the arch.     

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.     

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.     

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.     

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.     

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.     

RIGHT SINUS OF VALSALVA ANEURYSM RUPTURE INTO THE RIGHT ATRIUM: PRE- AND POSTOPERATIVE ECHOCARDIOGRAPHIC FINDINGS

Echocardiographic features of a patient with sinus of Valsalva aneurysm rupture into the right atrium are discussed. An abnormal echo representing the aneurysmal sac was seen in the anterior aortic root. Other abnormal echoes were demonstrated in the right atrium. Flutter was visible in systole and diastole in the tricuspid area. The pulmonic valve echogram was normal. After surgical repair, the above findings disappeared.     

Human sperm swimming in flow

To study the movement of human sperm, we have developed a microflow cell by miniaturizing our design for a preparative fractionation flow column. The microflow cell enabled us to view the movement of sperm over periods as long as 2 min. Sequential steps of filming, editing, and analysis revealed that the curved swimming patterns of sperm swimming in stagnant fluid become nearly straight tracks when the flow velocity is increased. However, the net swimming speed remained unchanged. Motile sperm accumulated near solid wall surfaces surrounding the fluid and oriented against the direction of the current; the velocity gradient was steepest in these regions. A laminar-flow preparative column separated motile sperm from dead sperm by carrying the nonmotile sperm and debris with the stream while leaving the motile sperm near the surrounding walls.     

Human sperm swimming in flow

Abstract. To study the movement of human sperm, we have developed a microflow cell by miniaturizing our design for a preparative fractionation flow column. The microflow cell enabled us to view the movement of sperm over periods as long as 2 min. Sequential steps of filming, editing, and analysis revealed that the curved swimming patterns of sperm swimming in stagnant fluid become nearly straight tracks when the flow velocity is increased. However, the net swimming speed remained unchanged. Motile sperm accumulated near solid wall surfaces surrounding the fluid and oriented against the direction of the current; the velocity gradient was steepest in these regions. A laminar-flow preparative column separated motile sperm from dead sperm by carrying the nonmotile sperm and debris with the stream while leaving the motile sperm near the surrounding walls.     

Regional epicardial strain in the embryonic chick heart during the early looping stages

Epicardial strains were measured in Hamburger-Hamilton stage 11 and 12 embryonic chick hearts (1.6-2.0 days of incubation). These stages include part of the early phase of cardiac looping, as the initially straight heart tube bends and twists to form a curved c-shaped tube. By analyzing the motion of microbeads placed on the myocardial surface, we measured strains near the outer curvature, in the central region, and near the inner curvature of the primitive ventricle. No significant differences in strain were found between stages. Relative to end diastole, all three regions shortened by about 10% during systole in the circumferential direction, and the outer curvature shortened longitudinally by about 5%. In contrast, and unlike strains in older hearts, the inner curvature and central regions elongated by approximately 5-10% in the longitudinal direction during systole. These results are consistent with microstructural data and suggest that the material properties of the outer curvature are relatively isotropic, whereas the properties of the central and inner curvature regions are orthotropic, with contractile stress exerted primarily in the circumferential direction.     

Helical flow around arterial bends for varying body mass

The three dimensionally curved aortic arch is modeled as a portion of a helical pipe. Pulsatile blood flow therein is calculated assuming helical symmetry and an experimentally measured pressure pulse. Appropriate values for the Womersley and Reynolds numbers are taken from allometric scaling relations for a variety of body masses. The flow structure is discussed with particular reference to the wall shear, which is believed to be important in the inhibition of atheroma. It is found that nonplanar curvature limits the severity of flow separation at the inner bend, and reduces spatial variation of wall shear.     

Sheep cardiac Purkinje fibers: Configurational changes during the cardiac cycle

Previous attempts to study the cytoarchitecture of cardiac Purkinje fibers with the scanning electron microscope (SEM) have been limited by the surrounding dense connective tissue. In this study the connective tissue was removed by treatment with 8N HCl, after adult sheep hearts were fixed in diastole or systole and tissue taken for SEM and transmission electron microscopy (TEM). In SEM, Purkinje fibers freely anastomosed in false tendons and formed a subendocardial plexus. In systole, medium and small-sized Purkinje fibers formed deep clefts not observed in diastole. The clefts are thought to be due to sarcolemmal folding and fiber buckling and may therefore affect conduction. The myofibrils beneath the laterally apposed sarcolemmas of adjacent Purkinje cells when fixed in systole were often observed as tightly curved arches in series. Similar configurations with expanded arches were observed in diastole. The formation of arches by myofibrils is unique to Purkinje fibers and is interpreted as the mechanism responsible for their compliance to stretch. The significance of contraction in producing the observed geometric changes in Purkinje fibers and the implications of their cytoarchitecture with respect to conduction are discussed.     

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.     

Flow patterns in three-dimensional left ventricular systolic and diastolic flows determined from computational fluid dynamics

A realistic model of the left ventricle of the heart was previously constructed, using a cast from a dog heart which was in diastole. Previous studies of the three-dimensional heart model were conducted in systole only. The purpose of this investigation was to extend the model to both systole and diastole, and to determine what the effect of a previous cardiac cycle was on the next cardiac cycle. The 25.8 cc ventricular volume was reduced by 40% in 0.25 seconds, then increased to the original volume in another 0.25 seconds and then allowed to rest for 0.25 seconds. Runs done with an ejection fraction of 60% showed little variation from one cardiac cycle to another after the third cardiac cycle was completed; the maximum velocity could vary by over 30% between the first and second cardiac cycles. In systole, centerline and cross-sectional velocity vectors greatly increased in magnitude at the aortic outlet. Most of the pressure drop occurred in the top 15% of the heart. The diastolic phase showed complex vortex formation not seen in the systolic contractions; these complex vortices could account for experimentally observed turbulent blood flow fluctuations in the aorta.