The effect of proximal artery flow on the hemodynamics at the distal anastomosis of a vascular bypass graft: computational study

The formation of distal anastomotic intimal hyperplasia (IH), one common mode of bypass graft failure, has been shown to occur in the areas of disturbed flow particular to this site. The nature of theflow in the segment of artery proximal to the distal anastomosis varies from case to case depending on the clinical situation presented. A partial stenosis of a bypassed arterial segment may allow residual prograde flow through the proximal artery entering the distal anastomosis of the graft. A complete stenosis may allow for zero flow in the proximal artery segment or retrograde flow due to the presence of small collateral vessels upstream. Although a number of investigations on the hemodynamics at the distal anastomosis of an end-to-side bypass graft have been conducted, there has not been a uniform treatment of the proximal artery flow condition. As a result, direct comparison of results from study to study may not be appropriate. The purpose of this work was to perform a three-dimensional computational investigation to study the effect of the proximal artery flow condition (i.e., prograde, zero, and retrograde flow) on the hemodynamics at the distal end-to-side anastomosis. We used the finite volume method to solve the full Navier-Stokes equations for steady flow through an idealized geometry of the distal anastomosis. We calculated the flow field and local wall shear stress (WSS) and WSS gradient (WSSG) everywhere in the domain. We also calculated the severity parameter (SP), a quantification of hemodynamic variation, at the anastomosis. Our model showed a marked difference in both the magnitude and spatial distribution of WSS and WSSG. For example, the maximum WSS magnitude on the floor of the artery proximal to the anastomosis for the prograde and zero flow cases is 1.8 and 3.9 dynes/cm2, respectively, while it is increased to 10.3 dynes/cm2 in the retrograde flow case. Similarly, the maximum value of WSSG magnitude on thefloor of the artery proximal to the anastomosis for the prograde flow case is 4.9 dynes/cm3, while it is increased to 13.6 and 24.2 dynes/cm3, respectively, in the zero and retrograde flow cases. The value of SP is highest for the retrograde flow case (13.7 dynes/cm3) and 8.1 and 12.1 percent lower than this for the prograde (12.6 dynes/cm3) and zero (12.0 dynes/cm3) flow cases, respectively. Our model results suggest that the flow condition in the proximal artery is an important determinant of the hemodynamics at the distal anastomosis of end-to-side vascular bypass grafts. Because hemodynamic forces affect the response of vascular endothelial cells, the flow situation in the proximal artery may affect IH formation and, therefore, long-term graft patency. Since surgeons have some control over the flow condition in the proximal artery, results from this study could help determine which flow condition is clinically optimal.     

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

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

Sequential venous anastomosis design to enhance patency of arterio-venous grafts for hemodialysis

Arterio-venous grafts (AVGs), the second best option as long-term vascular access for hemodialysis, face major issues of stenosis mainly due to development of intimal hyperplasia at the venous anastomosis which is linked to unfavorable hemodynamic conditions. We have investigated computationally the utility of a coupled sequential venous anastomotic design to replace conventional end-to-side (ETS) venous anastomosis, in order to improve the hemodynamic environment and consequently enhance the patency of AVGs. Two complete vascular access models with the conventional and the proposed venous anastomosis configurations were constructed. Three-dimensional, pulsatile blood flow through the models was simulated, and wall shear stress (WSS)-based hemodynamic parameters were calculated and compared between the two models. Simulation results demonstrated that the proposed anastomotic design provides: (i) a more uniform and smooth flow at the ETS anastomosis, without flow impingement and stagnation point on the artery bed and vortex formation in the heel region of the ETS anastomosis; (ii) more uniform distribution of WSS and substantially lower WSS gradients on the venous wall; and (iii) a spare route for the blood flow to the vein, to avoid re-operation in case of stenosis. The distinctive hemodynamic advantages observed in the proposed anastomotic design can enhance the patency of AVGs.     

A New Near-wall Residence Time Model Applied to Three Arterio-venous Graft End-to-side Anastomoses

