分享:連鑄結(jié)晶器內(nèi)高溫熔體多相流模型化研究進(jìn)展
1.
2.
連鑄結(jié)晶器內(nèi)高溫熔體多相流的運(yùn)動(dòng)狀態(tài)強(qiáng)烈地影響高品質(zhì)鋼的潔凈化、均勻化和精細(xì)化。該高溫熔體多相流為非穩(wěn)態(tài)湍流,同時(shí)耦合傳熱、傳質(zhì)、相變、反應(yīng)、電磁力等諸多過程,形成了十分復(fù)雜的非穩(wěn)態(tài)、非線性、非平衡多物理場,很難通過現(xiàn)場測試對其內(nèi)部的各種物理量進(jìn)行在線測量。在確保流動(dòng)和凝固過程相似的情況下,開展物理模擬實(shí)驗(yàn)以及對其多尺度傳輸現(xiàn)象的數(shù)值模擬,已成為闡述連鑄坯各類缺陷形成機(jī)理的首選。然而,與傳統(tǒng)流體力學(xué)問題相比,冶金高溫熔體多相流的突出特點(diǎn)是物性變化大、本構(gòu)關(guān)系式復(fù)雜、相界面影響因素多及邊界物理量梯度大等,且連鑄結(jié)晶器內(nèi)多物理場存在連續(xù)相大界面變形、離散相顆粒輸運(yùn)以及連續(xù)相-離散相轉(zhuǎn)變等復(fù)雜多變的多尺度界面現(xiàn)象,以及多尺度的湍流渦結(jié)構(gòu),這些都給高溫熔體多相流模型化研究帶來極大困難。相比于單相流,多相流具有相界面拓?fù)湫螤钭兓奶卣?本文從離散流界面尺度分布性、混合流界面跨尺度性、凝固界面多尺度性、以及湍流在揭示多尺度相界面結(jié)構(gòu)中的作用等4方面回顧了連鑄結(jié)晶器內(nèi)高溫熔體多相流模型化的研究進(jìn)展,并展望了未來可能的研究方向。
關(guān)鍵詞:
高品質(zhì)連鑄寬厚板和大圓坯是重要的戰(zhàn)略資源,主要用于國防大型艦艇、核電站安全殼、深海鉆井平臺、高壓鍋爐管、軸承套、高速列車輪等領(lǐng)域。探傷缺陷是連鑄寬厚板的主要質(zhì)量缺陷,形成原因包括氣泡/夾雜物、夾渣、氫致裂紋、中心偏析和疏松等[1~8],通過分析現(xiàn)場大量鋼板探傷檢測圖譜,發(fā)現(xiàn)目前最不受控的是非穩(wěn)態(tài)、非均勻的氣泡/夾雜物和夾渣缺陷[9~14],它們已成為制約連鑄寬厚板質(zhì)量進(jìn)一步提高的瓶頸。目前大圓坯主要采用模鑄的方法制備,但模鑄的生產(chǎn)效率、金屬收得率較低。用連鑄代替模鑄能夠?qū)崿F(xiàn)模具鋼制備的高效率和高收得率,但隨著連鑄坯截面擴(kuò)大,單位長度鑄坯熱容量增大、散熱面積減小,凝固方式由快速凝固向慢速凝固轉(zhuǎn)變,造成芯部鋼液的熱對流和溶質(zhì)元素的再分配過程加劇[15,16],一些由凝固速率慢導(dǎo)致的溶質(zhì)元素宏觀偏析及疏松縮孔等問題也越發(fā)突出。
向煉鋼反應(yīng)器中噴吹惰性氣體攪拌鋼液已成為現(xiàn)代冶金的重要技術(shù)手段。塞棒或上水口吹A(chǔ)r作為保障連鑄結(jié)晶器順行的關(guān)鍵技術(shù),其作用是防止水口堵塞、防止吸氣、促進(jìn)非金屬夾雜物上浮,從而提高鋼坯質(zhì)量。但小氣泡和黏附在其表面的非金屬夾雜物一旦被凝固坯殼捕捉,就會(huì)造成鑄坯“氣泡類”缺陷,惡化產(chǎn)品質(zhì)量;吹A(chǔ)r還會(huì)造成渣/金界面劇烈波動(dòng)并可能發(fā)生乳化,以致于液渣在上回流或漩渦的作用下被卷入鋼液中,形成卷渣。連鑄過程中高溫熔體(鋼液、熔渣、惰性氣體)多相流的運(yùn)動(dòng)狀態(tài)強(qiáng)烈地影響高品質(zhì)鋼的潔凈化、均勻化和精細(xì)化。結(jié)晶器內(nèi)的高溫熔體流動(dòng)為非穩(wěn)態(tài)湍流,同時(shí)耦合傳熱、傳質(zhì)、相變、電磁力、多相流等諸多過程,形成了十分復(fù)雜的非線性、非平衡多物理場,如圖1所示。該場存在連續(xù)相(渣/金界面、凝固前沿等)大界面變形、離散相(氣泡、夾雜物、渣滴、等軸晶粒等)顆粒輸運(yùn)以及連續(xù)相-離散相轉(zhuǎn)變(大氣團(tuán)的形成與破碎、渣滴的形成與再融合等)等復(fù)雜多變的多尺度界面現(xiàn)象,以及多尺度的湍流渦結(jié)構(gòu),常規(guī)的宏觀連續(xù)理論和方法對描述這類多尺度相界面共存及轉(zhuǎn)變問題面臨根本性困難。
圖1

