連鑄結晶器內鋼液流動的研究進展

唐海燕 李京社 吉傳波 梁永昌 吳拓 王建磊

( 北京科技大學)

0 引言

結晶器是連鑄機非常重要的部件，是一個強制水冷的無底鋼錠模，稱之為連鑄設備的“心臟”。其基本功能是使鋼液逐漸凝固成所需要規(guī)格、形狀的坯殼，并通過結晶器的振動，使坯殼脫離結晶器壁而不被拉漏以及不產(chǎn)生變形、裂紋等鑄坯缺陷［1］。因此，結晶器的性能對連鑄機的生產(chǎn)能力和鑄坯質量起著十分重要的作用。連鑄板坯的表面和內部缺陷與結晶器內鋼液的流動狀態(tài)密切相關。如從浸入式水口流出的流股沖向鑄坯窄面區(qū)域，由于沖刷作用，使窄面坯殼不均勻生長，影響坯殼的傳熱，從而導致角部縱裂紋產(chǎn)生，嚴重時甚至造成漏鋼;另外結晶器鋼水流動模式的改變引起過大的液面波動，破壞了彎月面的穩(wěn)定凝固，容易造成局部卷渣，從而引起皮下夾渣等表面缺陷。因此，筆者對國內外結晶器內鋼液流動行為進行了總結，并提出了未來的發(fā)展方向。

1 結晶器內鋼水流動模擬研究

結晶器內的流場屬于湍流流動。研究它的方法很多，主要有直接測量法、水模擬法和數(shù)值模擬法。由于鋼水的高溫( ～1800 K) 直接測量方法很難實現(xiàn)。Assar 等［2］用一種電磁感應器來測量彎月面附近鋼水的速度。也有部分學者用釘板試驗法研究結晶器表面液態(tài)渣層的形貌［3］。這些研究對實際連鑄生產(chǎn)過程中結晶器的流場提供了非常有價值的信息。然而，所有這些研究對結晶器內復雜流動現(xiàn)象僅僅提供有限的分析，而且介于現(xiàn)場技術原因有時這些測量結果并不是很準確。因而，水模型和數(shù)值模擬這兩種方法就成為研究連鑄結晶器內流場行為的主要方法。

對結晶器內鋼液流動行為的研究主要集中在鋼液的流場以及鋼/渣界面的物理行為。

1．1 結晶器鋼水流動的水模擬

液態(tài)水和鋼水的動力學粘度系數(shù)相差很小( ～20%) ，因而以相似原理為基礎的水模型實驗研究廣泛應用于模擬連鑄過程流動現(xiàn)象［4－5］。影響結晶器內鋼水流場的主要因素有: 浸入式水口結構參數(shù)如水口直徑、出口傾角、面積、水口形狀等以及連鑄工藝參數(shù)如水口浸入深度、拉速、吹氣量等，因此水模擬實驗也主要是圍繞這些方面展開。模型與原型比例常為1∶ 1、1∶ 2 或者1∶ 3。為便于研究，經(jīng)常使用示蹤法顯示流場。使用較多的示蹤法有: 粒子示蹤法、染色示蹤法及氣泡示蹤法。

對流場的定量測量，分為接觸式測量和非接觸式測量，主要對流場的速度、壓力以及液面波動等進行測量。Thomas 和Huang［6］在1∶ 1 結晶器水模型中用熱線式風速計測量了結晶器表面的水流速度。Honeyands 和Herberton［7］，在薄板坯水模型試驗中用超聲波流量傳感器研究結晶器內瞬態(tài)流動現(xiàn)象，超聲波流量傳感器安裝在水面上部測量表面波動。Gupta 和Lahiri［5］通過流量傳感器的研究表明:對結晶器底部封閉和開放式兩種結晶器，表面波動和形狀波動周期分別為5 s 和50 s，表面波動振幅隨著拉速的增加而增加。水模型試驗中同樣可以利用一些非侵入性測量儀器，如激光多普勒速度測量儀( LDV)［8］和粒子圖像測速儀( PIV)［9］。

關于水口形狀對流場的影響，R． Chaudhary等［10］研究了底部結構為山形和凹形水口對結晶器內流場的影響。結果表明使用山形水口，結晶器的表面流速比凹形水口高50%，因而更易造成表面卷渣。

N．A． Mcphereson［11］比較了橢圓形側孔和圓形側孔的不同之處，結果表明橢圓形側孔的較高部位處存在較大的壓力差，因此產(chǎn)生不穩(wěn)定的水流。而采用圓形出口，鋼液較為穩(wěn)定。

雷洪、朱苗勇等［12］研究了長方形水口和圓形水口對改善卷渣的影響，結果表明:前者比后者更為有效。

N．Tsukamoto 等［13］用1∶ 1 水模型實驗研究了浸入式水口的結構，發(fā)現(xiàn): 在相同出口面積下，方形側孔的流股噴出速度比圓形側孔大，流股沖擊深度也大。但方形側孔的流股散射程度小，對結晶器寬面的沖刷小，利于坯殼生長。

包燕平［14］針對薄板坯連鑄高拉速的需要，開發(fā)了一種耗散型浸入式水口，并通過與普通雙側孔水口的比較，得出了耗散型水口是一種適合薄板連鑄高拉速生產(chǎn)的新型水口。

關于連鑄和水口工藝參數(shù)對結晶器流場的影響，國內外的學者也做了大量的水模擬實驗研究。

