譚康伯 路宏敏 蘇濤

(西安電子科技大學,天線與微波國家重點實驗室,西安 710071)

基于系統(tǒng)性電磁兼容的考慮,對等離子環(huán)境中在軌航天器導體充放電現(xiàn)象中的能量特征進行了變分研究.通過電磁Collin原理,對等離子環(huán)境中導體系統(tǒng)幾何尺度與所帶電能的變分聯(lián)系進行了理論分析.在此基礎上,推廣了更具一般性的數(shù)值估值分析方法,并對復雜導體系統(tǒng)電磁參數(shù)、等離子環(huán)境特征與系統(tǒng)能量間的關系進行了實例分析.研究結(jié)果對于等離子環(huán)境中復雜帶電體的能量控制及相關的電磁環(huán)境效應與防護等研究具有積極意義.

1 引 言

等離子是物質(zhì)存在的一種特殊形態(tài)[1?7],對處于其中的導體系統(tǒng)形成影響.在等離子環(huán)境工程中,導體充放電現(xiàn)象將導致在軌故障,對在軌航天器的有效運行以及功能發(fā)揮構成嚴重威脅[8?16].等離子環(huán)境中航天器的帶電分析是系統(tǒng)性電磁兼容領域中需要持續(xù)研究的重點問題.

近年來,針對空間系統(tǒng)的靜電防護,電磁結(jié)合統(tǒng)計的等離子環(huán)境的物理模擬為相關研究的不斷深入提供了重要的基礎支撐.從動力學角度來看,等離子環(huán)境中的導體系統(tǒng)會受到來自周圍大量粒子的轟擊,進而在系統(tǒng)表面形成電荷積累[10?15].這就構成了由導體與其周圍等離子環(huán)境相關聯(lián)的特殊動力學共同體,并呈現(xiàn)特定的能量規(guī)律.

對于動力學系統(tǒng)能量特征的探討,變分方法原理簡潔[17?21],在分析力學中有著重要的物理意義,可以深入刻畫相關動力學的本質(zhì)規(guī)律.在電磁研究中,變分方法對于系統(tǒng)特征的評估具有應用價值.因此,等離子環(huán)境中帶電體能量的變分研究成為相關系統(tǒng)性電磁兼容理論研究中必須進一步深入探討的問題.

本文基于變分原理,通過電磁Collin分析對等離子環(huán)境中與導體目標相關的能量估計問題進行討論.

2 等離子中導體系統(tǒng)的Collin變分原理

在軌航天器與其周圍等離子環(huán)境構成了特殊的動力學系統(tǒng),其導體系統(tǒng)能量受兩方面共同影響.一方面,由于等離子環(huán)境中具有一定能量的電子與離子隨機運動[8?15],在導體表面上將產(chǎn)生帶電粒子的聚集作用.在高溫環(huán)境中,電子和離子有近似相同的熱力學溫度,即Te≈Ti.電子質(zhì)量me遠小于離子質(zhì)量mi,處于平衡態(tài)時電子速度遠大于離子速度,即ve?vi.由統(tǒng)計理論可知,到達導體目標表面的電子數(shù)目將大于離子數(shù)目.于是,導體表面形成大量負電荷聚集,構成負電能積累.另一方面,由電磁理論可知[22,23],導體的靜電能量為Ee=Q2/(2C),其中C是導體結(jié)構的電容,為導體系統(tǒng)自身屬性,導體電位U及電量Q的變化受其約束.當導體處于電位穩(wěn)定狀態(tài)后,其靜電特征滿足關系:Q=CU.將該關系與導體系統(tǒng)能量關系共同考慮,可以得到

從(1)式可以看出,當導體電位確定時,孤立導體電容的變化與系統(tǒng)靜電能量的變化間形成約束聯(lián)系.

為了確定導體電位,需要借助于對電子和離子的聚集作用的分析.當導體系統(tǒng)處于等離子環(huán)境中時,由于導體表面負場阻礙作用的影響,電子只有沿表面法向具有大到足夠克服相應負勢能的能量時,才能到達導體表面.于是有能量關系:?eU=me⊥min/2,其中e為電子電量,ve⊥min為電子到達導體所具有的最小速率.在此基礎上,依據(jù)等離子體中粒子的運動特性,從統(tǒng)計分析及麥克斯韋分布出發(fā),可得電子流密度Γe為

