×

ALE框架中可压缩流体的表面张力。 (英语) Zbl 07504716号

摘要:我们描述了一种求解具有毛细力的可压缩欧拉系统的任意拉格朗日-欧拉(ALE)方法。该算法分为两个步骤。首先,拉格朗日步骤基于以细胞为中心的方案[D.E.伯顿等,计算。流体83,33–47(2013;Zbl 1290.76095号);B.设计C.马泽兰,建筑。定额。机械。分析。178,第3期,327–372页(2005年;Zbl 1096.76046号),P.-H.梅尔等人,SIAM J.Sci。计算。291781-1824年第4期(2007年;兹比尔1251.76028)]. 将表面张力离散化,以便在离散水平上准确验证拉普拉斯定律。我们还提供了二阶空间扩展和低马赫数校正,这不会破坏该方案的良好平衡特性。对拉格朗日格式的几个问题进行了评估,特别是针对线性Richtmyer-Meshkov不稳定性,这是目标应用。第二步是重新分区和重新绘制,这得益于使用界面附近精确交点的扫描-区域方法。我们使用流体体积(VOF)方法跟踪界面。我们描述了混合细胞的处理,特别是热力学闭包和曲率计算。该新方案用于研究表面张力对非线性Richtmyer-Meshkov不稳定性的影响。

MSC公司:

