双扩散自然对流的格子Boltzmann模拟-中国科技论文在线.PDF

中国科技论文在线 第 19 卷 第 6 期 计 　　算 　　物 　　理 Vol. 19 ,No. 6 　 2002 年 11 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Nov. , 2002 ( ) [文章编号] 　1001246X 2002 双扩散 自然对流的格子 Boltzmann 模拟 郭照立 , 　李 　青 , 　郑楚光 (华中科技大学煤燃烧国家重点实验室 ,湖北 武汉 　430074) [摘 　要] 　建立了一个模拟双扩散 自然对流系统的格子 Boltzmann 模型 ,用此模型对受温度和浓度梯度驱动的方腔 流动进行了模拟 ,研究了浓度 Rayleigh 数对传热传质的影响 ,并与其它文献的结果进行了比较. [ 关键词] 　格子 Boltzmann 方法 ;双扩散对流 ;传热传质 [ 中图分类号] 　O35 [文献标识码] 　A 即多速模型和多分布模型. 前者是等温模型的推广 , 0 　引言 只用同一类分布函数 ,并用多个粒子速度. 此类模型 当溶质进入具有不同温度的溶液中时 , 由于存 的平衡态分布函数包含速度的高阶项 ,流体的温度 在浓度和温度梯度 ,混合流体中会有热量和质量的 由分布函数的二阶矩得到. 典型的多速模型有二维 [6 ] [7 ] 扩散 ,进而发生自然对流现象. 这种由浓度梯度和温 13 速模型 、三维 21 速和 25 速模型 等. 多速模型 度梯度的综合效应引起的自然对流称为双扩散 自然 的主要不足是 Prandtl 数是常数 , 且数值稳定性较 对流[1 ] . 双扩散 自然对流是一种常见的现象 ,在海洋 差. 此后人们针对这些 问题提 出了若干 改进方 学 、大气物理学 、地球物理学 、生物流体力学以及化 法[8～10 ] . 多分布模型的出发点是 ,如果粘性热耗散和 学动力学等领域中都会遇到这种流动. 由于热量和 压力功可以忽略 ,则温度可以看作是一个随流体运 质量的扩散速度一般不同 ,双扩散对流比单纯由温 动的被动标量 ,其运动可以用一个对流扩散方程描 度梯度引起的自然对流要复杂得多. 虽然人们很早 述 ,故可用一个独立于密度分布函数的新分布函数 就观察到了这种现象 ,但双扩散对流体运动的影响 来模拟温度场. 所以,在多分布模型中有两类分布函 机制至今仍然不很清楚. 已有研究表明 ,双扩散对流 数分别模拟速度场和温度场. 多分布模型的优点是 系统中包含许多令人感兴趣的非线性现象. 由于双 具有良好的数值稳定性. 典型的多分布模型是 1997 扩散对流的复杂性 ,理论分析方法受到很大的限制. 年 Shan 提出的基于多组分LB 方法的模型[11 ] . 随着计算机技术和数值计算方法的发展 ,数值模拟 上述热 LB 模型主要用于模拟单一流体 的流 方法已经成为研究此类问题的主要手段[2～4 ] .