圖1 結(jié)晶器內(nèi)多相、多物理場、多尺度特征
Fig.1 Multiphase, multi-physics, and multi-scale characteristics in the mold (SEN—submerged entry nozzle)
大尺寸的寬厚板和大圓坯造價(jià)高,加上連鑄過程復(fù)雜的高溫動(dòng)態(tài)環(huán)境,直接的實(shí)驗(yàn)研究既不經(jīng)濟(jì)、難度也極大。針對寬厚板和大圓坯研制長期以來過度依賴“制備-測試-制備”反復(fù)循環(huán)的“試錯(cuò)法”,致使研制周期長、耗費(fèi)高、效率低等問題,在確保凝固過程相似的情況下,開展凝固過程的熱模擬實(shí)驗(yàn),以及對其多尺度傳輸現(xiàn)象的數(shù)值模擬,已成為闡述寬厚板和大圓坯各類缺陷形成機(jī)理的首選。相比于單相流,多相流具有相界面拓?fù)湫螤钭兓奶卣?其復(fù)雜多變的界面結(jié)構(gòu),在瞬態(tài)時(shí)間上具有寬廣的空間尺度范圍,即表現(xiàn)出時(shí)間和空間2個(gè)方面的多尺度特征[17]。根據(jù)流動(dòng)形態(tài)與界面尺度的關(guān)系,連鑄結(jié)晶器內(nèi)的高溫熔體多相流的界面多尺度通常涉及幾個(gè)方面的問題:離散流界面尺度分布性、混合流界面跨尺度性、凝固界面多尺度性、以及湍流在揭示多尺度相界面結(jié)構(gòu)中的作用等。
1 離散流界面尺度分布性問題
離散流界面尺度分布性問題源于離散氣泡或非金屬夾雜物的聚并和破碎效應(yīng),造成離散氣泡或非金屬夾雜物具有較大的尺度范圍,并表現(xiàn)為小尺度相界面的不連續(xù)性。例如,初始均一粒徑(近似) Ar氣泡經(jīng)上水口彌散磚或塞棒進(jìn)入水口后,由于受到鋼液湍流的作用,氣泡間會(huì)發(fā)生強(qiáng)烈的聚并和破碎現(xiàn)象,造成進(jìn)入結(jié)晶器內(nèi)的Ar氣泡具有較廣的粒徑范圍,在結(jié)晶器內(nèi)形成復(fù)雜的多尺度泡狀流。冷態(tài)物理實(shí)驗(yàn)被廣泛用于研究結(jié)晶器內(nèi)的多尺度泡狀流[18~25],研究者[22~25]采用水模型實(shí)驗(yàn)分別測量了不同參數(shù)影響下的水口或結(jié)晶器內(nèi)的氣泡粒徑分布,揭示了吹A(chǔ)r連鑄結(jié)晶器中氣泡的分布及運(yùn)動(dòng)規(guī)律,發(fā)現(xiàn)水模型中可采集的氣泡平均粒徑范圍在0.5~5 mm之間。近期,Wu等[25]采用高速攝像機(jī)和ImageJ軟件識別并成功提取了不同粒徑氣泡的運(yùn)動(dòng)軌跡,如圖2[25]所示。采用低熔點(diǎn)合金和惰性氣體模擬連鑄結(jié)晶器中的多尺度泡狀流也是常用的研究方法之一[26~31],如陳芝會(huì)等[26]采用Pb-Sn-Bi合金研究了低頻行波磁場下磁感應(yīng)強(qiáng)度、拉速、吹A(chǔ)r量對氣泡運(yùn)動(dòng)行為及分布的影響;Timmel等[27,28]利用X射線技術(shù)測量了GaInSn合金熔液中Ar氣泡的數(shù)量和粒徑分布,發(fā)現(xiàn)氣泡粒徑在1~10 mm之間;Ren等[29]采用液態(tài)Hg和超聲Doppler測速儀研究了電磁場對結(jié)晶器流場及液面波動(dòng)的影響。然而這些實(shí)驗(yàn)均采用冷態(tài)水模型或低溫合金實(shí)驗(yàn),實(shí)際高溫鋼液內(nèi)的Ar氣泡運(yùn)動(dòng)狀態(tài)仍不可直接觀察和測量。目前,實(shí)際連鑄坯中的氣泡缺陷只能通過對鑄坯成品進(jìn)行檢測得到,部分學(xué)者[12,32~34]采用X射線技術(shù)、鑄坯刨層檢測等手段在凝固后的板坯中發(fā)現(xiàn)了不同尺寸的氣泡類缺陷,氣泡粒徑(Dm)在0.02~3 mm之間,并獲得了氣泡和非金屬夾雜物在鑄坯內(nèi)的三維空間分布,如圖3[12]所示。
圖2

圖2 不同粒徑氣泡的運(yùn)動(dòng)軌跡[25]
Fig.2 Trajectories of bubbles with different sizes[25]
(a) 1-1.5 mm (b) 1.5-2 mm (c) 2-2.5 mm (d) > 2.5 mm
圖3

