馮 娜，李得天，楊生勝，陳益峰，趙呈選，湯道坦

（蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室，蘭州 730000）

0 引言

近年來，電推進(jìn)系統(tǒng)因?yàn)槠涓咝?、高比沖，可以大幅提升衛(wèi)星平臺的有效載荷承載能力，延長衛(wèi)星的在軌使用壽命，降低發(fā)射成本成為未來長壽命、大功率地球同步軌道衛(wèi)星必將采用的技術(shù)[1-2]。與傳統(tǒng)化學(xué)推進(jìn)相比，電推進(jìn)羽流中的噴出物質(zhì)是高度電離的低能高密度的等離子體，這些帶電離子在運(yùn)動過程中受到航天器所處空間背景電磁環(huán)境及航天器自身電磁場的影響，又對環(huán)境電磁場產(chǎn)生影響，產(chǎn)生化學(xué)推進(jìn)所沒有的電磁環(huán)境和等離子體環(huán)境。而電推進(jìn)產(chǎn)生的等離子體又與地球同步軌道空間天然等離子體特性存在較大差異，導(dǎo)致衛(wèi)星表面帶電過程更為復(fù)雜。在航天器局域空間內(nèi)，電推進(jìn)羽流等離子體、空間等離子體和航天器三者間相互作用產(chǎn)生多種效應(yīng)，如干擾航天器表面帶電、降低太陽電池陣的輸出功率，影響衛(wèi)星通信信號等[6-7]。因此，需深入了解電推進(jìn)等離子體羽流與航天器表面帶電的影響，通過相關(guān)技術(shù)手段對電推進(jìn)系統(tǒng)進(jìn)行風(fēng)險評估，進(jìn)而采取相應(yīng)的防護(hù)措施是其在應(yīng)用衛(wèi)星上推廣的前提。

未來通信衛(wèi)星平臺均可能采用離子推進(jìn)技術(shù)，將不可避免地會帶來電推進(jìn)等離子體的帶電效應(yīng)。為此，主要介紹應(yīng)用電推進(jìn)系統(tǒng)后在航天器周圍產(chǎn)生的等離子體和電磁場等特殊環(huán)境，并探討這些特殊環(huán)境對航天器充放電效應(yīng)的影響，重點(diǎn)從理論研究、仿真模擬和地面試驗(yàn)三個方面對電推進(jìn)系統(tǒng)對航天器表面帶電的影響研究進(jìn)展進(jìn)行分析總結(jié)。最后，對如何有效分析電推進(jìn)等離子體羽流特性及其對航天器表面帶電影響的研究途徑進(jìn)行了展望。

1 電推進(jìn)等離子體對航天器的影響機(jī)理

地球同步軌道空間等離子體引起的衛(wèi)星表面充/放電效應(yīng)，對衛(wèi)星在軌正常工作和安全運(yùn)行具有重要影響。帶電效應(yīng)會導(dǎo)致電勢變化，同時會導(dǎo)致在航天器表面之間或者航天器表面與航天器結(jié)構(gòu)之間產(chǎn)生強(qiáng)電場，大于擊穿閾值時就會發(fā)生靜電放電現(xiàn)象。電推進(jìn)產(chǎn)生的電子擴(kuò)散速度快，離子擴(kuò)散慢，從而在電推進(jìn)產(chǎn)生的等離子體區(qū)域內(nèi)形成電子密度和離子密度不相等的現(xiàn)象，且由離子和電子組成了屏蔽層，產(chǎn)生電場E，該電場的形成將有助于GEO天然等離子體中的電子進(jìn)入衛(wèi)星表面區(qū)域，同時抑制離子進(jìn)入該區(qū)域，而該電場對電推進(jìn)等離子體中的電子和離子具有相反的作用，從而影響航天器帶電環(huán)境，改變原本的充放電過程，從而影響衛(wèi)星表面充電平衡及表面電位[3]，如圖1所示。

圖1 電推進(jìn)等離子體擴(kuò)散過程示意圖

