Measurement and prediction of flow through a replica segment of a mildly atherosclerotic coronary artery of man

Pressure distributions were measured along a hollow vascular axisymmetric replica of a segment of the left circumflex coronary artery of man with mildly atherosclerotic diffuse disease. A large range of physiological Reynolds numbers from about 60 to 500, including hyperemic response, was spanned in the flow investigation using a fluid simulating blood kinematic viscosity. Predicted pressure distributions from the numerical solution of the Navier-Stokes equations were similar in trend and magnitude to the measurements. Large variations in the predicted velocity profiles occurred along the lumen. The influence of the smaller scale multiple flow obstacles along the wall (lesion variations) led to sharp spikes in the predicted wall shear stresses. Reynolds number similarity was discussed, and estimates of what time averaged in vivo pressure drop and shear stress might be were given for a vessel segment.     

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

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

Integration of geometric separation mechanisms by implementing curved constrictions in a biochip microchannel fluidic separator

This paper investigates the effectiveness of using curved constrictions in the bifurcation region of T-type fluid separators for promoting flow development in the intervals between bifurcations. A design of biofluid separator is proposed and a mathematical analysis and a numerical simulation of the blood flow in microchannels are conducted. The design is based on a modification of an existing T-shaped biochip device which consists of a main channel and a series of perpendicularly positioned side channels. By means of bifurcation effect, the blood is separated into plasma concentration flow from the side channels and blood cell concentration flow from the main channel. In this design, curved constrictions are inserted between bifurcations to replace the original straight channel section, so that the constriction and curved channel effects can be induced apart from the existing bifurcation effect. The mathematical analysis is aimed to the flow field and shear stress of the blood fluid in the microchannel geometries employed in the current design, including bifurcation, constriction and curved channel. The numerical simulation and mathematical analysis result in agreed conclusions, giving some insights into the importance of the relevant geometries in promoting biofluid separation. The main results can be summarised as follows: (i) the constrictions can largely increase the shear stress by the ratio of square of the reduction of the sections between the constriction and parent main channel. (ii) The curved channel intervals can induce centrifugal force, smoothly transit the flow field and increase the chances depleting fluid from the cell-free layer. (iii) The thickness of the boundary layer skimmed into the side channels from the main channel is decreased in this design and can be controlled, falling into the cell-free layer region by adjusting the geometry of the side channels.     

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.     

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

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

A numerical study of flow in curved tubes simulating coronary arteries

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

Stent implantation alters coronary artery hemodynamics and wall shear stress during maximal vasodilation.

Coronary stents improve resting blood flow and flow reserve in the presence of stenoses, but the impact of these devices on fluid dynamics during profound vasodilation is largely unknown. We tested the hypothesis that stent implantation affects adenosine-induced alterations in coronary hemodynamics and wall shear stress in anesthetized dogs (n = 6) instrumented for measurement of left anterior descending coronary artery (LAD) blood flow, velocity, diameter, and radius of curvature. Indexes of fluid dynamics and shear stress were determined before and after placement of a slotted-tube stent in the absence and presence of an adenosine infusion (1.0 mg/min). Adenosine increased blood flow, Reynolds (Re) and Dean numbers (De), and regional and oscillatory shear stress concomitant with reductions in LAD vascular resistance and segmental compliance before stent implantation. Increases in LAD blood flow, Re, De, and indexes of shear stress were observed after stent deployment (P < 0.05). Stent implantation reduced LAD segmental compliance to zero and potentiated increases in segmental and coronary vascular resistance during adenosine. Adenosine-induced increases in coronary blood flow and reserve, Re, De, and regional and oscillatory shear stress were attenuated after the stent was implanted. The results indicate that stent implantation blunts alterations in fluid dynamics during coronary vasodilation in vivo.     

Oxidative stress-induced dysregulation of arteriolar wall shear stress and blood pressure in hyperhomocysteinemia is prevented by chronic vitamin C treatment

