伍飛飛 廖瑞金? 楊麗君 劉興華 汪 可 周 之

1)(輸配電裝備及系統(tǒng)安全與新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室(重慶大學(xué))，重慶 400030)

2)(淄博供電公司，山東省電力集團(tuán)公司，淄博 255000)

(2012年11月5日收到;2012年12月12日收到修改稿)

1 引言

電暈放電是一種非平衡態(tài)低溫等離子體放電過程，通常發(fā)生在曲率半徑很大的尖端電極附近[1].隨著我國電力工業(yè)的發(fā)展，高壓和超高壓輸電線路成為我國主要骨干網(wǎng)架并即將建立特高壓線路，電暈損失已不可忽視[2].同時(shí)，電暈放電帶來的可聽噪聲、電磁干擾、線路腐蝕、絕緣劣化以及擾民問題已日漸受到重視[3]，此外，電暈放電在靜電除塵、污水處理、空氣凈化、表面處理等工業(yè)設(shè)備中起著核心作用[4,5].特里切爾脈沖對于解釋電暈放電的眾多現(xiàn)象發(fā)揮著至關(guān)重要的作用，國內(nèi)外雖開展了大量關(guān)于電暈產(chǎn)生機(jī)理的研究，但至今未取得突破性進(jìn)展.對于電暈放電中特里切爾脈沖的形成機(jī)理，自1938年Trichel發(fā)現(xiàn)了特里切爾脈沖以來一直受到廣大研究者的青睞[6].特里切爾脈沖是只在電負(fù)性氣體中存在，通常表現(xiàn)為迅速上升(小至ns級)和緩慢衰減(數(shù)十ns)的一組規(guī)律脈沖電流信號[7,8].Trichel本人認(rèn)為間隙中負(fù)空間電荷的形成和消散導(dǎo)致了脈沖電流的重復(fù)出現(xiàn).Loeb等的試驗(yàn)[9]表明規(guī)律的脈沖電流只存在于電負(fù)性氣體之中，同時(shí)Loeb提出了負(fù)直流電暈特里切爾脈沖空間電荷運(yùn)動(dòng)過程的物理模型[10]，該模型考慮了電子崩過程，但無法解釋他試驗(yàn)中發(fā)現(xiàn)的1.5 ns快速上升沿.Zentner[7]開展了大量試驗(yàn)測量特里切爾脈沖的上升時(shí)間，并且發(fā)現(xiàn)上升沿可小至1.3 ns.Lama和Gallo[11]豐富了Trichel和Loeb提出的特里切爾脈沖的微觀過程，他們認(rèn)為放電起始階段棒尖端因湯森雪崩產(chǎn)生大量電子，接著，被中性分子吸附形成負(fù)離子使得運(yùn)動(dòng)速度減慢，形成的負(fù)離子云將棒尖附近電場降低至臨界值以下，導(dǎo)致新的電子崩無法形成，使得放電過程發(fā)生休克.只有當(dāng)負(fù)離子云遷移到距離陰極足夠遠(yuǎn)處，陰極附近場強(qiáng)恢復(fù)之后，新的特里切爾脈沖脈沖開始出現(xiàn).Lama和Gallo的研究無法對離子云內(nèi)部的粒子成分進(jìn)行分析，也無法對特里切爾脈沖過程中的其他重要微觀特征量進(jìn)行描述.