圖3 實(shí)際連鑄坯內(nèi)氣泡及夾雜物空間分布[12]
Fig.3 Spatial distributions of bubbles and inclusions in actual slab (Dm—bubble size)[12]
近期,本文作者通過水模型實(shí)驗(yàn)觀察到了結(jié)晶器內(nèi)離散氣泡的群體運(yùn)動(dòng)特征,如圖4所示,將氣泡的運(yùn)動(dòng)分布分成2部分:氣泡群和脫群氣泡。氣泡群內(nèi),兩相運(yùn)動(dòng)由氣泡群主導(dǎo),氣泡密度較大,易發(fā)生氣泡間的碰撞聚并;氣泡群外,脫群氣泡的運(yùn)動(dòng)由鋼液主導(dǎo),氣泡彌散分布,難以發(fā)生碰撞聚并。脫群氣泡的數(shù)量決定了被凝固坯殼捕捉的氣泡量,因此減少脫群氣泡的產(chǎn)生對鑄坯質(zhì)量的提升至關(guān)重要。氣泡的這種聚團(tuán)效應(yīng)實(shí)現(xiàn)了微觀顆粒向宏觀結(jié)構(gòu)的轉(zhuǎn)變,使系統(tǒng)行為發(fā)生質(zhì)的改變,其傳遞性能與離散體系截然不同。一般而言,界面現(xiàn)象在這一尺度上發(fā)揮了重要作用。目前結(jié)晶器內(nèi)兩相流實(shí)驗(yàn)研究多偏重于從冷態(tài)、等溫實(shí)驗(yàn)現(xiàn)象出發(fā)去探索泡狀流中的一些表觀參數(shù)變化規(guī)律,如流型圖、液面波動(dòng)、含氣率、離散氣泡的空間分布等;而對結(jié)晶器內(nèi)氣泡群存在的穩(wěn)定性條件以及脫群氣泡的產(chǎn)生機(jī)理仍缺乏深入的認(rèn)識,有很多基本理論問題亟待解決。
圖4

圖4 結(jié)晶器內(nèi)的氣泡群與脫群氣泡
Fig.4 Bubbles cluster and bubbles outside cluster
(a) macro-distribution of bubbles (b) 0 ms (c) 15 ms (d) 30 ms (e) 40 ms
由于多相流體系大多是非線性、非平衡的,從實(shí)驗(yàn)和理論的角度研究體系的時(shí)空多尺度相界面結(jié)構(gòu)還有很大困難,而且高溫鋼液不易直接觀察和測量,因此計(jì)算流體力學(xué)(computational fluid dynamics,CFD)成為研究高溫熔體多相流不可或缺的關(guān)鍵手段,有助于從物理層面理解流動(dòng)機(jī)制并為工程應(yīng)用提供相關(guān)的分布型數(shù)據(jù)。基于連續(xù)介質(zhì)假設(shè)的經(jīng)典力學(xué)方法,存在Euler和Lagrange 2種體系,其中Lagrange體系如離散相模型(discrete phase model,DPM)能夠相對容易地描述離散粒子的運(yùn)動(dòng)行為,因而在求解結(jié)晶器內(nèi)非金屬夾雜物、Ar氣泡的運(yùn)動(dòng)行為及分布的研究中得到廣泛應(yīng)用[35~45],研究者[9~11,35~41]建立了不同機(jī)理的夾雜物/氣泡捕捉模型,揭示了氣泡/夾雜物在結(jié)晶器內(nèi)的運(yùn)動(dòng)和捕捉位置,如圖5[10]所示。Zhang等[44]和Yang等[45]通過在DPM中引入氣泡的聚并和破碎模型,分析了結(jié)晶器內(nèi)的多尺寸泡狀流運(yùn)動(dòng)特征。但Lagrange體系存在2點(diǎn)不足:一是忽略了離散相顆粒體積,無法給出幾何意義上的相界面;二是涉及大量粒子時(shí)對計(jì)算機(jī)性能要求過高。針對小尺度相界面(小于網(wǎng)格尺寸)不連續(xù)性問題,一種基于平均場的Euler-Euler雙流體模型(two-fluid model,TFM)展示了較好的適用性,該模型簡化了小尺度相界面的幾何邊界描述,并在源項(xiàng)中添加了相間作用項(xiàng)。部分研究者[46~53]采用單一氣泡粒徑的TFM研究了結(jié)晶器內(nèi)氣液兩相流動(dòng)行為,得到了氣液兩相流動(dòng)的宏觀相分布規(guī)律,如Liu等[51]通過模擬預(yù)測得到了與實(shí)際現(xiàn)場發(fā)現(xiàn)的渣眼現(xiàn)象吻合較好的預(yù)測結(jié)果,見圖6[51]。
圖5

圖5 夾雜物在凝固坯殼上的捕捉位置[10]
Fig.5 Locations of captured inclusions on solidified shell[10]
(a) 1.5 s (b) 3.0 s (c) 10 s (d) 30 s (e) 100 s
圖6

圖6 結(jié)晶器上表面的渣眼分布[51]
Fig.6 Locations of exposed slag eyes on the top surface of the mold[51]
(a) industrial scene (b) numerical prediction (αg—gas volume fraction)
近年來,基于群數(shù)密度函數(shù)的群體平衡模型(population balance model,PBM)得到了快速發(fā)展,該方法根據(jù)彌散相氣泡的微觀行為機(jī)制建立氣泡相群數(shù)密度平衡方程,將氣泡的微觀現(xiàn)象(破碎、聚并、長大等)與其宏觀屬性(粒徑、表面積等)聯(lián)系起來,是目前研究彌散相系統(tǒng)的有效工具。該方法采用內(nèi)、外2個(gè)坐標(biāo)系描述與氣泡群相關(guān)的變量,內(nèi)坐標(biāo)系是指描述氣泡尺寸、表面積等內(nèi)部屬性的坐標(biāo)系,外坐標(biāo)系是指由流體對流和擴(kuò)散決定的氣泡時(shí)間、空間坐標(biāo)系。研究者將PBM與TFM相結(jié)合分別構(gòu)建了多氣泡組質(zhì)量傳遞(MUltiple SIze Group,MUSIG)模型[54~58]和平均氣泡數(shù)密度(average bubble number density,ABND)模型[59],考慮了不同的氣泡聚并、破碎機(jī)制。2種模型的預(yù)測結(jié)果均能夠揭示結(jié)晶器內(nèi)氣液兩相流的動(dòng)力學(xué)行為在空間、時(shí)間上的變化規(guī)律及氣泡尺寸變化規(guī)律。如圖7所示,分別給出了水模型實(shí)驗(yàn)及2種模型預(yù)測的浸入式水口內(nèi)的含氣率分布特征。通過定性比較,實(shí)驗(yàn)結(jié)果中的水口壁面低含氣率特征和水口出口處高含氣率特征能被較好地呈現(xiàn)在模擬預(yù)測結(jié)果中。
圖7

