炭/炭復合材料熱解炭基體微觀結構轉變動力學

黃清波,　張　丹，　白瑞成，　李愛軍，　孫晉良

(1.上海大學 材料復合及先進分散技術教育部工程中心， 上海200072)(2.上海大學 上海市應用數(shù)學及力學研究所， 上海200072)

?

炭/炭復合材料熱解炭基體微觀結構轉變動力學

黃清波1,張丹2，白瑞成1，李愛軍1，孫晉良1

(1.上海大學 材料復合及先進分散技術教育部工程中心， 上海200072)(2.上海大學 上海市應用數(shù)學及力學研究所， 上海200072)

摘要：對熱解炭沉積和織構形成過程進行動力學建模，重點分析C/C復合材料CVI制備工藝中基體炭形成時中織構/高織構(MT/HT)熱解炭之間輪廓分明的急劇轉變現(xiàn)象?；贚angmuir-Hinshelwood(L-H)理論和Particle-filler(P-F)概念模型，將MT和HT熱解炭作為炭的兩種亞穩(wěn)相，以氣相中占優(yōu)的兩種中間組分作為基體炭前驅體(線性小分子烴和小分子芳香烴)，考慮基體表面單分子沉積形成MT熱解炭和P-F雙分子反應形成HT熱解炭的過程，建立包含吸附/解吸附/脫氫的多步非均相熱解炭沉積和織構形成反應動力學模型，研究該動力學系統(tǒng)達到穩(wěn)態(tài)時熱解炭隨氣相組成變化的情況。結果表明，熱解炭沉積和織構形成過程曲線呈現(xiàn)“S”型特征，該曲線的線性穩(wěn)定性分析表明熱解炭沉積中的織構轉變是一個包含遲滯區(qū)間的雙穩(wěn)態(tài)過程，進一步的計算表明此遲滯區(qū)間的大小明顯受初始直鏈烴濃度以及沉積溫度的影響。

關鍵詞：炭/炭復合材料； 熱解炭； 微觀結構； 動力學

1前言

炭/炭(C/C)復合材料是一種高性能高溫結構材料，具有高比強度/比模量、良好的韌性，在高溫狀態(tài)下具有優(yōu)良的強度保持率，以及優(yōu)良的摩擦性能和較高的導電、導熱系數(shù)，廣泛應用于航空、航天等領域。按照制備工藝不同，炭/炭復合材料的基體炭可分為：瀝青炭、樹脂炭和氣相熱解炭，以后者為基體的炭/炭復合材料具有更好的綜合性能。等溫等壓化學氣相滲透(I-CVI)工藝是制備熱解炭基體最常用的工藝，其特征是采用氣相烴類為前驅氣體，經(jīng)復雜的氣相反應和表面沉積在炭纖維編制體內部得到熱解炭基體。根據(jù)工藝條件不同，烴類氣體的熱解可以生成具有不同織構的熱解炭基體，而不同的熱解炭織構直接影響到材料宏觀力學性能，盡管采用高分辨透射電子顯微鏡(HR-TEM)觀察發(fā)現(xiàn)高織構熱解炭(HT)內織構仍存在漸變的現(xiàn)象，但在偏光顯微鏡下觀察經(jīng)常發(fā)現(xiàn)不同織構間的變化是陡變的，其中最顯著的是中織構(MT)和高織構(HT)熱解炭之間的突變現(xiàn)象[1,2]。如圖1所示，相比子圖1(a，b)的單一織構占優(yōu)情況，圖1(c)中MT與HT熱解炭分別主導沉積前期和后期，且相互之間涇渭分明，并沒有顯著的過渡區(qū)。因此，探索C/C復合材料中MT/HT熱解炭突變動力學是深入理解熱解炭的表面沉積和織構形成過程的鑰匙，有助于對C/C復合材料熱解炭微觀結構的預測，并為最終實現(xiàn)熱解炭微觀結構的可控制備提供理論支持。

圖 1　偏光顯微鏡下不同熱解炭微觀結構的光學活性:(a)MT，(b)HT，(c)MT+HT

