殷 攀，趙洪亮, 2 ，張立峰,2, 王 森,，張建坤，范 巍

(1. 北京科技大學(xué)冶金與生態(tài)工程學(xué)院，北京100083；2. 稀貴金屬綠色回收與提取北京市重點(diǎn)實(shí)驗(yàn)室，北京100083； 3. 謙比希銅銅冶煉公司北京代表處，北京100029)

艾薩旋流頂吹熔煉過(guò)程數(shù)值模擬

殷 攀1，趙洪亮1, 2，張立峰1,2, 王 森1,3，張建坤3，范 巍3

(1. 北京科技大學(xué)冶金與生態(tài)工程學(xué)院，北京100083；2. 稀貴金屬綠色回收與提取北京市重點(diǎn)實(shí)驗(yàn)室，北京100083； 3. 謙比希銅銅冶煉公司北京代表處，北京100029)

以謙比希銅冶煉廠ISA熔煉爐為原型，通過(guò)水模擬和時(shí)下應(yīng)用廣泛的數(shù)值模擬方法對(duì)ISA爐頂吹熔煉過(guò)程進(jìn)行了初步研究.以噴槍為研究中心，探討了噴槍中是否加入旋流片對(duì)熔池流場(chǎng)以及熔煉過(guò)程中噴濺量的影響.結(jié)果表明：噴槍加入旋流片后有利于氣流沿徑向分布，形成鼓狀氣泡，氣流對(duì)液相的沖擊深度會(huì)減小，從而減輕對(duì)渣锍界面的干擾，有利于渣、锍分離，此外還能夠明顯減小爐內(nèi)的渣的噴濺量；不過(guò)噴槍加入旋流片后，對(duì)液相區(qū)域的攪拌能力會(huì)減弱.

ISA熔煉；旋流片；水模型；流場(chǎng)；數(shù)值模擬

相較于其它的熔池熔煉方法，艾薩(ISA)熔煉法具有熔煉速度快，生產(chǎn)率高；原料適應(yīng)性強(qiáng)；爐子結(jié)構(gòu)緊湊；工藝流程和操作相對(duì)簡(jiǎn)單等優(yōu)點(diǎn)[1]，在有色金屬冶煉工業(yè)逐漸得到人們的認(rèn)可.2002年云南銅業(yè)引進(jìn)國(guó)內(nèi)第一臺(tái)艾薩熔煉爐，2009年我國(guó)于贊比亞銅帶省建成海外投資的最大銅冶煉項(xiàng)目，其使用的也是艾薩熔煉法[2].因此近年來(lái)國(guó)內(nèi)關(guān)于艾薩爐的研究工作也逐漸開(kāi)展起來(lái)，其中浸沒(méi)式頂吹噴槍一直是學(xué)者研究的重點(diǎn).謙比希銅冶煉廠為了進(jìn)一步提高艾薩銅冶煉的生產(chǎn)能力，加大冶煉強(qiáng)度，需要對(duì)原有艾薩熔煉爐進(jìn)行擴(kuò)建，若單純依靠在工業(yè)生產(chǎn)中改變條件摸索出最佳的生產(chǎn)工藝將會(huì)花費(fèi)巨大的時(shí)間和經(jīng)濟(jì)成本.由于艾薩熔煉是在頂吹空氣的作用下對(duì)熔池內(nèi)熔體和爐料進(jìn)行劇烈的攪拌過(guò)程，涉及到熱量、質(zhì)量的傳遞和化學(xué)反應(yīng)速率等現(xiàn)象[3-4]，因此本文試圖用模擬的方法初步探討當(dāng)外部操作條件變化時(shí)會(huì)對(duì)熔煉過(guò)程中熔池流場(chǎng)產(chǎn)生什么影響，以期找到艾薩熔煉法的最佳工藝條件并為實(shí)際工業(yè)生產(chǎn)提供指導(dǎo).

1 研究方法

1.1 實(shí)驗(yàn)方法