T．honeyands 等［15］研究了薄板坯連鑄過程中浸入式水口的結構參數(shù)對射流充填率和傾角的影響。結果表明:在相同的條件下，大內徑的水口射流噴射角比小內徑的水口大。雷洪等［16］通過實驗說明，出口面積增大，明顯降低了水口出口處鋼水流股速度，一方面會減小流股沖擊窄面的速度、使表面波動降低;另一方面，由于鋼流速度偏角的減小，造成沖擊點上移引起更大的波動。

B．G．Thomas［17］等對不同水口傾角下結晶器流場進行了水模型實驗。結果表明: ①隨著水口傾角的下降，有效射流角度增加，射流與窄面的沖擊位置下移;②水口浸入深度小，液面波動劇烈; ③隨著浸入深度變大，射流與窄面的沖擊點下移，但流場的基本形狀不發(fā)生變化;④不同拉速下，結晶器內流場的基本形狀相似;但隨著拉速的增加，結晶器內的流速增大，沖擊深度變深，結晶器液面的波動加劇; ⑤吹入氬氣改變了結晶器上回流的流動方式，使注流與結晶器窄面的撞擊點上移，回流區(qū)的位置上移，當含氣鋼液流至上表面時，流向水口處的回流減緩，向結晶器寬面有輕微的漂移，其影響程度隨含氣率的增加和氣泡尺寸的降低而增加; 但吹入氬氣對下回流的影響很小。

Hua Bai 等［18］研究認為: 氬氣的吹入增大了鋼液的湍動，使射流傾角上傾，減小了回流區(qū)( 鋼液重新流回水口上方) 的大小，且隨著拉速的增加，吹氬的影響減弱。水口中的最小壓力受吹氬量的影響，吹氬量增加，水口中的壓降減少，最小壓力提高，從而減少了水口堵塞的可能性。

陸巧彤等［19］采用1∶ 1 的水模型研究了工藝參數(shù)對大板坯結晶器內流場表面波動的影響。結果表明:隨著拉速和吹氣量的增加，液面波動加劇。但結晶器寬度不同，吹氣量影響也不同，隨著結晶器寬度的增加，液面波動明顯增大。

GUPTA 和A． K． LAHIRI［5］研究了結晶器內液體流動的不對稱性，結果表明:當結晶器的厚寬比為1∶ 6．25 或更小時，水口出口處旋轉的渦流與結晶器寬面的碰撞會導致結晶器內的流動方式不對稱;射流的不穩(wěn)定與渦流方向的改變會引起流動方式的改變。

關于結晶器卷渣，He qinglin［20］利用物理模擬實驗研究了漩渦卷渣的機理，認為結晶器內漩渦卷渣主要是因為在水口附近水口兩側的流股作用形成的，偏流不是漩渦卷渣形成的主要原因，但它會增大漩渦卷渣的深度和發(fā)生頻率。通過優(yōu)化水口參數(shù)，改善偏流現(xiàn)象可以有效地減少或消除漩渦卷渣現(xiàn)象。朱苗勇［21－22］等利用物理模型系統(tǒng)研究了高拉速條件下結晶器內漩渦現(xiàn)象，分析了鋼渣界面行為以及液態(tài)渣的消耗機理。認為表面回流鋼渣卷混是引起結晶器卷渣的主要方式，拉速、水口的浸入深度、張角以及保護渣的粘度均對結晶器內的卷渣產(chǎn)生影響。通過增大水口向下的張角和水口浸入深度，可以防止結晶器內卷渣的發(fā)生。

在結晶器內鋼液流場研究中，最具代表性的水模研究成果是日本的NKK 鋼鐵公司［23］采用1∶ 3水模型研究結晶器液面波動、表面流速等對結晶器卷渣的影響，提出采用F 數(shù)評價結晶器鋼水卷渣情況，F(xiàn) 數(shù)的計算公式為:

式中:ρ——鋼水的密度，kg/m3;

QL——鋼水的流量，m3/s;

θ——注流撞擊窄面的角度;

V——注流的撞擊速度，m/s;

D——撞擊點距自由面之間的距離，m。

研究表明，F(xiàn) 數(shù)在2 ～4 時，結晶器鋼水卷渣最不容易發(fā)生，鑄坯及冷軋軋板的表面缺陷最小。

陸巧彤［24］，齊新霞［25］等用實驗驗證了以上的公式。他們的研究表明: 當F 數(shù)控制在3 ～5，即可以將液面波動控制在±(3 ～5) mm 的合理范圍內，卷渣的幾率最小。

1．2 結晶器內鋼水流動數(shù)值模擬

水模型在直觀認識結晶器內發(fā)生的基本現(xiàn)象和水口優(yōu)化方面發(fā)揮了很重要的作用，但畢竟結晶器內的現(xiàn)象比較復雜，進行結晶器內鋼流流動的計算機模擬很有必要，目前已成為研究的主要手段之一。

流體流動的數(shù)學模型主要涉及以下幾方面內容:①結晶器內流體流動情況、流動的非對稱性、合理的水口設計參數(shù)和工藝參數(shù)優(yōu)化［17，26］; ②研究結晶器自由液面現(xiàn)象，液渣層分布以及流動行為，結晶器液面波動大小以及卷渣情況分析［27－28］;③研究結晶器鋼液流動過程對氣泡運動軌跡［29］以及夾雜物上浮去除的影響［30］;④在電磁力作用下的鋼液流動情況，電磁制動對于結晶器內流動和液面波動的影響［31－32］。

