李若昕， 林柏洋， 顧浩， 林赫

(上海交通大學(xué)內(nèi)燃機(jī)所， 上海 200240)

?

摻混甲烷對乙烯層流預(yù)混火焰炭煙形貌的影響

李若昕， 林柏洋， 顧浩， 林赫

(上海交通大學(xué)內(nèi)燃機(jī)所， 上海 200240)

針對摻混甲烷的乙烯層流火焰中的炭煙顆粒從形貌角度進(jìn)行了研究，構(gòu)建了具有相同當(dāng)量比(2.07)、火焰溫度、氣體流速的層流預(yù)混乙烯火焰和摻混甲烷的乙烯火焰，并各自選取了3個高度進(jìn)行熱泳探針取樣，采用TEM和AFM方法對各個工況的炭煙顆粒粒徑分布、平均粒徑、平均體積當(dāng)量球徑進(jìn)行了對比。發(fā)現(xiàn)少量的甲烷摻混會令乙烯火焰中炭煙顆粒的尺寸變大，而當(dāng)甲烷摻混量變多時該現(xiàn)象消失。而通過對顆粒的圓度、球度和分形維數(shù)計(jì)算和分析，發(fā)現(xiàn)即使是很小的顆粒也具有相似的不規(guī)則形狀，這一現(xiàn)象廣泛存在于各種火焰和各種大小的顆粒中，不受摻混比的影響，并且越大的顆粒結(jié)構(gòu)越松散。

炭煙； 形貌； 熱泳探針； 取樣； 協(xié)同效應(yīng)

由化石燃料燃燒所產(chǎn)生的碳顆粒排放會對全球氣候、空氣質(zhì)量和人類健康產(chǎn)生嚴(yán)重危害[1-3]，20世紀(jì)90年代以來，燃燒過程炭煙顆粒的生成機(jī)理一直是研究的熱點(diǎn)。由Frenklach等人提出的Hydrogen-Abstraction-Carbon-Addition(HACA)機(jī)理[4-7]已被廣泛認(rèn)可，并成為炭煙生成模型的基礎(chǔ)理論之一。該理論認(rèn)為，大分子碳?xì)淙剂显谌紵龝r會在極短的時間內(nèi)裂解成小分子碳?xì)浠衔锶缫蚁?、丙烷等并進(jìn)行脫氫加碳反應(yīng)，然后不斷地凝并和聚合，最終形成以芳香族為核心的炭煙顆粒[8,19]。由此可見，對乙烯等小分子燃料的燃燒成煙過程進(jìn)行研究是研究大分子混合燃料燃燒過程中炭煙的生成機(jī)理的基礎(chǔ)。

在以前的小分子燃料研究中，有學(xué)者發(fā)現(xiàn)燃料摻混會對層流擴(kuò)散火焰的成煙過程造成影響[9-10]，也有人對苯、甲烷等燃料對PAHs生成所造成的影響進(jìn)行研究[7]。其中，Chung等人通過對炭煙體積分?jǐn)?shù)和數(shù)密度的測量，發(fā)現(xiàn)甲烷摻混對乙烯火焰中炭煙的生成有一定的促進(jìn)作用，并將這種自身成煙趨勢小、摻混后卻能夠促進(jìn)成煙的燃料摻混起到的促進(jìn)作用稱之為協(xié)同效應(yīng)(synergistic effect)[9-11]。Can Shao在隨后的研究中通過SMPS(掃描電遷移率顆粒物粒徑譜儀)對甲烷摻混量在協(xié)同效應(yīng)中的作用進(jìn)行了研究，并在模擬中探討了摻混對成煙的影響，認(rèn)為當(dāng)甲烷摻混量變大時反而會令協(xié)同效應(yīng)消失[12]。由于目前SMPS主要采用基于理想球形假設(shè)(假設(shè)炭煙顆粒為球狀)的模型計(jì)算方法計(jì)算顆粒的體積分?jǐn)?shù)，計(jì)算精確度與顆粒的真實(shí)形貌有很強(qiáng)的關(guān)聯(lián)性[16,20]。因此，本研究將從形貌角度對甲烷摻混對成煙的影響進(jìn)行分析，研究甲烷摻混所產(chǎn)生的協(xié)同效應(yīng)對形貌方面產(chǎn)生的影響。

