王雨瑋,任軍學(xué),*,吉林桔,湯海濱

1.北京航空航天大學(xué)宇航學(xué)院,北京100191 2.北京動(dòng)力機(jī)械研究所,北京100074

放電電壓和屏柵電壓對(duì)離子推力器性能的影響

王雨瑋1,任軍學(xué)1,*,吉林桔2,湯海濱1

1.北京航空航天大學(xué)宇航學(xué)院,北京100191 2.北京動(dòng)力機(jī)械研究所,北京100074

利用試驗(yàn)和數(shù)值模擬相結(jié)合的方法研究6 cm Kaufman離子推力器放電電壓和屏柵電壓的變化對(duì)其工作性能的影響。試驗(yàn)中,離子推力器使用氬氣作為推進(jìn)劑,測(cè)量了多組不同工況下的性能參數(shù)。此外,基于Goebel的理論模型模擬了放電電壓對(duì)束流電流和推進(jìn)劑利用率的影響;采用單元內(nèi)粒子-蒙特卡羅碰撞(PIC-MCC)方法模擬屏柵電壓對(duì)束流電流、推進(jìn)劑利用率和加速柵極電流的影響。試驗(yàn)和數(shù)值模擬結(jié)果一致,發(fā)現(xiàn)當(dāng)放電電壓逐漸增大時(shí),引出的束流電流和推進(jìn)劑利用率先增加然后趨于穩(wěn)定;當(dāng)屏柵電壓逐漸增大時(shí),引出的束流電流和推進(jìn)劑利用率先增加然后趨于穩(wěn)定,加速柵極電流先減小后趨于穩(wěn)定。研究可以為提高多模式離子推力器的性能提供參考。

離子推力器;多模式;放電電壓;屏柵電壓;氬氣

離子推力器具有比沖高、效率高、工作時(shí)間長(zhǎng)和安全環(huán)保等特點(diǎn),在空間任務(wù)中應(yīng)用廣泛。近年來(lái),隨著離子推力器的應(yīng)用逐漸由傳統(tǒng)的地球靜止軌道任務(wù)向深空探測(cè)和中低軌道高精度姿態(tài)控制等多元化任務(wù)發(fā)展,推力器在完成主推進(jìn)任務(wù)的同時(shí),有時(shí)還要兼顧航天器軌道保持和姿態(tài)控制的要求[1-3]。這就需要離子推力器可以在多模式下工作。目前,國(guó)外已經(jīng)成功研制了多種多模式離子推力器。例如應(yīng)用于DS-1航天器和DAWN航天器的30 cm NSTAR離子推力器,單臺(tái)推力器工作功率為0.49～2.31 k W,比沖為19 500～32 800 m/s,推力為19.5～92 m N,效率為38%～64%,共有5種推力模式[4-6]。XIPS-25推力器主要有2種工作模式,功率、比沖、推力、效率分別為4.5 k W、3.0 k W,35 000 m/s、34 000 m/s, 165 m N、79 m N,65%、63%[7]。NEXT推力器主要有2種工作模式,功率、比沖、推力、效率分別為6.9 k W、2.5 k W,41 900 m/s、40 000m/s,237 m N、80.4 m N,71%、64%[8]。

多模式離子推力器主要通過(guò)調(diào)節(jié)推力、比沖和效率來(lái)滿足不同任務(wù)的要求。加州理工大學(xué)的Goebel和Katz[9]及蘭州空間技術(shù)物理研究所的鄭茂繁和江豪成[10]認(rèn)為在不考慮束流發(fā)散和雙荷離子影響的情況下,推力和比沖的大小可以由束流大小反映;比沖的大小由屏柵電壓決定;利用束流電流、放電電壓和屏柵電壓可以計(jì)算放電損耗從而確定推力器的電效率。但是他們沒(méi)有從放電室的工作原理出發(fā)完整解釋放電電壓和屏柵電壓對(duì)束流電流的具體影響。北京航空航天大學(xué)的仇釬[11]、任軍學(xué)等[12-13]對(duì)屏柵的工作過(guò)程和束流的引出情況也進(jìn)行了仿真研究。但是他們的研究?jī)H限于對(duì)單個(gè)柵極孔的分析,沒(méi)有得到屏柵電壓與總的束流電流和總的加速柵極電流之間的關(guān)系曲線,也沒(méi)有進(jìn)行相關(guān)試驗(yàn)對(duì)結(jié)論進(jìn)行驗(yàn)證。