圖7 實(shí)驗(yàn)與模擬水口內(nèi)含氣率分布對比
Fig.7 Comparisons of gas void fraction profiles inside the SEN between experiment (a) and MUltiple SIze Group (MUSIG) model (b) and average bubble number density (ABND) model (c) simulations
PBM能否準(zhǔn)確地預(yù)測氣泡粒徑分布,關(guān)鍵在于是否能夠建立合理的氣泡聚并和破碎模型。然而氣泡的聚并和破碎機(jī)理相當(dāng)復(fù)雜[60~63],其中,氣泡聚并機(jī)理主要有:湍流渦機(jī)理和大氣泡尾渦機(jī)理;氣泡破碎機(jī)理主要有:湍流渦碰撞機(jī)理、液相剪切力拉伸機(jī)理和大氣泡表面不穩(wěn)定性機(jī)理。而且,以上這些機(jī)理模型均是基于常溫水模型推導(dǎo)而來,是否適用于高溫熔體當(dāng)中的Ar氣泡聚并和破碎行為,還有待于進(jìn)一步深入研究。
2 混合流界面跨尺度性問題
混合流界面跨尺度性現(xiàn)象源于流動(dòng)形態(tài)的轉(zhuǎn)化,即連續(xù)相大尺度界面與離散相小尺度界面的相互轉(zhuǎn)變。在連鑄結(jié)晶器內(nèi),當(dāng)渣/金界面鋼水流速過快、劇烈波動(dòng)或出現(xiàn)漩渦時(shí),保護(hù)渣會(huì)被卷入鋼液形成渣滴(即卷渣),造成了從連續(xù)相渣/金界面向離散相渣滴的轉(zhuǎn)變;反之,渣滴在鋼液流場的影響下會(huì)上浮,當(dāng)再次到達(dá)渣/金界面后很容易與原渣層融合,實(shí)現(xiàn)離散相渣滴到連續(xù)相渣/金界面的轉(zhuǎn)變。大量的冷態(tài)實(shí)驗(yàn)被用于研究結(jié)晶器內(nèi)的卷渣現(xiàn)象和卷渣機(jī)理。其中,結(jié)晶器液面波動(dòng)與卷渣行為聯(lián)系密切,可以較好地反映卷渣程度且易于觀測,因此關(guān)于卷渣的早期研究集中在液面波動(dòng)強(qiáng)度的量化。Teshima等[64]經(jīng)過大量實(shí)驗(yàn)數(shù)據(jù)分析,提出了液面波動(dòng)指數(shù)(F數(shù)),涉及鋼液密度、出鋼量、射流沖擊速率、射流傾角、沖擊深度、結(jié)晶器斷面尺寸等參數(shù),可間接評價(jià)卷渣程度。之后,研究者[65~67]對于不同的鋼種和操作工況歸納出了相應(yīng)的F數(shù)適用范圍。經(jīng)過大量的工業(yè)實(shí)踐和實(shí)驗(yàn)研究,發(fā)現(xiàn)卷渣機(jī)理主要分為3類[68~73]:剪切卷渣、漩渦卷渣和氣泡沖擊卷渣。其中,剪切卷渣是最為常見的卷渣形式,其實(shí)質(zhì)是渣層和鋼液間剪切流動(dòng)引起的相界面不穩(wěn)定性,因此對剪切卷渣影響最為直接的因素是鋼液表面流速,研究者[72,74~77]通過大量實(shí)驗(yàn)研究給出了引起剪切卷渣的臨界表面流速。漩渦卷渣是由結(jié)晶器兩側(cè)不對稱流動(dòng)引起的,在流動(dòng)較弱的一側(cè)會(huì)形成漩渦,漩渦強(qiáng)度較大時(shí),可將保護(hù)渣抽吸進(jìn)入鋼液,如圖8[78]所示,形成原因主要與滑動(dòng)水口、水口傾角、水口偏斜度、水口結(jié)瘤、塞棒位置、中間包底部非對稱流動(dòng)、水口兩側(cè)吹A(chǔ)r不均等因素有關(guān)[69,79~86]。氣泡沖擊卷渣是指當(dāng)氣泡進(jìn)入渣層并在渣/金界面或渣層內(nèi)破裂時(shí),產(chǎn)生的沖擊力致使保護(hù)渣劇烈波動(dòng)并進(jìn)入鋼液[87,88],這一作用與上升氣泡簇的流速相關(guān),而氣泡簇流速與吹氣量直接相關(guān),因此部分研究者將產(chǎn)生卷渣時(shí)的吹氣量定義為臨界吹A(chǔ)r量[88]。
圖8

