Numerical simulation of unsteady generalized Newtonian blood flow through differently shaped distensible arterial stenoses

An updated numerical simulation of unsteady generalized Newtonian blood flow through differently shaped distensible arterial stenoses is developed. A shear-thinning fluid modelling the deformation dependent viscosity of blood is considered for the characterization of generalized Newtonian behaviour of blood. The arterial model is treated as two-dimensional and axisymmetric with an outline of the stenosis obtained from a three-dimensional casting of a mildly stenosed artery. The full Navier-Stokes equations governing blood flow are written in the dimensionless form and the solution is accomplished by finite time-step advancement through their finite difference staggered grid representations. The marker and cell (MAC) method comprising the use of a set of marker particles moving with the fluid is used for the purpose. Results are obtained for three differently shaped stenoses – irregular, smooth and cosine curve representations. The present results do agree well with those of existing investigations in the steady state, but contrary to their conclusions the present findings demonstrate that the excess pressure drop across the cosine and the smooth stenoses is caused by neither their smoothness nor their higher degree of symmetry relative to the irregular stenosis, but is rather an effect of area cover with respect to the irregular stenosis. This effect clearly prevails throughout the entire physiological range of Reynolds numbers. Further the in-depth study in flow patterns reveals the development of flow separation zones in the diverging part of the stenosis towards the arterial wall, and they are influenced by non-Newtonian blood rheology, distensibility of the wall and flow unsteadiness in order to validate the applicability of the present model.