液橋內熱質耦合對流不穩(wěn)定性及旋轉磁場法控制1)

鄒 勇 朱桂平 李 來 黃護林,2)

*(南京航空航天大學能源與動力學院，南京210016)

?(安徽工業(yè)大學數(shù)理學院，安徽馬鞍山243032)

液橋內熱質耦合對流不穩(wěn)定性及旋轉磁場法控制1)

鄒 勇*,?朱桂平*李 來*黃護林*,2)

*(南京航空航天大學能源與動力學院，南京210016)

?(安徽工業(yè)大學數(shù)理學院，安徽馬鞍山243032)

浮區(qū)法因具有無坩堝接觸污染的生長優(yōu)點而成為生長高完整性和高均勻性單晶材料的重要技術.但熔體中存在的毛細對流會給浮區(qū)法晶體生長帶來極大挑戰(zhàn)，這是由于對流的不穩(wěn)定會導致晶體微觀瑕疵的產生和宏觀條紋等缺陷的形成.為了提高浮區(qū)法生長單晶材料的品質，研究浮區(qū)法晶體生長中毛細對流特性及如何控制其不穩(wěn)定性顯得尤為重要.本文采用數(shù)值模擬的方法對半浮區(qū)液橋內SixGe1?x體系中存在的熱質毛細對流展開研究并施加旋轉磁場對其進行控制.結果表明：純溶質毛細對流表現(xiàn)為二維軸對稱模式，溫度場主要由熱擴散作用決定，而濃度場則由對流和溶質擴散共同支配；純熱毛細對流呈現(xiàn)三維穩(wěn)態(tài)非軸對稱流動，濃度分布與熔體內熱毛細對流的流向密切相關，等溫線在對流較大的區(qū)域發(fā)生彎曲；耦合溶質與熱毛細對流則為三維周期性旋轉振蕩流.施加旋轉磁場后，熔體周向速度沿徑向向外增大，熔體內濃度場和流場均呈現(xiàn)二維軸對稱分布.

浮區(qū)法，毛細對流，表面張力，數(shù)值模擬，旋轉磁場

引言

浮區(qū)單晶生長技術[1-3]作為無坩堝法，是制備單晶的一種重要方法，廣泛用于高精度硅、高溫合金及其他半導體材料的生長[4-5].浮區(qū)法晶體生長是一個包含固體導熱、熔體對流、固液相變等多種傳熱流動方式的復雜熱--質輸運過程.生長過程中的傳熱流動特性直接影響功率消耗、溫度分布、界面形狀等宏觀晶體生長參數(shù)，并最終決定單晶內部的微觀結構和缺陷、雜質的分布.因此，有必要深入研究熔體中的熱質流動特性并對其進行有效控制，以改善單晶質量[6-9].由于浮區(qū)單晶生長中存在著自由表面，表面張力驅動的馬蘭哥尼(Marangoni)對流對熔體中雜質濃度的分布會產生極大的影響[10].單晶中雜質分布的不均勻性是制約晶體質量的主要因素之一，它所導致的材料光學、電學性質的不均勻性將最終損害半導體電子和光電子器件的性能.要想充分利用浮區(qū)結晶法的優(yōu)點生長出高品質的單晶，首先必須了解雜質分布不均勻性的起源和形成，即晶體生長中的雜質分凝現(xiàn)象.晶體生長過程中，生長界面處的雜質濃度分布除了受溫度場支配的生長界面的形狀、穩(wěn)定性的影響之外，還由質量輸運過程決定，即通過擴散過程和伴隨對流過程(動量輸運過程)實現(xiàn).擴散過程由熔區(qū)內的濃度梯度驅動，而動量輸運過程則緊密依賴于熔區(qū)內的流場.

目前，人們對于單純的熱毛細對流穩(wěn)定性從實驗[11-13]和理論[14-16]上已有較深入的研究.結果表明：隨著溫差增大，馬蘭哥尼對流將發(fā)生失穩(wěn)并出現(xiàn)熱流體波不穩(wěn)定性；進一步增大溫差，流場將出現(xiàn)更為復雜的振蕩流動，其振蕩特性與馬蘭哥尼數(shù)Ma有關.在地面條件下，熱毛細對流還將和重力產生的自然對流耦合在一起[17-19].吳勇強等[18]采用2cst硅油實驗研究了矮液橋浮力--熱毛細對流的起振及轉捩到混沌的過程，并分析了不同液橋的溫度振蕩頻率及相位變化情況.結果發(fā)現(xiàn)，液橋高徑比、體積比對起振溫差及振蕩頻率都有很大的影響.隨后，王佳等[19]對大尺寸液橋內部的浮力--熱毛細對流流場結構和流動模式展開研究.實驗發(fā)現(xiàn)在臨界Ma附近，流場內會出現(xiàn)行波現(xiàn)象，流動模式也會隨高徑比的變化而發(fā)生變化，繼續(xù)增大Ma，流動會進入混沌狀態(tài).