其次，電推進(jìn)系統(tǒng)產(chǎn)生的等離子體將以羽流的形式直接作用于衛(wèi)星表面，從而改變航天器表面電位。在GEO軌道等離子體作用下，衛(wèi)星表面通常帶較高的負(fù)電位，航天器帶電表面將從周圍濃稠等離子體中收集離子產(chǎn)生電流，大多數(shù)電荷交換（CEX）離子彌散于電推進(jìn)出口下游附近，受電推進(jìn)出口平面附近的電場作用，CEX離子會返流到航天器的帶電表面，如大面積的高壓偏置太陽帆板[4]。由于太陽帆板面積較大，不同電性能材料間的復(fù)合復(fù)雜，結(jié)構(gòu)精細(xì)，而電推進(jìn)等離子體的分布范圍較廣，離子數(shù)密度相差較大，電推進(jìn)等離子體與電池陣相互作用產(chǎn)生的寄生電流可達(dá)毫安量級，此效應(yīng)將改變太陽電池陣的有效工作點(diǎn)，并降低有效輸出功率[5-6]，如圖2所示。

圖2 電推進(jìn)羽流與航天器的相互作用示意圖

洛克希德·馬丁公司對多顆地球同步軌道衛(wèi)星的電位監(jiān)測顯示[7]：電推進(jìn)器羽流可降低航天器懸浮電位，但較低的航天器懸浮電位使得航天器表面濺射效應(yīng)增大，此效應(yīng)可能改變材料性能以及濺射腐蝕效應(yīng)，從而導(dǎo)致材料性能永久性退化?？偠灾?，航天器的懸浮電位需同時考慮航天器表面電流收集效應(yīng)以及材料濺射效應(yīng)的綜合作用[8]。

2 電推進(jìn)羽流對航天器表面帶電影響的國外研究

評估電推進(jìn)羽流對航天器的影響是一項(xiàng)復(fù)雜的工程，對于空間等離子體引起的航天器表面帶電的研究始于20世紀(jì)60年代，國外對于航天器表面充放電的物理機(jī)制進(jìn)行了大量的研究，開展了大量的地面及在軌飛行試驗(yàn)，研發(fā)了相應(yīng)的航天器表面充電效應(yīng)模擬軟件，且對于電推進(jìn)等離子體對于航天器的帶電影響也開展了一定的工作，以滿足電推進(jìn)系統(tǒng)與航天器的相容性和可靠性等要求。

2.1 理論分析研究

以霍爾推力器SPT100為例分析電推進(jìn)系統(tǒng)對航天器充放電效應(yīng)的影響，如圖3所示。

電推進(jìn)系統(tǒng)工作時，由空心陰極發(fā)射的電子在正交電磁場的作用下與從陽極氣體分配器出來的Xe原子產(chǎn)生碰撞產(chǎn)生帶正電的Xe+和Xe2+，Xe+在電場的作用下高速噴出產(chǎn)生推力，同時從空心陰極發(fā)射的另一部分電子與噴出的Xe+中和，使得束流保持電中性。實(shí)驗(yàn)測得距離電推進(jìn)出口30 cm處的等離子體的數(shù)密度約為1017m-3，而由高能Xe+和Xe原子交換電荷碰撞產(chǎn)生的低能CEX離子的密度相對較低（約為1012m-3），易擴(kuò)散至電推進(jìn)乃至整個航天器周圍。電推進(jìn)Xe離子能量為300 eV左右，電子能量約為4 eV，而電子的質(zhì)量約為Xe離子的1/200 000，分析可知電推進(jìn)產(chǎn)生等離子體中電子的擴(kuò)散速度是Xe離子的40倍。因此在電推進(jìn)產(chǎn)生等離子體的擴(kuò)散過程中，將形成如圖4所示的電子和離子分布[3]。

圖3 SPT100構(gòu)造示意圖

圖4 電推進(jìn)等離子體離子-電子層形成過程

