鋼液內(nèi)宏觀超聲空化特征的數(shù)值模擬

孔 為，蒼大強，王文波（北京科技大學冶金與生態(tài)工程學院，北京 100083）

鋼液內(nèi)宏觀超聲空化特征的數(shù)值模擬

孔 為，蒼大強，王文波（北京科技大學冶金與生態(tài)工程學院，北京 100083）

本文基于CFD軟件FLUENT，通過編寫UDF來模擬鋼液內(nèi)由于超聲波振動所產(chǎn)生的空化特征。結(jié)果表明：未加空化模型的鋼液流場內(nèi)存在正負壓，分別可達到5.5MPa和－5.4MPa；加入空化模型的鋼液流場內(nèi)只存在正壓，且壓力值顯著減小，最大值穩(wěn)定在0.10MPa～0.12MPa，最小值穩(wěn)定在0.08MPa～0.1MPa。兩者變化趨勢正相反。每一時刻都有一部分鋼液區(qū)域不析出氣體。加入空化模型后壓力積分值變化曲線位置隨時間小幅度增高，含氣率積分值曲線變化趨勢與此正相反。各亞松弛因子的變化對空化模擬計算結(jié)果的精確度產(chǎn)生的影響不大。湍流動能、湍流耗散率、密度和體積力等亞松弛因子對低壓流體比較敏感。改變湍流模型和壁面函數(shù)對空化模擬計算結(jié)果基本無影響。不可壓縮氣體的存在有利于空化作用的產(chǎn)生。

鋼液；超聲波；空化；數(shù)值模擬

超聲波能量足夠高時，就會產(chǎn)生“超聲空化”現(xiàn)象。對于空化現(xiàn)象的研究主要包括實驗研究[1]和數(shù)值模擬研究[2～5]。但主要集中在對水[1～5]及低熔點合金熔體[6]方面。關(guān)于空化對鋼液的方面，主要研究空化作用所產(chǎn)生的良好效果[7，8]，諸如改善鋼的凝固組織，均勻成份，改善夾雜，從而提高鋼的力學性能。對鋼液本身的空化特征研究尚屬空白。本文對鋼液內(nèi)宏觀超聲空化特征進行數(shù)值模擬研究，為進一步的試驗研究及超聲波引入到鋼的生產(chǎn)中提供一定的理論指導(dǎo)。

1 數(shù)學模型

模擬對象為剛玉坩堝中的鋼液在超聲波處理振動下的壓力場及空化場變化。坩堝尺寸為?55（內(nèi)徑）mm×120（內(nèi)高）mm，壁厚為3mm。工具頭材質(zhì)為Mo－Al2O3－ZrO2金屬陶瓷管[9]。超聲波處理采用上部導(dǎo)入法，工具頭從上部直接插入鋼液內(nèi)進行施振，液面高度為54mm。超聲波發(fā)生器工作頻率為20kHz。

1.1 計算假設(shè)

（1）假設(shè)整個模擬過程為絕熱的，不考慮熱對流的影響；

（2）假設(shè)鋼液為不可壓縮流體；

（3）忽略空化氣泡崩潰所產(chǎn)生的噴射流；

（4）假設(shè)金屬陶瓷管的物理屬性與剛玉相同。

1.2 主控方程

因為模型可簡化為二維軸對稱模型，于是相應(yīng)的主控方程分別為：

連續(xù)性方程

式中：u、v——分別為極坐標系下軸向（x）、徑向（r）2個方向上的速度分量；

ω——漩渦速度，m/s；

p——壓強，Pa；

Sx、Sr——極坐標軸軸向、徑向2個方向上的動量源分量（本研究中為 0），N/m3；

ρ——密度，kg/m3；

μ——動力粘性系數(shù)，單位為Pa·s。

1.3 定解條件

建立的計算區(qū)域模型如圖1所示。選用非定常求解器，計算0.2ms。湍流模型選擇標準k-ε模型；模型中的坩堝壁及工具頭壁（材質(zhì)均為默認的剛玉）均采用無滑移邊界條件，在近壁區(qū)采用標準壁面函數(shù)；出口條件為壓力出口；工具頭端面為動邊界，結(jié)合前述超聲波發(fā)生器工作參數(shù)，設(shè)置動邊界的位移方程為：

式中a=10μm，f=20kHz，t為超聲波發(fā)生器作用時間。將式（4）結(jié)合超聲工作參數(shù)通過編寫UDF程序加入到Fluent動網(wǎng)格設(shè)置選項中。加入空化模型的同時，采用多相流混合物模型。Fluent所采用的空化模型為Singhal等所發(fā)展的完全空化模型[10]。物性參數(shù)為：鋼液的密度為 7000kg/m3，動力粘性系數(shù)為 6.1·10-3Pa·s；鋼液內(nèi)氣體密度為1.225kg/m3，動力粘性系數(shù)為1.8·10-5Pa·s；氣化壓強為 6Pa。