本文采用了水模擬和數(shù)值模擬兩種方法對(duì)銅冶金ISA熔煉過(guò)程進(jìn)行研究，兩種方法各有自己的優(yōu)缺點(diǎn)，水模擬實(shí)驗(yàn)平臺(tái)容易搭建，研究?jī)?nèi)容比較直觀，但對(duì)于實(shí)際熔煉過(guò)程的復(fù)雜性卻很難模擬出熔煉過(guò)程中物質(zhì)傳輸、熱量傳輸以及化學(xué)反應(yīng)等現(xiàn)象并得出相關(guān)數(shù)據(jù)，只能定性分析.數(shù)值模擬則彌補(bǔ)了這方面的缺陷，拋開(kāi)模型計(jì)算的準(zhǔn)確性，數(shù)值模擬開(kāi)放的用戶自定義設(shè)置能夠很大程度地滿足各種研究問(wèn)題的需要，并且對(duì)于相關(guān)數(shù)據(jù)的獲取更方便.不過(guò)數(shù)值模擬并不能完全取代水模擬，兩者之間可以相互驗(yàn)證，更準(zhǔn)確的反映實(shí)際問(wèn)題[5-6].

水模擬實(shí)驗(yàn)以謙比希銅冶煉廠ISA爐為原型，依據(jù)相似原理按原型(Lp)與模型(Lm)比為10∶1 比例制作物理模型.實(shí)驗(yàn)時(shí)，用水代替冰銅，硅油代替爐渣，在此基礎(chǔ)上根據(jù)工業(yè)尺寸得到實(shí)驗(yàn)?zāi)Ｐ鸵合喔叨葹?8.5 cm，其中水相高度 8.5 cm，硅油高度10 cm.熔池流動(dòng)主要由氣泡浮力引起，為滿足動(dòng)力相似必須要求模型與原型的修正弗魯?shù)聹?zhǔn)數(shù)Fr(流動(dòng)由慣性力和重力主導(dǎo))相等，由此得到流量比：

其中，m代表模型，P代表原型，ρl,m為水的密度、ρl,p為冰銅的密度、ρg,m=ρg,p為空氣的密度.

工業(yè)通氣量(標(biāo)準(zhǔn)態(tài)下)為30 000～32 000 m3/h，由此得到實(shí)驗(yàn)通氣量為44.7～47.68 m3/h.具體實(shí)驗(yàn)時(shí)選取通氣量為45 m3/h.實(shí)驗(yàn)借助的測(cè)試工具有高速攝像機(jī)，每秒能拍攝500張圖片，用于捕捉熔池流場(chǎng)的變化.實(shí)驗(yàn)用到的通氣設(shè)備為廣渦渦漩空壓機(jī)，工作時(shí)通氣量和壓力能達(dá)到：2 m3/min，0.8 MPa.實(shí)驗(yàn)時(shí)頂吹浸沒(méi)式噴槍的浸沒(méi)深度可以控制在3、 4和 5 cm，噴槍直徑分別為3.5和4 cm.根據(jù)需要，實(shí)驗(yàn)側(cè)重設(shè)計(jì)了旋流片以考察旋流噴吹對(duì)流場(chǎng)的影響[7].

圖1為實(shí)驗(yàn)裝置簡(jiǎn)圖，ISA爐模型置于一矩形槽中，矩形水槽中充滿水，保持液位高度高于ISA模型中液位，由此可以消除高速相機(jī)拍攝時(shí)爐體曲面的折射，使熔池區(qū)域光線明暗一致.

圖1 實(shí)驗(yàn)裝置簡(jiǎn)圖Fig.1 Schematic view of the experimental device

圖2所示為實(shí)驗(yàn)所用的頂吹噴槍模型，為保證噴槍中氣體均勻分布，在噴槍頂部置有三個(gè)進(jìn)氣口，在噴槍內(nèi)部靠近下端部安裝一旋流片.調(diào)節(jié)噴槍插入深度為3 cm，調(diào)節(jié)流量計(jì)氣量為45 m3/h，實(shí)驗(yàn)通氣的同時(shí)高速相機(jī)開(kāi)始記錄爐內(nèi)液相的運(yùn)動(dòng)狀況.

圖2 帶旋流片噴槍模型Fig.2 Model of the gunjet with swirl plate

1.2 模擬方法

數(shù)值模擬方法也是基于艾薩熔煉爐，利用常用的建模軟件Gambit建立等比例的簡(jiǎn)易幾何模型，包括噴槍模型和艾薩爐體模型.噴槍模型是為了計(jì)算得到艾薩爐的速度入口信息，模擬過(guò)程主要在艾薩爐體中完成.噴槍模型和爐體幾何模型尺寸如圖3、圖4.圓柱形爐體的直徑為3.84 m，高度為7.173 m，球形底部的半徑為5.1 m，中心噴槍的直徑為0.4 m，噴槍距離圓柱爐體的底部為1.55 m.爐體中，冰銅層的高度為0.5 m，混合層的高度為0.85 m，渣層高度為1.85 m.