結晶器內的流動屬于湍流流動。研究湍流現(xiàn)象的數(shù)學模型一般包括平均雷諾數(shù)Navier －Stokes 模擬方法( RANS) 、大渦流模擬方法( LES) 和直接數(shù)值模擬方法( DNS)［33］。考慮到計算成本原因，RANS模型加上雙方程( k －ε) 模型廣泛應用于連鑄過程水口和結晶器流動現(xiàn)象的模擬。Huang 等［34］用非穩(wěn)態(tài)RANS 數(shù)學模型研究了水口條件的變化和鋼渣界面波動劇烈的情況下結晶器內瞬態(tài)流動現(xiàn)象。

在研究連鑄湍流現(xiàn)象中，更多的是用LES 和DNS 數(shù)學模型模擬湍流速度場。在雷諾數(shù)較高的情況下由于DNS 計算成本較大，所以LES 模型更加適合研究連鑄湍流現(xiàn)象。但是LES 模型只能解決大尺度渦流，一些耗散的小渦流通常會被過濾掉，這就導致產(chǎn)生殘余應力張量，這時可以用一種叫做次網(wǎng)格尺度的模型( SGS model) 來模擬。LES 模型最初是應用在一些簡單的結構例如在渠道或者管道中流體流動中湍流現(xiàn)象，在這些簡單的流體結構中也使用SGS 模型研究流體的湍流現(xiàn)象［35］。隨著計算機技術的飛速發(fā)展，在一些復雜的流動結構中也可以利用LES 模型來解決更復雜的問題，例如湍流射流現(xiàn)象［36］及鈍體周圍流體流動的湍流現(xiàn)象［37］。

B．G． Thomas［17，26，29］采用高雷諾數(shù)湍流模型對板坯結晶器內鋼液的流動進行了數(shù)值模擬，研究表明，水口的流股角度對結晶器內鋼液的流動有較大影響。

Creech［38］等采用有限差分模塊CFX4．2 對結晶器內穩(wěn)態(tài)流動和傳熱進行模擬，比較了標準K －ε模型、低雷諾數(shù)K－ε 模型和邊界修正K－ε 模型三種模型的仿真計算結果，指出邊界修正K －ε 模型與試驗測試結果吻合較好。而低雷諾數(shù)K －ε 模型對于網(wǎng)格劃分要求較高，只有當網(wǎng)格處理較好時才能保證得到理想的結果。采用Fortran 進行壁面修正的K－ε 模型可以預測結晶器內坯殼生長情況，有助于了解坯殼生長變化規(guī)律，可以避免其它模型出現(xiàn)的坯殼過薄的情況。對于結晶器液面波動的研究，基本上都采用N－S 方程和K －ε 湍流模型［39］，而模擬的難點主要是如何選擇合適的界面邊界條件以及確定界面形狀。對結晶器液面波動進行模擬研究時，主要研究對象局限于水油界面和自由表面。

朱苗勇等［40］用K －ε 模型對不同浸入式水口側孔傾角下結晶器內鋼液的流場進行了數(shù)值模擬研究，結果表明:當出口角度從向下到水平變化時，下部回流區(qū)變大，鋼液對結晶器窄邊的沖擊點上移，上部回流區(qū)變小，鋼流對熔池表面沖擊的強度增強，加劇了表面的波動;當水口出口角度從水平到向上變化時，上部回流區(qū)逐漸變小，直至消失，表面擾動和不穩(wěn)定加劇，易造成卷渣和鋼液面裸露，下部回流中心上移;當出口角度為水平時，射流到結晶器窄面的距離最短，沖擊速度最大。

Hua Bai 等［41］采用商業(yè)軟件CFX 對吹氬后浸入式水口內的鋼液行為進行了研究。結果表明: 大部分氬氣從水口出口的上半部分逸出，而出口下半部分的注流中含有很少的氣體; 氬氣的吹入使出口射流傾角上移，鋼液擾動增加，反流入水口的回流區(qū)減小;隨著拉速的增加，氬氣對流場的影響減小。與水模型的實驗結果比較吻合。

Kouji TAKATANI 等［42］采用LES( Large Eddy Simulation) 模擬法建立了結晶器內瞬變流場的數(shù)學模型，分析了向浸入式水口中吹入氬氣對鋼液和凝固坯殼的影響。他們把計算區(qū)域分為剛體區(qū)域和可移動區(qū)域兩部分，根據(jù)計算得出的壓力初步推測出了彎月面的形狀。

Zhang 等［43］對水口結瘤情況下的結晶器內夾雜物上浮情況進行了數(shù)值模擬研究，認為在水口結瘤情況下，結晶器內水口出口兩流股呈現(xiàn)明顯的不對稱性，未結瘤一側水口出口流股湍動能較強，引起較大的液面波動，容易造成小尺寸的夾雜物卷入鋼液內，同時由于水口結瘤，結晶器內溫度分布不均勻，造成出結晶器時坯殼過薄，嚴重時容易產(chǎn)生漏鋼。

流體流動模型對描述鋼液在結晶器內發(fā)生的湍流流動行為及其對液面卷渣夾雜上浮等現(xiàn)象本質起到了積極作用。但模型一般只考慮了坯殼較薄的結晶器內的流動情況，而且大多只考慮了鋼液流動的情況，未充分考慮兩相區(qū)之間的相互作用，尤其是傳熱和凝固過程對流動行為的影響。

2 結晶器內鋼水流動控制技術

結晶器內的鋼水流動控制除控制浸入式水口參數(shù)和連鑄工藝參數(shù)外，還有幾項很重要的技術就是電磁控制技術。用于板坯結晶器的電磁制動( EMBr) 、電磁攪拌( EMS) 和電磁流動控制( FC 結晶器)是結晶器鋼水流動控制技術的典型代表。