圖4中的a區(qū)域?yàn)殡娡七M(jìn)產(chǎn)生等離子體，c區(qū)域?yàn)镚EO軌道空間等離子體，b區(qū)域的形成主要是由于電推進(jìn)產(chǎn)生等離子體中的電子擴(kuò)散速度遠(yuǎn)大于Xe離子速度，導(dǎo)致該區(qū)域分別產(chǎn)生Xe離子層（正電荷）和電子層（負(fù)電荷），從而形成了一個由正電荷指向負(fù)電荷的電場。而電子與離子在此電場作用下的運(yùn)動機(jī)理類似于PN結(jié)，高密度等離子體區(qū)域相當(dāng)于N結(jié)，低密度區(qū)域?yàn)镻結(jié)。采用半導(dǎo)體載流子擴(kuò)散理論計算電推進(jìn)等離子體和空間天然等離子體之間的電流，利用漂移-擴(kuò)散方程計算發(fā)現(xiàn)電推進(jìn)等離子體擴(kuò)散過程中產(chǎn)生的勢壘為20～40 V。

研究表明，作為N結(jié)多數(shù)載流子的電子擴(kuò)散至其作為少數(shù)載流子的P結(jié)形成了電場Eint，而電推進(jìn)等離子體電勢Vps與空間等離子體電勢VpGEO差形成了外部電場Eext。對于帶負(fù)電的航天器而言，即自身的電勢低于空間等離子體電勢時，電推進(jìn)產(chǎn)生的電子由N結(jié)擴(kuò)散至P結(jié)，則外電場Eext會削減內(nèi)電場Eint，即由航天器擴(kuò)散至GEO空間等離子體環(huán)境的電子流則會將降低航天器的表面電位，如圖5所示。

圖5 電推進(jìn)等離子體形成的電場分布示意圖

當(dāng)航天器帶正電時，受到外電場的影響，電場Eint與屏蔽勢壘均會增大，電推進(jìn)羽流中電子能量較低而無法越過此勢壘，則此擴(kuò)散電流則可忽略。而Xe離子則可能流向空間等離子體區(qū)域，且GEO軌道電子受到電場的作用而向航天器漂移。但由于電推進(jìn)離子的運(yùn)動速度較慢，而空間電子流較小，為此航天器結(jié)構(gòu)電位的中和速度較慢。

航天器表面受到空間等離子體中帶電粒子及光子的撞擊形成了進(jìn)入航天器表面的電流，以及環(huán)境電子產(chǎn)生的二次電子電流Ise、背散射電子電流Isc、環(huán)境離子所產(chǎn)生的二次電子電流Isi[9]和光電子Ip[10]，并且綜合考慮電推進(jìn)等離子體中氙離子流Ixe以及中和器電子電流Ine，因此航天器任意表面上的總電流密度為：

式中：Φs為航天器表面電位，當(dāng)I（Φs）=0時所對應(yīng)的Φs即為航天器表面達(dá)到電流平衡時所具有的帶電電位。

當(dāng)電推進(jìn)SPT100不工作的情況下，環(huán)境電子的隨機(jī)通量（電子熱運(yùn)動所產(chǎn)生的通量）遠(yuǎn)大于離子隨機(jī)通量。因此，航天器經(jīng)常在月食期達(dá)到電流平衡時其表面帶負(fù)電[11-12]。國外同步軌道衛(wèi)星SCATHA上數(shù)據(jù)顯示平均入射到航天器表面的電子密度為1.0×109cm-3，且能量為2.49 keV[13]，最終在相對平靜期衛(wèi)星表面電位高至-4 000 V。當(dāng)電推進(jìn)SPT100工作時，由于Xe+帶正電，其對航天器表面的充電電流與Ie方向相反，當(dāng)SPT100的中和器與衛(wèi)星地相連時，中和器噴射的中和電子引起的電流Ine也與Ie方向相反。分析看出，空間電子電流Ie為航天器表面高負(fù)電位充電的主導(dǎo)因素，因此，與Ie反方向的Ine和Ine與對航天器表面充電有一定的緩解作用。另一方面，SPT100噴出的中和電子能量很低，極易受到航天器表面電勢的屏蔽，對航天器表面充電的貢獻(xiàn)非常有限，不會使航天器充電到很高的負(fù)電位。

