2 數(shù)值計(jì)算程序及網(wǎng)格的有效驗(yàn)證

本文對(duì)液橋的三維網(wǎng)格進(jìn)行了有效性驗(yàn)證，最終采用60×68×60的非均勻計(jì)算網(wǎng)格。為了驗(yàn)證計(jì)算程序的可靠性，本文對(duì)Vizman et al。文中的Rayleigh-Bénard結(jié)構(gòu)的算例進(jìn)行了重復(fù)。(采用導(dǎo)電的鎵熔體材料的相關(guān)參數(shù)為:密度ρ=6。1×10^3 kg/m^3，熱膨脹系數(shù)β=1。26×10^-4 /K，動(dòng)力學(xué)粘度ν=1。907×10^-7 kg/ms，擴(kuò)散系數(shù)κ_T=9。22×10^-6 kg/ms，電導(dǎo)率σ=3。87×10^6 A/Vm)。圖2為二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)下，本文計(jì)算得到的熔體在r=3R/4處周向方向的洛倫茲力的分布曲線(a)和文獻(xiàn)結(jié)果(b)。由圖可知，熔體在r=3R/4處熔體周向的洛倫茲力呈拋物線型分布，并在z=0。5H處達(dá)到最大值。本文的計(jì)算結(jié)果與Vizman et al。的結(jié)果較好地吻合，且計(jì)算所得的周向洛倫茲力的最大值與Vizman et al。通過(guò)數(shù)值計(jì)算和實(shí)驗(yàn)得到的最大值的相對(duì)誤差為1。8%。該結(jié)果驗(yàn)證了本文計(jì)算程序的可靠性。

(a) 本文的結(jié)果 (b) Vizman et al的結(jié)果

圖2 熔體在半徑為r=3R/4處沿其軸向方向的周向洛倫茲力分布

3 結(jié)果及分析

為了深入分析外加非均勻旋轉(zhuǎn)磁場(chǎng)對(duì)液橋熔體對(duì)流的控制作用，本文數(shù)值研究了微重力環(huán)境下，二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)對(duì)三維半浮區(qū)液橋表面張力對(duì)流的影響。文中相關(guān)參數(shù)的取值:旋轉(zhuǎn)磁場(chǎng)強(qiáng)度7 mT，旋轉(zhuǎn)頻率50 Hz;相應(yīng)的無(wú)量綱量參數(shù)Pr=0。01， Ta=3。72×10^5， Re=4。4×10^3;表面張力Ma數(shù)取15~40。

3.1 非均勻旋轉(zhuǎn)磁場(chǎng)對(duì)液橋熔體對(duì)流的影響

二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)作用下，液橋熔體的最大周向?qū)α魉俣鹊玫教岣撸畲筝S向?qū)α魉俣鹊玫捷^好的抑制。圖3描述了Ma數(shù)由15逐漸增加到40的過(guò)程中，二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)(7 mT、50 Hz)對(duì)液橋熔體對(duì)流速度的影響。與無(wú)磁場(chǎng)條件相比，二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)作用下，Ma=30時(shí)液橋熔體的最大周向速度增加了4。75倍，最大軸向速度降低了28。5%;Ma=35時(shí)液橋熔體的最大周向速度增加了3。48倍，最大軸向速度降低了25%;Ma=40時(shí)液橋熔體的最大周向速度增加了2。75倍，最大軸向速度降低了22%。

圖3 無(wú)磁場(chǎng)和非均勻旋轉(zhuǎn)磁場(chǎng)(7 mT、50 Hz)作用下，熔體最大周向速度和最大軸向速度曲線圖

本文所采用的旋轉(zhuǎn)磁場(chǎng)有限長(zhǎng)-α模型考慮到洛侖茲力的所有分量，其中周向洛倫茲力對(duì)熔體對(duì)流的影響最大，即外加旋轉(zhuǎn)磁場(chǎng)會(huì)對(duì)熔體產(chǎn)生周向攪拌作用。數(shù)值結(jié)果表明，微重力下二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)使液橋熔體周向流動(dòng)得到加強(qiáng)，該結(jié)論與理論分析是一致的。此外，熔體軸向?qū)α魉俣鹊玫接行б种疲芯勘砻鳎O對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)對(duì)液橋熔體兩方面的作用效果均有助于控制液橋熔體流，將其轉(zhuǎn)變?yōu)槎S軸對(duì)稱定常流動(dòng)。

