Haemodynamics of angioaccess venous anastomoses

A large percentage of arteriovenous haemodialysis angioaccess loop grafts (AVLG) fail within the first year after surgery, the occlusive lesions being found predominantly at the venous anastomosis site. This paper presents a detailed flow dynamic study of the AVLG system using three elastic, transparent bench-top flow models, which were based on the geometry of silicone rubber casts obtained at different times from a chronic animal model. Each model thus represented a different stage of the lesion development. Flow visualization and laser Doppler anemometer surveys of the flow field confirmed that the hydrodynamic factors favour lesion development near the stagnation point opposite the anastomotic toe, where the momentum of the impinging jet stream, combined with the oscillating wall shear stress generated in the vicinity of the stagnation point, acts in both directions. The accumulation of tracer particles in the region of flow separation is believed to be a combined contribution from the hydraulic forces and the inward motion of the vessel wall. As these hydrodynamic factors are enhanced upon further development of the occlusive lesion, a vicious cycle may be formed.     

Flow profiles and wall shear stress distribution at a hemodialysis venous anastomosis: preliminary study

The phasic velocity field in the vicinity of the venous anastomosis in a hemodialysis angioaccess arteriovenous fistula loop graft (AVLG) is investigated employing a laser Doppler anemometer (LDA) system. Detailed LDA velocity profiles are obtained by sectional survey performed in a transparent, elastic flow model which was fabricated to represent the geometry of the AVLG system under physiological pressure and flow waveforms. The geometry of the flow model was based on a silicone rubber cast obtained from an experimental dog model. In the present study, detailed distribution of velocity profiles is obtained. The distribution of wall shear stress in the model is computed from the slope of the local velocity profiles near the wall. The relationship between the results obtained by flow visualization and the LDA measurement is discussed.     

Computational model of blood flow in the aorto-coronary bypass graft

Background

Coronary artery bypass grafting surgery is an effective treatment modality for patients with severe coronary artery disease. The conduits used during the surgery include both the arterial and venous conduits. Long- term graft patency rate for the internal mammary arterial graft is superior, but the same is not true for the saphenous vein grafts. At 10 years, more than 50% of the vein grafts would have occluded and many of them are diseased. Why do the saphenous vein grafts fail the test of time? Many causes have been proposed for saphenous graft failure. Some are non-modifiable and the rest are modifiable. Non-modifiable causes include different histological structure of the vein compared to artery, size disparity between coronary artery and saphenous vein. However, researches are more interested in the modifiable causes, such as graft flow dynamics and wall shear stress distribution at the anastomotic sites. Formation of intimal hyperplasia at the anastomotic junction has been implicated as the root cause of long- term graft failure.Many researchers have analyzed the complex flow patterns in the distal sapheno-coronary anastomotic region, using various simulated model in an attempt to explain the site of preferential intimal hyperplasia based on the flow disturbances and differential wall stress distribution. In this paper, the geometrical bypass models (aorto-left coronary bypass graft model and aorto-right coronary bypass graft model) are based on real-life situations. In our models, the dimensions of the aorta, saphenous vein and the coronary artery simulate the actual dimensions at surgery. Both the proximal and distal anastomoses are considered at the same time, and we also take into the consideration the cross-sectional shape change of the venous conduit from circular to elliptical. Contrary to previous works, we have carried out computational fluid dynamics (CFD) study in the entire aorta-graft-perfused artery domain. The results reported here focus on (i) the complex flow patterns both at the proximal and distal anastomotic sites, and (ii) the wall shear stress distribution, which is an important factor that contributes to graft patency.

Methods

The three-dimensional coronary bypass models of the aorto-right coronary bypass and the aorto-left coronary bypass systems are constructed using computational fluid-dynamics software (Fluent 6.0.1). To have a better understanding of the flow dynamics at specific time instants of the cardiac cycle, quasi-steady flow simulations are performed, using a finite-volume approach. The data input to the models are the physiological measurements of flow-rates at (i) the aortic entrance, (ii) the ascending aorta, (iii) the left coronary artery, and (iv) the right coronary artery.

Results

