張歡雷 程禮盛 王 安 譚 晶 楊衛(wèi)民 安 瑛

(北京化工大學(xué) 機電工程學(xué)院, 北京 100029)

引 言

在碳纖維加工過程中,主流的石墨體間接加熱技術(shù)存在石墨體在高溫下易揮發(fā)、裝置難以長時間持續(xù)運行等問題,并且能量利用效率低、能耗高,因此克服高溫限制,開發(fā)高效、高質(zhì)量、節(jié)能環(huán)保的石墨化技術(shù)是未來的發(fā)展趨勢[1]。 近年來,微波加熱技術(shù)由于具有加熱速度快、便于控制、清潔衛(wèi)生等優(yōu)點,逐漸開始應(yīng)用于碳纖維的碳化和石墨化過程。Carrott等[2]利用微波爐對碳纖維進行熱處理,可在120 s 內(nèi)將碳纖維加熱至1 100 ℃左右。 微波爐用于加熱碳纖維存在的問題是微波腔體內(nèi)場強分布不均,導(dǎo)致碳纖維在走絲過程中不能持續(xù)高效升溫。因此,設(shè)計一種專用于碳纖維的微波加熱裝置具有重要的現(xiàn)實意義。 一些國內(nèi)外學(xué)者針對圓柱形微波諧振腔腔體設(shè)計進行了相關(guān)研究。 李濤等[3]針對輪胎硫化進行了圓柱諧振腔的設(shè)計及數(shù)值模擬。 侯影飛等[4]利用COMSOL 多物理場仿真軟件,模擬了負載半徑、負載位置、微波饋入功率對微波加熱均勻性以及加熱效率的影響。 曹湘琪等[5-6]對圓柱形微波加熱器的效率和均勻性進行了研究,并利用HFSS軟件分析了圓柱形微波加熱器饋口位置和長度、內(nèi)筒半徑和高度以及負載厚度對微波吸熱效率的影響。 Cha-um 等[7]則利用數(shù)值模擬和實驗研究了矩形微波諧振腔中樣品位置、樣品尺寸和微波功率對加熱效率的影響。

目前工業(yè)中主要采用矩形諧振腔作為碳纖維石墨化的加熱腔[8]。 矩形諧振腔在單位功率微波源下場強較低,并且場強分布不均勻,很難滿足碳纖維石墨化的加工要求。 槽波導(dǎo)具有結(jié)構(gòu)簡單、尺寸公差要求低、損耗低、色散低、功率容量高等優(yōu)點[9],相較于矩形諧振腔,使用槽波導(dǎo)諧振腔作為微波石墨化的加熱腔具有明顯的優(yōu)勢。 因此,本文首次提出使用矩形槽波導(dǎo)諧振腔作為碳纖維石墨化專用諧振腔進行優(yōu)化設(shè)計,并研究了走絲速度與保護氣體流速兩個工藝參數(shù)對加熱碳纖維絲束的影響。

1 矩形槽波導(dǎo)諧振腔尺寸優(yōu)化設(shè)計

矩形槽波導(dǎo)諧振腔的橫截面如圖1 所示,其中橫截面的高a、兩端凸起高a′、底部寬b、兩端凸起寬d 的大小均會對諧振腔內(nèi)電場分布、截止頻率、歸一化截止波長產(chǎn)生影響。 當(dāng)矩形槽波導(dǎo)諧振腔橫截面的尺寸參數(shù)滿足a′/a≤0.6、a/b≥0.8、d/a≥0.6時,其截止頻率趨于一個穩(wěn)定值[10]。 根據(jù)實際使用要求,選擇各橫截面尺寸為a=120 mm,b=160 mm,a′=72 mm,d=60 mm。

借助Ansys Electronics 仿真軟件構(gòu)建微波場,選擇微波頻率為2.45 GHz,波導(dǎo)口為BJ-22。 以電場強度為標(biāo)的,0.1 mm 為步進,在100 ～400 mm 之間優(yōu)化腔體長度L。 如圖2 所示,最終選擇腔體長度L為219.20 mm。

2 數(shù)值模型與物理模型

2.1 模型簡化

