1、高温蓄热式加热炉中的然气流流动分析李健 1，赵博宁 2摘要：本文运用 Fluent 软件对天然气在高温蓄热式加热炉中的燃烧过程进行了模拟研究。比较分析了影响炉膛内气流温度、流场以及烟气 浓度分布特征的主要因素。研究结果表明：当天然气和空气的预热温度都不变时，改变天然气和空气射流间的夹角，同 样会改变炉膛内的温度分布情况，炉内的最高温度几乎呈线性降低，炉内平均温度呈二次曲线变 化先增大后减小。当 烧嘴夹角为 35时，炉内和钢坯表面温度比较均匀，并且 NO 生成量少， 对钢坯的加热质量和降低环境污染非常有利。关键词：天然气；流场；温度场；浓度场引言天然气是一种优质的气体燃料。因其热值高、质量稳定、
2、设备简易、便于控制以及污染小等优点，在世界能耗结构中所占的比例越来越大。最近 10 年，很多企业已经对工业炉进行了燃天然气改造 12。本文以天然气在高温蓄热式锻造加热炉上的应用为背景，对其燃烧参数和炉膛结构参数进行了热态模拟。在本研究中，列微分方程组是研究的一个难点。再者，对微分方程组离散化求解也相当的困难 34。而这些数学上的困难可以直接借助于 Fluent 软件模型来求解。Fluent 求解的思路和步骤是首先根据实际情况抽象和建立物理模型，然后选择求解模型，设置边界条件和初始条件，设置迭代和控制参数并进行计算，最后进行后处理和分析得出结论，或者根据实验结果重新调整参数进行在计算直至达到可接
3、受的结果 56。1 物理模型和计算工况1.1 燃烧过程数值模拟的物理模型流体的流动，一般分为层流和湍流。由于炉膛的尺寸较大，燃料入口处流速较高，流动往往是湍流状态，所有的物理量都是空间和时间的随机变量，但是湍流流动仍遵循连续介质一般运动，并具有一定规律的统计学特征，其瞬时流动仍满足粘性流体流动方程，描述燃烧流场各瞬时的微分方程如下 8：连续性方程: （1）0iiutx动量方程： （2） ()jiiij efj ijjiupt xx组分方程： （3）sjsssj jjmumRtx能量方程： （4）()j jjP PjPhhCTCTTSt x k 方程: （5） 3/2()jjj jijiieft
4、 pkkCktxx方程: 21()jj jij jjjii efCut kkx （6）式中， 为流体密度,kg/m 3； 为时间，s；P 为压力,Pa; 代表 1，2，3， 为直角坐标系 j 方t ,ju向速度矢量，m/s； 为直角坐标系 j 方向的坐标； 为组分的质量分数，%； 是扩散项或源项的jxsmsR生成率，s -1； 为热焓源项，W/h 3; 分子粘性系数， ,； 为湍流动力粘性系数， ,hS0aPs1aPs其中 为系数， =0.09； 为湍流有效粘性系数，2/tuckucuef， ；k 为湍流动能， 为湍流耗散率， =1.44， =1.92，aP01efC21.2 计算工况本文共设
5、计以下几种气流夹角，来考察气流夹角对炉内温度场，流场分布及燃烧产物的影响。由于本文采用的是三喷口式烧嘴，中间喷口喷天然气，两边喷口喷空气，所以计算时天然气射流采取直喷方式，只改变两侧空气射流喷射方向，就可以实现改变两者射流夹角的目的。计算工况：烧嘴安装高度为650mm，烧嘴喷口间距为 450mm，空气预热温度为 800mm，燃气预热温度为 800mm，空气喷射方向与烧嘴轴线的夹角为 0 o 50 o之间分布。2 实体模型的建立以及简化2.1 实体模型的建立以 某 锻 造 厂 采 用 高 温 空 气 燃 烧 技 术 改 造 后 的 蓄 热 式 加 热 炉 为 研 究 对 象 ， 燃 烧 空 间
6、的 有 效 长 度 为4200mm， 宽 为 4200mm， 高 度 为 3600mm.， 炉 墙 两 侧 分 别 布 置 两 组 烧 嘴 ， 两 侧 烧 嘴 交 替 的 进 行 喷 气 和 排 烟 ，烧 嘴 所 用 的 天 然 气 总 量 为 759.02m3/h， 理 论 所 需 的 助 燃 空 气 量 为 7796.91m3/h。 烧 嘴 有 效 流 通 面 积 为 圆 截 面 ，燃 料 采 用 工 业 天 然 气 ， 主 要 成 分 为 CH4， 其 平 均 热 值 为 35000J/m3。 由 于 天 然 气 理 论 空 气 耗 量 为 9.4， 为 了 避免 空 气 喷 口 过 大
7、 ， 故 设 计 烧 嘴 有 三 个 喷 口 ， 中 间 喷 口 喷 天 然 气 ， 直 径 为 240mm,两 侧 喷 口 喷 空 气 ， 直 径 为360mm,空 气 喷 口 与 天 然 气 喷 口 可 成 一 定 夹 角 ， 相 邻 喷 口 组 成 一 个 燃 烧 单 元 ， 依 靠 横 向 动 量 形 成 湍 流 ， 以 便组 织 炉 膛 内 燃 料 燃 烧 。 同 时 ， 空 、 燃 气 烧 嘴 喷 口 水 平 安 装 ， 有 利 于 空 、 燃 气 的 充 分 混 合 ， 避 免 烧 嘴 火 焰 直烧 加 热 钢 坯 。 锻 造 加 热 炉 运 行 参 数 见 表 1。表 1 锻造
8、加热炉主要参数Table 1 The major of parameters of forge heating furnace名称 数据最高炉温 Tmax/ 1300最大生产能力/Kg/h 300燃料种类 天然气燃料低发热量/KJ/m 3 35炉底热强度/GJ/(m 2h) 1.506最大燃料消耗量/ Km 3/h 760空气过剩系数 1.022.2 建立实体模型和划分网格计算区域为三维区域，炉膛的结构如图 1 所示，结构不算复杂，但燃气入口相对炉子尺寸小，网格过大对入口空间内部速度和温度等计算不利，网格过小造成网格数目巨大，计算时间延长。本文采用固定正交的非均匀网格，在保证喷口出网格的墙体下
9、，网格的大小又靠近炉墙和喷口出的 30mm 渐渐把大为100mm 至 300mm 不等，共生成网格大约为 117891 个。网格生成后如图 2 所示，输出 3D 格式的 mesh 文件。2.34 边界条件的设定设定边界条件的类型，指定流动变量和热量在边界条件处的值，是 Fluent 使用过程中很关键的一部分，设定边界条件必须小心谨慎燃料和助燃空气入口定义设定为速度入口，各入口处射流的速度和方向依具体工况决定，表 2 给定射流入口的直径和射流组分的质量百分数。烟气出口定义为压力出口，所以烟气出口方向应与入口方向一致。表 2 射流入口边界条件（物质组分为质量百分比）Table 2 The boun
10、dary condition of the inlets入口 直径（ mm） 物质组分（ %）CH4 O2 N2gasinlet_1 240 1 0 0gasinlet_2 240 1 0 0airinlet-1 360 0 0.2201 0.7299airinlet-2 360 0 0.2201 0.7299airinlet-3 360 0 0.2201 0.7299airinlet-4 360 0 0.2201 0.72993炉内温度分布的分析图 1 至图 3 为炉内的平均温度，最高温度和平均温度的差，以及出口处平均温度随空气射流与天然气射流之间的夹角的变化。由图看到，当夹角由 5增加到
11、50时，炉内的最高温度几乎呈线性降低，炉内平均温度近似成二次曲线变化先增大后减小。当夹角处于 2035o时对应的炉内的平均温度较高；处于2540时最高温度与平均温度差较小。这主要是因为天然气和空气是由不同的喷口进入炉膛，在进入炉膛前没有混合，在进入炉膛后边扩散混合边燃烧。如果喷射角度过大，天然气和空气混合较早，燃烧太快，燃烧区靠近喷口附近，火焰高温不能很快传到炉膛中部或后区造成炉膛温度分布不均匀。如果夹角过小，天然气、空气在还未进行充分混合燃烧时亦被吸入对侧的蓄热室通道，所以夹角过小或过大都会造成炉内温度分布不均匀（注： Tmax代表最高温度； Tave代表平均温度） 。05105205305
12、4050245024750250最高温度/K空 气 喷 射 方 向 与 烧 嘴 轴 线 的 夹 角 /图1 炉内最高温度随的变化Fig.1 the Tmax of the furnace in every0510520530540502172025250平均温度/K空 气 喷 射 方 向 与 烧 嘴 轴 线 的 夹 角 /图2 炉内平均温度随的变化Fig.2 the Tave of the furnace in every051052053054050225027530325350最高温度与平均温度差/K空 气 喷 射 方 向 与 烧 嘴 轴 线 的 夹 角 /图 3 炉内最高温度和平均温度的
13、差随 的变化Figure 3 the Tmax - Tave of the furnace in every图 4 为各个工况下，沿天然气喷口中心轴线方向上的温度分布图。由图可看到，各工况均在距烧嘴喷口 1m 左右达到最高温度，随后温度变化不大，比较平缓。当采用蓄热式燃烧技术时，在前后两个换向周期，天然气、空气从相反方向的喷口喷入炉膛内混合燃烧，若控制好的话，后半个换向周期炉膛内的温度分布应与前半个换向周期的分布相反。即在后半个换向周期燃烧时，同样在距喷口 1m 左右，炉内的温度达到最高，随后变化不大，变得比较平稳。所以若给钢坯加热时，钢坯最好放置在距前后喷口 1mm 左右的位置。4 炉内烟气
14、的浓度分布特征天然气燃烧过程中生成的 NOx 中，NO 占 90%左右，其余为 NO27，因此本文中烟气的含量中只考虑了NO，而没有考虑 NO2。 图 6 至图 11 是烟气的浓度（燃烧产物的质纵量浓度）随空气射流与天然气射流夹角的变化图。从图看出，显然，在加热炉的实际生产调试中，到属于调试敏感区域。也就是说，控制如此小范围的夹角，对于伴随着环境的振动以及毛坯材料的碰撞的大尺寸加热炉，将会使设备性能的稳定性非常差。理论上，研究到这样质量浓度变化巨烈的区域是最有吸引力的，但是研究适于实际应用的加热炉空气射流与天然气射流的夹角对烟气的浓度的影响是本此工作的重点。本实验选择 从5到 50。每 10为
15、一个间隔。从这几个图上可以看出，有喷射角度和没有喷射角度时，除 NO 外，其-2.015.-.051.2.0546180246温度/K喷 口 几 何 轴 线 /m case1(0o)253 4(o)case6 7(5o)89 (o)图 4 天然气喷口中心轴线上得的温度分布图Figure 4 The temperature distribution along the fuel burners Geometry axis余烟气组份的浓度相差很大。当喷射角度为 5，CH 4的浓度和 O2的浓度是没有喷射角度的 1/12， CO2、H 2O 和 N2的浓度分别为没有喷射角度 1.56、1.51、1.
16、1 倍。这说明有喷射角度比没喷射角度时，天然气和空气混合充分，燃烧比较完全。并且有喷射角度时，随着喷射角度的增加，CH 4的浓度随之增加，O2和 N2的浓度反而减少，CO 2和 H2O 的浓度呈先增大，后减小的变化趋势，而 NO 的浓度出现先减小后增大的趋势。这说明喷射角度并非越大，天然气和空气混合就愈充分，对燃烧越有利。图5 出口处CH 4的质量浓度 图6 出口处O 2的质量浓度Fig.5 The average mass fraction Fig.6 The average mass fraction of CH4 in the burners outlet of O2 in the bu
17、rners outlet空气喷射方向与烧嘴轴线的夹角/ o空气喷射方向与烧嘴轴线的夹角/ oCH4的质量浓度/% 051052053054050.10.23.405.670.89 051052053054050.10.23.405.670.89O2的质量浓度/%空气喷射方向与烧嘴轴线的夹角/ oCO2的质量浓度/% 051052305450.9.10120.3.140.5NO2的质量浓度/%空气喷射方向与烧嘴轴线的夹角/ o051052305450.70.8.90.1.0.12图7 出口处CO 2的质量浓度 图8 出口处H 2O的质量浓度Fig.7 The average mass fract
18、ion Fig.8 The average mass fraction of CO2 in the burners outlet of H2O in the burners outlet图9 出口处N2的质量浓度 图10 出口处NO的质量浓度Fig.9 The average mass fraction Fig.10 The average mass fraction of N2 in the burners outlet of NO in the burners outlet051052053054050.6.70.68.90.7.10.72.30.740510520530540500.17
19、.20.5.20.75.30.25.30.75.4空气喷射方向与烧嘴轴线的夹角/ o 空气喷射方向与烧嘴轴线的夹角/ oNO的质量浓度/%N2的质量浓度/%观 察 图 10， NO 的 浓 度 随 喷 射 角 度 变 化 ， 当 喷 射 角 度 增 加 到 35之 前 ， NO 的 浓 度 在 逐 渐 减 少 。 当 由 35增 加 到 50时 ， 浓 度 又 开 始 增 加 。 NO 的 生 成 量 和 温 度 有 很 大 关 系 5， 当 炉 内 温 度 在 1500 附 近 变 化 时 ， 温 度每 增 加 100 ， NO 的 反 应 速 度 将 增 大 6 7 倍 。 当 温 度 为
20、 2000 时 ， NO 的 生 成 量 就 极 为 迅 速 。 结 合 炉 内 最 高温 度 随 喷 射 角 度 的 变 化 ， 看 到 炉 内 最 高 温 度 在 角 度 增 大 到 35之 前 ， 一 直 在 降 低 ， 但 是 过 了 35后 ， 开 始 缓慢 增 加 ， 所 以 NO 的 浓 度 也 开 始 增 加 。4 炉内流速的分布图 11-14 为喷口横截面的烟气速度分布图，对比这几种工况下的流速图，发现，当夹角为 2035时，在炉中心形成很大的漩涡，流速分布比较均匀的，因此温度比较均匀。而角度大于35湍流太大，形成了漩涡太大，不能使高温燃气与高温空气达到良好细化混合，结果致使
21、燃料燃烧图 12 20时，喷口横截面速度向量图 Figure12 The velocity distribution of flow in the burners sectional surface when=20o 图 11 5时，喷口横截面速度向量图Figure 11 The velocity distribution of flowin the burners sectional surface when=5o图 13 35时，喷口横截面速度向量图 Figure 13 The velocity distribution of flowin the burners sectional su
22、rface when=35o图 14 43时，喷口横截面速度向量图Figure 14 The velocity distribution of flowin the burners sectional surface when =43o不 充 分 ， 所 以 当 角 度 增 大 时 ， 出 口 处 的 天 然 气 的 浓 度 也 随 之 增 加 ， 而 出 口 温 度 随 之 减 小 。 角 度 小 于 35时 ，主 要 以 小 漩 涡 为 主 ， 虽 然 有 利 于 燃 料 和 空 气 的 混 合 ， 但 是 小 漩 涡 又 是 主 要 的 能 量 耗 散 区 ， 所 以 炉 内 最 高 温
23、 度比 较 高 ， 平 均 温 度 比 较 低 。5 结论本文运用 Fluent 对天然气在高温蓄热式加热炉中的燃烧过程进行研究。考察了天然气和空气的流速不变时，考察了改变射流间的夹角，炉膛内温度、流速以及烟气的浓度分布情况。本节得出以下结论： （1）空气射流喷射方向与天然气射流方向成一定的角度喷射时，不仅可以使天然气和空气射流剧烈混合，而且可以避免高温气流冲涮对面的炉墙，提高加热炉的寿命。（2）当夹角有 5增加到 50时，炉内的最高温度和出口处的平均温度几乎呈线性降低，炉内平均温度成二次曲线变化先增大后减小，并且炉膛内的高温区越靠近喷口处。（3）各工况均在距烧嘴喷口 1m 左右达到最高温度，
24、随后温度变化不大，比较平缓。所以若给钢坯加热时，钢坯最好放置在距前后喷口 1mm 左右的位置。（4）综合各方面的因素，当烧嘴夹角为 35时，炉内和钢坯表面温度比较均匀，并且 NO 生成量少，对钢坯的加热质量和降低环境污染非常有利。参考文献：1唐炼.世界能源供需现状与发展趋势.国际石油经济J.2005,1:3033.2Nabil Rafidi,Wlod zimierz Blasiak. Heat transfer characteristics of HTAC heating furnace using regenerative burners. Applied Thermal Engineer
25、ingJ.2006, 26:20272034.3IshiT, Zang C, Sugiyama S. The numerical and experimental study of non-premixed combustion flames in regenerative furnaces. Journal of Energy Resource TechnologyJ.2000,122:287293.4Gupta A K, Bolz S. Effect of air pre-heat temperature and oxygen concentration on the flame structure and emission. Journal of Energy Resources TechnologyJ.1999,121:209216.5 郝玉玲，邹宗树，蔡九菊等.燃气喷射方向对高温空气燃烧室内流动影响的数值模拟J.东北大学学报,2003,3: 276279.6曹诗倬.锻造加热设备M. 北京:机械工业出版社,1989:525,5159.7苏正川.高温空气燃烧的实验和数值模拟研究. 同济大学硕士学位 论文D,2003: 714.8任雁秋.加热炉热过程数学模型J. 包头钢铁学院学报。1999，18：424424（原刊于北京工业大学学报2009 年 8 月）