As a representative example of an end-to-side anastomosis, the distal portion of an arterio-venous hemodialysis graft with three graft-end designs and two outflow conditions has been considered. Arterio-venous graft failure for hemodialysis patients can be frequently related to stenotic developments, especially in the toe region. With critical blood particles such as LDL, platelets, and monocytes, being considered as a first-order approximation of; blood elements pathlines, near-wall residence times and wall stagnation points have been simulated to identify, in conjunction with the wall shear stress gradient, susceptible sites of stenotic developments. Assuming transient laminar three-dimensional non-Newtonian fluid flow, the transport equations were solved with validated finite-volume based CFD software and user-Fortran programs for the computation of fluid element pathlines and hemodynamic parameters. We conclude from this study that for the application of looped arterio-venous grafts, the graft-end geometry and the presence of retrograde outflow significantly influence the magnitude and location of the hemodynamic wall parameters considered. The new concept of near-wall blood “particle” behaviour, in conjunction with a wall shear stress based hemodynamic parameter, can be compared to similar histopathologic studies of intimal hyperplasia formation to understand better the role of hemodynamic factors in the biogenesis of the disease.     

Comparison of the hemodynamics in 6mm and 4-7 mm hemodialysis grafts by means of CFD

The aim of our study is to investigate with computational fluid dynamics (CFD) whether different arterial anastomotic geometries result in a different hemodynamics at the arterial (AA) and venous anastomosis (VA) of hemodialysis vascular access grafts. We have studied a 6mm graft (CD) and a 4-7 mm graft (TG). A validated three-dimensional CFD model is developed to simulate flow in the two graft types. Only the arterial anastomosis (AA) geometry differs. The boundary conditions applied are a periodic velocity signal at the arterial inlet and a periodic pressure wave at the venous outlet. Flow rate is set to 1,000 ml/min. The time dependent Navier-Stokes equations are solved. Wall shear stress (WSS), wall shear stress gradient (WSSG) and pressure gradient (PG) are calculated. Anastomotic flow is asymmetric although the anastomosis geometry is symmetric. The hemodynamic parameters, WSS, WSSG and PG, values at the suture line of the arterial anastomosis of the TG are at least twice as much as in the CD. Comparing the parameters at the two AA indicate that little flow rate increase introduces the risk of hemolysis in the TG whereas the CD is completely free of hemolysis. The hemodynamic parameter values at the venous anastomosis of the CD are 24 till 35% higher compared to the values of the TG. WSS values (> 3 Pa) in the VA are in the critical range for stenosis development in both graft geometries. The zones where the parameters reach extreme values correspond to the locations where intimal hyperplasia formation is reported in literature. In all anastomoses, the hemodynamic parameter levels are in the range where leucocytes and platelets get activated. Our simulations confirm clinical results where TG did not show a better outcome when compared to the CD.     

Simulation of particle-hemodynamics in a partially occluded artery segment with implications to the initiation of microemboli and secondary stenoses.

Computational results of laminar incompressible blood-particle flow analyses in an axisymmetric artery segment with a smooth local area constriction of 75 percent have been presented. The flow input waveform was sinusoidal with a nonzero average. The non-Newtonian behavior of blood was simulated with a modified Quemada model, platelet concentrations were calculated with a drift-flux model, and monocyte trajectories were described and compared for both Newtonian and Quemada rheologies. Indicators of "disturbed flow" included the time-averaged wall shear stress (WSS), the oscillatory shear index (OSI), and the wall shear stress gradient (WSSG). Implications of the vortical flow patterns behind the primary stenosis to the formation of microemboli and downstream stenoses are as follows. Elevated platelet concentrations due to accumulation in recirculation zones mixed with thrombin and ADP complexes assumed to be released upstream in high wall shear stress regions, could form microemboli, which are convected downstream. Distinct near-wall vortices causing a local increase in the WSSG and OSI as well as blood-particle entrainment with possible wall deposition, indicate sites susceptible to the onset of an additional stenosis proximal to the initial geometric disturbance.     

Numerical modeling of simulated blood flow in idealized composite arterial coronary grafts: steady state simulations

