李 陽(yáng)，陳 杰，王旭陽(yáng)，王 飛

1.陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西 渭南 714000；2.西安科技大學(xué)材料學(xué)院，陜西 西安 710054；3.西北工業(yè)大學(xué)凝固技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710072

碳熱還原合成碳化硅過(guò)程的數(shù)值模擬

李 陽(yáng)1，陳 杰2，王旭陽(yáng)3，王 飛1

1.陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西 渭南 714000；2.西安科技大學(xué)材料學(xué)院，陜西 西安 710054；3.西北工業(yè)大學(xué)凝固技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710072

為揭示碳化硅合成過(guò)程中能量及物質(zhì)擴(kuò)散機(jī)理，從而為碳化硅的提質(zhì)增產(chǎn)奠定理論基礎(chǔ)，采用數(shù)值模擬的方法對(duì)碳化硅合成過(guò)程中的溫度場(chǎng)、壓力場(chǎng)、氣體流動(dòng)規(guī)律進(jìn)行模擬研究。結(jié)果表明，隨著合成時(shí)間的延長(zhǎng)，爐內(nèi)熱量呈輻射狀向外擴(kuò)散，合成爐內(nèi)氣體呈現(xiàn)三維多向流動(dòng)特性，反應(yīng)進(jìn)行到24 h時(shí)CO氣體流量達(dá)到最大，而此時(shí)由于爐底透氣性差的原因，致使?fàn)t底部壓力高于其余位置，最大可達(dá)1.525×101 kPa，此時(shí)可加入少量木屑以增加爐底透氣性來(lái)改善因壓力過(guò)高所造成的噴爐事故。模擬結(jié)果得到了生產(chǎn)實(shí)踐驗(yàn)證。

碳熱還原 碳化硅 冶煉爐 數(shù)值模擬

碳化硅(SiC)具有高密度、高純度、高結(jié)晶性、高均勻性等優(yōu)越性能，目前依然沿用Acheson發(fā)明的爐型進(jìn)行生產(chǎn)，所采用的原理為碳熱還原法[1-4]。國(guó)內(nèi)外科學(xué)家一直致力于優(yōu)化碳化硅冶煉合成的方法，Mccolm[5]設(shè)計(jì)出一種新式β-SiC合成設(shè)備，該裝置中下部的SiO和Si在高溫下形成SiO2氣體，上升同碳粉層反應(yīng)生成SiC。王曉剛等[6]設(shè)計(jì)提出了多熱源的碳化硅合成設(shè)備，并投入生產(chǎn)。但是目前研究多圍繞工藝改進(jìn)及爐體的優(yōu)化設(shè)計(jì)，長(zhǎng)期以來(lái)由于合成爐內(nèi)溫度過(guò)高、壓力及氣體流動(dòng)變化較快，實(shí)驗(yàn)數(shù)據(jù)相對(duì)較少，這給開展碳化硅傳熱傳質(zhì)機(jī)理研究帶來(lái)很大困難，本工作利用計(jì)算機(jī)模擬的方法，模擬得到生產(chǎn)過(guò)程中溫度、壓力、氣體流動(dòng)速度變化規(guī)律，通過(guò)結(jié)合工業(yè)試驗(yàn)數(shù)據(jù)，揭示生產(chǎn)過(guò)程中能量及物質(zhì)傳遞的規(guī)律，以期為提高SiC成品質(zhì)量，改進(jìn)合成方案提供理論基礎(chǔ)。

1 實(shí)驗(yàn)部分

1.1 實(shí)驗(yàn)原料

作為硅質(zhì)原料的石英砂和碳質(zhì)原料的無(wú)煙煤均來(lái)自寧夏石嘴山，石英砂經(jīng)破碎過(guò)篩后使其粒度小于2 mm，純度為98.5%(w)，無(wú)煙煤經(jīng)破碎過(guò)篩后使其粒度小于1 mm。

1.2 實(shí)驗(yàn)設(shè)備及方法

碳化硅合成爐為自制設(shè)備，如圖1所示。爐墻一側(cè)開有石英玻璃窗口，窗口至爐芯間埋有中空的石墨碳管，碳管直徑3 cm。為避免爐芯氧化和煙氣干擾，在窗口上方開有通氣口，Ar氣流量為0.5 L/min。根據(jù)反應(yīng)不同階段給合成爐提供不同的供電功率，0～4 h為0.88×106W/m3、4～32 h為0.83×106W/m3、32～48 h為0.75×106W/m3。合成開始前將預(yù)先混合好的配合料裝入合成爐中，爐體頂部覆蓋碳粉，供電過(guò)程中通過(guò)紅外測(cè)溫儀對(duì)爐芯至爐壁測(cè)溫段進(jìn)行溫度采集，用于模擬結(jié)果的驗(yàn)證。