Benzinger和Hüttinger[3-5]證明了表面沉積反應是非均相反應，并與均相(熱解)反應相互競爭。其中非均相反應有兩種歷程[6,7]：生長歷程和形核歷程。生長歷程依靠基層表面的活性點化學吸附氣相中以線型小分子和小分子芳香烴為主的活性組分，并在石墨烯層邊緣生長，符合經(jīng)典的Langmuir-Hinshelwood動力學方程。形核歷程依靠基層表面物理吸附以多環(huán)芳烴(PAH)等大分子為主的組分形成新的石墨烯層，因此吸附平衡很重要，但不需要活性點的存在。Dong等[8]采用ICVI，沉積溫度為1 100 ℃時以甲烷為原料氣在較大壓力范圍內研究了熱解炭微觀結構的形成機理，提出了在生長階段熱解炭織構生成的Particle-filler模型。他們認為生成高織構熱解炭的必要條件是氣相中存在比例最優(yōu)化的芳香烴C6(Particle)與直鏈烴是C2(filler)，氣相中存在過多的芳香化合物，或者過多的直鏈烴，都將導致中低織構熱解炭的生成。對CVI過程來說，由于炭纖維預制體的比表面積AS/VR大，難以實現(xiàn)飽和吸附，氣相組分大量消耗在反應表面大的PAH分子難以在孔內形成，因此熱解炭沉積符合生長歷程[7]。Vignoles和Delhaes等[9,10]采用ICVI和丙烷為前驅氣體，認為以小分子線性烴類易沉積生成中織構熱解炭而小分子芳香烴類易生成高織構熱解炭，并且通過動力學分析，初步表明熱解炭織構轉變是屬于雙穩(wěn)態(tài)過程，兩個穩(wěn)態(tài)之間存在遲滯轉變區(qū)間。

基于Langmuir-Hinshelwood(L-H)理論，以氣相組分中占主要的小分子線性烴類C2(Filler)和小分子芳香烴類C6(Particle)作為基體炭前驅體，以MT和HT熱解炭作為兩種不同亞穩(wěn)相，假設在基體表面單分子前驅體沉積形成MT熱解炭而P-F雙分子反應形成HT熱解炭，即以Particle-filler概念模型為指導對熱解炭沉積過程進行動力學分析，研究氣相組分濃度比、初始C2濃度以及沉積溫度對熱解炭沉積和織構形成過程的影響。

2模型及參數(shù)設定

2.1熱解炭沉積和織構形成非均相表面反應動力學模型

考慮小分子線性烴類C2和小分子芳香烴類C6這兩種主要的氣相組分，將沉積成炭的過程按照L-H理論分為組分吸附、表面組分解吸附和表面組分脫氫3個過程。同時，基于P-F概念模型將MT和HT熱解炭作為炭的兩種亞穩(wěn)相，以氣相中占優(yōu)的兩種中間組分作為基體炭前驅體(線性小分子烴和小分子芳香烴)，考慮基體表面單分子沉積形成MT熱解炭和表面P-F雙分子反應形成HT熱解炭的過程。如圖2所示，反應機理中總共包含7種組分和9個動力學常數(shù)，不考慮初始氣相烴類原料氣對表面反應的影響，即假設入口原料氣在進入反應室并到達沉積表面時已經(jīng)完全轉化為線性小分子烴類C2和芳香小分子烴類C6，此時兩者的比值R為0。

圖 2　熱解炭沉積和織構形成非均相表面反應動力學模型

探索穩(wěn)態(tài)時基體近表面處不同的穩(wěn)態(tài)氣相組成對熱解炭微觀結構形成動力學的影響，因此不同氣相組分間轉化氣相反應動力學不在討論范圍內，即k0=0，k8=0。同時，根據(jù)HACA機制乙炔在苯環(huán)上的加成反應可以得到完善的芳香烴結構，但在表面反應中整個域內單個或數(shù)個此類反應并不能決定熱解炭的織構(圖3)，當表面上某兩相鄰域的邊界分別存Aarmchair以及Zig-zag位就可能在兩者之間會形成五元環(huán)以及七元環(huán)從而引起整個炭平面的上翹影響熱解炭的織構。由此可知，整個域內的熱解炭織構取決于整體的加成情況，即根據(jù)P-F概念模型只有F與P達到合適比例時，F(xiàn)才能完美填充P間的邊緣形成完美的炭層，否則就會形成圖3的非六元環(huán)。因此，僅以氣相中占優(yōu)的兩種中間組分作為基體炭前驅體(線性小分子烴和小分子芳香烴)研究其在表面沉積過程中的宏觀動力學特征。