This paper presents a comparative study of simulated blood flow in different configurations of simplified composite arterial coronary grafts (CACGs). Even though the composite arterial grafting is increasingly used in cardiac surgery, it is still questionable whether or not the blood flow in such grafts can adequately meet the demands of the native myocardial circulation. A computational fluid dynamics (CFD) model was developed to conduct computer-based studies of simulated blood flow in four different geometric configurations of CACGs, corresponding to routinely used networks in cardiac surgery coronary grafts (T, Y, Pi and sequential). The flow was assumed three-dimensional, laminar and steady and the fluid as Newtonian, while the vessel walls were considered as inelastic and impermeable. It was concluded that local haemodynamics, practically described by velocity, pressure drop, wall shear stress (WSS) and flow rates, may be strongly influenced by the local geometry, especially at the anastomotic sites. The computations were made at mean flow rates of 37.5, 75 and 150ml/min. The side-branch outflow rates, computed for each bypass graft, showed noticeable differences. The results, which were found both qualitatively and quantitatively consistent with other studies, indicate that the Pi-graft exhibits significantly less uniform distribution of outflow rates than the other geometric configurations. Moreover, prominent variations in WSS and velocity distribution among the assessed CACGs were predicted, showing remarkable flow interactions among the arterial branches. The lowest shear stress regions were found on the lateral walls of bifurcations, which are predominantly susceptible to the occurrence of coronary artery disease (CAD). In contrast, the highest WSS were observed at the turn of the arterial branches.     

Optimization of cardiovascular stent design using computational fluid dynamics

Coronary stent design affects the spatial distribution of wall shear stress (WSS), which can influence the progression of endothelialization, neointimal hyperplasia, and restenosis. Previous computational fluid dynamics (CFD) studies have only examined a small number of possible geometries to identify stent designs that reduce alterations in near-wall hemodynamics. Based on a previously described framework for optimizing cardiovascular geometries, we developed a methodology that couples CFD and three-dimensional shape-optimization for use in stent design. The optimization procedure was fully-automated, such that solid model construction, anisotropic mesh generation, CFD simulation, and WSS quantification did not require user intervention. We applied the method to determine the optimal number of circumferentially repeating stent cells (N(C)) for slotted-tube stents with various diameters and intrastrut areas. Optimal stent designs were defined as those minimizing the area of low intrastrut time-averaged WSS. Interestingly, we determined that the optimal value of N(C) was dependent on the intrastrut angle with respect to the primary flow direction. Further investigation indicated that stent designs with an intrastrut angle of approximately 40 deg minimized the area of low time-averaged WSS regardless of vessel size or intrastrut area. Future application of this optimization method to commercially available stent designs may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes.     

Large Negative Stress Phase Angle (SPA) attenuates nitric oxide production in bovine aortic endothelial cells

Hemodynamics plays an important role in cardiovascular physiology and pathology. Pulsatile flow (Q), pressure (P), and diameter (D) waveforms exert wall shear stress (WSS), normal stress, and circumferential strain (CS) on blood vessels. Most in vitro studies to date have focused on either WSS or CS but not their interaction. Recently, we have shown that concomitant WSS and CS affect EC biochemical response modulated by the temporal phase angle between WSS and CS (stress phase angle, SPA). Large negative SPA has been shown to occur in regions of the circulation where atherosclerosis and intimal hyperplasia are prevalent. Here, we report that nitric oxide (NO) biochemical secretion was significantly decreased in response to a large negative SPA of -180 deg with respect to an SPA of 0 degrees in bovine aortic endothelial cells (BAEC) at 5 h. A new hemodynamic simulator for the study of the physiologic SPA was used to provide the hemodynamic conditions of pro-atherogenic (SPA = -180 deg) and normopathic (SPA = 0 deg) states. The role of complex hemodynamics in vascular remodeling, homeostasis, and pathogenesis can be advanced by further assessment of the hypothesis that a large negative SPA is pro-atherogenic.     

PIV-validated numerical modeling of pulsatile flows in distal coronary end-to-side anastomoses

