×

基于通用曲线网格的流体体积与液位组耦合方法,具有精确的表面张力计算。 (英语) Zbl 1452.76111号

小结:本文结合VOF方法和LS方法的优点,提出了一种通用曲线网格上的流体体积与水平集耦合方法,用于界面流动模拟。在该方法中,引入了迭代几何运算来计算物理域和计算域上的水平集函数。特别是,在物理域上成功地实现了距离函数的直接构造,从而避免了计算域逆变换引起的数值误差。推导了基于水平集的曲线网格上的连续表面张力(CSF)模型,以精确地模拟表面张力,然后将该模型以及所提出的耦合体积跟踪方法纳入不可压缩Navier-Stokes解算器中。最后,在不规则区域上对时间反转单涡流、静态液滴、溃坝和单气泡上升四个基准问题进行了测试,以验证该方法处理复杂区域的准确性和鲁棒性。结果表明,该方法与前人的实验数据或数值计算结果吻合良好,并具有良好的质量守恒特性,表明该方法能够准确模拟不可压两相流。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
76T10型 液气两相流,气泡流
76D45型 不可压缩粘性流体的毛细管(表面张力)
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Kozyrakis,G。;Delis,A。;Kampanis,N.,复杂区域中不可压缩Navier-Stokes流动的有限差分求解器,Appl。数字。数学。,115, 275-298 (2017) ·Zbl 1422.65159号
[2] Quirk,J.J.,《计算任意复杂二维物体周围气体动力流动的非结构化网格的替代方法》,《计算》。流体,23,1,125-142(1994)·Zbl 0788.76067号
[3] 尼科尔斯,B。;希尔特,C。;Hotchkiss,R.,《Sola-vof:多自由边界瞬态流体流动的求解算法》(1980),洛斯阿拉莫斯科学实验室,技术代表。
[4] 科特,D。;Rider,W。;莫索,S。;布罗克·J。;Hochstein,J.,具有二维和三维表面张力的界面体积跟踪,(第34届航空航天科学会议和展览(1996)),859
[5] Nikseresht,A。;Alishahi,M。;Emdad,H.,使用计算域中拉格朗日传播的三维vof完成acv(气割车辆)周围的流场计算,Compute。结构。,86, 7-8, 627-641 (2008)
[6] López,J。;Hernández,J.,《通用网格中流体三维体积方法的分析和几何工具》,J.Compute。物理。,227, 12, 5939-5948 (2008) ·Zbl 1142.76042号
[7] 萨古,A。;文森特,S。;Caltagirone,J.-P.,浸没界面和边界的二阶曲线到笛卡尔变换。虚拟域和多相流的应用,计算。流体,46,1,422-428(2011)·Zbl 1431.76105号
[8] Nguyen,V.-T。;Park,W.-G.,两相不可压缩流动的流体体积(vof)界面锐化方法,计算。流体,152104-119(2017)·Zbl 1390.76498号
[9] Nguyen,V.-T。;Park,W.-G.,基于改进的流体体积法的Navier-Stokes解算器增强,用于复杂界面流动模拟,应用。海洋研究,72,92-109(2018)
[10] 奥舍,S。;Sethian,J.A.,《以曲率相关速度传播的前沿:基于哈密尔顿-雅可比公式的算法》,J.Compute。物理。,79, 1, 12-49 (1988) ·Zbl 0659.65132号
[11] Yue,W。;林,C.-L。;Patel,V.C.,用水平集方法对非定常多维自由表面运动进行数值模拟,国际期刊Numer。方法流体,42,853-884(2003)·Zbl 1078.76592号
[12] 王,Z。;邹强。;Reeve,D.,使用两相流cfd模型模拟溢出破碎波,计算。流体,38,10,1995-2005(2009)·Zbl 1242.76184号
[13] 萨斯曼,M。;Puckett,E.G.,用于计算三维和轴对称不可压缩两相流的耦合水平集和流体体积法,J.Compute。物理。,162, 2, 301-337 (2000) ·Zbl 0977.76071号
[14] 儿子,G。;Hur,N.,一种用于流体颗粒浮力驱动运动的耦合水平集和流体体积方法,Numer。热传输。,B部分,Fundam。,42, 6, 523-542 (2002)
[15] Haghshenas,M。;Wilson,J.A。;Kumar,R.,《表面张力主导两相流的代数耦合水平集-流体体积法》,国际J.Multiph。流量,90,13-28(2017)
[16] 斯科里什,M。;Garmory,A。;Dianat,M.,作为耦合水平集流体求解器体积的一部分,在一般凸网格中PLIC的迭代界面重建方法,J.Compute。物理。,368, 254-276 (2018) ·Zbl 1392.76053号
[17] 杨,X。;詹姆斯·A·J。;Lowengrub,J。;郑,X。;Cristini,V.,非结构化三角形网格的自适应耦合水平集/流体界面体积捕获方法,J.Comput。物理。,217, 2, 364-394 (2006) ·Zbl 1160.76377号
[18] 法拉利,A。;马格尼尼,M。;Thome,J.R.,在非均匀和非结构化网格中模拟微尺度两相流的柔性耦合液位集和流体体积(FLEXCLV)方法,国际J.Multiph。流量,91,276-295(2017)
[19] Sun,D。;Tao,W.,《计算不可压缩两相流的流体体积与水平集(VOSET)耦合方法》,《国际热质交换杂志》。,53, 4, 645-655 (2010) ·Zbl 1303.76126号
[20] Ling,K。;李,Z.-H。;Sun,D.-L。;何,Y.-L。;陶文清,不可压缩两相流的三维流体体积&液位集(VOSET)方法,计算。流体,118293-304(2015)·Zbl 1390.76679号
[21] 曹,Z。;Sun,D。;Yu,B。;Wei,J.,一种基于非结构化三角形网格重映射算法的流体体积与水平集耦合(VOSET)方法,国际热质传输杂志。,111232-245(2017)
[22] 曹,Z。;Sun,D。;魏杰。;Yu,B.,一种基于多维平流的非结构化三角网格流体体积与水平集耦合方法,Chem。工程科学。,176, 560-579 (2018)
[23] 曹,Z。;Sun,D。;Yu,B。;Wei,J.,基于非结构化四边形网格解析PLIC的流体体积与水平集耦合方法,数值。热传输。,B部分,Fundam。,73,4189-205(2018)
[24] 王,Z。;杨,J。;Stern,F.,在一般结构网格上构造距离函数的新流体体积法,J.Compute。物理。,231, 9, 3703-3722 (2012) ·Zbl 1402.65091号
[25] 王,Z。;杨,J。;顾,B。;Stern,F.,《用于模拟冲击破碎波锐界面的耦合水平集和流体体积法》,国际J·多波。流量,35,3227-246(2009)
[26] Brackbill,J。;科特,D.B。;Zemach,C.,《模拟表面张力的连续体方法》,J.Compute。物理。,100, 2, 335-354 (1992) ·Zbl 0775.76110号
[27] Yokoi,K.,基于CLSVOF方法、多矩方法和密度尺度CSF模型的自由表面流动实用数值框架:液滴飞溅的数值模拟,J.Compute。物理。,232, 1, 252-271 (2013)
[28] Shyy,W.,《复杂流体流动和传热的基于压力的计算算法要素》,(《传热进展》,第24卷(1994),爱思唯尔出版社),191-275
[29] 贝尔,J.B。;科尔拉,P。;Glaz,H.M.,《不可压缩Navier-Stokes方程的二阶投影法》,J.Compute。物理。,85, 2, 257-283 (1989) ·兹比尔0681.76030
[30] Wang,T。;李,H。;Feng,Y。;Shi,D.,动态自适应四叉树网格上的流体体积与水平集(VOSET)耦合方法,国际热质传输杂志。,67, 1, 70-73 (2013)
[31] J.C.马丁。;莫伊斯,W.J。;Penney,W.G。;价格,A。;Thornhill,C.,第四部分。刚性水平面上液柱坍塌的实验研究,Philos。事务处理。R.Soc.伦敦。A、 244882312-324(1952)
[32] 海辛,S.-R。;Turek,S。;库兹明,D。;北卡罗来纳州帕罗里尼。;伯曼,E。;Ganesan,S。;Tobiska,L.,《二维气泡动力学的定量基准计算》,国际期刊数值。《液体方法》,60,11,1259-1288(2009)·Zbl 1273.76276号
[33] Chiu,P.-H。;Lin,Y.-T.,求解不可压缩两相流的保守相场方法,J.Comput。物理。,230, 1, 185-204 (2011) ·Zbl 1427.76201号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。