Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood

Pulsatile flow in an axisymmetric rigid-walled model of an abdominal aorta aneurysm was analyzed numerically for various aneurysm dilations using physiologically realistic resting waveform at time-averaged Reynolds number of 300 and peak Reynolds number of 1607. Discretization of the governing equations was achieved using a finite element scheme based on the Galerkin method of weighted residuals. Comparisons with previously published work on the basis of special cases were performed and found to be in excellent agreement. Our findings indicate that the velocity fields are significantly affected by non-Newtonian properties in pathologically altered configurations. Non-Newtonian fluid shear stress is found to be greater than Newtonian fluid shear stress during peak systole. Further, the maximum shear stress is found to occur near the distal end of AAA during peak systole. The impact of non-Newtonian blood flow characteristics on pressure compared to Newtonian model is found insignificant under resting conditions. Viscous and inertial forces associated with blood flow are responsible for the changes in the wall that result in thrombus deposition and dilation while rupture of AAA is more likely determined by much larger mechanical stresses imposed by pulsatile pressure on the wall of AAA.     

Blood flow in abdominal aortic aneurysms: pulsatile flow hemodynamics

Numerical predictions of blood flow patterns and hemodynamic stresses in Abdominal Aortic Aneurysms (AAAs) are performed in a two-aneurysm, axisymmetric, rigid wall model using the spectral element method. Physiologically realistic aortic blood flow is simulated under pulsatile conditions for the range of time-averaged Reynolds numbers 50< or =Re(m)< or =300, corresponding to a range of peak Reynolds numbers 262.5< or =Re(peak) < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is characterized by a sequence of five different flow phases in one period of the flow cycle. Hemodynamic disturbance is evaluated for a modified set of indicator functions, which include wall pressure (p(w)), wall shear stress (tau(w)), and Wall Shear Stress Gradient (WSSG). At peak flow, the highest shear stress and WSSG levels are obtained downstream of both aneurysms, in a pattern similar to that of steady flow. Maximum values of wall shear stresses and wall shear stress gradients obtained at peak flow are evaluated as a function of the time-average Reynolds number resulting in a fourth order polynomial correlation. A comparison between predictions for steady and pulsatile flow is presented, illustrating the importance of considering time-dependent flow for the evaluation of hemodynamic indicators.     

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.     

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 simulation of sinusoidal flow in a straight elastic tube: effects of phase angles

Numerical simulations of flow in straight elastic (moving wall) tubes subjected to a sinusoidal pressure gradient were performed for conditions prevailing in large and medium sized arteries. The effects of varying the phase angle between the pressure gradient and the tube radius, the amplitude of wall motion, and the unsteadiness parameter (alpha) on flow rate and wall shear stress were investigated. Mean and peak flow rates and shear stresses were found to be strongly affected by the phase angle between the pressure gradient and the tube radius with greater sensitivity at higher diameter variation and higher alpha. In large artery simulations (alpha = 12), means flow rate was found to be 60% higher and peak flow rate to be 73% higher than corresponding rigid tube values for certain phase angles, while a threefold increase in mean wall shear stress and sevenfold increase in peak wall shear stress were observed in a sensitive phase angle range. Significant reversal in the wall shear stress direction occurred in the sensitive phase angle range even when there was negligible flow rate reversal. All effects were greatly diminished in simulations of medium sized vessels (alpha = 4). Some experimental evidence to support the predictions of a strong effect of phase angle on wall shear stress in large vessels is presented. Finally, physiological implications of the present work are discussed from a basis of aortic input impedance data, and a physical explanation for the extreme sensitivity of the flow field to small amplitude wall motion at high alpha is given.     

Velocity and wall shear stress patterns in the human right coronary artery.

Blood flow dynamics in the human right coronary artery have not been adequately quantified despite the clinical significance of coronary atherosclerosis. In this study, a technique was developed to construct a rigid flow model from a cast of a human right coronary artery. A laser photochromic method was used to characterize the velocity and wall shear stress patterns. The flow conditions include steady flow at Reynolds numbers of 500 and 1000 as well as unsteady flow with Womersley parameter and peak Reynolds number of 1.82 and 750, respectively. Characterization of the three-dimensional geometry of the artery revealed that the largest spatial variation in curvature occurred within the almost branch-free proximal region, with the greatest curvature existing along the acute margin of the heart. In the proximal segment, high shear stresses were observed on the outer wall and lower, but not negative, stresses along the inner wall. Low shear stress on the inner wall may be related to the preferential localization of atherosclerosis in the proximal segment of the right coronary artery. However, it is possible that the large difference between the outer and inner wall shear stresses may also be involved.     

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.     

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.     

A two-dimensional numerical analysis of unsteady flow in the carotid artery bifurcation. A comparison with three-dimensional in-vitro measurements and the influence of minor stenoses

