2. 中航商用航空发动机有限责任公司, 上海 200241
2. AVIC Commercial Aircraft Engine Co., Ltd., Shanghai 200241, China
焙烧竖炉是钒钛磁铁矿直接提钒工艺中的关键设备[1-2], 具有能耗低、炉内气氛可控等特点.焙烧过程在整个提钒工艺中至关重要, 其炉内的气固传热决定着钒钛磁铁矿焙烧的效果.
就床层本质而言, 钒钛磁铁矿焙烧竖炉属于大颗粒填充床的范畴.可借鉴大颗粒填充床内的气固流动与传热的相关研究, Al-Sumaily等[3]分别采用局部热力学平衡模型和非平衡模型分析了填充床内颗粒大小对流体流动和气固换热的影响.Pivem等[4]建立了逆流式移动床气固热力学非平衡传热模型, 讨论了基本物性参数、气固热当量等因素对传热过程的影响.Nijemeisland等[5]通过CFD模拟了圆柱形颗粒在竖管内的传热行为, 讨论了颗粒有无内热源对传热过程的影响.蒋鹭等[6]以欧拉多向流模型为基础, 建立了南钢球团竖炉的三维数理模型, 探讨了不同操作参数对竖炉炉温及焙烧产物Fe2O3产率的影响规律.Mei等[7]基于平均Re的N-S方程和κ-ε模型的双能量方程, 对球团竖炉整个均热段和冷却塔中共轭混合对流与三维湍流的热传递进行了数值模拟, 并研究了气体Re和Gr对炉内温度场分布的影响.
综上所述, 目前针对竖炉局部热力学非平衡的研究较少, 模拟过程中通常将气固温度看作一致; 而且气体动量方程阻力源项采用传统经验关联式, 模拟结果与实际情况偏差较大; 有关钒钛磁铁矿球团焙烧竖炉内气固传热过程及传热系数尚未进行针对性研究.针对以往研究的不足, 本文基于局部热力学非平衡模型, 通过软件二次开发功能UDF, 实现钒钛磁铁矿竖炉数学模型的2点改进:1)将竖炉内化学反应热以内热源的形式编译到能量方程中; 2)将在自制实验平台获得的料层压降公式定义到动量方程源项中.并在将球团矿下移速度编译到模型中的前提下, 建立了焙烧竖炉内气固逆流式气固传热过程三维稳态数值模型.基于此, 研究操作参数对传热过程的影响规律, 获得中试焙烧竖炉适宜的操作参数, 为后续的中试焙烧工艺的改进奠定基础.
1 模型的建立 1.1 物理模型本文以辽西某地区中试竖炉为参照对象, 建立竖炉的三维物理模型, 其结构如图 1所示.该竖炉为双炉膛对称结构, 因此探究单侧炉膛内的温度场和流场的分布就可以确定整座竖炉的热工行为.图 2为竖炉单侧的三维物理模型, 竖炉炉高6.2 m, 炉膛宽度1.1 m, 长4.2 m.
竖炉多进口、多出口的特性导致炉内气固流动及传热过程非常复杂, 在确保计算精度的要求下, 对焙烧竖炉的物理模型进行如下简化:
1) 钒钛磁铁矿焙烧竖炉运行工况稳定, 认为各参数在一定的范围内恒定;
2) 将焙烧竖炉内气体视为不可压缩流体, 忽略压力对其密度的影响;
3) 钒钛磁铁矿球团体积随温度变化很小, 忽略不计, 将整个床层看作各向同性的多孔介质.
1.2 数学模型本文所采用的气固传热与气体流动控制方程如下:
连续性方程为
(1) |
动量方程为
(2) |
式中:ρf为空气密度, kg · m-3; ui, uj分别为气体在i, j方向上的速度, m · s-1; Pij为表面的压力矢量, 包括静压力和气体黏性压力; gi为气体在i方向上的体积作用力, N · m-3; fi为作用在单位体积气体上的反方向阻力, N · m-3.
在动量守恒方程式(2)中增加动量损失源项Si, 用来描述气体流经竖炉内多孔介质时的动量输运过程:
(3) |
式中:1/α为黏性阻力系数; C2为惯性阻力系数; μ为气体动力黏度, Pa · s; u为气体流动速度, m · s-1.
通过实验获得描述焙烧竖炉内气流阻力特性的修正Ergun方程[8], 基于此得到多孔介质区域黏性阻力系数和惯性阻力系数.
(4) |
(5) |
(6) |
式中:ε为床层空隙率; dp为钒钛磁铁矿球团当量直径, m.
考虑到多孔介质计算域内, 流体和固体为不同的连续介质, 采用基于局部热力学非平衡状态的稳态双能量方程对其进行求解[9-10], 对竖炉内的气体和钒钛磁铁矿球团分别建立能量方程.
气相能量方程为
(7) |
固相能量方程为
(8) |
式中:ρs为钒钛磁铁矿球团的密度, kg · m-3; cs和cp分别为钒钛磁铁矿球团和空气的比热容, J · kg-1 · K-1; us和uf分别为颗粒下移速度和气体表观流速, m · s-1; Ts和Tf分别为钒钛磁铁矿球团和空气温度, K; λs和λf分别为钒钛磁铁矿球团和空气的导热系数, W · m-1 · K-1; hv为钒钛磁铁矿球团与气体间的对流体积换热系数[11], W · m-3 · K-1;Φ为钒钛磁铁矿球团焙烧化学反应内热源项, W · m-3.
1.3 边界条件竖炉下部冷却风入口采用速度入口边界条件, 将导风墙的下出口和竖炉上出口设为压力出口, 焙烧风入口采用速度入口边界条件.其余壁面设为绝热面, 不考虑炉壁的散热损失.
1.4 模型验证根据中试竖炉的产量计算得出, 单侧炉膛的球团下移速度为0.322 m/h, 此标准工况下, 焙烧风的流量和温度分别为3 456 m3/h和1 373 K, 冷却风流量和温度分别为3 024 m3/h和300 K.由于钒钛磁铁矿在竖炉内的各部分发生不同化学反应, 经计算得出该工况条件下炉内的化学反应热为665 W/m3[8].由于不同的焙烧温度所放出的反应热不同, 因此将化学反应热看作一个随高度变化的函数.
表 1为中试竖炉正常运行工况下不同高度位置空气温度测量值和数值模拟结果对比情况, 模拟结果的准确度是由实际测试结果与模拟结果的相对误差来衡量的.由表 1知, 测试值与模拟计算值存在一定的误差, 最大误差为7.91 %, 最小误差为3.38 %, 均在允许误差范围内.由此可知数值计算模型符合计算要求, 可用于探究炉内气固传热过程的数值研究.
利用已建立的模型, 分别探究影响竖炉内气固传热过程的3个主要参数:焙烧风流量(Qr)、冷却风流量(Qc)和球团矿下移速度(vs)对竖炉内气固传热过程的影响规律.
标准工况下竖炉内钒钛磁铁矿球团和气体温度分布如图 3所示.由图可知,球团矿在焙烧段的温度分布不均匀, 靠近焙烧风入口的位置, 球团的温度较高, 而远离焙烧风入口的位置, 球团的温度较低.
根据已有研究, 为了达到焙烧要求, 球团矿在竖炉内应达到的焙烧时间为4~6 h, 辽西新型钒钛磁铁矿的适宜焙烧温度为1 100~1 200 K[12].针对均热温度过低及焙烧时间过短的情况, 本文在此提出以焙烧时间为竖炉焙烧效果的评价指标, 通过研究竖炉操作参数变化对竖炉焙烧时间的影响, 确定适宜的操作参数.焙烧时间t的计算公式为
(9) |
式中:t为焙烧时间, h; H为达到焙烧温度料层高度, m; vs为钒钛磁铁矿球团下移速度,m/h.
2.1 焙烧风流量的影响当vs和Qc保持不变时, 焙烧风的通入量会影响t和温度.现有竖炉在实际生产中Qr为3 456 m3/h, 在Qc为1 512 m3/h的基础上逐渐改变Qr.考察Qr为2 523, 2 765, 2 995, 3 214, 3 456, 3 698, 3 905, 4 147, 4 389 m3/h时的温度分布, 并计算t, 如图 4所示.当Qr=2 523 m3/h时, t =3.65 h, 平均焙烧温度为1 010 K; 当Qr=2 995 m3/h时, t =4.6 h, 平均焙烧温度为1 087 K, 此时满足保温要求; 当Qr=4 389 m3/h时, 焙烧时间达到模拟最大值为6.7 h, 此时平均焙烧温度达到1 195 K.
随着Qr的增加, 焙烧时间和温度会随着增加, 后两者与前者呈非线性关系.Qr在2 523~3 214 m3/h时, t随着Qr的增加变化速率逐渐增大; 而当Qr在3 214~4 389 m3/h时, 此时t的变化速率逐渐减小, 原因在于当球团下移速度和冷却风流量保持不变时, 增大Qr会使钒钛磁铁矿球团温度升高, 竖炉沿高度方向达到保温温度要求的区间更大, t随之增长; 而随着Qr的增加, 球团温度升高, 焙烧风温度与球团温差减小, 换热速率降低, 球团焙烧温升速率降低.考虑到在达到保温要求的同时, 避免炉内温度过高引起的固结现象以及能源浪费, 应保持Qr在3 070~3 670 m3/h时比较适宜.
2.2 冷却风流量的影响在给定的初始冷却风流量条件下, 模拟得到的结果显示竖炉焙烧段温度偏低, t无法达到焙烧要求, 这是由于Qc过高造成的.针对此问题, 在实际生产Qc为3 024 m3/h的基础上, 逐渐减少Qc, 模拟Qc为3 024, 2 721.6, 2 419.2, 2 116.8, 1 814.4, 1 512, 1 210 m3/h的温度分布, 并计算t, 拟合曲线图如图 5所示.从图 5中可以看出,当Qc为1 210 m3/h时, t=5.3 h, 平均焙烧温度为1 150 K, 满足生产要求.图 5中, Qc与t和温度呈非线性关系, 因为随着Qc的增大, 流速逐渐增大, 换热速率加快, 单位时间内换热量增多, 钒钛磁铁矿球团沿高度方向的降温速率逐渐加快, t缩短且变化速率加快.而当温度降低时, 球团和空气比热减小, 温差减小导致换热速率降低, 焙烧温度降低速率逐渐减小.在满足保温要求的前提下, Qc应保持在1 210~1 430 m3/h比较合适.
为了使现有竖炉的传热效果得到最大化提升, 在确保钒钛磁铁矿球团在竖炉内的焙烧时间(4~6 h)充足和产量一定的前提下进行模拟, 模拟过程中, 通过改变竖炉的结构尺寸改变vs,竖炉原长度为4.2 m,表 2为改变竖炉长度计算得到的下移速度vs.
现有竖炉的vs为0.322 m · h-1, 逐渐改变vs, 模拟vs分别为0.193, 0.225, 0.258, 0.29, 0.354, 0.386, 0.419, 0.451 m · h-1时所得到的t.图 6为Qr=3 160 m3/h, Qc =1 468 m3/h时的t和焙烧温度随vs的变化关系曲线, 并与实际工况作对比分析, 从图中可以看出:随着vs的不断增大, t逐渐缩短.模拟结果显示t与vs呈非线性关系, 说明t不仅受到vs的影响, 还与球团自身氧化反应放出化学反应热量有关.Qr和Qc一定, vs越大, 单位质量的球团矿化学反应热总量越低, 达到焙烧要求的面积越小, 焙烧段的高度会相应减小.球团经预热到一定温度后才发生化学反应.随着球团下移速度的不断增加, 气固换热速率加快, 而焙烧风提供的热量一定.所以焙烧段球团温度的升高速率随vs的增加而降低, 进而使t随vs变化的变化量越来越小.
由图 6可知, 当vs=0.193 m/h时, 球团在炉内的t =8.02 h, 平均焙烧温度为1 170 K, 焙烧时间过长且vs过小, 不满足生产要求; 当vs=0.258 m/h时, t缩短至6 h, 平均焙烧温度降低至1 115 K, 此时已经满足保温要求, 说明t和vs都比较合适; 当vs=0.322 m/h时,t =4.9 h, 平均焙烧温度为1 078 K, 基本满足保温要求; vs逐渐增大至0.451 m/h时, t=3.43 h, 下移速度过大, t不满足生产要求.经计算, vs应控制在0.258~0.29 m/h较为适宜.
3 结论1) 以局部热力学非平衡模型为基础, 建立钒钛磁铁矿焙烧竖炉的三维稳态气固传热模型, 计算竖炉在实际工况条件下炉内各段气体温度值与实际测量值的最大误差小于7.91 %, 验证了数值模型的可靠性.
2) 随着焙烧风流量的增加, 钒钛磁铁矿的焙烧时间和温度逐渐增加, 且焙烧时间增长速率先增大后减小, 最后趋于定值; 随着冷却风流量的减少, 钒钛磁铁矿的焙烧时间和温度逐渐增加; 随着球团下移速度的增加, 竖炉内球团矿焙烧时间和温度逐渐减小, 且变化速率逐渐降低, 二者呈非线性关系.
3) 对于焙烧面积为6.8 m2的中试竖炉, 在要求焙烧时间为4~6 h的情况下, 其适宜的操作参数为:冷却风流量1 210~1 430 m3/h; 焙烧风流量3 070~3 670 m3/h; 竖炉球团下移速度0.258~0.29 m/h.
[1] |
张井凡, 董辉, 蔡九菊, 等.一种竖炉式氧化钠化焙烧方法及装置: 201110371297.9[P].2012-04-11.
( Zhang Jing-fan, Dong Hui, Cai Jiu-ju, et al.Method and device of oxidation-sodium roasting with shaft furnace type: 201110371297.9[P].2012-04-11. ) |
[2] |
董辉, 冯军胜, 李鹏, 等.
一种新型钒氧化钠化焙烧竖炉的研究[J]. 工业炉, 2012, 34(6): 5–8.
( Dong Hui, Feng Jun-sheng, Li Peng, et al. Research of a new-style roasting shaft furnace of vanadium oxide[J]. Industrial Furnace, 2012, 34(6): 5–8. DOI:10.3969/j.issn.1001-6988.2012.06.002 ) |
[3] |
Al-Sumaily G F, Nakayama A, Sheridan J, et al.
The effect of porous media particle size on forced convection from a circular cylinder without assuming local thermal equilibrium between phases[J]. International Journal of Heat and Mass Transfer, 2012, 55(13/14): 3366–3378.
|
[4] |
Pivem A C, Lemos M J S D.
Numerical simulation of a crossflow moving porous bed using a thermal non-equilibrium model[J]. International Journal of Heat and Mass Transfer, 2013, 67(4): 311–325.
|
[5] |
Nijemeisland M, Dixon A G, Stitt E H.
Catalyst design by CFD for heat transfer and reaction in steam reforming[J]. Chemical Engineering Science, 2004, 59(22): 5185–5191.
|
[6] |
蒋鹭, 黄山, 王天才, 等.
球团竖炉气固流动与焙烧过程耦合的三维数值模拟[J]. 东南大学学报(自然科学版), 2012, 42(2): 301–307.
( Jiang Lu, Huang Shan, Wang Tian-cai, et al. Three-dimensional numerical simulation of gas-solid flow and roasting process coupling in pelletizing shaft furnace[J]. Journal of Southeast University(Natural Science), 2012, 42(2): 301–307. ) |
[7] |
Mei D, Xing F, Wen M, et al.
Numerical simulation of mixed convection heat transfer of galvanized steel sheets in the vertical alloying furnace[J]. Applied Thermal Engineering, 2016, 93: 500–508.
DOI:10.1016/j.applthermaleng.2015.10.010 |
[8] |
孙用军.辽西钒钛磁铁矿直接提钒焙烧竖炉气固传热数值计算[D].沈阳: 东北大学, 2015.
( Sun Yong-jun.Numerical study of gas-solid heat transfer on roasting shaft furnace for vanadium titano-magnetite in western Liaoning province[D].Shenyang: Northeastern University, 2015. ) |
[9] |
Mahmoudi Y, Maerefat M.
Analytical investigation of heat transfer enhancement in a channel partially filled with a porous material under local thermal non-equilibrium condition[J]. International Journal of Thermal Sciences, 2011, 50(12): 2386–2401.
DOI:10.1016/j.ijthermalsci.2011.07.008 |
[10] |
Nakayama A, Kuwahara F, Sugiyama M, et al.
A two-energy equation model for conduction and convection in porous media[J]. International Journal of Heat and Mass Transfer, 2001, 44(22): 4375–4379.
DOI:10.1016/S0017-9310(01)00069-2 |
[11] |
王国胜, 董辉, 蔡九菊.
球团竖炉热工过程解析与模拟[J]. 钢铁, 2004, 39(4): 60–63.
( Wang Guo-sheng, Dong Hui, Cai Jiu-ju. Analysis and simulation of thermo process in pelletizing shaft furnace[J]. Iron and Steel, 2004, 39(4): 60–63. ) |
[12] |
邵颖聪.辽西钒钛磁铁矿直接提钒焙烧浸出工艺及竖炉热工过程研究[D].沈阳: 东北大学, 2015.
( Shao Ying-cong.Experimental study on roasting and leaching process of compact extracting vanadium process and analysis on thermo process in roasting shaft furnace[D].Shenyang: Northeastern University, 2015. ) |