2 計算結(jié)果與討論

2.1 流場計算結(jié)果與討論

為了研究超聲波振動產(chǎn)生的空化作用與不考慮空化作用對鋼液內(nèi)流場影響的不同。首先排除鋼液內(nèi)次相——氣體的影響，先不加入空化模型。步長為0.5μs，經(jīng)過400步迭代，得出鋼液內(nèi)絕對壓力最大值（如圖2）和最小值（如圖3）隨時間的變化情況?？梢钥闯觯駝舆^程中鋼液內(nèi)絕對壓力變化有很明顯的規(guī)律性。由于工具頭端面的振動頻率是20kHz，振動周期就是0.5×10-4s，所以在迭代400時間步數(shù)，0.2ms中運行了4個周期。鋼液內(nèi)壓力最大值和最小值變化也符合這個振動周期。絕對壓力出現(xiàn)了負壓，正負壓的幅值很大，達到了5.5MPa和-5.4MPa。當壓力最大值達到最高時，壓力最小值也為最高值。當壓力最小值達到最低時，壓力最大值也達到最低值。二者變化趨勢一致。負壓的出現(xiàn)使鋼液內(nèi)局部壓強低于液體的飽和蒸氣壓（氣化壓強），就可以產(chǎn)生空化現(xiàn)象。圖4為計算出的不同時刻壓力云圖?？砂l(fā)現(xiàn)迭代到25步時壓力最大值達到峰值，且分布在工具頭端面下一個很小的區(qū)域，壓力由這個小區(qū)域向外呈輻射狀減小。隨著迭代時間步數(shù)的增加，壓力整體變小，但原來相對高壓區(qū)域的壓力降低速率高于相對低壓區(qū)域壓力降低速率。以致迭代到75步時，原壓力最大值峰值區(qū)域變成壓力最小值峰值區(qū)域，同樣壓力由這個壓力最小值峰值區(qū)域向外呈輻射狀增大。迭代步數(shù)繼續(xù)增加，壓力整體變大，但原來相對低壓區(qū)域壓力增高速率高于相對高壓區(qū)域壓力增高速率。以致迭代到125步時，壓力分布情況重新回到迭代到25步時的情況。繼續(xù)迭代，壓力分布情況如前述循環(huán)變化。為了評價鋼液整體的壓力情況，計算出隨時間變化的鋼液內(nèi)絕對壓力積分值變化曲線，如圖5?？煽吹矫黠@的正弦曲線，完全符合前述的振動周期特征。壓力積分值最大值和最小值的出現(xiàn)時刻符合前述壓力最大值和壓力最小值變化趨勢。

加入空化模型后，同樣采用步長為0.5μs，先迭代400時間步數(shù)，計算出鋼液內(nèi)絕對壓力的最大值（如圖6）和最小值（如圖7）隨運行時間變化的曲線?？煽闯鰞汕€中的壓力全部為正壓，壓力幅值較未加空化模型時大幅降低，且曲線的周期與前述相同；曲線更加圓滑，這說明空化模型的加入可以優(yōu)化計算結(jié)果；壓力最大值變化曲線振蕩幅度隨時間減小。進一步延長迭代時間步數(shù)，發(fā)現(xiàn)曲線的振動幅度隨時間在減小到一定程度后，又逐漸增大，最后振動幅度穩(wěn)定到0.10MPa到0.12MPa之間。在400時間步數(shù)內(nèi)，壓力最小值變化曲線振蕩幅度隨時間逐步增大。進一步延長迭代時間步數(shù)，發(fā)現(xiàn)曲線的振動幅度隨時間在增大到一定程度后，又逐漸減小，最后振動幅度穩(wěn)定到0.08MPa到0.1MPa之間。圖8為迭代到400步時，隨時間變化鋼液內(nèi)含氣率最大值的變化曲線，可看到此曲線的變化規(guī)律與壓力最小值的變化規(guī)律完全相反，進一步延長迭代時間步數(shù)，可知最后曲線振動幅度穩(wěn)定在0.08到0.1。計算得出的各個時刻最小含氣率均為0，這說明每一時刻都有一部分鋼液區(qū)域不析出氣體。計算出不同時刻的壓力云圖和氣化率云圖，可發(fā)現(xiàn)，與圖4相似，也存在分布在工具頭端面下的壓力最大值峰值區(qū)域和壓力最小值峰值區(qū)域，整體壓力隨時間變化情況與未加空化模型時基本一致，但峰值區(qū)域外的其他外圍區(qū)域壓力分布較為均勻，不存在壓力變化的輻射狀分布?？栈瘏^(qū)域隨時間的分布變化情況與壓力分布變化情況正相反。為了評價加入空化模型后的鋼液整體的壓力及空化情況，計算出鋼液內(nèi)絕對壓力積分值變化曲線，如圖9，及含氣率積分值變化曲線，如圖10。與未加空化模型的壓力積分值變化曲線相比，加入空化模型后壓力積分值變化曲線位置會隨時間整體小幅度增高；含氣率積分值變化曲線與壓力積分值變化曲線的變化規(guī)律正相反，曲線整體位置隨時間小幅度降低。