基于反應機理(圖2)，具體得到如下多相化學反應：

(r1)

(r2)

(r3)

(r4)

(r5)

圖 3　表面上兩個基本結構單元之間

式中，S代沉積表面的活性位，C2H2*、C6H6*分別表示表面吸附的組分，主要是將P-F概念模型進行定量化，從而驗證該模型的正確性。在該模型中將氣相中所有的線性小分子烴以C2表示作為filler，將所有的大分子芳香烴以C6表示作為Particle，而對于各氣相組分在基體表面發(fā)生吸附時占有的活性位目前還沒有統(tǒng)一的定論，因此僅以含炭量來加以衡量，假設每個C2H2占有一個表面活性位，則每個C6H6占3個活性位，r1-r5均為基元反應，由此可建立表面吸附物濃度隨時間變化的動力學方程為：

(e1)

(e2)

(e3)

當生成的熱解炭織構和熱解炭沉積速率不再隨時間變化時，系統(tǒng)既達到穩(wěn)態(tài)，此時表面[C2H2*]和[C6H6*]保持恒定，即式(e1)和(e2)滿足下式：

(e4)

將上式e4代入(e1～e3)得：

(e5)

(e6)

(e7)

2.2非均相表面反應動力學模型設置

將表面組分覆蓋率θ作為變量引入方程(e5)～(e7)，并以C2H2的初始濃度[C2H2]0對系統(tǒng)進行無量綱化處理：

(e8)

(e9)

(e10)

式中：

a=[C2H2]/[C2H2]0,b=[C6H6]/[C2H2]0

(e11)

式中，反應速率常數(shù)ki與反應溫度T的關系滿足Arrhenius關系式：

(e12)

式中，Ai為指前因子，R為理想氣體常數(shù)，ΔEi為反應的激活能。

文獻[9]給出了沉積溫度T為1 373 K時各反應速率常數(shù)ki，由(e12)可知只要確定相關反應的激活能就能計算出相指前因子等相關參數(shù)的值。取文獻[11]的反應激活能作為相關反應的激活能值，見表1, 2。

表 1　文獻中的相關表面反應激活能[11]

Note： C.(S)—active site; C(B)—deposited carbon; C6H6#—benzene.

3模擬結果與討論

3.1氣相組成對熱解炭沉積的影響

Vignoles的熱解炭沉積模型和Hüttinger的Particle-filler概念模型均認同氣相組成對熱解炭織構形成過程具有決定性影響。此處以氣相中[C6]和[C2]比值R作為表征氣相組分成熟度的主要變量，假設沉積表面處線性小分子烴C2歸一化濃度a為1，沉積溫度T為1 373 K，計算了隨R值的變化(即隨小分子芳香化合物C6歸一化濃度的變化)熱解炭沉積系統(tǒng)達到穩(wěn)態(tài)時表面組分和熱解炭織構的演變規(guī)律。求解(e8～e10)方程組在不同R值下的正實數(shù)解，如圖4和5分別顯示了C2與C6組分的表面覆蓋率以及MT與HT熱解炭的沉積速率隨R的變化曲線，兩組曲線都一致呈現(xiàn)“S”型，存在著上下兩條實線分支和中間的虛線分支，為確定該“S”型曲線的三分支的狀態(tài)，需要進行系統(tǒng)穩(wěn)定性分析。所謂系統(tǒng)的穩(wěn)定性，就是指給一個已經(jīng)達到平衡狀態(tài)的系統(tǒng)施加一個隨機擾動，在擾動作用消失后經(jīng)過一段過渡過程后能否回復到原來的平衡狀態(tài)或足夠準確地回復到原來的平衡狀態(tài)的性能。若系統(tǒng)能恢復到原來的平衡狀態(tài)，則稱系統(tǒng)是穩(wěn)定的；若擾動消失后系統(tǒng)無法恢復到原來的平衡狀態(tài)，而是偏差越來越大，則系統(tǒng)不穩(wěn)定[12]。

表 2　熱解炭沉積和織構形成

Note：*activation energy from Table 1；#activation energy from gas-phase reaction[12]： A3+H=A3-4+H2and A3-4+C2H2=A3C2H2-4.