炭煙顆粒形貌研究的主要目的是從粒徑和相應(yīng)的幾何參數(shù)入手，為炭煙模型的建立和修改提供形貌上的數(shù)據(jù)和參考[16]。在常見的形貌研究中，TEM(透射電子顯微鏡)是廣為使用的方法之一[18，22-25]，可以有效地觀測納米級顆粒。AFM(原子力顯微鏡)由于能夠獲得顆粒的三維形貌，也成為了研究的重要手段[16,21,23]。本研究結(jié)合TEM和AFM技術(shù)對乙烯摻混甲烷火焰中的炭煙形貌進(jìn)行研究，探究甲烷摻混對乙烯火焰中炭煙形貌的影響規(guī)律。

1 試驗(yàn)設(shè)備和方法

試驗(yàn)系統(tǒng)見圖1，左圖為燃燒器示意，右圖為熱泳取樣裝置示意。試驗(yàn)采用的燃燒器為自主設(shè)計(jì)的類McKenna預(yù)混火焰燃燒器，燃料氣體、氧氣和氬氣組成的混合氣經(jīng)由下部進(jìn)入，由壓鑄合金多孔材料整流后燃燒產(chǎn)生均勻的一維火焰，外圈通過氮?dú)庾鞅Ｗo(hù)氣與空氣隔離，保證預(yù)混的當(dāng)量比不會受到外界干擾，火焰半徑為50 mm，并通過循環(huán)水冷卻器和鎧裝熱電偶監(jiān)控并維持燃燒器溫度在20～25 ℃。在火焰的上方10 mm處放置了一塊金屬滯止板，同樣通過水冷和熱電偶確保溫度在105～110 ℃，最終火焰和保護(hù)氣在此處滯止，兩個邊界溫度穩(wěn)定，進(jìn)而構(gòu)建成一個標(biāo)準(zhǔn)的BSS(Burner Stabilized Stagnation)火焰，該火焰的優(yōu)點(diǎn)在于穩(wěn)定可控，重復(fù)性好，便于分析和驗(yàn)證[7，12,17-18]。

在氣體流量控制方面使用了經(jīng)典的體積流量控制，利用恒定溫度下的音速噴嘴在一定范圍內(nèi)壓力與體積流量成正比的特性，配合數(shù)顯壓力表對氣體的體積流量進(jìn)行控制來調(diào)控當(dāng)量比和流速等參數(shù)。

本研究采用熱泳探針取樣(TSPD，Thermophoretic Sampling Particle Diagnostic)方法從火焰中獲得炭煙樣本。TSPD方法從20世紀(jì)就開始在顆粒研究中得到使用[13-14]，在炭煙顆粒形貌研究中也得到廣泛使用[18，22-25]。本研究中的熱泳探針取樣系統(tǒng)采用了電力驅(qū)動的方式，比傳統(tǒng)的壓縮空氣驅(qū)動更為精準(zhǔn)。電缸為FESTO DNCE系列直線電缸，在驅(qū)動端安裝了自鎖鑷子用以夾持取樣基底，以減少探針對火焰的干擾。除了電缸操作系統(tǒng)自帶的時間曲線以外，又配置了兩個高精度光電門，用以調(diào)試和校正探針在火焰中的停留時間。根據(jù)前人經(jīng)驗(yàn)，停留時間選為20 ms[13，16,20，22-25]。

圖1 整體試驗(yàn)設(shè)備示意

研究中，火焰溫度由熱電偶涂覆后修正的方法進(jìn)行測量[12]，由于成煙過程受火焰溫度等因素的影響很大，所以在保證火焰溫度、當(dāng)量比(2.07)和流速一致的情況下，改變乙烯和甲烷的比例，設(shè)計(jì)了4個工況(見表1)，表中C3為文獻(xiàn)中的標(biāo)準(zhǔn)火焰[7]，T1，T2，T3則分別為摻混10%，20%，30%甲烷的工況。在10 mm的BSS火焰中，每個工況分別選取了HAB(Height Above Burner)6 mm，7 mm，8 mm作為測量點(diǎn)。

表1 摻混火焰預(yù)設(shè)工況

