The effect of compliance on wall shear in casts of a human aortic bifurcation

Rigid and compliant casts of a human aortic bifurcation were subjected to physiologically realistic pulsatile fluid flows. At a number of sites near the wall in the approximate median plane of the bifurcation of these models, fluid velocity was measured with a laser Doppler velocimeter, and wall motion (in the case of the compliant cast) was determined with a Reticon linescan camera. The velocity and wall motion data were combined to estimate the instantaneous shear rates at the cast wall. Analysis showed that at the outer walls the cast compliance reduced shear rates, while at the walls of the flow divider the shear rate was increased.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

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.     

Pulsatile flow in a constricted channel.

A nonuniform channel is used as a simple model of a constricted arterial vessel. Flow patterns have been calculated for pulsatile flow with both sinusoidal and nonsinusoidal flow rates for a range of Reynolds number, Re, and Strouhal number, St. The results show that even for relatively low frequency flows a strong vortex wave will be generated with a complex wall shear stress distribution and peak values much greater than those found in steady or unsteady parallel flow. The vortex wave increases in strength with increasing Re and St, with its total length and wavelength independent of Re but inversely proportional to St. The form of the imposed flow rate is found to have an important effect on the flow and the shear stress distribution.     

Experimental and numerical study of pulsatile flows through stenosis: wall shear stress analysis.

Different shapes of pulsatile flows through a model of stenosis are experimentally and numerically modeled to validate both methods and to determine the wall shear stress temporal evolution downstream from the stenosis. Two-dimensional velocity measurements are performed in a 75% severity stenosis using a pulsed Doppler ultrasonic velocimeter. Finite element package is employed for the transient numerical simulations. Polynomial method, based on the experimental velocity values, is proposed to determine the wall shear stress temporal evolution. There is a good agreement between the numerical and experimental results. The wall shear stress temporal analysis shows oscillating wall shear stress values during the cycle with high wall shear stress values at the throat of about 120 dyn/cm2, and low values downstream from the stenosis of about - 2.5 dyn/cm2. The key result of the study is that the presence of the stenosis leads the artery to work in a direction which is opposite to the direction of a healthy artery.     

Visualization and finite element analysis of pulsatile flow in models of the abdominal aortic aneurysm

Pulsatile flows in glass models simulating fusiform and lateral saccular aneurysms were investigated by a flow visualization method. When resting fluid starts to flow, the initial fluid motion is practically irrotational. After a short period of time, the flow began to separate from the proximal wall of the aneurysm. Then the separation bubble or vortex grew rapidly in size and filled the whole area of the aneurysm circumferentially. During this period of time, the center of the vortex moved from the proximal end to the distal point of the aneurysm. The transient reversal flow, for instance, which may occur at the end of the ejection period, passed between the wall of the aneurysm and the centrally located vortex. When the rate and pulsatile frequency of flow were high, the vortex broke down into highly disturbed flow (or turbulence) at the distal portion of the aneurysm. The same effect was observed when the length of the aneurysm was increased. A reduction in pulsatile amplitude made the flow pattern close to that in steady flow. A finite element analysis was made to obtain velocity and pressure fields in pulsatile flow through a tube with an axisymmetric expansion. Calculations were performed with the pulsatile flows used in the visualization experiment in order to study the effects of change in the pulsatile wave form by keeping the time-mean Reynolds number and Womersley's parameter unchanged. Calculated instantaneous patterns of velocity field and stream lines agreed well with the experimental results. The appearance and disappearance of the vortex in the dilated portion and its development resulted in complex distributions of pressure and shear fields. Locally minimum and maximum values of wall shear stress occurred at points just upstream and downstream of the distal end of the expansion when the flow rate reached its peak.     

Comparison of steady and pulsatile flow in a double branching arterial model