實際上，由于半導體晶體的摻雜需要，尤其是合金材料，在浮區(qū)生長中，熔體中存在著很大的濃度梯度，這將引起強烈的宏觀分凝.表面處的濃度梯度誘導溶質毛細對流，這種毛細流和由溫度梯度產生的熱毛細對流耦合在一起，形成復雜的流動.Lyubimova等[20-21]利用線性穩(wěn)定性方法分析了浮區(qū)晶體生長中不同情況下三維熱--質耦合對流的發(fā)展過程，獲得了非軸對稱對流隨時間演化的特性和結構圖像.從他們得到的不同結晶速度和高徑比參數(shù)下熱馬蘭哥尼數(shù)MaT--溶質馬蘭哥尼數(shù)MaC的穩(wěn)定性平面圖可以看出：即使是弱的溶質對流擾動也會引起熱毛細流動穩(wěn)定性的大幅變化；弱的熱毛細效應的存在同樣會對溶質毛細流的穩(wěn)定性造成較大的影響.Minakuchi等[22]以全浮區(qū)液橋模型研究了零重力下溶質--熱毛細耦合對流，得到結論是盡管MaC比MaT大，但是自由表面上由濃度梯度引起的軸向對流強度反而比熱毛細對流小.在國內，游仁然和胡文瑞[23]首先采用了液橋模型開展耦合熱--質對流的研究，他們發(fā)現(xiàn)MaC對液橋中的流場和濃度場具有重要影響，但對溫度場分布的影響較弱.Zhou和Huai[24]對具有變界面的液橋中的熱毛細對流和溶質毛細對流進行了數(shù)值分析，研究發(fā)現(xiàn)，當MaT與MaC的比值為?1時，自由表面在冷、熱兩端凸出而中間收縮.目前，這些研究的重點主要集中在毛細對流的特征上，在毛細對流的控制方面的研究相對較少.

由于硅、鍺等熔體具有良好的導電性，通過外加磁場可有效地抑制熔體流動并影響熱質傳輸，從而減少晶體的生長條紋，控制晶體生長的界面、雜質的分離等.各種外部磁場[25-29]中，旋轉磁場具有低功耗的特點，其在導電熔體中產生洛倫茲力不僅能夠增加傳質速率，還能增強熔體流動的穩(wěn)定性，因此施加旋轉磁場對導電熔體進行主動控制是可行的.目前為止，旋轉磁場已用于Czochralski法[30-31]、液相擴散法[32]、Bridgman法[33]及浮區(qū)法[28-29]等晶體生長中，相應的理論分析及實驗研究已經(jīng)展開，取得了一定的成果.鑒于理論分析和實驗研究在高溫條件下存在的困難，目前僅有少量實驗報道溶質毛細對流對Si-Ge浮區(qū)晶體生長的不利影響[34]，利用旋轉磁場對Si-Ge體系中毛細對流進行控制的研究尚未見報道.本文采用半浮區(qū)液橋模型，通過數(shù)值模擬研究SixGe1?x體系中熱質毛細力對熔體對流的貢獻，對比分析各種毛細對流的特征以及旋轉磁場下熔體的傳熱傳質流動特性.

1 物理數(shù)學模型

本文計算采用了柱坐標系(r,θ,z)下的半浮區(qū)液橋模型，物理模型如圖1所示.在零重力條件下，不考慮自由面的變形，即認為液橋自由面為圓柱面，液橋高為L，懸浮在具有相同半徑R的圓盤之間，將固液界面簡化為無滑移平面，分別位于z=0和z=L處并保持溫差?T(?T=Th?Tc)不變.

假設熔體為不可壓縮的牛頓型黏性流體，表面張力作用在自由表面上，其大小可表示為溫度及濃度的線性函數(shù)[24]

式中σ0=σ(T0,C0)，γT=?σ/?T，γC=?σ/?C,C為濃度.只考慮表面張力隨溫度升高而減小，隨溶質濃度增大而增大這一特殊情況.假設流體的其他物性不隨溫度變化而僅隨濃度變化，且與濃度呈簡單的線性混合規(guī)律.施加的外部旋轉磁場強度和頻率分別為B0和ω，其分布可表示為[28-29]

式中，eθ和er分別是軸向和徑向單位矢量.