盡管國內(nèi)外眾多研究者開展了大量的試驗(yàn)研究，特里切爾脈沖的物理機(jī)理仍未明了.從試驗(yàn)方面推進(jìn)空氣放電機(jī)理的研究，需要花費(fèi)大量的時(shí)間精力，耗資巨大，且目前缺乏有效的關(guān)于大氣壓下空氣放電等離子體診斷手段[12]，放電過程中的很多特征量都還不能通過試驗(yàn)獲得.與此同時(shí)，國內(nèi)外的眾多學(xué)者在利用數(shù)學(xué)手段開展特里切爾脈沖微觀過程的研究中取得了大量成果.Kekez等[13]采用了一種改進(jìn)的電路方法去描述特里切爾脈沖的連續(xù)過程，然而他的模型無法對脈沖的形成機(jī)理進(jìn)行解釋.Tran等[14]利用ComsolMultiphysics軟件提出特里切爾脈沖的二維計(jì)算模型，然而他們的求解結(jié)果無法和相同電極結(jié)構(gòu)的試驗(yàn)數(shù)據(jù)相容，脈沖周期減小到試驗(yàn)數(shù)據(jù)的1/20.考慮粒子之間化學(xué)反應(yīng)和非平衡態(tài)的大氣壓空氣放電成為未來進(jìn)行科學(xué)研究的一個(gè)極其重要的方向[15].Nahomy等[16]在研究氣體放電過程中考慮了430種化學(xué)反應(yīng)過程，是目前考慮化學(xué)反應(yīng)類型最多的研究，為后續(xù)科研工作者奠定了良好的基礎(chǔ).Pancheshnyi等[17]在Nahomy基礎(chǔ)上選取了10種主要粒子提出N2：O2(9：1)混合氣體在壓強(qiáng)為760 Torr(1 Torr=1.33322×102Pa)時(shí)流注放電二維物理模型，利用該模型對電場分布、帶電粒子密度、反應(yīng)系數(shù)等放電參數(shù)進(jìn)行了詳細(xì)分析并且和試驗(yàn)數(shù)據(jù)具有較好的一致性.縱觀目前的研究成果可知，能夠全面反應(yīng)電暈放電特里切爾脈沖過程的的宏觀及微觀過程的物理模型尚欠缺，開展特里切爾脈沖的微觀機(jī)理研究具有重要的理論和實(shí)際價(jià)值.

本文基于二維流體動(dòng)理學(xué)模型，提出一種改進(jìn)的負(fù)電暈放電混合數(shù)值模型.該模型加入了12種粒子及它們之間的27種碰撞反應(yīng)，同時(shí)考慮光電離和二次電子發(fā)射作用.利用該模型對棒-板間距3.3 mm，施加電壓-5.0 kV情況下進(jìn)行數(shù)值計(jì)算，得到負(fù)電暈放電的特里切爾脈沖.重點(diǎn)分析了5個(gè)關(guān)鍵時(shí)刻的微觀特征量在特里切爾脈沖發(fā)展過程中的變化規(guī)律.

2 數(shù)值模型

2.1 控制方程

在流體動(dòng)力學(xué)模型中考慮碰撞反應(yīng)，并在混合數(shù)值模型中加入光電離項(xiàng)和二次電子發(fā)射過程來分析空氣放電的物理過程，其實(shí)質(zhì)就是將電子連續(xù)性方程、重粒子多組分?jǐn)U散輸運(yùn)方程、重粒子動(dòng)量方程、能量方程和泊松方程化為適當(dāng)?shù)钠⒎址匠探M，再將偏微分方程組歸一化后以離散的數(shù)值差分形式求解.

電子的連續(xù)控制方程[18]為

其中，Te是電子溫度，Γe是電子通量，nk是粒子k的密度，ηk是電子與重粒子k碰撞的能量損失系數(shù)，me是電子質(zhì)量，Mk是重粒子k質(zhì)量，T是空氣溫度，ve,n是電子動(dòng)量傳遞碰撞頻率.方程(3)右邊為源項(xiàng)：je·E表示電子由于焦耳加熱效應(yīng)獲得的能量，最后兩項(xiàng)分別表示非彈性碰撞和彈性碰撞的能量損失.

本模型中，對于重粒子(例如，正離子、負(fù)離子和中性粒子)，多組分輸運(yùn)方程組[22]表示為

其中，ρ為總的密度，u平均流體速度矢量，Rk為粒子k由于化學(xué)反應(yīng)的變化速率，jk為擴(kuò)散通量矢量[23]，

重粒子的動(dòng)量方程[25]表示為