通過對系統(tǒng)進行線性穩(wěn)定性分析[13]即可確定各分支的穩(wěn)定性。將e3帶入e1和e2，并進行無量綱化處理得到：

?tx=k1a[C2H2]0(1-x-3y)-(k2+k5)x-k7[S]0xy

-k7[S]0xy(e13)

取x=x0+Δx,y=y0+Δy代入方程(e13)中，即分別給曲線上的某個x0和y0分別加上擾動Δx和Δy得到一個線性化擾動方程(e14)：

(e14)

圖 4　C2和C6表面覆蓋率隨R值的變化

式中：

a1=-(k1a[C2H2]0+k2+k5+k7[S]0y0

a3=-3k1a[C2H2]0+k7[S]0x0

此式的特征方程為：

γ2+c1γ+c2=0

(e15)

式中：

c1=k2+k4+k5+k6+k1a[C2H2]0+

按照Routh-Hurwitz判據(jù)[14,15]，定態(tài)解x0和y0為穩(wěn)定解的充要條件是：

(c1>0)∩(c2>0)

(e16)

若上式中有任一條件不滿足，則系統(tǒng)的定態(tài)解失穩(wěn)。顯然c1中各項都大于零，所以針對圖中各段只需計算c2是否滿足(e16)即可。分別將圖4和5的上中下三個分支曲線上點對應的(R,x0,y0)值代入c2表達式，易驗證曲線的上下兩分支滿足(e16)，而中間虛線分支不滿足，因此,動力學上，上下兩分支穩(wěn)定而中間虛線分支是不穩(wěn)定的，實驗中該分支不存在，證明熱解炭沉積過程是一個存在遲滯回線的雙穩(wěn)態(tài)過程。

圖 5　HT和MT隨R值的(a)沉積速率及其(b)沉積速率比

可見，熱解炭沉積和織構形成過程是個含有遲滯回線的雙穩(wěn)態(tài)的動力學問題。當初始R較小時，隨著R的增加C2的表面覆蓋率沿著上分支緩慢減小，但此時表面主要還是由吸附的C2占據(jù)，結合圖5可知，HT熱解炭相對于MT熱解炭沉積速率優(yōu)勢很小，可判定此階段將導致整體的熱解炭織構呈現(xiàn)MT織構特征；當R持續(xù)增加到R2時，C2的表面覆蓋率將突變至下分支，之后沿著下分支變化，此時表面主要由大分子C6占據(jù)，主要沉積HT熱解炭且沉積速率較大，可判定此階段形成的熱解炭將呈現(xiàn)HT織構的特征。當R值從較大值減小至R1時，HT沉積機制將會突變?yōu)镸T沉積機制。在這過程中存在著如圖上所示的回形遲滯區(qū)間，對于遲滯現(xiàn)象的產生，Vignoles認為,由基體表面的雙分子反應引起的，而氣相擴散與表面多相反應的相互作用則改變產生遲滯現(xiàn)象的參數(shù)區(qū)間[9]。

3.2線性小分子烴對熱解炭沉積的影響

為分析線性小分子烴對熱解炭沉積的影響，研究不同的線性小分子烴濃度條件下系統(tǒng)達到穩(wěn)態(tài)時R值的變化對于熱解炭沉積的影響。設定沉積溫度T為1 373 K，線性小分子烴C2的歸一化濃度a分別為0.1、0.5、1，通過求解(e8～e10)方程組，獲得其正實數(shù)解，并對結果進行穩(wěn)定性分析，見圖6。

圖 6　不同初始小分子濃度C2(歸一化濃度)對于沉積反應的影響

從6圖可以看出當線性小分子烴C2初始歸一化濃度較小時，隨R的增加，表面將由MT沉積機制急劇轉變?yōu)镠T沉積機制，且在沉積過程中不存在遲滯回線；隨著初始線性小分子烴C2歸一化濃度的逐漸增加，(e8～e10)方程組系統(tǒng)開始出現(xiàn)非穩(wěn)態(tài)解，表面MT與HT熱解炭沉積機制的區(qū)域間開始出現(xiàn)遲滯現(xiàn)象，且遲滯回線區(qū)間逐漸增加，而MT/HT熱解炭沉積機制間的突變點也逐漸向右移動。

3.3沉積溫度對于熱解炭沉積的影響