3.2 非均勻旋轉(zhuǎn)磁場(chǎng)作用下液橋的對(duì)流特性分析

微重力環(huán)境下，無(wú)磁場(chǎng)時(shí)熔體表面張力對(duì)流隨著Ma數(shù)的增大將會(huì)產(chǎn)生失穩(wěn)，首先由二維軸對(duì)稱定常流轉(zhuǎn)變?yōu)槿S定常流。Ma=40時(shí)，熔體對(duì)流已經(jīng)經(jīng)歷過(guò)第一次失穩(wěn)。此時(shí)液橋熔體對(duì)流為定常三維流動(dòng)。圖4描述了無(wú)磁場(chǎng)條件下，液橋熔體Z=0。5H切面上的溫度等值線和流線圖。從圖4(a)中可以看出，熔體Z=0。5H切面上溫度分布具有明顯的非軸對(duì)稱性，且在自由表面附近分別形成了一對(duì)熱區(qū)和一對(duì)冷區(qū)。同時(shí)，從4(b)中看到，Z=0。5H切面上的流線亦為明顯的渦流結(jié)構(gòu):在自由表面附近，熔體從低溫區(qū)流向高溫區(qū)，這一流動(dòng)與自由表面上軸向表面張力驅(qū)動(dòng)流體從高溫區(qū)流向低溫區(qū)是相反的，結(jié)果表明半浮區(qū)液橋的周向流動(dòng)并不是表面張力驅(qū)動(dòng)的。表面張力只是作為周向流動(dòng)的微弱的反作用力而存在。

(a) 溫度等值線 (b) 流線圖

圖4 無(wú)磁場(chǎng)時(shí)， Ma=40的熔體Z=0。5H切面上的溫度等值線和流線圖

在外加二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)(7 mT、50 Hz)作用下，浮區(qū)液橋失穩(wěn)后的三維非軸對(duì)稱的對(duì)流特性得到有效地改善。圖5為二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)作用下，Ma=40時(shí)的半浮區(qū)液橋在Z=0。5H切面上的溫度等值線和流線圖。從5(a)中可以看出，二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)的周向攪拌作用下，熔體Z=0。5H切面上的溫度以中心軸為中心呈現(xiàn)出明顯的軸對(duì)稱分布特性;熔體在Z=0。5H切面上由液橋中心向自由表面流動(dòng)，形成周向渦流結(jié)構(gòu)，如圖5(b)所示。研究結(jié)果表明，二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)的周向攪拌作用下，Ma=40時(shí)熔體的三維非軸對(duì)稱表面張力對(duì)流得到有效地控制，此時(shí)二維軸對(duì)稱流動(dòng)特性明顯。

(a) 溫度等值線 (b) 流線圖

圖5 非均勻旋轉(zhuǎn)磁場(chǎng)作用下(7 mT， 50Hz)， Ma=40時(shí)熔體Z=0。5H切面上的溫度等值線和流線圖

4 結(jié)論

本文采用有限體積法數(shù)值研究了微重力環(huán)境下，外加二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)(7 mT、50 Hz)對(duì)半浮區(qū)液橋表面張力對(duì)流特性的影響。研究結(jié)果表明，無(wú)磁場(chǎng)作用下，半浮區(qū)液橋隨著Ma數(shù)(15—40)的逐漸增加將會(huì)產(chǎn)生首次失穩(wěn)由二維軸對(duì)稱流動(dòng)轉(zhuǎn)變?yōu)槿S定常流動(dòng)。外加二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)的周向攪拌作用下，液橋熔體的最大周向速度得到加強(qiáng)，同時(shí)熔體最大軸向?qū)α魉俣鹊玫接行б种疲撟饔眯Ч欣诳刂埔簶蛉垠w的三維表面張力對(duì)流并將其轉(zhuǎn)變?yōu)槎S軸對(duì)稱流動(dòng)。因此，外加二極對(duì)非均勻旋轉(zhuǎn)磁場(chǎng)可以作為一種控制液橋熔體表面張力對(duì)流的有效手段。