Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

Analysis of blood flow in an out-of-plane CABG model

Coronary artery bypass graft (CABG) is a routine surgical treatment for ischemic and infarcted myocardium. A large number of CABG fail postoperatively because of intimal hyperplasia within months or years. The cause of this failure is thought to be partly related to the flow patterns and shear stresses acting on the endothelial cells. An accurate representation of the flow field and associated wall shear stress (WSS) requires a detailed three-dimensional (3D) model of the CABG. The purpose of this study is to present a detailed analysis of blood flow in a 3D aorto/left CABG, bypassing the occluded left anterior descending coronary (LAD) artery. The analysis takes into account the influence of the out-of-plane geometry of the graft. The finite volume technique was employed to model the 3D blood flow pattern to determine the velocity and WSS distributions. This study presents the flow field distributions of the velocity and WSS at four instances of the cardiac cycle, two in systole and two in diastole. Our results reveal that the CABG geometry has a significant effect on the velocity distribution. The axial velocity profiles at different instances of the cardiac cycle exhibit strong skewing; significant secondary flow and vortex structures are seen in the in-plane velocity patterns. The maximum WSS on the bed of the occluded LAD artery opposite to the graft junction is 14 Pa in middiastole, whereas there is a significantly lower and more uniform distribution of WSS on the bed of the anastomosis. The present results indicate that nonplanarity of the blood vessel along with the inflow conditions has a substantial effect on the fluid mechanics of CABG that contribute to the patency of graft.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

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.     

Syringomyelia hydrodynamics: an in vitro study based on in vivo measurements

A simplified in vitro model of the spinal canal, based on in vivo magnetic resonance imaging, was used to examine the hydrodynamics of the human spinal cord and subarachnoid space with syringomyelia. In vivo magnetic resonance imaging (MRI) measurements of subarachnoid (SAS) geometry and cerebrospinal fluid velocity were acquired in a patient with syringomyelia and used to aid in the in vitro model design and experiment. The in vitro model contained a fluid-filled coaxial elastic tube to represent a syrinx. A computer controlled pulsatile pump was used to subject the in vitro model to a CSF flow waveform representative of that measured in vivo. Fluid velocity was measured at three axial locations within the in vitro model using the same MRI scanner as the patient study. Pressure and syrinx wall motion measurements were conducted external to the MR scanner using the same model and flow input. Transducers measured unsteady pressure both in the SAS and intra-syrinx at four axial locations in the model A laser Doppler vibrometer recorded the syrinx wall motion at 18 axial locations and three polar positions. Results indicated that the peak-to-peak amplitude of the SAS flow waveform in vivo was approximately tenfold that of the syrinx and in phase (SAS approximately 5.2 +/- 0.6 ml/s, syrinx approximately 0.5 +/- 0.3 ml/s). The in vitro flow waveform approximated the in vivo peak-to-peak magnitude (SAS approximately 4.6 +/- 0.2 ml/s, syrinx approximately 0.4 +/- 0.3 ml/s). Peak-to-peak in vitro pressure variation in both the SAS and syrinx was approximately 6 mm Hg. Syrinx pressure waveform lead the SAS pressure waveform by approximately 40 ms. Syrinx pressure was found to be less than the SAS for approximately 200 ms during the 860-ms flow cycle. Unsteady pulse wave velocity in the syrinx was computed to be a maximum of approximately 25 m/s. LDV measurements indicated that spinal cord wall motion was nonaxisymmetric with a maximum displacement of approximately 140 microm, which is below the resolution limit of MRI. Agreement between in vivo and in vitro MR measurements demonstrates that the hydrodynamics in the fluid filled coaxial elastic tube system are similar to those present in a single patient with syringomyelia. The presented in vitro study of spinal cord wall motion, and complex unsteady pressure and flow environment within the syrinx and SAS, provides insight into the complex biomechanical forces present in syringomyelia.     

The influence of inflow boundary conditions on intra left ventricle flow predictions

The combination of computational fluid dynamics (CFD) and magnetic resonance imaging (MRI) offers a promising tool that enables the prediction of blood flow patterns in subject-specific cardiovascular models. The influence of the model geometry on the accuracy of the simulation is well recognized. This paper addresses the impact of different boundary conditions on subject-specific simulations of left ventricular (LV) flow. A novel hybrid method for prescribing effective inflow boundary conditions in the mitral valve plane has been developed. The detailed quantitative results highlight the strengths as well as the potential pitfalls of the approach.     

A numerical calculation of flow in a curved tube model of the left main coronary artery