因此,本文首先從放電室微觀工作過(guò)程出發(fā),模擬分析6 cm Kaufman離子推力器放電電壓的變化對(duì)束流電流和推進(jìn)劑利用率的影響;屏柵電壓對(duì)束流電流、推進(jìn)劑利用率和加速柵極電流的影響。然后,進(jìn)行相關(guān)試驗(yàn)對(duì)分析結(jié)果進(jìn)行驗(yàn)證,得到參數(shù)之間的宏觀關(guān)系曲線。本文的研究可以為多模式離子推力器性能的提高和參數(shù)優(yōu)化設(shè)計(jì)提供參考,促進(jìn)其在工程上的應(yīng)用。

1 試驗(yàn)系統(tǒng)介紹

1.1 Kaufman離子推力器

典型的Kaufman離子推力器主要由主陰極、放電室、光學(xué)系統(tǒng)和中和陰極4部分組成。其工作原理是主陰極發(fā)射電子與放電室內(nèi)的中性氣體原子碰撞,使其電離產(chǎn)生離子;然后,屏柵極和加速柵極組成的光學(xué)系統(tǒng)將離子從放電室引出并利用柵極間的電壓將其加速噴出;最后中和陰極發(fā)射電子中和引出的離子束,使離開(kāi)推力器的粒子束保持電中性。典型的離子推力器結(jié)構(gòu)如圖1[14]所示。

本次研究使用6 cm Kaufman離子推力器,為便于試驗(yàn)操作,離子推力器的主陰極和中和陰極都使用直熱式鎢絲。主陰極燈絲位于放電室中央,采用雙陰極(一個(gè)陰極工作,另一個(gè)備用),可不開(kāi)真空罐切換陰極,能增加離子推力器的連續(xù)工作時(shí)間。放電室的結(jié)構(gòu)參數(shù)見(jiàn)表1。

光學(xué)引出系統(tǒng)采用雙柵極,屏柵極和加速柵極的孔徑和厚度都相等,都使用平面柵。柵極的結(jié)構(gòu)參數(shù)如表2所示。

圖1 Kaufman離子推力器的結(jié)構(gòu)[14]Fig.1 Structure of Kaufman ion thruster

表1 放電室的結(jié)構(gòu)參數(shù)Table 1 Structure parameters of discharge chamber

表2 柵極結(jié)構(gòu)參數(shù)Table 2 Structure parameters of screens

1.2 推進(jìn)劑儲(chǔ)供系統(tǒng)

推進(jìn)劑供給系統(tǒng)由氣瓶、壓力表、減壓閥、管路及質(zhì)量流量控制器等組成。試驗(yàn)采用氬氣作為推進(jìn)劑。使用S49 32/MT氣體質(zhì)量流量控制器,它使用氮?dú)鈽?biāo)定,本次試驗(yàn)使用的量程為5～500 m L/min。流量控制器的工作壓差范圍為50～300 k Pa,重復(fù)精度為±0.2%F.S。

1.3 真空系統(tǒng)

高空模擬真空系統(tǒng)主要包括模擬艙,真空泵和控制系統(tǒng)三個(gè)部分。真空艙的直徑為1.8 m。系統(tǒng)有三級(jí)泵,包括2臺(tái)高真空油擴(kuò)散泵,2臺(tái)羅茨泵和4臺(tái)機(jī)械泵。離子推力器點(diǎn)火前真空艙的本體壓力為4.7×10-3Pa,試驗(yàn)時(shí)模擬艙中的壓強(qiáng)不高于1.2×10-2Pa。