圖8 水模型實(shí)驗(yàn)中的漩渦卷渣現(xiàn)象[78]
Fig.8 Vortex slag entrapment phenomenon observed in water model experiment[78]
(a) near SEN (b) near the quarter width
另外,結(jié)晶器渣/金界面乳化現(xiàn)象將會(huì)顯著降低渣/金界面的界面性質(zhì)(浸潤性和表面張力等),從而加劇卷渣,惡化鑄坯質(zhì)量。Chung等[89]和Han等[90]采用X射線成像技術(shù)動(dòng)態(tài)觀察了氣泡在金屬液和熔渣中的傳輸過程,考察了氣泡、金屬液滴在渣相中的乳化現(xiàn)象。乳化行為與渣/金、渣/氣界面處的剪切力分布息息相關(guān)[91~94],剪切力引發(fā)的熔渣乳化與卷渣有3種不穩(wěn)定機(jī)制[68,73,74],即Kelvin-Helmholtz不穩(wěn)定性、Tylor-Saffman不穩(wěn)定性和Fluid流動(dòng)不穩(wěn)定性。部分研究者[92,93]通過理論和實(shí)驗(yàn)分析定量描述了結(jié)晶器內(nèi)的乳化和卷渣行為,但主要針對等溫的冷態(tài)水模型,而實(shí)際高溫熔體渣/金界面處存在較大的溫度梯度和濃度梯度,由此引發(fā)的Marangoni效應(yīng)[95] (因界面張力梯度引起的流動(dòng))在渣/金界面處變得突出,進(jìn)而影響乳化和卷渣行為。Zhou等[96]采用實(shí)際高溫熔體實(shí)驗(yàn)剖析了渣/金界面潤濕行為和界面張力之間的內(nèi)在相互影響機(jī)理,其結(jié)果也反映了Marangoni效應(yīng)在渣/金界面處的作用。Liu等[51,52]采用水模型實(shí)驗(yàn)觀察到了渣金界面附近的乳化現(xiàn)象,發(fā)現(xiàn)主要是由氣泡聚集造成的,或可稱為氣-渣乳化現(xiàn)象,如圖9[51]所示。
圖9

圖9 水模型實(shí)驗(yàn)中的乳化現(xiàn)象[51]
Fig.9 Emulsifying phenomenon observed in water model experiment[51]
(a) top surface of slag layer
(b) wide face of the mold
多尺度相界面共存及不連續(xù)性的數(shù)學(xué)描述一直是多相流數(shù)值模擬面臨的最大困難。針對大尺度相界面(遠(yuǎn)大于網(wǎng)格尺度)不連續(xù)性問題,一種基于標(biāo)識函數(shù)的界面捕捉模型(interface-capturing model,ICM)被廣泛采用,VOF和Level Set為當(dāng)前主流的ICM。ICM直接從原始的Navier-Stokes方程出發(fā),采用標(biāo)識函數(shù)統(tǒng)一系統(tǒng)方程的描述,避開了平均場的處理,沒有人為地丟失系統(tǒng)信息,被廣泛應(yīng)用于模擬結(jié)晶器內(nèi)氣-鋼-渣三相界面行為[97~108],目前研究成果主要集中于捕捉渣/金界面的波動(dòng)、渣眼分布、剪切卷渣等方面,如Liu等[51,52]和Li等[106~108]采用ICM方法計(jì)算了結(jié)晶器內(nèi)鋼-渣-氣多相流,復(fù)現(xiàn)了卷渣、渣眼等現(xiàn)象,如圖10[52]所示;并基于簡化的管流模型提出了新的液面波動(dòng)數(shù)——J指數(shù)。雖然采用ICM可以得到結(jié)晶器或鋼包內(nèi)的渣滴尺寸[106~108],但受到網(wǎng)格尺度的限制,無法捕捉小于網(wǎng)格尺度的渣滴;且界面以空間幾何尺度顯式的表征,無隱式的物理尺度表征。另外,ICM中多相速度共場的處理使其在相間作用的描述上凸顯不足,尤其是在界面相間速度差較大時(shí),需要補(bǔ)充界面摩擦力。如前面所述,TFM可描述小尺度相界面(小于網(wǎng)格尺度)的平均物理尺度,但對于大尺度相界面,由于數(shù)值擴(kuò)散導(dǎo)致其無明確的相間幾何邊界。因此,2種模型都無法同時(shí)有效地描述混合流中大小尺度共存的界面結(jié)構(gòu)。耦合ICM和TFM是解決混合流界面跨尺度性問題的思路,但關(guān)鍵是要澄清不同尺度間的相互作用關(guān)系和耦合的原則,建立控制系統(tǒng)的穩(wěn)定性條件及調(diào)控機(jī)理。近期,Sun等[109]發(fā)展了多尺度氣泡橋接模型,實(shí)現(xiàn)了ICM和DPM 2種方法的耦合,即Euler和Lagrange體系下的混合氣泡共存及相互轉(zhuǎn)變,并首先將其應(yīng)用到電解鋁領(lǐng)域,如圖11[109]所示,期待未來將其應(yīng)用到連鑄領(lǐng)域。
圖10

圖10 數(shù)值模擬捕捉到的漩渦卷渣現(xiàn)象[52]
Fig.10 Vortex slag entrapment phenomenon predicted by numerical simulation[52]
(a) interface of slag and steel
(b) 20 mm below the slag-steel interface
(c) flow field in the mold
圖11