In the present study a two-dimensional finite element model for incompressible Newtonian flow is applicated to the modelling of carotid artery flow. In earlier studies, the numerical model was validated experimentally for several flow configurations. In general the pulsatile flow is characterized by reversed flow regions at the non-divider side walls of both the internal and external carotid arteries. The unsteadiness of the flow is associated with rather complex spatial and temporal velocity distributions and leads to temporal variations of the location and length of the reversed flow regions. As a consequence, pronounced spatial and temporal variations in the wall shear stresses are found. At the non-divider side walls, wall shear stresses are relatively low and exhibits an oscillatory behaviour in space and time. At the divider side walls, wall shear stresses are relatively high and approximately follow the flow rate distribution in time. The aim of this study is not only to present two-dimensional calculations but also to compare the calculated two-dimensional velocity profiles with those from three-dimensional experiments. It is observed that in the common carotid artery and in the proximal parts of the internal and external carotid arteries, the two-dimensional numerical model provides valuable information with respect to the three-dimensional configuration. In the more distal parts of especially the internal carotid artery, deviations are found between the two-dimensional numerical and three-dimensional experimental model. These deviations can mainly be attributed to the neglect of the secondary velocity distribution in the two-dimensional model. In the two-dimensional numerical model the influence of a minor stenosis in the internal carotid artery is hardly distinguishable from a minor geometrical variation without stenosis. Full three-dimensional analyses of the influence of minor stenoses are needed to prove numerically whether in-vivo measurements of the axial velocity distribution are useful in the detection of minor stenoses.     

Intra-aneurysmal flow with helix and mesh stent placement across side-wall aneurysm pore of a straight parent vessel

Pulsatile flow fields in a cerebrovascular side-wall aneurysm model with a wide ostium after stenting are presented in terms of particle tracking velocimetry measurements and flow visualization. Among the stent parameters the shape, helix versus mesh, was selected to study its effect on the changes of intraaneurysmal hemodynamics for the reference of minimally invasive endovascular aneurysm treatment. The blocking ratio of the stents was fixed at 30%. The Womersley number was 3.9 and the mean, peak, and minimal Reynolds numbers based on the bulk average velocity and diameter of the parent vessel were 600, 850, and 300, respectively. Four consecutive flow-rate phases were selected to characterize the intra-aneurysmal flow. The results are characterized in terms of velocity vector field, regional average velocity, and intra-aneurysmal vorticity/circulation/wall shear stress. It is found that the hemodynamic features inside the aneurysm alter markedly with the shape of the stent and the size of the orifice. Both stents investigated induce favorable changes in the intra-aneurysmal flow stasis as well as direction and undulation of wall shear stresses. A comparison of the results of the helix to mesh stent shows that the former is more favorable for endovascular treatment.     

The fluid shear stress distribution on the membrane of leukocytes in the microcirculation

Recent in-vivo and in-vitro evidence indicates that fluid shear stress on the membrane of leukocytes has a powerful control over several aspects of their cell function. This evidence raises a question about the magnitude of the fluid shear stress on leukocytes in the circulation. The flow of plasma on the surface of a leukocyte at a very low Reynolds number is governed by the Stokes equation for the motion of a Newtonian fluid. We numerically estimated the distribution of fluid shear stress on a leukocyte membrane in a microvessel for the cases when the leukocyte is freely suspended, as well as rolling along or attached to a microvessel wall. The results indicate that the fluid shear stress distribution on the leukocyte membrane is nonuniform with a sharp increase when the leukocyte makes membrane attachment to the microvessel wall. In a microvessel (10 microns diameter), the fluid shear stress on the membrane of a freely suspended leukocyte (8 microns diameter) is estimated to be several times larger than the wall shear stress exerted by the undisturbed Poiseuille flow, and increases on an adherent leukocyte up to ten times. High temporal stress gradients are present in freely suspended leukocytes in shear flow due to cell rotation, which are proportional to the local shear rate. In comparison, the temporal stress gradients are reduced on the membrane of leukocytes that are rolling or firmly adhered to the endothelium. High temporal gradients of shear stress are also present on the endothelial wall. At a plasma viscosity of 1 cPoise, the peak shear stresses for suspended and adherent leukocytes are of the order of 10 dyn/cm2 and 100 dyn/cm2, respectively.     

Physical stresses at the air-wall interface of the human nasal cavity during breathing.

The nose is the front line defender of the respiratory system and is rich with mechanoreceptors, thermoreceptors, and nerve endings. A time-dependent computational model of transport through nasal models of a healthy human has been used to analyze the fields of physical stresses that may develop at the air-wall interface of the nasal mucosa. Simulations during quiet breathing revealed wall shear stresses as high as 0.3 Pa in the noselike model and 1.5 Pa in the anatomical model. These values are of the same order of those known to exist in uniform large arteries. The distribution of temperature near the nasal wall at peak inspiration is similar to that of wall shear stresses. The lowest temperatures occur in the vicinity of high stresses due to the narrow passageway in these locations. Time and spatial gradients of these stresses may have functional effects on nasal sensation of airflow and may play a role in the well-being of nasal breathing.     

Accuracy of mini peak flow meters in indicating changes in lung function in children with asthma.

