×

一种基于对称保护的二阶时间精确PISO方法。 (英语) Zbl 1521.76435号

摘要:本文提出了一种新的保守对称保二阶时间精确PISO压力-速度耦合方法,用于求解非结构并置网格上的不可压缩Navier-Stokes方程。这种隐式时间步进的新方法是Trias等人[35]针对显式时间步和非结构化并置网格的保守对称保增量压力投影方法的扩展。为了评估和比较这两种方法,我们在开源代码OpenFOAM中的一个统一求解器中实现了它们,在该代码中,我们使用Butcher数组来指定Runge-Kutta方法。因此,通过改变Butcher数组的条目,显式和对角隐式Runge-Kutta格式可以组合成一个求解器。我们使用泰勒-格林涡和盖驱动空腔流测试案例评估了所实现离散化方法的能量守恒特性和所选Runge-Kutta方案的时间一致性。最后,我们使用了一个更复杂的湍流通道流动测试案例,以进一步评估所提出的新的基于保守对称保持增量压力PISO的方法的性能。
虽然这两种实现的方法都基于对称保护离散化,但我们表明,当直接从泊松方程求解总压时,它们仍然会产生少量的数值耗散。当使用增量压力方法时,其中压力校正是由泊松方程求解的,这两种方法都是有效的完全保守方法。对于不可压缩湍流的高保真模拟,非常需要使用完全保守的方法。因此,对于此类模拟,所提出的数值方法预计具有较大的附加值,因为它们为在复杂几何体中执行真正节能的高保真模拟铺平了道路。此外,这两种方法都是在OpenFOAM中实现的,OpenFOAM在CFD社区中得到了广泛应用,因此该社区的很大一部分人可以从开发和实现的数值方法中受益。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
6500万08 含偏微分方程初值和初边值问题的有限体积法
65升06 常微分方程的多步、Runge-Kutta和外推方法
76D05型 不可压缩粘性流体的Navier-Stokes方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 解析-流露。用户指南,14.0版。ANSYS公司。;黎巴嫩;2011
[2] Archambeau,F。;梅奇托,N。;Sakiz,M.,Code saturne:计算湍流不可压缩流的有限体积代码:工业应用,国际有限体积杂志,1,1-62(2004)·Zbl 1490.76137号
[3] 阿舍尔,美国。;Ruuth,S。;Spiteri,R.,含时偏微分方程的隐式显式Runge-Kutta方法,应用数值数学,25,2-3,151-167(1997)·Zbl 0896.65061号
[4] 伯德·R。;斯图尔特,W。;Lightfoot,E.N.,《运输现象》(1960),John Wiley&Sons Inc。
[5] van der Blij,F.,《适合转变的细胞形状》(2007),荷兰格罗宁根国立大学
[6] Butcher,J.,《关于高阶Runge-Kutta过程》,《澳大利亚数学学会杂志》,4,2,179-194(1964)·兹比尔0244.65046
[7] 卡斯塔诺,S。;彼得罗尼奥,A。;Petris,G。;Armenio,V.,《使用OpenFOAM评估湍流大涡模拟的求解算法》,流体,4171,1-25(2019)
[8] 卡斯蒂格利奥尼,G。;Domaradzki,J.,使用低阶可压缩Navier-Stokes解算器进行真实几何形状的流动模拟中的数值耗散率和粘度,计算机与流体,119,37-46(2015)·Zbl 1390.76137号
[9] Chorin,A.,Navier-Stokes方程的数值解,计算物理杂志,22745-762(1968)·Zbl 0198.50103号
[10] D’Alessandro,V。;Binci,L。;Montelpare,S。;Ricci,R.,关于基于显式和隐式高阶Runge-Kutta格式的OpenFOAM求解器的开发,《计算物理通讯》,22214-30(2018)·Zbl 07693031号
[11] Felten,F。;Lund,T.,与不可压缩流的并置网格方案相关的动能守恒问题,《计算物理杂志》,215465-484(2006)·Zbl 1173.76382号
[12] Ferziger,J。;Perić,M.,《流体动力学的计算方法》(1997),柏林斯普林格出版社·Zbl 0869.76003号
[13] 海尔,E。;诺塞特,S。;Wanner,G.,《求解常微分方程》。I nonstiff problems(2008年),Springer-Verlag
[14] 火腿,F。;留置权,F。;Strong,A.,非均匀网格上不可压缩流动的完全保守二阶有限差分格式,《计算物理杂志》,177(1),117-133(2002)·兹比尔1066.76044
[15] 希肯,J。;火腿,F。;Militzer,J。;Koksal,M.,完全保守方法的偏移变换:复杂非结构网格上的湍流模拟,计算机物理,208704-734(2005)·Zbl 1329.76134号
[16] Issa,R.,用算子分裂法求解隐式离散流体流动方程,《计算物理杂志》,62,40-65(1985)·Zbl 0619.76024号
[17] van Kan,J.,粘性不可压缩流的二阶精确压力校正方法。,SIAM科学与统计计算杂志,7870-891(1986)·Zbl 0594.76023号
[18] 卡泽米·卡米亚布,V。;van Zuijlen,A。;Bijl,H.,高阶隐式Runge-Kutta格式在不可压Navier-Stokes方程并置有限体积离散化中的分析和应用。,计算机与流体,107-114(2015)·兹比尔1390.76455
[19] 科门,E。;卡米洛,L。;Shams,A。;Geurts,B。;Koren,B.,《准DNS和欠分辨率DNS中数值耗散的量化方法,以及湍流通道流准DNS和低分辨率DNS中数字耗散的影响》,《计算物理杂志》,345565-595(2017)·Zbl 1378.76035号
[20] 科门,E。;弗雷德里克斯,E。;科彭,T。;D’Alessandro,V。;Kuerten,J.,OpenFOAM中非结构化并置有限体积法中不同Runge-Kutta和速度插值方法的数值耗散率分析,计算物理通信,253,1-24(2020)·Zbl 07687597号
[21] 科门,E。;Shams,A。;卡米洛,L。;Koren,B.,OpenFOAM针对不同网格类型的准DNS功能,《计算机与流体》,96,87-104(2014)·兹比尔1391.76227
[22] 莫德斯蒂,D。;Pirozzoli,S.,非结构化网格上湍流可压缩流动的低耗散解算器,OpenFOAM实现,计算机与流体,152,14-23(2017)·Zbl 1390.76492号
[23] 蒙特基亚,M。;Brethouwer,G。;Wallin,S。;Johansson,A.V。;Knacke,T.,通过最小化数值耗散和使用显式代数SGS应力模型来改进OpenFOAM的LES,J Turbul,20:11-12,697-722(2019)
[24] Morinishi,Y。;Lund,T。;O.V.瓦西利耶夫。;Moin,P.,不可压缩流的完全保守高阶有限差分格式,J Comput Phys,14390-124(1998)·Zbl 0932.76054号
[25] Morinishi,Y。;O.瓦西利耶夫。;Ogi,T.,用于不可压缩流动模拟的柱坐标下完全保守有限差分格式,《计算物理杂志》,197,686-710(2004)·Zbl 1079.76602号
[26] Moser,R。;Kim,J。;Mansour,N.,《湍流通道流动的直接数值模拟》,Re{}_\tau=590,流体物理学,11,943-945(2008)·Zbl 1147.76463号
[27] Nieuwstadt,F。;Boersma,B。;J.Westerweel,《湍流》。《湍流理论与应用导论》(2016),斯普林格国际出版公司
[28] OpenFOAM®。用户指南,5.0版。OpenFOAM®基金会;2017
[29] Pascau,A.,非定常Navier-Stokes方程结构化共定位网格中的细胞面速度替代,国际数值方法流体杂志,65,812-833(2011)·Zbl 1428.76173号
[30] Rhie,C。;Chow,W.,《带后缘分离的孤立翼型湍流数值研究》。,AIAA J,211525-1532(1983)·Zbl 0528.76044号
[31] 桑德斯,B。;Koren,B.,应用于不可压缩Navier-Stokes方程的显式Runge-Kutta方法的精度分析,计算物理杂志,2313041-3063(2012)·Zbl 1402.65116号
[32] 沙申克;拉尔森,J。;Iacarino,G.,一个在无粘极限下具有精确质量、动量和动能守恒的共存不可压缩Navier-Stokes解算器,《计算物理杂志》,2294425-4430(2010)·兹比尔1305.76080
[33] 沈,C。;Sun,F。;Xia,X.,在openFOAM框架中实现基于密度的所有速度求解器,计算物理杂志,1852730-2741(2014)·Zbl 1360.76010号
[34] STAR-CCM+。用户指南,8.06版。CD-adapco;伦敦;2013
[35] Trias,F。;O.莱姆库尔。;奥利瓦,A。;佩雷兹·塞加拉,C。;Verstappen,R.,并置非结构化网格上navier-stokes方程的保对称离散化,J Comput Phys,258246-267(2014)·Zbl 1349.65386号
[36] Tuković,Z。;佩里奇,M。;Jasak,H.,一致二阶时间精确非迭代piso算法。,计算机与流体,166,78-85(2018)·Zbl 1390.76525号
[37] Verstapen,R。;Veldman,A.,湍流的对称性保持离散化。,《计算物理杂志》,187343-368(2003)·Zbl 1062.76542号
[38] Verstapen,R。;Wissink,J。;Veldman,A.,驱动空腔流动的直接数值模拟。,应用科学研究,51,377-381(1993)·Zbl 0800.76205号
[39] Versteeg,H。;Malalasekera,W.,《计算流体动力学导论:有限体积法》(2007),第二版,培生教育有限公司
[40] de Villiers,E.,《大涡模拟对有壁流动建模的潜力》(2006),帝国理工学院
[41] Vreman,A。;Kuerten,J.,《(R e\tau=180)下湍流通道流直接数值模拟数据库的比较》,流体物理学,26,015102,1-21(2014)
[42] Vreman,B.,可压缩湍流混合层的直接和大涡模拟(1995),特温特大学
[43] Vuorinen,V。;Keskinen,J。;杜维格,C。;Boersma,B.,《关于使用OpenFOAM®实现时间相关流的低耗散Runge-Kutta投影方法》,《计算机与流体》,93,153-163(2014)·Zbl 1391.76006号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。