The flow field and the wall shear stress are calculated throughout the cycle, but reported in this paper at two different instants of the cardiac cycle, one at the onset of ejection and the other during mid-diastole for both the right and left aorto-coronary bypass graft models. Plots of velocity-vector and the wall shear stress distributions are displayed in the aorto-graft-coronary arterial flow-field domain. We have shown (i) how the blocked coronary artery is being perfused in systole and diastole, (ii) the flow patterns at the two anastomotic junctions, proximal and distal anastomotic sites, and (iii) the shear stress distributions and their associations with arterial disease.

Conclusion

The computed results have revealed that (i) maximum perfusion of the occluded artery occurs during mid-diastole, and (ii) the maximum wall shear-stress variation is observed around the distal anastomotic region. These results can enable the clinicians to have a better understanding of vein graft disease, and hopefully we can offer a solution to alleviate or delay the occurrence of vein graft disease.

     

The effect of angle and flow rate upon hemodynamics in distal vascular graft anastomoses: an in vitro model study.

A steady flow, in vitro model of distal arterial bypass graft junctions was used to examine the effects of junction angle and flow rate on the local velocity field. Three test sections were fabricated from Plexiglas tubing having anastomotic junction angles of either 30, 45, or 60 deg. Flow visualization revealed velocity profiles skewed toward the outer wall with a flow split around a clear stagnation point along the outer wall. Laser Doppler anemometry [LDA] measurements confirmed a distinct stagnation point at the outer wall and both reverse and forward shear were detected immediately upstream and downstream, respectively, of this site. Axial velocities and shear rates along the outer wall were higher than along the inner wall and occurred in the junction angle order: 45, 60, and 30 deg. This study clearly identified changes in wall shear which varied with the anastomotic angle and flow rate.     

犬不同血管搭桥方法及搭桥血管内流场的计算机数值模拟

目的血管搭桥术后的内膜增生往往导致手术失败,而内膜增生与搭桥血管内的流场密切相关,为改善搭桥血管中的流场结构,作者设计了偏心搭桥手术方法,利用计算机数值模拟技术,探索偏心搭桥和传统搭桥血管中流场的变化,为血管搭桥方法提供优化设计方案。方法16只犬随机分为偏心搭桥组和传统搭桥组进行血管搭桥,测定搭桥前后血管几何数据,搭桥后近心端及远心端吻合口血流量和血压。按测定的血管几何数据,FLUENT 6.2模拟搭桥血管内的流场。结果偏心搭桥近心端和远心端吻合口不在同一平面。传统搭桥中,主体动脉远心端吻合口对应面处存在一个较低壁面剪切应力（WSS）区域及流体停滞点,离脚跟较近的一部分流体会形成涡漩,血流进入主体动脉后,还会表现出迪恩涡二次流;偏心搭桥中,主体动脉吻合口对应面上的低WSS区域和流体停滞点消失,血流接触到吻合口底面后,以切向旋转的方式改变其流动方向,不会形成涡漩,且当血流进入主体动脉后,立即发生螺旋流态且能持续很长一段。结论偏心搭桥能够产生血液旋动流,显著增加远心端血流量、提高WSS。     

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.     

The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis

We present an experimental and computational investigation of time-varying flow in an idealized fully occluded 45 degrees distal end-to-side anastomosis. Two geometric configurations are assessed, one where the centerlines of host and bypass vessels lie within a plane, and one where the bypass vessel is deformed out of the plane of symmetry, respectively, termed planar and non-planar. Flow experiments were conducted by magnetic resonance imaging in rigid wall models and computations were performed using a high order spectral/hp algorithm. Results indicate a significant change in the spatial distribution of wall shear stress and a reduction of the time-averaged peak wall shear stress magnitude by 10% in the non-planar model as compared to the planar configuration. In the planar geometry the stagnation point follows a straight-line path along the host artery bed with a path length of 0.8 diameters. By contrast in the non-planar case the stagnation point oscillates about a center that is located off the symmetry plane intersection with the host artery bed wall, and follows a parabolic path with a 0.7 diameter longitudinal and 0.5 diameter transverse excursion. A definition of the oscillatory shear index (OSI) is introduced that varies between 0 and 0.5 and that accounts for a continuous range of wall shear stress vector angles. In both models, regions of elevated oscillatory shear were spatially associated with regions of separated or oscillating stagnation point flow. The mean oscillatory shear magnitude (considering sites where OSI>0.1) in the non-planar geometry was reduced by 22% as compared to the planar configuration. These changes in the dynamic behavior of the stagnation point and the oscillatory shear distribution introduced by out-of-plane graft curvature may influence the localization of vessel wall sites exposed to physiologically unfavorable flow conditions.     

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.     

