電渣重熔工藝中渣金兩相的流動、傳熱及凝固

李寶寬，陸秋敏，陳 瑞，宋照偉

(1.東北大學 材料與冶金學院，沈陽 110004;2.沈陽鑄造研究所，沈陽 110021))

電渣重熔工藝中渣金兩相的流動、傳熱及凝固

李寶寬1，陸秋敏1，陳 瑞2，宋照偉2

(1.東北大學 材料與冶金學院，沈陽 110004;2.沈陽鑄造研究所，沈陽 110021))

本文利用計算流體力學方法對電渣重熔過程中的渣相和金屬相流動、傳熱以及相變進行了耦合數值模擬計算.通過分析結果，得到了電渣重熔過程中的溫度場分布、流場分布以及其中的凝固情況.并在此基礎之上，通過改變熔速來查看其對電渣重熔結果的影響，發現金屬熔速增大，會使金屬熔池內的溫度普遍提升，金屬熔池液變得更深更寬.

電渣重熔;溫度場;凝固場;流場

電渣重熔工藝越來越廣泛地應用于機械、交通、電力和化工等行業，用于制備模具鋼、工具鋼和機車曲軸、發電機組中轉子、定子、導葉等關鍵鑄件.電渣重熔過程是一種消耗金屬電極再重熔的過程.通過控制重熔過程中的精煉和凝固條件，制備出比原材料更高質量的鋼錠.電渣重熔的基本過程為:在電極和水冷壁間通過電流，由于電渣的電阻率較高產生大量焦耳熱，使得渣產生高溫，電極不斷地被熔化，熔化的液滴穿過渣層進入金屬熔池，從而凝固成擁有更好的組織結構和性能的新的鋼錠.電渣重熔過程包含電磁場、渣金兩相流動、傳熱和相變等多種復雜的物理過程.掌握這些復雜的物理過程對改善電渣冶金工藝、提高電渣鋼錠質量具有重要的意義.

關于鋼錠的電渣重熔工藝中傳輸與凝固方面，國內外已有較豐富的研究結果.早在1980年M.Choudhary 和 J.Szekely［1］就對小型鋼錠的電渣重熔結晶器內傳熱、凝固和渣/鋼熔體流動有了初步的研究，其方法和結果至今仍有借鑒作用.近年來，Viswanathan 等人［2，3］分析了單電極工藝參數對結晶器內溫度分布、流體流動及液池形狀影響.Kompan 等人［4］和 Bandyopadhyay 等人［5］研究了電渣熔煉高強度合金鋼的工藝特點及操作問題.Krane等人［6］分析了溶質元素偏析發展規律及相關參數影響.Mitchell等人［7，8］詳細研究了電渣重熔鎳基合金鋼錠內柱狀晶向等軸晶轉變與過渡行為.Kellar等人［9］和 Cefalu 等人［10］還發展了電渣重熔過程數學模型，給出了溫度場和液池形狀等工藝參量的詳細信息.

劉福斌等人［11］對電渣重熔過程渣池流場進行了數值模擬，考慮了電磁力和結晶器壁冷卻導致的浮力，并考察了重熔電流、填充比等一些參數對渣池流動的影響.

由于前人的模擬計算大多為單相的模擬，只對鋼錠進行計算沒有考慮渣的影響，僅把金屬鋼錠和電渣的關系簡化為一個渣溫的影響;或者只考慮液態鋼液的流動，而沒有對凝固的計算.本文在前人的基礎上同時考慮液態電渣和金屬鋼錠的影響，對電渣重熔過程進行了傳熱、流動和凝固三個方面的分析，并在此基礎上分析了電渣重熔過程中的操作參數對重熔本身的影響.

1 物理模型

1.1 幾何模型

由于實際的電渣重熔大多是用于鑄造圓形鋼錠，它的幾何形狀是軸對稱的，并且在豎直方向上的工藝條件是穩定的.所以計算區域可以簡化為軸對稱圖形(圖1).

圖1 鋼錠幾何模型Fig.1 Geometric simulation of steel ingot