2.2 數(shù)值仿真模擬

自從發(fā)現(xiàn)衛(wèi)星表面帶電引起衛(wèi)星故障以后，整星表面帶電分析技術(shù)就一直是國際宇航界研究的一項(xiàng)重要內(nèi)容，而數(shù)值模擬方法可以在航天器設(shè)計階段評估航天器的在軌充電效應(yīng)，為航天器防護(hù)設(shè)計和地面測試提供參考，并給出特定航天器失效是否由表面充放電造成的判斷依據(jù)。電推進(jìn)羽流在真空中的流動問題與傳統(tǒng)化學(xué)推進(jìn)羽流不同，是包含大量帶電粒子的等離子體。欲處理帶電離子在電磁場作用下的運(yùn)動，須引入PIC（Particle-In-Cell）方法求解羽流中自洽電磁場以及粒子的運(yùn)動，并引入MCC（Monte Carlo Collision）方法處理粒子碰撞問題[14]，其計算原理如圖6所示。

圖6 PIC-MCC/DSMC計算原理框圖

表面帶電是由于等離子體與衛(wèi)星表面材料相互作用，產(chǎn)生的凈電荷在衛(wèi)星表面積累的過程。其涉及了空間等離子體運(yùn)動、材料光電子和二次電子發(fā)射、材料的泄露電流等多種因素。如要獲得電推進(jìn)等離子體環(huán)境下航天器表面帶電的變化，需求解場和粒子之間的強(qiáng)耦合作用過程。推進(jìn)器等離子體具有與GEO軌道等離子體完全不同的分布特點(diǎn)，為此利用PIC方法進(jìn)行等離子體模擬時，需兼容兩類離子的分布特性，綜合考慮進(jìn)行空間網(wǎng)格劃分，選取不當(dāng)將會直接影響到模擬結(jié)果的正確性。

國際上已經(jīng)建有許多成熟的表面帶電分析軟件，如以NASCAP為代表的航天器帶電分析軟件的開發(fā)應(yīng)用。絕大部分軟件對于航天器表面帶電的求解大部分使用PIC方法來模擬航天器外的等離子體環(huán)境，且使用電流平衡方法來計算航天器表面的帶電情況，使國外已掌握整星的表面帶電分析仿真的技術(shù)，同時也節(jié)約試驗(yàn)費(fèi)用并加快了衛(wèi)星的研制過程。

美國航空航天局（NASA）與美國空軍（US Air Force）共同開發(fā)的NASCAP（NASA Charging Analyz?er Program）是最早的航天器帶電分析軟件[15]。是以等離子體的磁流體力學(xué)分析方法為理論基礎(chǔ)的模擬軟件，可用于航天器三維表面帶電分析[16]。NASCAP-2K采用了PIC方法模擬航天器周圍的等離子體環(huán)境，使用單麥克斯韋分布和雙麥克斯韋分布來模擬等離子體。在實(shí)際計算過程中，利用科學(xué)應(yīng)用國際公司（Science Applications International Corpora?tion，SAIC）開發(fā)的交互式仿真計算軟件“電推進(jìn)相互作用程序包（Electric Propulsion Interactions Code，EPIC）[17]可以在軸對稱條件下得到電推進(jìn)羽流的數(shù)值模型。隨之將得到的電推進(jìn)羽流模型及離子推力器作為一個特殊對象代入NASCAP的航天器構(gòu)型進(jìn)行自洽的等離子體和電磁場計算。電推進(jìn)羽流中的高能離子由EPIC羽流模型給出，結(jié)合航天器表面電位和鞘層，利用SOAP（Simple Object Access Protocol）結(jié)合NASCAP-2K來計算電推進(jìn)羽流與航天器對航天器表面電位的影響。對美國的NEXT推力器進(jìn)行的模擬中，假設(shè)航天器初始電位為-100 V，研究電推進(jìn)等離子體流與航天器表面電位之間的關(guān)系，且NASCAP得到的CEX離子回流電流與實(shí)際試驗(yàn)結(jié)果相當(dāng)一致，結(jié)果如圖7、8所示。