2．1 電磁制動

電磁制動技術( EMBr) 是通過對結晶器施加一個與鑄流方向垂直的靜態(tài)磁場而對流動的鋼液進行制動的技術。它可抑制水口射流速度，促進夾雜物和氣泡上浮。

電磁制動方式有很多種，分別有不同的用途。其中一種是局部電磁制動( local EMBr) ，它是在結晶器兩個寬面上分別施加兩塊磁鐵，在浸入水口出口附近產(chǎn)生一塊矩形區(qū)域磁場。這種類型的電磁制動能夠抑制并彌散水口射流、降低彎月面速度、減小結晶器表面波動［44］，如圖1 所示。

圖1 局部電磁制動工作原理圖解及對流場的作用效果

M．Y．Ha 等［45］用RANS 模型模擬了結晶器內施加電磁制動情況下流場、溫度場及坯殼凝固現(xiàn)象，結果表明由于外加磁場的作用有效抑制了水口射流，減少了射流對結晶器窄面的沖擊作用，降低了底部循環(huán)區(qū)域的穿透深度。Takatani 等［46］使用相似的模擬方法研究電磁制動對結晶器流場的影響。結果表明外加磁場的情況下能夠整體降低結晶器內鋼水的流動速度，并且發(fā)現(xiàn)如果施加的磁場強度過大、產(chǎn)生的洛倫磁力過強，會引起水口射流在沖擊結晶器窄面之前發(fā)生彎曲并消失。Kim 等［47］的研究表明，施加局部電磁制動會使結晶器內鋼水速度和動量有明顯的減小。

另外一種電磁制動是Ruler EMBr( 全幅一段電磁制動) 如圖2 所示。這種類型的EMBr 是在浸入水口出口處結晶器兩側安置兩塊矩形磁鐵薄片，每塊磁鐵都橫過整個結晶器寬面。和Local EMBr 一樣，這種制動方式也是為了穩(wěn)定結晶器彎月面速度和形狀。Harada 等［48］用數(shù)值模擬的方法比較了Local EMBr 和Ruler EMBr 對結晶器內流場影響的不同。研究結果表明盡管兩種電磁制動都能有效降低結晶器彎月面速度和結晶器下部循環(huán)區(qū)域的穿透深度，但Ruler EMBr 能夠更有效的穩(wěn)定彎月面鋼水流動。Zeze 等［49］比較了水銀模型和數(shù)學模型電磁制動對流場的影響，研究表明施加Ruler EMBr 更容易產(chǎn)生活塞流。

圖2 單條全幅EMBr

2．2 電磁攪拌

電磁攪拌是利用交流電產(chǎn)生連續(xù)變化的磁場來控制結晶器內鋼水的流動方式。圖3 為板坯電磁攪拌結晶器，在彎月面附近每個寬面上布置2 個攪拌器，這種電磁攪拌方式能夠使結晶器內鋼水順著結晶器四周有序流動，均勻彎月面處溫度進而改善鑄坯質量［50］。

圖3 板坯結晶器電磁攪拌

2．3 電磁流動控制

電磁流動控制( FC 結晶器) 含有兩個方向相反的制動磁場，第一個位于彎月面區(qū)域，另一個位于結晶器的下部，每一個磁場都覆蓋了板坯的整個寬度。FC 結晶器的上電磁場減小了結晶器彎月面湍流，可防止保護渣卷入凝固殼和角部橫裂; 下電磁場可減小鋼液向下的流速，有利于夾雜物和氣泡上浮［51］，如圖4 所示。

Idogawa et al［52］用數(shù)值模擬的方法研究FC 結晶器內鋼水的流動情況。研究表明FC 結晶器能夠有效降低結晶器內鋼水的整體流動速度。

3 本課題組的研究工作

圖4 FC 結晶器

林曉川等［53］以國內某鋼廠板坯連鑄結晶器為原型，采用1∶ 1的水模型，用OA型光纖式流速儀測定結晶器內流體流速，研究了拉速、浸入式水口出口角度、水口浸入深度、水口底面結構及結晶器斷面寬度等工藝參數(shù)對板坯結晶器內表面流速的影響。其實驗參數(shù)見表1。

表1 實驗工藝參數(shù)

結果表明: 拉速對表面流速的影響最大，隨著拉速的提高，結晶器內鋼液表面流速明顯增大，當斷面寬度為1650 mm，拉速由0．7 m/min 提高到1．4 m/min，表面流速由0．04 m/s 提高到0．1 m/s;波浪面結構的浸入式水口表面流速效果最優(yōu)。

此外，他們通過在結晶器表面不同位置處布置波高傳感器測定瞬時波高，研究了拉速、浸入式水口出口角度、水口浸入深度、水口底面結構等工藝參數(shù)對板坯結晶器內流場和液面波動行為的影響。結果表明: 拉速對板坯連鑄結晶器內液面波動的影響最大，水口傾角次之;各工藝參數(shù)之間對結晶器內部流場的影響是相互制約的;提高拉速時，應適當增大水口浸入深度來控制液面波動。在此實驗基礎上他們提出了優(yōu)化結晶器流場的工藝參數(shù)，并應用于現(xiàn)場取得了一定效果［54］。