The mechanism of mucus clearance in cough

An instability resembling an avalanche is proposed as the mechanism by which mucus is expelled from the respiratory tract during cough. The cough event was simulated in a model airway. In these experiments, air was forced through a channel whose walls were lined with a non-Newtonian material rheologically similar to tracheal mucus. Frames from high-speed cine photographs showed an unstable event which began as an undulation of the free surface and progressed to a catastrophic clearance of the channel. Measurements of the longitudinal pressure gradient support the hypothesis that the clearance event is initiated when the total stress applied to the mucus analog exceeds its finite yield stress. A continuum model predicts that yielding occurs within the bottom layers of the mucus analog. Calculations based upon estimates of tracheal geometry and air flow show that the clearance event studied here would be expected to occur during a cough but not during normal breathing. Experiments also show that a lubricant introduced between the channel walls and the mucus blanket can reduce the air flow rate required to precipitate the clearance.     

Hemodynamic factors at the distal end-to-side anastomosis of a bypass graft with different POS:DOS flow ratios

A pulsatile flow in vitro model of the distal end-to-side anastomosis of an arterial bypass graft was used to examine the effects that different flow ratios between the proximal outlet segment (POS) and the distal outlet segment (DOS) have on the flow patterns and the distributions of hemodynamic factors in the anastomosis. Amberlite particles were tracked by flow visualization to determine overall flow patterns and velocity measurements were made with Laser Doppler anemometry (LDA) to obtain detailed hemodynamic factors along the artery floor and the graft hood regions. These factors included wall shear stress (WSS), spatial wall shear stress gradient (WSSG), and oscillatory index (OSI). Statistical analysis was used to compare these hemodynamic factors between cases having different POS:DOS flow ratios (Case 1-0:100, Case 2-25:75, Case 3-50:50). The results showed that changes in POS:DOS flow ratios had a great influence on the flow patterns in the anastomosis. With an increase in proximal outlet flow, the range of location of the stagnation point along the artery floor decreased, while the extent of flow separation along the graft hood increased. The statistical results showed that there were significant differences (p<0.05) for the mean WSS between cases along the graft hood, but no significant differences were detected along the artery floor. There were no significant differences for the spatial WSSG along both the artery floor and the graft hood. However, there were significant differences (p<0.05) in the mean OSI between Cases 1 and 2 and between Cases 1 and 3 both along the artery floor and along the graft hood. Comparing these mechanical factors with histological findings of intimal hyperplasia formation obtained by previous canine studies, the results of the statistical analysis suggest that regions exposed to a combination of low mean WSS and high OSI may be most prone to the formation of intimal hyperplasia.     

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.     

Numerical study of hemodynamics and wall mechanics in distal end-to-side anastomoses of bypass grafts.

The development and progress of distal anastomotic intimal hyperplasia seems to be promoted by altered flow conditions and intramural stress distributions at the region of the artery-graft junction of vascular bypass configurations. From clinical observations, it is known that intimal hyperplasia preferentially occurs at outflow anastomoses of prosthetic bypass grafts. In order to gain a deeper insight into post-operative disease processes, and subsequently, to contribute to the development of improved vascular reconstructions with respect to long term patency rates, detailed studies are required. In context with in vivo experiments, this study was designed to analyze the flow dynamics and wall mechanics in anatomically correct bypass configurations related to two different surgical techniques and resulting geometries (conventional geometry and Miller-cuff). The influence of geometric conditions and of different compliance of synthetic graft, the host artery and the interposed venous cuff on the hemodynamic behavior and on the wall stresses are investigated. The flow studies apply the time-dependent, three-dimensional Navier-Stokes equations describing the motion of an incompressible Newtonian fluid. The vessel walls are described by a geometrically non-linear shell structure. In an iterative coupling procedure, the two problems are solved by means of the finite element method. The numerical results demonstrate non-physiological flow patterns in the anastomotic region. Strongly skewed axial velocity profiles and high secondary velocities occur downstream the artery-graft junction. On the artery floor opposite the junction, flow separation and zones of recirculation are found. The wall mechanical studies show that increased compliance mismatch leads to increased intramural stresses, and thus, may have a proliferative influence on suture line hyperplasia, as it is observed in the in vivo study.     