模型內(nèi)流場(chǎng)的數(shù)值計(jì)算使用的是時(shí)下流行的商用CFD軟件包Fluent.噴槍模型為單相流模型，采用standardk-ε模型計(jì)算，根據(jù)進(jìn)氣量(標(biāo)準(zhǔn)態(tài)下)Q=30 000 m3/h得到速度入口的進(jìn)氣速度大小，出口邊界設(shè)為壓力出口.求解器設(shè)置為穩(wěn)態(tài)求解，設(shè)置時(shí)間步長(zhǎng)為0.001 s，步數(shù)為2 000，計(jì)算得到出口信息.

圖3 ISA爐噴槍模型和旋流片 Fig.3 The model of ISA gunjet and swirl plate

圖4 ISA爐爐體模型(d=0.4 m，h=0.4 m) Fig.4 Model of ISA furnace (d=0.4 m, h=0.4 m)

艾薩爐體模型為氣渣锍的三相流動(dòng)過(guò)程，涉及到各項(xiàng)的動(dòng)量和能量交換，該模型計(jì)算包含如下特點(diǎn)：計(jì)算模型采用VOF多相流模型和realizablek-ε湍流模型；速度與壓力的耦合采用PISO算法，并對(duì)動(dòng)量、湍動(dòng)能和湍流耗散率選用二階迎風(fēng)格式離散；入口邊界條件為速度入口；出口邊界條件為壓力出口[8-9].關(guān)于該模型的控制方程如下：

(1)連續(xù)性方程：

(其中αq為相q(除主相)的體積分?jǐn)?shù)，ρq為相q的密度)

(2)動(dòng)量方程：

(3)realizablek-ε湍流方程:

爐體模型采用六面體網(wǎng)格進(jìn)行劃分，數(shù)目在28萬(wàn)左右，并對(duì)噴槍出口處速度變化較大的區(qū)域進(jìn)行了加密處理，求解器設(shè)置為瞬態(tài)求解，設(shè)置求解步長(zhǎng)為0.001s，總共計(jì)算時(shí)間為20s.由于湍流過(guò)程的相關(guān)變量存在湍流脈動(dòng)，所以對(duì)2s后的計(jì)算結(jié)果統(tǒng)計(jì)平均后再進(jìn)行計(jì)算.關(guān)于艾薩爐內(nèi)各項(xiàng)的物性如表1所示.

表1 冰銅，渣的實(shí)驗(yàn)參數(shù)

2 實(shí)驗(yàn)結(jié)果與分析

2.1 旋流片對(duì)噴槍速度分布的影響

圖5所示為噴槍直徑為0.4 m，不加入旋流片時(shí)出口處的速度矢量分布圖，圖6所示為噴槍直徑為0.4 m，噴槍中加入旋流片時(shí)出口處的速度矢量分布圖，兩種條件下保持入口通氣量均為 30 000 m3/h.很明顯，不加旋流片噴槍出口處速度分布中間略低，沿半徑方向增加后減小，速度變化梯度小，整體在55～88 m/s之間，各處速度方向一致，均是垂直于出口界面向外；加入旋流片的噴槍出口處速度分布也是中間較小，沿徑向先增大后減小，但速度梯度變化較大，較大的接近110 m/s，中心處速度甚至反向.各處速度方向也有一定夾角，總體呈現(xiàn)出螺旋的方向.由于在徑向有分速度，其有利于氣泡在徑向延展.

圖5 噴槍出口處速度分布(d=0.4 m，無(wú)旋流片)Fig.5 Velocity vector on the gunjet tip (d=0.4 m, without swirl plate)

圖6 噴槍出口處速度分布(d=0.4 m，含旋流片)Fig.6 Velocity vector on the gunjet tip (d=0.4 m, with swirl plate)

圖7 氣流沖擊最深時(shí)氣泡的形狀d=4 cm, h=5 cm, Q=45 m3/h)Fig.7 Shape of the bubble with deepest gas penetration (a)—不帶旋流片噴槍； (b)—帶旋流片噴槍

圖7所示為氣體噴吹達(dá)到最深處時(shí)刻的圖片.實(shí)驗(yàn)條件為：噴槍直徑為4 cm，噴槍插入深度為5 cm，頂吹噴槍通氣量為Q=45 m3/h，其中(a)所示為噴槍中不含旋流片，(b)所示噴槍中加入旋流片.