2.2 碰撞反應(yīng)類型

表1 模型中考慮的主要碰撞反應(yīng)

2.3 邊界條件

在混合模型中，邊界條件是數(shù)值計(jì)算過程非常重要的部分，電子通量在陰極和陽極的邊界條件[30]如下：

其中，Γe和Γi分別是電子密度通量和離子密度通量，υth,e是電子熱速率，γ為二次電子發(fā)射系數(shù)，在陰極設(shè)為0.004，在陽極設(shè)為0.

其中，υth,i和υth,s分別為離子和中性粒子的熱速率系數(shù).

圖1所示為負(fù)電暈放電外電路及計(jì)算區(qū)域示意.棒電極的曲率半徑為0.4 mm，板電極半徑為5.0 cm，棒-板間距3.3 mm.外電路是由DC源，電容C和保護(hù)電阻R組成.與直接連接DC電源相比，外電路的作用就是為極板之間提供穩(wěn)定放電電壓.在模型中，直流電壓取值范圍為-5.0 kV.電容C設(shè)為1.0 pF，并且保護(hù)電阻R為5.0 kΩ.

在棒尖端，假設(shè)空氣放電前存在的粒子(種電子和正離子)是最大值為1016m-3的高斯分布[14]：

其中，Nmax=1016m-3，r0=0μm，z0為棒尖的位置，s0=25μm.初始條件已經(jīng)被證明僅僅會加快脈沖的形成，不會改變放電特性[31].

圖1 負(fù)電暈放電外電路及計(jì)算區(qū)域示意圖

圖2 計(jì)算流程圖

2.4 計(jì)算實(shí)現(xiàn)

電暈放電物理計(jì)算模型通過軸對稱旋轉(zhuǎn)簡化為二維棒-板結(jié)構(gòu)(圖1虛線框內(nèi)部分)，該混合模型是通過基于有限元方法的ComsolMultiphysics軟件的等離子體模塊來求解實(shí)現(xiàn)的.環(huán)境條件設(shè)置為300 K和1.0 atm(1atm=1.01325×105Pa)，計(jì)算網(wǎng)格總共劃分為244884個(gè)單元.計(jì)算流程圖如圖2所示.

3 數(shù)值計(jì)算結(jié)果

3.1 特里切爾脈沖波形

圖3所示為特里切爾脈沖波形.結(jié)果表明：電流在t1-t2過程迅速上升到0.195 A，隨后開始緩慢下降，經(jīng)過60 ns的衰減之后維持在大約1 mA，波形特性與Zentner和Soria[7,8]的試驗(yàn)結(jié)果符合較好.本文選取特里切爾脈沖持續(xù)過程中的5個(gè)時(shí)間點(diǎn)(t1-t5)來分析負(fù)電暈特里切爾脈沖的微觀發(fā)展過程.

圖3 特里切爾脈沖波形

3.2 電場分布特性

在一個(gè)特里切爾脈沖持續(xù)過程中，電場強(qiáng)度分布對于分析放電的微觀物理過程具有相當(dāng)重要的意義.圖4表明了電場強(qiáng)度軸向分布的發(fā)展規(guī)律.在特里切爾脈沖起始階段，電場強(qiáng)度主要集中在陰極鞘附近，且迅速增大至最大值.隨著放電時(shí)間的發(fā)展，電場集中分布區(qū)域向陽極移動(dòng)且幅值變小.在特里切爾脈沖后期則電場分布重新集中至陰極附近.這是由于電子崩往陽極發(fā)展的過程中，因?yàn)榕鲎埠臀椒磻?yīng)產(chǎn)生大量帶電離子，離子云和電子形成的內(nèi)部場導(dǎo)致了棒-板間的電場畸變，當(dāng)負(fù)離子云遷移至陽極，間隙重新恢復(fù)電中性，一個(gè)特里切爾脈沖過程結(jié)束.陽極鞘附近的電場隨著放電時(shí)間的發(fā)展有略微變大的趨勢，這是因?yàn)檫_(dá)到陽極鞘的電子數(shù)量逐漸變多，陽極和鞘層電子的內(nèi)部場加強(qiáng)了陽極鞘附近的電場.由圖4可以發(fā)現(xiàn)，電子崩的發(fā)展方向就是電場變?nèi)醯姆较?，對電子崩的發(fā)展變得非常不利，這就是同樣棒-板結(jié)構(gòu)情況下，棒極為負(fù)極性時(shí)的擊穿電壓大于正極性的原因.