基于以上假設，外加旋轉磁場作用下，熔區(qū)內磁流體動力學方程組可表示為

式中，?是在柱坐標系下的算符，V是速度矢量，T為熔體溫度，C為熔體中Si的質量濃度.ρ為熔體密度，P為壓力，μ為動力黏度，κ為導熱系數(shù)，cp為比定壓熱容，D為Si在Ge中的擴散系數(shù).假設混合熔體的物性與Si濃度呈簡單的線性關系，即ρ=CρSi+(1 ?C)ρGe，μ=CμSi+(1 ?C)μGe，κ=CκSi+(1?C)κGe，cp=CcpSi+(1?C)cpGe.旋轉磁場在熔體內部產生的洛倫茲體積力可表示為[35]其中σe為電導率.

邊界條件和初始條件設定如下：

當z=L時

當z=0時

當r=R時

當t=0時

式中，n是自由表面的法向矢.

表1 SixGe1?x體系物性參數(shù)[34,36]及其他計算所需參數(shù)Table 1 Physical properties of the SixGe1?xmelt[34,36]and other parameters needed

2 計算方法

采用非結構化網(wǎng)格的有限體積法對控制方程進行空間離散，并在自由表面及交界面附近進行局部加密，網(wǎng)格獨立性驗證曲線如圖2所示，其對應的網(wǎng)格量如表2所示.由圖2可見，當網(wǎng)格數(shù)大于20萬時，液橋內監(jiān)測點P(r=R,θ=0,z=L/2)的無量綱溫度最大值變化較小，綜合考慮計算準確性和計算時間，選擇40×80×80(r向×θ向×z向)的網(wǎng)格進行數(shù)值計算.計算中采用時間步長為10?3s.對動量方程、能量方程中的對流項采用QUICK格式離散，擴散項都采用二階中心差分離散，時間項采用二階隱式推進法，壓力速度耦合采用PISO算法.算法的正確性同文獻[37]的結果進行了比較.表3顯示的是當Pr=0.01,Re=3500時，相同模型下熔體中最大速度值umax和周向最大速度值wmax同文獻[37]的對比結果，可以看出，兩者基本吻合.

圖2 網(wǎng)格獨立性驗證曲線Fig.2 Grid independence veri fi cation curve

表2 網(wǎng)格量Table 2 Grid numbers with di ff erent meshes

表3 本文和文獻[37]的計算結果對比Table 3 Comparison of the present result and the result of Ref.[37]

3 結果與分析

3.1 無磁場時熔區(qū)內毛細對流特征

零重力環(huán)境下，由重力引起的浮力對流消失，晶體生長過程中熔體質量的輸運，主要依賴擴散，因而可實現(xiàn)純擴散晶體生長過程以獲得較大的濃度梯度，同時又可以避免不穩(wěn)定浮力流對晶體質量的影響.但在空間浮區(qū)法晶體生長過程中，雖然克服了浮力對流的作用，但由表面張力梯度驅動的馬蘭哥尼對流在熔區(qū)中占據(jù)主導地位.本文中由于所研究的SixGe1?x體系中溶質Si濃度很大，自由表面溶質毛細力的影響亦不容忽視.數(shù)值結果表明，當僅考慮溶質毛細力時，由其驅動的對流沿自由表面從Si濃度較低的下底面流向濃度較高的上底面，在子午面內形成左右兩個對稱的對流渦，并在液橋中心區(qū)域形成逆向回流渦，結果如圖3(a)所示.在擴散作用下，同時受對流影響，子午面上半部Si濃度等值線(虛線所示)呈“W”型分布(由對稱性只顯示右半部)，但是，由于熔體的普朗特數(shù)遠小于施密特數(shù)，溫度等值線仍呈平行分布(圖中左半部點劃線所示)，這與游仁然和胡文瑞[23]得到的關于溶質毛細對流對溫度場影響較弱的結論是一致的.中截面z=L/2上速度軸對稱的分布特點(圖4(a))表明熔體中對流是二維定常流.

圖3 無磁場時子午面上流線(實線)、無量綱溫度(T*)(左邊點劃線)和Si濃度(右邊虛線)分布圖Fig.3 Streamline(solid line),dimensionless temperature(T*)(dash-dot line in the left)and Si concentration(dashed line in the right)distribution on the meridian plane without magnetic fi eld