結晶器的具體尺寸為:半徑0.07 m，高度為0.3 m.

1.2 物性參數

鑄坯在凝固過程中，鋼液的密度、比熱、黏性、導熱系數等物性參數會隨著溫度的變化而變化，但是這些變化對本模型的計算結果影響較小，故認為鋼液的物性參數為常數，如表1所示.

表1 鋼種參數Table 1 Parameters of steel

渣的物性參數也設為常數，如表2所示.

表2 電渣參數Table 2 Parameters of slag

1.3 基本假設

實際的電渣重熔過程具有復雜的邊界條件和工藝條件，幾乎不太可能對其進行精確的解析求解.在對電渣重熔進行計算之前，通常都必須針對具體的研究對象作如下假設:

(1)電渣重熔過程穩定，可將其簡化為二維計算;

(2)鋼液為不可壓縮牛頓流體，物性參數為常數;

(3)忽略因冷卻而引起的鋼錠收縮;

(4)電渣重熔過程的電磁現象相對于凝固和流動影響較小，因此不予考慮;

(5)將由電流產生的熱量簡化為有效渣溫.

1.4 流場控制方程

本文中流動為穩態，不可壓流體的連續方程、流體的動量方程、傳熱方程、描述湍流的κ－ε方程與湍動黏度μt方程如下:

在式(1)～(5)中:μeff為湍流有效粘度，σT為湍流普朗特數，Gk是由于平均速度梯度引起的湍動能k的產生項，C1ε和 C2ε為經驗常數，σk和 σε分別是與湍動能k和耗散率ε對應的Schmit數.

1.5 邊界條件

由于電渣重熔熱量主要來自于渣的焦耳熱，渣為發熱體且保持恒定溫度，故設置為第一類邊界條件，設渣溫為1 900 K，由渣溫關系換算熔速為0.57 kg/min.

鋼錠的熱量要通過渣皮、空隙、結晶器、冷卻水、以及外壁傳遞到系統之外，過程較為繁瑣，實際模擬較為復雜，為方便模擬，也采用第一類邊界條件，設壁面冷卻溫度為900 K.在對稱軸處為絕熱面，熱流密度為零.

2 結果分析

2.1 電渣重熔溫度場分布

圖2 電渣重熔中溫度等值線分布Fig.2 The isoline distribution of temperature in ESR system

圖2為整個鋼錠的數值模擬溫度分布場，從中可以看出，溫度由渣、金屬液池到凝固鋼錠逐漸降低，在凝固鋼錠區溫度自上而下均勻下降，在底部溫度基本一致.熔池中部溫度較為一致，邊緣部位溫度梯度較大.在鋼錠下部溫度基本一致.渣區溫度基本為1 600 K以上，液態金屬區約為1 300～1 600 K，凝固金屬區為1 300 K以下.布，即電渣重熔過程中渣金界面(鋼錠高度為0.245 m)、熔池內部(鋼錠高為0.21 m)以及鋼錠內部(鋼錠高為0.17 m)由鋼錠中心到邊緣的溫度分布圖.可見渣金界面和熔池內部在鋼錠中間與邊緣的溫差較大，這是因為結晶器受到內部冷卻水強制冷卻造成的.同時兩者在中間的溫度相差不大，在邊緣熔池內部的溫度明顯低于渣金界面，這是與熔池區域邊界受到冷卻形成凝固層溫度必然降低是相符的，在凝固鋼錠內部溫度分布較為平緩.

圖3 不同高度的溫度分布Fig.3 The temperature distribution of different heights

2.2 電渣重熔凝固場

圖4為電渣重熔凝固圖，在鋼錠中部形成一個碗狀的固液兩相的交界區，上部為液相包括液態電渣和未凝固的液態金屬，下部為凝固的金屬，交界區內為固液兩相混合區.可以清晰地看出熔池的形狀，最上面為金屬的液相線，下面為固相線，中間為兩相混合區.在底部出現的凸起應該是由于熔池內的回流造成的.

