多孔介质中变粘性Benard对流的格子Boltzmann模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083202 多孔介质中变粘性Benard 对流的格子 Boltzmann 模拟 1 2 1 荣伏梅 ，施保昌 ，郭照立 1 2 ( 华中科技大学煤燃烧国家重点实验室，华中科技大学数学与统计学院，武汉 430074) (Tel: 027Email: rongfm@126.com) 摘要：多孔介质中的流动与传热涉及许多工程应用，是当前工程热物理学科最为活跃的前沿研究领域 之一。本文用格子Boltzmann 方法模拟了填充多孔介质的方腔内的变粘性Benard 对流，分析了Pr 数的 变化和粘性随温度变化对流场和温度场的影响，进一步给出了相应的Nusselt 数。研究结果揭示了变粘 性情形下Benard 对流的一些特征。 关键词：Benard 对流、 变粘性、 多孔介质、 格子Boltzmann 方法 0 引言 两平行光滑平板之间充满 Boussinesq 流体，从底部加热而产生的对流称为 Benard 对流。这个模型最早是由Rayleigh （1916年）提出的，后来为许多人接受。Benard 对流 是研究稳定性和非线性动力学特性的经典问题之一，在工程上也有许多重要的应用，如 反应堆，石油开采，热交换等。在多孔介质条件下，Benard 对流的数值模拟已成为很多 学者研究的课题。许多文献研究的是粘性为常数的情形，然而温度的变化对粘性有强烈 [1] 的影响，比如水当温度从10度上升到50度粘性下降2.5倍 ，流体粘性随温度变化的 性质将直接改变常粘性情形下的流体传质传热性质。基于这一事实，我们对变粘性情形 下充满多孔介质的Benard 对流做了一些初步的研究。 1问题描述及控制方程 设区域为正方形的腔体充满各向同性，均质的饱和多孔介质 （如图 1 所示），左右 壁面绝热，上下壁面温度分别为T ，T （其中T T ）。当T ，T 的温度差不是太大 c h h c c h 时，流体只有热传导而无对流；当T ，T 的温差逐渐增大，差值达到某一值时，流体会 c h 产生对流，这样的问题就是 Rayleigh-Benard 对流问题。对于腔体充满多孔介质的二维 Rayleigh-Benard 对流问题，其无量纲形式的数学模型为： 资助项目：国家自然科学基金 （Nos.； 国家重大基础专项研究基金 （No. 2006CB705800 ） u 0 （1） u u   2 u  (p )   u+F （2 ）   e t    2  u   