熔體在熱毛細力的作用下，對流結構明顯不同于溶質毛細對流，由圖3(b)不難看出，雖然熱毛細對流在子午面內也形成左右兩個對稱的對流渦，但并未在液橋中心區(qū)域形成完整的回流渦.兩種對流渦心位置也有很大的不同，溶質毛細對流渦心在中截面下方，而熱毛細對流的渦心在中截面的上方.由于熱毛細對流比溶質毛細對流強度大，所以回流速度也較大，液橋底部Si濃度較低的熔體被輸運到液橋上方區(qū)域，在上底部附近形成較大的濃度梯度，其等值線(虛線所示)也由“W”型變成“π”型.等溫線(點劃線所示)在自由表面附近區(qū)域向上底面彎曲，此時熔體對流由二維軸對稱流動轉變?yōu)槿S穩(wěn)態(tài)非軸對稱流動，其周向波數(shù)m=4，如圖4(b)所示.

圖4 無磁場時中截面z=L/2速度(m/s)分布圖Fig.4 Velocity(m/s)distribution on z=L/2 section without magnetic fi eld

由圖3(c)可見當溶質毛細力和熱毛細力共同作用時，其液橋內對流結構與圖3(b)中所示熱毛細力單獨作用時的對流結構相近，此時對流渦心較熱毛細對流的渦心略有下移，所不同的是在液橋中心部區(qū)域形成了完整的逆向回流渦.此逆向回流渦結構在圖 3(a)所示的僅有溶質毛細力作用下的液橋對流結構中同樣存在，因此為溶質毛細對流結構所特有.然而，由監(jiān)測點P的速度及Si濃度隨時間變化曲線 (圖 5)可以看出：僅溶質毛細力驅動時，P點的速度很小，溶質毛細對流很弱，因此熱毛細對流的強度和耦合毛細對流強度相當.但是在溶質毛細力和熱毛細力共同驅動下，熔體對流呈現(xiàn)三維周期振蕩模式，由傅里葉頻譜分析可知其振蕩頻率f=0.0292Hz(τ=34.2s).溶質毛細流對表面處濃度影響很大，從圖5(b)可見，雖然熱毛細對流的強度比溶質毛細流的強度大，但是溶質毛細流驅動下P點的Si濃度比熱毛細流下P點的濃度要小很多，耦合后P點濃度介于兩者之間并隨時間呈周期性變化.圖6顯示了中截面z=L/2上Si濃度在一個周期內的振蕩圖樣，其周向波數(shù)m=4，與熱毛細對流的周向波數(shù)相同.由此可見，即使是弱的溶質對流也會引起熱毛細對流穩(wěn)定性的大幅變化，使穩(wěn)態(tài)的對流轉變?yōu)橹芷谛哉袷幜鳎@一點可以從文獻[20-21]的結論得以驗證，但與文獻[38]的結論不同，他們認為當熱毛細對流與溶質毛細對流耦合時，對流的穩(wěn)定性反而比僅有溶質毛細對流的穩(wěn)定性要高.

圖5 無磁場時監(jiān)測點P的速度及Si濃度隨時間變化曲線Fig.5 Curves of velocity and Si concentration of point P with time without magnetic fi eld

圖6 溶質與熱毛細力共同作用時中截面z=L/2上Si濃度周期性振蕩圖Fig.6 Periodic oscillation of Si concentration on z=L/2 section under the e ff ect of coupled solute and thermo-capillary force without magnetic fi eld

對文獻[20-21,38]中截然不同的論述，我們更傾向于認同文獻[20-21]的結果.就本文而言，溶質毛細對流和熱毛細對流都是沿自由表面從下底面向上底面流動，它們耦合時流動是加強的.分析監(jiān)測點P的數(shù)據(jù)可知，熔體僅由溶質毛細力驅動時的速度V1=0.21mm/s，對應溶質MaC=7.4；僅由熱毛細力驅動時的速度V2=6.75mm/s，對應熱MaT=202.7，兩者速度之比V2/V1=32.這是由于熱毛細力比溶質毛細力要大得多，兩者馬蘭哥尼數(shù)之比MaT/MaC=27.因而，當兩種不同形式的對流疊加時，流動穩(wěn)定性發(fā)生變化，速度呈現(xiàn)小幅振蕩模式.熔體內的濃度則發(fā)生較大幅度的振蕩.前文已經(jīng)分析，溶質熱毛細對流的渦心在熔體的下半部，所以處在中截面上監(jiān)測點P的濃度受下半部低濃度的熔體影響較大，C1=0.175；而純熱毛細對流的渦心在熔體的上半部，導致監(jiān)測點P受上半部影響大，C2=0.495.當兩種對流耦合時，高、低濃度的熔體相混合，不穩(wěn)定性增加，濃度產生較大的振蕩，振蕩頻率和速度振蕩頻率相同，其峰谷值為0.024.

3.2 旋轉磁場對熔區(qū)內毛細對流的影響