2 放電電壓對(duì)束流電流和推進(jìn)劑利用率的影響

2.1 理論計(jì)算方法

Goebel在2006年提出的會(huì)切磁場(chǎng)零維等離子體放電室模型得到了離子推力器放電損耗與推進(jìn)劑利用率之間的關(guān)系曲線。本節(jié)將基于Goebel的理論分析束流電流和推進(jìn)劑利用率隨放電電壓的變化情況。由于各個(gè)參數(shù)之間相互耦合,無(wú)法直接解出束流電流和推進(jìn)劑利用率隨放電電壓變化的關(guān)系式。因此,首先根據(jù)放電室內(nèi)中性氣體密度守恒和放電電流守恒求解出每個(gè)特定的放電電壓對(duì)應(yīng)的束流電流和推進(jìn)劑利用率的值。然后在2.2節(jié)中將其擬合,得到束流電流和推進(jìn)劑利用率隨放電電壓變化的關(guān)系曲線。

(1)中性氣體密度守恒公式

根據(jù)離開(kāi)放電室的中性氣體的流量Qout(未被電離的工質(zhì))等于進(jìn)入放電室的氣體流量Qin減去被電離并引出的離子束的氣體流量:

式中:Ib為引出的束流電流。離開(kāi)放電室的中性氣體的流量為

式中:n0為放電室內(nèi)中性氣體的密度;v0為中性原子速度;As為柵極面積;Ts為柵極透明度;ηc為克勞辛系數(shù)。推進(jìn)劑利用率ηm指的是產(chǎn)生單位束流電流需要消耗的推進(jìn)劑,是衡量放電室性能的關(guān)鍵參數(shù),其定義為

聯(lián)立式(1)～式(3)可以得到中性氣體密度的表達(dá)式為

此外,中性氣體的密度還可以由離子平衡方程得到,對(duì)于一個(gè)未磁化的等離子體來(lái)說(shuō),任意方向流出等離子體的電流Ii可以由波姆電流得到:

式中:ni為放電室內(nèi)離子的密度,約等于二次電子的密度ne;k為玻爾茲曼常數(shù);Te為電子溫度;mi為氬原子的質(zhì)量;Ai為所有離子吸收面積之和。電離室內(nèi)的離子由陰極發(fā)射的原初電子以及二次電子電離氣體產(chǎn)生,單位時(shí)間內(nèi)產(chǎn)生的離子數(shù)量Ip為

式中:np為原初電子密度;vp為原初電子速度;ve為二次電子的速度;σi為電離截面積;V為放電室的體積;〈σivp〉為原初電子的反應(yīng)速率;〈σive〉為二次電子的反應(yīng)速率?？偟碾x子產(chǎn)生速率等于總的離子吸收率:

聯(lián)立式(5)～式(7)可以得到:

聯(lián)立式(4)和式(8)可以得到中性氣體密度的守恒公式:

(2)放電電流守恒公式

首先,根據(jù)放電室內(nèi)輸入功率與輸出功率的守恒推導(dǎo)出與電子溫度和束流電流相關(guān)的放電損耗的表達(dá)式,具體過(guò)程如下。

求流動(dòng)到陽(yáng)極的二次電子電流Ia:式中:Id為放電電流;Iia為到達(dá)陽(yáng)極的離子電流;IL為原初電子電流。陰極發(fā)射的電流Ie為式中:Is為屏柵極電流;Ik為返回到陰極的離子電流。輸入放電室的功率Pin等于陰極發(fā)射電流乘以電子在電離室內(nèi)截獲的電壓Vk:

式中:Vd為放電電壓;Vc為陰極電壓;Vp為等離子體區(qū)電勢(shì);φ為相對(duì)于陽(yáng)極壁面的鞘層電勢(shì)。放電室消耗的功率Pout用于引出從陰極發(fā)射的原初電子以產(chǎn)生離子、激發(fā)態(tài)中性粒子和二次電子。離開(kāi)放電室到達(dá)電極的能量主要包括離子向陽(yáng)極、陰極和屏柵極的輸運(yùn),以及原初電子和二次電子與陽(yáng)極的碰撞。輸出功率等于以上各項(xiàng)之和,具體表達(dá)式為

式中:U+為氣體電離勢(shì)能;I*為激發(fā)態(tài)離子電流;U*為氣體激發(fā)勢(shì)能;TeV為電子溫度;ε為二次電子和陽(yáng)極壁面碰撞損失的能量。根據(jù)輸入功率等于輸出功率,聯(lián)立式(12)和式(13),同時(shí)利用式(10)和式(11)得到:

因此,放電損耗εe可以表示為然后,根據(jù)離子推力器的工作原理可得,引出的離子束流電流等于到達(dá)柵極的離子電流乘以柵極的透明度:

式中:va為離子聲速。將式(16)帶入式(15)(具體替換以及其他參數(shù)的帶入過(guò)程見(jiàn)文獻(xiàn)[9]),得到本節(jié)將使用的求解放電損耗公式:式中:Aas為陽(yáng)極壁面面積;fc為離子聲速和玻姆速度之比(Va/Vbohm);Ap為原初電子損失面積,根據(jù)式(17)得到放電電流的守恒公式:

放電電流守恒公式中的Ib由式(3)得到,ID是從試驗(yàn)中讀取的放電電流。

(3)特定放電電壓對(duì)應(yīng)的放電電流和推進(jìn)劑利用率求解方法

在放電室結(jié)構(gòu)參數(shù)和推力器工況確定的情況下,守恒式(9)和式(18)都只與電子溫度和束流電流相關(guān)。因此本節(jié)聯(lián)立式(9)和式(18)在每個(gè)放電電壓下計(jì)算都能得到一組電子溫度和束流電流的解滿足式(9)和式(18)。

2.2 結(jié)果與分析

從30 V開(kāi)始,以2.5 V為間隔,逐漸增加陽(yáng)極電壓,到60 V為止,求解每個(gè)特定的放電電壓對(duì)應(yīng)的束流電流和推進(jìn)劑利用率值。求解方法如2.1節(jié)所述。對(duì)數(shù)據(jù)進(jìn)行擬合,得到放電電壓對(duì)束流電流和推進(jìn)劑利用率之間的理論變化關(guān)系曲線。

同時(shí),利用試驗(yàn)得到每個(gè)放電電壓下,束流電流和推進(jìn)劑利用率的值,與理論結(jié)果進(jìn)行比較。試驗(yàn)工況如表3所示。

表3 改變放電電壓的試驗(yàn)工況Table 3 Conditions of discharge voltage changed experiments

圖2和圖3給出了束流電流和推進(jìn)劑利用率在相應(yīng)的陽(yáng)極電壓下理論計(jì)算和試驗(yàn)結(jié)果(放電電壓為42.5V時(shí),試驗(yàn)得到的束流電流和推進(jìn)劑利用率誤差較大,因此未予顯示)。

圖2 放電電壓對(duì)束流電流的影響Fig.2 Effect of discharge voltage on beam current

圖3 放電電壓對(duì)推進(jìn)劑利用率的影響Fig.3 Effect of discharge voltage on propellant utilization efficiency

由圖2和圖3可以看出,在30～60 V區(qū)間內(nèi),隨著放電電壓的增加,引出的束流電流和推進(jìn)劑利用率先增加然后趨于穩(wěn)定。這是因?yàn)橥屏ζ鞴ぷ鲿r(shí),電子是在磁場(chǎng)和電場(chǎng)的共同約束下轟擊推進(jìn)劑中性原子使其電離。只提高放電電壓時(shí),電場(chǎng)強(qiáng)度隨之增加,使得電磁場(chǎng)對(duì)電子的約束效果增強(qiáng)。與此同時(shí),電場(chǎng)強(qiáng)度增加也提高了電子運(yùn)動(dòng)速度,提高了電子的能量,增加了原子的電離率。更多的原子電離會(huì)使得放電室內(nèi)的離子密度增加,到達(dá)柵極的離子數(shù)會(huì)增加,所以其他條件不變的情況下,引出的離子束流電流和推進(jìn)劑利用率隨之增加。

