Abdominal aortic hemodynamics in young healthy adults at rest and during lower limb exercise: quantification using image-based computer modeling

Localization of atherosclerotic lesions in the abdominal aorta has been previously correlated to areas of adverse hemodynamic conditions, such as flow recirculation, low mean wall shear stress, and high temporal oscillations in shear. Along with its many systemic benefits, exercise is also proposed to have local benefits in the vasculature via the alteration of these regional flow patterns. In this work, subject-specific models of the human abdominal aorta were constructed from magnetic resonance angiograms of five young, healthy subjects, and computer simulations were performed under resting and exercise (50% increase in resting heart rate) pulsatile flow conditions. Velocity fields and spatial variations in mean wall shear stress (WSS) and oscillatory shear index (OSI) are presented. When averaged over all subjects, WSS increased from 4.8 +/- 0.6 to 31.6 +/- 5.7 dyn/cm2 and OSI decreased from 0.22 +/- 0.03 to 0.03 +/- 0.02 in the infrarenal aorta between rest and exercise. WSS significantly increased, whereas OSI decreased between rest and exercise at the supraceliac, infrarenal, and suprabifurcation levels, and significant differences in WSS were found between anterior and posterior sections. These results support the hypothesis that exercise provides localized benefits to the cardiovascular system through acute mechanical stimuli that trigger longer-term biological processes leading to protection against the development or progression of atherosclerosis.     

Inferior vena caval hemodynamics quantified in vivo at rest and during cycling exercise using magnetic resonance imaging

Compared with the abdominal aorta, the hemodynamic environment in the inferior vena cava (IVC) is not well described. With the use of cine phase-contrast magnetic resonance imaging (MRI) and a custom MRI-compatible cycle in an open magnet, we quantified mean blood flow rate, wall shear stress, and cross-sectional lumen area in 11 young normal subjects at the supraceliac and infrarenal levels of the aorta and IVC at rest and during dynamic cycling exercise. Similar to the aorta, the IVC experienced significant increases in blood flow and wall shear stress as a result of exercise, with greater increases in the infrarenal level compared with the supraceliac level. At the infrarenal level during resting conditions, the IVC experienced higher mean flow rate than the aorta (1.2 +/- 0.5 vs. 0.9 +/- 0.4 l/min, P < 0.01) and higher mean wall shear stress than the aorta (2.0 +/- 0.6 vs. 1.3 +/- 0.6 dyn/cm(2), P < 0.005). During exercise, wall shear stress remained higher in the IVC compared with the aorta, although not significantly. It was also observed that, whereas the aorta tapers inferiorly, the IVC tapers superiorly from the infrarenal to the supraceliac location. The hemodynamic and anatomic data of the IVC acquired in this study add to our understanding of the venous circulation and may be useful in a clinical setting.     

Numerical simulation of steady flow fields in a model of abdominal aorta with its peripheral branches

In the present study, a numerical calculation procedure based on a finite volume method was developed to simulate steady flow fields in a model of abdominal aorta with its peripheral branches. The study focused on the steady baseline flow fields and the wall shear stress (WSS) distribution as well as the localization of the reversed flow regions and results were compared to those obtained by other investigators. In the case of resting conditions, the existence of a region of reversed flow of about one to two diameters in size and next to the renal arteries and along the posterior wall as observed by other researchers was confirmed. However, under the exercise conditions this region could be wiped out. The flow reversal along the lateral walls proximal to the bifurcation persisted in both rest and exercise conditions. The WSS distribution and the wall shear stress gradient distribution were obtained. The lowest WSS occurred near the ostia of the renal arteries and the lateral walls of the iliac arteries. And the highest is always at the turn to the branch. The results were generally consistent with those obtained experimentally and numerically by other investigators. It was also shown that the steady flow might be used to depict the averaged behavior of pulsatile flow. The present computer code provides a platform for the future more realistic simulations.     

Hemodynamics simulation and identification of susceptible sites of atherosclerotic lesion formation in a model abdominal aorta