通過以上分析可知，熔區(qū)內流場結構對于溶質Si的輸運過程有直接而強烈的影響，因此均勻合理的流場結構有利于改善Si濃度分布的均勻性.在熔區(qū)內施加旋轉磁場，可以激發(fā)熔體運動產生強迫對流.旋轉磁場可以對熔體產生與磁場旋轉方向相同的周向攪拌作用，適當調節(jié)旋轉磁場的磁場強度，可以使該強迫對流成為熔體內的主導運動，從而在晶體生長過程中對傳熱傳質有較好的控制.

圖7顯示，選取的不同磁泰勒數(shù)(Ta在7500～17500之間)時子午面上的流動狀態(tài)均呈軸對稱分布，這是由于旋轉磁場產生的洛倫茲力對熔體攪拌作用的結果.從計算結果看，其溫度等值線為一系列平行線，此時熔體內對流結構與無磁場時各毛細力驅動的對流結構均不同.原來在液橋中心部區(qū)域的逆向回流渦被擠壓到液橋下部，并且隨著磁場的增加，因回流而削弱的表面毛細對流的程度更大，這樣在子午面上的4個頂點附近各形成一個對流渦胞，并且上下兩對對流渦胞關于z軸對稱.在對流作用下，除表面附近外，熔體中大部分區(qū)域Si濃度的等值線走向基本和流線一致，這樣，液橋上部分Si濃度較大的熔體向中截面方向流動，液橋下部分Si濃度較低的熔體也向中截面方向流動，因此在中截面附近形成較大的濃度梯度.

圖7 旋轉磁場下子午面上流線(實線)、無量綱溫度(T*)(左邊點劃線)和Si濃度(右邊虛線)分布圖Fig.7 Streamline(solid line),dimensionless temperature(T*)(dash-dot line in the left)and Si concentration(dashed line in the right)distribution on the meridian plane with rotating magnetic fi eld

圖8 不同磁泰勒數(shù)下中截面z=L/2速度(m/s)等值線圖Fig.8 Velocity(m/s)contour lines on z=L/2 section under di ff erent magnetic Taylor numbers

圖9 不同磁泰勒數(shù)下監(jiān)測點P的速度及Si濃度隨時間變化曲線Fig.9 Curves of velocity and Si concentration of point P with time under di ff erent magnetic Taylor numbers

圖8速度分布顯示周向波被完全抑制，在旋轉磁場洛倫茲力的攪拌下，周向速度關于液橋的中心軸對稱，其大小沿徑向向外增大并在自由表面處達到最大值.圖9顯示監(jiān)測點P的速度及濃度振蕩特征已經(jīng)消失，由此可見，施加3mT(對應Ta=7500)的旋轉磁場所產生的強迫對流已經(jīng)成為熔體內的主導運動，從而熔體由三維振蕩流轉變?yōu)槎S軸對稱流動；增大磁場強度，監(jiān)測點P的速度隨磁泰勒數(shù)近似線性增大.當B0=7mT(Ta=17500)時，測得P的速度為9cm/s，這和Dold等[28]施加7.5mT旋轉磁場生長硅單晶時得到的實驗數(shù)據(jù)一致，只是由于SixGe1?x熔體的物性參數(shù)不同于Si熔體的物性參數(shù)，這一速度略小于文獻[28]得到的速度值11cm/s.雖然監(jiān)測點P點的速度變化很大，但是該處的Si濃度基本相同，并且和純熱毛細流下P點的濃度差別不大.

圖10顯示中截面z=L/2上Si濃度也呈中心對稱分布，濃度等值線呈環(huán)狀結構，Si濃度周向分布的不均勻性被完全消除，表明旋轉磁場的作用可使?jié)舛瘸尸F(xiàn)軸對稱分布的特點.在自由表面附近，由于表面張力流的作用，使得下方的Si濃度低的熔體向上方輸運，所以中截面上的邊緣處Si濃度最小.從圖11中截面上Si濃度隨徑向變化曲線可以看出，濃度從中心處沿徑向向外逐漸降低，但在r/R=0.95處出現(xiàn)一個峰值.結合圖7不難發(fā)現(xiàn)，由于受到液橋下部逆向回流渦的擠壓，毛細對流渦胞的流線在r/R=0.90附近突然向下，然后出現(xiàn)一個峰值.在峰值兩側，軸向速度的方向經(jīng)歷了突變，沿半徑增加的方向，軸向速度由原來的負向變?yōu)檎?，原本中截面上方Si濃度高的熔體向下輸運變?yōu)橄路絊i濃度低的熔體向上輸運，因此在中截面上r/R=0.95處Si濃度出現(xiàn)一個峰值.對比不同磁泰勒數(shù)下中截面上Si濃度沿徑向變化曲線，我們發(fā)現(xiàn)Ta=17500時的濃度曲線較為平坦些，因此，針對文中設定的晶體生長條件，施加7mT(Ta=17500)的旋轉磁場可以對傳質進行更好的控制.

