余热锅炉进口烟道流场模拟报告(杭锅？).docVIP

余热锅炉进口烟道流场模拟报告 杭锅？ 1 前言 为了分析余热锅炉前部烟道（从炉前水平烟道进口到余热锅炉内高压过热器管束出口的这一段烟气流通区域）内的烟气流动和换热情况，在大型流体力学计算软件FLUENT-CFD开发平台上进行了数值模拟工作，得到了这段区域内的烟气速度和温度分布。 2 模拟对象的几何建模和网格划分 如前言中所述，数值模拟的对象为余热锅炉的进口烟道，包括燃机出口，锅炉进口烟道和高压过热器管束换热区这三个部分，其整体几何模型如图1所示。其中图中深绿色区域为高压过热器管束。 图1 模拟对象的整体几何模型 在建模过程中，对高压过热器换热模块进行了合理的简化，将原来的管束结构简化成相同外形尺寸的长方体，这主要是为了配合使用porous media model（多孔介质模型）。在完成几何模型的建立后，使用gambit软件对其进行网格划分，整个模型共划分网格数约为60万，网格划分情况如图2所示。 图2 模型网格划分情况 3 数值模拟方法 3.1 控制方程及求解方法 将烟气的流动视为三维稳态定常流动，采用标准k-ε双方程模型，该模型的控制方程组如下： 连续方程： 动量方程： 湍动能?k方程： 湍动能耗散率方程： 其中： ：烟气密度；,,分别为x，y，z坐标方向的速度分量； ：压力； ：湍流动能；：湍流动能耗散率；：有效粘性系数；：湍流动能剪力产生力， ：动力粘性系数, ； 方程中的常数=1.44，=1.92，=0.09 由于问题中还涉及烟气与过热器管束的换热，所以还要引入能量方程。 在数值计算过程中，上述方程的离散全部采用一阶上风格式，并选用SIMPLE算法来处理压力场和速度场的耦合。 3.2 选用的特殊模型 为了模拟高压过热器换热管束，引入porous media model，即多孔介质模型。 如前所述，在建模过程中将高压过热器换热管束简化成了长方体区域，现将该区域定义为多孔介质区，并设定内部阻力系数和单位体积能量源项。该模型不仅可以模拟烟气流经过热器管束时的压力损失和换热情况，还能通过设定多孔性（定义为多孔介质区流体的体积分数）参数来反映管束区流通截面缩小引起的烟气流速的增加。其中多孔区域的设置如表1所示。 表1 多孔区域参数表 模型进口烟温（k） 模型出口烟温(k) 烟气流量(kg/s) 高压过热器内换热量 （MW） 流经高压过热器后烟气压降 （Pa） 图5 Y＝0截面温度场 K 图6 Y＝0截面速度场 m/s 图7 Y＝0截面速度矢量 图8 Z＝0截面温度场 K 图9 Z＝0截面速度场 m/s 图10 Z＝0截面速度矢量 从图中可以看出，由于余热锅炉进口烟道中无吸热部件，从锅炉进口到进入高压过热器之前的这段区域内，烟气的温度没有发生变化。在进入换热区之后，由于烟气放热，其温度有了明显的下降。换热区内的温度分布层次鲜明，反映出了温度由高到低渐进变化的过程。另外，整个换热区内的温度分布情况都趋于一致，可见烟气与高压过热器管束之间的换热是比较均匀的。在出了换热区域后，烟气的温度基本不再变化。 在速度方面，由于燃机出口，既锅炉进口的烟气速度不均匀，使进口烟道中存在明显的速度梯度。在烟道高度方向，由于流通截面积迅速扩大，高速的烟气扩散能力差，所有存在上下两个回流区。在烟道水平方向，由于流通截面积未扩散，没有回流区。烟道进口中心的回流区是燃机公司提供的回流速度造成的。当烟气到达高压过热器（既多孔区域）时，由于烟气速度越高，所受阻力越大，反之阻力越小。根据力学平衡原理，高速烟气必将进行分流，使烟气在高压过热器中所受阻力趋向相同，图中模拟结果与此理论分析相吻合。 图11 X＝10截面温度场 K 图12 X＝10截面速度场 m/s 图13 X＝10截面速度矢量图 图11～13是高压过热器进口截面模拟结果。在高压过热器进口截面将有一部分高速烟气流向低速区域，所以此截面上的速度分布比较复杂，如图13所示。由于惯性作用，截面各处烟气速度数值仍存在差距，如图12所示。烟气速度分布的不均匀造成吸热不均匀，吸热不均匀又造成烟气温度的不均匀，如图11所示。 图14 出口截面温度场 K 图15 出口截面速度场 m/s 图16 出口截面速度矢量图 图14～16是模型出口截面模拟结果。烟气经过高压过热器时，由于前面所分析的力学平衡，速度不均匀性将受到改善。从模拟结果上看，当到达模型出口截面时，速度方向趋向一致，如图16所示；速度大小趋向相同，如图15所示。但由于在高压过热器中烟气速度一直存在不均匀，出口截面上的烟气速度也不均匀，如图14所示，但数字相差不大，不影响锅炉性能。 5 小结 通过对余热锅炉进口烟道的数值模拟，结果显示，余热锅炉螺旋鳍片管的均流效果好。本项目燃机出口烟气速度极不均