065通风空调百叶风口阻力特性与几何结构的关系正文.docx

通风空调百叶风口阻力特性与几何结构的关系西安建筑科技大学李丹李安桂摘要 由于百叶风口结构简单、安装方便等特点在通风空调工程中得到了广泛使用。但是，工程实践中发现：百叶倾斜角度，百叶宽度、厚度，百叶垂直间距甚至百叶本身的结构等对通风 效果、压力损失、阻力系数甚至噪声效果有重要的影响。本文基于理论分析，以数值研究方法对单侧百叶风口的阻力系数与几何结构的关系进行了探讨。数值模拟采用FLUENT 6.0软件，计算结果用Excel进行后处理。研究分析了百叶倾斜角度，百叶挡板宽度，百叶挡板垂直间距 及百叶挡板厚度对风口阻力系数的影响。关键词百叶风口数值模拟局部阻力系数Numerical Simulation of Shutter Resistance Character Used in Air-conditioning and VentilationBy Li Dan Li AnguiAbstract Based on the theoretical analyses, describe local resistance coefficient of the outlets through numerical simulation. Use the software of FLUENT 6.0 for numerical simulation and use Excel for processing the results. Analyze the different results of resistance coefficient brought by changing the obliquity of the boards, the width of the boards (b),the ration of width and upright distance between two boards (b/h) and the thickness (d) of boards.Keywords shutter, numerical simulation, local resistance coefficient 1引言侧送风方式具有布置简单，施工方便，能满足房间对射流扩散、速度衰减的要求，因而被 广泛用于空调房间送风。在侧送风方式中，百叶风口由于结构简单而被普遍应用，尤其是在自 然通风中，百叶风口更为常见。百叶风口的应用已有悠久的历史，但由于其结构的可变化范围 很大，实际上，工程经验表明人们至今对其和气流组织的关系尚未有清楚的认识。鉴于此，本 研究对单侧百叶风口的性能进行数值模拟，以便更深入地了解这一传统的风口动力特性，从而 在建筑空调的使用中更好地发挥其性能优势。2百叶风口的数值模拟模型的建立衡量风口性能优劣主要依据是在风压转换系数 K值相同的条件下，局部阻力系数小者为优，大者 为劣1，在百叶风口中，决定局部阻力系数的因素 主要有百叶挡板倾斜角度、百叶挡板宽度b、百叶挡板宽度与百叶挡板垂直间距的比值b/h及百叶挡板厚度di，因此，本研究从这几方面入手，进行 数值模拟。关于风口模型的建立，采用了以下几点假定 4 5 6 7:风口内空气流动为不可压缩气体紊流流动；由于本研究只关心百叶间的空气流动状况，且百叶风口在宽度方向（ X方向）对称，因此 将其简化为只具有厚度（Z方向）和高度（Y方向）的两维模型（即认为X方向无穷大）， 相当于用一个矩形截面来截此型体，从而大大减少了计算量，百叶风口的模型见图1 ；工程中使用标准的k-£方程已有足够的精度，本 文以稳态标准k-£方程来描述百叶风口内的空气 流动；由于百叶风口相比其他建筑空间较小，如果直接以百叶风口的进出口侧为进出口边界，可 能由于气流的未充分发展而使计算产生较大的误差，故本文通过对多个模型的模拟计算， 找出了两个假定的边界，距百叶的进口侧1000 mm，距出口侧1500 mm，在这两条边界处阻力系数的计算已基本不受所设边界距离的影响；数值计算选用控制容积法对控制方程的微分形式进行离散，并采用压力校正的SIMPLEC算法；对动量方程，紊流脉动动能方程（K方程）及紊流动能耗散率方程（£方程）采用QUICK 算法，对墙壁及百叶挡板均采用标准壁面函数法及无滑移壁面条件；百叶挡板的倾斜角度变化范围为 0°60 °，百叶挡板宽度b的变化范围为10 mm100 mm、 百叶挡板宽度与百叶挡板垂直间距的比值 b/h的变化范围为0.43.0，百叶挡板厚度d的 变化范围为0.6 mm2.1 mm；3数值模拟结果及分析1）阻力系数修正值a的确定本文模拟了当百叶入口侧风速从 0.5m/s变化到10m/s过程中的流动情况，计算了每种情况1)2)3)4)5)6)图1 百叶模型下流体流经风口时的风速v，根据公式a =虬，求得平均值作为a的确定值，其中，流过风口 V' 2时的v由V = 2；1 V计算，vi是第i个和第i+1个挡板间流体的速度，n为百叶风口的挡板”i=i数量。阻力系数数值计算的其他条件为：挡板 b/h=1.0，板长b=50 mm，板厚d=1.5 mm，空气密度取1.2 kg/m3，挡板倾斜角度为30。（假定挡板水平放置时为0°）。经过计算得出，工程上所应用的阻力系数值实际上比阻力系数的准确计算值大，二者的关一， 匚系为：°准确 =0.813。r工程2）百叶风口用于自然通风时，室外风速的变化对阻力系数的影响1933年，尼古拉兹对于六种不同人工粗糙度的圆管子的流动状况进行了大量的实验研究，结果将流动分成五个区域，即层流区域、紊流发生区、紊流光滑区、紊流过渡区和阻力平方区（自模区），得出当流动进入阻力平方区时，沿程阻力系数入仅与相对粗糙度K/d有关而与Re数无关。对于气流流过百叶这种板间流动，上述理论内容依然成立，本研究将通过数值模拟，确定流体进入阻力平方区时的风速。计算所得的阻力系数计算表如下所示：0.95 0.511.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10外界风速(m/s)图2外界风速变化对阻力系数的影响曲线由图2可以得出，百叶风口应用于自然通风时，当室外风速达到7m/s时，空气的紊流流动达到阻力平方区，此时阻力系数不再随外界风速的变化而改变。3）挡板倾斜角度的变化对阻力系数的影响阻力系数数值计算的其他条件为：挡板b/h=1.0，板长b=50 mm，板厚d=1.5 mm，空气密度取1.2 kg/m3,设定来流速度速为3m/s。5.9如心5.：3；沁 DD4.?8e+00i.l'jc+uO3.59.1+UU£.99.i+i：i0£.39.i+00l.Suc+uOLSu-i+uD6.i：n：i.i-U 11.70.-03S.U5.2+IJU?.63t+0i：i6.7沁 Mi6.36.i+i：n：i5.94c*005.51.1+LiU5.09.2*i：h：i4.24t+0i：i3.82.1+Oli3.40.2+002.9?i+Uii2.55c*002.13.1+Oli1.70.2+001.沧+口 0S.5S.-U 14.：Z4.2-i：H（c） A处放大图（a）挡板倾斜角度为30 °（ b）挡板倾斜角度为60 °图3不同挡板倾斜角度下的速度场通过图3、图4可看出：1）挡板倾斜角度越小则阻力损失越小。 阻挡来流的叶片后面，必然有涡存在，涡区 是低压和低速区，叶片倾斜角度大，其后面 的涡区域就较大，它的存在必然使叶片后形 成负压区，叶片前后形成压差，带来百叶风 口的阻力损失。并且挡板角度越小，气流遇 到的挡板阻挡越小，气流撞击挡板所形成的 撞击损失越小。虽然 0。倾角的挡板的阻力 系数最小，但其水平放置的挡板不具任何导 流作用，其作用与一空洞相似，故工程中一 般不予采用；大 国家自然科学基金资助项目（No.50278075,o o o o oL值确准数系力阻350478046)百叶挡板倾斜角度图4阻力系数随百叶挡板倾斜角度的变化曲2）随挡板倾斜角度的增大，挡板后涡区增大，并且，当挡板倾斜角度增大到60。时，挡板后形成了两个涡区，见图 3 （c），原因是当气流受到挡板阻挡，在流过挡板间隙时，速度增大， 如果挡板的倾斜角度过大，流体就不能很好的沿挡板倾斜方向流动，而与挡板脱离，使得挡板的导流效果减弱，两股高速气流互相干扰.于是，在挡板下又形成了一个涡流区，增加了涡流损 失（因形成涡流而消耗的能量）。3）挡板b/h的变化b与h如图5所示。作者依次模拟了 b/h值从0.43.0的变化。图5百叶的b和h图6阻力系数随百叶挡板b/h值的变化曲 线模拟条件为：来流速度为3 m/s，挡板倾斜角度为30 °，且板长b=50 mm，板厚d=1.5 mm， 空气密度取1.2 kg/m3；具体情况如图6和图7所示：通过模拟可以得出：1）随b/h值的增加，压力损失大体上逐渐增加。因为在板长b 一定的情况下，b/h越大，板间距h越小，由于它减小了气流弯曲半径，气流经过百叶挡板间隙绕过百叶挡板过程中由 于离心力的作用使气流弯曲所造成的阻力损失增大。同时通过图7可看到在b/h=1.2时，阻力系数出现了一个极小值，这和挡板下方形成的涡流区有关。因为挡板b/h值越小，涡流区越明显，当挡板b/h=0.4时，挡板下方出现了两个涡流区，同样也是受到挡板上方和 下方两股高速气流的影响，涡流的形成影响出口速度并增加了涡流损失，当b/h=1.2时，涡流区开始明显减弱，此时压力损失小，当b/h=2.2时，百叶风口下部挡板的涡流区开始消失（图7）；2）鉴于压力损失、出口速度及所形成的涡流区大小的综合考虑，作者推荐使用0.8 Wb/h W 2.0的百叶风口，鉴于在b/h=1.2时阻力系数最小，建议广泛使用。5.06«+004.56t+00：3.55.i*Oi：i：3.：3i：i1；+i：ii：i2.79.i*i：n：i2.54.i+i：i02.29c+001.7,9c*i：n：i1.53-i+i：n：i7.78c-015.27-i-Ol2.75-i-Ol5.91C+005.6Ut+i：iU5.29i+UU4.9如口口4.67t+i：iu4.35.2+UU4.u4t+i：iU3.73t+i：n：i3.4£t+i：ii：i3.11t+i：n：i2.80c+002.49c+001.87t+i：iU1.56t+UU9.招如口 16.22.-U 13.11.-U16.45t+ij06.11t+ijU5.77.;+i：ii：i543t+i：iO5.10t+i：i04.76.1+LiO4.42t+i：i04.U：:；-i+UU3.74.i+i：ii：i3.4ij.i+iju3.0600£.7£.i+i：i02.3Sc*002.u4.i+0u1.7i：i.i+i：i0L：3E，5Wi1.02006.804-01：=：=：r.-._ri.<大 国家自然科学基金资助项目（No.50278075，50478046）(a) b/h=0.4(b) b/h =1.2(c) b/h =2.2(f) C处的局部放大图(d) A处的局部放大图图7(e) B处的局部放大图不同挡板b/h值时的速度场对比3)挡板b值的变化作者依次模拟了 b从10 mm到拟条件为：挡板倾角为 30°，挡板b/h=1.0，板厚 d=1.5 mm,空气密度取1.2 kg/m3,具体情况见和图 8、9所示，可得出：随着百叶b值的增大，压力损失先减小后 增加，在b=40 mm时最小。因为流体经过的 挡板长度增大，为克服挡板表面摩擦力所 耗的能量就增多，压力损失将增大，但对 b/h值同为1.0的百叶，板间距将增大，气100 mm的变化，模1)图8阻力系数随b值的变化曲线2)3)流经过百叶挡板间隙绕过百叶挡板过程中由于离心力的作用使气流弯曲所造成的阻力 损失减小，二者的综合作用产生了一个板宽的最佳值，此时流体的压力损失最小； 由图10看出，挡板下面的涡流区随板宽的增加而增大，增大了涡流损失； 由表中数据可见，综合考虑到涡流区等的影响，建议选择 b=40 mm时具有最小的阻力系数值，建议普遍使用。10 mmW bW 50 mm的板宽，当(c) b=100 mm(a) b=10 mm?： S7.-H1(d) A处放大图(f) C处放大图图10 阻力系数随d值的变化曲线(e) B处放大图图9不同挡板b值时的速度场对比4) 挡板厚度d值的变化由于工程上普遍使用铝合金百叶风口，本文 模拟风口挡板厚度为0.6 mm2.1 mm,模拟条件 为：来流速度为 3m/s，挡板倾角为30。，挡板 b/h=1.0，板宽 b=50 mm,空气密度取 1.2 kg/m3；从图10可见，百叶挡板厚度对百叶风口的局 部阻力系数影响不大，当d=1.5 mm时局部阻力系 数最小，建议使用此值。5) 改变百叶板入口侧形状，如图11所示，由于 其良好的导流作用，将会减少由于板边缘垂直变 化引起涡流。4结论1)根据公式AP =。皓的物理意义，Z的准确计算值与工程应用值并不相等，二者的关系为匚刀 |> = 0.813 ；'工程2)根据紊流分区的理论，百叶风口用于自然通风时，当室外风速为7m/s时，空气流动进入阻力平方区(自模区)，百叶风口的阻力系数为一定值，不再随室外风速的改变而变化；3) 挡板倾斜角度越小，压力损失越小，因此在播音室、图书馆等要求进口风速比较小的场所， 宜采用挡板角度W 30。的百叶风口；在电影院、通用办公室等要求进口风速比较大的场所， 百叶挡板的角度应选W 45。，并可根据实际需要调节百叶倾斜角度来调节进出口风速；4) 百叶挡板的b/h值越大，压力损失越大，鉴于压力损失、板间形成涡流区的大小，建议百 叶挡板的b/h值取0.8到2.0之间，b/h=1.2为最佳值；5) 当百叶挡板宽度b=40 mm时，百叶挡板的阻力系数最小，考虑到涡流区等的影响，建议选 用10 mmW b W 50 mm的板宽；6) 百叶挡板的厚度可取为1.5 mm；7) 叶挡板的入口侧可选用具有倾斜度的楔形挡板来代替矩形挡板，利用其较好的导流作用， 可得到更好的出流效果。大 致谢：本研究得到国家自然科学基金资助(No.50278075，50478046 )。参考文献1 郭春信.地下空间自然通风.北京：中国建筑工业出版社，19932 陶文铨.数值传热学.西安：西安交通大学出版社，19883 卜增文.空调末端设备安装图集.北京：中国建筑工业出版社，20004 张祝新，程晓新.对工程流体力学节流孔口流量特性的讨论.机床与液压，1999，6: 68-695 邢春礼，秦裕琨.百叶窗煤粉浓缩器阻力特性研究.热能动力工程，1994，9(4)： 221-2276 E.P.,D.Sc Vishnevskii, Petrosoek Climate Maker Co., St. Pete rsburg, Russia,Aerodynamic Properties of Slot Outlets and Exhaust Hoods Engineering, 5thInternational Conference on Air Distribution in Rooms, ROOMVENT 96, July 17-19, 19967 赵彬，李先庭，彦启森.百叶风口送风射流的数值模拟.暖通空调，2001，31 (6)： 86-898 王致清，等.流体力学基础.北京：高等教育出版社，1991