等截面水泥烟道沿程阻力研究

1、工程建筑论文-等截面水泥烟道沿程阻力硏究摘耍：通过对等截面水泥烟道模型沿程阻力的实验硏究得到烟道内壁绝对粗糙 度k値，从而为烟道系统设计提供依据。运用fluent气流模拟，提出了烟道系 统沿程阻力估算公式，通过与理论计算结果的比对，证明了其在实践中应用的可 行性。关键词：住宅厨房烟道沿程阻力当量糙粒高度气流模拟沿程阻力估算目前我国大多数低、中、高层住宅的厨房采用了集中排烟系统，最高应用层数已 达48层，正在设计的有50层。集中排烟系统的流动主要由以下三个部分组成： 烟道内壁粗糙形成的烟气流动沿程阻力；主烟道与用户支管的烟气汇合产生的合 流阻力以及用户止逆阀玄因道排风帽等设备造成的局部阻力。对于

2、高层住宅而言， 数十米甚至超过一百米的内壁粗糙的烟道会造成很大的流动阻力，要求排油烟机 需具备较大的出口全压。由于烟道系统的设计通常采用估算法，即根据使用条件 及设计参数假设一个烟道的截面尺寸，然后过计算校该其尺寸是否能满足要求， 因此设计时合理地选用k値（管道内壁当量糙粒高度）成为设计有效烟道截面尺 寸的必要条件。廿前所使用的烟道一般均为玻纤网增强水泥混凝土管道，内壁的 光洁度较好。现有相关资料给出的混凝土管道内壁的绝对粗糙度一般为13mm， 相差范围较大。以前住宅层数不多时，沿程阻力对排油烟机的影响很小，但目前 高层住宅的层数越来越大，实际使用的烟道内气流流速最高可达十几米/秒，因 此，烟

3、道内壁k値的取値对烟道设计截面积的影响变得十分重要，有必要通过实 验硏究获得比较确切的数据。1等截面烟道沿程阻力系数入及k値的实验硏究烟气在等截面烟道内流动为克服沿程阻丿j引起的能量损失用压强损失表达：a（1）从公式（1）可以看出，沿程阻力计算的核心问题是各种流动条件下无因次系数乂 的取値。尼古拉兹实验比较完整地反映了沿程损失系数入的变化规律，并揭示了 影响a变化的主要因素有雷诺数re与相对粗糙度k/d。本课题最初的实验是采用原型烟道以空气为流动介质进行流动阻力测定，再依据 结果计算出相应的k値，但多次实验的结采均不理想。原因可能有几方面：测定 风量和风压的仪器精度不够，测定误差较大；烟道总长

4、度不够，总沿程损失较小， 使仪器精度误差的影响成为主要因素；风机性能冇限，使管内流动难以保持在阻 力平方区附近等。所以决定改用水为流动介质进行模型实验。1.1实验装置及实验原理实验系统如图1所示，它是由水箱1、出水阀门2 （调节流量用）、受试管段3（外截面尺寸为loommx 100mm，壁厚10mm）、测压管4（读取测定断面的静压差）、 三角堰5 （测量流量用）、回水沟6、水池7及水泵8等组成。根据相似理论，模型实验应与其原型之间有相似性。本实验的目的是测定烟道内 壁的k値。实验采用了实际烟道，所以模型与实型的k値是相等的，对实验结果无须修正。但是水和空气是两种性质的流体，两者的粘滞系数差别很

5、大,只有满足下述情况中的一种时从实验推导山的结果才能被采信:即模型的流动状态处 于不受模型律制约的范围内；或者是两种流动的同名准则数虽然不相等，但可以 认为两者的沿程阻力系数a都与雷诺数rc及管壁相对粗糙度k/d有关。所以， 在实验过程中，保证了管内高水流速度（通过减小烟道截面尺寸、提高水箱位置等实现），以使流动状态接近或进入阻力平方区；并保持测压管水位的高差在10cm以上，以使读数准确。由于水泥烟道是吸水材料，在实验之前对其进行相关处理。图1沿程阻力系数检测系统示意图12实验结果及分析烟道沿程阻力系数入及内壁绝对粗糙度k値的实验测定结果如表1所示。其中， 根据实验数据推导出的入値，k1是采用