圖7 航天器表面電位與電推進(jìn)離子電流隨時間的變化趨勢

美國的DAWN航天器上搭載了三部NSTAR離子推進(jìn)器，美國維吉尼亞理工學(xué)院采用COLISEUM[18]框架結(jié)合NASCAP-2K來模擬了NSTAR推力器CEX離子所形成的返流對太陽電池帆板的作用。由于CEX離子分布的發(fā)散性，且受到太陽電池板自身電勢的影響，不同密度的CEX離子撞擊在太陽電池帆板表面形成一個橫向的電勢梯度。

圖8 NEXT推力器CEX離子回流電流的模擬與實(shí)驗(yàn)值

歐洲航天局（ESA）利用航天器表面充電效應(yīng)的開源軟件 SPIS（Spacecraft Plasma Interactions Sys?tem），基于PIC方法模擬等離子體，并利用Vlasov-Poisson方程解算空間等離子體和電磁場分布。該軟件包含了龐大的航天材料庫，可提供電導(dǎo)率、二次發(fā)射率、場發(fā)射率及濺射系數(shù)等各類材料特性參數(shù)，模擬飛行器在等離子體環(huán)境影響的研究中，等離子體環(huán)境的研究，太陽電池陣在等離子體環(huán)境中相互影響的研究以及等離子體精確計算等方面具有很高的應(yīng)用前景[19]。為了適應(yīng)電推進(jìn)航天器表面充電計算的需要，CNES開發(fā)了電推進(jìn)羽流模塊，利用電推進(jìn)羽流的出流速度及羽流中的離子和中性粒子密度，計算電荷交換發(fā)生率，得到CEX離子的空間分布特征，最后與背景等離子體共同參與軟件的模擬計算過程。

奧地利國家技術(shù)研究院（Austrian Institute of Technology，AIT）的Tajmar等[6]利用等離子體模擬軟件SmartPIC模擬了SMART-1上搭載的PPS1350電推進(jìn)的整個羽流場的分布、航天器在電推進(jìn)等離子體下懸浮電位的變化趨勢以及交換工質(zhì)離子所形成的返流對太陽帆板的作用，搭載于SMART-1上的等離子體診斷包（Electric Propulsion Diagnostic Package，EPDP）所測得的在軌數(shù)據(jù)與SmartPIC模擬仿真結(jié)果對比，數(shù)據(jù)吻合良好[20]。

2.3 地面模擬實(shí)驗(yàn)

地面模擬實(shí)驗(yàn)一般在真空室內(nèi)進(jìn)行，真空室本底氣壓一般控制在1×10-4Pa以下，模擬電推進(jìn)羽流等離子體參數(shù)主要包括羽流的粒子成分、等離子體密度/溫度、中性氣體溫度、離子速度等。等離子體參數(shù)診斷方法一般有Langmuir探針、遲滯電位分析儀、發(fā)射光譜法等。如常采用的Langmuir探針實(shí)際上是一根金屬裸絲，探針本身需施加一個相對于等離子體偏置的電壓，并與某參考電極構(gòu)成回路。探針?biāo)┘拥钠秒妷号c回路電流構(gòu)成探針伏-安特性曲線，再依照理論假設(shè)對曲線進(jìn)行處理，可獲得探針附近的局域等離子體參數(shù)。但采用Langmuir探針作為等離子體診斷設(shè)備可能對羽流產(chǎn)生干擾。相對的發(fā)射光譜法作為一種非介入式診斷方法，其干擾較小，即利用發(fā)射光譜儀收集羽流等離子體的譜線，通過對線狀譜的分析，可以獲得羽流成分和相對含量，譜線的展寬可用于計算電子密度，再利用碰撞輻射模型和譜線的相對強(qiáng)度計算等離子體電子溫度。