假設初始線性小分子烴C2歸一化濃度a為1，在不同沉積溫度T為1 000、1 100、1 200 ℃條件下，通過求解(e8～e10)方程組，獲取其正實數(shù)解，并對結果進行穩(wěn)定性分析，見圖7。圖7與圖6規(guī)律相似，即在沉積溫度T較低時，隨R的增加，表面各組分覆蓋率會產生急劇變化，但并未出現(xiàn)遲滯的現(xiàn)象；隨著沉積溫度T的升高，表面MT與HT沉積機制的區(qū)域間開始出現(xiàn)遲滯且遲滯區(qū)間逐漸增加，而MT/HT沉積機制間的突變點也逐漸向左移動。即隨著沉積溫度的增加，體系由MT沉積機制突變?yōu)镠T沉積機制時所需達到的R值越小。說明在一定條件下升高沉積溫度T，可以使得系統(tǒng)由MT沉積機制轉變?yōu)镠T沉積機制。

圖 7　不同溫度下表面各組分的覆蓋率變化

固定C2的初始歸一化濃度a=1，在R分別為1、3、6時，求解(e7～e9)，取不同沉積溫度下的正實數(shù)解，并對結果進行穩(wěn)定性分析。得到不同R值下表面吸附的C2和 C6組分的濃度隨沉積溫度的變化曲線(圖8)，在CVI反應過程中，氣相烴類氣體的熱解以及熱解產物在基體表面的吸附、脫附、沉積以及表面雙分子反應本身都是吸熱與放熱的過程，因此即使在I-CVI過程中，體系內各處都會伴隨溫度的小范圍波動，當R很小時，體系內溫度的波動不會打破該處的反應平衡；只有當R大于某一臨界值時，體系內溫度的波動才有可能打破該處的反應平衡，如圖7，對應R為3和6時的表面各組分覆蓋率曲線，可以看出當體系處于C2下分支且遠離右端點或處于上分支且遠離左端點時(對應于C6上分支且遠離右端點或下分支且遠離左端點)，體系內微小的溫度波動不會打破該處的反應平衡。而當體系某處處于C2下分支的右端點(對應于C6上分支右端點)附近時，該處溫度略微的向下波動將會導致該處由HT沉積機制突變?yōu)镸T沉積機制；同樣地，而當體系某處處于C2上分支的左端點(對應于C6下分支左端點)附近時，該處溫度略微的向上波動將會導致該處由MT沉積機制突變?yōu)镠T沉積機制。

圖 8　表面吸附C2和C6組分的濃度隨溫度的變化

3.4討論

在CVI致密化過程，氣相中中大小分子濃度比R，初始小分子濃度以及沉積溫度T都對表面MT與HT熱解炭的形成有很大的影響。對于前兩種影響因素，可以通過原料氣的選擇以及初始濃度、預制體初始孔隙比表面積(A/V)、反應氣體在入口處的滯留時間來進行控制。

在低壓時，CVI系統(tǒng)中初始A/V相對較大，且隨著反應的進行A/V逐漸變的更大，氣相連續(xù)反應受阻，此時氣相組分主要由小分子直鏈烴組成，即R較小，由圖4和5可知此時表面主要由小分子C2占據(jù)，以MT熱解炭沉積為主且沉積速率較小。此時若想將體系由MT形成機制轉變?yōu)镠T形成機制，由圖7和8可知，提高體系的沉積溫度T可以使得MT沉積機制向HT沉積機制轉變，且R越小，需要的沉積溫度T越高；

在中壓時，甲烷分壓的增加，體系中大分子芳香烴C6持續(xù)增加，此時基體表面不易達到飽和吸附，當體系達到穩(wěn)定時，氣相中生成的C6等于基體表面因反應消耗的C6，從而阻止了C6在氣相中的持續(xù)長大。此時，R值較大，基體表面主要沉積HT熱解炭且沉積速率較快；隨著沉積的進行，A/V逐漸增加，氣相中C6的消耗速率逐漸增加，從而使R逐漸減小，由圖4可知如果R值降為R1，體系將會由HT形成機制突變?yōu)镸T形成機制；由圖7可知，若在此沉積過程中逐漸升高沉積溫度T，體系由HT沉積機制轉變?yōu)镸T沉積機制的突變點R1將會逐漸減小，從而避免因為沉積過程中由于A/V比的不斷增加導致系統(tǒng)由HT沉積機制向MT沉積機制轉變。

