Numerical simulation of laminar non-Newtonian blood flow under varying Reynolds numbers and geometric parameters using OpenFOAM

Year 2025, Volume: 3 Issue: 2, 150 - 163, 31.10.2025

Mehmet Numan Kaya

https://doi.org/10.59292/bulletinbiomath.1791463

Abstract

Blood flow within biomedical devices and vascular models is characterized by laminar dynamics at relatively low Reynolds numbers, where shear-dependent viscosity governs hemodynamic behavior. This behavior is commonly observed in biomedical applications, such as nozzles, and accurate modeling is essential. In this study, numerical simulations of laminar non-Newtonian blood flow are performed in OpenFOAM using the Bird-Carreau viscosity model to examine the influence of flow conditions and nozzle geometry. The FDA benchmark nozzle is employed as a reference geometry, and the computational setup is validated against available experimental data prior to the parametric study. Five throat Reynolds numbers, $Re = 100, 300, 500, 1000,$ and $1500$, are investigated together with two collector cone angles, $20^\circ$ and $40^\circ$, as well as three throat diameters, $D_t = 3, 4,$ and $5$ mm, to assess geometric effects. The results show that narrower throats and higher Reynolds numbers significantly increase both velocity and shear stress, highlighting the strong sensitivity of hemodynamics to geometric constriction. Pressure drop analysis further reveals that enlarging the throat diameter can substantially reduce losses; for example, at $Re=1000$ and $1500$, increasing the throat from 3 mm to 4 mm lowers the pressure drop by nearly 50%, while a further increase to 5 mm reduces it by about 40%. Overall, the study demonstrates that both geometric variations and flow conditions lead to significant changes in blood flow physics, underscoring their importance in hemodynamic applications.

Keywords

Laminar blood flow , Non-Newtonian fluid , Bird–Carreau model , CFD , OpenFOAM , blood , Reynolds number , Cone angle , Throat diameter

References