Combined MR imaging and numerical simulation of flow in realistic arterial bypass graft models

We report methods for (a) transforming a three-dimensional geometry acquired by magnetic resonance angiography (MRA) in vivo, or by imaging a model cast, into a computational surface representation, (b) use of this to construct a three dimensional numerical grid for computational fluid dynamic (CFD) studies, and (c) use of the surface representation to produce a stereo-lithographic replica of the real detailed geometry, at a scale convenient for detailed magnetic resonance imaging (MRI) flow studies. This is applied to assess the local flow field in realistic geometry arterial bypass grafts. Results from a parallel numerical simulation and MRI measurement of flow in an aorto-coronary bypass graft with various inlet flow conditions demonstrate the strong influence of the graft inlet waveform on the perianastomotic flow field. A sinusoidal and a multi harmonic coronary flow waveform both with a mean Reynolds number (Re) of 100 and a Womersley parameter of 2.7 were applied at the graft inlet. A weak axial flow separation region just distal to the toe was found in sinusoidal flow near end deceleration (Re = 25). At the same location and approximately the same point in the cycle (Re = 30) but in coronary flow, the axial flow separation was stronger and more spatially pronounced. No axial flow separation occurred in steady flow for Re = 100. Numerical predictions indicate a region in the vicinity of the suture line (where there is a local narrowing of the graft) with a wall shear magnitude in excess of five times that associated with fully developed flow at the graft inlet.     

Numerical 3D-simulation of pulsatile wall shear stress in an arterial T-bifurcation model

The structure of pulsatile blood flow and wall shear stress in a 90° T-bifurcation model is analysed numerically. The nonlinear Navier-Stokes equations for time-dependent incompressible Newtonian fluid flow are approximated using a newly developed pressure correction, finite element method. The wall shear stress is calculated from the finite element velocity field. The investigation shows viscous flow phenomena such as flow separation and stagnation and the distribution of high and low wall shear stress during the pulse cycle. Furthermore, the effect of a sharp corner the bifurcation edge on the wall shear stress is analysed. Detailed local flow investigation is required to examine fluid dynamic contribution to the development of arterial diseases such as atherosclerosis and thrombosis.     

SPY Imaging Assessment Correlates With Transesophageal Echocardiogram Assessment of Ventricular Function During Off-Pump Coronary Artery Bypass Grafting

Intraoperative assessment of graft anastomoses is commonly performed after off-pump coronary artery bypass grafting (OPCAB). The SPY imaging system allows intraoperative graft assessment. We document correlation between intraoperative SPY images and wall motion abnormality by transesophageal echocardiogram (TEE) during OPCAB. A 79-year-old female underwent OPCAB. Intraoperative graft patency assessment was performed with the SPY and left ventricular wall motion was assessed by TEE. SPY imaging demonstrated poor flow trough the distal vein graft anastomosis to the posterior descending artery, which correlated with a new posterior wall motion hypokinesis. After graft revision, SPY imaging demonstrated good distal flow and the TEE demonstrated normalization of the left ventricular posterior wall motion. SPY technology allows the surgeon to accurately assess graft patency intraoperatively and allows immediate correction of a technical problem.     

Multiphysics simulation of blood flow and LDL transport in a porohyperelastic arterial wall model

