Influence of Shear-Thinning Blood Rheology on the Laminar-Turbulent Transition over a Backward Facing Step

• Published in: Fluids (Basel) 2020-06, Vol.5 (2), p.57
• Main Authors: Kelly, Nathaniel S., Gill, Harinderjit S., Cookson, Andrew N., Fraser, Katharine H.
• Format: Article
• Language: English
• Subjects: Artificial organs; Backward facing steps; blood; Blood flow; Cardiovascular diseases; Computational fluid dynamics; Computer simulation; Coronary vessels; Eddies; Flow simulation; Fluid flow; Fluids; Geometry; Laminar flow; Large eddy simulation; Mathematical models; Mechanical failure; Newtonian; Newtonian fluids; non-Newtonian; Non-Newtonian fluids; Parameters; Proteins; Reynolds number; Rheological properties; Rheology; Shear flow; Shear thinning (liquids); shear-thinning; Simulation; Thinning; Thromboembolism; Thrombosis; transition; Turbulence models; Turbulent flow; Velocity; Viscosity; Yield stress
• ISSN: 2311-5521; 2311-5521
• DOI: 10.3390/fluids5020057

Cardiovascular diseases are the leading cause of death globally and there is an unmet need for effective, safer blood-contacting devices, including valves, stents and artificial hearts. In these, recirculation regions promote thrombosis, triggering mechanical failure, neurological dysfunction and infarctions. Transitional flow over a backward facing step is an idealised model of these flow conditions; the aim was to understand the impact of non-Newtonian blood rheology on modelling this flow. Flow simulations of shear-thinning and Newtonian fluids were compared for Reynolds numbers ( R e ) covering the comprehensive range of laminar, transitional and turbulent flow for the first time. Both unsteady Reynolds Averaged Navier–Stokes ( k − ω SST) and Smagorinsky Large Eddy Simulations (LES) were assessed; only LES correctly predicted trends in the recirculation zone length for all R e . Turbulent-transition was assessed by several criteria, revealing a complex picture. Instantaneous turbulent parameters, such as velocity, indicated delayed transition: R e = 1600 versus R e = 2000, for Newtonian and shear-thinning transitions, respectively. Conversely, when using a Re defined on spatially averaged viscosity, the shear-thinning model transitioned below the Newtonian. However, recirculation zone length, a mean flow parameter, did not indicate any difference in the transitional Re between the two. This work shows a shear-thinning rheology can explain the delayed transition for whole blood seen in published experimental data, but this delay is not the full story. The results show that, to accurately model transitional blood flow, and so enable the design of advanced cardiovascular devices, it is essential to incorporate the shear-thinning rheology, and to explicitly model the turbulent eddies.