连铸板坯结晶器内钢液吹氩行为的数值模拟.pdf_大学文库

D0I:10.13374/j.issn1001053x.2003.03.032 第25卷第3期北京科技大学学报 Vol.25 No.3 2003年6月 Journal of University of Science and Technology Beijing Jun.2003 连铸板坯结晶器内钢液吹氩行为的数值模拟于会香朱国森王新华张炯明王万军北京科技大学冶金学院，北京100083 摘要以宝钢一连铸板坯结晶器为研究对象，采用多相流模型计算了吹氩后结晶器内钢液的流场、温度场及氩气分布，结果表明：吹氩后，结晶器内上回流区钢液的流动强度增大，下回流区变小：钢液的温度梯度减小，温度分布均匀：氩气的分布不均匀，上回流区钢液的含气率大大高于下回流区关键词连铸；结晶器；吹氢：数值模拟分类号TF777.1 浸入式水口吹氩保护浇铸技术是指在浸人(4)不考虑结晶器振动及锥度等因素的影响；（⑤）式水口中，通过塞棒吹入氩气防止空气吸人的技氩气泡尺寸均匀一致，氩气温度与钢液的温度相术.该技术不仅能够有效防止水口堵塞，而且能同：(6)水口人口处氩气分布均匀. 强化搅拌，均匀钢水成分和温度，促进夹杂物上浮 1.2多相流模型以及改善铸坯质量.但是吹人氩气也会带来负面在多相流模型中，每一相有自己的一套求解影响，如导致保护渣卷入、钢水的裸露氧化及皮方程，各相的传输量通过相间传输项发生联系. 下缺陷等.因此，研究水口吹氩后结晶器内的钢对于气液两相流，气相和液相分别有自己的质量液的行为十分必要. 守恒方程、N—S方程、能量守恒方程以及x一ε双近年来许多学者采用商业软件CFX和数值方程，两相的速度和温度通过相间力和热传递项模拟法对这一问题进行了大量研究.对于结产生联系.其方程如下：晶器内的液相，建立质量守恒、动量守恒和K一￡ 0i(r.p.0.)+V-[r.(p.U..-r.Vo.)]- 双方程，在动量守恒方程中加人氩气泡的浮力项；对于结晶器内的气相，只求解质量守恒方程， rS.+cd0r-0.)+(m.D,-m.D.) 认为气泡一进入计算区域，便在垂直方向上达到其中，c(D。-中)为a,B相之间中的相间传输，c。= 0,Cog=Caa 了其稳态的最终速度，该方法与实际情况相比简 13计算条件化较大.本文则采用多相流模型对吹氩后钢液的 (1)边界条件.①模型的入口设在浸入式水口行为进行模拟，即为液相和气相分别建立一套求的人口处，其速度值根据入口与出口的质量守恒解方程，各相的传输量通过相间传输项发生联计算得出，温度T=1803K.②模型的出口设在结系. 晶器下方的出口处，采用质量边界条件，即出口 1数学模型处与入口处质量守恒，③模型的壁面设为达到但不包括糊状区的位置，即壁面处的温度设为恒定 1.1模型假设的温度1775K:壁面采用无滑移边界条件. 在建立多相流模型时本文采用了以下假设： (2)其他参数.计算所用工艺参数为：板坯断 (1)结晶器内钢液的流动为稳态粘性不可压缩流面尺寸1200mm×250mm,拉坯速度1.3m/min,水动；(2)结晶器液面设置为自由面，氩气完全由月口侧孔倾角-15°，水口浸人深度220mm,氩气泡由表面逸出：（③）不考虑结晶器内凝固壳的存在；直径3mm,水口人口处的含气率4%. 钢水和氩气的物理参数见表1· 收稿日期200206-21于会香女，24岁，助教