圖1 實(shí)驗(yàn)測(cè)溫裝置Fig.1 Schematic diagram of temperature measurement device

圖2 窯爐三維幾何模型Fig.2 Three-dimensional model of furnace

2 數(shù)值模擬

2.1 幾何模型

圖2為單熱源爐計(jì)算的幾何模型，其結(jié)構(gòu)尺寸采用工業(yè)實(shí)際尺寸，具體參數(shù)見表1。

表1 窯爐結(jié)構(gòu)尺寸Table 1 Structure of glass furnace size

2.2 簡(jiǎn)化假設(shè)及邊界條件設(shè)置

(1) 假定在整個(gè)合成過(guò)程中，爐內(nèi)物料混合充分均勻，并且比較密實(shí)，反應(yīng)料是各向同性的。

(2) 將合成爐內(nèi)的反應(yīng)物料等效為一種物質(zhì)，通過(guò)數(shù)學(xué)擬合方式得到其比熱容、導(dǎo)熱系數(shù)、密度隨時(shí)間變化的函數(shù)[7,8]。

(3) 激活能量方程選擇各向同性熱傳導(dǎo)模型，輻射模型選擇DO模型[9]。對(duì)熱源的四個(gè)邊界類型設(shè)定為HOTWALL（熱邊界），單熱源熱流密度加載12 W/cm2。

(4) 選擇完全預(yù)混的反應(yīng)模型，加載化學(xué)反應(yīng)，其中具體化學(xué)反應(yīng)方程式機(jī)理如表2所示[10]。

(5) 選擇爐內(nèi)反應(yīng)物為多孔介質(zhì)，并設(shè)定其氣孔率為30%，初始設(shè)定質(zhì)量比SiO2:C為4:3。

表2 碳熱還原合成碳化硅反應(yīng)機(jī)理Table 2 Structure of glass furnace size

2.3 數(shù)學(xué)模型

根據(jù)上述的假設(shè)，可以得到如下的模型方程[11]。

連續(xù)性方程：

動(dòng)量守恒方程：

能量守恒方程：

其中，ρ為密度，kg/m3；t為時(shí)間，s；u、v、w分別為x、y和z方向上的速度分量，m/s；Fx、Fy和Fz分別是單位質(zhì)量上的質(zhì)量力在三個(gè)方向上的分量；P為壓力，N/m2；μ為粘度，Pa·s；cp為比熱容，J/(kg·K)；λ為導(dǎo)熱系數(shù)，W/(m·K)；qv為單位時(shí)間、單位體積的化學(xué)反應(yīng)熱，W/m3。

2.4 網(wǎng)格劃分

圖3 冶煉爐的網(wǎng)格劃分Fig.3 The grid of furnace

在ICEM-CFD中選用以六面體為主的混合網(wǎng)格對(duì)合成爐主體結(jié)構(gòu)進(jìn)行網(wǎng)格劃分，網(wǎng)格數(shù)量為1 929 752，網(wǎng)格劃分的三維模型如圖3所示。

3 結(jié)果與討論

3.1 溫度場(chǎng)分析

圖4 合成爐內(nèi)不同時(shí)間的演變規(guī)律Fig.4 The evolution under different synthesis times in the furnace

圖8 不同時(shí)刻爐內(nèi)壓力分布Fig.8 Pressure distribution of furnace under different times

