圆形与圆柱形热源羽流计算方法的研究.pdfVIP

圆形和圆柱形热源羽流计算方法的研究 西安建筑科技大学张强朱树园周宇 摘 要 回顾与分析了圆形和圆柱形热源羽流流量、轴心速度：和轴心过余温度的三种预测 方法：虚拟点源法、积分法和数值模拟法．数值模拟时分别采用了四种k—g模型和MIT 零方程模型作为湍流模型，还讨论了不同的壁面函数和近壁处网格划分方法对模拟结果的 影响．模拟结果与虚拟点源法和积分法的半经验公式计算结果进行了比较，结果表明：三 种预测方法的准确程度大体上随离开热源的高度增加而改善．在离热源较远处三种方法的 结果基本一致，都较可靠．虚拟点源法和积分法均不能正确地预测靠近热源处轴心速度的 变化．数值模拟法如果采用合适的湍流模型和网格系统则可以获得较合理的预测结果，但 有些湍流模型如MIT零方程模型的偏差太大，不宜采用．在近壁处进行网格细分的情况下， 标准k一￡模型结合加强壁面函数法大体可以满足工程计算的要求． 关键词 自然对流羽流数值模拟湍流模型 1引言 热源上方因自然对流引起的空气流动称为热羽流或浮力羽流，简称羽流。羽流现象 是通风工程中经常遇到的现象，尤其在自然通风和置换通风中，它是空气流动的唯一或 主要动力。掌握羽流产生发展规律是通风工程基础研究和技术进步的核心问题之一。早 在上世纪五、六十年代前苏联M．A．谢比列夫n’21等许多学者就对羽流及其引起的室内热 分层现象进行了深入的理论分析和大量的实验研究工作，极大地推动了工业建筑自然通 系统而持续的研究，从理论和实验上深入揭示了羽流现象的物理及几何特征。七十年代 后期由北欧首先研发后来得以迅速推广的置换通风技术，以及由此形成的“北欧经验” 刚，既得益于这些基础研究的成果，也大大推动了羽流及其应用研究的开展。然而由于 实际热源上的湍流羽流流场和温度场无解析解，至今其计算都采用经验或半经验的公 式，因实验条件和拟合的方法不同，不同研究者得出的公式也不一致。另一方面，随着 计算流体力学(CFD)和数值传热学的迅速发展，数值模拟已经是室内气流组织研究的 主要手段之一，用数值模拟方法确定羽流流场和温度场是另一种选择。然而在数值模拟 时，首先要遇到湍流模型的选择问题。本文将针对圆形和圆柱形热源上的羽流计算问题 首先讨论现有经验和半经验公式，再比较不同湍流模型的数值计算结果，力求找出合适 的预测方法。 2点源羽流流量公式的比较 现有实际热源羽流的经验或半经验公式都是建立在点源公式的基础之上的，所以有 必要先讨论一下点源的计算公式。 应用热射流理论可以求得均一温度环境中点热源上湍流羽流流速分布、流量和过余 温度分布的半解析解。之所以称为半解析解足因为求解中引进了某些假定且公式中包含 必须由实验确定的常数。不同文献中实验常数的取值不尽相同导致推荐的公式(主要是 系数)略有差异。目前应用较多的点源羽流流量计算公式见表1。 表1公式中： qV．z——某一高度羽流体积流量 西一点热源对流换热量 2007年全国通风技术交流大会论文集 乒—计算断面至点源的垂直距离 正——羽流卷吸常量，仉=O．118 争—-羽流浮力通量，B=@g／(c．o刀 伊—-重力加速度 伊—·空气比定压热容 P-—咤气密度 产—空气热力学温度 为了对表1公式计算结果进行比较，取a=50W，环境温度300K计算羽流流量，结 果如图1所示。 表1点源羽流体积流量的计算公式 发表者(发表年) 计算公式 吼工：0．005中坍少 妣erzwinski和Popiolek瑚(1981) 图1五种点源羽流流量计算公式的比较 从图l可以看出随着高度的增加羽流流量增大，各公式计算差异也增大。超过2米后最 大(ASHRAE 1980推荐公式)最小(王炎生公式)值偏差约达10096。我国暖通专业的统 编教材工业通风》早期版本曾推荐ASHRAE1980公式，新版本则推荐王炎生公式Ⅱ3】， 显然是注意到了这一差别