陳永峰等［55］以某廠板坯結晶器為原型建立1∶ 1的物理模型，通過監(jiān)測液面波動，研究了浸入式水口結瘤物對結晶器流場造成的影響。實驗結果表明:水口結瘤物的存在會導致結晶器水口兩側流場的明顯不對稱，水口堵塞側鋼液更多集中于結晶器上部，下回流占據(jù)區(qū)域較小;未堵塞一側漩渦卷渣出現(xiàn)的頻率較高;較大拉速下，隨著水口傾角的增大，液面平均波動逐漸減弱; 拉速提高會造成水口兩側平均波高的差距變大，結瘤物脫離水口后液面恢復穩(wěn)定的時間延長;結瘤物脫離瞬間會導致流場迅速惡化，影響鑄坯質量。

付建勛等［56］研究了鑄坯展寬現(xiàn)象及其與拉速的關系。通過對5 個連鑄工廠結晶器設置的調查，對鑄坯展寬的普遍性進行了研究; 在提出衡量鑄坯展寬指標的基礎上，利用儀器在線監(jiān)測鑄坯的展寬。結果顯示: X60 鋼的目標展寬率范圍為1． 25% ～3．63%，平均為3．01%; Q235 鋼目標展寬率范圍為0．77% ～2．91%，平均為2．04%; 且鑄坯的展寬和拉速有明顯的一致性，拉速恒定，展寬也基本保持穩(wěn)定。拉速降低，展寬減小，展寬的變化略滯后于拉速的變化。

吉傳波等［57］采用數(shù)值模擬及鋼水表面插釘法，對電磁制動條件下薄板坯連鑄結晶器內鋼水流動進行研究，分析了不同磁場強度對結晶器內鋼水流動及彎月面波動的影響。模擬結果表明，電磁制動能夠明顯抑制高速流動的鋼水，減小對結晶器窄面的沖擊及鋼水表面卷渣的影響，磁場電流強度由180 A增加到284 A 時，鋼水表面速度最大值由無電磁制動的0．85 m/s 分別降低至0．50 m/s、0．16 m/s。插釘法可作為一種經(jīng)濟有效的方法檢測鋼水表面速度，且測量值與模擬結果比較吻合。通過合理控制磁場電流強度約260 A 可以優(yōu)化結晶器內鋼水流動，進而改善鑄坯質量。

他們也采用雷諾平均( RANS) 數(shù)學模擬方法，研究波浪形和山形水口底部結構對結晶器內鋼水湍流現(xiàn)象及表面流速的影響。表明波浪形水口可以抑制水口流出鋼水的射流，改善結晶器內鋼水流場，降低表面流速，穩(wěn)定液面，進而改善鑄坯表面質量［58］。

結晶器內鋼水渦流現(xiàn)象是連鑄過程中普遍存在的非定常現(xiàn)象。尤其是在高拉速條件下，容易導致連鑄結晶器內強烈的非穩(wěn)態(tài)湍流，而鋼液的湍流狀態(tài)又決定了夾雜物和氣泡的運動分布，進而最終影響鑄坯質量。因此，隨著對鑄坯質量要求的提高，冶金工作者越來越關注這一現(xiàn)象。課題組的吉傳波等人［59］用LES 大渦模擬的方法對水口和結晶器內湍流現(xiàn)象進行了分析，比較了雷諾應力對水口和結晶器內湍流脈動的影響，分析了結晶器寬面中心面內不同位置點的瞬時速度波動，以及采用本征正交分解方法對水口底部和結晶器內射流區(qū)域不同湍流尺度分離，分析不同模態(tài)下速度脈動和能量分布。

研究表明:

1) 水口底部的漩渦流呈現(xiàn)交替旋轉變化，水口出口區(qū)射流的上下波動，引起射流向結晶器內流動呈階梯狀分布;

2) 采用LES 對結晶器內速度模擬比實驗測量結果更能準確捕捉流體內的復雜湍流現(xiàn)象以及速度波動情況;

3) 采用本征正交分解方法對水口和結晶器寬面中心面處、水口底部及射流區(qū)域的速度隨時間和空間變化進行降維處理，能夠識別流動中的主要能量結構狀態(tài)，并通過奇異值分解方法對速度場進行重構，提取不同尺度湍流所積累能量的關系。

4) 控制結晶器內的湍流，最主要的是控制低階模態(tài)下含有高能量的大尺度湍流結構，由于低階模態(tài)下的流體包含主要的能量，通過控制低階模態(tài)下的湍流運動進而控制水口流出的流體向結晶器內擴散作用，為最終優(yōu)化鑄坯質量提供理論支持。

4 結晶器流場研究中國內外的差距及未來的發(fā)展方向

4．1 國內外差距

1) 目前國內所進行的水模擬實驗大多是用粒子、墨汁等作為示蹤劑，簡單地模擬流場情況，而國外已開始使用激光多普勒速度測量儀( LDV) 和粒子圖像測速儀( PIV) ，它們具有準確、清晰、多維等特點。

2) 國內水模擬實驗很少采用1∶ 1 的模型，主要原因是1∶ 1 模型體積太大，需要較大的實驗空間。根據(jù)相似原理模型可以按照比例縮小，但研究表明采用1:1 的模型實驗結果更準確。國外水模型通常按1∶ 1 制作。另外國內結晶器模擬還存在一個問題，就是由于空間限制，通常結晶器模型的長度不夠，造成實驗結果不盡準確。

3) 水模中經(jīng)常用到渣和夾雜物的模擬，通常采用一定配比的油來模擬渣層，利用聚乙烯粒子來模擬夾雜物，它們和水的密度比很難和鋼渣的密度比一致，到目前還沒有很好的辦法來模擬頂渣和夾雜物，尤其是液態(tài)夾雜物，所以水模得到的卷渣和夾雜物去除指標也只具有相對比較的價值，還不能和實際生產(chǎn)的情況完全對應起來。國外用玻璃粒子來模擬夾雜物。