Atherosclerosis localizes at a bend andor bifurcation of an artery, and low density lipoproteins (LDL) accumulate in the intima. Hemodynamic factors are known to affect this localization and LDL accumulation, but the details of the process remain unknown. It is thought that the LDL concentration will be affected by the filtration flow, and that the velocity of this flow will be affected by deformation of the arterial wall. Thus, a coupled model of a blood flow and a deformable arterial wall with filtration flow would be invaluable for simulation of the flow field and concentration field in sequence. However, this type of highly coupled interaction analysis has not yet been attempted. Therefore, we performed a coupled analysis of an artery with multiple bends in sequence. First, based on the theory of porous media, we modeled a deformable arterial wall using a porohyperelastic model (PHEM) that was able to express both the filtration flow and the viscoelastic behavior of the living tissue, and simulated a blood flow field in the arterial lumen, a filtration flow field and a displacement field in the arterial wall using a fluid-structure interaction (FSI) program code by the finite element method (FEM). Next, based on the obtained results, we further simulated LDL transport using a mass transfer analysis code by the FEM. We analyzed the PHEM in comparison with a rigid model. For the blood flow, stagnation was observed downward of the bends. The direction of the filtration flow was only from the lumen to the wall for the rigid model, while filtration flows from both the wall to the lumen and the lumen to the wall were observed for the PHEM. The LDL concentration was high at the lumenwall interface for both the PHEM and rigid model, and reached its maximum value at the stagnation area. For the PHEM, the maximum LDL concentration in the wall in the radial direction was observed at the position of 3% wall thickness from the lumenwall interface, while for the rigid model, it was observed just at the lumenwall interface. In addition, the peak LDL accumulation area of the PHEM moved about according to the pulsatile flow. These results demonstrate that the blood flow, arterial wall deformation, and filtration flow all affect the LDL concentration, and that LDL accumulation is due to stagnation and the presence of filtration flow. Thus, FSI analysis is indispensable.     

Tissue engineering skin flaps: which vascular carrier, arteriovenous shunt loop or arteriovenous bundle, has more potential for angiogenesis and tissue generation?

This study was designed to clarify which vascular carrier, the arteriovenous shunt loop or the arteriovenous bundle, has more potential as a vascular carrier for an artificial skin flap in rats. An arteriovenous shunt loop was constructed between the femoral artery and vein using an interpositional artery (group I) or vein (group II) graft. For arteriovenous bundle groups, the femoral artery and vein were used and subdivided into two groups: distal ligation type (group III) and flow-through type (group IV). The vascular pedicle was wrapped with an artificial dermis and implanted beneath the inguinal skin for 4 weeks. For the control group, a folded sheet of artificial dermis without any vascular carrier was embedded. In experiment 1, the volumes of generated tissue within the artificial dermis were measured in the experimental and control groups (n = 5 in each group). In experiment 2, the origin of new blood vessels sprouting from the arteriovenous shunt loop and arteriovenous bundle were evaluated histologically. The volume of generated tissue in the shunt groups was significantly greater than that in the bundle groups (p < 0.01). However, the bundle groups also showed a great potential for producing new tissue. Serial histological studies showed that new capillaries were derived not only from the vasa vasorum of the femoral vessels but directly from the femoral vein in both the shunt and the bundle groups. This "sprouting" was extensively exhibited in the group III. Although the arteriovenous shunt loop showed a greater potential for producing new tissue and capillaries, the distal ligation type of bundle was thought to be an effective and practical vascular carrier for producing a tissue-engineered skin flap.     

Flow pattern and shear stress distribution of distal end-to-side anastomoses. A comparison of the instantaneous velocity fields obtained by particle image velocimetry

