李志剛,陳宗勝
(國(guó)防科技大學(xué) 脈沖功率激光技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 安徽 合肥 230037)
等離子體是一種由大量自由運(yùn)動(dòng)的分子(或原子)、電子、離子構(gòu)成的特殊媒質(zhì),其電磁特性與內(nèi)部帶電粒子的變化密切相關(guān)[1]。當(dāng)高功率微波入射到等離子體中時(shí),兩者之間的相互作用明顯。一方面,入射波在等離子體區(qū)域產(chǎn)生強(qiáng)電磁場(chǎng),帶電粒子在強(qiáng)電磁場(chǎng)的作用下產(chǎn)生加速運(yùn)動(dòng),等離子體內(nèi)能提高,內(nèi)部碰撞反應(yīng)、非彈性碰撞加劇,等離子體狀態(tài)產(chǎn)生劇烈變化,內(nèi)部電子密度、電子能量、重粒子密度不斷增加;另一方面,等離子體的電磁特性也會(huì)發(fā)生較大改變,對(duì)入射波的損耗性能加劇,同時(shí),當(dāng)不斷增加的電子密度達(dá)到電磁截止密度時(shí),等離子體就會(huì)對(duì)入射波產(chǎn)生屏蔽作用,高功率微波無(wú)法透過(guò)等離子體區(qū)域[2-4]。
等離子體與高功率微波之間的相互作用可通過(guò)等離子體中帶電粒子的變化表現(xiàn)出來(lái),通常采用電子數(shù)密度進(jìn)行表征[5]。為深入研究等離子體與高功率微波之間的相互作用,本文采用等離子體流體近似方法,利用COMSOL軟件仿真分析研究了相互作用中的電子數(shù)密度和電子能的變化特性,著重分析了相互作用中產(chǎn)生的電子雪崩效應(yīng)的變化過(guò)程,探究等離子體產(chǎn)生電磁屏蔽的作用機(jī)理。
等離子體是由帶電粒子構(gòu)成的高度非線性組合體,想要精確描述等離子體的動(dòng)力學(xué)行為,必須要對(duì)等離子體中的每個(gè)粒子進(jìn)行軌跡跟蹤和分析,但是,由于等離子體內(nèi)部粒子數(shù)量實(shí)在太大,在實(shí)際研究過(guò)程中,這是無(wú)法實(shí)現(xiàn)的[6]。目前,只能根據(jù)不同條件和研究問(wèn)題,采用不同的近似方法,對(duì)研究問(wèn)題進(jìn)行簡(jiǎn)化分析,本文采用流體近似研究方法對(duì)相互作用中的電子分布特性進(jìn)行分析研究[7]。
流體近似方法是將等離子體等效為導(dǎo)電流體,采用經(jīng)典流體力學(xué)和電動(dòng)力學(xué)相結(jié)合的方法,研究等離子體與外界電磁場(chǎng)之間的相互作用,只用宏觀參量來(lái)描述等離子體的物理特性,著重分析等離子體的整體行為,不考慮其中單個(gè)粒子的運(yùn)動(dòng)狀態(tài)[8]。流體近似法主要用于描述等離子體的宏觀運(yùn)動(dòng),可通過(guò)求解流體力學(xué)方程組進(jìn)行求解[9]。在等離子體與高功率微波相互作用中,高功率微波的傳播過(guò)程可通過(guò)波動(dòng)方程進(jìn)行表征,等離子體內(nèi)部粒子的變化作用可通過(guò)電子漂移-擴(kuò)散方程和重物質(zhì)傳遞方程進(jìn)行計(jì)算分析。
將等離子體作為特殊電介質(zhì)進(jìn)行分析,等離子體中的麥克斯韋方程為:
▽×H=jωε0εr·E
(1)
▽×E=-jωμ0H
(2)
可得波動(dòng)方程為:
(3)
其中,E為電場(chǎng)強(qiáng)度,H為磁場(chǎng)強(qiáng)度,ω為入射波的角頻率,ε0、μ0、k0分別為真空中的介電常數(shù)、磁導(dǎo)率和波數(shù),εr為等離子體的相對(duì)介電常數(shù),可通過(guò)下式求解:
(4)
式中,ωp為等離子頻率,υ為電子碰撞頻率,i為復(fù)數(shù)標(biāo)號(hào)。
在相互作用中,電子數(shù)密度的變化可通過(guò)電子傳遞方程來(lái)進(jìn)行分析,可表示為:
(5)
Γe=-(μe·E)ne-De·▽ne
(6)
其中:ne為電子密度;μe、De分別為電子遷移率和電子擴(kuò)散率;Re為電子源項(xiàng),表征等離子體內(nèi)部碰撞反應(yīng)導(dǎo)致的電子的產(chǎn)生與消失,可通過(guò)下式求得。
(7)
式中,xj為j反應(yīng)中碰撞粒子的摩爾質(zhì)量分?jǐn)?shù),kj為j反應(yīng)的反應(yīng)速率,Nn為等離子體中總的粒子數(shù)密度。
在相互作用中,非電子粒子也會(huì)受到外加電磁場(chǎng)的影響,這些變化可通過(guò)重物質(zhì)傳遞方程來(lái)進(jìn)行分析,可表示為:
(8)
其中:ωk為第k種粒子的摩爾分?jǐn)?shù);ρ為氣體密度;μ為平均流體速度;jk為第k種粒子的擴(kuò)散通量,可通過(guò)下式表示。
jk=ρωkVk
(9)
(10)
圖1 等離子體與高功率微波相互作用計(jì)算模型Fig.1 Schematic diagram of the interactions between HPM and plasma
基于流體近似方法,本文采用COMSOL軟件對(duì)波動(dòng)方程、電子漂移-擴(kuò)散方程和重物質(zhì)傳遞方程進(jìn)行求解,對(duì)高功率微波入射到等離子體中產(chǎn)生的變化進(jìn)行分析[10]。計(jì)算模型如圖1所示,模擬分析高功率微波對(duì)等離子體管陣列進(jìn)行激勵(lì),等離子體層由左右延伸的等離子體陣列結(jié)構(gòu)組成,電磁波由上邊界入射傳播方向?yàn)?y方向,經(jīng)空氣層和等離子體層,傳輸至完美匹配層(Perfectly Matched Layer, PML)后被完全吸收。入射電磁波的功率大小采用電場(chǎng)分量來(lái)表述,假設(shè)入射電磁波場(chǎng)強(qiáng)幅值為E0,頻率為f,極化方向與等離子體管方向一致;然后,等離子體區(qū)域填充低氣壓氬氣作為等離子體源體,相互作用中的碰撞粒子和反應(yīng)在表1給出,等離子體區(qū)域氣體溫度設(shè)為T,氣體壓強(qiáng)為P。
表1 氬等離子體發(fā)生的碰撞反應(yīng)
本文重點(diǎn)研究了等離子體與高功率微波相互作用中電子分布特性產(chǎn)生的變化。初始時(shí)刻,內(nèi)部電子初始分布采用均勻分布,初始電子密度分別設(shè)置為1.0×1016m-3,初始電子能為3 eV,碰撞頻率為5×109Hz,高功率微波頻率選擇 6 GHz,電場(chǎng)強(qiáng)度為2×106V/m,氣體溫度設(shè)為300 K,壓強(qiáng)設(shè)為400 Pa,計(jì)算時(shí)間設(shè)置為0~0.01 s。
3.1.1 數(shù)值分析
圖2 等離子體與高功率微波相互作用中平均電子數(shù)密度的演變過(guò)程Fig. 2 Time dependence of the averaged electron density in the plasma of the interactions between HPM and plasma
圖3 等離子體與高功率微波相互作用中電子變化速率的演變過(guò)程Fig.3 Time dependence of electron change rate in the interactions between HPM and plasma
圖2給出了等離子體與強(qiáng)電磁脈沖相互作用中等離子體區(qū)域平均電子數(shù)密度的演化曲線,圖3為電子平均變化速率的演變過(guò)程,圖中橫縱坐標(biāo)均采用對(duì)數(shù)坐標(biāo)表示。從圖中可以看出,等離子體區(qū)域電子數(shù)密度發(fā)生劇烈變化,在計(jì)算時(shí)間(0~0.01 s)內(nèi),平均電子數(shù)密度從初始值1.0×1016m-3增長(zhǎng)至8.0×1022m-3,變化范圍近7個(gè)數(shù)量級(jí)。按照?qǐng)D中曲線變化趨勢(shì),電子數(shù)密度的變化過(guò)程在時(shí)間上可劃分為四個(gè)階段:
1)緩慢增長(zhǎng)階段(0 ~t1,這里t1=4.4×10-11s),平均電子數(shù)密度由1.0×1016m-3增加至1.1×1016m-3,電子變化速率在10 mol/(m-3·s)以下。產(chǎn)生這種現(xiàn)象的原因主要有兩個(gè)方面,一是這一階段高功率微波入射時(shí)間較短,并且等離子體內(nèi)電子數(shù)密度較小,等離子體對(duì)入射波的吸收能力較弱,因而只有較少的入射能量被等離子體吸收;二是等離子體具有很強(qiáng)的維持內(nèi)部粒子平衡的趨勢(shì),初始時(shí)刻等離子體內(nèi)部碰撞反應(yīng)以彈性碰撞(反應(yīng)1)為主,等離子體將吸收的電磁能量轉(zhuǎn)化為自身內(nèi)能,削弱外加電磁波對(duì)內(nèi)部平衡的干擾,同時(shí),吸收的電磁能量不足以激發(fā)更為頻繁的非彈性碰撞反應(yīng),新生電子產(chǎn)生速率相對(duì)較慢,和熄滅速度基本保持一致。
2)雪崩增長(zhǎng)階段(t1~t2,這里t2=1.2×10-9s),等離子體區(qū)域發(fā)生電子雪崩效應(yīng),在極短的時(shí)間內(nèi)(< 1 ns),平均電子數(shù)密度由1.1×1016m-3躍變至5.1×1018m-3,增長(zhǎng)2個(gè)數(shù)量級(jí),電子產(chǎn)生速率也由10 mol/(m-3·s)迅速增加至1.6×104mol/(m-3·s)。值得注意的是,在這一階段內(nèi),當(dāng)t=tc(這里,tc=5.0×10-10s)時(shí),平均電子數(shù)密度值為4.6×1017m-3,與入射波(6 GHz)的電磁截止密度相一致。產(chǎn)生這一現(xiàn)象的原因在于:初始階段,等離子體吸收電磁能量主要用于提高自身內(nèi)能,大量電子的運(yùn)動(dòng)速度增加,由低能電子轉(zhuǎn)變?yōu)楦吣茈娮?,具備大范圍發(fā)生非彈性碰撞反應(yīng)的能量基礎(chǔ),持續(xù)的電磁能量輸入,使得等離子體內(nèi)激發(fā)碰撞(反應(yīng)2)、電離碰撞(反應(yīng)3、反應(yīng)4)等反應(yīng)突顯、頻繁發(fā)生,電子產(chǎn)生速率遠(yuǎn)遠(yuǎn)大于電子熄滅速度,大量的電子和激發(fā)態(tài)氬原子不斷產(chǎn)生,同時(shí),新生電子和激發(fā)態(tài)氬原子也會(huì)參與到內(nèi)部非彈性碰撞反應(yīng)中,產(chǎn)生更多的電子和激發(fā)態(tài)氬原子,循環(huán)往復(fù),導(dǎo)致電子雪崩效應(yīng)產(chǎn)生,等離子體中電子數(shù)密度急劇增加,直至初始階段積累的高能電子基本消耗殆盡,電子雪崩增長(zhǎng)才會(huì)結(jié)束。
3)振蕩增長(zhǎng)階段(t2~t3,這里t3= 2.3×10-4s),平均電子數(shù)密度由5.1×1018m-3振蕩增加至7.6×1022m-3,電子產(chǎn)生速率有所減緩,由1.6×104mol/(m-3·s)逐漸減小為400 mol/(m-3·s)。這一階段電子數(shù)密度的變化較為復(fù)雜,主要是下面幾個(gè)因素導(dǎo)致的:一是在前一階段,高能電子被大量消耗,電離反應(yīng)4(電子與基態(tài)氬原子碰撞直接電離產(chǎn)生氬離子和新生電子)由于具有較高的能量壁壘(15.8 eV),反應(yīng)速率難以維持,導(dǎo)致等離子體電子產(chǎn)生速率減??;二是前一階段在產(chǎn)生大量電子的同時(shí),激發(fā)反應(yīng)也產(chǎn)生了大量的激發(fā)態(tài)氬原子,這些激發(fā)態(tài)氬原子參與到反應(yīng)5(電子與激發(fā)態(tài)氬原子碰撞電離產(chǎn)生氬離子和新生電子)中,電離反應(yīng)5由于所需能量(4.24 eV)較少,在這一階段更容易發(fā)生,這一反應(yīng)能夠在一定程度上維持電子的產(chǎn)生;三是電子數(shù)密度不斷增加,會(huì)導(dǎo)致碰撞反應(yīng)的發(fā)生概率也相應(yīng)增加,同時(shí)等離子體可以將更多的入射能量轉(zhuǎn)換為等離子體內(nèi)能,參與到內(nèi)部碰撞反應(yīng)中,為電子的產(chǎn)生持續(xù)提供能量。
4)平穩(wěn)保持階段(t3~ 0.01 s),等離子體區(qū)域平均電子數(shù)密度保持小幅增長(zhǎng),整體基本穩(wěn)定保持在8.0×1022m-3,電子產(chǎn)生速率基本維持在700 mol/(m-3·s) 左右。在這一階段,電子產(chǎn)生速率保持在較小的層次。產(chǎn)生這一現(xiàn)象的原因在于經(jīng)過(guò)前面階段的劇烈反應(yīng),高能電子和激發(fā)態(tài)氬原子被大量消耗,等離子體區(qū)域電子數(shù)密度大幅增長(zhǎng),單純的外加能量輸入難以維持頻繁、劇烈的非彈性碰撞反應(yīng),等離子體內(nèi)部逐漸達(dá)到新的平衡,內(nèi)部反應(yīng)趨于緩和,電子產(chǎn)生速率減緩。
3.1.2 空間分布
圖4給出了等離子體中的電子數(shù)密度空間分布隨作用時(shí)間的演變過(guò)程,整體來(lái)說(shuō),等離子體中電子數(shù)密度隨計(jì)算時(shí)間不斷增大。
(a) t=2.15×10-10 s (b) t=3.98×10-10 s (c) t=4.89×10-10 s
(d) t=5.99×10-10 s (e) t=7.35×10-10 s (f) t=9.03×10-10 s
(g) t=1.11×10-9 s (h) t=1.05×10-8 s (i) t=1×10-7 s
(j) t=1.1×10-5 s (k) t=1.05×10-4 s (l) t=0.001 s圖4 等離子體與高功率微波相互作用中電子數(shù)密度空間分布的演變過(guò)程Fig. 4 Time dependence of the spatial distribution of the electron number density in the interactions between HPM and plasma
初始階段(圖中t=2.15×10-10~4.89×10-10s)電子數(shù)密度變化最明顯的特征是整體均勻增長(zhǎng)。這是由于等離子體區(qū)域初始電子數(shù)密度值較小,入射波能夠在等離子體中傳輸,強(qiáng)電磁脈沖產(chǎn)生的外加電磁場(chǎng)將加速帶電粒子運(yùn)動(dòng),加劇內(nèi)部碰撞反應(yīng)的發(fā)生,尤其是非彈性碰撞反應(yīng),導(dǎo)致等離子體區(qū)域電子數(shù)密度急速增長(zhǎng)。另外,從單幅圖像對(duì)比來(lái)說(shuō),等離子體邊緣區(qū)域的電子數(shù)密度增長(zhǎng)速度要大于中心區(qū)域。產(chǎn)生這種現(xiàn)象的原因主要有兩個(gè)方面:一是等離子體具有一定的電子數(shù)密度,帶電粒子頻繁地進(jìn)行不同類型的碰撞反應(yīng),會(huì)對(duì)入射波產(chǎn)生一定的衰減,這樣導(dǎo)致入射電場(chǎng)周邊大中間??;二是柱狀結(jié)構(gòu)的等離子體類似于柱狀電介質(zhì),能夠改變?nèi)肷洳ǖ碾姶欧植?,?huì)形成環(huán)繞柱狀結(jié)構(gòu)的電磁分布,由于等離子體具有一定的衰減作用,入射到等離子體中的電磁波能量會(huì)受到衰減,圓柱周邊電磁作用強(qiáng)烈,電子碰撞反應(yīng)更為劇烈。
在電子數(shù)密度快速增長(zhǎng)階段(圖中t=5.99×10-10~1×10-7s)電子數(shù)密度峰值表現(xiàn)出向波源方向移動(dòng)的趨勢(shì),并且電子數(shù)密度峰值越來(lái)越尖銳。從t=5.99×10-10s開(kāi)始,等離子體區(qū)域下邊緣處電子數(shù)密度增長(zhǎng)減緩,與其他區(qū)域相比,電子密度值較低。這是由于等離子體兩側(cè)區(qū)域的電子數(shù)密度增加速率相對(duì)較大,電子數(shù)密度率先增長(zhǎng)至電磁截止密度。這樣,就會(huì)使得入射到等離子體下邊緣區(qū)域的電磁能量大幅度減小,電子數(shù)密度增加減緩,出現(xiàn)下邊緣區(qū)域電子數(shù)密度相比于其他區(qū)域較低的現(xiàn)象,電子數(shù)密度峰值開(kāi)始向等離子體上邊緣區(qū)域移動(dòng)。另外,隨著等離子體上邊緣區(qū)域電子數(shù)密度進(jìn)一步增加,透射到下邊緣等離子體區(qū)域的電磁能量會(huì)進(jìn)一步減小,碰撞反應(yīng)速率進(jìn)一步降低,下邊緣區(qū)域的電子數(shù)密度增長(zhǎng)進(jìn)一步減緩,等離子體區(qū)域電子數(shù)密度峰值進(jìn)一步上移。
在小幅增長(zhǎng)、穩(wěn)定保持階段(圖中t=1.10×10-5~0.001 s),等離子體區(qū)域電子分布開(kāi)始趨于均勻,電子數(shù)密度峰值降低。在前一階段中,峰值位置大量的高能電子和激發(fā)態(tài)氬原子被消耗,等離子體區(qū)域頻繁進(jìn)行的非彈性碰撞反應(yīng)趨于緩和,電子數(shù)密度增長(zhǎng)減緩,同時(shí)等離子體區(qū)域存在較大的濃度梯度,上邊緣區(qū)域的高濃度電子就會(huì)向下邊緣區(qū)域漂移擴(kuò)散,在下邊緣區(qū)域碰撞反應(yīng),等離子體區(qū)域電子密度逐漸趨于均勻,最終達(dá)到新的內(nèi)部平衡,入射波無(wú)法維持等離子體內(nèi)電子快速增長(zhǎng),也無(wú)法對(duì)新產(chǎn)生的內(nèi)部平衡產(chǎn)生破壞,等離子體吸收入射波提供的電磁能量,來(lái)維持正常的內(nèi)部反應(yīng)。
3.2.1 數(shù)值分析
圖5給出了等離子體區(qū)域平均電子能量隨時(shí)間的演變曲線。從圖中可以看出:在0 ~t1時(shí)間內(nèi)(這里,t1=5.0×10-10s),平均電子能由3 eV快速增加到112 eV,在t1~t2時(shí)間內(nèi)(這里,t2=1.2×10-9s),平均電子能急劇減小,由112 eV減小至8 eV,在t2~t3時(shí)間內(nèi)(這里,t3=2.3×10-4s),平均電子振蕩減小至2 eV,而后保持穩(wěn)定。值得注意的是,平均電子能達(dá)到峰值點(diǎn)的時(shí)間與電子數(shù)密度達(dá)到截止密度的時(shí)間相一致;平均電子能急速衰減結(jié)束的時(shí)間與電子數(shù)密度雪崩增長(zhǎng)階段的結(jié)束時(shí)間相一致。這一現(xiàn)象說(shuō)明,直接電離碰撞是電子雪崩效應(yīng)產(chǎn)生的主要反應(yīng)。
圖5 等離子體與強(qiáng)電磁脈沖相互作用平均電子能的演變過(guò)程Fig. 5 Time dependence of the average electron energy in the interactions between HPM and plasma
產(chǎn)生這一變化的原因?yàn)椋撼跏茧A段,入射波會(huì)在等離子體區(qū)域產(chǎn)生很強(qiáng)的電磁場(chǎng),內(nèi)部電子受到外加電磁場(chǎng)的作用產(chǎn)生加速運(yùn)動(dòng),電子內(nèi)能迅速增加,由低能電子轉(zhuǎn)化為高能電子,平均電子能增加迅速。而后,等離子體內(nèi)碰撞反應(yīng)開(kāi)始加劇,電子數(shù)密度也開(kāi)始急劇增加,高能電子開(kāi)始被消耗,當(dāng)?shù)入x子體區(qū)域平均電子數(shù)密度達(dá)到截止密度時(shí),等離子體開(kāi)始對(duì)入射波產(chǎn)生屏蔽作用,電磁能量傳輸受限,平均電子能增加受到限制;與此同時(shí),前期產(chǎn)生的高能電子會(huì)在劇烈的非彈性碰撞中被大量消耗,等離子體區(qū)域平均電子能急劇下降。當(dāng)內(nèi)部碰撞反應(yīng)趨于緩和后,高能電子消耗減小,等離子體區(qū)域逐漸建立新的粒子平衡,電子數(shù)密度增加減緩,電子分布開(kāi)始趨于均勻,平均電子能也逐漸趨于穩(wěn)定。
3.2.2 空間分布
圖6給出了等離子體區(qū)域電子能的空間分布隨作用時(shí)間的演變過(guò)程。從圖中可以看出,初始階段,等離子體區(qū)域電子能表現(xiàn)為整體增長(zhǎng),由于外加電磁場(chǎng)分布的影響,電子能會(huì)出現(xiàn)邊緣大中間小的現(xiàn)象,同時(shí),隨著電子數(shù)密度不斷增加,等離子體上游區(qū)域的電子能增加速度要快于下游區(qū)域,電子能峰值向波源方向移動(dòng),開(kāi)始時(shí)會(huì)出現(xiàn)兩個(gè)峰值點(diǎn),這是計(jì)算模型的對(duì)稱結(jié)構(gòu)產(chǎn)生的。
(a) t=2.15×10-10 s (b) t=3.98×10-10 s (c) t=4.89×10-10 s
(d) t=5.99×10-10 s (e) t=7.35×10-10 s (f) t=9.03×10-10 s
(g) t=1.11×10-9 s (h) t=1.05×10-8 s (i) t=1×10-7 s
(j) t=1.1×10-5 s (k) t=1.05×10-4 s (l) t=0.001 s圖6 等離子體與強(qiáng)電磁脈沖相互作用中電子能空間分布的演變過(guò)程Fig.6 Time dependence of the spatial distribution of the electron energy in the interactions between HPM and plasma
隨著作用時(shí)間的增加,電子能峰值繼續(xù)朝向等離子體上邊緣區(qū)域移動(dòng),兩峰值點(diǎn)也逐漸匯合,在等離子體上邊緣區(qū)域產(chǎn)生單一的峰值點(diǎn)。而后,等離子體中雪崩效應(yīng)產(chǎn)生,大量的高能電子被消耗用于產(chǎn)生新的粒子,等離子體區(qū)域的電子能急速下降,但峰值位置仍在等離子體上邊緣區(qū)域。最后,等離子體內(nèi)部非彈性碰撞反應(yīng)變緩,等離子體內(nèi)部達(dá)到新的平衡,電子數(shù)密度峰值開(kāi)始減小,等離子體區(qū)域平均電子能開(kāi)始變得均勻。
電子數(shù)密度是表征等離子體物理特性的一項(xiàng)重要因素,當(dāng)高功率微波入射到等離子體中時(shí),等離子體區(qū)域的電子分布特性會(huì)產(chǎn)生劇烈變化。本文采用COMSOL軟件仿真分析了相互作用中等離子體區(qū)域電子分布特性,計(jì)算結(jié)果表明,等離子體內(nèi)部激烈的碰撞反應(yīng)是導(dǎo)致電子分布特性變化的重要原因,在高功率微波作用下,等離子體區(qū)域電子數(shù)密度在數(shù)值上會(huì)產(chǎn)生顯著的雪崩增長(zhǎng)變化,空間分布上電子數(shù)密度峰值趨于入射方向移動(dòng);電子能隨入射激勵(lì)增加呈增長(zhǎng)趨勢(shì),隨電子數(shù)密度的增加而減小。
國(guó)防科技大學(xué)學(xué)報(bào)2020年1期