當(dāng)放電電壓為45 V時(shí),理論和試驗(yàn)對(duì)應(yīng)的束流電流大小分別為110 mA和106 mA,推進(jìn)劑利用率大小分別為10.0%和9.6%。當(dāng)放電電壓大于45 V后,束流電流和推進(jìn)劑利用率的增速逐漸變緩,然后趨于穩(wěn)定。放電電壓過(guò)高,會(huì)產(chǎn)生大量雙核離子,導(dǎo)致離子的生成成本增加。此外,計(jì)算原初電子損失面積的表達(dá)式為[9]

式中:B為環(huán)尖處磁場(chǎng)強(qiáng)度;m為電子質(zhì)量;LC為磁尖總長(zhǎng)。從式(19)也可以發(fā)現(xiàn)當(dāng)放電電壓增加時(shí),磁尖處原初電子的損失面積增大,直接損失的原初電子增加,會(huì)影響電離率,使束流電流和推進(jìn)劑利用率的增速減慢。

3 屏柵電壓對(duì)束流電流、推進(jìn)劑利用率和加速柵極電流的影響

3.1 理論計(jì)算方法

采用單元內(nèi)粒子-蒙特卡羅碰撞(PICMCC)方法,模擬束流離子及CEX離子在柵極間的分布情況,可以得到在其他條件不變的情況下,每個(gè)特定屏柵電壓對(duì)應(yīng)的單個(gè)小孔的束流電流、加速柵極電流、離子的運(yùn)動(dòng)軌跡和分布情況。計(jì)算過(guò)程中,認(rèn)為電子是符合玻爾茲曼分布,采用的粒子比例是1 000,計(jì)算8 000步后,離子數(shù)穩(wěn)定。計(jì)算區(qū)域和等離子體參數(shù)的選擇見(jiàn)文獻(xiàn)[12]。

此外,總的束流電流和總的加速柵極電流由上述單個(gè)小孔的電流乘以小孔數(shù)量得到。推進(jìn)劑利用率利用式(3)計(jì)算得到。

3.2 結(jié)果與分析

從300 V開(kāi)始,以50 V為間隔,逐漸增加屏柵電壓,到600 V為止。首先根據(jù)3.1節(jié)的方法計(jì)算得到各個(gè)屏柵電壓對(duì)應(yīng)的束流電流和加速柵極電流的值。然后根據(jù)式(3)計(jì)算每個(gè)束流電流對(duì)應(yīng)的推進(jìn)劑利用率。將得到的束流電流、推進(jìn)劑利用率和加速柵極電流值進(jìn)行擬合,得到它們與屏柵電壓之間的關(guān)系曲線。同時(shí),利用試驗(yàn)得到每個(gè)屏柵電壓下,束流電流和推進(jìn)劑利用率的值,與理論結(jié)果進(jìn)行比較。試驗(yàn)工況如表4所示。

表4 改變柵極電壓的試驗(yàn)工況Table 4 Condition of screen grid voltage changed experiments

圖4和圖5給出了束流電流和推進(jìn)劑利用率在相應(yīng)的屏柵電壓下理論計(jì)算和試驗(yàn)結(jié)果。圖6為加速柵極電流隨屏柵電壓變化的試驗(yàn)結(jié)果和理論結(jié)果。

圖4 屏柵電壓對(duì)束流電流的影響Fig.4 Effect of screen grid voltage on beam current

圖5 屏柵電壓對(duì)推進(jìn)劑利用率的影響Fig.5 Effect of screen grid voltage on propellant utilization efficiency