76倍 流体力学
65-XX岁 数值分析
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 贝尔,M。;Nunziato,J.,《反应性颗粒材料中爆燃-爆轰转变(DDT)的两相混合物理论》,国际J.Multiph。流量,12,6,861-889(1986)·Zbl 0609.76114号
[2] Barlow,A.J。;Maire,P.-H。;Rider,W.J。;Rieben,R.N。;Shashkov,M.J.,《模拟高速可压缩多材料流动的任意拉格朗日-欧拉方法》,J.Compute。物理。,322, 603-665 (2016) ·Zbl 1351.76093号
[3] Barth,T。;Jespersen,D.,非结构网格上迎风格式的设计和应用,(第27届航空航天科学会议(1989)),366
[4] 贝茨,J.W。;施密特,A.J。;卡拉西克,M。;Zalesak,S.T.,使用FastRad3D代码对平面惯性约束聚变目标中烧蚀瑞利-泰勒不稳定性的数值模拟,Phys。Plasmas,23,12,第122701条,pp.(2016)
[5] Belytschko,T.等人。;刘伟凯。;莫兰,B。;Elkhodary,K.,《连续统和结构的非线性有限元》(2013),John Wiley&Sons·Zbl 1279.74002号
[6] Benson,D.J.,《拉格朗日和欧拉水文代码中的计算方法》,《计算》。方法应用。机械。工程,99,2,235-394(1992)·Zbl 0763.73052号
[7] 伯恩特,M。;布雷尔,J。;加莱拉,S。;Kucharik,M。;Maire,P.-H。;Shashkov,M.,多材料任意拉格朗日-欧拉方法的两步混合保守重映射,J.Compute。物理。,230, 17, 6664-6687 (2011) ·Zbl 1408.65077号
[8] Brackbill,J.U。;科特,D.B。;Zemach,C.,《模拟表面张力的连续体方法》,J.Compute。物理。,100, 2, 335-354 (1992) ·Zbl 0775.76110号
[9] 伯顿,D.E。;卡尼,T.C。;摩根,N.R。;Sambasivan,S.K。;Shashkov,M.J.,固体动力学的以细胞为中心的拉格朗日-戈杜诺夫类方法,计算。流体,83,33-47(2013)·Zbl 1290.76095号
[10] Carré,G。;德尔·皮诺,S。;Després,B。;Labourasse,E.,基于任意维一般非结构网格的以单元为中心的拉格朗日流体动力学格式,J.Compute。物理。,228, 14, 5160-5183 (2009) ·Zbl 1168.76029号
[11] Chang,C.-H。;邓,X。;Theofanous,T.,《界面不稳定性的直接数值模拟:一致、保守、全速、尖锐界面方法》,J.Compute。物理。,242, 946-990 (2013) ·Zbl 1299.76097号
[12] Chauveheid,D.,FVCF-ENIP方法框架下的表面张力新算法,Eur.J.Mech。B、 流体,50,175-186(2015)·Zbl 1408.76378号
[13] Chiapolino,A。;索雷尔,R。;Nkonga,B.,《在非结构化网格上锐化可压缩流体的扩散界面》,J.Compute。物理。,340, 389-417 (2017) ·Zbl 1376.76068号
[14] 克拉利斯,J.-M。;Jaouen,S。;Raviart,P.-A.,计算线性扰动平面对称气体动力学流动的拉格朗日坐标下的Godunov型方法,J.Compute。物理。,198, 1, 80-105 (2004) ·Zbl 1107.76367号
[15] Cole,R.H.,《水下爆炸》(1965年),多佛出版社
[16] Coquerelle,M。;Glockner,S.,《承受表面张力的两相流在水平集框架中的四阶准确曲率计算》,J.Compute。物理。,305838-876(2016)·Zbl 1349.76053号
[17] 德尔·皮诺,S。;拉博拉塞,E。;Morel,G.,两个可压缩流与摩擦耦合系统的渐近保持多维ALE方法,J.Compute。物理。,363, 268-301 (2016) ·Zbl 1392.76036号
[18] Després,B.,以单元为中心的拉格朗日GLACE格式在任意维一般网格上的弱一致性,计算。方法应用。机械。工程,199,41,2669-2679(2010)·Zbl 1231.76177号
[19] Després,B。;拉古提埃,F。;拉博拉塞,E。;Marmajou,I.,《非结构网格上多组分流动的反耗散传输方案》,《国际有限期刊》第7卷,第30-65页(2010年)
[20] Després,B。;Mazeran,C.,《二维拉格朗日气体动力学和拉格朗夫系统》,Arch。定额。机械。分析。,178,3327-372(2005年)·Zbl 1096.76046号
[21] Dirac,P.,《量子力学原理》。牛津大学出版社,第27期(1981年)·Zbl 0477.46058号
[22] 杜兰德,O。;Jaouen,S。;Soulard,L。;O.休泽。;Colombet,L.,在存在粘度和表面张力的情况下,使用原子和连续方法进行微喷射的比较模拟,J.Appl。物理。,第122、13条,第135107页(2017年)
[23] 杜兰德,O。;Soulard,L.,用混合方法对冲击加载铜射流破碎和喷出物产生的大规模分子动力学研究,J.Appl。物理。,第111、4条,第044901页(2012年)
[24] 费希特,S。;蒙兹,C.-D。;罗德,C。;Zeiler,C.,具有相变和表面张力的可压缩液体蒸汽流动的近似黎曼解算器,计算。流体,169169-185(2018)·Zbl 1410.76171号
[25] Fedkiw,R.P。;Aslam,T。;梅里曼,B。;Osher,S.,《多材料流动界面的非振荡欧拉方法》(鬼流体方法),J.Compute。物理。,152, 2, 457-492 (1999) ·Zbl 0957.76052号
[26] 保险丝,D。;Popinet,S.,《模拟表面张力下气泡动力学问题的全马赫方法》,J.Compute。物理。,374, 752-768 (2018) ·Zbl 1416.76310号
[27] Fyfe,D.E。;奥兰,E.S。;Fritts,M.J.,三角网格上拉格朗日流体动力学的表面张力和粘度,J.Compute。物理。,76, 2, 349-384 (1988) ·兹伯利0639.76043
[28] 加莱拉,S。;Maire,P.-H。;Breil,J.,使用VOF界面重建的二维非结构化以细胞为中心的多材料ALE方案,J.Compute。物理。,229, 16, 5755-5787 (2010) ·Zbl 1346.76105号
[29] Garrick,D.P。;Owkes,M。;Regele,J.D.,基于有限体积HLLC的表面张力可压缩界面流动方案,J.Compute。物理。,339, 46-67 (2017) ·Zbl 1375.76099号
[30] Ghidaglia,J-M.,《毛细管力:体积公式》,《欧洲力学杂志》。B、 流体,59,86-89(2016)·Zbl 1408.76381号
[31] Goldman,R.,隐式曲线和曲面的曲率公式,计算。辅助Geom。设计。,22、7、632-658(2005),几何建模和微分几何·Zbl 1084.53004号
[32] Hoch,P.,《求解可压缩流动的任意拉格朗日-欧拉策略》(2009),CEA,技术报告HAL:HAL-00366858
[33] Hoch,P。;Labourasse,E.,基于局部凸壳保持的以细胞为中心的ALE方案的框架不变量和最大值原理的二阶扩展,国际期刊Numer。《液体方法》,761043-1063(2014)
[34] 伊藤,K。;库努吉,T。;Ohno,S。;卡米德,H。;Ohshima,H.,非结构化网格上界面法线和曲率的高精度计算方法,J.Compute。物理。,273, 38-53 (2014) ·兹比尔1351.76287
[35] Kataoka,I.,两相流的局部瞬时公式,国际J.Multiph。流量,12,5,745-758(1986)·Zbl 0613.76114号
[36] Kucharik,M。;Shashkov,M.,多材料任意拉格朗日-欧拉方法的一步混合重映射算法,J.Compute。物理。,231, 7, 2851-2864 (2012) ·Zbl 1323.74108号
[37] Labourasse,E.,多维有限体积激波捕获方案的低马赫数修正及其在拉格朗日框架中的应用,计算。流体,179372-393(2019)·Zbl 1411.76090号
[38] 拉法里,B。;纳尔多,C。;斯卡多维利,R。;Zaleski,S。;Zanetti,G.,《利用SURFER,J.Compute对多相流中的合并和破碎进行建模》。物理。,113, 1, 134-147 (1994) ·兹比尔0809.76064
[39] Lee,J。;林,K-C。;Eklund,D.,《高速推进系统燃油喷射的挑战》,AIAA J.,53,6,1405-1423(2015)
[40] Legendre,D。;Magnaudet,J.,粘性线性剪切流中球形气泡的升力,J.流体力学。,368, 81-126 (1998) ·Zbl 0927.76020号
[41] Legendre,D。;Magnaudet,J。;Mougin,G.,粘性液体中并排上升的两个球形气泡之间的流体动力学相互作用,J.流体力学。,497, 133-166 (2003) ·Zbl 1060.76026号
[42] 雷,X。;杜,Z。;Li,J.,用基于GRP的能量分裂方法模拟可压缩多流体流动,计算。流体,181416-428(2019)·Zbl 1410.76244号
[43] 卢特瓦克,G。;Falcovitz,J.,《向量的斜率限制:一种新的向量限制算法》,Int.J.Numer。《液体方法》,65,11-12,1365-1375(2011)·Zbl 1453.76100号
[44] Maire,P.-H.,非结构网格上二维可压缩流体流动的高阶以单元为中心的拉格朗日格式,J.Compute。物理。,228, 7, 2391-2425 (2009) ·Zbl 1156.76434号
[45] Maire,P.-H.,《聚变约束惯性的模型贡献》(2011),波尔多第一大学,博士论文
[46] Maire,P.-H。;Abgrall,R。;布雷尔,J。;Ovadia,J.,《二维可压缩流动问题的以单元为中心的拉格朗日格式》,SIAM J.Sci。计算。,29, 4, 1781-1824 (2007) ·Zbl 1251.76028号
[47] Mazeran,C.,《结构数学与流体力学拉格朗日二维近似》(2007年11月),波尔多大学I期博士论文
[48] Meshkov,E.,由冲击波加速的两种气体界面的不稳定性,流体动力学。,4, 5, 101-104 (1969)
[49] Morgan,N.,《水平集方法的一种新的液-气相变技术》(2005),佐治亚理工学院,博士论文
[50] Owkes,M。;考布尔,E。;塞内卡尔,J。;Currie,R.,动态气液界面预测模拟曲率评估尺度的重要性,J.Compute。物理。,365, 37-55 (2018) ·Zbl 1395.76101号
[51] 佩里高德,G。;Saurel,R.,《具有毛细管效应的可压缩流动模型》,J.Compute。物理。,209, 1, 139-178 (2005) ·Zbl 1329.76301号
[52] Popinet,S.,《表面张力驱动界面流动的精确自适应解算器》,J.Compute。物理。,228, 16, 5838-5866 (2009) ·Zbl 1280.76020号
[53] Popinet,S.,表面张力的数值模型,年。Rev.流体机械。,50, 49-75 (2018) ·Zbl 1384.76016号
[54] 波皮内特,S。;Zaleski,S.,《准确表示表面张力的前向跟踪算法》,国际期刊编号。液体方法,30,6,775-793(1999)·Zbl 0940.76047号
[55] 雷纳迪,Y。;Prost,M.Renardy,《流体体积法表面张力的抛物线重建》,J.Compute。物理。,183, 2, 400-421 (2002) ·Zbl 1057.76569号
[56] Richtmyer,R.D.,可压缩流体冲击加速度中的泰勒不稳定性,Commun。纯应用程序。数学。,13, 2, 297-319 (1960)
[57] 罗德,C。;Zeiler,C.,具有相变和表面张力的可压缩两相流的松弛黎曼解算器,Appl。数字。数学。,95, 267-279 (2015) ·Zbl 1320.76095号
[58] 索雷尔,R。;Gavrilyuk,S。;Renaud,F.,《具有内部自由度的多相模型:冲击气泡相互作用的应用》,《流体力学杂志》。,495, 283-321 (2003) ·Zbl 1080.76062号
[59] 施密德迈尔,K。;佩蒂帕斯,F。;丹尼尔,E。;Favrie,N。;Gavrilyuk,S.,《毛细管效应可压缩流动的模型和数值方法》,J.Compute。物理。,334, 468-496 (2017) ·Zbl 1375.76198号
[60] 施瓦茨,L。;数学研究所(斯特拉斯堡),《发行理论》,第2卷(1957年),赫尔曼:赫尔曼巴黎·兹伯利0078.11003
[61] Shin,S。;Juric,D.,使用水平轮廓重建方法模拟三维多相流,用于无连接的前方跟踪,J.Compute。物理。,180, 2, 427-470 (2002) ·Zbl 1143.76595号
[62] 萨斯曼,M。;Puckett,E.G.,用于计算三维和轴对称不可压缩两相流的耦合水平集和流体体积法,J.Compute。物理。,162, 2, 301-337 (2000) ·Zbl 0977.76071号
[63] 萨斯曼,M。;斯梅雷卡,P。;Osher,S.,《计算不可压缩两相流解的水平集方法》,J.Compute。物理。,114, 1, 146-159 (1994) ·兹比尔0808.76077
[64] Theofanous,T.,《牛顿和粘弹性液体的空气裂解》,年。Rev.流体机械。,43, 661-690 (2011) ·Zbl 1299.76217号
[65] Tryggvason,G。;邦纳,B。;Esmaeeli,A。;Juric,D。;Al-Rawahi,N。;陶伯,W。;Han,J。;Nas,S。;Jan,Y.-J.,多相流计算的前置跟踪方法,J.Compute。物理。,169, 2, 708-759 (2001) ·Zbl 1047.76574号
[66] S.O.Unverdi。;Tryggvason,G.,《粘性、不可压缩、多流体流动的前跟踪方法》,J.Compute。物理。,100, 1, 25-37 (1992) ·Zbl 0758.76047号
[67] Yang,Y。;张,Q。;夏普,D.H.,《Richtmyer-Meshkov不稳定性的小振幅理论》,物理学。流体,6,5,1856-1873(1994)·Zbl 0831.76040号
[68] Youngs,D.,具有大流体畸变的时间依赖性多材料流动,(流体动力学数值方法(1982))·Zbl 0537.76071号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。