Abstract. The patency and overall performance of implanted bypass grafts is closely related to hemodynamics and its influence on vessel remodelling. In this regard, numerical investigation of blood flow in models reconstructed from clinical data may, next to clinical research, provide a valuable insight into the problem of graft failures, which are usually associated with restenosis and/or occlusive intimal hyperplasia. In this study, numerical results of pulsatile non-Newtonian blood flow in three realistic aorto-coronary bypass models are presented and further discussed with emphasis placed on the distribution of wall shear stress (WSS) and oscillatory shear index (OSI). Blood's shear-thinning behaviour is described by the CarreauYasuda model. Assuming all model walls to be impermeable and inelastic, the numerical solution of the mathematical model, which has the form of time-dependent non-linear system of Navier-Stokes (NS) equations, is carried out on the basis of the three-stage fractional step method and cell-centred finite volume method formulated for hybrid unstructured tetrahedral grids. The viscous terms of the NS equations are time discretised implicitly using the CrankNicolson scheme. The convective terms are solved explicitly and their computation utilises a local time-stepping technique in order to improve the overall computational efficiency of the developed CFD code.