圖5為渣金界面、熔池內部和鋼錠內部的由中心到邊緣的液相率圖，由圖可見渣金界面基本為液相，僅在邊緣處有稍微波動.而在熔池內部，中間部位仍為液相，在邊緣液相率則出現明顯下降并且到達零，說明在邊緣已經形成凝固殼.在鋼錠內部液相率則完全是零，說明已經凝固，與前面的場圖也是相符的.

2.3 電渣重熔速度場

圖6分別為電渣重熔的速度矢量圖和速度等值線圖，可見速度區域分布在渣相區和熔池內，在鋼錠凝固區速度為零.在渣相中部速度較為集中且方向向下，其他區域則較為稀疏，在中部左右兩側形成兩個較小的回流.在金屬熔池內也存在兩個較為明顯的回流，而熔池邊緣也存在微小的回流.速度最大處大約在渣相的中部，在整個流場區域有漩流處相比其他區域速度較大.

圖4 電渣重熔中液相率等值線分布Fig.4 The isoline distribution of liquid fraction in ESR system

圖5 不同高度的液相率分布Fig.5 The liquid fraction distribution of different heights

圖7為渣金界面、熔池內部和凝固鋼錠內部的速度曲線圖.可見，渣相的中部的速度值最大，然后越往邊緣速度越小并在邊緣處有微小波動.在金屬熔池內，速度變化較為頻繁，出現兩個較大的峰值，這與熔池內形成的兩個回流是相符的.而在凝固鋼錠區速度則為零.

2.4 熔速對電渣重熔的影響

圖6 電渣重熔中速度分布Fig.6 The distribution of velocity in ESR system

圖7 不同高度的速度大小分布Fig.7 The velocity magnitude distribution of different heights

為了考察熔速對電渣重熔的影響，以熔速0.57 kg/min為標準.為使得結果具有可比性，保持冷卻條件等其他條件不變，改變熔速，得到一個在0.35 kg/min和0.79 kg/min的電渣重熔情況，然后三者進行對比.

為了更準確地說明3種熔速情況下溫度和凝固的變化情況，下面從縱向和橫向兩個方向描述其變化，縱向選擇在金屬相區的中軸線上，橫向仍為金屬熔池內部，即是鋼錠高度為0.21 m處.

圖8是在鋼錠中軸線上的溫度分布.由圖可知在鋼錠底部三者的溫度基本相同，沿軸向越靠近熔池溫度差別越大，且變化較為均勻，在達到0.19 m時溫度變化較為平緩.熔速最大的熔池溫度約為1 800 K，中速為1 550 K，熔速最低的為1 400 K.

圖8 不同熔速對鋼錠中心線溫度的影響Fig.8 Effect of different melting speeds on temperature of ingot centerline

圖9 不同熔速對鋼錠中心線液相率的影響Fig.9 Effect of different melting speeds on liquid fraction of ingot centerline

圖9為不同熔速下在金屬相區沿中軸線的液相率變化，可見熔速最大的0.79 kg/min最先由固相變為液相，并看出熔速為0.79、0.57和0.35 kg/min的熔池深度分別約為9、8和7.5 cm.

圖10為在3種熔速下金屬熔池內的橫向溫度分布圖.隨著熔速的降低即渣溫的降低，內部溫度也逐漸降低，三者的溫度變化趨勢較為一致，在熔池內溫度基本保持一致，在邊緣處迅速降低.

圖11為金屬熔池內的液相率變化曲線.金屬熔池內液相率均為1，靠近邊緣處發生變化.由液態最先變化的是熔速最低的0.35 kg/min，最晚變化的為0.79 kg/min.可知熔速越大金屬熔池也越寬.

由以上的分析可總結出:當熔速增加，渣溫升高時，會使得熔池內部的溫度也相應升高，并且熔池的寬度和深度都會有所增加.熔速的增大或減小改變對熔池形狀和溫度造成的變化較為均勻.

圖10 不同熔速對鋼錠高0.21 m溫度影響Fig.10 Effect of different melting speeds on temperature of ingot height at 0.21 m

