韓永輝，蒼大強，高 海，劉福海，張興利

（1.北京科技大學冶金與生態(tài)工程學院，北京 100083；2.河北鋼鐵股份有限公司承德分公司，河北承德 067002）

在以釩鈦磁鐵礦為原料的承德鋼廠（以下簡稱承鋼），鋼包作為轉(zhuǎn)爐提釩和轉(zhuǎn)爐精煉的中間容器，是煉鋼過程中必不可少的設(shè)備［1-2］。在不改變工藝的前提下，在鋼包中進行脫磷是一種經(jīng)濟實用、簡便易行的新方法［3］。但由于轉(zhuǎn)爐提釩后出半鋼時間較短，動力條件成為影響脫磷效率的主要瓶頸［4-11］。為此，本文利用商業(yè)軟件Fluent對鋼包內(nèi)的底吹模式及側(cè)吹氧槍的噴吹射流性質(zhì)進行數(shù)值模擬計算，研究不同側(cè)吹條件下熔池的流動特性，以期為現(xiàn)場脫磷生產(chǎn)工藝提供理論依據(jù)。

1 模型建立

1.1 模擬方案

根據(jù)水模型實驗結(jié)果，底吹最佳位置是0.7R-120°，最佳的底吹流量是100 L／min。最佳的頂吹方案為頂吹流量3 000 Nm3／h、頂槍角度13°、頂槍高度500 mm。對承鋼鋼包采用3種噴吹方式進行模擬，并且模型中將底部透氣磚簡化為4根管狀噴嘴，且考慮渣層的影響，設(shè)計的模擬方案如下。

方案1：鋼包氧槍側(cè)吹氣量為2 500 Nm3／h；

方案2：鋼包氧槍側(cè)吹氣量為2 800 Nm3／h；

方案3：鋼包氧槍側(cè)吹氣量為3 000 Nm3／h。

1.2 幾何模型及網(wǎng)格劃分

1.2.1 幾何模型

采用1∶1比例進行建模計算，模型計算區(qū)域包括底吹噴管、鋼包內(nèi)熔池以及上部凈空，鋼包模型幾何尺寸參照承鋼100 t鋼包設(shè)計圖紙尺寸。鋼包物理模型的主要參數(shù)詳見表1。鋼包主體為類圓錐結(jié)構(gòu)，如圖1所示。

表1 鋼包物理模型的主要參數(shù)Tab.1 Main parameters of physical model of ladle

圖1 鋼包模型及底吹點示意圖Fig.1 Schematic diagram of ladle model and bottom blowing point

1.2.2 網(wǎng)格模型

建模時，利用Fluent軟件包中的Gambit模塊對建立好的幾何模型進行網(wǎng)格劃分［12］?？臻g域網(wǎng)格的劃分采用內(nèi)密外疏、下密上疏劃分法。最小網(wǎng)格尺寸位于氣體入口處，計算域中最大網(wǎng)格尺寸為18 mm，網(wǎng)格數(shù)約為13萬，網(wǎng)格質(zhì)量為0.4以上。本模擬以整個實際鋼包尺寸為計算域進行計算，包括頂部的汽液兩相區(qū)，鋼包網(wǎng)格形式見圖2。

圖2 鋼包模型網(wǎng)格圖Fig.2 Grid of ladle model

1.3 控制方程及邊界條件

假設(shè)底吹噴管與鋼包內(nèi)部所有連接處都很光滑，忽略摩擦作用，氧槍壁面為絕熱面。射流運動過程中的控制方程［13-16］如下。

連續(xù)性方程：

動量守恒方程：

能量守恒方程：

對于湍流流動，引入標準的κ-ε方程。

κ方程：

ε方程：

為研究鋼包復吹攪拌效果是否能夠達到設(shè)計標準，研究利用Fluent模擬軟件采用1∶1比例模型對不同底吹情況下鋼包流場的影響進行模擬，利用VOF模型模擬多相流的相互影響，采用標準k-ε湍流模型、無滑移壁面和標準壁面函數(shù)。計算過程中，壓力和速度采用PISO算法耦合，進行非穩(wěn)態(tài)計算。壓力采用Body Force Weighted方法進行差分，體積分數(shù)采用Ger-Reconstruct進行差分，其余變量采用一階迎風格式差分。