本研究采用的TEM為Tecnai G2 Spirit Biotwin，加載電壓為120 kV，基底選擇300目超薄碳膜。AFM為美國DI公司的MOLTIMODE NANOSCOPE，基底為直徑5 mm的云母片。在圖像分析中，使用的軟件為ImageJ和AFM自帶的分析軟件Nanoscope。

2 結(jié)果和討論

根據(jù)前人研究，C3火焰中從HAB 4 mm開始有炭煙顆粒產(chǎn)生[7，17,18]，因此選擇6 mm，7 mm，8 mm作為觀測點(diǎn)，并從粒徑、圓度、球度、分形維數(shù)等參數(shù)對顆粒進(jìn)行研究。

在模型的預(yù)測和SMPS分析中認(rèn)為，少量的甲烷摻混會增加成煙趨勢，產(chǎn)生協(xié)同效應(yīng)，而過量的甲烷摻混將會使其消失[12]。在顆粒的成長過程中，顆粒粒徑越大，說明反應(yīng)越充分、時間越充足。由于研究的對象是一維火焰，所以當(dāng)工況不同時，相同高度下的粒徑越大也就意味著成煙的趨勢越大，因此比較顆粒粒徑的變化將很直觀地反映出協(xié)同效應(yīng)的有無。

圖2示出典型的炭煙在TEM下的圖像，而由于炭煙顆粒為不規(guī)則形狀而非球形，這里使用形貌研究中比較通用的方法來計(jì)算其粒徑，也就是用一個最小的圓將炭煙顆粒包圍，并將其直徑作為該炭煙顆粒的粒徑[16，20]，圖2中黑色的圓圈直徑代表顆粒直徑。

圖2 TEM圖像(Flame C3,HAB 8 mm)

為了讓結(jié)果具有代表性，在每個工況點(diǎn)選取了20多張TEM照片近200多個顆粒進(jìn)行統(tǒng)計(jì)分析。圖3示出C3火焰中炭煙顆粒的粒徑分布情況，從分布圖中可以明顯看出，隨著高度的增加，大的顆粒開始變多，所占的比例逐漸增加。這是因?yàn)?，在同一火焰中，由于火焰高度的增加，生成的炭煙顆粒有著更多的時間和空間進(jìn)行長大、凝并和團(tuán)聚，這也就意味著顆粒的粒徑會隨著高度的增加而增加，另外，這個過程是持續(xù)而同時進(jìn)行的，所以小顆粒數(shù)量也不會明顯減少。

圖3 C3火焰各個高度下的粒徑分布

圖4示出摻混了甲烷后的T系列火焰中各個工況的粒徑分布。從圖中可以看到，在T1火焰中，相同高度下的炭煙粒徑分布明顯比C3火焰更加偏向于大顆粒一邊，也就是說，在T1火焰中，相同高度下成煙速度更快，具有更大的成煙趨勢。而在T2，T3火焰中，大顆粒的炭煙粒徑分布就沒有T1火焰那么明顯，甚至比C3還要低， T2，T3火焰中成煙的趨勢和C3火焰相差不多，甚至更小。

圖4 T系列火焰粒徑分布

為了更明顯地表示各個火焰中粒徑的變化，選取每個工況點(diǎn)的粒徑平均值作為該工況點(diǎn)的平均粒徑，并將它們繪制在一張圖表中。圖5示出4種火焰中每個高度下的平均粒徑，發(fā)現(xiàn)T1火焰明顯在各個高度上的粒徑都高于C3火焰，而T2，T3火焰則相差不多且低于C3火焰。

圖5 各工況下的平均粒徑對比

可以看出，在小量的甲烷摻混所產(chǎn)生的火焰中(10%T1),炭煙粒徑會明顯高于不摻混的C3火焰，而隨著甲烷摻混的比例加大(20%T2，30%T3)，炭煙粒徑則不會像少量摻混的工況一樣具有那么明顯的變化，甚至?xí)葥交熘斑€要更小。

