湍流热对流DNS多分辨率并行直接求解方法

doi:10.19596/j.cnki.1001-246x.7974

计算物理 ›› 2019, Vol. 36 ›› Issue (6): 641-647.DOI: 10.19596/j.cnki.1001-246x.7974

湍流热对流DNS多分辨率并行直接求解方法

包芸, 何建超, 方明卫

中山大学航空航天学院, 广东广州 510275

收稿日期:2018-10-08 修回日期:2018-12-25 出版日期:2019-11-25 发布日期:2019-11-25
作者简介:包芸(1960-),女,教授,从事计算流体力学研究,E-mail:stsby@mail.sysu.edu.cn
基金资助:
国家自然科学基金（11772326）资助项目

Multiple-resolution Parallel Direct Method of DNS for Turbulence Thermal Convection

BAO Yun, HE Jianchao, FANG Mingwei

School of Aeronautics and Astronautics, Sun Yat-sen University, Guangzhou, Guangdong 510275, China

Received:2018-10-08 Revised:2018-12-25 Online:2019-11-25 Published:2019-11-25

摘要/Abstract

摘要： 依据温度标量场与动量计算的空间和时间计算分辨率不同的特点，采用两套网格，建立多分辨率双网格并行直接求解方法，用以解决极高Ra数湍流热对流DNS模拟巨大计算工作量的难题.在两套网格的数据交换上，根据每个细网格都满足连续方程，设计了速度的守恒平移插值方法.二维极高Ra数湍流热对流的计算结果表明，采用多分辨率双网格并行直接求解方法的DNS计算，可以使计算工作量降低近一个量级.瞬时温度场显示，双网格方法的计算结果可以很好地描述极高Ra数下快速运动的小尺寸漩涡团状羽流，得到的结果与原网格一致，不同方法计算得到的传热Nu数误差不超过1%.

关键词: Rayleigh-Bénard热对流, 湍流, 极高Ra数, 多分辨率, MrPDM-DNS

Abstract: Spatial and time resolution of thermal scalar field and velocity fields are different in turbulence thermal convection. A multiple-resolution parallel direct method for DNS with two different girds is presented to solve the problem of huge computational cost in simulation of turbulence thermal convection with very high Rayleigh numbers. A constant translation interpolation method of speed was designed to meet continuous equations at each finer gird in data transformation between coarser gird and finer gird. Simulations of 2D turbulence thermal convection with very high Rayleigh numbers show that computational cost is reduced an order of magnitude by using this method. Small-scale eddy-like plumes, that move rapidly in thermal convection, are well performed in instantaneous temperature fields. It is consistent with results with a single gird. Difference of Nusselt number obtained by two methods is below 1%.

Key words: Rayleigh-Bénard convection, turbulence, high Rayleigh number, multiple-resolution, MrPDM-DNS

中图分类号:

O358

包芸, 何建超, 方明卫. 湍流热对流DNS多分辨率并行直接求解方法[J]. 计算物理, 2019, 36(6): 641-647.