圖10 不同磁泰勒數(shù)下中截面z=L/2上Si濃度分布Fig.10 Si concentration contour lines on z=L/2 section under di ff erent magnetic Taylor numbers

圖11 不同磁泰勒數(shù)下中截面z=L/2上Si濃度隨徑向變化曲線Fig.11 Si concentration on z=L/2 section versus radial direction curves under di ff erent magnetic Taylor numbers

4 結論

本文采用半浮區(qū)液橋模型，用數(shù)值方法研究了溶質毛細對流及熱毛細對流的流動特征，并分析了旋轉磁場對耦合溶質--熱毛細對流流場及濃度場的影響.得到以下結論：

(1)僅受溶質毛細力作用時，熔體在自由表面附近形成低濃度區(qū)域，熔體中部區(qū)域濃度由擴散和對流共同支配，子午面上濃度等值線呈“W”型；溫度場主要由擴散作用決定，呈軸對稱分布；溶質毛細對流為二維軸對稱流.

(2)在純熱毛細力作用下，熔體內濃度場受熱毛細對流的影響很大，濃度分布與熔體內對流流向相關，子午面上濃度等值線呈“π”型分布；溫度等值線在自由表面附近區(qū)域向上底面彎曲，熱毛細對流呈現(xiàn)三維穩(wěn)態(tài)非軸對稱流動，其周向波數(shù)m=4.

(3)耦合溶質--熱毛細對流為三維旋轉振蕩流，振蕩頻率f=0.0292Hz，其周向波數(shù)與熱毛細對流的周向波數(shù)相同，m=4.施加旋轉磁場后，在洛倫茲力的攪拌下，熔體周向速度沿徑向向外增大.在選取的3種不同磁泰勒數(shù)(Ta在7500～17500之間)旋轉磁場作用下，熔體內的周向波均被完全抑制，濃度場和流場均呈現(xiàn)二維軸對稱分布.因此，施加旋轉磁場有助于熔體流動的穩(wěn)定性和濃度分布、溫度分布的均勻性，有利于溶質具有周向對稱分布特點的合金晶體的生長.

1 Otani S,Hirata K,Adachi Y,et al.Floating zone growth and magnetic properties of Y2C two-dimensional electride.Journal of Crystal Growth,2016,454:15-18

2 Hossain MM,Watauchi S,Nagao M,et al.E ff ects of growth parameters on silicon molten zone formed by infrared convergent-heating floating zone method.Journal of Crystal Growth,2017,459:105-111

3 Sabanskis A,Surovovs K,Virbulis J.3D modeling of doping from the atmosphere in fl oating zone silicon crystal growth.Journal of Crystal Growth,2017,457:65-71

4 Muiznieks A,Virbulis J,Ldge A,et al.Floating zone growth of silicon//Peter Rudolph ed.Handbook of Crystal Growth(Second Edition).Boston:Elsevier,2015.241-279

5 Da?bkowska HA,Da?bkowski AB,Hermann R,et al.Floating zone growth of oxides and metallic alloys//Peter Rudolph ed.Handbook of Crystal Growth(Second Edition).Boston:Elsevier,2015.281-329

6 Lyubimova TP,Skuridyn RV.The in fl uence of vibrations on the stability of thermocapillary flow in liquid zone.International Journal of Heat and Mass Transfer,2014,69(2):191-202

7 Lin K,Dold P,Benz KW,Numerical study of in fl uences of buoyancy and solutal Marangoni convection on flow structures in a germanium-silicon fl oating zone.Crystal Research and Technology,2005,40(6):550-556

8 Surovovs K,Muiznieks A,Sabanskis A,et al.Hydrodynamical aspects of the fl oating zone silicon crystal growth process.Journal of Crystal Growth,2014,401:120-123

9 Capper P,Zharikov E.Oscillatory-driven fluid flow control during crystal growth from the melt//Peter Rudolph ed.Handbook of Crystal Growth(Second Edition).Boston:Elsevier,2015.951-993

10 Tsukada T.The role of marangoni convection in crystal Growth//Peter Rudolph ed.Handbook of Crystal Growth(Second Edition).Boston:Elsevier,2015.871-907

11 Schwabe D,Benz S.Thermocapillary flow instabilities in an annulus under microgravity—results of the experiment magia.Advances in Space Research,2002,29:629-638

