×

用于湍流充满气泡的通道流的高度可扩展DNS解算器。 (英语) Zbl 1410.76082号

摘要:我们提出了一种高效且高度可扩展的离散气液流动直接数值模拟(DNS)求解器,其中包含大量可变形气泡。我们将此应用于湍流中的(O(10^4))气泡。这是通过最先进的数值方法以及起源于单相和两相流的快速可扩展数值算法的精心组合实现的。介绍了曲率计算、气泡碰撞和变效率泊松方程求解器算法的特点和关键要素。通过与Hysing基准进行比较,确定了上升气泡准确平流的分辨率要求[S.Hysing公司等,《国际期刊数字》。方法流体60,No.11,1259–1288(2009;Zbl 1273.76276号)]. 曲率计算采用广义高度函数法。与有限差分方法相比,使用GHF时,我们观察到与理论收敛速度一致,并且伪速度显著降低。分析了气泡与区域壁之间的主要相互作用机制,表明了基本数值方法的二阶收敛性。对Navier-Stokes(NS)求解器和气体体积分数求解器的并行性能进行了详细分析。分析表明,在计算网格上,NS解算器具有10亿个单元,且气体体积分数解算器的理想缩放比例高达\(O(10^3)\)个气泡,近似线性缩放到18000个核心。除此之外,还发现了高达\(O(10^4)\)气泡的可接受开销。一个装有10000个气泡的湍流通道的下流结构模拟说明了计算能力。计算速度场的一阶和二阶统计以及统计稳态下平均气体体积分数场的分布。在所考虑的情况下,观察到壁剪切应力增加了47%,这是由嵌入的气泡引起的湍流修正引起的。在高体积分数下,湍流和气泡之间的相互作用导致通道核心宽区域内速度波动的均方根值强烈衰减。

MSC公司:

76层65 湍流的直接数值模拟和大涡模拟
76M20码 有限差分方法在流体力学问题中的应用
76T10型 液气两相流,气泡流
PDF格式BibTeX公司 XML格式引用
全文: DOI程序 链接

参考文献:

[1] Pope,S.B.,《湍流》(2001),IOP出版社
[2] Tryggvason,G。;达比里,S。;Aboulhasanzadeh,B。;Lu,J.,多相流直接数值模拟中的多尺度考虑a,Phys Fluids,25,3,031302,(2013)
[3] 卢,J。;Tryggvason,G.,垂直通道中湍流气泡下流的数值研究,《物理流体》,18,10,103302,(2006)
[4] 卢,J。;Tryggvason,G.,垂直通道中湍流气泡下流中气泡尺寸的影响,化学工程科学,62,11,3008-3018,(2007)
[5] 卢,J。;Tryggvason,G.,垂直通道中湍流气泡上流中气泡变形能力的影响,《物理流体》,20,4,040701,(2008)·Zbl 1182.76476号
[6] 博洛特诺夫,I.A。;Jansen,K.E。;Drew,D.A。;Oberai,A.A。;拉希,R.T。;Podowski,M.Z.,湍流两相泡状通道流动的分离直接数值模拟,国际J多道流,37,6,647-659,(2011)
[7] Bois,G.,《垂直通道中湍流泡状流的直接数值模拟:朝向改进的二阶雷诺应力模型》,Nucl Eng-Des,321,92-103,(2017)
[8] 桑塔雷利,C。;Fröhlich,J.,垂直湍流通道中球形气泡的直接数值模拟。气泡尺寸和双分散性的影响,《国际多相流杂志》,81,27-45,(2016)
[9] 斯潘丹,V。;梅斯奇尼,V。;Ostilla-Mónico,R。;Lohse,D。;Querzoli,G。;de Tullio,M.D.,《用于变形界面和膜的完全解析模拟的并行相互作用势方法与浸没边界法》,《计算物理杂志》,348567-590,(2017)
[10] Prosperetti,A。;Tryggvason,G.,多相流的计算方法,(2007),剑桥大学出版社·Zbl 1166.76003号
[11] 希特,C.W。;Nichols,B.D.,自由边界动力学的流体体积(VOF)方法,计算物理杂志,39,1,201-225,(1981)·Zbl 0462.76020号
[12] S.O.Unverdi。;Tryggvason,G.,粘性、不可压缩、多流体流动的前跟踪方法,《计算物理杂志》,100,1,25-37,(1992)·兹比尔0758.76047
[13] 萨斯曼,M。;斯梅雷卡,P。;Osher,S.,《计算不可压缩两相流解的水平集方法》,《计算物理杂志》,114,146-159,(1994)·Zbl 0808.76077号
[14] Shin,S。;Yoon,我。;Juric,D.,直接模拟二维和三维多相流的局部前沿重建方法,《计算物理杂志》,230,17,6605-6646,(2011)·Zbl 1408.76549号
[15] 萨斯曼,M。;Puckett,E.G.,计算三维和轴对称不可压缩两相流的耦合水平集和流体体积方法,J Comput Phys,162,2301-337,(2000)·Zbl 0977.76071号
[16] Rider,W.J。;Kothe,D.B.,重建体积跟踪,计算机物理杂志,141,2112-152,(1998)·Zbl 0933.76069号
[17] 斯卡多韦利,R。;Zaleski,S.,用最小二乘拟合和分裂欧拉-拉格朗日平流进行界面重建,国际数值方法流体,41,3,251-274,(2003)·Zbl 1047.76080号
[18] 利奥维奇,P。;Rudman,M。;Liow,J.-L。;莱克哈尔,D。;Kothe,D.,一种具有平面保持界面重建的三维非平面投影体积跟踪算法,计算流体,35,10,1011-1032,(2006)·Zbl 1177.76292号
[19] Cifani,P。;Michalek,W.R。;Priems,G.J.M。;Kuerten,J.G.M。;范德盖尔德,C.W.M。;Geurts,B.J.,VOF方法中表面压缩方法和界面重建方法的比较,计算流体,136,421-435,(2016)·Zbl 1390.76414号
[20] Noh,W.F。;Woodward,P.,Slic(简单线界面计算),第五届流体动力学数值方法国际会议论文集,恩舍德,330-340,(1976),施普林格·Zbl 0382.76084号
[21] Puckett,E.G。;Almgren,A.S。;贝尔,J.B。;马库斯,D.L。;Rider,W.J.,《用于跟踪变密度不可压缩流中流体界面的高阶投影方法》,《计算物理杂志》,130,2,269-282,(1997)·Zbl 0872.76065号
[22] 雷纳迪,Y。;Renardy,M.,Prost:流体体积法表面张力的抛物线重建,《计算物理杂志》,183,2400-421,(2002)·Zbl 1057.76569号
[23] 伦纳德,B.P。;Lock,A.P。;MacVean,M.K.,《保守的显式无限制时间步长多维常保平流方案》,《Mon Weather Rev》,124,11,2588-2606,(1996)
[24] 康明斯,S.J。;弗朗索瓦,M.M。;Kothe,D.B.,从体积分数估算曲率,计算结构,83,6,425-434,(2005)
[25] 巴尔图森,M.W。;Kuipers,J.A.M。;Deen,N.G.,《流体体积法表面张力模型的临界比较》,化学工程科学,109,65-74,(2014)
[26] Popinet,S.,表面张力数值模型,《流体力学年鉴》,50,1-28,(2018)
[27] 科亚吉,E。;Boersma,B.J.,用多标记正面捕获方法对液-液界面上液滴撞击的数值模拟,计算物理杂志,228,12,4444-4467,(2009)·Zbl 1395.76056号
[28] Cifani P.TBF求解器。10.5281/泽诺多1222098。;Cifani P.TBF求解器。10.5281/天诺多.1222098。
[29] Kim,J。;莫因,P。;Moser,R.,低雷诺数下充分发展的通道流中的湍流统计,《流体力学杂志》,177133-166,(1987)·Zbl 0616.76071号
[30] 莫瑟,R.D。;Kim,J。;Mansour,N.N.,湍流通道流动的直接数值模拟τ=590,流体物理学,11,4,943-945,(1999)·Zbl 1147.76463号
[31] Vreman,A.W。;Kuerten,J.G.M.,雷诺数下湍流通道流直接数值模拟数据库的比较τ=180,流体物理学,26,1015102,(2014)
[32] Vreman,A.W。;Kuerten,J.G.M.,湍流通道中速度和压力空间导数的统计,流体物理学,26,8,085103,(2014)
[33] Vreman,A.W.,《不可压缩Navier-Stokes方程的投影方法:无滑移壁附近的压力》,《计算物理杂志》,263,353-374,(2014)·Zbl 1349.76547号
[34] Dong,S。;Shen,J.,《大密度比两相不可压缩流相场模拟中包含常数矩阵的时间步进格式》,《计算物理杂志》,231,17,5788-5804,(2012)·Zbl 1277.76118号
[35] Popinet,S.,《表面张力驱动界面流的精确自适应解算器》,《计算物理杂志》,228,16,5838-5866,(2009)·Zbl 1280.76020号
[36] 海辛,S。;Turek,S。;库兹明,D。;北卡罗来纳州帕罗里尼。;伯曼,E。;Ganesan,S.,二维气泡动力学的定量基准计算,国际数值方法流体,60,11,1259-1288,(2009)·Zbl 1273.76276号
[37] 范德波尔,E.P。;Ostilla-Mónico,R。;Donners,J。;Verzicco,R.,《强湍流壁边界流动的铅笔分布有限差分代码》,计算流体,116,10-16,(2015)·Zbl 1390.76250号
[38] Kwakkel,M。;布鲁根,W.-P。;Boersma,B.J.,《利用聚结/破碎模型对液滴流的CLSVOF方法的扩展》,《计算物理杂志》,253,166-188,(2013)·Zbl 1349.76121号
[39] 达贡,L。;Menon,R.,Openmp:共享内存编程的行业标准API,IEEE计算科学工程,5,1,46-55,(1998)
[40] 格罗普,W。;Lusk,E。;剂量,N。;Skjellum,A.,MPI消息传递接口标准的高性能、可移植实现,并行计算,22,6,789-828,(1996)·Zbl 0875.68206号
[41] Lee,M.J。;做哦,B。;Kim,Y.B.,不可压缩Navier-Stokes方程的正则分步方法和一致边界条件,计算物理杂志,168,1,73-100,(2001)·Zbl 1074.76585号
[42] Wesseling,P.,《计算流体动力学原理》,29,(2009),Springer Science&Business Media·Zbl 1185.76005号
[43] Rai,M.M。;Moin,P.,《使用有限差分格式直接模拟湍流》,《计算物理杂志》,96,1,15-53,(1991)·Zbl 0726.76072号
[44] Verzicco,R。;Orlandi,P.,柱坐标下三维不可压缩流动的有限差分格式,《计算物理杂志》,123,2,402-414,(1996)·兹伯利0849.76055
[45] Chan,T.F。;Wan,W.L.,非光滑系数椭圆线性系统的鲁棒多重网格方法,计算应用数学杂志,123,1,323-352,(2000)·兹伯利0966.65099
[46] Leonard,B.P.,基于二次上游插值的稳定和精确对流建模程序,Comput Methods Appl Mech Eng,19,1,59-98,(1979)·Zbl 0423.76070号
[47] 库兰特,R。;弗里德里希斯(Friedrichs,K.)。;Lewy,H.,《数学年鉴》,100,1,32-74,(1928)
[48] Brackbill,J.U。;科特,D.B。;Zemach,C.,《模拟表面张力的连续介质方法》,《计算物理杂志》,100,2,335-354,(1992)·Zbl 0775.76110号
[49] 丹尼·F。;van Wachem,B.G.M.,毛细波对数值时间步长的限制,《计算物理杂志》,285,24-40,(2015)·Zbl 1351.76028号
[50] Wesseling P.1995年。多重网格方法简介(编号ICASE-95-11)。弗吉尼亚州汉普顿科学与工程计算机应用研究所。;Wesseling P.1995年。多重网格方法简介(编号ICASE-95-11)。弗吉尼亚州汉普顿科学与工程计算机应用研究所。
[51] 多德,M.S。;Ferrante,A.,《不可压缩双流体流动的快速压力校正方法》,《计算物理杂志》,273,416-434,(2014)·Zbl 1351.76161号
[52] 卡努托,C。;侯赛尼,M.Y。;Quarteroni,A.M。;Thomas Jr,A.,《流体动力学中的光谱方法》(2012),Springer Science&Business Media
[53] 弗朗索瓦,M.M。;康明斯,S.J。;Dendy,E.D。;科特,D.B。;西西里岛,J.M。;Williams,M.W.,体积跟踪框架内连续和尖锐界面表面张力模型的平衡力算法,《计算物理杂志》,213,1,141-173,(2006)·Zbl 1137.76465号
[54] 医学博士托里。;克劳特曼,L.D。;Mjolsness,R.C。;Hirt,C.,Nasa-vof2d:自由表面不可压缩流动的计算机程序,技术报告,(1985),洛斯阿拉莫斯国家实验室,新墨西哥州(美国)
[55] Harvie,D.J.E。;Davidson,M.R。;Rudman,M.,流体体积模拟中寄生电流产生的分析,应用数学模型,30,10,1056-1066,(2006)·Zbl 1163.76401号
[56] 潘,K.-L。;法律,C.K。;Zhou,B.,迎面双液滴碰撞中合并和反弹的实验和机理描述,应用物理学杂志,103,6,064901,(2008)
[57] 阿夫哈米,S。;Bussmann,M.,将接触角应用于三维VOF模拟的高度函数,国际数值方法流体,61,8,827-847,(2009)·Zbl 1252.76057号
[58] 卢,J。;比斯瓦斯,S。;Tryggvason,G.,《垂直通道中层流泡状流的DNS研究》,《国际多道流杂志》,32,6,643-660,(2006)·Zbl 1136.76567号
[59] 弗里戈,M。;Johnson,S.G.,《西方最快的傅里叶变换》,技术报告,(1997),马萨诸塞州理工学院
[60] 弗里戈,M。;Johnson,S.G.,《FFTW3的设计与实现》,美国电气与电子工程师协会(Proc IEEE),93,2,216-231,(2005)
[61] 安德森,E。;Bai,Z。;比肖夫,C。;南卡罗来纳州布莱克福德。;德梅尔,J。;Dongarra,J.,《LAPACK用户指南》(1999),宾夕法尼亚州费城工业和应用数学学会·Zbl 0934.65030号
[62] 王世凯。;Lee,S.J。;琼斯,O.C。;Lahey,R.T.,泡状两相流中的三维湍流结构和相位分布测量,国际多相流杂志,13,3,327-343,(1987)
[63] 卡辛斯基,O.N。;Randin,V.V.,《垂直管道中的向下气泡气液流》,《国际多相流杂志》,25,1,109-138,(1999)·Zbl 1137.76628号
[64] 斯卡多韦利,R。;Zaleski,S.,《矩形网格中连接线性界面和体积分数的分析关系》,《计算物理杂志》,164,1,228-237,(2000)·Zbl 0993.76067号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。