高炉风口回旋区化学反应过程的数值模拟研究

《高炉风口回旋区化学反应过程的数值模拟研究》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高炉风口回旋区化学反应过程的数值模拟研究梁小平1，赵欣1，王雨1，八木顺一郎2（1.重庆大学材料与工程学院，重庆400044；2.东北大学多元物质科学研究所，日本仙台980-8577）摘要：通过建立高炉风口回旋区内气固两相的流动、传热、传质及化学反应过程的数学模型，对风口回旋区的化学反应过程进行了数值模拟研究，预测分析了鼓风条件对回旋区内气相成分分布的影响，以期为实际高炉鼓风操作的合理控制提供理论依据。关键词：高炉，风口，回旋区，化学反应过程，数学模型中国分类号：TF5361前言在炼铁高炉中，高温高速

2、的空气通过风口鼓入充满焦炭的固定床中，由于鼓风的作用，在风口前沿附近存在一个焦炭在其内作循环运动的区域即为风口回旋区。高炉回旋区内的化学反应情况,将直接影响着高炉下部煤气的分布、上部炉料的均衡下降、以及整个高炉内的传热传质过程[1-5]，因此，了解回旋区内的化学反应过程以及气体成分的变化规律对于炼铁生产实践具有指导意义。本文以攀钢2000m3高炉为研究对象，通过建立数学模型，对风口回旋区的化学反应过程进行数值模拟，预测回旋区气相成分的分布，研究鼓风条件的变化对气体成分分布的影响。2数学模型的建立如图1

3、所示的高炉下部区域即为本文数学模型的计算区域。本文采用Euler-Euler方法建立了风口回旋区的数学模型，模型中考虑了气固两相的流动、传热、传质及焦炭燃烧等化学反应对高炉下部区域气相成分的影响，模型的基本方程如下：(1)式中角标i分别代表气相g及固相s；为各相的体积分数，且；为待求变量，包括各相的速度、热焓（或温度）及各组分的质量分数（即浓度）；变量及为对应于待求变量的有效扩散系数及源项，其中源项包括各相之间的相互作用和化学反应等。将方程（1）展开，分别应用于气相及固相的流动、传热及传质过程即得到气

4、固两相的动量守恒方程、质量守恒方程及能量守恒方程；在一定的边界条件下，解出方程（1）中的待求变量即可获得风口回旋区的各相的速度、温度以及气相各组成的浓度分布，由此即可对回旋区气相成分分布进行预测计算。考虑到风口回旋区内所发生的化学反应，模型中气相的组成包括O2、CO、CO2、N2、H2O及H2这六种组分。图1风口回旋区数学模型的计算区域Fig.1Computationzonefornumericalmodelofraceway1.1化学反应模型回旋区内发生的化学反应主要包括鼓风中的氧气与焦炭发生的碳的

5、不完全燃烧反应、碳的完全燃烧反应反应、碳的熔损反应以及鼓风湿分的高温反应。在考虑化学反应影响的基本方程的源项时，需要根据反应动力学原理来确定化学反应速度，上述四个化学反应的数学模型方程如下。①碳的不完全燃烧反应C+O2＝CO（2）式（2）中，R1：化学反应①的反应速率，kg-C/s；：反应面积，m2；：O2的浓度；：化学反应的速率常数；：系数；kc1：扩散速率常数。②碳的完全燃烧反应C+O2＝CO2（3）式（3）中，R2：化学反应②的反应速率，kg-C/s；A：系数。③碳的熔损反应C+CO2＝2CO（

6、4）式（4）中，R3：化学反应③的反应速率速度，kg-C/s；：化学反应的速率常数；：CO2的浓度；：系数；kc3：扩散速率常数。④鼓风湿分的高温反应C+H2O＝CO+H2（5）式（5）中，R4：化学反应④的反应速率度，kg-C/s；：化学反应的速率常数；：H2O的浓度kg/mol；：系数；kc4：扩散速率常数。1.2异相间的热交换由于高炉下部的计算区域内同时有气固两相存在，各相之间有传热，在建立模型时必须计入这种相与相之间传热的影响，故在气固两相的传热模型方程的源项中涉及到各相之间的传热系数（g、s

7、代表气固两相）。气固间的传热系数采用经秋山修正的Ranz式计算[6]：（6）式（6）中，Re：雷诺准数；Pr：普朗特准数；：气体的导热系数；：焦炭颗粒的直径1.3异相间的动量交换气固两相的同时存在使得其运动亦相互影响，在各相的流动模型方程中这种影响效果已计入方程右端的源项中，表现为作用力（g、s分别代表气固两种相态的物质）。气固间的相互作用力采用经过修正后的Erugun方程[7]：（7）其中异相间的接触面积按下列公式计算：式（7）中，：粘度；：孔隙度；：速度；：形状因子；下标g、s、l：分别表示气、固

8、、液相。1.4边界条件上部入口处和风口入口处的参数作为已知值按攀钢实际生产条件确定，对称轴及下部出口处按处理，壁面处为无滑移条件，壁面传热按第一类边界条件处理。3计算及结果分析以SIMPLE算法编制数值求解计算程序，程序用FORTRAN语言编制。对于计算区域采用非均匀的正交交错网格进行网格划分，共划分为31×47＝1517个网格。由于风口直径远远小于炉体直径及炉高，为了提高计算精度，在网格划分时以风口直径为最小间距、按等比级数对其附近区域网格进行加密划分