6、柯列勃洛克公式（式2）计算而得；k2是 采用阿里特苏里公式（式3）计算得出的；k3则是采用粗糙区的希弗林松公式（式 4）算出来的。平均値kpj是不同情况下k1値的平均。21凶(3)(2)a (4)表1等截面烟道沿程阻力系数测定当量直径(m)长度(m)流量 (1/s)流速 (m/s)沿程 损失 (pa)沿程阻 力系数a雷诺 数rek(mm)kik2k30. 08367.4850. 004210. 680. 0580. 0274434000. 190. 190. 320. 005000.810. 0840. 0282515440.230. 250. 360. 005680. 920. 1060.

7、0276585540. 220.230. 330. 006841. 100. 1530. 0275705120. 230. 250. 330. 007201. 160. 1650. 0268742230.210. 220. 29kpj0. 22表1中的沿程阻力系数是根据公式(1)计算得出得。根据表1中所列不同流速 时的沿程阻力系数入及雷诺数rc，对照莫迪图发现，此次实验过程中烟道内的 水流的流动状态处于紊流过渡区内，因此入値既与烟道内壁绝对粗糙度k値有 关，也与雷诺数re有关。对比ki、k2、k3的计算値可以看出，采用柯列勃洛克 公式算得的k1值相对比较精确;而用粗糙区的希弗林松公式算出的k3

8、値因为没 有考虑re对a的影响，导致结果偏大且误差较大。由表1可知，受测烟道的内壁绝对粗糙度k约为0.22mm o但在烟道设计中还不 能直接采用该値。因为批量生产的烟道的均匀性比较差(日前烟道生产主要采用 手工方式)，而且在施工安装时各节烟道的接口处会有水泥凸出，另外还要考虑 到使用后烟道内的积油污等。所以，笔者认为在系统设计时，1値可根据实际工 程屮烟道产品的选用状况在0.5mm1.5mm的范围内进行选择。烟道产品制作以 及系统施工时，采用精良的制作工艺，减小k値，是减小烟道截面积冋时保证系统有良好运行效果的一个重要前提。有些企业已经开始改用机械化生产，烟道质 量冇望提高。1.3实际设计时沿

9、程阻力系数入的确定当k値确定以后，不同应用条件下的沿程阻力系数入也就可以通过计算确定了， 进而就可以计算沿程阻力。大部分情况下烟道内烟气的流动状态处于紊流光滑区 和紊流过渡区，此时粗糙区的希弗林松公式显然不适用，而柯列勃洛克公式又过 丁繁琐，不适合工程计算使用，那么，选用阿里特苏里公式计算入値是否能够得 到合理的结果呢？表2屮的数据是笔者根据表1屮列出的k1及相应的re采用阿 里特苏里公式反算入（表中定义为;i a）値，并将计算结果与实验所得的入値进 行比较。表2采用不同公式计算入値所得结果比较k10. 190.230. 220.230.21re434005154458554705127422

10、3a0. 02740. 02820. 02760. 02750. 0268人a0. 02740. 02780. 02730. 02720. 0266谋差 e = | x-xa /xx100%01.421.091.090. 75由表2可知，在k値与re値都相等的情况下，分别采用柯列勃洛克公式和阿里 特苏里公式算出的入値其误差不大于1.5% o此外，笔者亦曾用阿氏公式计算过 在不同k/d的情况下，烟道内气体流速较高时的沿程阻力系数入値，并将其与有 关资料提供的相应条件下的入値进行对比，发现两者的误差不大于4% （通常稍 偏小）。所以，烟道系统设计时，在已确定k値与re値的情况下，采用阿里特 苏里公

11、式计算沿程阻力系数入是可行的。2气流模拟分析与沿程阻力简化计算方法通过上一节中对沿程阻力实验结果的分析可知：实际工程烟道设计时取k=1.5， 并由此得出k/d的値，再根据相应的re选取z进行沿程阻力计算是安全的。然而，经由莫迪图查得入値再通过公式(1)计算得到气流流过某一定距离所要克 服的沿程损失的値毕竟是比较麻烦的。所以笔者结合对多个实际烟道系统的测试 结果并采用气流模拟计算后'提出了一种简使的估算方法。表3烟道截面尺寸表(烟道壁厚10mm)适用层数同吋使用系数同吋开机数烟道尺寸(mm)670%5250x30071813400x50019 3060%18500x50031 50306

