張欣利，張文懋，喬通，張世雄，趙定國
( 華北理工大學冶金與能源學院，河北唐山063210)
摘要: 高壓底吹條件下鋼液流動運動行為復雜，采用數值模擬分析了壓力、底吹流量、底吹位置及鋼液面高度對熔池內鋼液流動的影響。模擬結果表明: 高壓條件下，隨著壓力的增加，熔池內的速度分布均勻，梯度減小程度小。鋼液的平均速度隨底吹流量的增加而增加。隨底吹噴孔位置向熔池壁靠近，主循環回路變長。隨著熔池液面高度的上升，氣液兩相區的徑向范圍變廣，速度梯度減小，熔池內死區比例降低; 達到一定高度后，氮氣到達熔池液面時速度變小。研究高壓底吹條件下鋼液的流動行為，對鋼包內鋼液底吹精煉具有指導意義。
關鍵詞: 高壓底吹氮氣; 鋼液流動; 數值模擬
0 引言
常壓下鋼包底吹氣體的數值模擬和物理模擬已有大量的研究，但對高壓條件下鋼包內流場的數理模擬較少，高壓底吹氣體對鋼液流場、湍動能等的影響規律還未見定論。
煉鋼過程中，噴吹氣體攪拌鋼液可以均勻成分和溫度，加快傳質。因此，冶金專家學者對中間包內的流場［1 － 3］、RH 真空精煉過程的流場［4 － 6］、結晶器內的流場［7 － 9］等進行了大量的研究。對鋼液流場的研究一般是處于常壓和真空狀態下的，對高壓條件下的流場的研究較少。研究高壓條件下流場的變化對高壓冶金的發展具有指導意義。
1 研究方案設計
采用Fluent 數值模擬軟件對高壓底吹氮氣時鋼液的流場進行模擬，設計選用L9( 34 ) 的正交表進行流場模擬實驗，正交實驗的因素與水平如表1 所示，模擬條件如圖1 所示。
按照設計方案進行正交實驗，高壓底吹氮氣鋼液流場影響因素的正交實驗設計如表2 所示。
2 模型的建立
對熔池內高壓底吹氮流場的數學模擬有以下假設: 高壓底吹氮過程中，驅使鋼液循環流動的主動力為氣泡浮力; 流體設定為不可壓縮的粘性流體，鋼液密度為常數; 流動為定常流動; 熔池壁面、底面光滑，頂部液面為光滑的自由面; 氣泡是大小均勻，具有同一直徑的球體; 運動過程中不發生化學反應; 忽略溫度對鋼液流場的影響。
用Gambit 建立熔池模型，劃分網格，如圖2( a) 、2( b) 所示，該偏心底吹模型共67 465 個網格，將Gambit 劃分的網格導入Fluent 中，按正交實驗設計進行運算。底吹噴嘴位于X 軸距熔池中心a mm處，鋼液面距離熔池底部b mm。取熔池的3 個主要截面進行分析，分別為x = a 的A 截面，y = 0 的B截面，z = b 的C 截面，如圖2( c) 所示。
3 高壓底吹鋼液流場模擬結果與分析
按表2 方案進行正交模擬實驗，得到高壓底吹條件下B 截面速度矢量圖，如圖3 所示。
由圖3 可知，高壓底吹條件下熔池底吹氮氣時，氮氣自熔池底噴嘴進入熔池內，帶動熔池內鋼液向上流動形成上升流股。當氣液兩相區到達熔池液面時，氣液分離，氮氣進入熔池上部的高壓氮氣氛，鋼液在慣性力和重力的共同作用下，形成流向熔池壁的表面流。鋼液形成的表面流到達熔池壁后，在高壓作用和熔池壁的阻擋下轉而沿側壁向熔池底部流動。在氣液兩相區的卷吸作用下，再次流入氣液兩相區，形成循環流。此流場形成的主要驅動力是慣性力和氣泡浮力，在其驅使下氣液兩相區的中心速度最大，速度由中心沿徑向方向逐漸減小，循環流中心速度接近為零。高壓底吹條件下的B 截面速度分布云圖，如圖4 所示。
由圖4 可知，隨著底吹速度的增加，氣液兩相區的上升速度增加，增強了對氣液兩相區附近鋼液的卷吸能力，氣液兩相區附近同一位置的鋼液向上流動的速度增加，使得氣液兩相區在同一高度的橫截面積增加，帶動周圍鋼液上升的范圍增加，熔池頂部及底部的鋼液速度也明顯增大，使得熔池內部整體的速度隨之增加，速度分布更均勻，速度梯度減小，有利于氮氣對鋼液的攪拌和混勻效果，減小死區。同時，瞬時吹入熔池內的氮氣泡數量隨底吹速度的增加而增加，氮氣泡與熔池內鋼液的總接觸面積增加，有利于熔池內氮氣向鋼液中傳質的效果。
對比圖3( c) 和圖4( c) 可知，當底吹噴孔位置位于偏心21 mm 時，上升流股很快與熔池壁發生接觸，使得大型氣液兩相區的橫截面積減小，靠近熔池壁一側的小循環流消失。結合圖3( a) 、3( b) 、3( c)和圖4( a) 、4( b) 、4( c) 可知，隨著底吹噴孔位置向熔池壁的偏移，流場中氣液兩相區遠離熔池壁一側主循環流的回路變長，速度分布均勻，速度梯度減小。由此表明，底吹噴孔位置的偏移，使得小循環流消失，主循環回路變長，有利于氮氣對鋼液的傳質和混勻效果，減小了鋼液死區比例。
4 結論
(1) 高壓條件下，鋼液的平均速度隨底吹流量的增加而增加。高壓條件下，隨底吹噴孔位置向熔池壁靠近，使得小循環流消失，主循環回路變長，有利于鋼液的混勻和傳質; 過于偏心底吹，氣液兩相區與熔池壁的碰撞，造成能量損失和熔池壁材料的損耗。
(2) 高壓條件下，隨著熔池液面高度的上升，氣液兩相區的徑向范圍變廣，速度梯度減小，熔池內死區比例降低。達到一定高度后，氮氣到達熔池液面時速度變小，動能變少。
參考文獻
［1］張美杰，汪厚植，顧華志，等． 氣幕擋墻對中間包內鋼液流場影響的數值模擬［J］． 鋼鐵研究學報，2006，18( 6) : 17 ～ 21．
［2］ 韓建剛． 宣鋼十二流中間包流場水模實驗研究［J］． 河北冶金，2016，( 1) : 19 ～ 23．
［3］ 張美杰，黃奧，汪厚植，等． 底吹中間包的流場及夾雜物運動軌跡的數值模擬研究［C］/ / 2005 中國鋼鐵年會． 2005．
［4］ 孫亮，艾立群，趙俊花，等． RH 精煉過程鋼液流動行為與循環流量數值模擬［J］． 鋼鐵釩鈦，2009，30( 2) : 28 ～ 32．
［5］Miki Y，Shimada Y，Thomas B G，et al． Model of Inclusion Removal during RH Degassing of Steel［J］． Iron ＆ Steelmaker，1997．
［6］ 張琳，孫彥輝，朱進鋒，等． ＲH 精煉過程鋼液流動數值模擬和應用［J］． 北京科技大學學報，2009，31( 7) : 821 ～ 825．
［7］鄭麗艷，耿賽曉． 連鑄坯結晶器內鋼液流動數學模擬研究［J］． 河北冶金，2014，( 6) : 1 ～ 6．
［8］ 竇冠琪，王寶峰，李建超，等． 寬厚板坯連鑄結晶器內的鋼液流場的數值模擬［J］． 連鑄，2013，1( 1) : 27 ～ 32．
［9］ 陳永，張長利，沈厚發，等． 連鑄結晶器電磁攪拌磁場及鋼液流場模擬［J］． 鋼鐵研究學報，2008，20( 7) : 11 ～ 14．