We aimed to test the hypothesis that an enhanced level of reactive oxygen species (ROS) is primarily responsible for the impairment of nitric oxide (NO)-mediated regulation of arteriolar wall shear stress (WSS) in hyperhomocysteinemia (HHcy). Thus flow/WSS-induced dilations of pressurized gracilis muscle arterioles (basal diameter: approximately 170 microm) isolated from control (serum Hcy: 6 +/- 1 microM), methionine diet-induced HHcy rats (4 wk, serum Hcy: 30 +/- 6 microM), and HHcy rats treated with vitamin C, a known antioxidant (4 wk, 150 mg. kg body wt-1.day-1; serum Hcy: 32 +/- 10 microM), were investigated. In vessels of HHcy rats, increases in intraluminal flow/WSS-induced dilations were converted to constrictions. Constrictions were unaffected by inhibition of NO synthesis by N omega-nitro-L-arginine methyl ester (L-NAME). Vitamin C treatment of HHcy rats reversed the WSS-induced arteriolar constrictions to L-NAME-sensitive dilations but did not affect control responses. Similar changes in responses were obtained for the calcium ionophore A-23187. In addition, diastolic and mean arterial blood pressure and serum 8-isoprostane levels (a marker of in vivo oxidative stress) were significantly elevated in rats with HHcy, changes that were normalized by vitamin C treatment. Taken together, our data show that in chronic HHcy long-term vitamin C treatment, by decreasing oxidative stress in vivo, enhanced NO bioavailability, restored the regulation of shear stress in arterioles, and normalized systemic blood pressure. Thus our study provides evidence that oxidative stress is an important in vivo mechanism that is primarily responsible for the development of endothelial dysregulation of WSS in HHcy.     

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient ▿P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The ▿P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published ▿P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion.Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%.In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.     

Velocity measurements in steady flow through axisymmetric stenoses at moderate Reynolds numbers

The velocity field in the neighborhood of axisymmetric constrictions in rigid tubes was investigated using laser Doppler anemometry and flow visualization. Upstream flow conditions were steady; and Reynolds numbers were in the range 500-2000, values which are representative of the larger arteries in humans. Stenoses of 25, 50 and 75% area reduction were studied. Velocity profiles are presented in sufficient detail to allow comparison with computational biofluid dynamics models. Wall shear stresses were estimated from the near wall velocity gradient, and the nature of observed poststenotic flow disturbances is discussed. Results indicate that flow disturbances of discrete oscillation frequency may be more valuable than turbulence as an indicator of early stages of stenosis development. Additionally, despite the fact that poststenotic turbulence exists for the higher degrees of stenosis and Reynolds numbers, the resulting wall shear stresses are only three to four times greater than the Poiseuille value and are considerably less than the wall shear stress within the stenosis itself.     

Haemodynamic assessment of human coronary arteries is affected by degree of freedom of artery movement

Abnormal haemodynamic parameters are associated with atheroma plaque progression and instability in coronary arteries. Flow recirculation, shear stress and pressure gradient are understood to be important pathogenic mediators in coronary disease. The effect of freedom of coronary artery movement on these parameters is still unknown. Fluid–structure interaction (FSI) simulations were carried out in 25 coronary artery models derived from authentic human coronaries in order to investigate the effect of degree of freedom of movement of the coronary arteries on flow recirculation, wall shear stress (WSS) and wall pressure gradient (WPG). Each FSI model had distinctive supports placed upon it. The quantitative and qualitative differences in flow recirculation, maximum wall shear stress (MWSS), areas of low wall shear stress (ALWSS) and maximum wall pressure gradient (MWPG) for each model were determined. The results showed that greater freedom of movement was associated with lower MWSS, smaller ALWSS, smaller flow recirculation zones and lower MWPG. With increasing percentage diameter stenosis (%DS), the effect of degree of freedom on flow recirculation and WSS diminished. Freedom of movement is an important variable to be considered for computational modelling of human coronary arteries, especially in the setting of mild to moderate stenosis.

Abbreviations: 3D: Three-dimensional; 3DR: Three-dimensional Reconstruction; 3D-QCA: Three-dimensional quantitative coronary angiography; ALWSS: Areas of low wall shear stress; CAD: Coronary artery disease; CFD: Computational fluid dynamics; %DS: Diameter stenosis percentage; EPCS: End point of counter-rotating streamlines; FSI: Fluid–structure interaction; IVUS: Intravascular ultrasound; LAD: Left anterior descending; MWSS: Maximum wall shear stress; SST: Shear stress transport; TAWSS: Time-averaged wall shear stress; WSS: wall shear stress; WPG: Wall pressure gradient; MWPG: Maximum wall pressure gradient; FFR: Fractional flow reserve; iFR: Instantaneous wave-free ratio     