圖11 底部視角下的氣泡聚并過程[109]
Fig.11 Transient dumbbell coalescence process from a bottom view[109]
(a) experiment (b) numerical simulation (DPM—discrete phase model, ICM—interface-capturing model)
3 凝固界面多尺度性問題
合金凝固過程是相界面熱和溶質(zhì)傳輸?shù)倪^程,也是一個(gè)耦合擴(kuò)散、熱力學(xué)非平衡、結(jié)晶各向異性和對流多物理場作用的問題,涉及宏觀(熱溶質(zhì)對流、外場強(qiáng)迫對流、夾雜物漂浮等)和微觀(溶質(zhì)析出/擴(kuò)散、局部溶質(zhì)再分配、晶粒形核和生長等)不同尺度間的傳輸,直接導(dǎo)致連鑄坯內(nèi)部缺陷的形成。在確保凝固過程相似的情況下,開展凝固過程的熱模擬實(shí)驗(yàn),是闡述鑄坯凝固組織及內(nèi)部缺陷形成機(jī)理的有效辦法。但金屬的不透明性阻礙了對凝固過程中液相流動(dòng)的實(shí)時(shí)觀測。研究[110~128]表明,某些水溶液(如NH4Cl)和有機(jī)透明物質(zhì)(如丁二腈)的凝固過程與金屬凝固過程非常類似,可以再現(xiàn)金屬凝固過程中許多重要的物理現(xiàn)象。由于它們具有透明性并在接近室溫的溫度區(qū)間內(nèi)凝固,因而便于實(shí)驗(yàn)條件的控制和對凝固過程的直接動(dòng)態(tài)觀測。Tsukamoto等[110]運(yùn)用實(shí)時(shí)觀察法對熔融過程中晶體的生長和形貌變化進(jìn)行了研究。介萬奇和周堯和[111]通過開展NH4Cl水溶液凝固過程的模擬實(shí)驗(yàn),研究了柱狀晶向等軸晶轉(zhuǎn)變(columnar-equiaxed transformation,CET)的條件及規(guī)律。田陸等[126]通過開展有機(jī)透明混合溶液的凝固模型實(shí)驗(yàn),發(fā)現(xiàn)凝固過程中存在著由凝固前沿的下降流和中心上升流組成的蝶形對流。還有部分學(xué)者采用某些低熔點(diǎn)金屬(如Al-Cu合金)并輔助先進(jìn)的測量技術(shù)進(jìn)行了熱態(tài)凝固實(shí)驗(yàn)研究。Li等[127]采用X射線成像技術(shù)實(shí)時(shí)觀測了在磁場作用下定向凝固過程中枝晶的生長行為,發(fā)現(xiàn)熱電磁流動(dòng)將顯著地影響凝固過程中的微觀和宏觀偏析。Niu等[128]和Liu等[129]搭建了一套水冷凝固可視化實(shí)驗(yàn)系統(tǒng),采用NH4Cl、Na2S2O3等有機(jī)溶液考察了結(jié)晶器內(nèi)晶粒微觀結(jié)構(gòu)及宏觀凝固組織演化,發(fā)現(xiàn)中心加入冷源(模擬喂鋼帶)能顯著降低結(jié)晶器內(nèi)部溫度和橫向溫度梯度,從中心發(fā)展的結(jié)晶雨促發(fā)柱狀晶向等軸晶轉(zhuǎn)變,有效阻礙冷卻壁面柱狀晶的生長,顯著提升最終凝固組織中等軸晶比例,如圖12[128,129]所示。
圖12

圖12 凝固實(shí)驗(yàn)捕捉到的宏觀凝固組織演化[128,129]
Fig.12 Evolution of macroscopic solidification structure observed in solidification experiment
(a) device for solidification[129] (b) without inner cooler[128] (c) with inner cooler[128]
由于合金凝固過程的復(fù)雜性,對其多尺度傳輸現(xiàn)象的模型化已經(jīng)成為研究凝固傳輸過程的主要手段。20世紀(jì)60年代末,Flemings等[130]基于固相無擴(kuò)散和液相處于平衡狀態(tài)的假設(shè)導(dǎo)出了著名的局部溶質(zhì)再分配方程,標(biāo)志著凝固過程溶質(zhì)傳輸模型化的開始,該模型的貢獻(xiàn)在于它首次給出了逆偏析、負(fù)偏析、正偏析和通道偏析的統(tǒng)一數(shù)學(xué)模型。20世紀(jì)80年代,Bennon和Incropera[131]采用經(jīng)典混合物理論建立了描述二元系固-液相變過程動(dòng)量、熱量和溶質(zhì)傳輸?shù)倪B續(xù)介質(zhì)模型,該模型假定凝固中的固相和液相充分混合,相與相之間不存在微觀界面,整個(gè)凝固體系視為一連續(xù)介質(zhì)。最近,Li等[132]利用連續(xù)焓-多孔介質(zhì)模型,并考慮固相有限擴(kuò)散,研究了3.3 t鋼錠的凝固及偏析的形成過程。Beckermann和Viskanta[112]采用體積平均方法建立了描述二元合金凝固的體積平均模型(volume-averaged model,VAM)。隨后,Ni和Beckermann[133]提出了兩相VAM,標(biāo)志著凝固傳輸模型的又一重大進(jìn)展。該模型對凝固體系中的固相和液相分別采用一套微分方程描述其傳輸過程,其特點(diǎn)主要在于把微觀組織和界面?zhèn)鬏敩F(xiàn)象與宏觀傳輸方程聯(lián)系起來。Li等[134]采用兩相VAM預(yù)測了鋼鑄錠凝固過程中的宏觀偏析。但該類方法對固相的處理較為簡化,為克服該缺點(diǎn),Wu和Ludwig[135]將固相細(xì)分為柱狀晶和等軸晶,提出了可描述柱狀晶、等軸晶及液相混合凝固的三相VAM,但模型中假定柱狀晶為柱狀,等軸晶為球狀,忽略了晶粒形貌的影響。為此,Wu等[136]又提出了考慮柱狀晶、等軸晶、柱狀晶間液相、等軸晶間液相及枝晶外液相等五相的VAM。該模型可用來預(yù)測鑄錠的宏觀結(jié)構(gòu),包括CET形成、各相體積分?jǐn)?shù)及分布、宏觀偏析等。但該模型引入了一些微觀尺度上的參變量(凝固微觀組織假定),其準(zhǔn)確性還有待確定。凝固微觀組織通常直接決定材料的最終力學(xué)性能和使用性能。通過對鑄坯凝固微觀組織演變的模擬[137~141],能很好地預(yù)測材料的性能。元胞自動(dòng)機(jī)有限元模型(cellular automat finite element,CAFE)是首個(gè)結(jié)合了宏觀熱流計(jì)算和微觀晶粒生長的計(jì)算模型,Li等[140]和Wang等[141]采用CAFE模型對H13模具鋼的電渣重熔過程進(jìn)行了模擬,研究了電磁攪拌、晶粒形核密度和渣池溫度對組織形貌的影響,結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻和較好。近期,Liu等[129]和Niu等[142,143]利用三相VAM等模型研究了有/無內(nèi)部冷源情況下的凝固組織演變行為,發(fā)現(xiàn)內(nèi)部冷源通過加速形核、增強(qiáng)對流和促進(jìn)柱狀晶向等軸晶轉(zhuǎn)變等作用機(jī)制可提高凝固組織中等軸晶的比率,并建立了喂入冷源參數(shù)與等軸晶率的數(shù)學(xué)關(guān)系式,如圖13[129]所示,該研究為連鑄結(jié)晶器喂鋼帶技術(shù)的推廣應(yīng)用奠定了理論基礎(chǔ)。目前研究人員對各個(gè)尺度上凝固現(xiàn)象的研究已相對成熟,但是如何實(shí)現(xiàn)宏-微觀跨尺度凝固傳輸現(xiàn)象的耦合仍是一個(gè)難點(diǎn)。
圖13