Employing the rabbit's abdominal aorta as a suitable atherosclerotic model, transient three-dimensional blood flow simulations and monocyte deposition patterns were used to evaluate the following hypotheses: (i) simulation of monocyte transport through a model of the rabbit abdominal aorta yields cell deposition patterns similar to those seen in vivo, and (ii) those deposition patterns are correlated with hemodynamic wall parameters related to atherosclerosis. The deposition pattern traces a helical shape down the aorta with local elevation in monocyte adhesion around vessel branches. The cell deposition pattern was altered by an exercise waveform with fewer cells attaching in the upper abdominal aorta but more attaching around the renal orifices. Monocyte deposition was correlated with the wall shear stress gradient and the wall shear stress angle gradient. The wall stress gradient, the wall shear stress angle gradient and the normalized monocyte deposition fraction were correlated with the distribution of monocytes along the abdominal aorta and monocyte deposition is correlated with the measured distribution of monocytes around the major abdominal branches in the cholesterol-fed rabbit. These results suggest that the transport and deposition pattern of monocytes to arterial endothelium plays a significant role in the localization of lesions.     

Does Lower Limb Exercise Worsen Renal Artery Hemodynamics in Patients with Abdominal Aortic Aneurysm?

Renal artery stenosis (RAS) and renal complications emerge in some patients after endovascular aneurysm repair (EVAR) to treat abdominal aorta aneurysm (AAA). The mechanisms for the causes of these problems are not clear. We hypothesized that for EVAR patients, lower limb exercise could negatively influence the physiology of the renal artery and the renal function, by decreasing the blood flow velocity and changing the hemodynamics in the renal arteries. To evaluate this hypothesis, pre- and post-operative models of the abdominal aorta were reconstructed based on CT images. The hemodynamic environment was numerically simulated under rest and lower limb exercise conditions. The results revealed that in the renal arteries, lower limb exercise decreased the wall shear stress (WSS), increased the oscillatory shear index (OSI) and increased the relative residence time (RRT). EVAR further enhanced these effects. Because these parameters are related to artery stenosis and atherosclerosis, this preliminary study concluded that lower limb exercise may increase the potential risk of inducing renal artery stenosis and renal complications for AAA patients. This finding could help elucidate the mechanism of renal artery stenosis and renal complications after EVAR and warn us to reconsider the management and nursing care of AAA patients.     

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.     

A numerical study of flow in curved tubes simulating coronary arteries

Numerical simulations of pulsatile flow in coronary arteries which take into account the curvature associated with the bending of arteries over the surface of the heart are presented for resting, excited and drug induced states. The study was motivated by reported observations of atherosclerotic plaque localization on the inner curvature of coronary arteries. The simulated flow field appears quasi-steady under resting conditions with wall shear stress always highest on the outside wall and only a single secondary flow vortex in the half tube. However, reversal of wall shear stress direction at the inside wall does occur under resting flow conditions and this is not a quasi-steady characteristic. The flow field is markedly unsteady under excited conditions with wall shear stress sometimes peaking on the inside wall and an increase in the magnitude of wall shear stress reversal on the inside wall. However, only a single secondary flow vortex in the half tube is observed. Implications of the simulations for the role of fluid mechanics in coronary artery atherosclerosis are also discussed.     

Blood flow patterns in the proximal human coronary arteries: relationship to atherosclerotic plaque occurrence

Atherosclerotic plaques in human coronary arteries are focal manifestations of systemic disease, and biomechanical factors have been hypothesized to contribute to plaque genesis and localization. We developed a computational fluid dynamics (CFD) model of the ascending aorta and proximal sections of the right and left coronary arteries of a normal human subject using computed tomography (CT) and magnetic resonance imaging (MRI) and determined the pulsatile flow field. Results demonstrate that flow patterns in the ascending aorta contribute to a pro-atherosclerotic flow environment, specifically through localization of low and oscillatory wall shear stress in the neighborhood of coronary orifices. Furthermore, these patterns differ in their spatial distribution between right and left coronary arteries. Entrance effects of aortic flow diminish within two vessel diameters. We examined relationships between spatial distributions of wall shear stress and reports of plaque occurrence in the literature. Results indicate low wall shear stress is co-located with increased incidence of lesions, and higher wall shear stresses are associated with lesion-resistant areas. This investigation does not consider plaque progression or advanced lesions, inasmuch as the CFD model was developed from a normal individual and the clinical data used for comparisons were obtained from autopsy specimens of subjects who died from non-cardiovascular causes. The data reported are consistent with the hypothesis that low wall shear stress is associated with the localization of atherosclerotic lesions, and the results demonstrate the importance of aortic flow on flow patterns in the proximal segments of the coronary arteries.     

