Numerical Study of Bifurcated Blood Flow in Three Different Blood Viscosity Models

doi:10.1007/s12204-022-2512-8

Abstract

Abstract: This study aimed to analyze the hemodynamic effects of bifurcated vessels using different bloodviscosity models. Three-dimensional models of bifurcated vessels in the popliteal artery were constructed basedon CT images, and hemodynamic parameters of the Newtonian, Casson, and two-phase models were calculatedby the computational fluid dynamics method. Blood flowed through the popliteal artery. Blood flow velocitychanged after the bifurcated vessel, with accelerated blood flow velocity in the anterior tibial artery. A lowvelocity vortex region with a region of low wall shear stress (WSS) was generated outside the bifurcated vessel.Local non-Newtonian importance factors of great than 1 (i.e., IL > 1) occurred during the cardiac cycle, andIL > 1.75 occurred at the beginning and end of the cycle. Compared with the Casson and two-phase models, theNewtonian model has a larger vortex region and lower WSS. Low-velocity vortex regions and low WSS regionsin the bifurcated vessels may contribute to the development of atherosclerosis. Blood exhibited non-Newtonianfluid properties in bifurcated vessels (IL > 1), and the effect of non-Newtonian properties was more pronouncedat the beginning and end of heartbeats (IL > 1.75). The Newtonian model predicts a higher risk of atherosclerosisformation and the effect of non-Newtonian properties of blood should be considered in hemodynamic studies.

Key words: hemodynamics, computational fluid dynamics, non-Newtonian fluid, non-Newtonian importance factor

摘要： 本研究旨在分析不同血液黏度模型分叉血管的血流动力学影响。基于CT图像构建腘动脉分叉血管三维模型，通过计算流体力学（CFD）方法计算牛顿模型、Casson模型和两相模型的血流动力学参数。血液流经腘动脉分叉血管后血液流速发生改变，胫前动脉血液流速加快；分叉血管外侧产生低速涡流区并伴随低壁面切应力(WSS)区域。在心动周期内局部非牛顿重要性系数IL > 1，心动初期和心动末期 IL > 1.75。与Casson模型和两相模型相比，牛顿模型的涡流区更大，WSS更低。分叉血管的低速涡流区和低壁面切应力区可能会导致动脉粥样硬化的产生。血液在分叉血管表现为非牛顿流体特性（IL > 1），在心动初期和心动末期非牛顿特性影响显著（IL > 1.75）。牛顿模型预测的动脉粥样硬化形成的风险更高；血流动力学研究中，应考虑血液非牛顿特性的影响。

关键词: 血流动力学，计算流体力学，非牛顿流体，非牛顿重要性系数

CLC Number:

Q 66
R 318

WU Hui(吴辉)， FU Rongchang* (富荣昌)， YANG Xiaoyu (杨晓玉)， LI Xianzheng (李现政)， WANG Zhaoyao (王召耀). Numerical Study of Bifurcated Blood Flow in Three Different Blood Viscosity Models[J]. J Shanghai Jiaotong Univ Sci, 2023, 28(4): 450-.

References

