基于格子Boltzmann方法的封闭空腔内自然对流数值模拟研究.pdf.pdfVIP

第 9 卷第 5 期 过 程 工 程 学 报 Vol.9 No.5 2009 年 10 月 The Chinese Journal of Process Engineering Oct. 2009 基于格子 Boltzmann 方法的封闭三角腔自然对流的数值模拟 1 1,2 1 雍玉梅 ， 杨 超 ， 毛在砂 (1. 中国科学院过程工程研究所绿色过程与工程重点实验室，北京 100190 ；2. 江苏省海洋资源开发研究院，江苏 连云港 222005) 摘 要：建立了二维不可压缩 D2G9 格子 Boltzmann 模型，耦合二维 TD2Q5 热格子 Boltzmann 模型，在非平衡态外 推的边界条件下，首先对不同 Eckert 数(Ec)和 Prandtl 数(Pr)时 Couette 流的温度场进行数值模拟，计算结果与解析解 吻合良好，且在 Ec 变化很大的条件下，计算结果仍与解析解相符，验证了模型的准确性和稳定性. 然后对封闭三角 空腔内不同 Rayleigh 数(Ra)下的自然对流流场和温度场进行了数值模拟，结果与文献计算值吻合良好，说明格子 Boltzmann 方法的 TD2Q5 热模型可用于高 Ra 时的空腔热流动模拟. 关键词：格子 Boltzmann 方法；不可压缩；热模型；数值模拟；空腔流 中图分类号：TQ021.3 文献标识码：A 文章编号：1009−606X(2009)05−0841−07 1 前 言 增加了计算量. He 等[7]直接在平衡态分布函数中忽略了 密度变化引起的高阶马赫数项，独立动力学变量是压力 空腔流动是航空工程中常遇到的问题，如飞行器武 而不是密度，但需增加额外的限制条件. Guo 等[8]建立的 器舱等，具有很高的工程实用价值. 另外由于空腔流动 D2G9 不可压缩 LBM 模型的平衡态函数之和保证为一 涉及非定常流、流体动力不稳定性、声−涡相互作用、 个常数，在一定程度上克服了以上模型的不足. 自然对流等许多基本的前沿问题，因而也具有广泛的理 LBM 中的反弹格式边界条件[9]非常适合模拟复杂 论研究意义[1]. 封闭空腔内流动和换热在电子器件冷 几何尺寸的流体流动，但只能在边界处给出一阶精度， 却、太阳能采集器、建筑物绝热、动力电站封闭母线及 会降低 LBM 在流体内部的二阶精度. 为了弥补格式精 旋转电机散热等领域有广泛的应用[2]. 封闭腔内的自然 度不高带来的缺点，陆续提出了动力学格式、外推方法 对流换热是计算流体力学与数值传热学研究的经典课 和非平衡态外推方法等具有二阶精度的无滑移速度边 题之一，人们用实验和传统的数值方法研究了不同几何 界条件. 动力学格式利用宏观物理量，在边界上直接求 区域内的热传导问题[3] ，但高 Rayleigh 数(Ra)条件下， 解未知分布函数的方程组，增加了计算量[10] ；外推方法 由于现象的复杂性，用传统的方法并不能很好解决，因 的数值稳定性较差[11] ；非平衡态外推具有二阶精度，数 此一些学者开始用热格子 Boltzmann 模型来模拟空腔内 值稳定性好，且计算简单、容易实现[12]. [4] 自然对流的速度和温度场分布 . 热格子 Boltzmann 模型可分为两种