Influence of the renal artery ostium flow diverter on hemodynamics and atherogenesis

The structure and function of the renal artery ostium flow diverter on the caudal side of the renal branch point were previously reported; in this study, we further evaluate the diverter?s possible functions. The protrusion of this structure into the abdominal aorta suggests that the diverter may preferentially direct blood flow to the renal arteries, and that it may also influence flow patterns and recirculation known to be involved in atherogenesis. Three-dimensional computational fluid dynamics (CFD) simulations of steady and pulsatile blood flow are performed to investigate the influence of diverter size and position, and vascular geometry, on the flow patterns and fluid mechanical forces in the neighborhood of the diverter. CFD results show that the flow diverter does affect the blood distribution; depending on the diverter?s position, the flow to the renal arteries may be increased or reduced. Calculated results also demonstrate the diverter?s effect on the wall shear stress (WSS) distribution, and suggest that the diverter contributes to an atherogenic environment in the abdominal aorta, while being atheroprotective in the renal arteries themselves. These results support previous clinical findings, and suggest directions for further clinical study. The results of this work have direct implications in understanding the physiological significance of the diverter, and its potential role in the pathophysiological development of atherosclerosis.     

Spatial comparison between wall shear stress measures and porcine arterial endothelial permeability

A better understanding of how hemodynamic factors affect the integrity and function of the vascular endothelium is necessary to appreciate more fully how atherosclerosis is initiated and promoted. A novel technique is presented to assess the relation between fluid dynamic variables and the permeability of the endothelium to macromolecules. Fully anesthetized, domestic swine were intravenously injected with the albumin marker Evans blue dye, which was allowed to circulate for 90 min. After the animals were euthanized, silicone casts were made of the abdominal aorta and its iliac branches. Pulsatile flow calculations were subsequently made in computational regions derived from the casts. The distribution of the calculated time-dependent wall shear stress in the external iliac branches was directly compared on a point-by-point basis with the spatially varying in vivo uptake of Evans blue dye in the same arteries. The results indicate that in vivo endothelial permeability to albumin decreases with increasing time-average shear stress over the normal range. Additionally, endothelial permeability increases slightly with oscillatory shear index.     

Effect of alterations in femoral artery flow on abdominal vessel hemodynamics in swine

In support of an in vivo investigation in swine of the influence of changes in fluid dynamic wall shear on arterial macromolecular permeability, a procedure has been developed to alter the flows in the porcine posterior arterial vasculature by opening and closing a reversible arteriovenous shunt placed on one of the femoral arteries. Laparoscopic techniques were used to place appropriately modified Transonic Systems ultrasonic flow probes on both external and circumflex iliac arteries, and on the terminal aorta. Flow measurements were made prior to shunt placement, and with the shunt open and closed, to measure the influence of altered external iliac artery flow on the distribution to the infrarenal abdominal vessels. Similar experiments were carried out to relate the flow rates in the external iliac arteries to those in the femoral arteries, which are more accessible. Based on the relationships among the measured flow rates, rules have been developed to estimate the major infrarenal flows in the pig, at baseline and with the shunt opened and closed, from only the flow rates measured at the two femoral arteries.     

Comparison of physiological and post-endovascular aneurysm repair infrarenal blood flow

Endovascular aneurysm repair (EVAR) of abdominal aortic aneurysms results in redirection of blood through the deployed endograft (EG). Even though EVAR is clinically effective, the absolute flow restoration is not warranted. Our purpose was to compare the physiological with the post-EVAR infrarenal flow conditions. We developed patient-specific models based on computed tomography data of five healthy volunteers and ten patients treated with the Endurant^® stent-graft system. Wall shear stress (WSS), helicity, pressure and velocity fields were calculated using computational fluid dynamics. The results showed a decrease of peak WSS on the part of the EG that resides in the iliac arteries, compared to the physiological value (p = 0.01). At the abdominal part, the average helicity seems to increase after EVAR, while at the iliac arteries part, the intensity of helical flow seems physiological. Pressure drop and peak velocity in the iliac arteries part are lower than the physiological values (p = 0.04). The comparison revealed that most hemodynamic properties converge to normal levels at the abdominal part whereas statistically significant variations were observed in the iliac arteries part. The delineation of the differences between physiological and postoperative flow data could pave the way for the improvement of EG designs.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