• [1] Onder, A., Incebay, O. and Yapici, R. Computational fluid dynamics simulating of the FDA benchmark blood pump with different coefficient sets and scaler shear stress models used in the power-law hemolysis model. Journal of Artificial Organs, 28, 184-191, (2025).
• [2] Mustafaoğlu, M., Kotçioğlu, I. and Yeşilyurt, M.K. Numerical analysis of the three-dimensional model of pulsatile and non-Newtonian blood flow in a carotid artery with local occlusion. Mathematical Modelling and Numerical Simulation With Applications, 5(1), 97-116, (2025).
• [3] Avci, M. CFD-Based extensional stress and hemolysis risk evaluation in the U.S. food and drug administration (FDA) benchmark nozzle configurations. Fluids, 10(9), 224, (2025).
• [4] Liu, H., Lan, L., Abrigo, J., Ip, H.L., Soo, Y., Zheng, D. et al. Comparison of Newtonian and non-Newtonian fluid models in blood flow simulation in patients with intracranial arterial stenosis. Frontiers in Physiology, 12, 718540, (2021).
• [5] Mamun, K., Rahman, M.M., Akhter, M.N. and Ali, M. Physiological non-Newtonian blood flow through single stenosed artery. AIP Conference Proceedings 1754(1), 040001, (2016).
• [6] Mehrabi, M. and Setayeshi, S. Computational fluid dynamics analysis of pulsatile blood flow behavior in modelled stenosed vessels with different severities. Mathematical Problems in Engineering, 2012(1), 804765, (2012).
• [7] Çutay, A., Çermik, Ö. and Kaya, A. Evaluation of the vessel with different stenosis structures using CFD approach. Kahramanmaraş Sütçü Imam Üniversitesi Mühendislik Bilimleri Dergisi, 28(1), 245-257, (2025).
• [8] Gijsen, F.J.H., van de Vosse, F.N. and Janssen, J.D. The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model. Journal of Biomechanics, 32(6), 601-608, (1999).
• [9] Lakzian, E. and Akbarzadeh, P. Numerical investigation of unsteady pulsatile Newtonian/non-Newtonian blood flow through curved stenosed arteries. Bio-Medical Materials and Engineering, 30(5–6), 525-540, (2020).
• [10] Fuchs, A., Berg, N., Fuchs, L. and Prahl Wittberg, L. Assessment of rheological models applied to blood flow in human thoracic aorta. Bioengineering, 10(11), 1240, (2023).
• [11] Faraji, A., Sahebi, M. and SalavatiDezfouli, S. Numerical investigation of different viscosity models on pulsatile blood flow of thoracic aortic aneurysm (TAA) in a patient-specific model. Computer Methods in Biomechanics and Biomedical Engineering, 26(8), 986-998, (2022).
• [12] Bilgi, C. and Atalık, K. Numerical investigation of the effects of blood rheology and wall elasticity in abdominal aortic aneurysm under pulsatile flow conditions. Biorheology, 56(1), 51-71, (2019).
• [13] Panchal, P.M., Hathi, D.S., Shah, N.K. and Lakdawala, A.M. A numerical study on pulsatile non-Newtonian hemodynamics in double-fusiform abdominal aortic aneurysms. Physics of Fluids, 34(3), 031908, (2022).
• [14] Johnston, B.M., Johnston, P.R., Corney, S. and Kilpatrick, D. Non-Newtonian blood flow in human right coronary arteries: steady state simulations. Journal of Biomechanics, 37(5), 709-720, (2004).
• [15] Johnston, B.M., Johnston, P.R., Corney, S. and Kilpatrick, D. Non-Newtonian blood flow in human right coronary arteries: transient simulations. Journal of Biomechanics, 39(6), 1116-1128, (2006).
• [16] Kandangwa, P., Torii, R., Gatehouse, P.D., Sherwin, S.J. and Weinberg, P.D. Influence of right coronary artery motion, flow pulsatility and non-Newtonian rheology on wall shear stress metrics. Frontiers in Bioengineering and Biotechnology, 10, 962687, (2022).
• [17] Bernabeu, M.O., Nash, R.W., Groen, D., Carver, H.B., Hetherington, J., Krüger, T. and Coveney, P.V. Impact of blood rheology on wall shear stress in a model of the middle cerebral artery. Interface Focus, 3(2), 20120094, (2012).
• [18] Oliveira, I., Santos, G.B., Gasche, J.L., Militzer, J. and Baccin, C.E. Non-Newtonian blood modeling in intracranial aneurysm hemodynamics: impact on the WSS and OSI metrics for ruptured and unruptured cases. Journal of Biomechanical Engineering, (2021).
• [19] Singh, D. and Singh, S. Computational analysis of patient-specific pulsatile blood flow: The influence of non-Newtonian models on wall shear stress assessment. Physics of Fluids, 36(1), 013123, (2024).
• [20] Liu, S., Wang, S., Tian, H., Xue, J., Guo, Y., Yang, J. et al. Comparison of blood viscosity models in different degrees of carotid artery stenosis. PeerJ, 13, e19336, (2025).
• [21] Peng, C., Liu, J., He, W., Qin, W., Yuan, T., Kan, Y. et al. Numerical simulation in the abdominal aorta and the visceral arteries with or without stenosis based on 2D PCMRI. International Journal for Numerical Methods in Biomedical Engineering, 38(3), e3569, (2021).
• [22] Khanafer, K.M., Gadhoke, P., Berguer, R. and Bull, J.L. Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood. Biorheology, 43(5), 661-679, (2006).
• [23] Manolesos, M., Celik, Y., Ramsay, H., Karande, R., Wood, B., Dinwoodie, I. et al. Performance improvement of a Vestas V52 850kW wind turbine by retrofitting passive flow control devices. Journal of Physics: Conference Series, 2767(2), 022027, (2024).
• [24] Kaya, M.N., Satcunanathan, S., Meinke, M. and Schröder, W. Leading-edge noise mitigation on a rod–airfoil configuration using regular and irregular leading-edge serrations. Applied Sciences, 15(14), 7822, (2025).
• [25] Kaya, M.N. Numerical Investigation of Hydrofoil Cavitation Using OpenFOAM: Effect of Thickness and Camber. Black Sea Journal of Engineering and Science, 8(5), 1585-1594, (2025).
• [26] Kocaaslan, O., Güleren, K.M., Saracoğlu, B.H. and Yasa, T. Dönen patlama motorlarındaki oluşan dalga yapısının sayısal olarak incelenmesi. Journal of Thermal Science and Technology, 44(1), 33-45, (2024).
• [27] Carreau, P.J. Rheological equations from molecular network theories. Transactions of the Society of Rheology, 16(1), 99-127, (1972).
• [28] Cho, Y.I. and Kensey, K.R. Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows. Biorheology, 28(3-4), 241-262, (1991).
• [29] GitHub Repository, NextFOAM: BARAM. [https://github.com/nextfoam/baram], [Accessed: 14/09/2024]
• [30] Hariharan, P., Giarra, M., Reddy, V., Day, S.W., Manning, K.B., Deutsch, S. et al. Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations. Journal of Biomechanical Engineering, 133(4), 041002, (2011).
• [31] Trias, M., Arbona, A., Massó, J., Miñano, B. and Bona, C. FDA’s nozzle numerical simulation challenge: non-Newtonian fluid effects and blood damage. PLoS One, 9(3), e92638, (2014).
• [32] Krishnamoorthy, G. and Gholizadeh, N. Blood damage analysis within the FDA benchmark nozzle geometry at laminar conditions: Prediction sensitivities to software and non-Newtonian viscosity models. Symmetry, 16(9), 1165, (2024).
• [33] Fang, P., Wu, P., Feng, H. and Huang, H. Enhancing thrombosis prevention in medical devices: The role of turbulence in washout performance using FDA benchmark nozzle model. Computer Methods and Programs in Biomedicine, 261, 108647, (2025).
• [34] Malinauskas, R.A., Hariharan, P., Day, S.W., Herbertson, L.H., Buesen, M., Steinseifer, U. et al. FDA benchmark medical device flow models for CFD validation. In American Society for Artificial Internal Organs (ASAIO) Platinum 70th Anniversary Special Edition (pp. 317-327). Boca Raton, USA: CRC Press, (2024).