圖4 電場強(qiáng)度軸向分布

3.3 凈空間電荷分布發(fā)展規(guī)律

圖5是一個(gè)特里切爾脈沖過程中5個(gè)不同時(shí)間的凈電荷分布.由圖可知，在整個(gè)棒-板間隙，大部分放電區(qū)域都是電中性的，只有在陰極鞘和陽極鞘附近有帶正電的等離子體特性.在特里切爾脈沖起始階段，離子云集中在陰極附近，電荷密度因強(qiáng)場電離和電子碰撞反應(yīng)迅速增加，快速增加的正離子電荷被棒極吸收而使得放電電流具有快速的上升沿.隨著時(shí)間的發(fā)展，帶負(fù)電性的離子云，向陽極發(fā)散似移動(dòng)，當(dāng)離子云發(fā)展至陽極時(shí)，間隙的帶電粒子迅速被吸收，間隙重新恢復(fù)電中性.

3.4 電子密度分布特性

電子和空氣中中性分子的碰撞電離是空氣放電過程中最為主要的電離過程，此外，電子是外加電場與重離子之間傳遞能量的主要載體，因此研究特里切爾脈沖過程中電子密度分布規(guī)律對于進(jìn)一步揭示電暈放電的微觀物理機(jī)理具有重要意義.圖6為一個(gè)特里切爾脈沖過程中5個(gè)不同時(shí)刻的軸向電子密度分布.由圖可以發(fā)現(xiàn)，在陰極鞘內(nèi)，電子密度近似為0，這是因?yàn)殡娮颖慌懦鉄o法進(jìn)入陰極[32].電子密度在陰極鞘外層具有最大值，這主要是由于棒極端附近電場最為集中，場致發(fā)射效應(yīng)最強(qiáng).隨著時(shí)間增加，體等離子體內(nèi)電子密度繼續(xù)增大.這是因?yàn)槌蜿枠O的基本電子崩建立起來，激勵(lì)反應(yīng)和電離反應(yīng)也迅速增強(qiáng)的緣故.

圖5 凈電荷密度二維分布(C/m3) (a)t1(1.0 ns);(b)t2(2.6 ns);(c)t3(20 ns);(d)t4(55 ns);(e)t5(62 ns)

圖6 軸向電子密度分布

從t5時(shí)刻的電子密度可以看出，靠近陽極時(shí)電子密度迅速增大，這種現(xiàn)象是因?yàn)閹щ娏Ｗ幼矒絷枠O表面發(fā)生二次電子發(fā)射導(dǎo)致的，此時(shí)電暈放電轉(zhuǎn)入自持階段[33].陰極鞘附近的電子在陰極強(qiáng)電場下開始加速并達(dá)到空氣分子的激發(fā)能，此時(shí)和中性粒子發(fā)生碰撞形成激發(fā)態(tài)，而處于激發(fā)態(tài)的粒子因容易發(fā)生輻射復(fù)合而使得電子密度有個(gè)減少的過程.隨后，電子在電場加速下能量達(dá)到電離能，使得電離程度大大提高，電子密度逐漸升高.此時(shí)電子主要有兩類，一類是快電子，這部分電子從陰極表面產(chǎn)生后，一直被電場加速到該處，另一類是慢電子，這部分電子從陰極表面發(fā)射出來后雖然經(jīng)過電場加速，但是同時(shí)經(jīng)歷了多次非彈性碰撞，電子能量小于電離能(接近激發(fā)能)，所以這個(gè)分布帶具有強(qiáng)烈的電離和激發(fā)過程.62 ns時(shí)刻整個(gè)放電區(qū)域內(nèi)電場強(qiáng)度變化不大，且低于大氣壓空氣電離強(qiáng)度，高能電子的能量以彈性碰撞和非彈性碰撞的形式消耗掉，電離作用減弱，所以靠近陽極表面電子密度會出現(xiàn)一個(gè)下降過程.隨后電子又重新獲得加速，使得電離過程變強(qiáng)，且電子崩通過陽極發(fā)生中和反應(yīng)，使得陽極附近的電子密度在整個(gè)放電區(qū)域是最大的.空氣放電過程中通常會伴隨著大量發(fā)光輻射現(xiàn)象，發(fā)光強(qiáng)度取決于電子和其他粒子間碰撞能量轉(zhuǎn)移過程.圖7所示為t5時(shí)刻電子密度分布的二維分布圖，和Antao等[34]所得到的負(fù)電暈放電成像圖片有著極其相近的分布.