Steady flow visualization in a rigid canine aortic cast

Steady flow studies were conducted in a transparent canine aortic cast. The cast segment stretched from the aortic valve to beyond the renal arteries and included all major branches. Flow was visualized by analysis of dye streaklines. Flow rates for basal and exercising cardiovascular states were simulated. The Reynolds numbers in the ascending aorta for basal and exercising conditions were 900 and 1587 respectively. Aortic core flow was laminar in basal simulations. Disturbed flow commenced in the upper descending aorta with exercising flow rates. Separation zones existed along the inner curvature of the aortic arch and the proximal walls of the brachiocephalic, left subclavian, and coeliac arteries. Such zones may exist over a portion of the cardiac cycle. If either renal artery was occluded, then a vortex formed. This vortex is associated with high shear regions which correlate well with sites where sudanophilic lesions have been reported in cholesterol-fed nephrectomized rabbits.     

Experimental determination of velocity profiles and wall shear rate along the rabbit aortoiliac bifurcation: relationship to vessel wall low-density lipoprotein (LDL) metabolism.

We have determined the velocity profiles and wall shear rates along the New Zealand White (NZW) rabbit aortoiliac bifurcation. A pulsatile perfusion apparatus was used to impose physiologic pressure and flow waveforms on nine freshly excised NZW bifurcation segments. Pulsed Doppler velocimetry (PDV) was utilized to construct velocity profiles at five measurement sites: within the infrarenal aorta; immediately distal to the apex of the bifurcation; and, more distally along the iliac arteries. Wall shear rate was derived from a numerical differentiation of the experimental velocity profiles. The results of this study indicate that the average shear rate was lower along the lateral (approximately 40 s-1) vs medial (approximately 240 s-1) wall of the proximal iliac branch. The degree of flow reversal along the proximal lateral walls (20 +/- 2%) exceeded that along the proximal flow divider wall (1 +/- 1%). Flow at the distal iliac measurement sites and within the infrarenal aorta was approximately symmetric. These findings complement our companion in vivo study [Berceli et al., Arteriosclerosis 10, 688-694 (1990)] wherein we determined the rates of low-density lipoprotein (LDL) incorporation and catabolism along this symmetrically bifurcating conduit. Taken together, these studies provide original information regarding the effects of hemodynamics on one presumed atherogenic risk factor, namely, LDL metabolism.     

Laser Doppler anemometer measurements of pulsatile flow in a model carotid bifurcation

Hemodynamics at the human carotid bifurcation is important to the understanding of atherosclerotic plaque initiation and progression as well as to the diagnosis of clinically important disease. Laser Doppler anemometry was performed in a large scale model of an average human carotid. Pulsatile waveforms and physiologic flow divisions were incorporated. Disturbance levels and shear stresses were computed from ensemble averages of the velocity waveform measurements. Flow in the common carotid was laminar and symmetric. Flow patterns in the sinus, however, were complex and varied considerably during the cycle. Strong helical patterns and outer wall flow separation waxed and waned during each systole. The changing flow patterns resulted in an oscillatory shear stress at the outer wall ranging from -13 to 9 dyn cm-2 during systole with a time-averaged mean of only -0.5 dyn cm-2. This contrasts markedly with an inner wall shear stress range of 17-50, (mean 26) dyn cm-2. The region of transient separation was confined to the carotid sinus outer wall with no reverse velocities detected in the distal internal carotid. Notable disturbance velocities were also time-dependent, occurring only during the deceleration phase of systole and the beginning of diastole. The present pulsatile flow studies have aided in identifying hemodynamic conditions which correlate with early intimal thickening and predict the physiologic level of flow disturbances in the bulb of undiseased internal carotid arteries.     