Influence of stenosis on hemodynamic parameters in the realistic left coronary artery under hyperemic conditions

The current study investigates the hyperemic flow effects on heamodynamics parameters such as velocity, wall shear stress in 3D coronary artery models with and without stenosis. The hyperemic flow is used to evaluate the functional significance of stenosis in the current era. Patients CT scan data of having healthy and coronary artery disease was chosen for the reconstruction of 3D coronary artery models. The diseased 3D models of coronary artery shows a narrowing of >50% lumen area. Computational fluid dynamics was performed to simulate the hyperemic flow condition. The results showed that the recirculation zone was observed immediate to the stenosis and highest wall shear stress was observed across the stenosis. The decrease in pressure was found downstream to the stenosis as compared to the coronary artery without stenosis. Our analysis provides an insight into the distribution of wall shear stress and pressure drop, thus improving our understanding of hyperemic flow effect under both conditions.     

Effect of geometrical assumptions on numerical modeling of coronary blood flow under normal and disease conditions

Shear stress plays a pivotal role in pathogenesis of coronary heart disease. The spatial and temporal variation in hemodynamics of blood flow, especially shear stress, is dominated by the vessel geometry. The goal of the present study was to investigate the effect of 2D and 3D geometries on the numerical modeling of coronary blood flow and shear stress distribution. We developed physiologically realistic 2D and 3D models (with similar geometries) of the human left coronary artery under normal and stenosis conditions (30%, 60%, and 80%) using PROE (WF 3). Transient blood flows in these models were solved using laminar and turbulent (k-ω) models using a computational fluid dynamics solver, FLUENT (v6.3.26). As the stenosis severity increased, both models predicted a similar pattern of increased shear stress at the stenosis throat, and in recirculation zones formed downstream of the stenosis. The 2D model estimated a peak shear stress value of 0.91, 2.58, 5.21, and 10.09 Pa at the throat location under normal, 30%, 60%, and 80% stenosis severity. The peak shear stress values at the same location estimated by the 3D model were 1.41, 2.56, 3.15, and 13.31 Pa, respectively. The 2D model underestimated the shear stress distribution inside the recirculation zone compared with that of 3D model. The shear stress estimation between the models diverged as the stenosis severity increased. Hence, the 2D model could be sufficient for analyzing coronary blood flow under normal conditions, but under disease conditions (especially 80% stenosis) the 3D model was more suitable.     

The effect of shear on the desorption of liposomes adsorbed to bacterial biofilms

With the aid of a flow cell assembly the desorption of cationic liposomes prepared from mixtures of dipalmitoylphoshatidylcholine (DDPC), cholesterol, and either dimethyldioctadecylammonium bromide (DDAB) or 3,beta[N-(N1,N-dimethylethylenediamine)-carbamoyl]cholesterol (DC-chol) from immobilized biofilms of Staphylococcus aureus has been studied as a function of shear stress by confocal microscopy. A shear stress theory has been adapted from fluid mechanics of laminar flow between parallel plates and used to determine the critical shear stress for liposome desorption. The critical shear stress for both DDAB and DC-chol liposomes has been determined as a function of cationic lipid content and hence surface charge as reflected in their zeta potentials. The critical shear stress has been used to obtain the potential energy of liposome-biofilm interaction which together with the electrostatic interaction energy has enabled estimates of the London-Hamaker constants to be made. The values of the London-Hamaker constants at small liposome-bacterial cell separation were found to be independent of liposome composition.     

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

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

Effects of luminal shear stress on cerebral arteries and arterioles