2.2 計算參數(shù)的影響與討論

通過改變算法及物性參數(shù)，可以研究其對加入空化模型的數(shù)值模擬計算結(jié)果的影響。本文主要研究4種參數(shù)的改變對計算結(jié)果的影響：①亞松弛因子；②湍流模型；③壁面函數(shù)；④不可壓縮氣體的質(zhì)量分數(shù)。

首先考慮亞松弛因子的改變，不同的亞松弛因子組合列入表1。

表1 亞松弛因子的不同組合

圖11為迭代到125步時，不同亞松弛因子組合的絕對壓力云圖?？汕逦乜吹?組壓力分布均較為合理，但2、3、4組圖中的中心靠左下位置均出現(xiàn)兩塊三角形白色區(qū)域，這說明此處計算失真?？梢姼鱽喫沙谝蜃拥淖兓鶗τ嬎憬Y(jié)果的精確度產(chǎn)生一定的影響，但影響不大。系統(tǒng)默認的亞松弛因子組合，即第一組為最佳組合。另外各組壓力最低值大小比較情況為：1=4>2>3。這說明2組（即湍流動能、湍流耗散率等亞松弛因子）及3組（即密度、體積力等亞松弛因子）對低壓流體比較敏感，更易產(chǎn)生空化氣泡。

考慮3種不同湍流模型對計算結(jié)果的影響：（1）Standard k-ε 模型；（2）RNG k-ε 模型；（3）Realizable k-ε模型。其中Standard k-ε模型只適合完全湍流的流動過程模擬；RNG k-ε模型是針對高雷諾數(shù)流動問題的；Realizable k-ε與標準k-ε模型相比有2個主要的不同點，即Realizable k-ε模型增加了湍流粘性和為耗散率增加了傳輸方程，這個方程用于處理層流的速度波動。另外由于空化在工具頭端面附近區(qū)域有較強的集中性，因此選擇何種壁面函數(shù)，也可能會影響空化模擬的計算結(jié)果精度。同樣考慮三種不同壁面函數(shù)：（1）Standard；（2）Non-equilibrium；（3）Enhanced。為總體反映各模型及各壁面函數(shù)對計算結(jié)果的影響，均采用計算迭代到125步時各種模型的含氣率積分值作為評價標準。圖12為3種模型的隨時間變化的含氣率積分值變化曲線。圖13為3種壁面函數(shù)隨時間變化的含氣率積分值變化曲線?？煽闯鰞蓤D中的三條曲線均完全重合，這說明這三種模型及三種壁面函數(shù)對計算結(jié)果基本無影響，對一般的鋼液內(nèi)空化模擬問題計算都有較好適應(yīng)性。

液體中往往存在不可壓縮氣體。即使不可壓縮氣體質(zhì)量分數(shù)的微小變化也會對模擬的結(jié)果產(chǎn)生很大的影響。先迭代400步，考慮在三組不同不可壓縮氣體質(zhì)量分數(shù)情況下，對鋼液內(nèi)計算結(jié)果的影響。三組不可壓縮氣體質(zhì)量分數(shù)取值分別為 10-2、10-5、10-8。圖14為這三種不同情況所計算出的含氣率積分值曲線。圖14（1）為三種情況的曲線合圖，可看出隨著不可壓縮氣體質(zhì)量分數(shù)的增加，含氣率也有明顯的提高，而含氣率的增高比例遠大于不可壓縮氣體質(zhì)量分數(shù)的增加值，這說明不可壓縮氣體的存在有利于空化作用的產(chǎn)生，提高空化比例。從圖14（2）、圖 14（3）中，可看出：當不可壓縮氣體質(zhì)量分數(shù)為10-2和10-5時，這兩種的情況的含氣率變化曲線整體形狀改變不大；然而當不可壓縮氣體質(zhì)量分數(shù)降低為10-8時，如圖14（4），曲線形狀發(fā)生很大變化，含氣率變化不再有正弦周期，而是隨時間逐漸變大；增加迭代步數(shù)后，發(fā)現(xiàn)含氣率不再一直增大，也無明顯變化規(guī)律，但最高為0.001。

3 結(jié)論

（1）未加空化模型的鋼液流場內(nèi)存在正負壓，正負壓分別可達到5.5MPa和-5.4MPa，且壓力變化帶有一定的周期性，周期性符合超聲波工作參數(shù)。