由于TEM只能觀測到顆粒的二維圖像，為了觀察到顆粒的三維圖像，本研究使用了AFM進(jìn)行輔助分析。圖6a示出AFM顯示的典型顆粒三維形貌，圖6b示出一個球缺狀顆粒的側(cè)視圖及相關(guān)尺寸。從圖6a可以看出，顆粒一般呈現(xiàn)為球缺形或團(tuán)簇狀，如同液體滴在基底上一樣，這是由于初生炭煙呈半液態(tài)沉積在云母片上的緣故[16]；同時可以觀察到多個顆粒團(tuán)簇在一起的形態(tài)，這種形態(tài)是初生顆粒團(tuán)聚在一起所形成的特殊形態(tài)。對于圖6b中所示的球缺型顆粒，可以使用球缺體積公式V=πh(3a2+h2)/6計(jì)算其體積，其中a為底面半徑，h為球缺體高，再根據(jù)其體積計(jì)算體積當(dāng)量球徑Ds，也就是等同于其體積的標(biāo)準(zhǔn)球體的直徑。圖7示出了各個工況下平均體積當(dāng)量球徑的比較，可以看出和TEM分析中的平均粒徑趨勢相似，并略小于平均粒徑，這是因?yàn)門EM顯示的是顆粒的二維投影，并且會有一些重疊的顆粒存在，會導(dǎo)致兩個粒徑稍有差異。

從結(jié)果上來看，不管是從TEM的二維圖上，還是AFM的三維圖上，都能觀察到類似的協(xié)同效應(yīng)，但只在甲烷摻混量很小(10%)的時候出現(xiàn)。文獻(xiàn)中的計(jì)算表明，隨著甲烷的摻入，C2H4摩爾分?jǐn)?shù)下降，并導(dǎo)致其脫氫產(chǎn)物C2H2的摩爾分?jǐn)?shù)也下降，同時CH3和C3H3的摩爾分?jǐn)?shù)則會相應(yīng)增加[12]；C2H4，C2H2和CH3，C3H3都會對PAHs中的重要產(chǎn)物A4的生成造成影響，如HACA機(jī)理中常見的增加C2H2的反應(yīng)，兩個C3H3的成環(huán)反應(yīng)(ring formation)等等；因此，在少量的甲烷摻混(10%，T1)中，CH3和C3H3的增加非常顯著，C2H4和C2H2則只是少量的減少，從而促進(jìn)了PAHs的生成，使得成煙趨勢變大，顆粒成長速度變快；而當(dāng)甲烷摻混量加大的時候(20%，T2和30%，T3)，C2H4和C2H2顯著減少，相比之下CH3和C3H3的增長則不明顯[12]，使得成煙趨勢變小、顆粒成長速度降低。

圖6 AFM下的炭煙顆粒及球缺體參數(shù)

圖7 各個工況下當(dāng)量體積球徑的比較

圓度、球度和分形維數(shù)在形貌研究中常用于表示實(shí)際顆粒形態(tài)與球形假設(shè)的差異，以及顆粒結(jié)構(gòu)的疏密程度。圓度C=Pp2/Ap，式中：Pp為顆粒的周長；Ap為顆粒的投影面積。圓度C意味著顆粒與標(biāo)準(zhǔn)圓形的幾何對比，越接近圓形，圓度越小，而標(biāo)準(zhǔn)圓具有最小的圓度。球度S=Da/Dc，式中：Da為和投影面積相等的圓的直徑(圓面積當(dāng)量直徑)，Dc為圍住整個顆粒的最小圓的直徑，球度意味著顆粒的幾何形狀的尖銳程度，一條直線的球度為0，而標(biāo)準(zhǔn)圓的球度為1[16]。分形維數(shù)Df則采用了經(jīng)驗(yàn)公式的方法進(jìn)行計(jì)算[21-23]，顆粒的分形維數(shù)越小，意味著其結(jié)構(gòu)越疏散。圖8示出各個工況下的形貌參數(shù)的相應(yīng)比較。從圖中可以看出，各個工況下，彼此的圓度和球度相差不大，分別在20～60和0.2～0.6之間波動，并且可以看出，大部分顆粒的分形維數(shù)都處在相近的區(qū)間內(nèi)，而面積最大的一些顆粒的分形維數(shù)都比較小，這意味著它們的結(jié)構(gòu)更加松散。

圖8 各個工況下圓度、球度和分形維數(shù)的比較