OBJECTIVE: To describe the local hemodynamics and pressure losses of crural bypass anastomoses using instantaneous velocity fields acquired by particle image velocimetry (PIV). METHODS: Silastic models of a Taylor patch, a Miller cuff and a femoro-crural patch prosthesis (FCPP) were attached to a circuit driven by a Berlin Heart, providing a pulsatile flow with an amplitude of 450 to 25 ml/min (mean 200 ml/min). An outflow resistance of 0.5 mmHg/ml/min (peripheral resistance units, PRU) was modeled using small silastic tubes providing a phase shift of -12 degrees between flow and pressure curves. The working fluid consisted of a glycerine/water mixture with a viscosity of 4 mPas. Hollow glass spheres with a mean size of 9-13 microm were used as tracer particles. Instantaneous velocity fields were obtained by means of PIV and shear rates as well as shear stresses were calculated. Triggered by the flowmeter signal, 10 measurements at 100 ms intervals per cardiac cycle were obtained. The pressures were measured on the inflow and at both distal outflows. The resulting mean pressure losses due to flow separation and distal fluid acceleration were calculated. RESULTS: Inside the Taylor patch anastomosis a large flow separation at the hood containing a clockwise rotating vortex was found. Additionally a smaller flow separation at the heel and a flow stagnation zone on the floor of the recipient artery were observed. Conversely, inside the Miller cuff a counterclockwise rotating vortex was seen inside a large heel flow separation. The FCPP also showed typical separation areas at the hood and heel of the anastomosis, although these were smaller compared to the other anastomoses. Inside the FCPP anastomosis no vortex creation was observed throughout the cardiac cycle. The mainstream velocities at the inlet levels were comparable for the three anastomoses. A significant fluid acceleration was present at the antegrade as well as the retrograde outlets of the Taylor and Miller cuff, while the fluid acceleration at the antegrade outflow of the FCPP was small, which was attributed to the end-to-end configuration of the antegrade FCPP leg. The calculated normalized antegrade and retrograde pressure losses for the Taylor form were 0.90 and 0.88, for the Miller cuff 0.89 and 0.86 and for the FCPP 0.94 and 0.86, respectively. The shear stresses inside the flow separations of the three anastomoses were significantly lower than normal wall shear stresses. High shear stress levels were found inside the transition zones between flow separation and high velocity mainstream. CONCLUSIONS: The flow pattern inside cuffed or funnel shaped anastomoses consists of large flow separation zones, which are thought to be associated with intimal hyperplasia development. In addition, fluid accelerations at the distal outlets result in pressure losses, which may contribute to impaired crural perfusion.     

Physiological flow analysis in significant human coronary artery stenoses

To evaluate the local hemodynamics in flow limiting coronary lesions, computational hemodynamics was applied to a group of patients previously reported by Wilson et al. (1988) with representative pre-angioplasty stenosis geometry (minimal lesion size d(m)=0.95 mm; 68% mean diameter stenosis) and with measured values of coronary flow reserve (CFR) in the abnormal range (2.3+/-0.1). The computations were at mean flow rates (Q) of 50, 75 and 100 ml/min (the limit of our converged calculations). Computed mean pressure drops Deltap were approximately 9 mmHg for basal flow (50 ml/min), approximately 27 mmHg for elevated flow (100 ml/min) and increased to an extrapolated value of approximately 34 mmHg for hyperemic flow (115 ml/min), which led to a distal mean coronary pressure p(rh) of approximately 55 mmHg, a level known to cause ischemia in the subendocardium (Brown et al., 1984), and consistent with the occurrence of angina in the patients. Relatively high levels of wall shear stress were computed in the narrow throat region and ranged from about 600 to 1500 dyn/cm(2), with periodic (phase shifted) peak systolic values of about 3500 dyn/cm(2). In the distal vessel, the interaction between the separated shear layer wave, convected downstream by the core flow, and the wall shear layer flow, led to the formation of vortical flow cells along the distal vessel wall during the systolic phase where Reynolds numbers Re(e)(t) were higher. During the phasic vortical mode observed at both basal and elevated mean flow rates, wide variations in distal wall shear stress occurred, distal transmural pressures were depressed below throat levels, and pressure recovery was larger farther along the distal vessel. Along the constriction (convergent) and throat segments of the lesion the pulsatile flow field was principally quasi-steady before flow separation occurred. The flow regimes were complex in the narrow mean flow Reynolds number range Re(e)=100-230 and a frequency parameter of alphae=2.25. The shear layer flow disturbances diminished in strength due to viscous damping along the distal vessel at these relatively low values of Re(e), typical of flow through diseased epicardial coronary vessels. The distal hyperemic flow field was likely to be in an early stage of turbulent flow development during the peak systolic phase.