從這組對(duì)比實(shí)驗(yàn)可以看出，噴槍加入旋流片后形成的氣泡圓滑，幾何外形好，氣泡整體呈現(xiàn)鼓狀并且氣泡底部中間往里凹陷，有利于氣泡沿徑向擴(kuò)展分布；噴槍不加旋流片形成的氣泡表面不規(guī)則有很多凸起并且氣泡底部往外凸，相比之下有更大的透入深度.除此之外噴槍加入旋流片后，在相同氣量下，氣流對(duì)底部銅渣界面的擾動(dòng)較小，這有利于冰銅與渣的分離與澄清.

圖8 艾薩爐內(nèi)監(jiān)控點(diǎn)位置Fig.8 Monitoring points in ISA furnace

2.2 旋流片對(duì)熔池流場(chǎng)和噴濺量的影響

此外對(duì)熔池中氣液界面、渣锍界面和底部這三個(gè)位置附近的點(diǎn)進(jìn)行了速度監(jiān)控，圖8所示為監(jiān)控點(diǎn)的位置.圖10(a) 、(b) 和(c)所示為有無(wú)旋流片條件下各監(jiān)控點(diǎn)處平均速度的對(duì)比圖.

圖9 噴濺量隨時(shí)間變化曲線Fig.9 Splashing with time

圖10 各監(jiān)控點(diǎn)平均速度隨距離的變化(d=0.4 m， h=0.3 m)Fig.10 Relationship between mean velocity and distance from the axis at different monitoring points(d=0.4 m, h=0.3 m) (a)—z=1.8 m; (b)—z=0.75 m; (c)—z=0 m)

從CFD的后處理中可以很直接的得到爐體流場(chǎng)速度變化云圖，不過(guò)對(duì)于某些位置數(shù)據(jù)的微小變化從云圖中并不能很好地分辨，所以企圖建立監(jiān)控點(diǎn)，將云圖的信息轉(zhuǎn)化為直接的數(shù)據(jù).從 ISA爐體簡(jiǎn)易模型可以看出整個(gè)爐體是呈圓柱狀的，計(jì)算區(qū)域是對(duì)稱的，只要清楚爐體1/4區(qū)域的流場(chǎng)信息，整個(gè)流場(chǎng)的信息都可以知道.所以在選取監(jiān)控點(diǎn)的時(shí)候只選取噴槍某一側(cè)的位置，其他對(duì)稱區(qū)域的信息也是一樣的.

在實(shí)際生產(chǎn)過(guò)程中，ISA爐下料口容易黏結(jié)，這極有可能是由于爐內(nèi)渣的噴濺過(guò)于嚴(yán)重造成的，因此運(yùn)用Fluent自帶的UDF編程功能對(duì)頂吹過(guò)程中爐內(nèi)渣的噴濺量進(jìn)行了考察.實(shí)驗(yàn)時(shí)，計(jì)算了噴濺高度超過(guò)液面1.25 m的渣量，渣量用渣密度、渣體積分?jǐn)?shù)和網(wǎng)格體積這三者的乘積表示.圖9中曲線1為噴槍中不加旋流片的情況，曲線2為加入旋流片的情況，可以明顯的看到噴槍加入旋流片后噴濺量明顯減小，一方面這有可能是加入旋流片后，液相區(qū)域徑向速度增大，在豎直方向上速度減??；而不加旋流片，液相在豎直方向上速度波動(dòng)很大，因此噴濺量也較多.

圖10(a)、(b)和(c)所示為噴槍直徑為d=0.4 m，噴槍浸沒(méi)深度h=0.3 m時(shí)，在10 s到20 s這段時(shí)間內(nèi)各監(jiān)控點(diǎn)處速度變化的平均值，圖10(a)所示為1、2、3、4四個(gè)監(jiān)控點(diǎn)，z=1.8 m；圖10(b)所示為5、6、7、8、9五個(gè)監(jiān)控點(diǎn)，z=0.75 m；圖10(c)所示為10、11、12、13、14五個(gè)監(jiān)控點(diǎn)，z=0 m.從三個(gè)圖中都可以看出，噴槍加入旋流片后各監(jiān)控點(diǎn)的平均速度都比不加旋流片的時(shí)候要低，也就是說(shuō)加入旋流片后液相區(qū)域的攪拌減弱了.液相區(qū)域越靠近爐底的部分，和側(cè)邊壁位置速度會(huì)大幅減小.

3 結(jié) 論

(1)在ISA爐頂吹噴槍中加入旋流片有利于形成鼓狀氣泡，擴(kuò)展氣流在徑向的分布，并且氣流對(duì)熔池的沖擊深度減小，對(duì)渣、锍界面擾動(dòng)減小，有利于渣、锍分離.