4) 在實驗設備方面目前國內主要應用中國水利水電科學院開發(fā)的DJ800 型水工多功能檢測系統(tǒng)，該系統(tǒng)配有電導率儀、波高傳感器和壓力傳感器，能夠分別檢測溶液的電導率、液面的波高和某一位置的壓力。對流體流速的測量僅局限于對結晶器表面局部點速度的測量。國外的一些水模擬實驗室用到了超聲波多普勒測速儀( UDV) 和粒子圖像測速儀( PIV) ，它們具有分辨率高、測速更準確、多分量同時測量等特點，能夠更全面地反映流場全貌。

4．2 未來的發(fā)展方向

在目前結晶器內鋼液流動的研究中，對于瞬態(tài)的、脈動的行為還不能有效預測，對于高溫液態(tài)鋼水的一些測量工作還不能如愿展開，還需要加大科技的投入。

1) 對連鑄結晶器內鋼水流場進行數(shù)值模擬，大部分研究者采用的都是標準k－ε雷諾平均( RANS) 數(shù)學模型，模型對鋼水流場進行時均化處理，總的來說可以預測流場的平均運動狀態(tài)，但是對鋼水湍流中的小尺度脈動不能有效進行計算，因此無法預測到流體隨時間變化的脈動行為;

2) 對縮小比例的結晶器內水模型采用大渦模擬( LES) ，并通過PIV 實驗室測量結果與模擬結果對比，發(fā)現(xiàn)LES 模型可以準確預測結晶器內鋼水流場，并能夠捕捉到小尺度脈動湍流，但是并未對LES模擬結果的具體湍流信息進行分解描述，只是簡單地預測了鋼水的瞬態(tài)流動形態(tài);

3) 對結晶器內流場測量，多是針對水模型中局部速度及結晶器表面波動，對現(xiàn)場高溫條件下鋼水表面速度的測量工作較少;

4) 對連鑄結晶器內施加電磁制動的研究，由于實驗條件所限大多是采用數(shù)值模擬方法預測其效果，而且對薄板坯內鋼水流場模擬多忽略了凝固坯殼的影響，事實上在凝固坯殼界面的質量和動量損失對結晶器內鋼水流場影響較大。

鑒于冶金工業(yè)中連鑄結晶器工作環(huán)境的特殊性，因此，未來的發(fā)展方向應該還是通過物理和數(shù)值模擬的方法研究結晶器內鋼液的流場。為了準確顯示結晶器內流場運動情況，物理模擬力求采用1∶ 1的模型，并配合多種測量手段，以獲得結晶器內三維空間上的流速分布。數(shù)學模擬則應開發(fā)計算速度更快、能對結晶器不同尺度的湍流信息進行分解描述的模型，如由本課題組開發(fā)的Filter URANS 模型，采用了本征正交分析( POD) 方法，對水口和結晶器內流場進行降維重構，提取不同尺度湍流所積累能量的關系，并且識別流動中的主要能量結構狀態(tài)。

5 結束語

結晶器內鋼液的流動狀態(tài)與連鑄坯的表面和內部缺陷密切相關。國內外主要通過水模擬配合多種測量手段以及數(shù)值模擬研究結晶器內鋼液的流動狀態(tài)，考察各種因素對流動的影響。水模擬具有直觀、可動態(tài)實時地進行控制測量等特點，但受模型比例及水的物性限制，且冷態(tài)和熱態(tài)對流場也有一定的影響，尤其是在吹氣條件下。數(shù)值模擬具有速度快、成本低等特點，但其中網(wǎng)格的劃分、邊界條件和殘差的設置、求解方法的選擇對計算結果影響很大。如何進一步提高物理模擬和數(shù)學模擬的準確性是未來的主要發(fā)展方向。

［1］陳登福，張大江．板坯連鑄結晶器內液面波動數(shù)值模擬［J］．過程工程學報，2010，10(1) :206 －210．

［2］M． B． Assar，P． H． Dauby，and G． D． Lawson． Opening the Black Box: PIV and MFC Measurements in a Continuous Caster Mold［C］．The 83rd Steelmaking Conference Proceedings，397 －411．

［3］H． Nam，H． S． Park and J． K． Yoon． Numerical Analysis of Fluid Flow and Heat Transfer in the Funnel Type Mold of a Thin Slab Caster［J］．ISIJ International，2000，40(9) :886 －892．

［4］Gupta，D．a(chǎn)nd Lahiri，A．K．．Water Modeling Study of the Jet Characteristics in a Continuous Casting Mold［J］． Steel Research，1992，63(5) :201 －204．

［5］Gupta，D． and Lahiri，A． K． A Water Model Study of the Flow Asymmetry inside a Continuous Slab Casting Mold［J］． Metallurgical and Materials Transactions B，1996，27(5) :757 －764．

［6］B．G．Thomas and X．Huang．Effect of Argon Gas on Fluid Flow in a Continuous Slab Casting Mold［C］． The 76th Steelmaking Conference Proceedings，Iron and Steel Society，Warrendale，PA，1993:273 －289．

［7］T． Honeyands and J． Herbertson． Oscillations in Thin Slab Caster Mold Flows［C］．127th ISIJ Meeting，ISIJ，Tokyo，Japan，1994．

［8］Adrian，R．J．Chapter 4．Laser Velocimetry Fluid Mechanics Measurements［M］．Taylor ＆ Francis，1996:175 －299．