其中k為玻爾茲曼常數(shù),ne和ni為電子和離子密度.在粒子流作用下,等離子環(huán)境中孤立導體的電位平衡特征可由理論確定[13?15].負場的增強對電子向?qū)w表面聚集形成阻礙作用,進而使流向?qū)w表面的電子流密度逐漸減少,其中二次電子效應對表面電子數(shù)目也將產(chǎn)生影響.當?shù)竭_孤立導體表面的電子流、離子流以及二次電子流作用平衡時,表面電位將達到穩(wěn)定狀態(tài).此時有流平衡關系:Γe= Γi+Γse,其中,Γse為二次電子流密度.Γse= δeΓe+ δiΓi,其中δe為目標表面材料由電子引起的二次電子系數(shù),δi則為由離子引起的二次電子系數(shù).結(jié)合(2)式、(3)式,將其代入流平衡關系,可得(1)式中所需的平衡電位U.

由上面分析可知,在等離子環(huán)境中的帶電導體系統(tǒng)呈現(xiàn)一定的能量特征,對其進行研究需要考慮電容的實際情況.Collin原理是變分在電磁研究中的重要延拓,可以對導體系統(tǒng)電容進行實用分析[24,25].為了探討問題的內(nèi)在聯(lián)系,本文引入變分方法進行研究.由該原理可知,較大的孤立導體有著較大的電容.

圖1中所示的導體結(jié)構S1?S??S2,C1

在相同電位和相同平均電荷密度條件下,該形式表明相同的表面積有相同的總電荷,于是對應相同的電容.在此理論原理基礎上,當電位得到一致性的確定后,由幾何尺度所確定的電容將形成對于平衡態(tài)等離子系統(tǒng)中靜電能的直接變分估計:

圖1 Collin原理圖Fig.1.Collin’s principle.

3 實例分析和一般性數(shù)值推廣

首先通過簡單實例,對Collin變分分析的應用進行探討.考慮空間尺度為a×a×a的方形導體目標.根據(jù)上一節(jié)中的變分條件,對該目標做相同表面積估值,所對應的等效半徑為則可得到其相應的電容估值為

將(6)式代入(1)式,并結(jié)合所得平衡電位,可得平衡態(tài)下孤立方形導體電能的理論估值√為

從上面的結(jié)果可以看到,對于方形目標,電能與幾何尺度在變分條件下成線性比例關系.該結(jié)果和文獻[15]所得的電能與面積關系在物理上保持一致.對于電子所引起的二次效應,類似于[14,15,26,27],取δ= δi≈ δe且ni≈ ne,可得到二次電子效應對于帶電體能量的影響,其結(jié)果如圖2所示.顯然,上面方法在對一些相對簡單結(jié)構的電能特征的分析應用方面具有理論優(yōu)勢.

作為一種較好的實用分析方法,Collin變分原理為處理等離子環(huán)境中更復雜的導體系統(tǒng)提供了理論工具.為了有利于更普遍的應用,以面對復雜結(jié)構所提出的一般性分析要求,下面將基于變分原理(5)式,同時考慮等離子平衡態(tài)分析理論中的統(tǒng)計平均性,引入數(shù)值離散技術,對分析理論方法做進一步的推廣.

圖2 二次電子效應對等離子環(huán)境中帶電立方體的影響Fig.2.Ef f ect of secondary electron emission on cube in plasma environment.

在等離子環(huán)境中,電子和離子流作用于導體系統(tǒng)表面的法向方向,形成電荷及能量積累.考慮到疊加特性,可將目標表面的待求區(qū)域進行N次離散剖分,以離散面片近似積分微元.通過離散,使得原有導體目標的幾何構型S?近似為離散剖分幾何構型S,其過程如圖3所示.平衡態(tài)時,設ρ表示目標上的電荷密度函數(shù),則由空間電位關系可得如下積分算子方程:

其中L表示靜場積分算子,R=|r?r′|,待求的未知函數(shù)是產(chǎn)生目標電位的電荷密度ρ.

圖3 導體構型的空間離散Fig.3.Spatial discreteness of conductor con fi guration.

電子和離子流作用于各導體面片法向.根據(jù)幾何構型進行離散,采用基函數(shù)fn將待求電荷密度函數(shù)展開,有從該式可以看出,只要求得系數(shù)an,就可進一步得到電荷密∫ 度 的分布情況.定義該問題的內(nèi)積為將前面電荷密度的展開式代入內(nèi)積表達式中,并且依據(jù)典型的矩陣向量形式[28],取檢驗采樣函數(shù)為wm(r?rm),其中rm為每個離散面片的中點,以此標定法向流作用的空間特征.至此,(8)式轉(zhuǎn)化為相應的矩陣代數(shù)方程:

其中l(wèi)m,n= ?wm,Lfn?,bm?wm,U?. 通過求解上面矩陣代數(shù)方程,得到平衡態(tài)電荷密度分布函數(shù)的系數(shù)[an]=[lm,n]?1[bm],進而得到電荷密度函數(shù).于是,在前面數(shù)值離散分析的基礎上,即可得到等離子環(huán)境中相應導體目標的平衡態(tài)系統(tǒng)的能量估值:

通過理論分析可以看到,(10)式是(1)式在Collin原理下的離散近似,其結(jié)果與實際目標的離散度相關.

下面將基于Collin變分的離散方法應用于具體的等離子環(huán)境中導體目標的靜電能量分析中,并進行討論.分析中計及電子所引起的二次效應,并取T=Te≈Ti.首先,考慮處于等離子環(huán)境中a=0.5 m的立方導體,通過(7)式和(10)式,可以分別得到其靜電能量理論分析和數(shù)值分析估值(圖4).

圖4 等離子環(huán)境中方形導體能量特征的估值對比Fig.4.Comparison for estimated energy of cube in plasma environment.

由上面結(jié)果可以看到,對于該方形導體目標,理論與數(shù)值分析具有一定的符合度.在等離子環(huán)境中,較高的等離子體溫度對應于更為激烈的粒子運動,于是在粒子流作用下,導體結(jié)構上積累了更多能量.其次,考慮到數(shù)值途徑為研究更復雜的導體系統(tǒng)提供了便利,下面給出具體實例分析.考慮半徑為0.25 m球體子部分和邊長為0.5 m立方體子部分的組合結(jié)構,其中兩個子部分的體心間距為0.43 m,通過分析可以得到該帶電系統(tǒng)在等離子環(huán)境中的能量特性(圖5).由結(jié)果可知,在相同的等離子溫度時,該組合系統(tǒng)較之單一立方體的幾何形體更大,自身的能量更高.由于該組合結(jié)構較之方形于幾何構型上的變化,在等離子溫度為2 keV時,系統(tǒng)所帶能量提升了約15.3%.通過分析可以看到,數(shù)值估值能夠較為普遍地處理等離子環(huán)境中的復雜導體系統(tǒng),在應用方面具有一定的一般性.值得注意的是,對于復雜系統(tǒng),在研究其能量特性時,需具體考慮其整體幾何構型特征.

圖5 等離子環(huán)境中復雜系統(tǒng)的數(shù)值估值Fig.5.Numerical estimation of the composited system in plasma environment.

4 結(jié) 論

本文基于變分原理,探討了導體目標幾何構型與電能間的電磁Collin聯(lián)系,對等離子環(huán)境中導體的電能進行了較為簡潔的變分估值分析.出于Collin原理在處理復雜導體系統(tǒng)時實用性的考慮,在理論研究的基礎上,本文結(jié)合數(shù)值分析技術,為等離子環(huán)境中導體系統(tǒng)電能的估值提供了更具實用價值的一般性分析方法.本文研究有利于等離子環(huán)境中復雜帶電體的能量控制及相關的電磁環(huán)境效應與防護技術的深入.