利用TS-6電推進(jìn)地面實(shí)驗(yàn)設(shè)備模擬測試20 cm氙離子推力器工作時，衛(wèi)星表面相對于周圍等離子體（假設(shè)周圍等離子體電勢為參考零電勢）的電位。工作時其真空度約為6×10-3Pa，且認(rèn)為真空艙殼體與束流等離子體等電位，電推進(jìn)系統(tǒng)懸浮在等離子體中，電推進(jìn)殼體電位對應(yīng)于空間應(yīng)用時航天器表面的電位。因此，航天器表面相對于周圍等離子體的電勢，就可以通過測量電推進(jìn)外殼對真空艙殼體之間的電位（即懸浮電位）獲得，其實(shí)驗(yàn)結(jié)果如表1所列。

測試結(jié)果可看出，懸浮電位約為-10 V，該結(jié)果與DS-1及NSTAR-30在軌運(yùn)行測量數(shù)據(jù)-6～-10 V也相當(dāng)。即可認(rèn)為，SJ-9A衛(wèi)星應(yīng)用離子電推進(jìn)后，電推進(jìn)造成的高密度人工等離子體對航天器表面充電起到了保護(hù)作用，航天器懸浮電位被鉗制在一個較低的負(fù)電位。

3 總結(jié)及建議

美國、歐空局和俄羅斯等國得到的經(jīng)驗(yàn)表明：了解電推進(jìn)等離子體對航天器表面帶電的影響，并對其帶來的效應(yīng)進(jìn)行全面分析，是確保電推進(jìn)器在地球同步軌道衛(wèi)星上成功應(yīng)用的前提。國外除了開展電推進(jìn)羽流理論分析之外，還開發(fā)了相應(yīng)的軟件仿真工具，具體分析推進(jìn)器工作對衛(wèi)星表面電位的影響，為衛(wèi)星設(shè)計者進(jìn)行電推進(jìn)器的設(shè)計和布局以及控制航天器羽流效應(yīng)風(fēng)險提供支持?？紤]到目前研究現(xiàn)狀，電推進(jìn)等離子體對航天器帶電的影響需要開展以下幾個方面的研究工作：

（1）進(jìn)行比較充分的電推進(jìn)特殊環(huán)境特性及其對航天器表面帶電效應(yīng)的地面試驗(yàn)。環(huán)境特性的重點(diǎn)為束流等離子體和交換電荷等離子體；帶電效應(yīng)的重點(diǎn)為交換電荷等離子體在太陽電池陣上所引起的帶電效應(yīng)和推進(jìn)器羽流對航天器結(jié)構(gòu)電位的影響；

（2）國內(nèi)電推進(jìn)系統(tǒng)在軌飛行主要用于驗(yàn)證電推進(jìn)系統(tǒng)的技術(shù)可行性，而對電推進(jìn)引入的環(huán)境效應(yīng)關(guān)注較少。對比國外衛(wèi)星搭載的表面電位監(jiān)測儀等監(jiān)測設(shè)備，實(shí)時監(jiān)測電推進(jìn)器工作對衛(wèi)星表面的影響，如有可能應(yīng)在之后的飛行試驗(yàn)中增加電推進(jìn)系統(tǒng)空間環(huán)境效應(yīng)的在軌監(jiān)測手段和儀器設(shè)備。具體監(jiān)測儀器的設(shè)計制造工作可參考SMART-1進(jìn)行；

（3）配合地面模擬試驗(yàn)和空間飛行試驗(yàn)，建立電推進(jìn)羽流與航天器相互作用仿真分析軟件。根據(jù)美國NASA、歐空局的研發(fā)經(jīng)驗(yàn)，軟件研發(fā)過程中可以廣泛地通過國內(nèi)、國外合作，借鑒國外比較成熟的羽流等離子體模型，建立基于PIC/MCC方法的等離子體數(shù)據(jù)分析模型，開發(fā)羽流和航天器相互作用分析軟件。結(jié)合空間飛行試驗(yàn)數(shù)據(jù)，對分析模型進(jìn)行修正和驗(yàn)證，使其更加準(zhǔn)確和實(shí)用。

