×

一种泊松方程方法,用于规定任意横截面模型中计算流体动力学模拟的充分发展的非牛顿入口条件。 (英语) 兹比尔07478812

摘要:在复杂几何形状(如解剖血管模型)的内部流动的计算流体动力学(CFD)模拟中,规定入口边界条件具有挑战性。在缺乏患者特定的入口速度数据的情况下,长血管的一种常见方法是假设入口流量充分发展。然而,在不规则横截面的容器中,由于缺乏通用的封闭式分析解,规定充分发展的条件是复杂的。在本研究中,我们开发了一种简单的泊松方程方法,用于规定任意横截面CFD模型中牛顿流体或非牛顿流体充分发展的入口条件。我们首先推导了非牛顿流体充分发展流动的广义泊松方程,然后发展并验证了在任何平面边界域上数值计算解的方法。此外,我们对该方法进行了简单扩展,以规定非正交入口速度,该速度表示上游管道以非正交角度连接到CFD入口的充分流动。这可用于调查内部流动CFD模拟中常见的不确定性来源,这是由于缺乏有关进水口处准确流向的信息。与几种牛顿和非牛顿基准验证解决方案的比较表明,该方法非常准确。作为一个实际的演示案例,我们使用该方法为患者特定的下腔静脉(IVC)模型中的非牛顿血液流动在多个非圆形入口上规定充分发展的条件。最后,我们通过使用患者特定的下腔静脉模型进行敏感性研究,进一步证明了该方法的实用性,其中我们研究了入口流速流向对非牛顿下腔静脉血流动力学的影响。考虑到其简单性和计算效率,该方法被证明远远优于在如此复杂的几何形状中规定充分发展的入口条件的替代方法。为了便于采用我们的泊松方程方法,我们将OpenFOAM源代码和本研究中的相关测试用例作为开源软件分发。

MSC公司:

76倍 流体力学
92年XX月 生物学和其他自然科学

软件:

开放式泡沫
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 艾科克,K.I。;坎贝尔,R.L。;林奇,F.C。;曼宁,K.B。;Craven,B.A.,《血液流变学和患者解剖对下腔静脉血流动力学的重要性》,Ann.Biomed。工程,44,3568-3582(2016)
[2] 艾科克,K.I。;坎贝尔,R.L。;林奇,F.C。;曼宁,K.B。;Craven,B.A.,下腔静脉滤器在特定患者和理想下腔静脉几何形状中栓塞捕捉性能的计算预测,生物医学。模型。机械双醇。,16, 1957-1969 (2017)
[3] 艾科克,K.I。;坎贝尔,R.L。;曼宁,K.B。;萨斯特里,S.P。;肖恩茨,S.M。;林奇,F.C。;Craven,B.A.,《根据患者的具体情况预测下腔静脉滤器性能的计算方法》,J.Biomech。工程,136,第081003条pp.(2014)
[4] 艾科克,K.I。;哈里哈兰,P。;Craven,B.A.,使用喷墨3D打印制作的解剖血管模型中的粒子图像测速测量,实验流体,58,154(2017)
[5] 艾科克,K.I。;Rebelo,N。;Craven,B.A.,《制造解决方案的方法》,《商用有限元解算器中的弹性静态固体力学问题的代码验证》,计算。结构。,229,第106175条pp.(2020)
[6] 伯德·R·B。;R.C.阿姆斯特朗。;Hassager,O.,《聚合物液体动力学》。第1卷:流体力学(1987),John Wiley&Sons公司:John Willey&Sons,Inc.多伦多
[7] 伯德·R·B。;Stewart,W.E。;Lightfoot,E.N.,《运输现象》(2002),John Wiley&Sons,Inc.:纽约John Willey&Sons公司
[8] 波恩,D。;罗德斯,S。;Groenik,M。;拉斐,S。;Shahidzadeh-Bonn,N。;Coussot,P.,磁共振成像在流体力学中的一些应用:复杂流动和复杂流体,阿努。流体力学版次。,40, 209-233 (2008) ·Zbl 1136.76006号
[9] 负荷,R.L。;Faires,J.D.,《数值分析》(2011),Brooks/Cole:Brooks/Cole Boston,MA
[10] Capaldi,F.M.,《连续介质力学:结构和生物材料的本构建模》(2012),剑桥大学出版社:剑桥大学出版社,纽约·Zbl 1345.74001号
[11] 卡斯特罗,M。;普特曼,C.M。;《颅内动脉瘤的计算流体动力学建模:母动脉分割对动脉瘤内血流动力学的影响》,美国神经放射学杂志。,27, 1703-1709 (2006)
[12] Cebral,J.R。;Castro,医学硕士。;伯吉斯·J·E。;Pergolizzi,R.S。;谢里丹,M.J。;Putman,C.M.,《使用患者特定的计算血流动力学模型评估脑动脉瘤破裂风险的特征描述》,美国神经放射杂志。,26, 2550-2559 (2005)
[13] 程,C.P。;Herfkens,R.J。;Taylor,C.A.,《利用磁共振成像对在体休息和自行车运动期间下腔静脉血流动力学进行量化》,美国生理学杂志。心脏循环。生理学。,284,H1161-H1167(2003)
[14] Cho,Y.I。;Kensey,K.R.,血液的非牛顿粘度对患病动脉血管内流动的影响。第1部分:稳定流动,生物流变学,28,241-262(1991)
[15] 克雷文,B.A。;艾科克,K.I。;Herbertson,L.H。;Malinauskas,R.A.,一种基于CFD的克里金替代建模方法,用于预测血液控制医疗设备中特定设备的溶血幂律系数,Biomech。模型。机械双醇。,18, 1005-1030 (2019)
[16] 克雷文,B.A。;艾科克,K.I。;Manning,K.B.,患者平均下腔静脉中的稳定流动:第二部分-计算流体动力学验证,心血管。工程技术。,9, 654-673 (2018)
[17] 克雷文,B。;帕特森,E。;Settles,G.,《犬嗅觉的流体动力学:独特的鼻腔气流模式作为巨大畸形的解释》,J.R.Soc.Interface,7933-943(2010)
[18] 克雷文,B。;帕特森,E。;定居,G。;Lawson,M.,犬鼻腔气流高保真计算流体动力学模型的开发和验证,J.Biotech。工程,131,文章091002 pp.(2009)
[19] De Stefano,G。;Denaro,F.M。;Riccardi,G.,通过基于局部平均值的数值程序分析三维后向台阶不可压缩流,国际。J.数字。《液体方法》,28,1073-1091(1998)·Zbl 0931.76057号
[20] E.W.,Chaves,《连续介质力学笔记》(2013),施普林格科学与商业媒体:西班牙巴塞罗那施普林格科技与商业媒体·Zbl 1277.74001号
[21] Formaggia,L。;Vergara,C.,《流体动力学中一般缺陷边界条件的规定》,米兰J.数学。,80, 333-350 (2012) ·Zbl 1342.76040号
[22] Formaggia,L。;Vergara,C.,《PDE的缺陷边界条件及其在血流动力学中的应用》,(Pietro,D.A.D.;Ern,A.;Formaggia,L.,PDE的数值方法(2018年),施普林格国际出版公司:瑞士施普林格国际出版公司Cham),285-312·Zbl 1416.65444号
[23] 加拉赫,M.B。;艾科克,K.I。;克雷文,B.A。;Manning,K.B.,患者平均下腔静脉内的稳定流动第i部分:休息和运动条件下的粒子图像测速测量,心血管。工程技术。,9, 641-653 (2018)
[24] 戈麦斯,A。;Pushparajah,K。;辛普森,J.M。;吉斯,D。;谢夫特,T。;Penney,G.,从多个记录的回波多普勒视图进行三维速度重建的敏感性分析,医学图像分析。,17, 616-631 (2013)
[25] 戈麦斯,A。;de Vecchi,A。;Pushparajah,K。;辛普森,J。;吉斯,D。;谢夫特,T。;Penney,G.,彩色多普勒图像和心室壁运动的三维心室内血流图,(Mori,K.;Sakuma,I.;Sato,Y.;Barillot,C.;Navab,N.,《医学图像计算和计算机辅助干预》,医学图像计算与计算机辅助干预,MICCAI 2013(2013),柏林-海德堡:柏林-斯普林格,海德堡),476-483
[26] Hackbusch,W.,张量空间和数值张量微积分(2012),Springer-Verlag:Springer-Verlag纽约·Zbl 1244.65061号
[27] 哈里哈兰,P。;Giarra,M。;雷迪,V。;戴,S.W。;曼宁,K.B。;德国,S。;斯图尔特,S.F。;迈尔斯,M.R。;伯曼,M.R。;Burgreen,G.W。;帕特森,E.G。;Malinauskas,R.A.,FDA基准喷嘴模型的多实验室粒子图像测速分析,以支持计算流体动力学模拟的验证,J.Biomech。Eng.,133,文章041002 pp.(2011)
[28] Hund,S.J。;安塔基,J.F。;Massoudi,M.,《计算血液损伤的湍流应力表示法》,国际。工程科学杂志。,48, 1325-1331 (2010) ·Zbl 1231.76367号
[29] H.Jasak,将有限面积法移植到OpenFOAM Plus,载于:《第五届OpenFOAM用户会议论文集》,德国威斯巴登,2017年。
[30] 姜杰。;Strother,C.,《不同心率下颅内动脉瘤的计算流体动力学模拟:一项针对患者的研究》,J.Biomech。工程,131,第091001条pp.(2009)
[31] 约翰斯顿,B.M。;约翰斯顿,P.R。;科尼,S。;Kilpatrick,D.,《人体右冠状动脉中的非牛顿血流:稳态模拟》,J.Biomech。,37, 709-720 (2004)
[32] Joseph,A.A。;沃伊特·D。;Frahm,J.,《再次访问下腔静脉-呼吸动作的实时流MRI》,NMR Biomed。(2020)
[33] 昆杜,P。;科恩,I。;Dowling,D.,《流体力学》(2012),爱思唯尔学术出版社·Zbl 1302.76001号
[34] Kyriakou,A。;Neufeld,E。;Szczerba,D。;Kainz,W。;卢辛格,R。;科泽克,S。;McGregor,R。;Kuster,N.,《人类初级血管中磁流体动力学效应的特定患者模拟和测量》,Physiol。测量。,33, 117-130 (2012)
[35] 洛茨,J。;梅耶,C。;Leppert,A。;Galanski,M.,《用相控MR成像进行心血管流量测量:基本事实和实施》,《放射学》,22,651-671(2002)
[36] Malecha,Z。;Mirosław,Ł。;托姆扎克,T。;科扎,Z。;马蒂卡,M。;Tarnawski,W。;Szczerba,D.,使用OpenFOAM对腹主动脉中的3D血流进行基于GPU的模拟,Arch。机械。,63, 137-161 (2011) ·Zbl 1301.76086号
[37] Malinauskas,R.A。;哈里哈兰,P。;戴,S.W。;Herbertson,L.H。;Buesen,M。;斯坦西弗,美国。;艾科克,K.I。;好,不列颠哥伦比亚省。;德国,S。;曼宁,K.B。;Craven,B.A.,FDA CFD验证医疗器械基准流模型,ASAIO J.,63,150-160(2017)
[38] 马克尔,M。;Chan,F。;马里兰州艾利。;婚礼,K。;Draney,M。;埃尔金斯,C。;帕克,D。;Wicker,R。;泰勒,C。;Herfkens,R.,《时间分辨三维相控MRI》,J.Magn。Reson公司。成像,17499-506(2003)
[39] 马克尔,M。;Kilner,P.J。;Ebbers,T.,通过心血管磁共振绘制心脏和大血管的综合4d速度图,心血管学杂志。Magn.公司。资源。,13, 7 (2011)
[40] Milner,J.S。;摩尔,J.A。;Rutt,B.K。;Steinman,D.A.,《人类颈动脉分叉的血流动力学:利用正常受试者磁共振成像重建模型的计算研究》,J.Vasc。外科学,28,143-156(1998)
[41] 莫里森,F.A.,《理解流变学》(2001),牛津大学出版社:牛津大学出版社,纽约·Zbl 1012.76500号
[42] 莫伊尔,K.R。;安提瓜,L。;Steinman,D.A.,基于图像的颈动脉分叉CFD模型的入口条件:假设充分流动合理吗?,生物计量学杂志。工程,128,371-379(2006)
[43] Nithiarasu,P。;萨佐诺夫一世。;Yeo,S.Y.,《明日医学中的患者特异性建模》,241-280(2012年),德国柏林施普林格大学
[44] Rygg,A.D。;Van Valkenburgh,B。;克雷文,B.A.,《嗅觉对犬鼻腔气流和气味沉积的影响》,《化学》。Senses,42683-698(2017)
[45] 图科维奇。;Jasak,H.,表面张力主导界面流体流动的移动网格有限体积界面跟踪方法,计算与流体,55,70-84(2012)·Zbl 1291.76223号
[46] 王,C。;佩坎,K。;de Zélicourt,D。;霍纳,M。;A.帕里哈尔。;Kulkarni,A。;Yoganathan,A.P.,《解剖性全腔静脉肺连接中流动不稳定性的CFD建模进展》,《生物医学年鉴》。工程师,35,1840-1856(2007)
[47] White,F.M.,《粘性流体流动》(2006),McGraw-Hill:纽约州纽约市McGraw-Hill
[48] White,F.M.,《流体力学》(2011),The McGraw-Hill Companies,Inc。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。