The flow pattern in the left main coronary artery has been calculated using an idealized geometry and by numerically solving the full Navier-Stokes equations for a Newtonian fluid. Two different forms for the entrance velocity profile were used, one a time-varying, flat profile and the other a time-varying, less flat velocity profile. The results obtained demonstrate the presence of secondary motions for conditions simulating flow in the left main coronary artery, with maximum secondary flow velocities being on the order of three to four percent of the maximum axial velocity. This secondary flow phenomenon has an important influence on the wall shear stress distribution, in spite of the fact that there is virtually no alteration in the axial velocity profile. The maximum ratio of the outer wall shear stress to that on the inner wall is 1.4 at a Reynolds number of Re = 270, and it increases with increasing Reynolds number, reaching a value of 1.7 at Re = 810. Although there are significant differences in the results in the immediate vicinity of the inlet for the two different forms of the entrance velocity profile used, this difference does not persist far into the tube. Independent of the choice of the entrance velocity profile, it appears that there will be significant secondary flow effects on the wall shear stress.     

A numerical simulation of steady flow fields in a bypass tube.

Steady flow in a complete by-pass tube was simulated numerically. The study was to consider a complete flow field, which included both the by-pass and the host tubes. The changes of the hemodynamics were investigated with three parameters: the inlet flow Reynolds number (Re), anastomotic angle (alpha) and the position of the occlusion in the host tube. The baseline flow field was set up with Re=200, alpha=45 degrees and the centered position of occlusion. The parametric study was then conducted on combination of Re=100, 200, 400, alpha=35 degrees, 45 degrees, 60 degrees, 75 degrees, 90 degrees and three occlusion positions: left, center and right. It was found that in the baseline case, large slow/recirculation flows could be seen in the host tube both upstream and downstream of the occlusion. The separation points were on the opposite walls to the junctions. Recirculation zones were also found near the toe and in the proximal outer wall of the by-pass tube. Their sizes were about one diameter of the tube or smaller. In some cases, pairing vortices could be seen in the host tube upstream of the occlusion. The shear rate distribution associated with the flow fields was presented. The flow pattern obtained was agreeable to those observed experimentally by other investigators. The difference of the flow fields between a complete bypass and simple anastomosis was discussed. The present numerical code provides a preliminary simulation/design tool for bypass graft flows.     

Computational investigation of subject-specific cerebrospinal fluid flow in the third ventricle and aqueduct of Sylvius

The cerebrospinal fluid flow in the third ventricle of the brain and the aqueduct of Sylvius was studied using computational fluid dynamics (CFD) based on subject-specific boundary conditions derived from magnetic resonance imaging (MRI) scans. The flow domain geometry was reconstructed from anatomical MRI scans by manual image segmentation. The movement of the domain boundary was derived from MRI brain motion scans. Velocimetric MRI scans were used to reconstruct the velocity field at the inferior end of the aqueduct of Sylvius based on the theory of pulsatile flow in pipes. A constant pressure boundary condition was assigned at the foramina of Monro. Three main flow features were observed: a fluid jet emerging from the aqueduct of Sylvius, a moderately mobile recirculation zone above the jet and a mobile recirculation below the jet. The flow in the entire domain was laminar with a maximum Reynolds number of 340 in the aqueduct. The findings demonstrate that by combining MRI scans and CFD simulations, subject-specific detailed quantitative information of the flow field in the third ventricle and the aqueduct of Sylvius can be obtained.     

Transitional flow at the venous anastomosis of an arteriovenous graft: potential activation of the ERK1/2 mechanotransduction pathway

We present experimental and computational results that describe the level, distribution, and importance of velocity fluctuations within the venous anastomosis of an arteriovenous graft. The motivation of this work is to understand better the importance of biomechanical forces in the development of intimal hyperplasia within these grafts. Steady-flow in vitro studies (Re = 1060 and 1820) were conducted within a graft model that represents the venous anastomosis to measure velocity by means of laser Doppler anemometry. Numerical simulations with the same geometry and flow conditions were conducted by employing the spectral element technique. As flow enters the vein from the graft, the velocity field exhibits flow separation and coherent structures (weak turbulence) that originate from the separation shear layer. We also report results of a porcine animal study in which the distribution and magnitude of vein-wall vibration on the venous anastomosis were measured at the time of graft construction. Preliminary molecular biology studies indicate elevated activity levels of the extracellular regulatory kinase ERK1/2, a mitogen-activated protein kinase involved in mechanotransduction, at regions of increased vein-wall vibration. These findings suggest a potential relationship between the associated turbulence-induced vein-wall vibration and the development of intimal hyperplasia in arteriovenous grafts. Further research is necessary, however, in order to determine if a correlation exists and to differentiate the vibration effect from that of flow related effects.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Hydrodynamic modeling of cerebrospinal fluid motion within the spinal cavity