12 Schwabe D,Cramer A,Schneider J,et al.Experiments on the multiroll-structure of thermocapillary fl ow in side-heated thin liquid layers.Advances in Space Research,1999,24:1367-1373

13 Schwabe D,Zebib A,Sim BC.Oscillatory thermocapillary convection in open cylindrical annuli.Part 1.Experiments under microgravity.Journal of Fluid Mechanics,2003,491:239-258

14 Chang CE,Wilcox WR.Inhomogeneities due to thermocapillary flow in floating zone melting.Journal of Crystal Growth,1975,28:8-12

15 Smith MK,Davis SH.Instabilities of dynamic thermocapillary liquid layers:Part 1.Convective instabilities.Journal of Fluid Mechanics,1983,132:119-144

16 Li K,Matsumoto S,Imaishi N,et al.Marangoni fl ow in fl oating half zone of molten tin.International Journal of Heat and Mass Transfer,2015,83:575-585

17 姜歡,段俐,康琦.矩形液池熱毛細對流轉捩途徑研究.力學學報,2015,47(3):422-429(Jiang Huan,Duan Li,Kang Qi.Study on transition to chaos of thermocapillary convection in a rectangular liquid pool.Chinese Journal of Theoretical and Applied Mechanics,2015,47(3):422-429(in Chinese))

18 吳勇強,段俐,李永強等.大普朗特數(shù)大液橋浮力--熱毛細對流地面實驗.力學學報,2012,44(6):981-989(Wu Yongqiang,Duan Li,Li Yongqiang,et al.Ground experiments of bouyant thermocapillary convection of large scale liquid bridge with large Prandtl number.Chinese Journal of Theoretical and Applied Mechanics,2012,44(6):981-989(in Chinese))

19 王佳,吳笛,段俐等.大尺寸液橋熱毛細對流失穩(wěn)性地面實驗研究.力學學報,2015,47(2):580-586(Wang Jia,Wu Di,Duan Li,et al.Ground experiments of bouyant thermocapillary convection of large scale liquid bridge with large Prandtl number.Chinese Journal of Theoretical and Applied Mechanics,2015,47(2):580-586(in Chinese))

20 Lyubimova TP,Skuridyn RV.Numerical modelling of threedimensional thermo-and solutocapillary-induced fl ows in a fl oating zone during crystal growth.The European Physical Journal Special Topics,2011,192(1):41-46

21 Lyubimova TP,Skuridin RV,Faizrakhmanova IS.Thermo-and soluto-capillary convection in the fl oating zone process in zero gravity conditions.Journal of Crystal Growth,2007,303:274-278

22 Minakuchi H,Okano Y,Dost S.A three-dimensional numerical simulation study of the Marangoni convection occurring in the crystal growthofSixGe1?xbythe fl oat-zonetechniqueinzerogravity.Journal of Crystal Growth,2004,266:140-144

23 游仁然,胡文瑞.浮區(qū)中熱和溶質的毛細對流.半導體學報,1992,13(4):209-216(You Renran,Hu Wenrui.Analysis of thermo-solutal-capillary convection in fl oating zone.Chinese Journal of Semiconductors,1992,13(4):209-216(in Chinese))

24 Zhou XM,Huai XL.Free surface deformation of thermosolutocapillary convection in axisymmetric liquid bridge.Microgravity Science and Technology,2015,27(1):39-47

25 Robertson GD,Oconnor JD.Magnetic fi eld e ff ects on fl oat-zone Si crystal growth:Strong axial fi elds.Journal of Crystal Growth,1986,72:110-122

26 Hirata H,Hoshikawa K,Silicon crystal growth in a cusp magnetic field.Journal of Crystal Growth,1989,96:747-755

27 鄒勇,張銀,唐碩捷等.微重力下靜態(tài)磁場對浮區(qū)法硅單晶生長的影響.人工晶體學報,2015,44(9):2343-2349(Zou Yong,Zhang Yin,Tang Shuojie,et al.E ff ect of static magnetic fi elds on growth of fl oating-zone silicon single crystal under microgravity.Journal of Synthetic Crystals,2015,44(9):2343-2349(in Chinese))

28 Dold P,Cr?ll A,Lichtensteige M,et al.Floating zone growth of silicon in magnetic fi elds:IV.Rotating magnetic fi elds.Journal of Crystal Growth,2001,231:95-106

29 YaoL,ZengZ,LiX,etal.E ff ectsofrotatingmagnetic fi eldsonthermocapillary fl ow in a fl oating half-zone.Journal of Crystal Growth,2011,316:177-184

