×

两相不可压缩流动的流体体积(VOF)界面锐化方法。 (英语) 兹比尔1390.76498

摘要:针对两相不可压缩流动,在通用曲线网格上发展了一种三维流体体积(VOF)界面锐化方法。在该方法中,针对双流体界面的平流,建立了VOF离散格式。为了保持界面清晰度,在每个平流时间步长后通过求解界面锐化方程进行处理,从而减少离散化方案求解中的数值扩散误差。为了证明平流格式的准确性和能力,进行了几次数值实验,包括三次纯平流基准试验。结果表明,该方法能够可靠有效地实现尖锐的界面,并获得了合理的质量守恒。对于粘性不可压缩流动的流场,采用双重时间预处理方法求解Navier-Stokes方程。采用全隐式方法和高效的上下对称高斯-赛德尔(LU-SGS)算法,将双时间步进技术作为子迭代方案,以及时推进求解。为了验证所提出的计算不可压缩自由表面流的方法,模拟了干式水平床溃坝水流和一个自由度半球的进水情况。本文给出了预测结果与现有实验数据的比较。

MSC公司:

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

参考文献:

[1] Gingold,R.A。;Monaghan,J.J.,《平滑粒子流体动力学:非球形恒星的理论和应用》,Mon Not R Astron Soc,181,375-389,(1977)·Zbl 0421.76032号
[2] 桑德斯,K。;普拉卡什,M。;克利里,P.W。;Cordell,M.,《光滑粒子流体动力学在闸门溢洪道水流建模中的应用》,应用数学模型,38,4308-4322,(2014)·Zbl 1428.76159号
[3] Oñate,E。;Idelsohn,S.R。;Celigueta,硕士。;Rossi,R.,《自由表面流中流体-多体相互作用和河床侵蚀分析的粒子有限元方法进展》,计算方法应用机械工程,1971777-1800,(2008)·Zbl 1194.74460号
[4] 奥尔森,E。;Kreiss,G.,两相流的保守水平集方法,《计算物理杂志》,210,225-246,(2005)·Zbl 1154.76368号
[5] 萨斯曼,M。;斯梅雷卡,P。;Osher,S.,《计算不可压缩两相流解的水平集方法》,《计算物理杂志》,114146-159,(1994)·Zbl 0808.76077号
[6] Osher,S。;Fedkiw,R.P.,水平集方法:概述和一些最新结果,《计算物理杂志》,169463-502,(2001)·兹伯利0988.65093
[7] 希特,C.W。;Nichols,B.D.,自由边界动力学的流体体积(VOF)方法,计算物理杂志,39,201-225,(1981)·Zbl 0462.76020号
[8] Noh W.F.、Woodward P.SLIC(简单线路接口计算)。收件人:A.Vooren,P.Zandbergen,编辑。第五届流体动力学数值方法国际会议论文集,1976年6月28日至7月2日,第二十大学,恩斯切德:施普林格-柏林-海德堡;1976年,第330-40页。;Noh W.F.、Woodward P.SLIC(简单线路接口计算)。收件人:A.Vooren,P.Zandbergen,编辑。第五届流体动力学数值方法国际会议论文集,1976年6月28日至7月2日,第二十大学,恩斯切德:施普林格-柏林-海德堡;1976年,第330-40页。
[9] Youngs,D.L.,具有大流体畸变的时间相关多材料流动,流体动力学中的数值方法,273-285,(1982),纽约学术出版社·Zbl 0537.76071号
[10] 所以,K.K。;胡晓云。;Adams,N.A.,两相可压缩流模拟的反扩散界面锐化技术,《计算物理杂志》,2314304-3323,(2012)·Zbl 1426.76428号
[11] 所以,K.K。;胡晓云。;Adams,N.A.,两相不可压缩流中界面变陡的反扩散方法,《计算物理杂志》,230,5155-5177,(2011)·Zbl 1416.76334号
[12] 萨斯曼,M。;Puckett,E.G.,用于计算三维和轴对称不可压缩两相流的耦合水平集和流体体积法,《计算物理杂志》,162,301-337,(2000)·Zbl 0977.76071号
[13] 王,Z。;杨,J。;顾,B。;Stern,F.,一种用于跌落破碎波锐界面模拟的耦合水平集和流体体积法,国际J多相流,35,227-246,(2009)
[14] Nguyen,V.-T。;Vu,D.-T。;公园,W.-G。;Jung,Y.-R.,使用双重时间伪压缩方法和流体界面体积跟踪算法对水冲击力进行数值分析,计算流体,103,18-33,(2014)·Zbl 1391.76486号
[15] Nguyen,V.-T。;Park,W.-G.,基于VOF方法的复杂三维水冲击问题的自由表面流求解器,国际J数值方法流体,82,3-34,(2015)
[16] Pilliod,J.E。;Puckett,E.G.,追踪材料界面的二阶精确流体体积算法,计算机物理杂志,199465-502,(2004)·兹比尔1126.76347
[17] 沙巴齐,K。;Paraschivoiu,M。;Mostaghimi,J.,基于三角网格重映射的二阶精确体积跟踪,《计算物理杂志》,188100-122,(2003)·Zbl 1127.76337号
[18] Dyadechko,V。;Shashkov,M.,《流体界面力矩重建》,(2005),洛斯阿拉莫斯国家实验室,技术代表LA-UR-05-7571
[19] 安,H.T。;Shashkov,M.,广义多面体网格上的多材料界面重建,计算物理杂志,2262096-2132,(2007)·Zbl 1388.76232号
[20] Owkes,M。;Desjardins,O.,《保守、三维、非分裂、几何输运的计算框架及其在流体体积(VOF)方法中的应用》,《计算物理杂志》,270,587-612,(2014)·Zbl 1349.76636号
[21] O.乌宾克。;Issa,R.I.,《捕捉任意网格上尖锐流体界面的方法》,《计算物理杂志》,153,26-50,(1999)·Zbl 0955.76058号
[22] Leonard,B.P.,《不连续性对流建模的简单高精度分辨率程序》,《国际数值方法流体》,81291-1318,(1988)·Zbl 0667.76125号
[23] Heyns,J.A。;马兰,A.G。;Harms,T.M。;Oxtoby,O.F.,《使用流体体积法模拟自由表面流动的压缩表面捕获公式的开发》,《国际数值方法流体》,71,788-804,(2013)·Zbl 1430.76471号
[24] Rusche,H.,高相分数分散两相流的计算流体动力学,(2002),帝国理工学院,博士论文
[25] Nguyen,V.-T。;Vu,D.-T。;Park,W.-G.,使用Navier-Stokes计算和移动Chimera网格方案对自由下落物体的入水进行三维数值模拟,(第十届ACFD会议,韩国大酒店,(2014),济州),85-90
[26] Nguyen,V.-T。;Vu,D.-T。;公园,W.-G。;Jung,Y.-R.,使用双重时间伪压缩方法和流体界面体积跟踪算法对水冲击力进行数值分析,计算流体,103,18-33,(2014)·Zbl 1391.76486号
[27] Zalesak,S.T.,流体的全多维通量修正传输算法,《计算物理杂志》,31335-362,(1979)·Zbl 0416.76002号
[28] Desjardins,O。;穆罗,V。;Pitsch,H.,《模拟湍流雾化的精确保守水平集/虚流体方法》,《计算物理杂志》,2278395-8416,(2008)·Zbl 1256.76051号
[29] Nguyen,V.-T。;Ha,C.-T。;Park,W.-G.,轴对称物体进出水的多相流模拟,(ASME 2013年会议记录,IMECE2013,美国加利福尼亚州,(2013),圣地亚哥)
[30] Nguyen,V.-T。;Vu,D.-T。;公园,W.-G。;Jung,C.-M.,在6自由度运动中使用移动Chimera网格法求解进水体的Navier-Stokes解算器,计算流体,140,19-38,(2016)·Zbl 1390.76602号
[31] Toro,E.F.,Riemann解算器和流体动力学数值方法。实用简介,(2009),施普林格出版社,柏林-海德堡·Zbl 1227.76006号
[32] 洛博夫斯克。;Botia-Vera,E。;卡斯特拉纳,F。;Mas-Soler,J。;Souto-Iglesias,A.,溃坝期间动压荷载的实验研究,《流体结构杂志》,48,407-434,(2014)
[33] Dressler,R.F.,水力溃坝波的理论和实验比较,国际科学与水文协会,319-328,(1954)
[34] De Backer,G。;万托雷,M。;比尔斯,C。;De Pré,J。;维克多·S。;De Rouck,J.,水对轴对称物体影响的实验研究,《应用海洋研究》,31143-156,(2009)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。