（2）加入空化模型的鋼液流場內(nèi)只存在正壓，且壓力值較未加空化模型時顯著減小。壓力最大值穩(wěn)定到0.10MPa到0.12MPa之間，壓力最小值穩(wěn)定在0.08MPa到0.1MPa之間；兩者的變化趨勢正相反。每一時刻都有一部分鋼液區(qū)域不析出氣體。壓力積分值變化曲線位置隨時間整體小幅度增高，含氣率積分值曲線變化趨勢與此正相反。

（3）各亞松弛因子的變化會對空化模擬計算結(jié)果的精確度產(chǎn)生一定的影響，但影響不大。湍流動能、湍流耗散率、密度和體積力等亞松弛因子對低壓流體比較敏感。

（4）改變湍流模型和改變壁面函數(shù)對空化模擬計算結(jié)果基本無影響。

（5）不可壓縮氣體的存在有利于空化作用的產(chǎn)生，提高空化比例。

[1]王國剛，孫冬柏，張秀麗，等.空泡潰滅過程中力學損傷行為[J].北京科技大學學報，2007，29（5）:483-485.

[2]Kim J O，Choi S，Kim J H.Vibroacoustic Characteristics of Ultrasonic Cleaners[J].Applied Acoustics，1999，58（2）:211-228.

[3]Dubus B，Granger C，Mosbah P，et al.Numerical Modeling of the Ultrasonic Cavitation Field and Experimental Evaluation of Bubble Density[A].Proceedings of the 3rd Conference on the Applications of Power Ultrasound in Physical and Chemical Processes[C].Sevilla:Sociedad Espaňola de Acústica，2001:169-174.

[4]Kaltenbacher M，Land-H，Hoffelner J，et al.Use of Modern Simulation for Industrial Applications of High Power Ultrasonics[A].Proceedings of the IEEE Ultrasonics Symposium[C].München:IEEE，2002:673-678.

[5]Lewis J P，Gardner S，Corp I.3D.Simulation of an Ultrasonic Cleaning Vessel-Modelling Methodology and Practical Comparisons[A].In Anglo-French Physical Acoustics Conference [C].Kent:Institute of Physics，2006:1-4.

[6]劉榮光，李曉謙，胡仕成.鋁合金熔體中超聲空化效應(yīng)和氣泡運動過程模擬[J].開發(fā)研究，2007（6），16.

[7]李杰，陳偉慶.超聲波熔體處理對高碳鋼凝固組織和性能的影響.鑄造，2008，57（10）：1009-1011.

[8]李杰，陳偉慶，王曉峰.超聲波功率對高碳鋼中夾雜物的影響.北京科技大學學報，2009，31（9）：1112-1121.

[9]李杰，陳偉慶，何北星，等.超聲波處理高溫鋼液工具頭材質(zhì)的選擇研究[J].北京科技大學學報，2007，29(12)：1246-1249.

[10]A K Singha，l H Y L，i M M Athavale，et a.l Mathematical Basis and Validation of the Full Cavitation Model [J].Journal of Fluid Engineering，2002，124：617-624.

Numerical Simulation of Macroscopic Ultrasonic-Cavitation Features in M olten Steel

KONG Wei，CANG DaQiang，Wang WenBo
（School of Metallurg ical and Ecolog ical Eng ineering，University of Sience and Technology Beijing，Beijing 100083，Beijing China）

The research has been based on CFD software Fluent and cavitation features generated by ultrasonic vib ration in m olten steel have been sim ulated through p rog ramm ing UDF.Results showed that w ithout cavitation model,positive p ressure and negative p ressure were p resent in molten steel and positive p ressure reached 5.5Mpa while negative p ressure-5.4Mpa；w ith cavitation model,there was only positive p ressure and the p ressure value decreased obviously.Maximum value of p ressure stabilized the value between 0.10Mpa and 0.12Mpa.Minimum value of p ressure stab ilized the value between 0.08Mpa and 0.1Mpa.Changes of maximum value of p ressure had a contrary trend w ith changes of m inimum value of p ressure.In every moment,there were som e reg ions w ithout separating out gas in m olten steel.With cavitation model,curve of variations of integral value of absolute p ressure increased slightly,and variations of gas volum e frac tion had a contrary trend w ith variations of integ ral value of absolute p ressure.Underrelaxation factors had some effects on accuracy of calculation results of cavitation simulation though the effec ts were small.Turbulent kinetic energy,turbulent d issipation rate,density and body forces were sensitive to low-p ressure fluid.Changes of turbulent model and wall function had no effects on calculation results of cavitation simulation.No-condensab le gas in existence could be infavor of p roducing cavitation.

Molten steel；Ultrasonic；Cavitation；Numerical simulation

T P391.9;

A；

1006－9658（2011）02－5

2010－10－15

2010－152

孔為（1983－），男，博士研究生