從以上4幅圖可以看出,隨著屏柵電壓的增加,引出的束流電流和推進(jìn)劑利用率先增加然后趨于穩(wěn)定,加速柵極電流先減小后趨于穩(wěn)定。這種現(xiàn)象主要是導(dǎo)流系數(shù)變化引起的[15]。

在其他影響參數(shù)(柵極結(jié)構(gòu)、放電室等離子體密度、加速柵極電壓等)確定的情況下,屏柵電壓決定導(dǎo)流系數(shù)的大小。利用3.1節(jié)的方法,計(jì)算屏柵電壓為300 V、450 V和600 V三個(gè)工況下單個(gè)孔離子數(shù)量分布情況,分析在300～600 V范圍內(nèi)屏柵電壓對(duì)導(dǎo)流系數(shù)的影響,如圖7所示。

從圖7可以得到,屏柵電壓從300～600 V變化時(shí),隨著屏柵電壓的增加,導(dǎo)流系數(shù)不斷減小,放電等離子體鞘層弧度增加,離子束的入射角度發(fā)生變化,被聚焦程度越來(lái)越大。

因此,當(dāng)柵極電壓為300 V時(shí),導(dǎo)流系數(shù)過(guò)大,很多離子直接轟擊在加速柵極上游表面,導(dǎo)致引出的束流電流較小,加速柵極電流較大,產(chǎn)生濺射,影響柵極壽命。隨著屏柵電壓的增大,導(dǎo)流系數(shù)逐漸減小,束流電流和推進(jìn)劑利用率隨之增加,加速柵極電流隨之減小。從圖4～圖6綜合來(lái)看,離子推力器使用氬氣作推進(jìn)劑時(shí),在本節(jié)試驗(yàn)工況下,最佳導(dǎo)流系數(shù)大概出現(xiàn)在550 V附近。當(dāng)屏柵電壓超過(guò)550 V后,隨著柵極電壓的升高,導(dǎo)流系數(shù)變得過(guò)小,束流被過(guò)度聚焦,在柵極孔內(nèi)發(fā)生交叉,離子會(huì)轟擊加速柵極孔內(nèi)壁和表面,增加損耗。因此束流電流的增速逐漸減小,電流值逐漸趨于穩(wěn)定;加速柵極電流減少的速率也變慢,電流值逐漸趨于穩(wěn)定。

此外,圖6中加速柵極電流的試驗(yàn)數(shù)據(jù)要比理論數(shù)據(jù)值高。將圖6與NSTAR-30離子推力器在功率2.3 k W、推力92.7 m N,比沖31 270 m/s、效率61.8%工況下模擬得到的加速柵極電流隨屏柵電壓的理論變化情況(如圖8[9]所示)對(duì)比,可以得到圖7的理論數(shù)據(jù)可靠性較高。

圖6 加速柵極電流隨屏柵電壓變化的試驗(yàn)結(jié)果和理論結(jié)果Fig.6 Experiment and theory result of acceleration grid electric current’s changing with screen grid voltage

圖7 屏柵電壓為300 V、450 V和600 V時(shí)離子數(shù)量分布Fig.7 Distribution of ion numbers at the screen grid voltage of 300 V,450 V and 600 V

圖8 NSTAR-30推力器加速柵極電流隨屏柵電壓的變化Fig.8 Changing of accelerating grid current with screen grid voltage’s changing of NSTAR-30 thruster

而造成試驗(yàn)中加速柵極電流過(guò)高的原因主要有兩個(gè):一是本次試驗(yàn)使用的推力器是用于高校教學(xué)研究的產(chǎn)品,推力器的加工工藝有局限性,尤其是屏柵極和加速柵極的平行度、小孔的對(duì)準(zhǔn)精度不夠高,會(huì)導(dǎo)致加速柵極電流增加;二是試驗(yàn)過(guò)程中,屏柵極和加速柵極的熱變形以及點(diǎn)火后真空艙內(nèi)真空度的提高也會(huì)造成的加速柵極電流的增加。