根據(jù)實驗情況給定邊界條件，入口邊界為氣體質(zhì)量入口，出口邊界為鋼包頂端氣體壓力出口，其他壁面為絕熱面。初始條件中設(shè)定合理的氣液兩相體積比例和初始速度。時間步長為非穩(wěn)態(tài)步長（采用非穩(wěn)態(tài)計算方式），迭代時間為10-5s，均方根殘差收斂標準為10-4。

將Gambit模塊生成的網(wǎng)格文件讀入Fluent中，設(shè)置其材料特性、初始條件與邊界條件。各部分的邊界條件設(shè)置如下：

1）底吹入口為inlet，采用流量入口，流量值的設(shè)定需根據(jù)流量進行計算，入口氣體溫度取環(huán)境溫度，298 K；

2）鋼包壁面設(shè)為wall，采用無滑移邊界條件，近壁面流場采用標準壁面函數(shù)計算；

3）鋼包上層空間出口為outlet，采用壓力出口，取值為環(huán)境壓力101 325 Pa，溫度為環(huán)境溫度1 300 K。

殘差收斂條件為除能量取10-6外，其他均為10-3?？紤]到氧槍為四面體網(wǎng)格，為得到正確計算結(jié)果，采用低殘差因子進行計算。模擬采用利用DELL M910，計算7 d可以計算約10 s左右的數(shù)據(jù)。

參考鋼包內(nèi)渣層、鋼液及底吹氣體的實際情況，得出物性參數(shù)（如表2所示）。

表2 鋼包內(nèi)物性參數(shù)表Tab.2 Physical parameters in ladle

2 結(jié)果分析

鋼包內(nèi)，不同的側(cè)吹流量對熔池攪拌效果不同，熔池在得到底吹氣的動能后進行傳輸運動（包括物質(zhì)，速度和溫度的傳輸運動），且熔池獲得的動能越大，傳輸運動越明顯。但在運動過程中產(chǎn)生形態(tài)不同的流股漩渦和相互干擾的熔池流股，傳輸運動受漩渦及流股的影響，使得不同區(qū)域內(nèi)熔池流股的速度和傳質(zhì)能力產(chǎn)生不同的改變，所以不同區(qū)域內(nèi)的速度和湍動能不一定呈線性分布。因此根據(jù)承鋼100 t鋼包進行1∶1模擬得到不同區(qū)域內(nèi)的速度和湍動能分布情況。

用底吹的鋼包，氣體從底部噴孔噴出，在鋼包中產(chǎn)生氣泡，氣泡在鋼液浮力作用下上升，同時帶動周圍鋼液向上運動，處于噴嘴正上方的液體形成強烈的向上股流。當氣液兩相區(qū)形成的上升流達到熔池液面后，氣體溢出熔池，而達到液面的鋼水被線面驅(qū)動流向四周，然后向下流動，從而形成循環(huán)流動。

同時，側(cè)吹氧槍產(chǎn)生的超音速氧氣射流穿過爐氣抵達渣層表面，因超音速氧氣沖力大于渣層表面張力，氧氣穿透渣層，并在熔池上方形成沖擊深坑。被熔池彈回的氧氣，部分直接進入爐氣，另一部分沿熔池表面運動，遇到轉(zhuǎn)爐爐壁阻礙向上流動形成循環(huán)流動。

2.1 側(cè)吹流量對熔池速度的影響

鋼液中速度死區(qū)集中在鋼包底吹流股之間的區(qū)域及熔池上部遠離脫磷氧槍一帶，如圖3所示。

圖3 鋼包不同截面處速度場分布Fig.3 Field distribution at different section rate of ladle

鋼液流速隨不同的側(cè)吹流量的增大而增大，熔池中速度死區(qū)逐漸降低。不同頂吹流量條件下，鋼包內(nèi)速度分別為0.016，0.017和0.018 m／s；速度死區(qū)分別占熔池體積的6.3%，7.6%和8.0%。不同流量下鋼包橫截面平均速度見表3。

表3 不同流量下鋼包橫截面平均速度Tab.3 Different flow cross-section average velocity of ladle

2.2 側(cè)吹流量對熔池湍動能的影響