從圖8可以看出，在各個火焰中生成的不同高度不同大小炭煙顆粒都具有明顯的不規(guī)則形狀，即使是很小的顆粒的形狀也不是圓形，并且在各個工況都存在這個現(xiàn)象，而圓度、球度和分形維數(shù)的變化都不是特別明顯，說明摻混對其影響很小。這是由于在HACA機(jī)理中，炭煙是由PAHs長大而生成的，以芳香族為核心，脂肪族為外圍不斷長大[19]，在這個過程中，外圍碳鏈不斷增大，這種增長是不規(guī)則的，而隨后的凝并過程則是將其結(jié)合成大的具有不規(guī)則形狀的顆粒，另外，顆粒團(tuán)聚的過程也是不規(guī)則的，這些綜合因素造成了不管是在任何一種摻混火焰中，任何大小的顆粒都不是圓形，而是呈現(xiàn)不規(guī)則碎片形(Fractal)。

3 結(jié)束語

本研究構(gòu)建了4種具有不同甲烷摻混量的火焰(0%，C3；10%，T1；20%，T2；30%，T4)，并分別在3個典型高度進(jìn)行了熱泳探針取樣和相應(yīng)的TEM與AFM觀測分析，根據(jù)試驗(yàn)和分析結(jié)果可知，從形貌角度看，當(dāng)少量甲烷摻混時，乙烯火焰中的炭煙顆粒成長速度變快，相同火焰高度下產(chǎn)生的炭煙顆粒粒徑明顯變大，而隨著更多量的甲烷摻混，乙烯火焰中的炭煙顆粒平均粒徑反而比原來更小，分布上也趨向于小顆粒。表明小劑量甲烷摻混對乙烯火焰的成煙具有明顯的增強(qiáng)作用，也就是會產(chǎn)生協(xié)同效應(yīng)；而摻混更多的甲烷反而會抑制和減少乙烯火焰的成煙趨勢，協(xié)同效應(yīng)消失。這個現(xiàn)象是甲烷摻混對乙烯火焰中PAHs生成造成的雙重作用而產(chǎn)生的綜合結(jié)果，和模型預(yù)測的結(jié)果以及SMPS試驗(yàn)所得到的結(jié)果相吻合。此外，在各個火焰中，顆粒從小到大都一致具有相似特征的不規(guī)則形狀，即使是很小的顆粒也不是圓形，摻混和協(xié)同效應(yīng)對其的影響微乎其微。

[1] Frances C Moore.Climate change and air pollution:exploring the synergies and potential for mitigation in industrializing countries[J].Sustainability,2009(1)：43-54.

[2] Bond T C,Doherty S J,Fahey D W,et al.Bounding the role of black carbon in the climate system:A scientific assessment[J].J.Geophys. Res. Atmos.，2013，118：5380-5552.

[3] Frank B,Schuster M E,Schl?gl R，et al.Emission of Highly Activated Soot Particulate：The Other Side of the Coin with Modern Diesel Engines[J].Angew. Chem. Int. Ed.，2013，52：2673-2677.

[4] Frenklach M. Reaction mechanism of soot formation in flames[J].Phys. Chem. Chem. Phys.，2002(4)：2028-2037.

[5] Frenklach M,Wang H.Detailed modeling of soot par-ticle nucleation and growth[J].Proc. Combust. Inst，1990，23：1559-1566.

[6] Frenklach M,Wang H.Soot Formation in Combustion:Mechanisms and Models of Soot Formation[J].Springer Series in Chemical Physics,1994，59：162-190.

[7] Abid A D,Heinz N,Tolmachoff E D,et al.On evol-ution of particle size distribution functions of incipient soot in premixed ethylene-oxygen-argon flames[J].Combust. Flame，2008，154：775-788.

[8] Yoon S,Lee S,Chung S.Effect of mixing methane,ethane,propane,and propene on the synergistic effect of PAH and soot formation in ethylene-base counterflow diffusion flames[J].Proc Combust Inst，2005，30:1417-1424.

[9] HWANG J Y,LEE W,KANG H G,et al.Synergistic Effect of Ethylene-Propane Mixture on Soot Formation in Laminar Diffusion Flames[J].Combust. Flame，1998，114:370-380.

[10] Lee S M,Yoon S S,Chung S H.Synergistic effect on soot formation in counterflow diffusion flames of ethylene-propane mixtures with benzene addition[J].Combustion and Flame,2004，136:493-500.