4 結(jié)束語(yǔ)

本文通過(guò)試驗(yàn)和理論模型相結(jié)合的方法,分析了6 cm Kaufman離子推力器在使用氬氣作為推進(jìn)劑時(shí),放電電壓對(duì)束流電流、推進(jìn)劑利用率的影響;以及屏柵電壓對(duì)束流電流、推進(jìn)劑利用率、加速柵極電流和導(dǎo)流系數(shù)的影響,并給出了參數(shù)之間宏觀的變化關(guān)系曲線。理論模擬和試驗(yàn)結(jié)果相符合,結(jié)論如下:

1)只改變放電電壓時(shí),隨著電壓的增加,引出的束流電流和推進(jìn)劑利用率先增加然后趨于穩(wěn)定。

2)只改變屏柵電壓時(shí),隨著屏柵電壓的增加,導(dǎo)流系數(shù)不斷減小,離子束被聚焦的程度越來(lái)越高;引出的束流電流和推進(jìn)劑利用率先增加,到達(dá)最佳導(dǎo)流系數(shù)狀態(tài)后趨于穩(wěn)定;加速柵極電流逐漸下降,到達(dá)最佳導(dǎo)流系數(shù)狀態(tài)后趨于穩(wěn)定。對(duì)于多模式的離子推力器,設(shè)定工作模式時(shí),應(yīng)盡量避免使用過(guò)高的放電電壓,減少雙核離子的產(chǎn)生電功率的損失;而且應(yīng)盡量使推力器工作在最佳導(dǎo)流系數(shù)狀態(tài)附近。

References)

[1] 鄭茂繁,耿海,梁凱,等.用于小行星探測(cè)的離子推力器技術(shù)研究[J].深空探測(cè)學(xué)報(bào),2015,2(3):236-240. ZHENG M F,GENG H,LIANG K,et al.Research on ion thruster technology for asteroid exploration[J]. Journal of Deep Space Exploration,2015,2(3):236-240(in Chinese).

[2] 溫正,王敏,仲小清.多任務(wù)模式電推進(jìn)技術(shù)[J].航天器工程,2014,23(1):118-123. WEN Z,WANG M,ZHONG X Q.Multitask mode electric propulsion technologies[J].Space Engineering, 2014,23(1):118-123(in Chinese).

[3] OH D Y,SANTIAGO G.Analytic optimization of mixed chemical-electric orbit raising missions[C]∥Proceedings of the 26th International Electric Propulsion Conference,2001.

[4] RAWLIN V K,SOVEY J S,HAMLEY J A.An ion propulsion system for NASA’s deep space missions, AIAA-1999-4612[R].Reston:AIAA,1999.

[5] CARDELL G,UILOA-SEVERINO A,GROSS M.The design and operation of the DAWN power system, AIAA-2012-3898[R].Reston:AIAA,2012.

[6] BROPHY J R,MARCUCCI M G,GANAPATHI C B,et al.The ion propulsion system for DAWN, AIAA-2003-4542[R].Reston:AIAA,2003.

[7] GOEBEL D M,MARTINEZ-LAVIN M,BOND T A,et al.Performance of XIPS electric propulsion in on-orbit station keeping of the Boeing 702 spacecraft, AIAA-2002-4348[R].Reston:AIAA,2002.

[8] HERMAN D A,SOULAS G C,PATTERSON M J. Performance evaluation of the prototype model NEXT ion thruster,AIAA-2007-5212[R].Reston:AIAA, 2007.

[9] GOEBEL D M,KATZ I.Fundamentals of electric propulsion:ion and hall thrusters[M].New Jersey: John Wiley&Sons Inc,2008:91-232.

[10] 鄭茂繁,江豪成.離子推力性能評(píng)價(jià)方法[J].真空與低溫,2012,18(4):223-227. ZHENG M F,JIANG H C.Method of performance evaluation for ion thruster[J].Vacuum&Cryogenics, 2012,18(4):223-227(in Chinese).