This study employed particle image velocimetry (PIV) to validate a numerical model in a complementary approach to quantify hemodynamic factors in distal coronary anastomoses and to gain more insights on their relationship with anastomotic geometry. Instantaneous flow fields and wall shear stresses (WSS) were obtained from PIV measurement in a modified life-size silastic anastomosis model adapted from a conventional geometry by incorporating a smooth graft-artery transition. The results were compared with those predicted by a concurrent numerical model. The numerical method was then used to calculate cycle-averaged WSS (WSS(cyc)) and spatial wall shear stress gradient (SWSSG), two critical hemodynamic factors in the pathogenesis of intimal thickening (IT), to compare the conventional and modified geometries. Excellent qualitative agreement and satisfactory quantitative agreement with averaged normalized error in WSS between 0.8% and 8.9% were achieved between the PIV experiment and numerical model. Compared to the conventional geometry, the modified geometry produces a more uniform WSS(cyc) distribution eliminating both high and low WSS(cyc) around the toe, critical in avoiding IT. Peak SWSSG on the artery floor of the modified model is less than one-half that in the conventional case, and high SWSSG at the toe is eliminated. The validated numerical model is useful for modeling unsteady coronary anastomotic flows and elucidating the significance of geometry regulated hemodynamics. The results suggest the clinical relevance of constructing smooth graft-artery transition in distal coronary anastomoses to improve their hemodynamic performance.     

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.     

Numerical investigation of oxygen mass transfer in a helical-type artery bypass graft

Local oxygen lack in arterial walls (hypoxia) plays a very important role in the initiation, progression and development of intimal hyperplasia (IH) and thrombosis. Aiming to find out whether a helical-type artery bypass graft (ABG) is hypoxia beneficial, a numerical study was carried out to compare oxygen transport between a helical-type ABG and a conventional-type ABG. The dimensionless mass transfer coefficient (Sherwood number) was introduced to evaluate the oxygen mass transfer distribution and detailed oxygen wall flux was computed. The results show that the intrinsic geometry of a helical-type ABG resulted in improved hypoxia and the oxygen-depleted fluid located proximally to the occluded section as compared with that of a conventional-type ABG. However, benefits aside, distinct double low regions (low wall shear stress (WSS) and hypoxia) which might be most prone to IH and more localised and thicker boundary layer of oxygen-depleted fluid were observed at the helical-type ABG. This may explain why the helical flow plays a detrimental role at some locations in the human body. In addition, it was observed that although low WSS region was always accompanied with low oxygen supply, the oxygen transport rate did not adjust simultaneously with flow. The change in oxygen distribution usually lagged behind the flow change. A physiological WSS region may be associated with hypoxia condition. This study captured the qualitative trend of oxygen distribution in ABGs and the effect of helical geometry on reducing hypoxia, which is useful in the structural design of swirling flow vascular devices.     

Non-Newtonian effects of blood flow on hemodynamics in distal vascular graft anastomoses

Non-Newtonian fluid flow in a stenosed coronary bypass is investigated numerically using the Carreau-Yasuda model for the shear thinning behavior of the blood. End-to-side coronary bypass anastomosis is considered in a simplified model geometry where the host coronary artery has a 75% severity stenosis. Different locations of the bypass graft to the stenosis and different flow rates in the graft and in the host artery are studied. Particular attention is given to the non-Newtonian effect of the blood on the primary and secondary flow patterns in the host coronary artery and the wall shear stress (WSS) distribution there. Interaction between the jet flow from the stenosed artery and the flow from the graft is simulated by solving the three-dimensional Navier-Stokes equation coupled with the non-Newtonian constitutive model. Results for the non-Newtonian flow, the Newtonian flow and the rescaled Newtonian flow are presented. Significant differences in axial velocity profiles, secondary flow streamlines and WSS between the non-Newtonian and Newtonian fluid flows are revealed. However, reasonable agreement between the non-Newtonian and the rescaled Newtonian flows is found. Results from this study support the view that the residual flow in a partially occluded coronary artery interacts with flow in the bypass graft and may have significant hemodynamic effects in the host vessel downstream of the graft. Non-Newtonian property of the blood alters the flow pattern and WSS distribution and is an important factor to be considered in simulating hemodynamic effects of blood flow in arterial bypass grafts.     

Hemodynamic Analysis in an Idealized Artery Tree: Differences in Wall Shear Stress between Newtonian and Non-Newtonian Blood Models