圖4為合成爐內(nèi)溫度分布隨時(shí)間的演變趨勢(shì)圖，從圖中可以看出以熱源為中心熱量均勻向外傳遞，形成了一個(gè)近似圓形的均勻溫度場(chǎng)。當(dāng)合成時(shí)間12 h時(shí)，熱源溫度已經(jīng)達(dá)到1 800 ℃，此時(shí)靠近熱源處的爐料已經(jīng)發(fā)生反應(yīng)生成SiC，即發(fā)生化學(xué)反應(yīng)。從圖中可以看出實(shí)驗(yàn)裝置24 h爐體頂部以發(fā)生小規(guī)模噴爐。36 h此時(shí)已經(jīng)行到反應(yīng)后期爐內(nèi)壓力明顯下降。圖9為距離熱源X為0.1，0.4，0.6 m處停爐時(shí)壓力分布的散點(diǎn)圖，從圖中可以看出，在窯爐豎直方向壓力呈現(xiàn)明顯的梯度變化，爐底至熱源底部壓力下降幅度相對(duì)較小，而在熱源及周圍附近由于反應(yīng)劇烈且離頂部出口越來(lái)越近，此范圍內(nèi)壓力下降梯度較為明顯，在靠近爐頂附近由于直接與周圍大氣相連，壓力變化較小基本與周圍大氣壓一致。

圖9 距離熱源不同位置Y方向壓力分布散點(diǎn)Fig.9 Scatter diagram of pressure distribution under different distance of heat source along Y direction

3.3 氣體流動(dòng)特性分析

圖10為不同時(shí)刻氣體流動(dòng)軌跡圖，從圖中可以看出隨著溫度的升高壓力的增大，大量氣體由熱源部向周圍呈放射狀擴(kuò)散。反應(yīng)進(jìn)行到在20 h時(shí)，由于外層的爐料也達(dá)到反應(yīng)溫度1 500 ℃，反應(yīng)生成的氣體（主要為CO）一部分向爐外流動(dòng)，一部分向著熱源方向流動(dòng)，此時(shí)與熱源處反應(yīng)生成的氣相物質(zhì)交匯產(chǎn)生一個(gè)氣相回流圈。圖11為SiC結(jié)晶筒橫截面及頂部放大圖可以看出由于氣相回流圈的存在，碳化硅結(jié)晶筒中形成一條明顯的分界線，在分界線內(nèi)圈SiC明顯較為致密，氣孔相對(duì)比較小且分布均勻，而在分界區(qū)域上部氣孔大小不一且分布雜亂。反應(yīng)后期熱源附近分解的碳化硅在結(jié)晶筒中產(chǎn)生較大的孔洞。

圖10 不同時(shí)刻氣體流動(dòng)軌跡Fig.10 Gas flow trajectory on the cross section of furnace at different times

圖11 SiC結(jié)晶筒橫截面及頂部放大圖譜Fig.11 Cross section of SiC crystal tube and the top of enlargement pattern

圖12為三維氣體流動(dòng)軌跡圖及工業(yè)實(shí)際圖，可以看出反應(yīng)生成氣體由熱源處向四周擴(kuò)散，在爐壁遇到阻礙后聚集在爐壁四周，聚集著的氣體通過(guò)保溫墻的縫隙向外部擴(kuò)散，在工業(yè)生產(chǎn)的爐體也可以看到，爐體四周傾斜著滲出的火苗，火焰顏色基本為亮黃色和淡藍(lán)色，說(shuō)明滲出的氣體中基本為CO氣體和含氫類物質(zhì)。模擬的結(jié)果和工業(yè)實(shí)際相符合。

圖12 三維氣體流動(dòng)軌跡圖及工業(yè)實(shí)際圖譜Fig.12 Three-dimensional trajectory diagram of gas flow and industrial production scene diagram

4 結(jié) 論

a）合成過(guò)程中熱量由熱源中心部位成放射狀向四周擴(kuò)散，隨著合成時(shí)間的延長(zhǎng)碳化硅的結(jié)晶筒逐漸長(zhǎng)大，由于合成爐內(nèi)溫度梯度的存在，使得結(jié)晶筒不同位置成品含量、致密度不盡相同；36 h爐芯表面溫度已高達(dá)2 600 ℃，此時(shí)爐芯周圍SiC發(fā)生分解，適當(dāng)延長(zhǎng)反應(yīng)時(shí)間促使少量SiC分解可以提高一級(jí)品的產(chǎn)量。

b）由于爐底透氣性差，產(chǎn)生的氣體難以擴(kuò)散，造成此處壓力局部過(guò)高，最高處可達(dá)1.525×101 kPa，生產(chǎn)時(shí)極易發(fā)生噴爐事故；生產(chǎn)時(shí)可在爐底底增加透氣材料，爐壁四周加裝氣體收集裝置，從而可以有效釋放爐底壓力、收集有害氣體。