[11] 仇釬.離子發(fā)動(dòng)機(jī)交換電荷離子分布的數(shù)值模擬[D].北京:北京航空航天大學(xué),2010. QIU Q.Numerical simulation of charge-exchange propagation in the plume of ion thruster[D].Beijing: Beihang University,2010(in Chinese).

[12] 任軍學(xué),謝侃,湯海濱,等.離子發(fā)動(dòng)機(jī)加速柵極孔擴(kuò)大腐蝕的粒子模擬[J].推進(jìn)技術(shù),2013,34(10): 1432-1440. REN J X,XIE K,TANG H B,et al.Particle simulation of the acceleration grid aperture erosion in ion thruster[J].Journal of Propulsion Technology, 2013,34(10):1432-1440(in Chinese).

[13] 任軍學(xué),顧左,郭寧,等.離子發(fā)動(dòng)機(jī)羽流特征的數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào),2013,28(6):1372-1378. REN J X,GU Z,GUO N,et al.Numerical simulation of characteristics of ion thruster plume[J]. Journal of Aerospace Power,2013,28(6):1372-1378(in Chinese).

[14] WIRZ R E.Discharge plasma process of ring-cusp ion thrusters[D].Pasadena:California Institute of Technology,2005.

[15] 于達(dá)仁,劉輝,丁永杰,等.空間電推進(jìn)原理[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2012:166-167. YU D R,LIU H,DING Y J,et al.Fundamentals of space electric propulsion[M].Harbin:Harbin Institute of Technology Press,2012:166-167(in Chinese).

(編輯:范真真)

Effects of discharge voltage and screen grid voltage on performance of ion thruster

WANG Yuwei1,REN Junxue1,*,JI Linju2,TANG Haibin1
1
.School of Astronautics,Beihang University,Beijing 100191,China 2.Beijing Power Machinery Institute,Beijing 100074,China

The effects of discharge voltage and screen grid voltage changing to working performance of a 6 cm Kaufman ion thruster were studied by combination of experiments and numerical simulation.In these experiments,the argon was used as propellant for the ion thruster,and many groups of data were gathered.Furthermore,based on Goebel's theory model,the effects of discharge voltage changing on beam current and propellant efficiency were simulated.And the particle in cell Monte-Carlo collide(PIC-MCC)theory was used tosimulate the effects of screen voltage on beam current,propellant efficiency and acceleration grid current.Experiments and numerical simulation drew the same conclusions.When discharge voltage increases,ion beam current and propellant utilization efficiency tend to be stable after a corresponding increase.A similar beam electric current and propellant utilization efficiency changing pattern is observed when the screen grid voltage is increased,but acceleration grid voltage tends to be stable after a corresponding decrease.This study can provide a reference for increasing multi-modes ion thrusters'performance.

ion thruster;multi-modes;discharge voltage;screen grid voltage;argon

V439+.1

:A

10.3780/j.issn.1000-758X.2016.0014

2015-11-12;

:2015-12-10;錄用日期:2015-12-30;< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間

時(shí)間:2016-02-24 13:31:27

http:∥www.cnki.net/kcms/detail/11.1859.V.20160224.1331.005.html

北京市高等教育“青年英才”計(jì)劃(YETP1129)

王雨瑋(1992-),男,碩士研究生,412544372@qq.com

*通訊作者:任軍學(xué)(1980-),男,副教授,rjx_buaa@163.com,主要研究方向?yàn)榛鸺l(fā)動(dòng)機(jī)及空間電推進(jìn)

王雨瑋,任軍學(xué),吉林桔,等.放電電壓和屏柵電壓對(duì)離子推力器性能的影響[J].中國(guó)空間科學(xué)技術(shù),2016,36(1):77-84.WANG Y W,REN J X,JI L J,et al.Effects of discharge voltage and screen grid voltage on performance of ion thruster[J].Chinese Space Science and Technology,2016,36(1):77-84(in Chinese).

http:∥zgkj.cast.cn