[1] ABUGATTAS C, AGUIRRE A, CASTILLO E, et al.Numerical study of bifurcation blood flows using threedifferent non-Newtonian constitutive models [J]. Applied Mathematical Modelling, 2020, 88: 529-549.
[2] LIU Y, SUN A. Influence of forward acceleration onhemodynamic characteristics of carotid arteries: A numerical simulation [J]. Journal of Beijing Universityof Aeronautics and Astronautics, 2019, 45(3): 493-498(in Chinese).
[3] LU Y, ZHANG S S. Rheological characteristics analysis of blood flow at the site of vascular lesion based onCFD [J]. Journal of Northeastern University (NaturalScience), 2020, 41(1): 148-152 (in Chinese).
[4] MORALES H G, LARRABIDE I, GEERS A J, et al.Newtonian and non-Newtonian blood flow in coiledcerebral aneurysms [J]. Journal of Biomechanics, 2013,46(13): 2158-2164.
[5] SUZUKI T, TAKAO H, SUZUKI T, et al. Variabilityof hemodynamic parameters using the common viscosity assumption in a computational fluid dynamicsanalysis of intracranial aneurysms [J]. Technology andHealth Care, 2017, 25(1): 37-47.
[6] ABBASIAN M, SHAMS M, VALIZADEH Z, et al. Effects of different non-Newtonian models on unsteadyblood flow hemodynamics in patient-specific arterialand Programs in Biomedicine, 2020, 186: 105185.
[7] OU C B, HUANG W, YUEN M M F, et al. Hemodynamic modeling of leukocyte and erythrocyte transport and interactions in intracranial aneurysms bya multiphase approach [J]. Journal of Biomechanics,2016, 49(14): 3476-3484.
[8] OSTROWSKI Z, MELKA B, ADAMCZYK W, et al.CFD analysis of multiphase blood flow within aortaand its thoracic branches of patient with coarctationof aorta using multiphase Euler - Euler approach [J].Journal of Physics: Conference Series, 2016, 745:032112.
[9] ARIBAS E, CELEBI M S. A thermal based RBC aggregation model for two-phase blood flow [J]. KoreaAustralia Rheology Journal, 2020, 32(2): 121-136.
[10] XIE H W, ZHANG Y. The effect of red blood cells onblood heat transfer [J]. International Journal of Heatand Mass Transfer, 2017, 113: 840-849.
[11] ZHANG W M, LIU J L, YAN Q, et al. Computationalhaemodynamic analysis of left pulmonary artery angulation effects on pulmonary blood flow [J]. InteractiveCardioVascular and Thoracic Surgery, 2016, 23(4):519-525.
[12] CONG M, ZHAO H, XU X, et al. Hemodynamic analysis on anomalous origin of the right coronary arteryfrom the left coronary artery sinus [J]. Journal of Medical Biomechanics, 2020, 35(3): 284-288 (in Chinese).
[13] QIAO Y Q, FU R C, ZHANG L H. Effects of sandtherapy on hemodynamics in the femoral artery bifurcation at different degrees of stenosis [J]. ChineseJournal of Tissue Engineering Research, 2020, 24(8):1218-1224 (in Chinese).
[14] STARY H C, CHANDLER A B, DINSMORE R E, etal. A definition of advanced types of atherosclerotic lesions and a histological classification of atherosclerosis.A report from the Committee on Vascular Lesions ofthe Council on Arteriosclerosis, American Heart Association [J]. Arteriosclerosis, Thrombosis, and VascularBiology, 1995, 15(9): 1512-1531.
[15] XIONG J, LUO W, LI D, et al. Hemodynamic studyon non-Newtonian properties of Fontan procedure [J].Journal of Medical Biomechanics, 2021, 36(6): 862-868 (in Chinese).
[16] MENG H, TUTINO V M, XIANG J, et al. HighWSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth,and rupture: Toward a unifying hypothesis [J]. AJNRAmerican Journal of Neuroradiology, 2014, 35(7):1254-1262.
[17] TANG D L, YANG C, MONDAL S, et al. A negative correlation between human carotid atherosclerotic plaque progression and plaque wall stress: In vivoMRI-based 2D/3D FSI models [J]. Journal of Biomechanics, 2008, 41(4): 727-736.
[18] GALLO D, BIJARI P B, MORBIDUCCI U, et al.Segment-specific associations between local haemodynamic and imaging markers of early atherosclerosisat the carotid artery: An in vivo human study [J].Journal of the Royal Society, Interface, 2018, 15(147):20180352.
[19] CHAPPELL D C, VARNER S E, NEREM R M, et al.Oscillatory shear stress stimulates adhesion moleculeexpression in cultured human endothelium [J]. Circulation Research, 1998, 82(5): 532-539.
[20] RIEDER M J, CARMONA R, KRIEGER J E, et al.Suppression of angiotensin-converting enzyme expression and activity by shear stress [J]. Circulation Research, 1997, 80(3): 312-319.
[21] BALLYK P D, STEINMAN D A, ETHIER C R. Simulation of non-Newtonian blood flow in an end-to-sideanastomosis [J]. Biorheology, 1994, 31(5): 565-586.
[22] FROLOV S V, SINDEEV S V, LIEPSCH D, et al.Newtonian and non-Newtonian blood flow at a 90?-bifurcation of the cerebral artery: A comparative studyof fluid viscosity models [J]. Journal of Mechanics inMedicine and Biology, 2018, 18(5): 1850043.
[23] AMIRI M H, KESHAVARZI A, KARIMIPOURA, et al. A 3-D numerical simulation of nonNewtonian blood flow through femoral artery bifurcation with a moderate arteriosclerosis: Investigating Newtonian/non-Newtonian flow and its effects onelastic vessel walls [J]. Heat and Mass Transfer, 2019,55(7): 2037-2047.
[24] MORAVIA A, SIMO?NS S, EL HAJEM M, et al. Invitro flow study in a compliant abdominal aorta phantom with a non-Newtonian blood-mimicking fluid [J].Journal of Biomechanics, 2022, 130: 110899.
[25] WEDDELL J C, KWACK J, IMOUKHUEDE P I, etal. Hemodynamic analysis in an idealized artery tree:Differences in wall shear stress between Newtonian andnon-Newtonian blood models [J]. PLoS ONE, 2015,10(4): e0124575.