湍動能越大，攪拌效果越好。圖4為鋼包不同的側(cè)吹流量時，不同截面處湍流動能云圖。從圖中可以看出，鋼包中下部離噴嘴越遠，湍動能越小。這是因為在底部吹氣量不大的情況下，上升氣流帶動的上升液體量相對較小。

圖4 鋼包不同截面處湍流動能云圖Fig.4 Contours of turbulent kinetic energy of ladle at different section

表4 不同流量下鋼包橫截面平均湍動能Tab.4 Average turbulent kinetic energy under different flow rate of ladle cross section

如表4所示，鋼液湍動能隨不同的側(cè)吹流量的增大而增大，同時氣體對熔池的攪拌能力也隨之提高。鋼包中湍動能分布由上到下，由底吹氣柱到四周降低。不同側(cè)吹條件下，鋼包內(nèi)湍動能分別為0.026，0.027和0.028 m2／s2。

2.3 側(cè)吹流量對熔池沖擊凹坑的影響

圖5顯示了不同側(cè)吹流量下的鋼包頂部鋼液沖擊面積，表5給出了不同頂吹流量下的熔池沖擊面積。

圖5 不同側(cè)吹流量下鋼包頂部鋼液沖擊面積Fig.5 Impact area at different side flow of molten steel ladle

表5 不同頂吹流量下熔池沖擊面積Tab.5 Impact area at different side flow of molten steel ladle

由圖5所示，不同側(cè)吹流量下鋼包頂處鋼液沖擊面積分別為4.68，4.72和5.03 m2。鋼液沖擊面積隨側(cè)吹氣量的增大而增大，因此加大側(cè)吹氣量對鋼包渣層有攪拌效果，促進鋼液的脫磷能力。

3 結(jié) 論

鋼包底吹在0.7R-120°位置上，底吹氣量為100 L／min時，經(jīng)模擬計算得出以下結(jié)論。

1）鋼液中速度死區(qū)集中在鋼包底吹流股之間的區(qū)域及熔池上部遠離脫磷氧槍一帶。鋼液流速隨不同的側(cè)吹流量的增大而增大，熔池中速度死區(qū)逐漸降低。不同頂吹流量條件下，鋼包內(nèi)速度分別為0.016，0.017和0.018 m／s；速度死區(qū)分別占熔池體積的6.3%，7.6%和8.0%。

2）鋼液湍動能隨不同的側(cè)吹流量的增大而增大，同時氣體對熔池的攪拌能力也隨之提高。鋼包中湍動能分布由上到下，由底吹氣柱到四周降低。不同側(cè)吹條件下，鋼包內(nèi)湍動能分別為0.026，0.027和0.028 m2／s2。

3）不同側(cè)吹流量下鋼包頂處鋼液沖擊面積分別為4.68，4.72和5.03 m2。鋼液沖擊面積隨側(cè)吹氣量的增大而增大，因此加大側(cè)吹氣量對鋼包渣層有攪拌效果，促進鋼液的脫磷能力。

4）本次模擬結(jié)果與水模實驗結(jié)果吻合，鋼液湍動能隨不同的側(cè)吹流量的增大而增大，最佳側(cè)吹流量為3 000 Nm3／h，速度死區(qū)集中在鋼包底吹流股之間的區(qū)域及遠離熔池上部氧槍一帶。

／References：

［1］ 馮 捷，張紅文.煉鋼基礎(chǔ)知識［M］.北京：冶金工業(yè)出版社，2005.FENG Jie，ZHANG Hongwen.The Steel Basic Knowledge［M］.Beijing：Metallurgical Industry Press，2005.

［2］ 雷 亞，楊治立，任正德，等.煉鋼學［M］.北京：冶金工業(yè)出版社，2010.LEI Ya，YANG Zhili，REN Zhengde，et al.Steelmaking Learn［M］.Beijing：Metallurgical Industry Press，2010.

［3］ 杜建良，王寶華，李敬賓，等.承鋼150噸轉(zhuǎn)爐低磷鋼冶煉工藝研究［J］.四川冶金，2011，33（5）：25-28.DU Jianliang，WANG Baohua，LI Jingbin，et al.The research of low phosphorus steel smelting technology of 150 tons converter in Chengde iron &steel group Co.Ltd.［J］.Sichuan Metallurgy，2011，33（5）：25-28.