[11] Yoon S,Lee S,Chung S.Effect of mixing methane,ethane,propane,and propene on the synergistic effect of PAH and soot formation in ethylene-base counterflow diffusion flames[J].Proc Combust Inst,2005，30:1417-1424.

[12] Can Shao,Bin Guan,Baiyang Lin,et al.Effect of methane doping on nascent soot formation in ethylene-based laminar premixed flames[J].Fuel,2016，186：422-429.

[13] Dobbins R A,Megaridis C M.Morphology of Flame-Generated Soot As Determined by Thermophore tic Sampling[J].Langmuir,1987(3):254-259.

[14] Robert P Lucht.Applied combustion diagnostics[J].Combustion and Flame,2002,133:507-508.

[15] Jeonghoon Lee,Igor Altman,Mansoo Choi.Design of thermophoretic probe for precise particle sampling[J].Aerosol Science,2008(39):418-431.

[16] Marina Schenk,Sydnie Lieb,Henning Vieker,et al.Morphology of nascent soot in ethylene flames[J].Proceedings of the Combustion Institute,2015(35)：1879-1886.

[17] Gu C,Lin H,Camacho J,et al.Particle size distri-bution of nascent soot in lightly and heavily sooting premixed ethylene flames[J].Combust Flame,2016，165:177-187.

[18] Hai Wang.Formation of nascent soot and other condensed-phase materials in flames[J].Proceedings of the Combustion Institute，2011，33：41-67.

[19] Jeremy P Cain,Joaquin Camacho,Denis J Phares,et al.Evidence of aliphatics in nascent soot particles in premixed ethylene flames[J].Proceedings of the Combustion Institute,2011，33：533-540.

[20] Marina Schenk,Sydnie Lieb.Imaging Nanocarbon Materials:Soot Particles in Flames are Not Structurally Homogeneous[J].Chem Phys Chem,2013，14:3248-3254.

[21] Nan Chen,Chonglin Song,Gang Lv,et al.Atom force microscopy analysis of the morphology,attractive force,adhesive force and Young's modulus of diesel in-cylinder soot particles[J].Combust Flame，2015，162:4649-4659.

[22] 李政.預(yù)混火焰分析系統(tǒng)開發(fā)及反應(yīng)工況對炭煙結(jié)構(gòu)特性的影響[D].天津:天津大學(xué),2011.

[23] 洪亮.預(yù)混火焰炭煙顆粒的形貌結(jié)構(gòu)特征及團(tuán)聚力的研究[D].天津:天津大學(xué),2012.

[24] 李朝.基于熱泳取樣的乙烯富氧擴(kuò)散火焰中碳黑生成及形貌檢測[D].武漢:華中科技大學(xué),2012.

[25] 梁軍輝.氧濃度對擴(kuò)散火焰中煙黑的生成及分布特性的影響研究[D].杭州:浙江大學(xué),2012.

[編輯： 潘麗麗]

Influence of Methane Doping on Soot Morphology of Ethylene Laminar Premixed Flame

LI Ruoxin, LIN Boyang, GU Hao, LIN He

(Internal Combustion Engine Research Institute, Shanghai Jiao Tong University, Shanghai 200240, China)

Soot of ethylene laminar premixed flame mixed with methane was researched from the view of morphology. Ethylene laminar premixed flames and ethylene flames mixed with methane which had the same equivalent ratio (2.07), flame temperature and gas velocity were first produced. Three samples were sampled from each flame at different heights by TSPD, and then the distributions of soot diameters, mean diameters and mean equivalent spherical diameters of different soot samples were compared through the TEM and AFM analysis methods. It was noticed that low proportion of methane could make the soot size become large while more methane would led to the disappearance of this phenomenon. Through further calculation and analysis of the circularity, sphericity and fractal dimension, the very small particles were found to be not circular but fractal, which existed widely in different flames and different size of particles. Besides, the phenomenon was immune to the mixing ratio and the larger particles had looser shapes.

soot； morphology； TSPD； sampling； synergistic effect

2016-12-13；

2017-03-10

李若昕(1990—)，男，碩士，主要從事火焰中炭煙顆粒的形貌分析的相關(guān)工作；liruoxin1234@sina.cn。

10.3969/j.issn.1001-2222.2017.03.003

TK411.2

B

1001-2222(2017)03-0014-06