圖7 t5(62 ns)時(shí)的電子密度二維分布

3.5 特里切爾脈沖后期的帶電粒子特性

圖8和圖9分別是負(fù)電暈放電過程中特里切爾脈沖后期電子的產(chǎn)生和消失速率，該模型中考慮的與電子產(chǎn)生和消失有關(guān)的反應(yīng)有8種.從圖中可以發(fā)現(xiàn)，對于電子產(chǎn)生速率影響最大的兩種碰撞反應(yīng)是R1和R2.R1和R2分別為與N2和O2有關(guān)的碰撞電離反應(yīng)，兩者分布曲線大致相同，但由于O2的電離閾值能量(12.06 eV)比N2的電離閾值能量(15.6 eV)要低，使得放電間隙內(nèi)R2的反應(yīng)速率略大于R1.電子消失的主要反應(yīng)包括R12，R13，R14，R22，R23和R27，其中R22是涉及N2和N+2和電子的三體復(fù)合反應(yīng)，在電子消失過程中占絕對優(yōu)勢.R12，R13，R14，R23和 R27對電子消失的作用很小.

圖8 電子的產(chǎn)生速率

圖9 電子的消失速率

圖9是特里切爾脈沖后期正負(fù)離子密度軸向分布曲線.由圖中可以發(fā)現(xiàn)，正離子密度分布曲線與電子密度分布曲線類似，這是因?yàn)檎x子主要是由電子的碰撞電離產(chǎn)生.與正離子密度相比，負(fù)離子密度的明顯特征是從陽極到陰極過程中緩慢增加，這是因?yàn)樨?fù)離子主要是由電子和中性分子(原子)的附著反應(yīng)而產(chǎn)生，在靠近陰極鞘的區(qū)域電子密度較大，另外由于在陰極鞘內(nèi)存在較大電場強(qiáng)度的緣故，使得負(fù)離子快速向陽極運(yùn)動(dòng)，使得陰極附近密度較低.由圖示坐標(biāo)可以發(fā)現(xiàn)，負(fù)離子密度比陽離子且整體要小1個(gè)數(shù)量級左右.這是因?yàn)樵趶?qiáng)場作用下，主要發(fā)生的是電離反應(yīng)，產(chǎn)生大量正離子;O2分子難以捕捉到運(yùn)動(dòng)速度很快的電子發(fā)生吸附反應(yīng)形成負(fù)離子;同時(shí)，負(fù)離子容易和周圍大量的正離子之間發(fā)生碰撞反應(yīng)而形成中性粒子而降低負(fù)離子密度.

圖10 正負(fù)離子密度軸向分布規(guī)律

4 結(jié)論