Data are presented to compare fluid flow parameters for steady flow with those for time-varying flow in a simplified two branch model which simulates the region of the abdominal aorta near the celiac and superior mesenteric branches of the dog. Measurements in the model included laser doppler anemometry velocity profiles during steady flow, sinusoidal flow with a superimposed mean flow (referred to as simple oscillatory flow) and arterial pulsatile flow. Shear rate measurements were made by an electrochemical technique during steady flow. Flow visualization studies were done during steady and pulsatile flow. Fluid flow effects in the simplified model during steady flow showed many similarities to the results from previous steady flow studies in a canine aortic cast. Shear rates in the region of the proximal (first, or celiac) branch were independent of flow rates in the distal (second, or mesenteric) branch, but the shear pattern within the proximal branch changed significantly as flow in the proximal branch increased. Shear rates on the proximal flow divider (leading edge into the distal branch) depended primarily on the flow rate to the proximal branch, but not on flow to the distal branch. At certain daughter branch flow ratios (approximately 2:1, proximal to distal), flow separation was promoted at the outer wall of the second branch, but flow separation did not occur in the first branch. In contrast to the canine aortic case results, flow separation was never detected on the distal (mesenteric) flow divider of the simplified model. This observation reflects the subtle effects of geometry on flow since the mesenteric flow divider in the canine cast protrudes into the main flow whereas the distal flow divider in the simplified model does not. There were distinct differences in the flow phenomena between steady, simple oscillatory and arterial pulsatile flow. Peak shear rates during pulsatile flow were as much as 10--100 times greater than steady flow shear rates at comparable mean flow rates. Particularly noteworthy for the pulsatile flow with a Womersley parameter of sixteen were very blunt velocity profiles throughout systole, and the absence of flow separation or reversal in those regions of the model that exhibited flow separation during steady flow. The shape of the waveform influences the nature of the flow during time-varying flows. Future studies of fluid dynamics in model systems must consider the pulsatile nature of the flow if a true interpretation of arterial flow phenomena is to be made.     

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.     

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.     

Pulsatile flow of Casson's fluid through stenosed arteries with applications to blood flow

The effects of non-Newtonian nature of blood and pulsatility on flow through a stenosed tube have been investigated. A perturbation method is used to analyse the flow. It is of interest to note that the thickness of the viscous flow region is non-uniform (changing with axial distance). An analytic relation between viscous flow region thickness and red cell concentration has been obtained. It is important to mention that some researchers have obtained an approximate solution for the flow rate-pressure gradient equation (assuming the ratio between the yield stress and the wall shear to be very small in comparison to unity); in the present analysis, we have obtained an exact solution for this non-linear equation without making that assumption. The approximate and exact solutions compare well with one of the exact solutions. Another important result is that the mean and steady flow rates decrease as the yield stress theta increases. For the low values of the yield stress, the mean flow rate is higher than the steady flow rate, but for high values of the yield stress, the mean flow rate behaviour is of opposite nature. The critical value of the yield stress at which the flow rate behaviour changes from one type to another has been determined. Further, it seems that there exists a value of the yield stress at which flow stops for both the flows (steady and pulsatile). It is observed that the flow stop yield value for pulsatile flow is lower than the steady flow. The most notable result of pulsatility is the phase lag between the pressure gradient and flow rate, which is further influenced by the yield stress and stenosis. Another important result of pulsatility is the mean resistance to flow is greater than its steady flow value, whereas the mean value of the wall shear for pulsatile flow is equal to steady wall shear. Many standard results regarding Casson and Newtonian fluids flow, uniform tube flow and steady flow can be obtained as the special cases of the present analysis. Finally, some applications of this theoretical analysis have been cited.     

Post-stenotic core flow behavior in pulsatile flow and its effects on wall shear stress

Arteries of several species, including man, tend to adjust their diameters such that the mean wall shear stress is in the range of 10-20 dynes cm-2. Additionally, intimal thickening in the human carotid bifurcation correlates well with the reciprocal of wall shear stress as determined in model studies. The correlation indicates that wherever the local mean wall shear stress exceeds approximately 10 dynes cm-2, the artery tends to be spared from intimal thickening. However, it is not known whether mean shear stress, i.e. the time-averaged value, or the instantaneous shear stress is the appropriate correlative variable. Each of these variables suggests different mechanisms for the reaction of the artery wall to its hemodynamic environment. It is therefore important to devise means by which the effects of mean shear and pulsatile shear can be separated in the study of atherogenesis. The present investigation examines the post-stenotic flow field in Plexiglas models under pulsatile conditions approximating those in the aortas of the cynomolgus monkey, an animal often employed in atherogenesis research. Behavior of the core flow and its effects on wall shear stress are studied for stenoses of 75 and 90% area reductions using laser velocimetry. The results show that the post-stenotic field contains regions in which the mean wall shear stress is low, but the pulsatile excursions are large.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Numerical simulations of pulsatile flow in an end-to-side anastomosis model