［9］Adrian，R．J．Particle－Imaging Techniques for Experimental Fluid Mechanics［J］．Annual．Rev．Fluid Mech．，1991，23:261 －304．

［10］R．Chaudhary，Go－Gi Lee，B．G．Thomas and Seon－Hoyo Kim．Transient Mold Flow with Well －and Mountain －Bottom Nozzles in Continuous Casting of Steel［J］． Metallurgical and Materials Transaction B，2008，139:870 －884．

［11］N．A．Mcpherson．Continuously Cast Clean Steel［C］．Steelmaking Conference Proceedings，1985:13 －25．

［12］雷洪，朱苗勇，邱同榜．板坯連鑄結晶器流場優(yōu)化［J］． 煉鋼，2000，16(3) :29 －31．

［13］N． Tsukamoto，K． Ichikawa，E． Iida． Improvement of Submerged Nozzle Design Based on Water modeling Examination of Tundish Slide Gate［C］． Steelmaking Conference Proceedings，1991:803－808．

［14］Bao Yanping． Experimental Study of Fluid Flow in Thin Slab Continuous Caster Mould with Water－Model［J］． Journal of University of Science and Technology Beijing．1999，6(1) :15 －18

［15］T． Honeyands，J． Lucas，J． Chambers． Preliminary Modeling of Steel Delivery to Thin Slab Caster Moulds［J］． Steelmaking Conference Proceedings．1992:451 －459 ．

［16］雷洪，朱苗勇，汪濕泉，許海虹，王文忠，樊俊飛，弁濟寧．水口吹氬對結晶器彎月面波動的影響［J］． 中國有色金屬學報，1998，8( 增刊2) :468 －471．

［17］Thomas B G，Mika L J，Najjar F M．Simulation of fluid flow inside a continuous slab casting machine［J］． Metallurgical and Materials Transactions B，1990，21B:387 －400．

［18］Hua Bai，Thomas B G．Turbulent Flow of Liquid Steel and Argon Bubbles Slide－Gate Tundish Nozzles: PartⅡEffect of Operation Conditions and Nozzle Design［J］． Metallurgical and Materials Transactions B，2001，32B:269 －284．

［19］陸巧彤，楊榮光，王新華，等． 板坯連鑄結晶器內液面波動的水模型研究［J］．包頭鋼鐵學院學報，2006，25(1) :13 －17．

［20］Qinglin HE． Observations of Vortex Formation in the Mould of a Continuous Slab Caster［J］．ISIJ International，1993，33(2) :343－345．

［21］于海歧，朱苗勇．板坯連鑄結晶器電磁制動和吹氬過程的多相流動現(xiàn)象［J］．金屬學報，2008(5) :619 －625．

［22］孟祥寧，朱苗勇．高拉速板坯連鑄結晶器液態(tài)渣消耗機理分析［J］．金屬學報，2009(3) :485 －489．

［23］手嶋俊雄，久保田淳，鈴木干雄等．スラブ高速鑄造時の連鑄鑄型內溶鋼流動にぉよぼす鑄造條件の影響［J］． 鐵と鋼，1993，79(5) :576 －582．

［24］陸巧彤，王新華，于會香，張炯明，王萬軍．F 數(shù)計算及其與板坯連鑄結晶器內鋼水卷渣的關系［J］． 北京科技大學學報．2007，29(8) :811 －815．

［25］齊新霞，包燕平．結晶器鋼液卷渣指數(shù)的討論［J］．鋼鐵研究，2005(3) :17 －20．

［26］B．G．Thomas，F(xiàn)．M．Najjar．Finite －Element Modeling of Turbulent Fluid Flow and Heat Transfer in Continuous Casting［J］．Applied Mathematical Modeling，1991，15(5) :226 －243．

［27］G． A．，Panaras，A．，Theodorakakos and G． Bergeles． Numerical Investigation of the Free Surface in a Continuous Casting［J］．Metallurgical and Materials Transactions B． 1998，29B: 1117 －1126．

［28］R．Mcdavid and B．G．Thomas．Flow and Thermal Behavior of the Top Surface Flux/Powder Layers in Continuous Casting Molds［J］．Metallurgical and Materials Transactions B，1996，27B:672－685．

［29］B．G．Thomas，X．Huang and R．C．Sussman．Simulation of Argon Gas Flow Effects in a Continuous Slab Caster［J］． Metallurgical and Materials Transactions B．1994，25B:527 －547．

［30］Lkuo Sawada，Hiroyuki Tanaka and Lemitsu Takigawa．Numerical Study of Fluid Flow in the Continuous Slab Caster［C］．Proceedings of the Sixth International Iron and Steel Congress，1990:334－339．

［31］K．Takatani．Analysis of Heat Transfer and Fluid Flow in the Continuous Casting Mold with Electromagnetic Brake［J］． ISIJ International，1989，29:1063 －1068．

［32］X．Huang and B．G．Thomas．Modeling of steel grade transition in continuous slab casting process［J］． Metallurgical and Materials Transactions B．1993，25B:379 －390．

［33］Pope，S． B． Turbulent Flows［M］． Cambridge University Press，Cambridge UK，2000:771．

［34］Huang，X． and Thomas，B． G． Modeling of Transient Flow Phenomena in Continuous Casting of Steel［J］． Canadian Metall．Quart．，1998，37(304) :197 －212．

［35］Schumann，U．Subgrid Length－Scales for Large Eddy Simulations of Stratified Turbulence［J］．Theoretical and Computational Fluid Dynamics，1991，2:279 －290．