為了驗證第1 節(jié)所設(shè)計的矩形槽波導(dǎo)諧振腔應(yīng)用于微波加熱碳纖維進行石墨化的有效性,利用Ansys Electronics 對微波加熱碳纖維的過程進行數(shù)值模擬。 為簡化計算過程和節(jié)省計算時間,對模型進行以下合理假設(shè):1)碳纖維絲束均勻且各向同性;2)碳纖維在加熱過程中介電性保持不變;3)碳纖維絲束和空氣的初始溫度均勻;4)空氣與碳纖維之間的化學(xué)反應(yīng)忽略不計;5)腔體邊界為完美邊界;6)端口由y 方向的橫向電磁駐波場激勵。

2.2 物理模型

選擇微波頻率為2.45 GHz,波導(dǎo)口為BJ-22,矩形槽波導(dǎo)長度L 為219.20 mm。 所選用的負載樣品為直徑1 mm 的碳纖維絲束,其位置如圖3 所示,走絲方向為x 軸正向。 絲束的相對介電常數(shù)為12-4j[11],導(dǎo)熱系數(shù)為10.465 W/(m·K),密度為1 700 kg/m3,恒壓熱容為795.5 J/(kg·K),絲束與腔體初始溫度均設(shè)為20 ℃。

2.3 控制方程

利用Ansys 仿真軟件構(gòu)建3 組雙向耦合。 首先建立微波-傳熱雙向耦合,采用Maxwell 方程計算碳纖維石墨化過程中的電磁傳播,電場波控制方程及耦合方式如下。

式中,μr為碳纖維相對磁導(dǎo)率;k0為自由空間的波數(shù);εr為碳纖維相對介電常數(shù);σ 為電導(dǎo)率;ω 為角頻率;ε0為真空介電常數(shù),8.85 ×10-12F/m;E 為電場強度,V/m。

微波場傳熱方程用傅里葉能量平衡公式來耦合,表達式為

式中,ρ 為碳纖維密度,kg/m3;cp為碳纖維常壓比熱容,J/(kg·K);T 為熱力學(xué)溫度,K;k 為碳纖維熱導(dǎo)率,W/(m·K);f 為微波頻率;δ 為碳纖維介電損耗角;Q 為體積熱源。

碳纖維絲束被電磁加熱之后,成為流場傳熱與輻射傳熱的熱源。 將氣體流場與傳熱雙向耦合,用以計算碳纖維絲束加熱過程中保護氣體所帶走的熱量。 將式(4)的Navier-Stokes 方程與式(5)的傳熱方程聯(lián)立,如下所示。

式中,F 為單位體積流體受的外力。

將流場中計算得出的速度u 與壓強p 導(dǎo)入傳熱方程中,可求得傳熱方程所需的溫度T。 溫度T 改變反向引起流體方程中與溫度有關(guān)的材料屬性的變化,繼而改變速度u 與壓強p。

最后構(gòu)建輻射-傳熱之間的耦合關(guān)系,用以計算加熱過程中絲束與腔體所發(fā)出的輻射引起的熱量傳遞。 使用能量守恒(廣義的傳熱)方程聯(lián)合求解有源邊界條件,如式(6)所示。

2.4 網(wǎng)格劃分

物理模型設(shè)置為兩個域,一個為氬氣域,另一個為絲束域。 全部網(wǎng)格均劃分為自由四面體網(wǎng)格,最大網(wǎng)格單元尺寸精細到微波波長的1/10。 絲束部分網(wǎng)格進行特別細化, 最大網(wǎng)格尺寸設(shè)為7.065 mm,最小尺寸設(shè)為0.211 9 mm。 網(wǎng)格數(shù)量為40 萬。

首先進行網(wǎng)格獨立性驗證。 針對網(wǎng)格劃分,分別進行網(wǎng)格數(shù)量為20 萬、30 萬、40 萬的3 次計算,以絲束所能達到的最高溫度與最高溫度所在位置作為標(biāo)準,結(jié)果如表1 所示,可以看出三者之間的溫度偏差較小,所在位置變化也很小,因此認為網(wǎng)格數(shù)量符合精度要求。 下文采用40 萬網(wǎng)格數(shù)量進行仿真模擬。

3 仿真結(jié)果

3.1 腔體內(nèi)電場與溫度分布