Development of many conditions and disorders, such as atherosclerosis and stroke, are dependent upon hemodynamic forces. To accurately predict and prevent these conditions and disorders hemodynamic forces must be properly mapped. Here we compare a shear-rate dependent fluid (SDF) constitutive model, based on the works by Yasuda et al in 1981, against a Newtonian model of blood. We verify our stabilized finite element numerical method with the benchmark lid-driven cavity flow problem. Numerical simulations show that the Newtonian model gives similar velocity profiles in the 2-dimensional cavity given different height and width dimensions, given the same Reynolds number. Conversely, the SDF model gave dissimilar velocity profiles, differing from the Newtonian velocity profiles by up to 25% in velocity magnitudes. This difference can affect estimation in platelet distribution within blood vessels or magnetic nanoparticle delivery. Wall shear stress (WSS) is an important quantity involved in vascular remodeling through integrin and adhesion molecule mechanotransduction. The SDF model gave a 7.3-fold greater WSS than the Newtonian model at the top of the 3-dimensional cavity. The SDF model gave a 37.7-fold greater WSS than the Newtonian model at artery walls located immediately after bifurcations in the idealized femoral artery tree. The pressure drop across arteries reveals arterial sections highly resistive to flow which correlates with stenosis formation. Numerical simulations give the pressure drop across the idealized femoral artery tree with the SDF model which is approximately 2.3-fold higher than with the Newtonian model. In atherosclerotic lesion models, the SDF model gives over 1 Pa higher WSS than the Newtonian model, a difference correlated with over twice as many adherent monocytes to endothelial cells from the Newtonian model compared to the SDF model.     

Effects of different stent designs on local hemodynamics in stented arteries

Following the deployment of a coronary stent and disruption of an atheromatous plaque, the deformation of the arterial wall and the presence of the stent struts create a new fluid dynamic field, which can cause an abnormal biological response. In this study 3D computational models were used to analyze the fluid dynamic disturbances induced by the placement of a stent inside a coronary artery. Stents models were first expanded against a simplified arterial plaque, with a solid mechanics analysis, and then subjected to a fluid flow simulation under pulsatile physiological conditions. Spatial and temporal distribution of arterial wall shear stress (WSS) was investigated after the expansion of stents of different designs and different strut thicknesses. Common oscillatory WSS behavior was detected in all stent models. Comparing stent and vessel wall surfaces, maximum WSS values (in the order of 1Pa) were located on the stent surface area. WSS spatial distribution on the vascular wall surface showed decreasing values from the center of the vessel wall portion delimited by the stent struts to the wall regions close to the struts. The hemodynamic effects induced by two different thickness values for the same stent design were investigated, too, and a reduced extension of low WSS region (<0.5Pa) was observed for the model with a thicker strut.     

Procedure for Human Saphenous Veins Ex Vivo Perfusion and External Reinforcement

The mainstay of contemporary therapies for extensive occlusive arterial disease is venous bypass graft. However, its durability is threatened by intimal hyperplasia (IH) that eventually leads to vessel occlusion and graft failure. Mechanical forces, particularly low shear stress and high wall tension, are thought to initiate and to sustain these cellular and molecular changes, but their exact contribution remains to be unraveled. To selectively evaluate the role of pressure and shear stress on the biology of IH, an ex vivo perfusion system (EVPS) was created to perfuse segments of human saphenous veins under arterial regimen (high shear stress and high pressure). Further technical innovations allowed the simultaneous perfusion of two segments from the same vein, one reinforced with an external mesh. Veins were harvested using a no-touch technique and immediately transferred to the laboratory for assembly in the EVPS. One segment of the freshly isolated vein was not perfused (control, day 0). The two others segments were perfused for up to 7 days, one being completely sheltered with a 4 mm (diameter) external mesh. The pressure, flow velocity, and pulse rate were continuously monitored and adjusted to mimic the hemodynamic conditions prevailing in the femoral artery. Upon completion of the perfusion, veins were dismounted and used for histological and molecular analysis. Under ex vivo conditions, high pressure perfusion (arterial, mean = 100 mm Hg) is sufficient to generate IH and remodeling of human veins. These alterations are reduced in the presence of an external polyester mesh.