30 Grants I,Gerbeth G.The suppression of temperature fl uctuations by a rotating magnetic fi eld in a high aspect ratio Czochralski con fi guration.Journal of Crystal Growth,2007,308:290-296

31 Zhou X,Huang H.Numerical simulation of Cz crystal growth in rotating magnetic fi eld with crystal and crucible rotations.Journal of Crystal Growth,2012,340:166-170

32 Yildiz E,Dost S.A numerical simulation study for the combined e ff ect of static and rotating magnetic fi elds in liquid phase di ff usion growth of SiGe.Journal of Crystal Growth,2007,303:279-283

33 Wang L,Shen J,Shang Z,et al.Preparation of gradient material in Sn-Cd peritectic alloy using rotating magnetic fi eld.Journal of Crystal Growth,2013,375:32-38

34 Campbell TA,Schweizer M,Dold P,et al.Float zone growth and characterization of Ge1?xSix(x6 10 at%)single crystals.Journal of Crystal Growth,2001,226:231-239

35 Rakoczy R.Enhancement of solid dissolution process under the influence of rotating magnetic fi eld.Chemical Engineering and Processing:Process Intensi fi cation,2010,49(1):42-50

36 AbbasogluS,SezaiI.Three-dimensionalmodellingofmelt fl owand segregation during Czochralski growth of GexSi1?xsingle crystals.International Journal of Thermal Sciences,2007,46:561-572

37 Levenstam M,Amberg G.Hydrodynamical instabilities of thermocapillary fl ow in a half-zone.Journal of Fluid Mechanics,1995,297:357-372

38 Minakuchi H,Takagi Y,Okano Y,et al.The relative contributions of thermo-solutal Marangoni convections on fl ow patterns in a liquid bridge.Journal of Crystal Growth,2014,385:61-65

INSTABILITY OF COUPLED THERMO-SOLUTE CAPILLARY CONVECTION IN LIQUID BRIDGE AND CONTROL BY ROTATING MAGNETIC FIELD1)

Zou Yong*,?Zhu Guiping*Li Lai*Huang Hulin*,2)
*(College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics，Nanjing210016，China)
?(School of Mathematics and Physics，Anhui University of Technology，Ma’anshan243032，Anhui，China)

Floating zone method is an important technology for growth of high-integrity and high-uniformity single crystal materials due to its free of crucible contamination.However,the capillary convection in the melt brings a great challenge to the fl oating zone crystal growth.This is because the instability of convection will cause the formation of some crystal defects such as microscopic imperfections and macroscopic stripes.Therefore,it is very important to investigate the behaviors of the capillary fl ow and control its instability in order to improve the quality of the produced single crystal materials.In this paper,numerical simulations are performed to investigate all kind of the capillary convection in the half floating liquid bridge on the SixGe1?xsystem.And the impact of the external rotating magnetic fi eld is also investigated on the stability of capillary convection.The results show that the purely solute capillary convection is a two-dimensional axisymmetric model，and the temperature fi eld is mainly determined by thermal di ff usion while the concentration fi eld is dominated by convection and solute di ff usion together.On the other hand,the purely thermo-capillary convection presents three-dimensional unsteady axisymmetric fl ow.The concentration distribution is closely related to the fl ow direction of thermo-capillary convection.The isotherms bend in the region with strong convection.The coupled solute and thermocapillary convection is a three-dimensional periodic rotating oscillatory fl ow.When the rotating magnetic fi eld is applied,the circumferential velocity of the melt increases with increasing radius.Both the concentration fi eld and the fl ow fi eld in the melt show a two-dimensional axisymmetric distribution.

fl oating zone，capillary convection，surface tension，numerical simulation，rotating magnetic fi eld

O363.2,O782+.6

A doi：10.6052/0459-1879-17-102

2017–03–27 收稿，2017–09–26 錄用,2017–09–26 網(wǎng)絡版發(fā)表.

1)國家自然科學基金(51276089)、江蘇省“六大人才高峰”(2015-XNY-003)和安徽工業(yè)大學青年基金(QZ201517)資助項目.

2)黃護林，教授，博士，主要研究方向：磁流體流動與傳熱.E-mail:hlhuang@nuaa.edu.cn

鄒勇,朱桂平,李來,黃護林.液橋內熱質耦合對流不穩(wěn)定性及旋轉磁場法控制.力學學報,2017,49(6):1280-1289

Zou Yong,Zhu Guiping,Li Lai,Huang Hulin.Instability of coupled thermo-solute capillary convection in liquid bridge and control by rotating magnetic fi eld.Chinese Journal of Theoretical and Applied Mechanics,2017,49(6):1280-1289