低雷诺数下矩形竖腔内混合对流的格子Boltzmann方法模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083165 低雷诺数下矩形竖腔内混合对流的格子 Boltzmann方法模拟 郭亚丽 沈胜强 李维仲 （大连理工大学海洋能源利用与节能教育部重点实验室，大连116024） （Tel: Email: ylguo@dlut.edu.cn） 摘要 文中采用格子Boltzmann力矩模型模拟了不同纵横比 （A）、瑞利数 （Ra）、雷诺数 （Re）下两维矩 形腔内的混合流动。用Richardson数 （Ri）衡量自然对流和强制对流模式对换热影响的相对程度。通过对 模拟结果的分析，表明在低雷诺数下，当流动是强制对流占主导时，换热可看作矩形腔内上部强制对流换 热Nu cv 加权该部分尺寸与下部导热Nu cd 加权下部导热尺寸的综合作用结果，得到整体换热强度随纵横 比变化的经验公式，并给出其适用范围。 关键词 格子Boltzmann方法；混合对流；纵横比；矩形腔 0 前言 近些年来格子Boltzmann方法 （LBM）因其程序设计简单、易处理复杂边界问题、适于 大规模并行计算等特点，受到越来越多的关注。该方法基于对流场的介观层次的描述，它将 流场离散为规则的格子，假设流体仅沿格线运动，并在到达格点后与当地粒子相互作用后继 续运动，从而完成质量、动量及能量的输运。它已被广泛地用于解决多相流、多组分流、带 相变的传热问题等多种领域的研究。目前LBM对换热问题的模拟主要有多速模型、双分布函 数模型和有限差分模型。本文采用双分布函数模型，为了节省计算时间，动量方程采用D2Q9 模型而能量方程采用D2Q5模型进行计算。 顶盖驱动的方腔流常常用来作为验证数值方法可靠性的计算对象。在顶盖驱动的方腔流 中，若腔壁存在温差，由于粘性的存在，流体受浮升力、惯性力与粘性力的共同作用，流动 变得更为复杂。不少文献研究了方腔内的流动，但对具有传热的顶盖驱动的方腔流所作的研 [1-2] 究很少。同样，对矩形腔内流动的研究也较少，一些文献 数值模拟了纵横比不为 1 时的 矩形腔内的顶盖驱动流，考察了不同雷诺数下矩形腔内涡旋的变化情况，依然没有涉及到同 时存在换热时流动与换热相互作用的分析。一般而言，在矩形腔尤其是纵横比 （高度与宽度 之比 A）较大 （A2）的深矩形腔内，由于向深度方向的扰动较慢且所需的格点数较多，若 没有有效的数值方法，实现准确模拟非常困难，如文献[3]所言，即使对Stokes流动而言， 数值模拟深矩形腔内的流动也很难。非绝热条件下由于温差产生浮升力使流体运动受到影 响，流体速度反过来又会影响传热过程，这使模拟的难度进一步增加。 本文采用格子Boltzmann力矩模型，模拟了具有不同纵横比的矩形腔内的混合对流。比 较了不同雷诺数、不同瑞利数及不同纵横比下的流动和换热特点，分析了Richardson数对 混合对流的影响，得出低雷诺数下纵横比对平均Nusselt数影响的分析拟合公式及其适用范 围。 1 Lattice Boltzmann方法 1.1 D2Q9力矩模型 本文采用的D2Q9模型基于向量空间和力矩空间的结合，即用9维向量定义流体在x 位置t时刻的状态，则可构建基于离散速度{c }的N 9 维向量空间F R 9N ，同时构建基 i 于力矩{f }的标量空间 R 9N ，其中Ｎ为空间节点数。显然，在每一格子节点处相应离 i 散速度有9个独立的力矩。力矩空间与向量空间的关系如下[4] m a f （1） j ij i