本文詳細(xì)介紹了一種基于流體動(dòng)力學(xué)理論改進(jìn)得到的適用于棒-板電極負(fù)電暈放電微觀過程研究的混合數(shù)值模型.利用該模型，本文對棒-板間距3.3 mm，施加電壓-5.0 kV進(jìn)行數(shù)值計(jì)算，得到負(fù)電暈放電的特里切爾脈沖.通過分析一個(gè)特里切爾脈沖持續(xù)過程中5個(gè)時(shí)間點(diǎn)的電場分布、凈空間電荷分布、電子密度分布及特里切爾脈沖后期的帶電粒子等特征量描述特里切爾脈沖微觀發(fā)展過程.

圖11 帶電離子密度軸向分布

主要結(jié)論如下：

1.豐富并量化描述了Lama和Gallo提出的特里切爾脈沖的重復(fù)機(jī)理：放電起始階段棒尖端因湯森雪崩產(chǎn)生大量電子和正離子，正離子在從陰極迅速吸收電子恢復(fù)中性，這個(gè)過程導(dǎo)致了特里切爾脈沖的快速上升沿;電子在電場力作用下向陽極移動(dòng)時(shí)，被中性分子吸附變成負(fù)離子，負(fù)離子云降低了棒尖附近電場至臨界值以下，導(dǎo)致新的電子崩無法形成，使得放電過程發(fā)生休克;雪崩產(chǎn)生的電子向陽極移動(dòng)的過程中不斷發(fā)生碰撞反應(yīng)和吸附反應(yīng)，當(dāng)負(fù)離子云緩慢遷移到陽極過程中間隙電流隨之緩慢衰減，當(dāng)負(fù)離子云遷移至陽極時(shí)，陰極附近場強(qiáng)恢復(fù)至放電起始階段，新的特里切爾脈沖開始.

2.在一個(gè)特里切爾脈沖持續(xù)時(shí)間內(nèi)，隨著放電時(shí)間的增加，電場集中分布區(qū)域向陽極發(fā)展，幅值逐漸減小.電子崩的發(fā)展方向就是電場變?nèi)醯姆较?，對電子崩的發(fā)展非常不利.

3.在特里切爾脈沖持續(xù)過程中，大部分放電區(qū)域都是電中性的，只有在陰極鞘和陽極鞘附近有帶正電的等離子體特性，帶負(fù)電的離子云隨著放電時(shí)間的發(fā)展緩慢向陽極發(fā)散式移動(dòng).

4.在一個(gè)特里切爾脈沖持續(xù)時(shí)間內(nèi)，陰極鞘內(nèi)電子密度幾乎為0;特里切爾脈沖起始階段陰極鞘附近電子密度迅速增加至最大值并保持基本不變;隨著時(shí)間增加，放電間隙內(nèi)電子密度整體增加，并且有往陽極發(fā)展的趨勢;當(dāng)電子崩發(fā)展至陽極，陽極鞘內(nèi)，電子密度出現(xiàn)最大值，此時(shí)放電間隙迅速恢復(fù)電中性.

[1]Liu Z Y 2005 Ultra-hig grid(Beijing：China Economic Publishig)(in Chinese)[劉振亞2005特高壓電網(wǎng)(北京：中國經(jīng)濟(jì)出版社)]

[2]Shu Y B，Hu Y 2007 Proceedings of the CSEE 27 1(in Chinese)[舒印彪，胡毅2007中國電機(jī)工程學(xué)報(bào)27 1]

[3]Zheng Y S，He J L，Zhang B 2011 High Voltage Engineering 37 752(in Chinese)[鄭躍勝，何金良，張波2011高電壓技術(shù)37 752]

[4]In L，WU B，Zhang P，Wang Y Q 2004 Chin.Phys.Lett.21 1993

[5]Stoffels E，F(xiàn)likweert A J，Stoffels W W and Kroesen G M W 2002 Plasma Sources Sci.Technol.11 383

[6]G W Trichel1938 Phys.Rev.54 1078

[7]Zentner R 1970 Z.Angew.Physik 29 294

[8]Soria H C，Pontiga F，Castellanos A 2007 J.Phys.D：Appl.Phys.40 4552

