268 送风孔板阻力系数模拟研究.doc

送风孔板阻力系数模拟研究 同济大学 朱 春 张 旭摘要 孔板作为一种风口形式，不仅可以用于洁净室送风末端，也常常用于地铁站台通风、列车车厢送风等人员密集或是空间相对狭小等场合。本文首先理论分析影响孔板阻力特性的参数，接着利用CFD模拟了不同开孔率下孔板的阻力特性，得出了孔板阻力系数的计算关联式，以便于工程实际设计与计算。关键词 孔板 开孔率 阻力特性 0引言 孔板作为一种风口形式，不仅可以用于洁净室送风末端，也常常用于地铁站台通风、列车车厢内送风等人员密集或是空间相对狭小等场合，起到均匀送风的作用。对航空航天等某些特殊通风领域，当送风管尺寸不能变化时，孔板可作为末端局部阻力件来调节各送风支路阻力平衡，此时需要已知孔板的阻力系数。常用的通风空调手册中1，没有述及多孔孔板的局部阻力计算，因此有必要对孔板的阻力特性进行研究。孔板风口一般为平板上均匀开设多个圆孔，风流通过孔板时形成多股小孔射流，然后混合汇聚，均匀送出。为计算孔板的阻力，需要给出每一个小孔的入流条件，分析各个小孔的流动特征，这给理论计算带来许多麻烦和不确定性，因此必须简化模型。简单的办法是将孔板风口等效为一个简单开口，其面积与孔板风口的有效通过面积相等，这样可以确保入流的动量流量和质量流量与实际一致2，简化模型如图1。根据分析，把气流流过孔板过程分解为突缩和突扩两个过程。对于突扩过程，由2-2断面、3-3断面应用贝努利方程，并代入动量方程，可得鲍尔德-卡尔诺公式： ( 1 ) Fc为自由出流收缩断面2-2的面积，引入收缩系数= Fc / F2，和开孔率（也称为扩张比）n=F2 / F3，公式（2）可写为 （2）当Re 104时，收缩系数可认为定值。图1 简化的孔板出风示意图对于突缩过程，由于收缩过程产生漩涡的复杂性，只给出经验公式3： （3）因此，气流通过孔板的阻力系数为： （4）理论分析中，只考虑了突扩、突缩过程的“冲击”损失，忽略了两个过程的沿程阻力损失，仅适用于Re104的紊流情况。简化的模型中，同时忽略了多孔射流的混合损失，而该损失与开孔孔径和开孔率有关。本文通过CFD模拟方法，分别研究影响多孔孔板阻力系数的因素，包括：孔板厚度l, 小孔孔径d，孔板开孔率n和孔板入流平均Re数的函数。2 模拟计算为分析以上各因素对孔板阻力特性的影响，研究思路为只改变单个影响参数，其它参数取为定值，采用CFD软件模拟方法，分别对各因素进行分析研究。建立如图2的四方体几何模型，模型尺寸为0.30.111m，模型中间截面设孔板，上、下底面分别设送风、出风口，两风口与四方体截面等尺寸，即为0.30.11m。四方体各侧面默认为对称边界条件。根据达西-维伊斯巴赫公式： ( 6 ) 给定送风入流速度，采用k-计算模型，可模拟计算出入口和出口的压力差，因而得出不同几何条件下孔板的阻力系数。图2 孔板几何模型2.1 孔板开孔尺寸d设定入口流速0.202m/s，孔板开孔率n为13.4%，孔板厚度l为3mm，分别取孔板开孔直径为4mm、6mm和8mm，开孔均匀布置。图3给出孔径4的截面模拟结果： 图3 孔径4的压力场和速度场 根据以上CFD计算结果，整理各孔径下的阻力系数见表1。表1 不同孔径的阻力系数（n=13.4%）开孔直径d压力差p（Pa）阻力系数4mm2.901186mm3.141288 mm2.99122由表1可见，相同开孔率条件下，不同孔径的阻力系数最大相差7.5%，计算阻力系数与孔板孔径大小不呈同一趋势变化关系，因此可认为孔板阻力与开孔大小无关。2.2 孔板厚度l设定入口流速0.202m/s，孔板开孔率13.8%，开孔直径6mm，分别取不同开孔厚度，开孔均匀布置。5mm厚孔板模拟结果见图4。 图4 5mm厚孔板压力场和速度场根据CFD处理结果，计算各孔板厚度下的阻力系数见图5。图5 不同孔板厚度的阻力系数（n=13.8%）由图5，相同开孔率条件下，1mm与10mm孔板厚度的阻力系数相差8.2%，15mm孔板厚度的阻力系数相差3.5%，可认为孔板阻力与孔板厚度无关，。2.3 孔板开孔率为研究开孔率与阻力系数的关系，设定入口流速0.202m/s，取孔板厚度3mm，开孔直径6mm，开孔均匀布置。模拟结果见图6，图7为n取10%40%的部分结果。图6 n-曲线关系 图7 n取10%40%的阻力系数根据图6、图7结果，孔板开孔率与阻力系数表现出明显的基本函数关系，根据曲线回归分析，最后拟合n-的幂函数关系如下： R2 = 0.9931 ( 7 )为工程简化计算，可如下拟合公式： R2 = 0.9885 （8）2.4 Re数取孔板开孔率13.8%，开孔直径6mm，孔板厚度3mm，设定不同速度，模拟Re数与关系如表2。（其中运动粘性系数取20空气值：15.710-6m2/s）入口速度v1Re103阻力系数0.2020.2431272.022.4312320.224.3120可见Re在25024000内，阻力系数相差5.5%。认为Re30时，射流进入湍流区4，阻力系数与Re数无关。3 结果检验 为检验模拟结果的可靠性，根据文献4的圆管内设多孔孔板的一组实验数据，与模拟公式（7）曲线比较结果如图7，计算R2=0.9738，可见模拟曲线与试验结果强相关。 图7 n-模拟曲线与试验结果比较4 结论综合上述分析，当Re30，孔板阻力特性仅与开孔率有关，与孔板厚度、孔板孔径大小及Re数均无关。在开孔均匀布置情况下，开孔率与阻力系数（- n）关系式为； R2 = 0.9931上式中开孔率n取值范围为150%。参考文献1 陆耀庆实用供热空调设计手册M北京：中国建筑工业出版社，19932 Nielsen P VDescription of supply openings in numerical models for room air distributionJ ASHRAE Transactions1992，98(1)：963-971.3 华绍曾，杨学宁 等编译实用流体阻力手册M北京：国防工业出版社，19854 余常昭紊动射流M北京：高等教育出版社，19935 赵彬，李先庭，