在高壓時，此時基體表面容易達到過飽和吸附，導致氣相中烴類氣體的不斷長大，此時，氣相主要由多環(huán)芳香烴與PAH組成，熱解炭生成向形核機制轉變，基體表面主要沉積MT熱解炭。

此結論很好的解釋了Rezinik等采用熱壁反應器在1 100 ℃下甲烷的ICVI實驗結果[16]，即當甲烷分壓為5 kPa時，纖維表面依次向外沉積了中織構與低織構炭：此時R較小，表面主要由小分子C2占據(jù)，隨著沉積的進行，A/V逐漸增大，R持續(xù)減小，因此整個過程以MT熱解炭沉積為主；當甲烷分壓為10 kPa時，纖維表面沉積的是純高織構炭：此時R較大，符合HT沉積機制，隨著沉積的進行A/V逐漸減小，R持續(xù)減小但未降至R1以下，因此整個沉積過程都處于HT沉積機制；而當甲烷分壓為20 kPa時，纖維表面先形成中織構炭后再形成高織構炭：初始沉積時基體表面已達到過飽和吸附，導致氣相中烴類氣體的不斷長大，導致MT熱解炭的沉積，隨著基體表面不斷地沉積A/V逐漸增大，使得整個反應系統(tǒng)脫離成核機制重新回歸生長機制開始沉積HT熱解炭。

Pauw在1 100 ℃，總壓為100 kPa(10%甲烷，90%氬氣)下研究了熱壁反應器中滯留時間τ對于在硅表面沉積的熱解炭織構演化的影響[17]，見圖9，當τ較小時主要沉積MT且沉積的炭顆粒平均直徑較小，只有當τ大于某一臨界值時才會沉積HT熱解炭且沉積的炭顆粒平均直徑較大(在該實驗條件下，臨界τ為1.4 s)。這主要是由于在τ較小時，氣相組分以小分子直鏈烴為主，基體表面主要是小分子的沉積(MT沉積機制)，所以沉積的炭顆粒平均直徑較??；當τ較大時，氣相組分不再以小分子為主，相當一部分已由大分子占據(jù)，此時表面主要發(fā)生的是雙分子反應(HT沉積機制)，此時沉積的炭顆粒平均直徑較大。

圖 9　在無氧條件下沉積1.5 h，Si表面沉積炭顆粒的

Hu等[6]利用熱壁反應器研究了甲烷壓強為25 kPa時，在不同滯留時間以及A/V值下董青石表面沉積的熱解炭織構隨沉積溫度的變化。如圖10所示，在較短的滯留時間(0.125 s)和高的A/V值條件下，熱解炭氣相沉積受化學生長模式控制，所得材料的織構在1 100～1 200 ℃存在最大值(OA最小值)，說明該條件下對高織構的生成最優(yōu)[18]。在此只關注溫度對于熱解炭沉積織構的影響，由圖8可知，在沉積溫度低于1 060 ℃左右時，不論R值多大，表面都是沉積MT熱解炭；隨著沉積溫度的增加，R值大的體系將先于R值小的體系由MT沉積機制突變?yōu)镠T沉積機制，這在圖10中得到了很好的驗證，即在A/V值為0.79 mm-1滯留時間為1 s時，體系由MT沉積突變?yōu)镠T沉積的臨界溫度為1 040 ℃左右，而滯留時間為0.125 s時，臨界溫度為1 080 ℃度左右；同樣在A/V值為3.2 mm-1也表現(xiàn)出了相似的趨勢。由圖7發(fā)現(xiàn)，隨著沉積溫度的升高，體系由MT沉積突變?yōu)镠T沉積的所需達到的R值越小，說明在生長機制階段沉積溫度越高，體系越容易沉積HT熱解炭，而在實驗中我們發(fā)現(xiàn)在A/V值為0.79 mm-1滯留時間為1 s時，當沉積溫度達到1 240 ℃時，體系將由HT沉積突變?yōu)镸T沉積，與模擬結果不符合，這主要是由于在的模型中未考慮氣相反應對于表面沉積的影響，隨著沉積溫度的升高，氣相烴類氣體的裂解速率逐漸大于氣相反應中間產物在基體表面的吸附速率而導致基體表面的飽和吸附，從而使得氣相組分逐漸長大，體系由生長機制轉變?yōu)槌珊藱C制，從而導致由HT沉積機制轉變?yōu)镸T沉積機制。Guellali等研究得到如果給于充足的滯留時間讓多環(huán)芳香烴分子內重排，形核機制也可以導致高織構熱解炭的生成[19]。