[9]Loeb L B，Kip A F，Hudson G G 1941 Phys.Rev.60 714

[10]Loeb L B 1965 Electricalcoronas:their basic physicalmechanisms(Berkeley：University of California Press)p226—230

[11]Lama W L，GalloC F 1974 J.Appl.Phys.45 103—13

[12]MichaelA L，Allan J L 2007 Plasma discharge principle and materials processing(Beijing：Science Press)(in Chinese)[邁克爾 A 力伯曼，阿倫J里登伯格2007等離子體放電原理與材料處理(北京：科學(xué)出版社)]

[13]Kekez M M，Savic P，Lougheed G D 1982 J.Phys.D：Appl.Phys.15 1963

[14]Tran T N，Golosnov I O，Levin P L，Georghiou G E 2009 IEEE Conf.on ElectricalInsulation and Dielectric Phenomena，CEIDP ’09(Virginia Beach，VA，18—21 October 2009)p 592—5

[15]AgostinoR D，F(xiàn)avia P，Oehr C，Wertheimer M R 2005 Plasma Processes and Polymers 2 7

[16]Nahomy J，F(xiàn)erreira C M，Gordiets B，Pagnon D，Touzeau M，Vialle M，2010 J.Phys.D：Appl.Phys.107 093304

[17]PancheshnyiS V，StarikovskiiA Y 2003 J.Phys.D：Appl.Phys.36 268

[18]Hagelaar G J M，Pitchford L C 2005 Plasma Sources Sci.Technol.14 722

[19]Zheleznyak M D，Mnattskanyan A K 1977 ZhurnalTekhnicheskoiFiziki47 2497

[20]Yu V S，Larsson A，GubanskiS M，Akyuz M 2001 J.Phys.D：Appl.Phys.34 614

[21]Liu X H，He W，Yang F，Wang H Y，LiaoR J，XiaoH G 2012 Chin.Phys.B 21 75201

[22]Philip D N，Janzen A R，Aziz R A 1972 J.Chem.Phys.57 1100

[23]Brokaw R S 1969 Ind.Eng.Chem.Process Des.8 240

[24]Bird R B，Stewart W E，Lightfoot E N 1960 Transport Phenomena(Madison：Madison Press)

[25]Farouk T，F(xiàn)arouk B，GutsolA，F(xiàn)ridman A 2008 Plasma Sources Sci.Technol.17 035015

[26]Curtiss C F，Bird R B 1999 Ind.Eng.Chem.Res.38 2515

[27]Xu X J，Zhu D C 1996 Air discharge physical(Shanghai：Fudan University Press)(in Chinese)[徐學(xué)基，諸定昌1996空氣放電物理(上海：復(fù)旦大學(xué)出版社)]

[28]Liu X H，He W，Yang F，XiaoH G，Ma J 2011 High Voltage Engineering 37 1614(in Chinese)[劉興華，何為，楊帆，肖漢光，馬俊2011高電壓技術(shù)37 1614]

[29]Du H L，He L M，Lan Y u，Wang F 2011 Acta Phys.Sin.60 115201(in Chinese)[杜宏亮，何立明，蘭宇丹，王峰2011物理學(xué)報(bào)60 115201]

[30]He W，Liu X H，Yang F，Wang H u，LiaoR J，XiaoH G 2012 Jpn.J.Appl.Phys.51 026001

[31]Tran T N，Golosnoy I O，Lewin P L，Georghiou G E 2011 J.Phys.D：Appl.Phys.44 015203

[32]Yu D R，Qing S W，Yan G J，Duan P 2011 Chin.Phys.B 20 65204

[33]Wang P X，F(xiàn)an F G，ZirilliF，Chen J H 2012 IEEE T.Plasma Sci.40 421

[34]AntaoD S，Staack D A，F(xiàn)ridman A，F(xiàn)arouk B 2009 Plasma Sources Sci.Technol.18 035016

[35]Mahadev S，Raja L L 2010 J.Appl.Phys.107 093304