Dynamic curvature strongly affects wall shear rates in a coronary artery bifurcation model

This study was motivated by the need for a better understanding of coronary artery blood flow patterns and their possible role in atherosclerosis formation. Of particular interest in this study was the effects of the dynamic deformation due to myocardial contraction on wall shear rate patterns in the coronary arteries. A better understanding of these effects on wall shear rate in a bifurcation geometry and an evaluation of the importance of these effects was desired. A three-dimensional computer model of a bifurcation lying on the surface of a sphere with time-varying radius of curvature was employed to simulate the motion and deformation of the arteries. The results indicated low mean shear rates along the myocardial wall and very high shear rate variations (over 100% of the static mean shear rate) along the outer wall. The results obtained using a quasi-static analysis were found to underestimate the dynamic wall shear rate variation along the myocardial and outer walls. It was concluded that dynamic geometry effects are important in determining sites of low mean and oscillating wall shear that have been associated with atherogenesis in curved, bifurcating arteries.     

Allometric scaling of wall shear stress from mice to humans: quantification using cine phase-contrast MRI and computational fluid dynamics

Allometric scaling laws relate structure or function between species of vastly different sizes. They have rarely been derived for hemodynamic parameters known to affect the cardiovascular system, e.g., wall shear stress (WSS). This work describes noninvasive methods to quantify and determine a scaling law for WSS. Geometry and blood flow velocities in the infrarenal aorta of mice and rats under isoflurane anesthesia were quantified using two-dimensional magnetic resonance angiography and phase-contrast magnetic resonance imaging at 4.7 tesla. Three-dimensional models constructed from anatomic data were discretized and used for computational fluid dynamic simulations using phase-contrast velocity imaging data as inlet boundary conditions. WSS was calculated along the infrarenal aorta and compared between species to formulate an allometric equation for WSS. Mean WSS along the infrarenal aorta was significantly greater in mice and rats compared with humans (87.6, 70.5, and 4.8 dyn/cm(2), P < 0.01), and a scaling exponent of -0.38 (R(2) = 0.92) was determined. Manipulation of the murine genome has made small animal models standard surrogates for better understanding the healthy and diseased human cardiovascular system. It has therefore become increasingly important to understand how results scale from mouse to human. This noninvasive methodology provides the opportunity to serially quantify changes in WSS during disease progression and/or therapeutic intervention.     

Image-based modeling of hemodynamics in coronary artery aneurysms caused by Kawasaki disease

Kawasaki Disease (KD) is the leading cause of acquired pediatric heart disease. A subset of KD patients develops aneurysms in the coronary arteries, leading to increased risk of thrombosis and myocardial infarction. Currently, there are limited clinical data to guide the management of these patients, and the hemodynamic effects of these aneurysms are unknown. We applied patient-specific modeling to systematically quantify hemodynamics and wall shear stress in coronary arteries with aneurysms caused by KD. We modeled the hemodynamics in the aneurysms using anatomic data obtained by multi-detector computed tomography (CT) in a 10-year-old male subject who suffered KD at age 3?years. The altered hemodynamics were compared to that of a reconstructed normal coronary anatomy using our subject as the model. Computer simulations using a robust finite element framework were used to quantify time-varying shear stresses and particle trajectories in the coronary arteries. We accounted for the cardiac contractility and the microcirculation using physiologic downstream boundary conditions. The presence of aneurysms in the proximal coronary artery leads to flow recirculation, reduced wall shear stress within the aneurysm, and high wall shear stress gradients at the neck of the aneurysm. The wall shear stress in the KD subject (2.95-3.81 dynes/sq cm) was an order of magnitude lower than the normal control model (17.10-27.15 dynes/sq cm). Particle residence times were significantly higher, taking 5 cardiac cycles to fully clear from the aneurysmal regions in the KD subject compared to only 1.3 cardiac cycles from the corresponding regions of the normal model. In this novel quantitative study of hemodynamics in coronary aneurysms caused by KD, we documented markedly abnormal flow patterns that are associated with increased risk of thrombosis. This methodology has the potential to provide further insights into the effects of aneurysms in KD and to help risk stratify patients for appropriate medical and surgical interventions.