［4］ 董香山，苑世林.承德地區(qū)資源狀況對承鋼持續(xù)發(fā)展的影響分析［J］.承鋼技術(shù)，2005（3）：1-2.DONG Xiangshan，YUAN Shilin.Effect of resources in Chengde area on the sustainable development of Chengde steel［J］.Chengde Iron and Steel Technology，2005（3）：1-2.

［5］ 曾建華，何 為，陳 永.轉(zhuǎn)爐半鋼冶煉脫磷技術(shù)研究［J］.鋼鐵釩鈦，2012，33（4）：68-72.ZENG Jianhua，HE Wei，CHEN Yong.Study on technologies of dephosphorization during semisteel smelting in converter［J］.Iron Steel Vanadium Titanium，2012，33（4）：68-72.

［6］ 潘秀蘭，王艷紅，梁慧智，等.世界鐵水預處理技術(shù)發(fā)展現(xiàn)狀［N］.中國冶金報，2011-02-01（C02）.PAN Xiulan，WANG Yanhong，LIANG Huizhi，et al.Current situation of the development of hot metal pretreatment technology all over the world［N］.China Metallurgical News，2011-02-01（C02）.

［7］ 劉 瀏，曾加慶.純凈鋼及其生產(chǎn)工藝的發(fā)展［J］.鋼鐵，2000，35（3）：68-72.LIU Liu，ZENG Jiaqing.Development of clean steel and its production process［J］.Iron &Steel，2000，35（3）：68-72.

［8］ 黃希祜.鋼鐵冶金原理［M］.北京：冶金工業(yè)出版社，1981.HUANG Xihu.Iron and Steel Metallurgy Principle［M］.Beijing：Metallurgical Industry Press，1981.

［9］ 郭漢杰.冶金物理化學教程［M］.北京：冶金工業(yè)出版社，2006.GUO Hanjie.Metallurgical Physical Chemistry Tutorial［M］.Beijing：Metallurgical Industry Press，2006.

［10］ 韓其勇.冶金過程動力學［M］.北京：冶金工業(yè)出版社，1983.HAN Qiyong.Metallurgical Process Dynamics［M］.Beijing：Metallurgical Industry Press，1983.

［11］ 張家蕓.冶金物理化學［M］.北京：冶金工業(yè)出版社，2004.ZHANG Jiayun.Metallurgical Physical Chemistry［M］.Beijing：Metallurgical Industry Press，2004.

［12］ 于 勇.FLUENT入門與進階教程［M］.北京：北京理工大學出版社，2008.YU Yong.FLUENT Introductory and Advanced Tutorial［M］.Beijing：Beijing Institute of Technology Press，2008.

［13］ 朱苗勇，井本健夫，肖澤強.多孔噴吹鋼包內(nèi)流動和混合過程的數(shù)學物理模擬［J］.金屬學報，1995，31（10）：435-439.ZHU Miaoyong，INOMOTO Takeo，XIAO Zeqiang.Mathematical and physical modeling of flow and mixing in ladles stirred by argon through multi-tuyere［J］.Acta Metallurgica Sinica，1995，31（10）：435-439.

［14］ 沈巧珍，吳曌環(huán)，陽 方，等.100 t鋼包底吹氬數(shù)值的模擬研究［J］.武漢科技大學學報，2010，33（6）：566-569.SHEN Qiaozhen，WU Zhaohuan，YANG Fang，et al.Mathematical simulation of 100 t bottom argon blowing ladle［J］.Journal of Wuhan University of Science and Technology，2010，33（6）：566-569.

［15］ 幸 偉，沈巧珍，王曉紅，等.鋼包底吹氬過程數(shù)學物理模擬研究［J］.煉鋼，2005，21（6）：33-36.XING Wei，SHEN Qiaozhen，WANG Xiaohong，et al.Water model and numerical simulation on argon blowing from bottom of ladle［J］.Steelmaking，2005，21（6）：33-36.

［16］ 任三兵，陳義勝，黃宗澤，等.大型鋼包雙孔吹氬最佳位置的探討［J］.包頭鋼鐵學院學報，2003，22（3）：193-197.REN Sanbing，CHEN Yisheng，HUANG Zongze，et al.The discussion of optimized location of two holes in large ladle furnace during argon blowing［J］.Journal of Baotou University of Iron and Steel Technology，2003，22（3）：193-197.