The fluid that resides within cranial and spinal cavities, cerebrospinal fluid (CSF), moves in a pulsatile fashion to and from the cranial cavity. This motion can be measured hy magnetic resonance imaging (MRI) and may he of clinical importance in the diagnosis of several brain and spinal cord disorders such as hydrocephalus, Chiari malformation, and syringomyelia. In the present work, a geometric and hydrodynamic characterization of an anatomically relevant spinal canal model is presented. We found that inertial effects dominate the flow field under normal physiological flow rates. Along the length of the spinal canal, hydraulic diameter was found to vary significantly from 5 to 15 mm. The instantaneous Reynolds number at peak flow rate ranged from 150 to 450, and the Womersle number ranged from 5 to 17. Pulsatile flow calculations are presented for an idealized geometric representation of the spinal cavity. A linearized Navier-Stokes model of the pulsatile CSF flow was constructed based on MRI flow rate measurements taken on a healthy volunteer. The numerical model was employed to investigate effects of cross-sectional geometry and spinal cord motion on unsteady velocity, shear stress, and pressure gradientfields. The velocity field was shown to be blunt, due to the inertial character of the flow, with velocity peaks located near the boundaries of the spinal canal rather than at the midpoint between boundaries. The pressure gradient waveform was found to be almost exclusively dependent on the flow waveform and cross-sectional area. Characterization of the CSF dynamics in normal and diseased states may be important in understanding the pathophysiology of CSF related disorders. Flow models coupled with MRI flow measurements mnay become a noninvasive tool to explain the abnormal dynamics of CSF in related brain disorders as well as to determine concentration and local distribution of drugs delivered into the CSF space.     

Inlet conditions for image-based CFD models of the carotid bifurcation: is it reasonable to assume fully developed flow?

BACKGROUND: Computational fluid dynamics tools are useful for their ability to model patient specific data relevant to the genesis and progression of atherosclerosis, but unavailable to measurement tools. The sensitivity of the physiologically relevant parameters of wall shear stress (WSS) and the oscillatory shear index (OSI) to secondary flow in the inlet velocity profiles was investigated in three realistic models of the carotid bifurcation. METHOD OF APPROACH: Secondary flow profiles were generated using sufficiently long entrance lengths, to which curvature and helical pitch were added. The differences observed were contextualized with respect to effect of the uncertainty of the models' geometry on the same parameters. RESULTS: The effects of secondary velocities in the inlet profile on WSS and OSI break down within a few diameters of the inlet. Overall, the effect of secondary inlet flow on these models was on average more than 3.5 times smaller than the effect of geometric variability, with 13% and 48% WSS variability induced by inlet secondary flow and geometric differences, respectively. CONCLUSIONS: The degree of variation is demonstrated to be within the range of the other computational assumptions, and we conclude that given a sufficient entrance length of realistic geometry, simplification to fully developed axial (i.e., Womersley) flow may be made without penalty. Thus, given a choice between measuring three components of inlet velocity or a greater geometric extent, we recommend effort be given to more accurate and detailed geometric reconstructions, as being of primary influence on physiologically significant indicators.     

MRI-based modeling of spleno-mesenteric confluence flow

Characterization of hepatic blood flow magnitude and distribution can lead to a better understanding of the pathophysiology of liver disease. However, the underlying patterns and dynamics of hepatic flow, such as the helical flow structure that often develops following the spleno-mesenteric confluence (SMC) of the hepatic portal vein, have not yet been comprehensively studied. In this study, we used magnetic resonance image (MRI)-based computational models to study the effects of the helical flow structure and SMC geometry on portal blood flow distribution. Additionally, we examined these flow dynamics with four-dimensional (4D) flow MRI in a group of 12 cirrhotic patients and healthy subjects. A validation model was also created to compare computational data to particle image velocimetry (PIV) data. We found significant correlations between flow structure development, vessel geometry, and blood flow distribution in both virtually modified models and in healthy and cirrhotic subjects. However, the direction of these correlations varied among vessel configuration types. Nonetheless, validation model results displayed good qualitative agreement with computational model data.     

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.     

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.     

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.

     

Helical flow as fluid dynamic signature for atherogenesis risk in aortocoronary bypass. A numeric study

The main purpose of the study was to verify if helical flow, widely observed in several vessels, might be a signature of the blood dynamics of vein graft anastomosis. We investigated the existence of a relationship between helical flow structures and vascular wall indexes of atherogenesis in aortocoronary bypass models with different geometric features. In particular, we checked for the existence of a relationship between the degree of helical motion and the magnitude of oscillating shear stress in conventional hand-sewn proximal anastomosis. The study is based on the numerical evaluation of four bypass geometries that are attached to a simplified computer representation of the ascending aorta with different angulations relative to aortic outflow. The finite volume technique was used to simulate realistic graft fluid dynamics, including aortic compliance and proper aortic and graft flow rates. A quantitative method was applied to evaluate the level of helicity in the flow field associated with the four bypass models under investigation. A linear inverse relationship (R = -0.97) was found between the oscillating shear index and the helical flow index for the models under investigation. The results obtained support the hypothesis that an arrangement of the flow field in helical patterns may elicit damping in wall shear stress temporal gradients at the proximal graft. Accordingly, helical flow might play a significant role in preventing plaque deposition or in tuning the mechanotransduction pathways of cells. Therefore, results confirm that helical flow constitutes an important flow signature in vessels, and its strength as a fluid dynamic index (for instance in combination with magnetic resonance imaging flow visualization techniques) for risk stratification, in the activation of both mechanical and biological pathways leading to fibrointimal hyperplasia.