c）合成爐內(nèi)在反應(yīng)中期會(huì)產(chǎn)生明顯的氣相回流圈，在分界線內(nèi)圈SiC明顯較為致密，氣孔相對(duì)比較小且分布均勻。

[1] Casady Ssdy J B, Johnson R W. Status of silicon carbide as a wide-band gap semiconductor for high-temperature applications: a review[J]. Solid State Electron, 1996, 39(10):1409-1422.

[2] ‘Katoh Y, Snead L L, Henager C H, et al. Current status and recent research achievements in SiC/SiC composites [J]. Journal of Nuclear Materials, 2014, 455(1):387-397.

[3] 趙 佶. 全球半導(dǎo)體產(chǎn)業(yè)掀起研發(fā)碳化硅芯片潮流 [J]. 半導(dǎo)體信息, 2014, (4):36-41.Zhao Ji. Semiconductor industry inspired a trend of study on silicon carbide chip [J]. Semiconductor Information, 2014, (4):36-41.

[4] Kong Y, Shen X D, Cui S. Effect of carbothermal reduction temperature on microstructure of fiber reinforced silicon carbide porous monoliths with high thermal resistance [J]. Chinese Journal of Inorganic Chemistry, 2014, 30(12):2825-2831.

[5] Macolm I J. Informing, shaping, and working of high-performance ceramics [M]. New York: Chapman and Hall, 1988:31-40.

[6] 王曉剛, 李曉池. 冶煉碳化硅的多爐芯爐及其生產(chǎn)碳化硅的方法:中國(guó). 1364996 [P]. 2002-08-21.

[7] 王曉剛. 碳化硅合成理論與技術(shù) [M]. 陜西科學(xué)技術(shù)出版社, 2001:12-19.

[8] 楊振明,田沖,矯義來(lái)等.泡沫碳化硅的制備及應(yīng)用[J]化學(xué)反應(yīng)工程與工藝,2013,29(3):269-275. Yang Zhen-ming,Tian Chong,Jiao Yi-lai etal.Preparation and Applications of Foam SiC. Chemical Reaction Engineering and Technology,2013,29(3):269-275.

[9] Chen J, Wang X G, Li Y. Diffusion mechanism of energy flow in multi-heat-source synthesis of SiC [J]. Journal of Wuhan University of Technology(Materials Science Edition), 2015, 2:266-270.

[10] 陳 杰. 多熱源多向流體系碳化硅材料制備理論及應(yīng)用 [D]. 西安:西安科技大學(xué), 2011.

[11] John D, Anderson Jr. 計(jì)算流體力學(xué)入門 [M]. 清華大學(xué)出版社, 2002:39-42.

Numerical Simulation of Synthesis Silicon Carbide by Carbothermal Reduction

Li Yang1, Chen Jie2, Wang Xuyang3, Wang Fei1
1. Shaanxi Railway Institute, Weinan 714000, China; 2. School of Material Science and Engineering, Xi’an University of Science and Technology, Xi’an 710054, China; 3. State Key Laboratory of Solidification Processing, Northwestern Polytechnical University, Xi’an 710072, China

In order to reveal the mechanism of energy and material diffusion in the process of silicon carbide synthesis so as to lay the theoretical foundation to improve the yield and quality of silicon carbide, the numerical simulation method was used to simulate the temperature field, pressure field and gas flow in the process of silicon carbide synthesis. The results showed that with the increase of the synthesis time, the heat in the furnace was radiated outward and the gas in the synthesis furnace exhibited three-dimensional multi-directional flow characteristic. When the reaction progressed to 24 h, the gas flow rate reached the maximum. At this moment, due to the poor ventilation at the bottom of the furnace, the pressure at this point was higher than the remaining position and the maximum was 1.525×101 kPa. At this time, a small amount of wood chips can be added to increase the permeability of the furnace to improve the pressure to avoid the spouting of furnace by the over high pressure. The simulation results have been verified by the production practice.

carbothermal reduction; silicon carbide; smelting furnace; numerical simulation

TQ173

A

1001—7631 ( 2017 ) 01—0090—07

10.11730/j.issn.1001-7631.2017.01.0090.07

2016-12-07;

2017-01-13。

李 陽(yáng)（1990—），男，助教。E-mail:554636457@qq.com。

陜西省教育廳科研計(jì)劃項(xiàng)目（12JK0785）；陜西鐵路工程職業(yè)技術(shù)學(xué)院2016第二批科研基金項(xiàng)目（KY2016-34）