颅内动脉瘤两相血流动力学分析及PIV实验研究
|
摘要
为了探索血液动力学因素对颅内动脉瘤的影响,分析动脉瘤的形成、生长和破裂与血流动力学的关系,为治疗颅内动脉瘤提供技术参考和研究方法.数值模拟部分:依据个体化脑部动脉瘤患者的CT三维重建模型,从计算流体力学的角度出发,将血液假设为血浆和红细胞组成的两相流流体进行瞬态数值模拟.并且与经典牛顿流体与非牛顿流体进行对比.实验部分:利用与数值模拟相同的边界条件进行PIV流场可视化分析.得到一个心跳周期内的血液动力学参数分布情况和相同条件下实验结果.与牛顿流体和非牛顿流体相比,两相流在特征点的血液流速波动更为复杂;在壁面剪切力分布更加不"平滑",红细胞的粒子速率和沉积对颅内动脉瘤有着至关重要的影响,两相流更符合真实血液流动.实验结果证明血液流动会在动脉瘤内部形成低速旋涡区域,并且随着时间的变化旋涡中心不断转移,与数值模拟结果相吻合.
Objective:To explore the influence of hemodynamic factors on intracranial aneurysms,and to analyze the relationship between the formation,growth and rupture of aneurysms and hemodynaamics,and to provide technical references and research methods for the treatment of intracranial aneurysms.Methods:Numerical simulation:According to the CT three-dimensional reconstruction model of individual brain aneurysm patients,from the perspective of computational fluid dynamics,the blood was assumed to be a transient twophase fluid composed of plasma and red blood cells.And compared with classic Newtonian fluids and non-Newtonian fluids.Experimental section:The PIV flow field visualization analysis was performed using the same boundary conditions as the numerical simulation.RESULTS:A distribution of hemodynamic parameters during the heartbeat period and experimental results under the same conditions were obtained.Conclusion:Compared with Newtonian fluids and non-Newtonian fluids,the fluctuation of blood flow velocity at the characteristic points of the two-phase flow is more complicated;the shear stress distribution on the wall is less "smooth";the particle velocity and sedimentation of red blood cells have an effect on intracranial aneurysms.Of critical importance,the two-phase flow is more in line with the true blood flow.The experimental results show that the blood flow will form a low velocity vortex area inside the aneurysm,and the vortex center will continue to shift with time,which is consistent with the numerical simulation results.
Objective:To explore the influence of hemodynamic factors on intracranial aneurysms,and to analyze the relationship between the formation,growth and rupture of aneurysms and hemodynaamics,and to provide technical references and research methods for the treatment of intracranial aneurysms.Methods:Numerical simulation:According to the CT three-dimensional reconstruction model of individual brain aneurysm patients,from the perspective of computational fluid dynamics,the blood was assumed to be a transient twophase fluid composed of plasma and red blood cells.And compared with classic Newtonian fluids and non-Newtonian fluids.Experimental section:The PIV flow field visualization analysis was performed using the same boundary conditions as the numerical simulation.RESULTS:A distribution of hemodynamic parameters during the heartbeat period and experimental results under the same conditions were obtained.Conclusion:Compared with Newtonian fluids and non-Newtonian fluids,the fluctuation of blood flow velocity at the characteristic points of the two-phase flow is more complicated;the shear stress distribution on the wall is less "smooth";the particle velocity and sedimentation of red blood cells have an effect on intracranial aneurysms.Of critical importance,the two-phase flow is more in line with the true blood flow.The experimental results show that the blood flow will form a low velocity vortex area inside the aneurysm,and the vortex center will continue to shift with time,which is consistent with the numerical simulation results.
引文
[1]Nieuwkamp D J,SETA L E,ALGRA A,et al.Changes in case fatality of aneurysmal subarachnoid haemorrha,ge over time,according to age,sex,and region:a metaranalysis[J].Lancet Neurology,2009,8(7):598-599.
[2]Borghia A,Wooda N B,Mohiaddinb R H,et al.Fluid-solid interaction simulation of flow and stress pattern in thoracoabdominal aneurysms:A patient-specific study[J].J Fluid Structure,2008,24:270-280
[3]CONNOLLY E S Jr,RABINSTEIN A A,CARHUAPOMA J R,et al.Guidelines for the management of aneurysmal subarachnoid h emorrhage:A guideline for healthcare professionals from the American Heart Association/American Stroke Association[J].Stroke,2012,43(6):1711-1737.
[4]GONZALEZ C F,CHO Y I,MORET J,et al.Intracranial aneurysms:Flow analysis of their origin and progression[J].Am J Neuroradiol,1992,13(1):181-188.
[5]Cheng C,Tempel D,Van Haperen R,et al.Atherosclerotic lesion size and vulnerability are determined by patterns of fluid shear stress[J].Circulation,2006,113(22):2744-2753..
[6]刘莹,章德发,毕勇强等.主动脉弓及分支血管内非稳态血流分析[J].应用数学和力学,2015,36(4):432-439.Liu Y,Zhang DF,Bi YQ,et al.The analysis of unsteady blood flowin aortic bifurcation[J].Applied Mathematics and Mechanics,2015,36(4):432-439.
[7]Perktold K,Rappitsch G.Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model[J].J Biomech,1995,28(7):845-856.
[8]Liu Y,Lai Y,Nagaraj A,et al.Pulsatile flow simulation in arterial vascular segments with intravascular ultrasound images[J].Med Eng Phys,2001,23(8):583-595.
[9]Tang D,Yang C,Ku D.A 3D thin-wall model with fluid structure interactions for blood flow in carotid arteries with symmetric and asymmetric stenoses[J].Comput Struct,1999,72(1-3):357-377
[10]Cebral J R,Lohner R.From medical images.to anatomically accurate finite element grids[J].International Journal for Numerical Meth,ods in Engineering,2001,51(51):985-1008.
[11]Liu B.The influences of stenosis on the downstream flow pattern in curved arteries[J].Med Eng Phys,2007,29(8):868-876.
[12]Vimmr J.Numerical analysis of non-Newtonian blood flow and wall shear stress in realistic single,double an d triple aorto-coronary bypasses.[J].Int J Numer Method Biomed Eng,2013,29(10):1057-1081.
[13]王质刚.血液净化学[M].北京:科学技术出版社,2010.
[14]Wu W T,Aubry N,Massoudi M.On the coefficients of the interaction forces in a two-phase flow of a fluid infused with particles
[15]KrishnanB.Chandran,StanleyE.Rittgers,AjitP.Yoganathan.生物流体力学[M].机械工业出版社,2015.
[16]Davies P F.Hemodynamic shear stress and the endothelium in cardiovascular pathophysiology[J].Nat Clin Pract Card,2009,6(1):16-26.
[17]SFORZA D M,PUTMAN C M,SCRIVANO E,et al.Bloodflow characteristics in a terminal basilar tip aneurysm prior to its fatal rupture[J].Am J Neuroradiol,2010,31(6):1127-1131.
[18]JING L,FAN J,WANG Y,et al.Morphologic and hemodynamic analysis in the patients with multiple intracranial aneurysms:Ruptured versus unruptured[J].PLoS One,2015,10(7):e0132494.
[19]XIANG J,NATARAJAN S K,TREMME M,et al.Hemodynamic-morphologic discriminants for intracranial aneurysm rupture[J].Stroke,2011,42(1):144-152.
[20]CEBRAL J R,MUT F,WEIR J,et al.Quantitative characterization of the h emodynamic environment in ruptured and unruptured brain aneurysms[J].Am J Neuroradiol,2011,32(1):145-151.
[21]KENICHI K,TOMOAKI T.Treatment strategy and followup evaluation for an unruptured anterior communicating artery aneurysm associated with pseudo-occlusion of the internal carotid artery using computational fluid dynamics simulations[J].Turk Neurosurg,2014,24(1):111-116
[2]Borghia A,Wooda N B,Mohiaddinb R H,et al.Fluid-solid interaction simulation of flow and stress pattern in thoracoabdominal aneurysms:A patient-specific study[J].J Fluid Structure,2008,24:270-280
[3]CONNOLLY E S Jr,RABINSTEIN A A,CARHUAPOMA J R,et al.Guidelines for the management of aneurysmal subarachnoid h emorrhage:A guideline for healthcare professionals from the American Heart Association/American Stroke Association[J].Stroke,2012,43(6):1711-1737.
[4]GONZALEZ C F,CHO Y I,MORET J,et al.Intracranial aneurysms:Flow analysis of their origin and progression[J].Am J Neuroradiol,1992,13(1):181-188.
[5]Cheng C,Tempel D,Van Haperen R,et al.Atherosclerotic lesion size and vulnerability are determined by patterns of fluid shear stress[J].Circulation,2006,113(22):2744-2753..
[6]刘莹,章德发,毕勇强等.主动脉弓及分支血管内非稳态血流分析[J].应用数学和力学,2015,36(4):432-439.Liu Y,Zhang DF,Bi YQ,et al.The analysis of unsteady blood flowin aortic bifurcation[J].Applied Mathematics and Mechanics,2015,36(4):432-439.
[7]Perktold K,Rappitsch G.Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model[J].J Biomech,1995,28(7):845-856.
[8]Liu Y,Lai Y,Nagaraj A,et al.Pulsatile flow simulation in arterial vascular segments with intravascular ultrasound images[J].Med Eng Phys,2001,23(8):583-595.
[9]Tang D,Yang C,Ku D.A 3D thin-wall model with fluid structure interactions for blood flow in carotid arteries with symmetric and asymmetric stenoses[J].Comput Struct,1999,72(1-3):357-377
[10]Cebral J R,Lohner R.From medical images.to anatomically accurate finite element grids[J].International Journal for Numerical Meth,ods in Engineering,2001,51(51):985-1008.
[11]Liu B.The influences of stenosis on the downstream flow pattern in curved arteries[J].Med Eng Phys,2007,29(8):868-876.
[12]Vimmr J.Numerical analysis of non-Newtonian blood flow and wall shear stress in realistic single,double an d triple aorto-coronary bypasses.[J].Int J Numer Method Biomed Eng,2013,29(10):1057-1081.
[13]王质刚.血液净化学[M].北京:科学技术出版社,2010.
[14]Wu W T,Aubry N,Massoudi M.On the coefficients of the interaction forces in a two-phase flow of a fluid infused with particles
[15]KrishnanB.Chandran,StanleyE.Rittgers,AjitP.Yoganathan.生物流体力学[M].机械工业出版社,2015.
[16]Davies P F.Hemodynamic shear stress and the endothelium in cardiovascular pathophysiology[J].Nat Clin Pract Card,2009,6(1):16-26.
[17]SFORZA D M,PUTMAN C M,SCRIVANO E,et al.Bloodflow characteristics in a terminal basilar tip aneurysm prior to its fatal rupture[J].Am J Neuroradiol,2010,31(6):1127-1131.
[18]JING L,FAN J,WANG Y,et al.Morphologic and hemodynamic analysis in the patients with multiple intracranial aneurysms:Ruptured versus unruptured[J].PLoS One,2015,10(7):e0132494.
[19]XIANG J,NATARAJAN S K,TREMME M,et al.Hemodynamic-morphologic discriminants for intracranial aneurysm rupture[J].Stroke,2011,42(1):144-152.
[20]CEBRAL J R,MUT F,WEIR J,et al.Quantitative characterization of the h emodynamic environment in ruptured and unruptured brain aneurysms[J].Am J Neuroradiol,2011,32(1):145-151.
[21]KENICHI K,TOMOAKI T.Treatment strategy and followup evaluation for an unruptured anterior communicating artery aneurysm associated with pseudo-occlusion of the internal carotid artery using computational fluid dynamics simulations[J].Turk Neurosurg,2014,24(1):111-116