(2)噴槍加入旋流片后對(duì)于熔池內(nèi)液相區(qū)域的攪拌是減弱的.

(3)噴槍加入旋流片后可以明顯減小爐內(nèi)的噴濺量.

[1]朱祖澤, 賀家齊. 現(xiàn)代銅冶金學(xué)[M]. 北京:科學(xué)出版社, 2003: 273-291. (Zhu Z Z, He J Q. Modern copper metallurgy[M]. Beijing: Science Press, 2003: 273-291.)

[2]黃善富, 楊新國(guó). 贊比亞謙比希銅冶煉艾薩爐技術(shù)集成化創(chuàng)新項(xiàng)目與成果運(yùn)用[J]. 中國(guó)有色金屬, 2010(S1): 242-251. (Huang S F, Yang X G. Integrated innovation project and application of ISA copper smelting technology in chambishi zambia[J]. China Nonferrous Metals, 2010(S1): 242-251.)

[3]張波, 史誼峰, 楊曉琴, 等. 艾薩爐水模擬氣量與混合均勻時(shí)間研究[J]. 云南冶金, 2007, 36(4): 31-33. (Zhang P, Shi Y F, Yang X Q,etal. Study on hydraulic modeling of the airflow rate and uniform mixing time for ISA smelting process[J]. Yunnan Metallurgy, 2007, 36(4): 31-33. )

[4]張波, 洪新, 陳朝軼, 等. 艾薩爐水模擬研究[C]//全國(guó)冶金物理化學(xué)學(xué)術(shù)會(huì)議. 2008. (Zhang P, Hong X, Chen C Y,etal. ISA smelting process hydraulic modeling research[C]//National Conference on Metallurgical Physical Chemistry. 2008.)

[5]韓麗輝, 劉云, 朱榮. CFX數(shù)值模擬在冶金工程專業(yè)中的應(yīng)用[J]. 實(shí)驗(yàn)室研究與探索, 2008, 27(10): 52-56. (Han L H, Liu Y, Zhu R. Application of CFX numerical simulation in metallurgical engineering specialty[J]. Research and Exploration in Laboratory, 2008, 27(10): 52-56.)

[6]張大江, 陳登福, 王翠娜, 等. 不同寬度板坯結(jié)晶器內(nèi)流場(chǎng)的水模型和數(shù)值模擬研究[J]. 特殊鋼, 2010, 31(1): 14-17. (Zhang D J, Chen D F, Wang C N,etal. A study on water model and numerical simulation of flow field in slab casting mold with different width[J]. Special steel, 2010, 31(1): 14-17).

[7]Huda N, Naser J, Brooks G,etal. CFD modeling of swirl and nonswirl gas injections into liquid baths using top submerged lances[J]. Metallurgical and Materials Transactions B, 2010, 41(1): 35-50.

[8]Huda N, Naser J, Brooks G,etal. CFD modeling of gas injection in top submerged lance smelting[C]//The Minerals, Metals & Materials Society. 2009.

[9]Chibwe D K, Akdogan G, Taskinen P,etal. Modelling of fluid flow phenomena in peirce-smith copper converters and analysis of combined blowing concept[C]//The Southern African Institute of Mining and Metallurgy: Pyrometallurgical Modelling. 2014.

Numerical simulation of ISA smelting process with top swirl blowing

Yin Pan1, Zhao Hongliang1, 2, Zhang Lifeng1, 2, Wang Sen1,3, Zhang Jiankun3, Fan Wei3

(1. School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, China; 2. Beijing Key Lab of Rare & Precious Metals Green Recycling and Extraction, Beijing 100083, China; 3. Chambishi Copper Smelter LTD, Beijing Representative Office, Beijing 100029, China)

Taken ISA smelting furnace in Chambish Copper Smelter Plant as a prototype,through water simulation and numerical simulation used widely nowadays,the ISA smelting process was studied preliminarity.Taken the gunjet as a focus,the effect of swirl plate in the gunjet on bath flow field and splashing in melting process was discussed. The result showed that the top submerged gunjet with swirl plate is beneficial for radical distribution of gas flow and to form drum-shaped bubbles. These the impact depth of the gas to the bath decreases,so it is benefical for the matte-slag separation.However,the gunjet with swirl plate can not give a strong stirring to the bath.

ISA smelting; swirler; water model; flow field;numerical simulation

10.14186/j.cnki.1671-6620.2017.02.005

TF 806

A

1671-6620(2017)02-0104-06