12、00x600注:表中的同吋开机数是笔者对每一类中楼层数最多的住宅楼按同吋使用系数为60%计算。表3列出的是日前有关设计图集规定的不同层数住宅厨房适用的烟道截面尺寸。 对应每层用户的排风量为3oom3/h。笔者采用fluent软件按照表中的条件分别 进行气流模拟，以期得山一定流速的气流流经单位长度特定截面尺寸的烟道所需 克服的沿程阻力。2.1模型的建立由于实际空气是粘性流体，流动基本为湍流流动，因此fluent计算选择k £两方程模型作为湍流模型。将空气设定为理想气体。由k方程、e方程、动量方程、能量方程、连续方程一起构成了室内空气流动与换热的基本控制方程。其通 用形式为：在fluen

13、t的离散选项中对离散方法做出如下选择：对压力p选用标准方法，压 力速度耦合采用simple算法；湍流动能k和湍流耗散率能量和动量选用一 阶迎风格式；收敛准则：各流动项残差小于10-3，能量项残差小于10-6 ；选用 标准k-6方程。计算时选用六面体结构网格，并假设颗粒均匀分布，即默认c = 0.5取値。气流模拟分別对烟道内壁绝对粗糙度k二1.5和k二1.0两种情况进行了计算。之所 以取k=1.0主耍是考虑到日前的烟道制作加工工艺相对于以前已经有了很大的 改进与提高，在某些系统施工较精良的场合，不使沿程阻力的估算値与实际相比 偏大太多，造成不必要的设备选型的浪费。2.2气流模拟结果分析气流模拟计

14、算结果如表4与表5所示。将其与通风管道单位长度摩擦阻力线解图 对比后发现，当k= 15和k= 1.0时，四种情况下所求得的比摩阻与线解图提供 的数据基本相符。因此'可以采用表中数据对系统的沿程阻力进行估算。表4 k二1.5时烟道内的比摩阻适用层 数同吋使用 系数同时开 机数末段烟道内流 速(m/s)烟道尺寸 (mm)比摩阻 (pa/m)6100%66. 67250x3003. 0271870%135. 42400x5001.0619 3060%186500x5001. 1131 50306. 94600x6001. 14表5 k= 1.0时烟道内的比摩阻适用层 数同时使用 系数同时开

15、机数末段烟道内流 速(m/s)烟道尺寸 (mm)比摩阻(pa/m)6100%66. 67250x3002. 6871870%135. 42400x5000. 9019 3060%186500x5000. 9631 50306. 94600x6000. 982.3沿程阻力简化计算方法由于不同层数的用户向同一烟道内排烟，造成了等截面烟道内各段的流量、流速 都不同，沿程阻力计算虽然不难，但却非常烦琐。烟气在系统内流动所要克服的 沿程阻力可由公式(6)计算得到。(6)式中总沿程阻力 paai 第i段烟道的沿程阻力系数，对应于该层烟道内的流速vi1 每层烟道的长度，m ；n 总层数；n同时开机数。因沿程

16、阻力损失占i，所以，旦启j，其中，01=0 ；同时，由于入1的数量级与回的数量级相比非常小，因此(7)由式7可以推出烟气流经烟道系统产生的沿程损失的估算公式8 :式中，rmn对应于系统总风量的烟道比摩阻，pa/m。从上述公式推导可知，o表6列出了当k= 15，排风量为300m3/h、400m3/h、500m3/h及600m3/h时分别采用公式(6)、(8)计算6层、18层、30层 及50层烟道系统总沿程阻力的结果。表中数据显示，烟道系统总沿程阻力的理 论计算结果通常约为估算结果的80%90%，具体偏差値的大小和系统内的排风 量有关。因此，对式(8)乘以一个修正系数a，即可得到较为精确的系统总沿

17、程 阻力的估算値，不同排风量下a的取値如表7所示。(9)表6两种方法计算k二1.5时烟道系统总沿程阻力结果比较适 用时用 司使同i1j-排风量(m3/h)300400500600层数数 系开机数a亠凶(pa)a亠凶(pa)a亠凶(pa)>(pa)60% o11721a4005.48110%773.30%66071.503080 o 4 312&24表7不同排风量下修正系数a的取値排风量(m3/h)300400500600修正系数a0. 880. 880. 840. 82采用式9来计算烟道的总沿程阻力只要根据每户额定排烟量、层高、总层数、冋 时开机数等设计参数，再计算一个最大比摩阻即可，显然比逐层计算简单得多， 特别是在修改假设烟道截面积需要重新计算时，优越性更明显。需要说明的是， 笔者同时还在编制烟道设计的软件，将为设计提供更大的方便。3结论3.1设计高层住宅厨房集中排风系统时内壁绝对粗糙度k取为1.5mm是安全的， 若系统制作工艺精良，k可以取为1.0，甚至更小。3.2在已确定k値与