×

用连续强迫浸没边界法施加自由滑移条件。 (英语) Zbl 1351.76111号

摘要:净化流体中球形和椭球形气泡的数值模拟需要在气泡表面施加自由滑移边界条件。本文描述了一种在浸没边界法中实现自由滑移边界条件的数值方法。与多相流的其他数值方法相比,实现起来并不简单。原因是浸没边界法将液相和气相视为密度和粘度恒定的场,在气泡内有一个虚拟流体。离散相的运动通过求解其每个元素的动量平衡来明确计算,并通过Navier-Stokes方程中的附加源项与连续相耦合。本文首先说明特别的方法不成功。在此基础上,提出了在气泡表面施加适当的直接力的新方法。一个中心发现是,在网格步长和气泡直径之间的共同比率下,需要考虑曲率项才能获得令人满意的结果。新方法首先针对球面对象开发,然后扩展到一般曲面界面。这是通过引入一个局部坐标系来实现的,该坐标系借助于该点处曲面的两个主曲率来近似拉格朗日标记附近的曲面。然后,对具有或不具有规定恒定角速度的球形和椭球形物体验证了该数值方案。结果表明,该方法与无滑移边界方法具有相似的收敛特性。将计算结果与绕球体蠕动流动的解析解和在贴体网格上获得的数值参考数据进行了比较。数值试验证实了该方法的优良性能。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
6500万08 含偏微分方程初值和初边值问题的有限体积法
76T10型 液气两相流,气泡流
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Aland,S。;施瓦兹,S。;弗罗里奇,J。;Voigt,A.,《液态金属中气泡的建模和数值近似》,《欧洲物理》。J.规格顶部。,220, 185-194 (2013)
[2] 艾伦,M.P。;Tildesley,D.J.,《液体的计算机模拟》(1989),牛津大学出版社·Zbl 0703.68099号
[3] 阿尔维斯,S.S。;Orvalho,S.P。;Vasconcelos,J.M.T.,《气泡污染对上升速度和传质的影响》,《化学》。工程科学。,60, 1-9 (2005)
[4] 安德森,E。;Bai,Z。;比肖夫,C。;南卡罗来纳州布莱克福德。;德米尔,J。;Dongarra,J。;杜克罗兹,J。;格林鲍姆,A。;Hammarling,S。;麦肯尼,A。;Sorensen,D.,《LAPACK用户指南》(1987),工业和应用数学学会·Zbl 0934.65030号
[5] Balaras,E.,在大规模模拟中使用固定笛卡尔网格上的外力场建模复杂边界,计算。流体,33,375-404(2004)·Zbl 1088.76018号
[6] 博拉兹贾尼,I。;Ge,L。;Sotiropulos,F.,模拟流体-结构与复杂三维刚体相互作用的Curvililinear浸没边界方法,J.Comput。物理。,227, 7587-7620 (2008) ·兹比尔1213.76129
[7] Chan-Braun,C。;García-Villalba,M。;Uhlmann,M.,《过渡粗糙明渠流中作用于颗粒的力和力矩》,J.流体力学。,684, 441 (2011) ·Zbl 1241.76393号
[8] 克利夫特,R。;格雷斯,J.R。;Weber,M.E.,《气泡、液滴和粒子》(1978),多佛出版社
[9] Crowe,C.,《多相流手册》(2006),CRC出版社·Zbl 1098.76001号
[10] Cuenot,B。;Magnaudet,J。;Spennato,B.,《微溶性表面活性剂对球形气泡周围流动的影响》,J.流体力学。,339, 25-35 (1997) ·Zbl 0892.76011号
[11] Diebel,J.,《表示姿态:欧拉角、单位四元数和旋转向量》,《矩阵》,58,15-16(2006)
[12] Duineveld,P.C.,《高雷诺数纯水中气泡的上升速度和形状》,《流体力学杂志》。,292, 325-332 (1995)
[13] 贝尔·菲拉(Bel Fdhila),R。;Duineveld,P.C.,《表面活性剂在高雷诺数和Peclet数下对球形气泡上升的影响》,Phys。流体,8310(1996)·Zbl 1027.76673号
[14] Ferzinger,J.H。;Peric,M.,流体动力学计算方法(2002),Springer-Verlag·Zbl 0998.76001号
[15] 弗朗索瓦,M。;E.乌兹戈伦。;J·杰克逊。;Shyy,W.,《多相流浸没边界技术的多重网格计算》,《国际数值杂志》。热流体流动方法,14,98-115(2004)·Zbl 1064.76081号
[16] 弗罗里奇,J。;施瓦兹,S。;海特卡姆,S。;桑塔雷利,C。;张,C。;沃格特,T。;博登,S。;Andruszkiewicz,A。;埃克特,K。;Odenbach,S.,磁场对液态金属中气泡行为的影响,《欧洲物理学》。J.规格顶部。,220, 167-183 (2013)
[17] Hadamard,J.S.,Mouvement permanent lend d’une sphere liquide et visqueuse dans un liquide visqueux,C.R.学院。科学。,1521735年(1911年)
[18] 哈佩尔,J。;Brenner,H.,《低雷诺数流体动力学》(1983),施普林格:荷兰施普林格
[19] 哈洛,F.H。;Welch,J.E.,含自由表面流体随时间变化的粘性不可压缩流动的数值计算,物理。流体,82182-2189(1965)·Zbl 1180.76043号
[20] 海特卡姆,S。;西德伦克汉。;Fröhlich,J.,《紧密包装球体:机械稳定性对封闭球体结构的影响》,Phys。修订稿。,108, 148302 (2012)
[21] Hua,J。;Stene,J.F。;Lin,P.,使用前沿跟踪方法对粘性液体中三维气泡上升的数值模拟,J.Compute。物理。,227, 3358-3382 (2008) ·Zbl 1329.76262号
[22] Kajishima,T。;Takiguchi,S.,粒子团簇和粒子诱导湍流之间的相互作用,《国际热流杂志》,23,639-646(2002)
[23] Kajishima,T。;高口,S。;滨崎步,H。;Miyake,Y.,由于旋涡脱落,垂直平面通道中颗粒流的湍流结构,JSME Int.J.Ser。B、 44526-535(2001)
[24] Kempe,T.,《带碰撞的粒子流界面再求解模拟的数值方法》(2013),德累斯顿大学流体力学研究所,德国TUD出版社Verlag der Wissenschaften GmbH,博士论文
[25] Kempe,T。;Fröhlich,J.,粘性流体中球形粒子界面重解模拟的碰撞建模,J.流体力学。,709, 445-489 (2012) ·Zbl 1275.76206号
[26] Kempe,T。;Fröhlich,J.,一种改进的浸没边界法,用于模拟颗粒载流,J.Compute。物理。,2313663-3684(2012年)·Zbl 1402.76143号
[27] Kempe,T。;沃温克尔,B。;Fröhlich,J.,《关于明渠水流泥沙输移界面再溶解模拟中碰撞建模的相关性》,国际J.Multiph。流量,58214-235(2014)
[28] Kidaneariam,A。;Chan-Braun,C。;Doychev,T。;Uhlmann,M.,有限尺寸水平明渠流的直接数值模拟,低固体体积分数的重颗粒,New J.Phys。,2015年5月15日(2013年)
[29] Kim,D。;Choi,H.,任意运动物体周围流动的浸没边界法,J.Compute。物理。,212662-680(2006年)·Zbl 1161.76520号
[30] Kim,J。;Kim,D。;Choi,H.,《模拟复杂几何形状流动的浸没边界有限体积法》,J.Compute。物理。,171, 132-150 (2001) ·Zbl 1057.76039号
[31] Kim,S。;Karrila,S.J.,《微流体动力学》(2005),多佛出版社
[32] Kühnel,W。;亨特,B.,《微分几何:曲线-曲面-流形》(2005),美国数学学会·Zbl 1082.53001号
[33] Lamb,H.,《流体动力学》(1993),多佛出版社·Zbl 0828.01012号
[34] 李,Z。;Ito,K.,《浸没界面法:涉及界面和不规则域的偏微分方程的数值解》(2006),工业数学学会·Zbl 1122.65096号
[35] Loth,E.,《分散颗粒、液滴和气泡运动的数值方法》,Prog。能源燃烧。科学。,26, 3, 161-223 (2000)
[36] Magnaudet,J。;Eames,I.,《非均匀流动中高雷诺数气泡的运动》,年。流体力学版次。,32, 659-708 (2000) ·Zbl 0989.76082号
[37] 马朱姆达尔,S。;艾卡里诺,G。;Durbin,P.,《具有自适应结构边界非协调网格的RANS解算器》,353-366(2001),斯坦福大学湍流研究中心,年度研究简报
[38] Malysa,K。;克拉索夫斯卡,M。;Krzan,M.,表面活性物质对气泡运动和与各种界面碰撞的影响,高级胶体界面科学。,114, 205-225 (2005)
[39] McLaughlin,J.B.,水中气泡运动的数值模拟,J.胶体界面科学。,184, 614-625 (1996)
[40] 梅尔,M。;Yadigaroglu,G。;Smith,B.L.,《PLIC-VOF方法中包含表面张力的新技术》,《欧洲力学杂志》。B流体,21,61-73(2002)·Zbl 1064.76084号
[41] Michaelides,E.E.,《粒子、气泡和液滴:它们的运动、热量和质量传递》(2006),世界科学出版社
[42] 米塔尔·R。;Dong,H。;博兹库塔斯,M。;Najjar,F.M。;瓦尔加斯,A。;von Loebbecke,A.,《复杂边界不可压缩流动的通用锐界面浸没边界法》,J.Compute。物理。,227, 4825-4852 (2008) ·Zbl 1388.76263号
[43] 米塔尔,R。;Iaccarino,G.,《浸入式边界方法》,Annu。流体力学版次。,37, 239-261 (2005) ·Zbl 1117.76049号
[44] Mohd-Yusof,J.,《复杂几何形状流动模拟的浸没边界/B样条组合方法》,317-327(1997),美国宇航局埃姆斯/斯坦福大学湍流研究中心,年度研究简报
[45] Osher,S。;Fedkiw,R.P.,《水平集方法:概述和一些最新结果》,J.Compute。物理。,169, 463-502 (2001) ·Zbl 0988.65093号
[46] Peskin,Ch.S.,心脏血流的数值分析,J.Compute。物理。,25, 220-252 (1977) ·Zbl 0403.76100号
[47] Peskin,Ch.S.,浸没边界法,数值学报。,11, 479-517 (2003) ·Zbl 1123.74309号
[48] 彼得斯,F。;Els,C.,《自来水中缓慢和快速气泡的实验研究》,《化学》。工程科学。,82, 194-199 (2012)
[49] Pöschel,T。;施瓦格,T.,《计算颗粒动力学:模型和算法》(2005),施普林格出版社
[50] Prosperetti,A。;Tryggvason,G.,《多相流计算方法》(2011),剑桥大学出版社
[51] 雷德菲尔德,J.A。;Houghton,G.,雷诺数为0.02-5000时单个气泡的传质和阻力系数,化学。工程科学。,20, 131-139 (1965)
[52] 罗黑尔,I。;刘永明。;迪恩,N.G。;造渣机,H.M。;巴尔图森,M.W。;Van Sint Annaland,M。;Kuipers,J.A.M.,《关于中等和高雷诺数下气泡群中气泡的阻力》,化学。工程科学。,66, 3204-3211 (2011)
[53] 罗马,A.M。;Peskin,C.S。;Berger,M.J.,浸没边界法的自适应版本,J.Compute。物理。,153, 509-534 (1999) ·Zbl 0953.76069号
[54] 罗曼,F。;那不勒斯,E。;米利奇,B。;Armenio,V.,一种改进的曲线网格浸没边界法,计算。流体,38,8,1510-1527(2009)·Zbl 1242.76210号
[55] Rybczynski,W.,将在公牛的einem zähen Medium中为Bewegung einer flüssigen Kugel提供支持。阿卡德。科学。克拉科夫。,40-46(1911),(德语)
[56] 桑塔雷利,C。;Fröhlich,J.,关于泡沫群中的配对相关函数,Kerntechnik,150-51(2013)
[57] 桑塔雷利,C。;施瓦兹,S。;Fröhlich,J.,垂直湍流道中轻椭球粒子的数值模拟,(第八届ICMF多相流国际会议,第八届国际ICMF多相流国际大会,韩国济州(2013))
[58] 斯卡多维利,R。;Zaleski,S.,自由表面和界面流动的直接数值模拟,年度。流体力学版次。,31, 567-603 (1999)
[59] Schwarz,S.,《磁流体力学流中粒子和气泡的浸没边界法》(2014年),德累斯顿大学流体力学研究所,德国TUDpress Verlag der Wissenschaften GmbH,博士论文
[60] 施瓦兹,S。;Fröhlich,J.,通过浸没边界法中的分析定义形状表示可变形气泡,(Simos,T.E.;Psihoyios,G.;Tsitouras,C.;Anastasi,Z.,《ICNAAM会议录》,ICNAAM,美国物理学会会议系列,第1479卷(2012)),104-108
[61] 施瓦兹,S。;Fröhlich,J.,《暴露于磁场中的高纵横比容器中气泡链的模拟》,《欧洲物理学》。J.规格顶部。,220, 195-205 (2013)
[62] 施瓦兹,S。;Fröhlich,J.,在纵向磁场作用下液态金属中单个气泡运动的数值研究,国际期刊Multiph。流量,62134-151(2014)
[63] 施瓦兹,S。;Kempe,T。;Fröhlich,J.,用浸没边界法计算粘性流中轻粒子运动的时间离散方案,J.compute。物理。,281, 591-613 (2015) ·兹比尔1351.76223
[64] 斯派克,J.H。;Aksel,N.,流体力学(2007),Springer·Zbl 1140.76001号
[65] Sridhar,G。;Katz,J.,涡旋夹带的微观气泡上的阻力和升力,Phys。流体,7389(1995)
[66] 萨斯曼,M。;法特米,E。;斯梅雷卡,P。;Osher,S.,不可压缩两相流的一种改进的水平集方法,计算。流体,27663-680(1998)·Zbl 0967.76078号
[67] 唐,Y。;Kriebitzsch,S.H.L。;彼得斯,E.A.J.F。;van der Hoef,医学硕士。;Kuipers,J.A.M.,《直接数值模拟高精度结果的方法:中等雷诺数下稠密气固流动中的阻力》,国际J.Multiph。流量,6273-86(2014)
[68] Tomiyama,A。;加藤岛,I。;尊,I。;Sakaguchi,T.,《正常重力和微重力条件下单个气泡的阻力系数》,JSME Int.J.Ser。B、 41472-479(1998)
[69] Tryggvason,G。;邦纳,B。;Esmaeeli,A。;Juric,D。;Al-Rawahi,N。;Tauber,W。;Han,J。;Nas,S。;Jan,Y.-J.,多相流计算的前置跟踪方法,J.Compute。物理。,169708-759(2001年)·Zbl 1047.76574号
[70] Tryggvason,G。;斯卡多维利,R。;Zaleski,S.,《气液多相流的直接数值模拟》(2011),剑桥大学出版社·Zbl 1226.76001号
[71] Uhlmann,M.,《模拟颗粒流的直接强迫浸没边界法》,J.Compute。物理。,209, 448-476 (2005) ·Zbl 1138.76398号
[72] Uhlmann,M.,湍流状态下垂直颗粒通道流动的界面重解析直接数值模拟,Phys。流体,20553305(2008)·Zbl 1182.76785号
[73] S.O.Unverdi。;Tryggvason,G.,《粘性、不可压缩、多流体流动的前跟踪方法》,J.Compute。物理。,100, 25-37 (1992) ·Zbl 0758.76047号
[74] 沃温克尔,B。;Kempe,T。;Fröhlich,J.,《泥沙输移的颗粒再溶解模拟》,(第八届多相流国际会议,第八届国际多相流会议,韩国济州(2013))
[75] 沃温克尔,B。;肯普,T。;Fröhlich,J.,《带全分辨率移动床的明渠紊流中的流体-颗粒相互作用》,《高级水资源》。(2014)
[76] 沃温克尔,B。;Kempe,T。;弗罗里奇,J。;Nikora,V.I.,明渠水流泥沙输移的数值模拟,(Murillo,R.,《河流流动》(2012)),507-514
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。