2019, Vol. 40
铁水预处理脱硫是降低钢中有害元素硫含量的主要途径, 目前最常用的脱硫方法为KR法和喷吹法[1-2].镁蒸气铁水脱硫是气液反应, 目前应用镁脱硫方法的脱硫效率都较低, 很重要的原因是铁水中产生的镁蒸气气泡较大[3]; 因此, 提高脱硫效率的关键是将产生的镁蒸气大气泡破碎为小气泡, 使其均匀分散在铁水中, 通过增加气液接触面积来提高脱硫效率[4].针对这一问题, 结合KR法和喷吹法的优点, 张廷安等提出了新型脱硫技术——“机械搅拌原位脱硫”的新工艺[5-6]; 该工艺利用载气携带镁基脱硫剂并将其喷吹进入铁水中, 在铁水中原位反应产生镁蒸气.经过搅拌桨搅拌, 镁蒸气的大气泡被破碎为小气泡, 并均匀广泛地分布在铁水中, 与铁水中的硫反应, 达到脱硫目的.刘燕等[7-8]通过水模型实验, 研究了桨型、搅拌转速、搅拌桨浸入深度等因素对铁水脱硫气泡行为特性的影响, 得到了最佳工艺条件.
为了研究气泡直径与操作条件的准数方程, 刘燕等[8]研究了因次分析法在水模型实验中的应用.运用因次分析得到了氧气底吹炼铜过程中气泡直径与操作参数之间的准数方程[9]和加压搅拌浸出体系下氧气的含气率与操作参数之间的准数方程[10].本文通过水模型实验研究了铁水脱硫过程中气泡直径与各影响因素之间的关系, 得到了相关准数方程.
1 实验本实验采用相似原理建立了某钢厂120t铁水包的有机玻璃圆柱筒水模型实验设备.设备直径为435mm, 高为515mm, 液面高度为350mm.由于水与铁水的动力黏度相似, 因此用水代替铁水, 用空气代替铁水中的镁蒸气, 以此模拟镁蒸气在铁水中的分布状态.实验过程中采用高速照相机(500帧/s)采集熔池中的气泡分布瞬时状态, 用图像处理软件Image Pro-Plus对采集的照片进行统计处理, 从而得到气泡的平均直径.将得到的气泡平均直径代入准数方程, 得到准数方程中的相关系数, 最终确定气泡平均直径的准数方程.研究气泡平均直径与设备参数、操作条件之间的关系, 对指导放大实验过程、进行工业生产具有重要意义.
2 气泡平均直径准数方程的建立及系数确定 2.1 气泡平均直径准数方程的建立由实验可知, 影响气泡平均直径dB的主要因素有搅拌转速n、通气流量Q、搅拌桨浸入深度h、偏心度e等.由文献[9]可知气泡平均直径dB也与容器直径D、气体密度ρg、气体出口速度v、气体黏度μg、液体密度ρl、液体黏度μl、液面高度H、气液表面张力σg-l、重力加速度g、入射气体压强p等因素有关.由因次分析原理可以得出一般的函数形式为
|
(1) |
上述变量的因次如表 1所示.
|
表 1 变量因次表 Table 1 Variables in dimensional analysis |
总变量数N=15, 独立变量数k=3, 由π定理可建立N-k=12个无因次组合.选取ρl, D, v为独立变量.对于只含长度因次的变量dB, h和H, 在构造无因次∏时可以直接用独立变量D表示.e为无因次量, 可直接写出.经过计算, 各个∏分别表示为
|
(2) |
|
(3) |
|
(4) |
|
(5) |
|
(6) |
|
(7) |
|
(8) |
|
(9) |
|
(10) |
|
(11) |
|
(12) |
|
(13) |
对于∏1, 代入相关变量的因次, 可得关系式:
|
(14) |
由此可得方程组:
|
(15) |
解得a1=0, b1=1, c1=-1.
故
|
(16) |
|
(17) |
|
(18) |
式中Re′为相对雷诺数.
|
(19) |
式中Re为雷诺数, 可用来表征流体流动情况的无量纲数.
|
(20) |
式中We为韦伯数, 代表惯性力和表面张力效应之比.
|
(21) |
式中Eu为欧拉数, 描述动量传递的特征数, 体现了在流动过程中动量损失率的相对大小.
|
(22) |
式中Fr为弗劳德数, 表示重力对流动影响的准数.
将各因次代入式(1)可得
|
(23) |
通过分析本文的研究内容, 可知主要影响气泡的是We数, 因此可以忽略Re′, Re, Eu和Fr.式(23)可以简化为
|
(24) |
此外, 实验中ρg, ρl和H/D为定量, 于是又可得到
|
(25) |
式(25)中
根据因次分析的特点, 可认为不同搅拌桨桨型和喷嘴结构下因次公式的形式是相同的, 不同的是具体经验公式的结果.各种现象的准则关系式在自变量介于某一范围内往往都可以采用幂函数的形式表述, 因次可以拟合出经验准则关系式:
|
(26) |
式中A, B, C, D′, E和F为拟合系数.
将式(26)两边取对数得
|
(27) |
依据上述线性关系对实验数据进行处理, 可以得到拟合系数A, B, C, D′, E和F.
对于偏心搅拌, 因为We=ρlDv2/σg-l, 式中各变量都恒定, 因此We为定值, 所以式(27)中lnA, FlnWe为定值, 这些定值可以合并为一项G=lnA+FlnWe, G也是恒定的.
式(27)可以改写为
|
(28) |
为了求拟合系数B, 在搅拌转速n变化的同时必须保证其他量不变, 因此式(28)可简化为
|
(29) |
将镁蒸气脱硫的水模型实验数据线性拟合, 可以得到B=－0.07, H=－5.676, 拟合曲线如图 1所示.
|
图 1 ln(dB/D)与ln(Dn/v)的拟合曲线 Fig.1 Fitting curve between ln(dB/D) and ln(Dn/v) |
同理, 由图 2~图 4依次可以得到C=0.045, D′=-0.243, E=-0.358, G=-5.481, 则经验公式为
|
(30) |
|
图 2 ln(dB/D)与ln(Q/vD2)的拟合曲线 Fig.2 Fitting curve between ln(dB/D) and ln(Q/vD2) |
|
图 3 ln(dB/D)与ln(h/D)的拟合曲线 Fig.3 Fitting curve between ln(dB/D) and ln(h/D) |
|
图 4 ln(dB/D)与lne的拟合曲线 Fig.4 Fitting curve between ln(dB/D) and lne |
取搅拌桨偏心度e=0.4, 搅拌桨浸入深度h=250mm, 通气流量分别为1.0, 1.5, 2.0, 2.5m3/h时, dB与搅拌转速n之间的关系可以表示为
|
(31) |
将Q值代入式(31)得到不同通气流量下, 搅拌转速对气泡平均直径的影响结果如图 5所示.
|
图 5 不同通气流量和搅拌转速下的气泡平均直径 Fig.5 Effect of rotation speed on average bubble diameter under different gas flow rates |
从图 5中可以看出, 在通气流量相同时, 随着搅拌转速的增加, 气泡的平均直径减小.这是因为在通气流量一定时, 增大搅拌转速, 单位时间内熔池的搅拌次数增多, 形成的剪应力也就越大, 对气泡的破坏也越强, 因此气泡的平均直径越小.从图 5还可以看出, 当搅拌转速一定时, 随着通气流量增加, 熔池内的气泡平均直径也增大.
3.2 不同偏心度下搅拌转速对气泡平均直径的影响取Q=1.0m3/h, 搅拌桨浸入深度h=250mm, 搅拌桨偏心度分别为0.2, 0.3和0.4时, dB与搅拌转速n之间的关系可以表示为
|
(32) |
将e值代入式(32)得到不同偏心度下, 搅拌转速对气泡平均直径的影响如图 6所示.可以看出, 偏心度相同时, 随着搅拌转速的增加, 气泡平均直径减小.在相同的搅拌转速下, 随着偏心度的增加, 气泡的平均直径逐渐减小.因为当搅拌转速一定时, 随着搅拌桨偏心度增大, 搅拌桨产生的漩涡远离了桨轴, 形成的切向流很不稳定, 在熔池中形成了很强的紊流状态, 阻止了气泡的聚并, 同时又进一步促进了气泡的破碎.
|
图 6 不同偏心度和搅拌转速下的气泡平均直径 Fig.6 Effect of rotation speed on average bubble diameter under different impeller eccentricities |
取Q=1.0m3/h, 偏心度e=0.4, 搅拌桨浸入深度分别为220, 250和280mm时, dB与搅拌转速n之间的关系可以表示为
|
(33) |
将h值代入式(33)得到不同搅拌桨浸入深度下, 搅拌转速对气泡平均直径的影响如图 7所示.
|
图 7 不同搅拌桨浸入深度和搅拌转速下的气泡平均直径 Fig.7 Effect of rotation speed on average bubble diameter under different impeller immersion depths |
从图 7中可以看出, 在搅拌桨浸入深度相同时, 随着搅拌转速的增加, 气泡平均直径减小.在相同的搅拌转速下, 随着搅拌桨浸入深度的增加, 气泡的平均直径也减小.这是通过因次分析, 只计算了气泡的平均直径, 忽略了随着搅拌桨浸入深度的增加, 熔池内部分气泡聚合的情况.
4 结论1) 增大通气流量后, 由于气泡之间的聚合作用, 气泡的平均直径随之增大, 降低了气泡的比表面积和气液反应的接触面积, 不利于提高脱硫效率.
2) 增加搅拌转速、搅拌桨偏心度和搅拌桨浸入深度均有利于减小气泡平均直径，并有利于气泡在熔池中的均匀分布, 从而提高脱硫效率.
3) 通过对水模型实验数据的拟合, 得出气泡平均直径与操作参数之间的准数方程为
|
| [1] |
Liu Y, Zhang T A, Sano M, et al.
Mechanical stirring for highly efficient gas injection refining[J]. Transactions of Nonferrous Metals Society of China, 2011, 21(8): 1896–1904.
DOI:10.1016/S1003-6326(11)60947-3 |
| [2] |
Liu Y, Sano M, Zhang T A, et al.
Intensification of bubble disintegration and dispersion by mechanical stirring in gas injection refining[J]. ISIJ International, 2009, 49(1): 17–23.
DOI:10.2355/isijinternational.49.17 |
| [3] |
Liu Y, Zhang Z M, Liu J N, et al.
Cold model study on Mg desulfurization of hot metal under mechanical stirring[J]. Journal of Iron and Steel Research, International, 2013, 20(12): 1–6.
DOI:10.1016/S1006-706X(13)60208-3 |
| [4] |
Liu Y, Zhang T A, Sano M, et al.
Study on absorption rate by eccentric mechanical stirring in gas injection refining for iron and steel making[J]. Journal of Iron and Steel Research, International, 2011, 18(sup2): 166–171.
|
| [5] |
张廷安, 刘燕, 赫冀成, 等.原位机械搅拌炉外脱硫方法及装置: ZL20081001186303.3[P]. 2008-10-03.
( Zhang Ting-an, Liu Yan, He Ji-cheng, et al.External furnace desulfurization method and device with in-situ mechanical stirring: ZL20081001186303.3[P]. 2008-10-03. ) |
| [6] |
Liu Y, Shao P, Zhao H L, et al.
Simulation on gas injection refining process with mechanical stirring[J]. Telkomnika Indonesian Journal of Electrical Engineering, 2012, 10(8): 2194–2201.
|
| [7] |
Liu Y, Zhang Z M, Sano M, et al.
Improvement of impeller blade structure for gas injection refining under mechanical stirring[J]. Journal of Iron and Steel Research, International, 2014, 21(2): 135–143.
DOI:10.1016/S1006-706X(14)60022-4 |
| [8] |
刘燕, 张廷安, 王强, 等.
因次分析法在水模型实验中的应用[J]. 工业炉, 2007, 29(6): 14–17.
( Liu Yan, Zhang Ting-an, Wang Qiang, et al. Dimension analysis theory in the application of fluid mechanics[J]. Industrial Furnace, 2007, 29(6): 14–17. ) |
| [9] |
Wang D X, Liu Y, Zhang Z M, et al.
Dimensional analysis of average diameter of bubbles for bottom blown oxygen copper furnace[J]. Mathematical Problems in Engineering, 2016, 2016: 1–8.
|
| [10] |
田磊, 刘燕, 吕国志, 等.
加压搅拌浸出体系下氧气的气含率[J]. 中国有色金属学报, 2017, 27(3): 655–665.
( Tian Lei, Liu Yan, Lyu Guo-zhi, et al. Gas holdup in pressure agitation leaching system[J]. The Chinese Journal of Nonferrous Metals, 2017, 27(3): 655–665. ) |