金 焱，袁 辉，祝俊俊，王 炜，薛正良，程常桂
（武汉科技大学钢铁冶金及资源利用省部共建教育部重点实验室，湖北 武汉，４３００８１）
摘 要：基于Ｆｌｕｅｎｔ软件，采用预混燃烧模型对高炉风口回旋区内温度场进行模拟，分析风量、风压、喷煤量等参数对高炉风口回旋区内温度分布的影响。结果表明，随着风量、风压、喷煤量的增加，风口回旋区内温度最高处离风口端部的距离逐渐增大，风口回旋区内最高温度逐渐降低；风口的堵塞会使风口回旋区内温度最高处与风口端部的距离缩短，使风口回旋区内最高温度升高。
关 键 词：高炉；风口回旋区；温度分布；风量；风压；喷煤量；温度模拟
高炉风口是高炉炼铁送风的重要部位，承受炉内高速煤粉的磨蚀、高温炉气的冲刷以及炉料的撞击，极易损坏［１－４］。风口前端部的熔蚀会导致风口失效［５］，其寿命的长短直接影响到高炉能否顺利运行、获得高产及降低炼铁成本，因此风口回旋区的温度分布对风口的受损状况有很重要影响。有研究人员对高炉风口理论燃烧温度进行分析，得到计算风口理论燃烧温度的经验公式，以及煤质、富氧量、喷煤量、鼓风湿度、热风温度等因素对风口理论燃烧温度的影响［６－７］。张丽丽等［８］应用气固两相流动理论和煤粉燃烧数学模型对煤粉燃烧过程进行了数值模拟，能准确预测风口回旋区的燃烧情况。本研究基于Ｆｌｕｅｎｔ软件，采用预混燃烧模型对高炉风口回旋区内温度场进行模拟，分析风量、风压、喷煤量等参数对风口回旋区内温度分布的影响，以期为优化高炉风口工艺参数提供参考。
１ 模型的建立
１．１ 控制方程
本文选用Ｆｌｕｅｎｔ软件中的预混燃烧模型对风口回旋区内煤粉燃烧和煤气温度进行模拟。风口回旋区内温度场模拟的基本方程为：
能量守恒方程
式中：ρ为流体密度，ｋｇ／ｍ３；ｕｉ（ｉ＝１，２，３）为沿ｉ方向的速度分量，ｍ／ｓ；ｘｉ（ｉ＝１，２，３）为坐标分量；Ｔ 为燃烧初始温度，Ｋ；Ｋｅｆｆ为有效传热系数，Ｗ／（ｍ２·Ｋ）。
质量守恒方程
式中：ｃ为无因次反应进程变量；Ｓｃｔ为湍流施密特准数；
为平均反应速率；μｔ为湍流黏度，Ｐａ·ｓ；Ｓｃ为反应进程源项。
１．２ 三维模型
运用Ｆｌｕｅｎｔ软件对高炉送风系统进行建模，选定建模尺寸（单位：ｍ），设定建模比例为１∶１，建立的高炉送风系统三维模型如图１所示。采用自由划分的方式对模型进行网格划分，单元尺寸设定为０．３，共得到２０７９５０６个网格和３５９６７３个网格节点。
１．３ 假设条件及边界条件
假设风口回旋区内煤气为不可压缩流体；忽略高炉炉缸耐火材料与煤气、焦炭的物理化学反应及高温等条件产生的副反应；入流热风垂直于风口方向；忽略壁面效应对流动区域的影响。
边界条件为：
（１）入口边界条件：以热风流量作为入口条件。
（２）炉缸壁面：流体采用壁面函数法计算，高炉壁面采用无滑移边界条件，高炉内壁面附近黏性底层中的流体采用标准壁面函数处理，默认为耦合条件。对壁面附近的边界层进行网格划分和局部加密。
（３）对称面：为简化计算，两风口间的对称面设为对称边界，即垂直于对称面的速度分量及所有其它变量沿对称面法线方向的梯度为零。
１．４ 计算方法
通过ＩＣＥＭ 软件划分模型的离散网格，并将网格导入Ｆｌｕｅｎｔ软件中，采用３Ｄ 隐式差分分离求解的方式进行计算。
２ 结果与分析
为了分析风口参数对风口回旋区内温度分布的影响，模型计算时，利用某钢铁厂３２００ｍ３ 高炉的实际生产数据，选取了５月中某３个工作日的风口数据，这３个工作日中没有风口堵塞；另又选取了１月中某３个工作日的风口数据，其中２个工作日有风口堵塞。
２．１ 风量对风口回旋区内温度分布的影响
风量对风口回旋区内温度分布的影响如图２所示。从图２（ａ）中可以看出，随着风量的增加，风口回旋区内温度最高处离风口端部的距离逐渐增大，其线性相关系数为０．８６，表明风量对回旋区内温度最高处离风口端部的距离影响较大，这是因为，增加风量可以增加高炉的鼓风动能，而鼓风动能越大，会使风口回旋区和燃烧带扩大，使回旋区的深度增加，因此风口回旋区内温度最高处离风口端部的距离也就越大。从图２（ｂ）中可以看出，随着风量的增加，风口回旋区内最高温度呈下降趋势，这是由于风量增加加快了炉缸内流体的对流换热，有利于炉缸内温度趋于均匀，炉缸内煤气分布也更加均匀，导致风口回旋区内温差下降。
２．２ 喷煤量对风口回旋区内温度分布的影响
喷煤量对风口回旋区内温度分布的影响如图３所示。从图３（ａ）中可以看出，随着喷煤量的增加，风口回旋区内温度最高处离风口端部的距离逐渐增大，其线性相关系数为０．７７，这是因为，增加喷煤量可以增加高炉的鼓风动能，鼓风动能越大，使回旋区和燃烧带扩大，使回旋区的深度增加，因此风口回旋区内温度最高处离风口端部的距离也就越大。从图３（ｂ）中可以看出，风口回旋区内最高温度随喷煤量的增加而呈下降趋势，这是因为，喷煤时煤粉中挥发分的分解需要吸收较多的热量，喷煤量越大，吸收的热量越多，使得风口温度降低。
２．３ 风压对风口回旋区内温度分布的影响
风压对风口回旋区内温度分布的影响如图４所示。从图４（ａ）中可以看出，随着风压的增加，风口回旋区内温度最高处离风口端部的距离逐渐增大，其线性相关系数为０．８４，表明风压对回旋区内最高温度处离风口端部的距离影响较大，这是因为，增加风压会使风口回旋区的深度增加，从而导致回旋区内温度最高处离风口端部的距离越远。从图４（ｂ）中可以看出，风口回旋区内最高温度随风压的增加呈下降趋势，且线性相关系数达到０．８１，表明风压对风口回旋区内最高温度影响较大，这是因为，风压增加使高炉内热风体积减小，风口回旋区内流体对流换热加强，使其温度分布更加均匀，导致风口回旋区内最高温度降低。
２．４ 风口堵塞状况对风口回旋区内温度分布的影响
风口堵塞状况对风口回旋区内温度分布的影响如图５所示。从图５（ａ）中可以看出，当风口堵塞５个时，其回旋区内最高温度最高，风口堵塞１个时，其回旋区最高温度和无堵塞时的最高温度差别不大，这是因为，有一定数量的风口被堵塞时，会使其它未堵塞的风口流体流速变大，使反应物混合更加充分，燃烧反应加强，热量来源充足，使风口回旋区最高温度升高。从图５（ｂ）中可以看出，堵塞５个风口时，回旋区内温度最高处距离风口端部最近，风口堵塞１个时回旋区内最高温度的位置和无堵塞时差别不大，这是因为，风口堵塞较多时，会使部分风口回旋区的深度减小，从而导致风口回旋区的温度最高点与风口端部的距离减小。由此可见，风口的堵塞会使风口回旋区内温度最高处与风口端部的距离缩短，使风口回旋区内最高温度升高。
３ 结论
（１）风口回旋区的温度分布与风口工艺参数的变化相关，随着风量、风压、喷煤量的增加，风口回旋区内温度最高处离风口端部的距离逐渐增大，风口回旋区内最高温度逐渐降低。
（２）风口的堵塞会使风口回旋区内温度最高处与风口端部的距离缩短，使风口回旋区内最高温度升高。
参 考 文 献
［１］Ｊｏｈｎ Ｇ Ｍａｔｈｉｅｓｏｎ，Ｊｏｈｎ Ｓ Ｔｒｕｅｌｏｖｅ，Ｈａｒｏｌｄ Ｒｏｇｅｒｓ．Ｔｏｗａｒｄ ａｎ ｕｎｄｅｒｓｔａｎｄｉｎｇ ｏｆ ｃｏａｌ ｃｏｍｂｕｓｔｉｏｎ ｉｎｂｌａｓｔ ｆｕｒｎａｃｅ ｔｕｙｅｒｅ ｉｎｊｅｃｔｉｏｎ［Ｊ］．Ｆｕｅｌ，２００５，８４：１２２９－１２３７．
［２］樊勇保，李晓桥，李玲，等．基于Ｆｌｕｅｎｔ的高炉风口流场和温度场的模拟［Ｊ］．特种铸造及有色合金，２００９，２９（４）：３２４－３２６．
［３］Ｃｈｅｎ Ｃｈｉｎｇｗｅｎ．Ｎｕｍｅｒｉｃａｌ ａｎａｌｙｓｉｓ ｆｏｒ ｔｈｅ ｍｕｌｔｉｐｈａｓｅ ｆｌｏｗ ｏｆ ｐｕｌｖｅｒｉｓｅｄ ｃｏａｌ ｉｎｊｅｃｔｉｏｎ ｉｎｓｉｄｅ ｂｌａｓｔｆｕｒｎａｃｅ ｔｕｙｅｒｅ［Ｊ］．Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃａｌ Ｍｏｄｅｌｉｎｇ，２００５，２９：８７１－８８４．
［４］李洋龙，程树森，陈川．高炉风口风量分配数学模［Ｊ］．东北大学学报：自然科学版，２０１６，３７（３）：３５７－３６２．
［５］李炎．基于ＦＬＵＥＮＴ的高炉风口温度场和流场模拟及结构优化［Ｄ］．武汉：武汉科技大学，２０１５：１－５．
［６］李肇毅．高炉风口理论燃烧温度分析［Ｊ］．宝钢技术，２０１１（４）：５－７．
［７］王颖生，竺维春．关于首钢高炉风口理论燃烧温度计算经验公式的探讨［Ｊ］．首钢科技，２００２（５）：２５－２７．
［８］张丽丽，史岩彬，陈举华．高炉风口回旋区煤粉燃烧过程三维数值模拟［Ｊ］．青岛大学学报：工程技术版，２００５，２０（４）：４４－４８．