圖13 鋼帶初始溫度(
Fig.13 Relation between ratio of equiaxed zone and initial temperature of steel strip[129] (d0—initial thickness of steel strip, D—diffusion coefficient of solute,
4 湍流在揭示多尺度界面結(jié)構(gòu)中的作用
連鑄結(jié)晶器內(nèi)的高溫熔體流動(dòng)處于強(qiáng)烈的非穩(wěn)態(tài)湍動(dòng)狀態(tài),相界面結(jié)構(gòu)受各相物理量的湍流脈動(dòng)及多尺度渦的控制,因此湍流模型的有效性也是多尺度相界面描述的關(guān)鍵。根據(jù)湍流數(shù)值模擬對不同尺度渦的處理方法,相界面可以分為3類:確定型、統(tǒng)計(jì)型和過濾型。表1給出了3種典型的湍流數(shù)值模擬方法。基于直接數(shù)值模擬(direct numerical simulation,DNS)的相界面屬于確定型界面,但DNS除了受計(jì)算機(jī)容量的限制外,對于具有復(fù)雜邊界問題的處理方法還未從理論上完全解決,目前難以應(yīng)用于冶金工程尺度。基于Reynolds平均模擬(Reynolds-averaged Navier-Stokes,RANS)的相界面為統(tǒng)計(jì)型界面,它仍是目前結(jié)晶器內(nèi)多相流研究的主流選擇[40~44,48~51,144~151],但需要強(qiáng)調(diào)的是,此時(shí)相界面的微細(xì)結(jié)構(gòu)以及湍流渦的多尺度效應(yīng)被抹平,得到的是統(tǒng)計(jì)平均意義上的相界面,無法反映受不同尺度渦影響的界面運(yùn)動(dòng)和變形,例如,RANS k-ε (k為湍動(dòng)能,ε為湍動(dòng)能耗散率)模型預(yù)測的薄板坯連鑄結(jié)晶器內(nèi)的對稱渦結(jié)構(gòu),與水模型實(shí)驗(yàn)的非對稱結(jié)果不符[152]。大渦模擬(large eddy simulation,LES)通過過濾操作將湍流運(yùn)動(dòng)分解成大尺度渦和小尺度渦2部分,大尺度渦通過直接求解N-S方程獲得,而小尺度渦則通過建立亞格子模型求解,因此得到的相界面屬于過濾型界面。LES已成功地應(yīng)用于結(jié)晶器內(nèi)的單相流動(dòng)計(jì)算[9~11,36,40,41,152~159],獲得了鋼液的瞬態(tài)非對稱流場特征。研究者[22,35,36,50,100,107,108,160~164]分別將LES模型與ICM或TFM耦合,獲得了結(jié)晶器內(nèi)的氣液兩相瞬態(tài)流動(dòng)特征,如圖14[50]所示,解決了以往湍流模型無法捕捉氣相對液相造成的湍流脈動(dòng)壓力的難題,改善了局部湍流速度參數(shù)、氣體流動(dòng)速率、液體運(yùn)動(dòng)參數(shù)的預(yù)測精度,與實(shí)驗(yàn)測量結(jié)果吻合較好。但目前發(fā)展的所有亞格子模型均假定流動(dòng)為湍流、充分發(fā)展和各向同性,不能適應(yīng)強(qiáng)各向異性湍流的計(jì)算。而且在結(jié)晶器泡狀流中,氣泡的脈動(dòng)運(yùn)動(dòng)會(huì)對高溫熔體的湍流脈動(dòng)產(chǎn)生影響(如氣泡尾跡增大流體小尺度湍流),目前的亞格子模型中未考慮該項(xiàng)的作用。因此,在LES的應(yīng)用中,適用于復(fù)雜各向異性的多相湍流的亞格子模型還有待于從機(jī)理上進(jìn)行更深入的研究。
表1 3種湍流數(shù)值模擬方法的基本方程和特點(diǎn)
Table 1
Method | Continuity | Momentum equation | Solvable scale | Model | Computation |
---|---|---|---|---|---|
equation | requirement | ||||
DNS |
|
|
Eddies of all scales | No | Huge |
|
|||||
LES |
|
|
Eddies of large scales | Subgrid scale model | High |
|
|||||
RANS |
|
|
Eddies of average scale | Time-averaged turbulence model | Low |
|
圖14

圖14 結(jié)晶器內(nèi)鋼液-Ar氣兩相瞬態(tài)流場特征[50]
Fig.14 Characteristics of molten steel-argon gas two-phase transient flow field in the mold[50]
Color online
(a) RANS (b) LES
(c) comparison of numerical simulation with experiment (
5 結(jié)論與展望
連鑄結(jié)晶器是控制鋼水潔凈度的最后環(huán)節(jié),是連鑄設(shè)備的“心臟”。結(jié)晶器內(nèi)的多物理場存在連續(xù)相大界面變形、離散相顆粒輸運(yùn)以及連續(xù)相-離散相轉(zhuǎn)變等復(fù)雜多變的多尺度界面現(xiàn)象,以及多尺度的湍流渦結(jié)構(gòu),形成了復(fù)雜的高溫熔體多相流。連鑄坯的各類質(zhì)量缺陷與結(jié)晶器內(nèi)高溫熔體多相流傳輸現(xiàn)象密切相關(guān)。因此,研究結(jié)晶器內(nèi)高溫熔體多相流傳輸現(xiàn)象對于促進(jìn)連鑄坯質(zhì)量的改進(jìn)具有重要意義,是連鑄領(lǐng)域一直以來的研究重點(diǎn),受到廣泛關(guān)注。然而,面對連鑄過程復(fù)雜的高溫動(dòng)態(tài)環(huán)境,采用直接的實(shí)驗(yàn)研究既不經(jīng)濟(jì)、難度又極大。因此, 計(jì)算流體力學(xué)方法就顯得極為重要。多年來,國內(nèi)外研究人員圍繞結(jié)晶器內(nèi)的高溫熔體多相流模型化開展了大量的研究工作,推動(dòng)了連鑄技術(shù)的快速進(jìn)步。本文針對結(jié)晶器內(nèi)高溫熔體多相流的界面多尺度性,分別從離散流界面尺度分布性、混合流界面跨尺度性、凝固界面多尺度性以及湍流在揭示多尺度相界面結(jié)構(gòu)中的作用等4方面回顧了近年來在該領(lǐng)域的研究現(xiàn)狀。
隨著新型高端裝備向著“大型化、一體化、高性能化”發(fā)展,大截面連鑄坯制造面臨更復(fù)雜的合金體系和更苛刻的組織/缺陷控制要求。然而,隨著連鑄坯的截面擴(kuò)大和高合金化,凝固方式由快速凝固向慢速凝固轉(zhuǎn)變,造成芯部鋼液的熱對流和溶質(zhì)元素的再分配過程加劇,進(jìn)而導(dǎo)致連鑄坯組織不均勻、宏觀偏析、疏松/縮孔、熱裂和雜質(zhì)富集等非均質(zhì)問題日益突出,已成為制約大截面高合金鋼連鑄坯質(zhì)量進(jìn)一步提高的瓶頸。為了滿足冶金產(chǎn)品日益嚴(yán)格的潔凈化、均勻化和精細(xì)化要求,亟需從基礎(chǔ)層面繼續(xù)深入開展連鑄過程高溫熔體多相流模型化研究,主要包括以下幾個(gè)方面。
(1) 氣泡聚并和破碎模型的研究。描述氣泡聚并、破碎微觀機(jī)制的數(shù)學(xué)模型是雙流體群體平衡模型的關(guān)鍵。然而氣泡的聚并和破碎機(jī)理非常復(fù)雜,目前氣泡聚并和破碎模型中仍有許多不確定的可調(diào)參數(shù),尤其是較寬的氣泡尺寸分布范圍以及不同氣泡形狀的影響很難進(jìn)行定量描述。
(2) 混合流界面跨尺度性問題的研究。連鑄過程涉及復(fù)雜的混合流界面跨尺度性問題,目前的數(shù)學(xué)模型尚不能較好地再現(xiàn)該過程。耦合ICM和TFM提供了解決該問題的思路,但關(guān)鍵是要澄清不同尺度間的相互作用關(guān)系和耦合的原則,建立控制系統(tǒng)的穩(wěn)定性條件及調(diào)控機(jī)理。
(3) 凝固過程多相宏-微觀跨尺度傳輸模型的研究。目前對宏觀、微觀上的凝固現(xiàn)象的研究已相對成熟,但是如何將不同尺度上的凝固傳輸現(xiàn)象耦合起來卻是一個(gè)難點(diǎn)。需要建立能夠準(zhǔn)確表征相界面微觀傳輸量的界面?zhèn)鬏斣错?xiàng),進(jìn)而實(shí)現(xiàn)宏-微觀跨尺度凝固傳輸現(xiàn)象的雙向耦合。
(4) 大渦模擬亞格子模型的研究。在氣液兩相流中,氣泡的脈動(dòng)運(yùn)動(dòng)會(huì)對液相的湍流脈動(dòng)產(chǎn)生影響(如氣泡尾跡增大流體小尺度湍流),目前的亞格子模型中未考慮該項(xiàng)的作用。因此,在大渦模擬的實(shí)際應(yīng)用中,適用于復(fù)雜各向異性湍流及多相湍流的亞格子模型還有待于從機(jī)理上進(jìn)行更深入的研究。
來源-金屬學(xué)報(bào)