The effect of luminal shear stress was studied in cerebral arteries and arterioles. Middle cerebral arteries (MCA) and penetrating arterioles (PA) were isolated from male Long-Evans rats, mounted in a tissue bath, and pressurized. After the development of spontaneous tone, inside diameters were 186 +/- 5 microm (n = 28) for MCA and 65 +/- 3 microm (n = 37) for PA. MCA and PA constricted approximately 20% with increasing flow. Flow-induced constriction persisted in MCA and PA after removal of the endothelium. After removal of the endothelium, the luminal application of a polypeptide containing the Arg-Gly-Asp amino acid sequence (inhibitor of integrin attachment) abolished the flow-induced constriction. Similarly, an antibody specific for the beta(3)-chain of the integrin complex significantly inhibited the flow-induced constriction. The shear stress-induced constriction was accompanied by an increase in vascular smooth muscle Ca(2+). For example, a shear stress of 20 dyn/cm(2) constricted MCA 8% (n = 5) and increased Ca(2+) from 209 +/- 17 to 262 +/- 29 nM (n = 5). We conclude that isolated cerebral arteries and arterioles from the rat constrict to increased shear stress. Because the endothelium is not necessary for the response, the shear forces must be transmitted across the endothelium, presumably by the cytoskeletal matrix, to elicit constriction. Integrins containing the beta(3)-chain are involved with the shear stress-induced constrictions.     

Computation of hemodynamics in the left coronary artery with variable angulations

The purpose of this study was to investigate the hemodynamic effect of variations in the angulations of the left coronary artery, based on simulated and realistic coronary artery models. Twelve models consisting of four realistic and eight simulated coronary artery geometries were generated with the inclusion of left main stem, left anterior descending and left circumflex branches. The simulated models included various coronary artery angulations, namely, 15°, 30°, 45°, 60°, 75°, 90°, 105° and 120°. The realistic coronary angulations were based on selected patient's data with angles ranging from narrow angles of 58° and 73° to wide angles of 110° and 120°. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics. The wall shear stress, wall shear stress gradient, velocity flow patterns and wall pressure were measured in simulated and realistic models during the cardiac cycle. Our results showed that a disturbed flow pattern was observed in models with wider angulations, and wall pressure was found to reduce when the flow changed from the left main stem to the bifurcated regions, based on simulated and realistic models. A low wall shear stress gradient was demonstrated at left bifurcations with wide angles. There is a direct correlation between coronary angulations and subsequent hemodynamic changes, based on realistic and simulated models. Further studies based on patients with different severities of coronary artery disease are required to verify our results.     

Reduced NO-dependent arteriolar dilation during the development of cardiomyopathy

Our previous studies have suggested that there is reduced nitric oxide (NO) production in canine coronary blood vessels after the development of pacing-induced heart failure. The goal of these studies was to determine whether flow-induced NO-mediated dilation is altered in coronary arterioles during the development of heart failure. Subepicardial coronary arterioles (basal diameter 80 microm) were isolated from normal canine hearts, from hearts with dysfunction but no heart failure, and from hearts with severe cardiac decompensation. Arterioles were perfused at increasing flow or administered agonists with no flow in vitro. In arterioles from normal hearts, flow increased arteriolar diameter, with one-half of the response being NO dependent and one-half prostaglandin dependent. Shear stress-induced dilation was eliminated by removing the endothelium. Arterioles from normal hearts and hearts with dysfunction but no failure responded to increasing shear stress with dilation that reached a maximum at a shear stress of 20 dyn/cm(2). In contrast, arterioles from failing hearts showed a reduced dilation, reaching only 55% of the dilation seen in vessels of normal hearts at a shear stress of 100 dyn/cm(2). This remaining dilation was eliminated by indomethacin, suggesting that the NO-dependent component was absent in coronary microvessels after the development of heart failure. Similarly, agonist-induced NO-dependent coronary arteriolar dilation was markedly attenuated after the development of heart failure. After the development of severe dilated cardiomyopathy and heart failure, the NO-dependent component of both shear stress- and agonist-induced arteriolar dilation is reduced or entirely absent.     

Three-dimensional numerical simulations of flow through a stenosed coronary bypass

This work analyzes the flow patterns at the anastomosis of a stenosed coronary bypass. Three-dimensional numerical simulations are performed using a finite elements method. We consider a geometrical model of the host coronary artery with and without a 75% severity stenosis for three different locations from the anastomosis. The flow features - velocity profiles, secondary motions and wall shear stresses - are compared for different configurations of the flow rate and of the distance of the anastomosis from the site of occlusion (called distance of grafting). The combination of the junction flow effects - counter rotating vortices - with the stenosis effects - confined jet flow - is particularly important when the distance of grafting is short. Given that the residual flow issued from the pathologic stenosis being non-negligible after two weeks grafting, models without stenosis cannot predict the evolution of the wall shear stress in the vicinity of the anastomosis.     

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.