利用时间推进法与QUICK格式研究矩形方腔自然对流.pdfVIP

中围T程热物理学会 传热传赝学 第卜一厢年会 编号：053180 方腔自然对流 陈晓春。， 朱颖心。，王元。 (①清华大学建筑技术科学系．北京100084．②西安交通大学流体工程系，西安710049) 法，对矩形方腔自然对流问题进行了数值模拟研究。计算结果表明：随着时间的推进，流场逐渐稳定 下米．最后趋于定常；将最终稳定的速度场与实验和大涡模拟结果进行比较，发现数值模拟能准确的 预测Ⅲ速度场、温度场以及揣动能场。 关键词：矩形方腔自然对流低雷诺数湍流模型数值模拟时间推进法 1． 引言 封闭腔体内的湍流自然对流是流动与传热研究领域中的一个基础问题，它在工程上 有着十分广泛的应1L|j背景，如太刚能集热器、地板辐射采暖、双层玻璃窗、建筑物中的 能量蓄积和电子冷却装置等等，所有这些工程上的应用都显示出对其进行深入研究的必 要性。 近年来，对于封闭方腔的湍流自然对流研究主要是采用低雷诺数湍流模型进行【1“， 尽管有些学者提出了自然对流边界层中的速度分布律，但使用壁面函数像通常湍流场计 算那样求解有限空间自然对流并无多少成功的实例【”。在对流项差分格式上，大多数学 scheme)【124】，主要考虑到混和格式是绝对 者则采用具有一阶精度的混和格式(hvbrid 稳定的n其在计算过程中收敛特性较好。然而实际计算中，采用这种方式处理，获得 的结果并不能令人满意，对流项差分格式精度不够是造成这种结果的一个重要原因。 本文利用Lam．Bremhorst低雷诺数湍流模型对矩形方腔的自然对流问题进行计算， 在对流项差分格式上采用了0u；ck格式，为了获得收敛结果，在计算方法上采用了时间 推进法。从整个计算过程看，采用上述处理方式获得相当满意的结果。 2．数值方法 2．1控制方程 如图1所示的矩形方腔，重力作用方向为一y，由于方腔的两个侧壁面存在温差， 行计算，控制方程可写为： 连续性方程： a(pu。1／苏，=0 基金项m剧家自然科学基金硕Ef资助(NN 152 ㈤ 言(雕)+毒(以¨一考+毒卜考一p瓦j+助(丁一兀k 能量方程： 塑地塑：旦f—生翌一万1 ㈩ i 8t j瓠|。x?＼、pPr白xj3 一p弼=“[考+善J-詈p螈 (4) 一p可=鲁篆 (5) k方程： O西k+u．酝c3k=毒[(y+毒]善]+最+G；一占 (6) s方程： --如+ui (7) tOx．丝o 良．Il 盯．J缸．1 +舶+刳考l 湍流粘性系数采_I{j如r公式计算： =一T VI=C。九kT (8) 方稗中阻尼函数f，『，厶干¨，：，分别表示如r【6】： 153 ‘=[】一exp(一0．0165Re。)]2(1+20，5／Re，)