第卷第期年月北京科技大学学报连铸板坯结晶器内钢液吹氨行为的数值模拟于会香朱国森王新华张炯明王万军北京科技大学冶金学院，北京摘要以宝钢一连铸板坯结晶器为研究对象，采用多相流模型计算了吹氨后结晶器内钢液的流场、温度场及氢气分布，结果表明吹氢后，结晶器内上回流区钢液的流动强度增大，下回流区变小钢液的温度梯度减小，温度分布均匀氢气的分布不均匀，上回流区钢液的含气率大大高于下回流区关键词连铸结晶器吹氢数值模拟分类号浸人式水口吹氢保护浇铸技术是指在浸人式水口中，通过塞棒吹人氢气防止空气吸人的技术该技术不仅能够有效防止水口堵塞，而且能强化搅拌，均匀钢水成分和温度，促进夹杂物上浮以及改善铸坯质量但是吹人氢气也会带来负面影响，如导致保护渣卷人、钢水的裸露氧化及皮下缺陷等因此，研究水口吹氢后结晶器内的钢液的行为十分必要近年来许多学者采用商业软件和数值模拟法对这一问题进行了大量研究卜对于结晶器内的液相，建立质量守恒、动量守恒和双方程，在动量守恒方程中加人氢气泡的浮力项对于结晶器内的气相，只求解质量守恒方程，认为气泡一进人计算区域，便在垂直方向上达到了其稳态的最终速度该方法与实际情况相比简化较大本文则采用多相流模型对吹氢后钢液的行为进行模拟，即为液相和气相分别建立一套求解方程，各相的传输量通过相间传输项发生联系不考虑结晶器振动及锥度等因素的影响氢气泡尺寸均匀一致，氢气温度与钢液的温度相同水口人口处氢气分布均匀多相流模型在多相流模型中，每一相有自己的一套求解方程，各相的传输量通过相间传输项发生联系对于气液两相流，气相和液相分别有自己的质量守恒方程、萍一方程、能量守恒方程以及双方程，两相的速度和温度通过相间力和热传递项产生联系其方程如下景。，。。 · 。。二成一甲。数学模型模型假设在建立多相流模型时本文采用了以下假设结晶器内钢液的流动为稳态粘性不可压缩流动结晶器液面设置为自由面，氢气完全由自由表面逸出不考虑结晶器内凝固壳的存在收稿日期刁于会香女，岁，助教呱嗯以呜一叫噜鸿一，叫其中，试妈一叭为渭相之间必的相间传输，。，叨 · 计算条件边界条件 ①模型的人口设在浸人式水口的人口处，其速度值根据人口与出口的质量守恒计算得出，温度 ②模型的出口设在结晶器下方的出口处，采用质量边界条件，即出口处与人口处质量守恒 ③模型的壁面设为达到但不包括糊状区的位置，即壁面处的温度设为恒定的温度壁面采用无滑移边界条件其他参数计算所用工艺参数为板坯断面尺寸，拉坯速度，水口侧孔倾角一，水口浸人深度，氢气泡直径，水口人口处的含气率钢水和氢气的物理参数见表 DOI ：10．13374／j ．issn1001－053x．2003．03．032