圖11 不同熔速對鋼錠高度0.21 m處液相率影響Fig.11 Effect of different melting speeds on liquid fraction of ingot height at 0.21 m

3 結論

(1)通過計算得到了電渣重熔內部溫度、流動和相變的情況.從上到下溫度逐漸下降，從中心向邊緣溫度也逐漸降低;在渣相的中部形成一個集中的向下的流動，且速度較大，在金屬熔池內，左右兩邊有兩個金屬液的回流;在金屬熔池底部形成一個碗狀的凝固邊界.

(2)通過3種不同金屬熔速的計算比較，可以看出金屬熔速的改變可以影響電渣重熔的溫度分布和凝固，金屬熔速增大，使得金屬熔池內的溫度普遍提升，金屬熔池也變得更深更寬，反之亦然.

［1］Choudhary M，Szekely J.Modeling of fluid flow and heat transfer in industrial－ scale ESR system［J］.Ironmaking and Steelmaking，1981(5):225－231.

［2］Viswanathan S，Melgaard D K，Patel D K，et al.Effect of processing parameters on temperature profiles，fluid flow，and pool shape in the ESR process［C］//LMPC 2005 －Proceedings of the 2005 International Symposium on Liquid Metal Processing and Casting，2005:145－154.

［3］Viswanathan S，Patel A D，Melgaard D K，et al.2D steadystate analysis of temperature profiles，fluid flow，and pool shape in the ESR process，Modeling of Casting［J］.Welding and Advanced Solidification Processes － XI，2006(2):977－984.

［4］Kompan Ya Yu，Zhadkevich M L，Protokovilov I V，et al.High strength alloys by magnetically controlled electroslag melting［J］.Materials Technology，2003，18(4):214 － 217.

［5］Bandyopadhyay，RaoT R，Krishna P.Development of ultrahigh strength steel through electroslag remelting with inoculation［J］.Ironmaking ＆ Steelmaking，2006，33(4):337－343.

［6］Krane，Matthew J M，Cefalu Shawn，et al.Segregation development in electroslag remelting processes，Modeling of Casting［J］.Welding and Advanced Solidification Processes －XI，2006:969 －976.

［7］Mitchell A.Solidification in remelting processes［J］.Materials science and engineering A ，2005，413:10－18.

［8］Mitchell A，Cockcroft S L，Wang T.Some observations on the columnair－to－equiaxed transition in nickel－base superalloy［J］.High temperature materials and processes，2005，24(2):131－138.

［9］Kelkar K M，Mok J，Patankar SV，et al.Computational modeling of electroslag remelting processes［J］.Journal de physique IV，2004，120:421 －428.

［10］Cefalu S A，VanEvery K J，Krane M J M.Modeling of electroslag remelting of Ni－ Cr － Mo alloys［J］.Multiphase phenomena and CFD modeling and simulation in materials processes，2004:279 －288.

［11］劉福斌，姜周華，臧喜民，等.電渣重熔過程渣池流場的數值模擬［J］.東北大學學報，2009，30(7):1013－1017.

The flow field，heat transfer and solidifying of slag and ingot in ESR system

LI Bao-kuan1，LU Qiu-min1，CHEN Rui2，SONG Zhao-wei2，

(1.School of Materials＆Metallurgy，Northeastern University，Shenyang 110004，China;2.Shenyang Research Institute of Foundry，Shenyang 110021，China)

This paper presents a numerical analysis on slag，steel fluid and ingot of ESR using computational fluid dynamics method.According to the results，we get the flow field，heat transfer and solidifying in ESR system.And on this basis，we change melting speeds to see the results of the impact of ESR .With the metal melting speed increases，the temperature of the metal molten pool is rised，and the metal pool became wider.

ESR;temperature field;solidifying field;flow field

TF 744

A

1671-6620(2011)S1-0120-06

2010-10-15.

國家自然科學基金重點項目資助 (50934008).

李寶寬 (1963—)，男，遼寧遼陽人，東北大學教授，E－mail:libk@smm.neu.edu.cn.

book=134,ebook=188