對腔體中間切面進行電場強度與溫度分析,仿真結(jié)果如圖4 和圖5 所示。 圖4 中有3 塊場強較強且均勻的區(qū)域,絲束處于中間區(qū)域。 仿真結(jié)果顯示,其輸入反射系數(shù)S11為-17.4,證明此時微波加熱效率 較 高。 根 據(jù) 圖 5 結(jié) 果, 絲 束 最 高 溫 度 為1 575.40 ℃,且溫度較高區(qū)域較為集中,證明其升溫效果較均勻。

3.2 保護氣體流速對絲束溫度的影響

實際生產(chǎn)過程中,保護氣體流動具有帶走碳纖維石墨化過程中產(chǎn)生的反應(yīng)廢氣、保持反應(yīng)腔內(nèi)正壓以及降低腔體溫度等重要作用,因此保護氣體需要有一定的流速。 為了分析碳纖維絲束在不同氣體流速下的升溫特性,固定走絲速度和微波功率分別為1 mm/s 和700 W,在0.01 ～0.10 m/s 范圍內(nèi)改變保護氣流速進行穩(wěn)態(tài)模擬,結(jié)果如圖6 所示。 圖6數(shù)據(jù)表明,絲束所能達到的最高溫度隨著保護氣流速的增加呈非線性降低。 氣流速度對最高溫度出現(xiàn)的軸向位置(最高溫度x 值)也有一定影響,隨著氣流速度的不斷增大,最高溫度點到絲束入口位置的距離不斷增大。 從圖中結(jié)果來看,兩條曲線隨氣體流速的變化均較為明顯:氣流速度從0.01 m/s 增大到0.10 m/s,最高溫度下降了近200 ℃,最高溫度所在位置移動了12 mm 左右。 這可能是受到對流換熱的影響,即隨著氣體流速的增大,單位時間內(nèi)氣體可以帶走更多的熱量。

3.3 走絲速度對絲束溫度的影響

為了分析碳纖維絲束在不同走絲速度下的升溫特性,固定保護氣體流速和微波功率分別為0.1 m/s和700 W,在0.5 ～3.0 mm/s 范圍內(nèi)改變走絲速度,進行穩(wěn)態(tài)模擬,結(jié)果如圖7 所示。

從圖7 可以看出,絲束所能達到的最高溫度隨著走絲速度的增大呈非線性降低;最高溫度所在位置到絲束入口的距離隨著走絲速度的增大幾乎呈線性增大。 總的來說,走絲速度從0.5 mm/s 增大到3.0 mm/s,最高溫度下降20 ℃左右,變化幅度不是特別大。 最高溫度x 值變化有些波動,這可能是由于絲束在運動過程中會引起氣流變化。 實際生產(chǎn)過程中,更高的走絲速度意味著更高的加工效率,因此在纖維于微波場中停留時間有保證的情況下,可以盡量提高走絲速度。

4 結(jié)論

(1)本文所優(yōu)化設(shè)計的矩形槽波導(dǎo)諧振腔腔體電場場強分布較均勻,輸入反射系數(shù)S11可以達到-17.4,證明加熱效率較高;通過溫度分布圖可以得出,絲束升溫區(qū)域較為集中,絲束可以穩(wěn)定升溫。 以上結(jié)果驗證了本文腔體結(jié)構(gòu)設(shè)計的合理性。

(2)在使用微波反應(yīng)器加熱碳纖維過程中,碳纖維溫度分布與保護氣體流速、走絲速度緊密相關(guān)。隨著流體速度的增大,絲束能達到的最高溫度降低較快。 而在實際生產(chǎn)過程中,由于保護氣體需要一定流速才能起到帶走反應(yīng)廢氣的作用,所以在實際應(yīng)用中應(yīng)注意保護氣體流速的選擇。

(3)走絲速度從0.5 mm/s 增大到3.0 mm/s,最高溫度下降20 ℃左右,證明絲束能達到的最高溫度隨走絲速度的增大略有降低,但幅度不大,因此在絲束于腔體中停留時間有保證的情況下,可以大幅提高走絲速度。

(4)保護氣流速與走絲速度均導(dǎo)致絲束最高溫度所在位置發(fā)生變化,因此在石墨化過程中,絲束測溫點位置應(yīng)根據(jù)保護氣流速與走絲速度進行綜合計算得到。