OBJECTIVE--To assess whether mini flow meters used to measure peak expiratory flow can track changes in lung function and indicate clinically important changes. DESIGN--Comparison of measurements with a spirometer and different brands of mini flow meter; the meters were allocated to subjects haphazardly. SUBJECTS--12 boys with asthma aged 11 to 17 attending boarding school. MAIN OUTCOME MEASURES--Peak expiratory flow measured twice daily for three months with a spirometer and at least one of four brands of mini flow meter. RESULTS--The relation between changes in lung function measured with the spirometer and those measured with the mini flow meters was generally poor. In all, 26 episodes (range 1-3 in an individual child) of clinically important deterioration in lung function were detected from the records obtained with the spirometer. One mini flow meter detected six of 19 episodes, one detected six of 15, one detected six of 18, and one detected three of 21. CONCLUSIONS--Not only are the absolute values of peak expiratory flow obtained with mini flow meters inaccurate but the clinical message may also be incorrect. These findings do not imply that home monitoring of peak expiratory flow has no place in the management of childhood asthma but that the values obtained should be interpreted cautiously.     

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.     

Microvascular Endothelial Cells Migrate Upstream and Align Against the Shear Stress Field Created by Impinging Flow

At present, little is known about how endothelial cells respond to spatial variations in fluid shear stress such as those that occur locally during embryonic development, at heart valve leaflets, and at sites of aneurysm formation. We built an impinging flow device that exposes endothelial cells to gradients of shear stress. Using this device, we investigated the response of microvascular endothelial cells to shear-stress gradients that ranged from 0 to a peak shear stress of 9–210 dyn/cm². We observe that at high confluency, these cells migrate against the direction of fluid flow and concentrate in the region of maximum wall shear stress, whereas low-density microvascular endothelial cells that lack cell-cell contacts migrate in the flow direction. In addition, the cells align parallel to the flow at low wall shear stresses but orient perpendicularly to the flow direction above a critical threshold in local wall shear stress. Our observations suggest that endothelial cells are exquisitely sensitive to both magnitude and spatial gradients in wall shear stress. The impinging flow device provides a, to our knowledge, novel means to study endothelial cell migration and polarization in response to gradients in physical forces such as wall shear stress.     

Model of forced expiratory flows and airway geometry in infants.

Early measurements of autopsied lungs from infants, children, and adults suggested that the ratio of peripheral to central airway resistance was higher in infants than older children and adults. Recent measurements of forced expiration suggest that infants have high flows relative to lung volume. We employed a computational model of forced expiratory flow along with physiological and anatomic data to evaluate whether the infant lung is a uniformly scaled-down version of the adult lung. First, we uniformly scaled an existing computational model of adult forced expiration to estimate forced expiratory flows (FEF) and density dependence for an 18-mo-old infant. The values obtained for FEF and density dependence were significantly lower than those reported for healthy 18-mo-old infants. Next, we modified the model for the infant lung to reproduce standard indexes of expiratory flow [forced expiratory volume in 0.5 s (FEV(0.5)), FEFs after exhalation of 50 and 75% forced vital capacity, FEF between 25 and 75% expired volume] for this age group. The airway sizes obtained for the infant lung model that produced accurate physiological measurements were similar to anatomic data available for this age and larger than those in the scaled model. Our findings indicate that the airways in the infant lung model differ from those in the scaled model, i.e., middle and peripheral airway sizes are larger than result from uniform downscaling of the adult lung model. We show that the infant lung model can be made to reproduce individual flow-volume curves by adjusting lumen area generation by generation.     

Microvascular Endothelial Cells Migrate Upstream and Align Against the Shear Stress Field Created by Impinging Flow

At present, little is known about how endothelial cells respond to spatial variations in fluid shear stress such as those that occur locally during embryonic development, at heart valve leaflets, and at sites of aneurysm formation. We built an impinging flow device that exposes endothelial cells to gradients of shear stress. Using this device, we investigated the response of microvascular endothelial cells to shear-stress gradients that ranged from 0 to a peak shear stress of 9–210 dyn/cm². We observe that at high confluency, these cells migrate against the direction of fluid flow and concentrate in the region of maximum wall shear stress, whereas low-density microvascular endothelial cells that lack cell-cell contacts migrate in the flow direction. In addition, the cells align parallel to the flow at low wall shear stresses but orient perpendicularly to the flow direction above a critical threshold in local wall shear stress. Our observations suggest that endothelial cells are exquisitely sensitive to both magnitude and spatial gradients in wall shear stress. The impinging flow device provides a, to our knowledge, novel means to study endothelial cell migration and polarization in response to gradients in physical forces such as wall shear stress.     

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.     

Estimation of disruption of animal cells by laminar shear stress

Using mechanical cell properties measured by micromanipulation, and a model of cell distortion in laminar flow fields, a method has been developed for predicting disruption of animal cells by laminar shear stresses. Predictions of the model were compared with measured losses of cell number and viability of TB/C3 murine hybridomas sheared in a cone and plate viscometer at shear rates up to 3950 s(-1), and shear stresses up to 600 Nm(-2), achieved by enhancement of viscosity with dextran. In all cases, the experimental, results and predictions were within 30%. Such excellent agreement suggests it might be possible to use micromanipulation measurements of animal cell mechanical properties to predict cell damage in more complex flow fields, such as those in bioreactors. (c) 1992 John Wiley & Sons, Inc.