BAO Yun, HE Jianchao, FANG Mingwei. Multiple-resolution Parallel Direct Method of DNS for Turbulence Thermal Convection[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(6): 641-647.

参考文献

[1] AHLERS G, GROSSMANN S, LOHSE D. Heat transfer and large scale dynamics in turbulent Rayleigh-Bénard convection[J]. Reviews of Modern Physics, 2009, 81(2):503-537.
[2] STEVENS R J A M, ROBERTO V, DETLEF L. Radial boundary layer structure and Nusselt number in Rayleigh-Bénard convection[J]. Journal of Fluid Mechanics, 2010, 643(3):495-507.
[3] VAN DER POEL E P, STEVENS R J A M, LOHSE D. Comparison between two-and three-dimensional Rayleigh-Bénard convection[J]. Journal of Fluid Mechanics, 2013, 736:177-194.
[4] 包芸, 宁浩, 徐炜. 湍流热对流大尺度环流反转时的角涡特性[J]. 物理学报, 2014, 63(15):154703.
[5] CHONG K L, HUANG S D, KACZOROWSKI M, et al. Condensation of coherent structures in turbulent flows[J]. Physical Review Letters, 2015, 115(26):264503.
[6] BAO Y, CHEN J, LIU B F, et al. Enhanced heat transport in partitioned thermal convection[J]. Journal of Fluid Mechanics, 2015, 784, R5.
[7] HE P, BAO Y. Effect of free-slip boundaries on flow and heat transport characteristics in two-dimensional Rayleigh-Bénard convection[J]. Chinese Journal of Conputational Physics, 2019, 36(5):542-550.
[8] LI Beibei, WANG Tingting, CHEN Jian,et al. Lattice Boltzmann study of nonlinear characteristics of double diffusive mixed convection in an enclosure[J]. Chinese Journal of Computational Physics, 2016, 33(2):156-162.
[9] XIE Y C, XIA K Q. Turbulent thermal convection over rough plates with varying roughness geometries[J]. Journal of Fluid Mechanics, 2017, 825:573-599.
[10] ZHU X, RJAM S, VERZICCO R, et al. Roughness-facilitated local 1/2 scaling does not imply the onset of the ultimate regime of thermal convection[J]. Physical Review Letters, 2017, 119(15):154501.
[11] ZHANG Y Z, SUN C, BAO Y, et al. How surface roughness reduces heat transport for small roughness heights in turbulent Rayleigh-Bénard convection[J]. Journal of Fluid Mechanics, 2018, 836.
[12] ZHU X, MATHAI V, STEVENS R, et al. Transition to the ultimate regime in two-dimensional Rayleigh-Bénard convection[J]. Physical Review Letters, 2018, 120(14):144502.
[13] OSTILLA-MONICO R, YANG Y, POEL E P V D, et al. A multiple-resolution strategy for direct numerical simulation of scalar turbulence[J]. Journal of Computational Physics, 2015, 301(C):308-321.
[14] BAO Y, LUO J, YE M. Parallel direct method of DNS for two-dimensional turbulent Rayleigh-Bénard convection[J]. Journal of Mechanics, 2018,34(2):159-166.
[15] BAO Yun, YE Mengxiang, LUO Jiahui.An efficient parallel direct method for turbulent thermal convection[J]. Chinese Journal of Conputational Physics, 2017, 34(6):651-656.
[16] 何鹏, 黄茂静, 包芸. 二维湍流Rayleigh-Bénard对流温度边界层脉动影响特性研究[J]. 中国科学G:物理学力学天文学, 2018.
[17] 黄茂静, 包芸. 湍流热对流近底板流态与温度边界层特性[J]. 物理学报, 2016, 65(20):149-156.
[18] 包芸, 高振源, 叶孟翔. 湍流热对流Prandtl数效应的数值研究[J]. 物理学报, 2018,67(1):158-164.
[19] 黄茂静, 包芸. 湍流Rayleigh-Bénard热对流温度剖面中对数律研究[J]. 中国科学G:物理学力学天文学, 2017, 47(6):064701.
[20] SHARIMAN B I, SIGGIA E D. Heat transport in high-Rayleigh-number convection[J]. Phys Rev A, 1990, 42(6):3650.

湍流热对流DNS多分辨率并行直接求解方法

Multiple-resolution Parallel Direct Method of DNS for Turbulence Thermal Convection

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	牛霄, 倪国喜, 马文铧. 自适应多分辨率方法在反应多相流数值模拟中的应用[J]. 计算物理, 2020, 37(6): 639-652.
[2]	周志超, 许凌飞, 任天荣, 顾村锋. 基于GCV-FFT方法的超声速压缩拐角流场气动光学效应计算[J]. 计算物理, 2020, 37(3): 284-298.
[3]	何鹏, 包芸. 滑移边界对二维Rayleigh-Bénard热对流流态和传热特性的影响[J]. 计算物理, 2019, 36(5): 542-550.
[4]	蔡林峰, 李昊歌, 陈伟芳. 无限展长后掠翼边界层定常横流不稳定转捩分析[J]. 计算物理, 2019, 36(2): 175-181.
[5]	李兆辉, 时钟. 稳定分层二层流非湍流/湍流层无平均剪切密度界面处的湍流[J]. 计算物理, 2018, 35(6): 631-648.
[6]	张学良, 谭惠丽, 唐国宁, 邓敏艺. 心肌组织机械形变的元胞自动机模拟[J]. 计算物理, 2018, 35(3): 294-302.
[7]	潘宏禄, 李俊红, 程晓丽, 马汉东. 气动光学效应湍流脉动模型系数修正[J]. 计算物理, 2018, 35(2): 194-204.
[8]	包芸, 叶孟翔, 罗嘉辉. 湍流热对流的高效并行直接求解方法[J]. 计算物理, 2017, 34(6): 651-656.
[9]	李俊红, 张亮, 俞继军, 程晓丽. 高超声速可压缩流中粗糙壁热流研究[J]. 计算物理, 2017, 34(2): 165-174.
[10]	赵亚洲, 马智博. SPH核函数光滑长度最优选取和自适应准则研究[J]. 计算物理, 2017, 34(1): 29-38.
[11]	谷建法, 戴振生, 古培俊, 叶文华, 郑无敌, 邹士阳. 点火靶驱动不对称性与表面粗糙度的模耦合模拟[J]. 计算物理, 2016, 33(6): 645-651.
[12]	潘宏禄, 关发明, 袁湘江, 卜俊辉. 高超声速平板边界层/圆柱粗糙元非定常干扰[J]. 计算物理, 2015, 32(5): 537-544.
[13]	甘文彪, 周洲. 基于层流动能湍流模型的数值模拟方法[J]. 计算物理, 2013, 30(2): 169-179.
[14]	聂德明, 林建忠, 黄利忠. 旋转体系中湍流的LBM模拟[J]. 计算物理, 2013, 30(1): 11-18.
[15]	王浩明, 赵海波, 郑楚光. 考虑双向耦合效应的格子Boltzmann气固两相湍流模型[J]. 计算物理, 2013, 30(1): 19-26.