A potential interaction between the local hemodynamics and the artery wall response has been suggested for vascular graft failure by intimal hyperplasia (IH). Among the various hemodynamic factors, wall shear stress has been implicated as the primary factor responsible for the development of IH. In order to explore the role of hemodynamics in the formation of IH in end-to-side anastomosis, computational fluid dynamics is employed. To validate the numerical simulations, comparisons with existing experimental data are performed for both steady and pulsatile flows. Generally, good agreement is observed with the velocity profiles whereas some discrepancies are found in wall shear stress (WSS) distributions. Using the same end-to-side anastomosis geometry, numerical simulations are extended using a femoral artery waveform to identify the possible role of unsteady hemodynamics. In the current simulations, Carreau-Yasuda model is used to account for the non-Newtonian nature of blood. Computations indicated a disturbed flow field at the artery-graft junction leading to locally elevated shear stresses on the vascular wall. Furthermore, the shear stress distribution followed the same behavior with oscillating magnitude over the entire flow cycle. Thus, distal IH observed in end-to-side artery-graft models may be caused by the fluctuations in WSS's along the wall.     

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.     

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.     

Numerical investigation of the non-Newtonian pulsatile blood flow in a bifurcation model with a non-planar branch

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Pulsatile poststenotic flow studies with laser Doppler anemometry

The pulsatile flow field distal to axisymmetric constrictions in a straight tube was studied using laser Doppler anemometry. The upstream centerline velocity waveform was sinusoidal at a frequency parameter of 7.5 and mean Reynolds number of 600. Stenosis models of 25, 50 and 75% area reduction were employed and velocity data were derived by ensemble averaging methods. Extensive measurements of the pulsatile velocity profiles are reported, and wall shear rates were computed from the near wall velocity profile gradients. The experiments indicate that a permanent region of poststenotic flow separation does not exist even for the severest constriction, in contrast to results for steady flow. Values of wall shear stress were greatest near the throat of the constriction and were relatively low in the poststenotic region, including the region of most intense flow disturbance. Turbulence was found only for the 75% stenosis model and was created only during a segment of the cycle. Although much emphasis has been placed upon turbulence in the detection of arterial stenoses, particularly as identified by Doppler ultrasound spectral broadening, the present study implies that identification of flow disturbances of an organized nature may be more fundamental in recognizing mild to moderate disease. Additionally, the relationship of these flow field results to the animal aortic coarctation model often employed in atherogenesis studies is discussed.     

Unsteady and three-dimensional simulation of blood flow in the human aortic arch

A three-dimensional and pulsatile blood flow in a human aortic arch and its three major branches has been studied numerically for a peak Reynolds number of 2500 and a frequency (or Womersley) parameter of 10. The simulation geometry was derived from the three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo using CAT scan imaging on a human aorta. The numerical simulations were obtained using a projection method, and a finite-volume formulation of the Navier-Stokes equations was used on a system of overset grids. Our results demonstrate that the primary flow velocity is skewed towards the inner aortic wall in the ascending aorta, but this skewness shifts to the outer wall in the descending thoracic aorta. Within the arch branches, the flow velocities were skewed to the distal walls with flow reversal along the proximal walls. Extensive secondary flow motion was observed in the aorta, and the structure of these secondary flows was influenced considerably by the presence of the branches. Within the aorta, wall shear stresses were highly dynamic, but were generally high along the outer wall in the vicinity of the branches and low along the inner wall, particularly in the descending thoracic aorta. Within the branches, the shear stresses were considerably higher along the distal walls than along the proximal walls. Wall pressure was low along the inner aortic wall and high around the branches and along the outer wall in the ascending thoracic aorta. Comparison of our numerical results with the localization of early atherosclerotic lesions broadly suggests preferential development of these lesions in regions of extrema (either maxima or minima) in wall shear stress and pressure.     

Numerical analysis of flow in an elastic artery model.

Oscillatory and pulsatile flows of Newtonian fluids in straight elastic tubes are simulated numerically with the aid of Ling and Atabek's "local flow" assumption for the nonlinear convective acceleration terms. For the first time, a theoretical assessment of the local flow assumption is presented, and the range of validity of the assumption is estimated by comparison with perturbation solutions of the complete flow problem. Subsequent simulations with the local flow model indicate that the flow field and associated wall shear stress are extremely sensitive to the phase angle between oscillatory pressure and flow waves (impedance phase angle). This phase angle, which is a measure of the wave reflection present in the system, is known to be altered by arterial disease (e.g., hypertension) and vasoactive drugs. Thus, the paper elucidates a mechanism by which subtle changes in systemic hemodynamics (i.e., phase angles) can markedly influence local wall shear stress values.