［36］Jahnke，S．，Kornev，N．，Leder，A． and Hassel，E． Les Simulation of Jet Mixing Processes with Heat Transfer in Turbulent Pipe Flow［J］． Chemical Engineering Technology，2004，27( 3) : 243 －248．

［37］Michelassi，V．，Wissink，J．G．，F(xiàn)rohlich，J．a(chǎn)nd Rodi，W．Large－Eddy Simulationof Flow around Low － Pressure Turbine Blade with Incoming Wakes［J］．AIAA Journal，2003，41( 11) :2143 －2156．

［38］D． T． Creech and B． G． Thomas． 3 － D Turbulent Multiphase Modeling of Molten Steel Flow and Heat Transfer in a Continuous Slab Caster［C］．CFX User＇s Conference，Wilmington，DE，1998．

［39］劉和平，王忠英．板坯結晶器液面波動的數(shù)學物理模擬及其特點［J］．鋼鐵研究，2002，2:47 －54．

［40］朱苗勇，劉家奇，肖澤強．板坯連鑄結晶器內鋼液流動過程的模擬仿真［J］．鋼鐵，1996，31(8) :23 －27．

［41］Bai，H．a(chǎn)nd Thomas，B．G．Turbulent Flow of Liquid Steel and Argon Bubbles in Slide － Gate Tundish Nozzles: Part I，Model Development and Validation［J］．Metallurgical and Materials Transactions B，2001，32B(2) :253 －267．

［42］Kouji TAKATANI，Yoshinori TANIZAWA． Mathematical Model for Transient Fluid Flow in a Continuous Casting Mold［J］． ISU International．2001，41(10) :1252 －1261．

［43］LiFeng Zhang，YuFeng Wang and XiangJun Zuo． Flow Transport and Inclusion Motion in Steel Continuous Casting Mold under Submerged Entry Nozzle Clogging Condition［J］． Metallurgical and Materials Transactions B，2008，39B:534 －550．

［44］蔡開科．連鑄技術的進展．鋼鐵工業(yè)的前沿技術－冶金工業(yè)前沿科技信息研究班專集［C］．北京:中國金屬學會，2000:180．

［45］Ha，M． Y．，H． G． Lee and S． H． Seong． Numerical Simulation of Three － Dimensional Flow，Heat Transfer，and Solidification of Steel in Continuous Casting Mold with Electromagnetic Brake［J］． Journal of Materials Processing Technology，2003: 322 －339．

［46］Takatani，K．，K．Nakai，N． Kasai，T．Watanabe，and H． Nakajima． Analysis of Heat Transfer and Fluid Flow in the Continuous Casting Mold with Electromagnetic Brake［J］．ISIJ International，1989，29(12) :1063 －1068．

［47］Kim，D． S．，W． S． Kim，and K． H． Cho． Numerical Simulation of the Coupled Turbulent Flow and Macroscopic Solidification in Continuous Casting with Electromagnetic Brake［J］ISIJ International，2000，40(7) :670 －676．

［48］H．Harada，T．Toh，T．Ishii，K．Kaneko and E． Takeuchi． Effect of Magnetic Field Conditions on the Electromagnetic Braking Efficiency［J］．ISIJ International，2001，41(10) :1236 －1244．

［49］Zeze，M．，H． Harada，E． Takeuchi，and T． Ishii． Application of DC Magnetic Field for the Control of Flow in the Continuous Casting Strand［C］． The 76th Steelmaking Conference，1993，20(11) :267 －272．

［50］Hackl，H． R．，A． F． Lehman，J． E． A． Eriksson，and S． G． Kollberg．Superior Slab Casting Using Electromagnetics in AISTech Preliminary Program．2006．

［51］陳登福，張大江． 板坯連鑄結晶器內液面波動數(shù)值模擬［J］．過程工程學報，2010，10(1) :206 －210．

［52］A．Idogawa，M．Sugizawa，S．Takeuchi，K． Sorimachi and T． Fujii． Control of Molten Steel Flow in Continuous Casting Mold by Two Static Magnetic Fields Imposed on Whole Width［J］．Materials Science and Engineering: A，1993，173 (1 －2) :293 －297．

［53］林曉川，李京社，張興利，梁新維． 板坯連鑄結晶器內鋼液表面流速的水模型研究［J］．中國冶金，2010，20(12) :15 －20．

［54］林曉川，李京社，楊樹峰，孫麗媛，李嬌．板坯連鑄結晶器內液面波動的水模型研究［J］．煉鋼，2011，27(2) :56 －60．

［55］陳永峰，李京社，安麗娜，孫麗媛，劉潤藻． 水口結瘤對結晶器流場影響的物理模擬［J］．工業(yè)加熱，2012，41(5) :38 －41．

［56］付建勛，李京社，張慧，張興中． 鑄坯展寬現(xiàn)象及其與拉速的關系［J］．鋼鐵研究學報，2012，22(2) :9 －12．

［57］C．Ji，J．Li，H．Tang，S．Yang．Effect of EMBr on flow in slab continuous casting mold and evaluation using nail dipping measurement［J］．Steel Research International［J］，2012，83( 6) :546 －553．

［58］李京社，程愛民，高雅巍，吉傳波，孫麗媛． 水口結構對連鑄結晶器內鋼水流動的影響［J］．河南冶金，2012，20(2) :1 －3．

［59］Chuanbo Ji，Jingshe Li，Shufeng yang，Liyuan Sun． Large eddy stimulation of turbulent fluid flow in a liquid metal model of continuous casting［J］． Journal of Iron and Steel Research International［J］．2012，19 (12) :35 －41．