圖 10　甲烷壓強25 kPa、A/V值分別為0.79 mm-1(○/●)和

4結論

為理解CVI沉積過程中熱解炭微觀織構的突變現(xiàn)象，根據(jù)Langmuir-Hinshelwood理論并量化Particle-filler模型，建立了一個簡化的熱解炭沉積和織構形成模型，證實了在一定條件下熱解炭的沉積是一個雙穩(wěn)態(tài)過程(分別是以MT和HT沉積機制主導的穩(wěn)態(tài)過程)并存在一個遲滯區(qū)間，且此遲滯區(qū)間受氣相組分濃度比，初始直鏈烴C2濃度以及沉積溫度T的影響。量化的P-F模型能很好地解釋熱解炭結構的突變現(xiàn)象以及多種沉積條件下以甲烷為原料氣的熱解炭織構的演化過程。由于此模型將整個CVI過程中參與表面沉積的組分簡化成了小分子直鏈烴C2類與大分子C6類并未考慮其它組分的沉積，也未考慮氣相熱解反應對于表面吸附組分的影響，因此，不能完全解釋個別的工藝參數(shù)設置下CVI過程中熱解炭沉積過程。

下一步將適當?shù)脑黾訁⑴c表面沉積的氣相組分，引入Monto Carlo方法研究熱解炭的沉積過程，結合表面沉積機制研究整個CVI系統(tǒng)中熱解炭沉積過程，從而模擬整個CVI過程中熱解炭沉積演變過程。

參考文獻

[1]Fitzer E, Manocha L M. Carbon Reinforcements and Carbon-carbon Composites[M]. Berlin, Heidelberg: Springer, 1998: 190-226.

[2]Beznik B, Guellali M, Gerthsen D, et al. Microstructure and mechanical properties of carbon/carbon composites with multilayered pyrocarbon matrix[J]. Mater Lett, 2002, 52: 14-19.

[3]Hüttinger K J. CVD in hot wall reactors-The interaction between homogeneous gas-phase and heterogeneous surface reactions[J]. Chem Vap Deposition, 1998, 4: 151-158.

[4]Benzinger W, Becker A, Hüttinger K J. Chemistry and kinetics of chemical vapor deposition of pyrocarbon: I. Fundamentals of kinetics and chemical reaction engineering[J]. Carbon, 1996, 34: 957-966.

[5]Benzinger W, Hüttinger K J. Chemical vapor infiltration of pyrocarbon: I. Some kinetic considerations[J]. Carbon, 1996, 34: 1465-1471.

[6]Hu Z, Hüttinger K J. Influence of pressure, temperature and surface area/volume ratio on the texture of pyrolytic carbon deposited from methane[J]. Carbon, 2003, 41: 749-758.

[7]Hu Z, Hüttinger K J. Mechanisms of carbon deposition-a kinetic approach[J]. Carbon, 2002, 40: 617-636.

[8]Dong G, Hüttinger K J. Consideration of reaction mechanisms leading to pyrolytic carbon of different textures[J]. Carbon, 2002, 40: 2515-2528.

[9]Vignoles G L, Baconneau O, Brauner C M. Interaction between gas diffusion and multistable heterogeneous chemical kinetics in C/C composite processing[J]. Fundamentals of gas-phase and surface chemistry of vapor phase deposition II/Process control, diagnostics and modeling in semiconductor manufacturing IV, 2001, 13: 237-244.

[10]Delhaes P. Attempts to chemical vapour infiltrate pyrocarbons: evidence for a spatial bistability[J]. Carbon, 2003, 41: 1093-1095.

[11]Lacroix R, Fournet R, Ziegler D I, et al. Kinetic modeling of surface reactions involved in CVI of pyrocarbon obtained by propane pyrolysis[J]. Carbon, 2010, 132-144.

[12]徐偉, 張中偉, 白瑞成， 等. 丙烷化學氣相沉積均相熱解反應動力學模擬[J]. 新型炭材料, 2014, 29: 67-77.