[1]張?zhí)炱?國外離子和霍爾電推進(jìn)技術(shù)最新進(jìn)展[J].真空與低溫，2006，12（4）:187-193.

[2]Dudzinski L A，Pen Pencil E J.Electric propulsion require?ments and mission analysis under NASA’s in-space propul?sion technology project[C]//30thInternational Electric Propul?sionConference，Italy，2007.

[3]Brosse S，Clerc S.Modeling of the plasma thruster impact on spacecraft charging[C]//The 29thInternational Electric Propul?sionConference，2005.

[4]劉磊，張慶祥，王立，等.電推進(jìn)羽流與航天器相互作用的研究現(xiàn)狀和建議[J].航天器環(huán)境工程，2011，28（5）:440-445.

[5]楊生勝，田愷，秦曉剛，等.電推進(jìn)羽流效應(yīng)分析[C]//電推進(jìn)系統(tǒng)環(huán)境影響研討會論文集，北京，2012:45-54.

[6]Tajmar M，Sedmik R，Scharlemann C.Numerical simulation of smart-1 hall-thruster plasma interactions[J].Journal of Pro?pulsionandPower，2009，25（6）:1178-1188.

[7]LikarJJ，BogoradAL，MalkoTR，etal.Interactionofcharged spacecraft with electric propulsion plume:On orbit data and ground test results[J].IEEE transactions on nuclear science，2006，53（6）:3602-3606.

[8]Lee J，Reissner A，Tajmar M，et al.Simulation of the Space?craftElectricPropulsionInteractiononDubaiSat-2usingSPIS[C]//The 32ndinteraction electric propulsion conference，Wies?baden，2011.

[9]KatzI，MandellM，JongewardG，etal.Theimportanceofaccu?rate secondary electron yields in modeling spacecraft charging[J].JournalofGeophysicalResearch:SpacePhysics，1986，91（A12）:13739-13744.

[10]Hastings D，Cho M，Chang P.Interactions between a plasma flowand a highly biased solar array[C]//In 25thJoint Propul?sionConference，1989.

[11]DeForest S E.Spacecraft charging at geosynchronous orbit[J].JGeophysRes，1972，77（4）:661-659

[12]Whipple E C.Potential of surface in space[J].Prog Phys，1981，44:1197

[13]Garrett H B，Whittlesey A C，Stevens N J.Design guidelines for assessing and controlling spacecraft charging effects[M].National Aeronautics and Space Administration，Scientific andTechnicalInformationBranch，1984.

[14]邵福球.等離子體離子模擬[M].北京：科學(xué)出版社，2002.

[15]Mandell M J，Davis V A，Mikellides I G.NASCAP-2K Pre?liminary Documentation[R].Science Applications Interna?tionalCorpsanDiegoCA，2002.

[16]Mandell M J，Stannard P R，Katz I.NASCAP Programmer’s ReferenceManual[M].NASA，1993.

[17]Mikellides I G，Mandell M J，Kuharski R A，et al.The Elec?tric PropulsionInteractions Code（EPIC）[C]//The 8thSpace?craftChargingTechnologyConference，AL，2004.

[18]InstitutionsFAQ.Thedevelopmentofaflexible，usableplas?mainteractionmodelingsystem[M].DefenseTechnicalInfor?mationCenter，2002.

[19]Roussel J F，Rogier F，Dufour G，et al.SPIS modelling capa?bilities，achievements and prospects[C]//Proceedings of 10thSpacecraftChargingTechnologyConference，2007.

[20]Arkhipov B A，Vertakov N M，Didenko B E，et al.Modelling and Experimental Verification of Hall and Ion Thrusters at ESTEC[C]//SpacecraftPropulsion，2000，465:683-690.