·216· 北京科技大学学报 2003年第3期表1物性参数 Table 1 Physical parameters 材料密度/(kgm) 分子粘性系数/(kgm's)导热系数(WmK)恒压热容/kgK) 钢水 7020 5.5×103 34 680 氩气 0.56 7.42×10-5 0.017 523.38 1.4数值求解表面上温度很高，自由表面活跃；温度场的不对采用英国AEA公司提供的CFX4.3商业软称性有所增加. 件，对上述模型进行了计算.求解在PentiumⅢ 23水口吹氩后结晶器内氩气分布的基本特征 -550计算机上进行，计算时间约为300min,迭代图3为吹氩时水口与结晶器内钢液含气率的 1.6×10次以后，计算达到收敛，质量残差为10 分布情况.可见，水口内钢液含气率的分布情况为：水口内钢液含气率较高，在4.4%以上；水口 2计算结果与讨论出口上半部分钢液的含气率最高，最大值达 13.27%;水口底部的“凹”形区内，含气率几乎为 21水口吹氩后结晶器内钢液流场的基本特征图1(a,b)分别为不吹氩和吹氩时结晶器内 0.结晶器内钢液含气率的分布情况为：氩气主要集中在水口出口上半部分，上表面靠近水口壁面钢液的流场.可见，吹氩前后结晶器内钢液的流附近，水口底部外壁，上下回流区的涡心处；结晶动相似，基本特征为：(1)钢液从水口流出后，形器下方钢液的含气率很低，在10数量级. 成一股很强的射流，撞击到结晶器窄面后，分成向上和向下两股钢流.(2)向上的流股在自由表面附近形成两个较强回流，向下的流股也形成与上部回流方向相反、范围更大的两个回流，其强度在向下延伸过程中逐渐减弱.但是，吹氩后结晶器内钢液的上回流强度增强，下回流强度明显减噻弱，下回流的涡心明显上移，下回流区减小.主要原因是氩气在吹入结晶器后对钢液有提升作用，致使钢液向下的速度减小，向上的速度增加. yu=1.3236 Ven=1.3376 (a)不吹氩(b)水口入口处含气率为4% 图2不吹氩与吹氮时结晶器内钢液的温度场图 Fig.2 Simulated temperature fields without and with argon injection 662是 (a)不吹氩 b)水口入口处含气率为4% 图1不吹氩与吹氩时结晶器内钢液的流场图 774720 Fig.1 Simulated velocity fields without and with argon in- jection 22水口吹盒后结晶器内钢液温度场的基本特征图2(a,b)分别为不吹氩四和吹氩时结晶器内钢液的温度场.可见，与不吹氩时相比，吹氩后结图3吹氩时漫入式水口与结晶器内含气率的分布图晶器内钢液的温度场发生了较大的变化：钢液的 Fig.3 Gas distribution in the submerged entry nozzle and 温度普遍增高，温度梯度小，分布趋于均匀；自由 mold when argon injection

一北京科技大学学报年第期表物性参数材料钢水氢气密度伙 · 一分子粘性系数 · 一 ’ 一 ’ 导热系数 · 一 ’ 一 ‘ 恒压热容 · 一 ’ 一，一，数值求解采用英国公司提供的商业软件，对上述模型进行了计算求解在一计算机上进行，计算时间约为，迭代、次以后，计算达到收敛，质量残差为一计算结果与讨论水口吹盔后结晶器内钢液流场的基本特征图分别为不吹氢 ‘，” 和吹氢时结晶器内钢液的流场可见，吹氢前后结晶器内钢液的流动相似，基本特征为钢液从水口流出后，形成一股很强的射流，撞击到结晶器窄面后，分成向上和向下两股钢流向上的流股在自由表面附近形成两个较强回流，向下的流股也形成与上部回流方向相反、范围更大的两个回流，其强度在向下延伸过程中逐渐减弱但是，吹氢后结晶器内钢液的上回流强度增强，下回流强度明显减弱，下回流的涡心明显上移，下回流区减小主要原因是氢气在吹人结晶器后对钢液有提升作用，致使钢液向下的速度减小，向上的速度增加表面上温度很高，自由表面活跃温度场的不对称性有所增加水口吹氛后结晶器内氨气分布的基本特征图为吹氢时水口与结晶器内钢液含气率的分布情况可见，水口内钢液含气率的分布情况为水口内钢液含气率较高，在以上水口出口上半部分钢液的含气率最高，最大值达水口底部的 “ 凹 ” 形区内，含气率几乎为结晶器内钢液含气率的分布情况为氢气主要集中在水口出口上半部分，上表面靠近水口壁面附近，水口底部外壁，上下回流区的涡心处结晶器下方钢液的含气率很低，在一数量级卜卜泪刀力潜吸肠脚缓脚朋加用蒯蓄象蓄备誉哥蜜誉聋考鑫象参象吞戴苏爹潇蓄澡健落滋罐通准通灌健落姜蓄略噜滚展傲蓉赘探灌壤灌落噶礴鹰凄着蘑孰、、、乐一介少价、阵 ‘ 厂一雄摧黔 “ 丫愧衡一不吹氢水口人口处含气率为图不吹鼠与吹盗时结晶器内钢液的温度场图一门了〔曰临盯之曰亡，了巴鬓麒一攫创 … 」勺图矽交矛甲甲，岁女，冷亨亏个兮七令丫‘亏卜告今杏，老俐今卜，护，犷、，宁心‘才、资平，‘ ‘石咨，孕念，勺俄斤啥李，牛才心不于今今奋甲十卜令‘ ‘飞，冲︺，平侮琳喻丫凌哈甲叹去壬毒，小守争争之奋少启七今心几母瑞﹂卜不吹氢水口人口处含气率为图不吹盘与吹盘时结晶器内钢液的流场图啥 · 。水口吹盆后结晶器内钢液温度场的基本特征图，分别为不吹氢 ‘川和吹氢时结晶器内钢液的温度场可见，与不吹氢时相比，吹氢后结晶器内钢液的温度场发生了较大的变化钢液的温度普遍增高，温度梯度小，分布趋于均匀自由吹益时浸入式水口与结晶器内含气率的分布图卜

Vol.25 No.3 于会香等：连铸板坏结晶器内钢液吹氩行为的数值模拟 ·217· 3结论 in a continuous slab casting mold [A].Steelmaking Con- ference Proceedings [C].America,1993.273 (1)将大型商业软件C℉X4.3应用于冶金领 4 Bai Hua,Thomas B G.Effect of clogging,argon injection 域，采用多相流模型(multi-phase flow model)模拟 and casting conditions on flow rate and air aspiration in 计算了水口吹氩时结晶器内钢液的流场、温度场 submerged entry nozzles [A].Steelmaking Conference 及氩气分布. Proceedings [C].America,2000.183 (2)吹氩后结晶器内钢液流场与温度场的基 5 Li Baokuan,Okane Toshimitsu,Umeda Takateru.Model- 本特征为：流场由上下四个回流区组成，上回流 ing of molten metal flow in a continuous casting process considering the effects of argon gas injection and static 的流动强度增强，下回流强度减弱，下回流区域 magnetic-field application[J].Metall Mater Trans B,2000, 减小.温度梯度减小，分布均匀. 31B:1491 (3)吹氩后含气率分布的基本特征为：氩气分 6 Zhang Jiongming,Wang Wanjun,He Jicheng.Effect of 布不均匀，主要分布在水口出口上半部分、壁面 gas injection on fluid flow in mold [J],J Univ Sci Technol 附近、底部壁面附近和结晶器上下回流区的涡心 Beijing.I998,5(1):16 处；水口底部“凹”槽中和结晶器下方钢液的含气 7 Xu D,Jones WK Jr,Evans J W.Experimental investiga- 率很低 tion of the argon bubble behavior in the mold during con- tinuous casting of steel-A water model study [A].Steel- 参考文献 making Conference Proceedings [C].America,1999.155 1 Bai Hua,Thomas B G.Turbulent flow of liquid steel and 8于会香，张炯明，王万军，等.板坯连铸浸入式水口 argon bubbles slide-gate tundish nozzles:Part I Model de- 出口速度对结晶器流场的数值模拟[J刀.北京科技大 velopment and validation [J].Metall Mater Trans B. 学学报，2002,245：492 2001,32B:253 9马范军，文光华，李刚.板坯连铸结晶器内吹入气体 2 Bai Hua,Thomas B G.Turbulent flow of liquid steel and 对钢液行为的影响[.炼钢，2000,16(3)：42 argon bubbles slide-gate tundish nozzles:Part II Effect of 10王妍，郑炜，朱苗勇.宝钢一炼钢厂连铸机结晶器浸 operation conditions and nozzle design [J].Metall Mater 人式水口吹氩的水模研究[].宝钢技术，2000,6：26 Trans B,2001,32B:269 11于会香.连铸板坯结晶器内浸入式水口的优化研究 3 Thomas B G,Huang X.Effect of argon gas on fluid flow [D1北京：北京科技大学，2002 Numerical Simulation of Fluid Behavior in Continuous Slab-Casting Mold with Argon Injection YU Huixiang,ZHU guosen,WANG Xinhua.ZHANG Jiongming,WANG Wanjun Metallurgical Engineering School,University of Science and Technology Beijing,Beijing 100083,China ABSTRACT Taking the mold of No.I continuous slab-casting machine in Baosteel as the object,a multi-phase model was applied to simulate the velocity,temperature fields and gas distribution of molten steel with argon injec- tion.The results indicate that the velocity of upward flow increases and the lower re-circulation region appears small,the gradient of temperature decreases.Argon bubbles distribute unevenly and the gas holdup of liquid steel in the upper region is far higher than that of the lower region. KEY WORDS continuous casting;mold;argon injection;numerical simulation