(Xu W, Zhang Z, Bai R, et al. Kinetic modeling of gas-phase reactions for CVD from propane[J]. New Carbon Mater, 2014, 29: 67-77.)

[13]胡壽松.自動控制原理[M]. 北京: 科學出版社, 2008: 1-632.

(Hu S. Principle of Automatic Control[M]. Beijing: Science Press, 2008: 1-632.)

[14]Lv L. Modulation effect of control parameter on chemical bistable system[J]. Journal of Shenyang power higher school, 2003, 3: 3-4.

[15]Xu D. Routh table and the stability of linear time-invariant system[J]. Chinese Science Bulletin, 1984, 8: 454-456.

[16]Reznik B, Gerthsen D, Zhang W, et al. Texture changes in the matrix of an infiltrated carbon fiber felt studied by polarized light microscopy and selected area electron diffraction[J]. Carbon, 2003, 41: 369-384.

[17]Pauw V D, Collin A, Send W, et al. Deposition rates during the early stages of pyrolytic carbon deposition in a hot-wall reactor and the development of texture[J]. Carbon, 2006, 44: 3091-3101.

[18]張偉剛. 化學氣相沉積-從烴類氣體到固體炭[M]. 北京:科學出版社, 2007: 12-128.

(Zhang W. Chemical Vapor Deposition from Hydrocarbon Gas to Solid Carbon[M]. Beijing: Science Press, 2007: 1-40.)

[19]Guellali M. Textures of pyrolytic carbon formed in the chemical vapor infiltration of capillaries[J]. Carbon, 2003, 41: 97-104.

Foundationitem: Doctoral Fund of Ministry of Education(20113108120019); National Natural Science Foundation of China(11202124); Aeronautical Science Foundation of China(2013ZF6001); Fund of Shanghai Municipality(13521101202).

Authorintroduction: HUANG Qing-bo, Master Student. E-mail: hqb_0222@shu.edu.cn

Simulation of the kinetics of pyrolytic carbon deposition in C/C composites

HUANG Qing-bo1,ZHANG Dan2,BAI Rui-cheng1,LI Ai-jun1,SUN Jin-liang1

(1.ResearchCenterforCompositeMaterials,ShanghaiUniversity,Shanghai200072,China)(2.ShanghaiInstituteofAppliedMathematicsandMechanics,ShanghaiUniversity,Shanghai200072,China)

Abstract:A multistep heterogeneous reaction kinetic model for pyrocarbon deposition is proposed to investigate the sharp and clear transition between the high-texture (HT) and medium-texture (MT) pyrocarbons in C/C composites synthesized by chemical vapor infiltration (CVI). The model is based on the Langmuir-Hinshelwood mechanism and a particle-filler conceptual model, which models both the pyrocarbon deposition and the texture formation. The model assumes that adsorption, desorption and dehydrogenation reactions are involved. Unimolecular dehydrogenation reactions of either light linear hydrocarbons as the source of fillers (F) or light aromatic species as the source of particles (P) result in the formation of MT pyrocarbon, while a bimolecular dehydrogenation reaction between P and F species leads to the formation of HT pyrocarbon. The relationship between the types of pyrocarbons and gas-phase compositions is simulated under steady state. It is found that MT and HT pyrocarbon formation are two dominant stable processes with a hysteresis interzone that is affected by gas composition, initial linear hydrocarbon concentration and deposition temperature. Simulated results account for the sharp and clear transition between MT and HT pyrocarbon, and agree well with most pyrocarbon evolution studies under various conditions during CVI with only a few exceptions that may be caused by simplification in constructing the model.

Keywords:Carbon/carbon composites; Pyrolytic carbon; Texture; Kinetics

文章編號：1007-8827(2016)02-0167-09

中圖分類號:TQ342+.74

文獻標識碼:A

收稿日期：2016-01-05;修回日期：2016-04-02

基金項目：教育部博士點基金(20113108120019)；國家自然科學基金(11202124)； 航空科學基金(2013ZF6001)； 上海市科委基金(13521101202).

通訊作者：李愛軍，教授. E-mail：aijun.li@shu.edu.cn

作者簡介：黃清波，碩士研究生. E-mail: hqb_0222@shu.edu.cn

Corresponding author:LI Ai-jun, Professor. E-mail: aijun.li@shu.edu.cn