離子霍爾推力器束流中和與耦合研究進展

王同宇，張天平，陳娟娟，楊 樂

（蘭州空間技術物理研究所 真空技術與物理重點實驗室，蘭州 730000）

0 引言

離子推力器和霍爾推力器是目前航天器上使用最廣泛的動力裝置之一，相比而言，在同樣功率下，離子推力器比沖高、推力小，霍爾推力器比沖低、推力大。為有效發揮離子推力器高比沖、霍爾推力器大推力的優勢，解決傳統離子和霍爾電推力器高比沖和大推力無法兼具的難題，提出了離子霍爾混合推力器概念。

離子霍爾混合推力器由環型離子推力器[1]和霍爾推力器[2]同軸嵌套組成。在兩種推力器同時工作情況下，由于束流等離子體電勢、密度、速度等方面的顯著差別，同時共用中和器，使得混合推力器束流特性與單一電推力器工作情況存在很大不同。離子霍爾混合電推力器的近羽流場等離子耦合是一個復雜的物理過程，分析離子推力器束流離子與霍爾推力器束流離子的耦合過程以及束流離子和電子中和反應過程，可為離子霍爾混合推力器參數優化提供技術支持。

在調研了國內外離子霍爾混合推力器研究的基礎上，重點關注了推力器束流耦合研究進展，并初步分析了離子霍爾混合推力器束流耦合的特殊問題。

1 離子霍爾混合推力器

美國NASA的Dankanich等[3]在2006年對航天器搭載多個離子推力器或者霍爾推力器，以及混合搭載離子推力器和霍爾推力器的方案進行了分析，研究發現合理的搭配使用離子推力器和霍爾推力器可以實現在保證可靠壽命的前提下提高功推比和功率。通過混合使用推力器一方面需要額外設計鑒定試驗及技術分析；另一方面飛行器搭載多個推力器會占用太多的空間，增加的推進劑補給線路及其他額外線路設計使航天器設計更加笨重，使得熱條件和操作環境的約束更多。所以相較單一的電推力器，使用多種不同的電推力器組合，有利于簡化電推進系統設計。

2011年NASA格倫研究中心的Patterson等[4]研究提出了一種離子霍爾混合推力器的設計概念，可以滿足未來電推進技術高比沖、大推力的技術要求。圖1為美國格林中心離子霍爾混合推力器結構示意圖。

圖1 離子霍爾混合推力器結構示意圖Fig.1 Theoretical model of ion and Hall hybrid-thruster

離子霍爾混合推力器由新型環型離子推力器和傳統霍爾推力器同軸嵌套組成，外環為環型離子推力器，內環為傳統型霍爾推力器。采用環型離子推力器有利于實現混合推力器一體化結構實現，由于環型離子推力器的放電室和陽極呈環狀結構，可以增大放電室內部放電陽極面積，提高放電電流和等離子體分布均勻性，保持更小間距下的較高柵極離子透明度，解決大口徑柵極加工制造、熱變形及柵極熱穩定性等技術難題[5-6]。

離子霍爾混合推力器中和器設計位于推力器外圍中心軸線上，既作為空心陰極為霍爾推力器等離子體產生提供電子源，又作為中和器為混合推力器提供用于束流中和的電子，該結構可減小推力器外部輪廓，大幅減小推力器幾何結構尺寸，便于加工制造組裝，具備應用范圍更廣的優勢。

國內開展離子霍爾混合推力器研究的機構是蘭州物理空間技術研究所[7-8]，該所在完成環型離子推力器研制的基礎上，設計了離子外環霍爾內環的混合推力器結構，并設計了環型離子推力器和霍爾推力器組合時的磁路組件，以滿足推力器磁場設計準則，同時設計了環型套筒磁屏蔽結構來隔離霍爾內環與離子外環間的磁場相互干擾。目前離子霍爾混合推力器原理樣機的設計為：整體外圓直徑約57 cm、功率1.5～6.5 kW、推力80～210 mN、比沖1500～5000 s。

2 束流中和分析

通常認為離子束中和過程可以分為三個步驟：電子-離子摻混、電子-離子耦合及電中性等離子體束推進。電子-離子混摻過程中的主要影響因素是陰極的安裝位置，推力器的電中性等離子體主要是通過電子-離子的耦合產生的，因此對束流中和的研究也主要集中在電子-離子耦合方面。

2.1 離子推力器束流中和

離子推力器主要由空心陰極、放電室、離子光學系統和中和器四個部分組成。通過這些組件的協調工作，離子推力器可以順利完成預期的使命。中和器的作用是發射電子中和推力器產生的離子束流使其達到宏觀電中性。整個中和過程機制復雜，主要涉及粒子間相互作用、粒子與場之間相互耦合，這些物理機制將導致非磁化等離子體發生振蕩，而碰撞過程會導致等離子體振蕩的幅度隨時間衰減。振蕩過程使等離子體逐漸從非平衡態向近平衡態演化。當離子束不能發生有效中和時，很大一部分離子會返回到航天器表面。這一問題首先是由Staff[9]在1960年關于靜電推力器的論文中提出的。

圖2 離子推力器羽流示意圖Fig.2 Schematic diagram of ion thruster plume

在人們認識到離子推力器和霍爾推力器束流需要電子中和后，許多研究人員對離子束流中和現象進行了分析。French[10]發現在中和過程中電子電流密度存在著振蕩現象。Mirels[11]在1961年發現中和器位置和熱電子的運動可能對離子束中和影響不大。Marshall等[12]指出如果電子從束流區外部固定位置通過漂移運動進入束流區，在離子束的中心仍會產生虛陰極并使推力降低。

為了解決離子束中和過程的相關問題，研究者采用數值方法進行研究。1965年，Brauch等[13]采用二維數值模擬方法讓電子通過在束流區的振蕩來實現中和，同時理論描述了振蕩并非中和過程的唯一機制。之后還有研究表明空間電荷場的波動使得熵增加從而使電子混合，但是在中和器下游并未發現有波動現象存在。

1979年Parks等[14]利用理論分析的方法研究了羽流區等離子體的演化過程。假設束流離子和中和電子之間存在著一個靜電勢阱，呈麥克斯韋分布的電子在其約束下與束流離子不斷發生碰撞，另外，通過試驗測量束流中和過程中電子的溫度變化，發現在整個束流區域內電子的溫度始終保持不變，且其與空間的電位之間滿足barometric定律。因此，在理論分析時假設空心陰極發射的電子呈熱平衡狀態，該分析得到了束流的密度和電位，以及目標粒子與推力器出口位置之間的距離滿足高斯定律。該模型最大的缺點在于假設束流離子被噴出后立刻達到了穩態，整個分析過程忽略了動態束流的存在。

近幾年，許多關于離子推力器離子束中和過程的數值模擬研究陸續出現。

Othmer等[15-17]采用三維PIC模型進行中和過程的模擬，重新對中和機理進行了更深入的研究，揭示了“噴注速度比”這一重要的物理量。

噴注速度比為：

Brieda等[18]使用了名為DRACO的三維等離子體模擬軟件，采用全粒子PIC方法對單離子推力器近羽流區中和過程進行了模擬。通過“多域模型”獲得了中和電子的速度分布，所謂“多域模型”是將近陰極區與近羽流區進行區分，這種模型使計算具有很高的真實度。將該模型與其他預混合束流模型和浮動陰極電位模型進行了比較，近羽流數值模擬研究結果表明陰極電子的引入導致了束流核心的振動。要獲得類似于玻爾茲曼分布的密度分布須要實現振動的衰減，但是在模型中沒有得到相應的結果。

Wang等[19]研究發現束流中和與傳播過程是兩個耦合的過程，熱電子在冷離子中最初的擴散過程會在束流中形成一個正電勢阱。這個伴隨束流離子傳播的勢阱會束縛隨后進入的熱電子，隨著沿束流方向電勢阱與被束縛電子之間的作用電子離子的耦合得以實現。

2.2 霍爾推力器束流中和

不同于離子推力器，霍爾推力器沒有柵極結構，霍爾推力器由陶瓷加速通道和空心陰極兩部分構成，空心陰極不僅給推力器提供放電室電離電子，又作為中和器中和離子束流。在空心陰極和加速通道之間存在一個區域，叫做耦合區，研究認為束流中和在這個區域已經發生。Alebarede等[20]的研究甚至懷疑霍爾推力器使工質發生電離的電子和中和電子是否存在交界，試驗發現在陰極出口電子會呈現線性分布如圖3所示。

圖3 耦合放電近陰極區光強分布圖Fig.3 The light distribution of coupling discharge near the cathode

Sommerville等[21]通過實驗發現，陰極的位置對推力器近場束流中和有顯著影響，當陰極距離加速通道越遠時，電子則須要克服磁場區域作用才能到達離子束流。實驗證明可以通過降低陰極電勢同時提高近場束流電勢保證電子能到達束流離子區域。此外，發現霍爾推力器束流中和電子傳導過程與背壓有關，隨著背壓增大，碰撞頻率增大，傳導阻抗減小，有利于減小束流耦合電壓[20]。

Hofer等[22]的試驗表明，將陰極放置在推力器的中央可以減小束流發散角并且可以產生軸對稱的羽流，另一方面中央放置的陰極會消除外加能量產生的偏航扭矩從而可以帶來推力的微小提升。

Smith等[23]對近場區電子行為進行了蒙特卡洛模擬，研究表明，電子能量分布從通道口往束流下游遷移過程中，電子溫度會逐漸升高，其能量分布也從麥克斯韋分布逐漸變為非麥克斯韋分布。Knoll等[24]實驗測量到近場區存在多種振蕩，這些振蕩會影響電導特性，影響束流及電子的特性，從而影響中和過程。

2.3 混合推力器束流中和基本問題

蘭州空間技術物理研究所研制的離子霍爾混合推力器具有離子外環、霍爾內環的結構，共用一個中和陰極。根據其結構特點，有霍爾工作模式、離子工作模式和混合工作模式。存在的問題有：其中一種推力器工作時啟動另一推力器的過程對混合推力器束流中和的影響以及對中和器產生的影響；一種推力器工作時關閉另一推力器的過程對混合推力器束流中和的影響以及對中和器的影響；一種推力器工作時改變另外一個推力器的工況對推力器束流中和的影響以及對中和器的影響。由此可見，中和器的適應能力決定推力器能否可靠工作。

陰極的穩定性及位置對推力器工作至關重要，混合推力器的多種工作模式啟動過程，可能會導致陰極不穩定，造成推力器束流中和發生震蕩進而造成推力器放電不穩定，甚至造成推力器突然熄火的危險。推力器工作的瞬態變化甚至會使中和電子與沿程束流離子碰撞引起的振蕩最終耦合進電源中，從而影響推力器性能，因此這對陰極技術提出了很高的要求。

3 推力器組束流耦合試驗

3.1 離子推力器組束流耦合

美國NASA格林研究中心[25-26]對離子推力器組工作情況進行了一系列研究，將4個7 kW的離子推力器組成如圖4擺放的組合，比較了不同推力器組合工作在高中低三個功率等級情形下的束流特性與單獨一個推力器工作時的束流特性，研究發現：相鄰兩個或者三個推力器以相同功率同時工作時，其比沖、推力以及效率與單一推力器工作時相比基本不變，但是推力器柵極沖擊電流增大。通過分析實驗數據認為當多個離子推力器同時工作時束流環境中中性氣體密度增大，導致束流中CEX碰撞頻率增大，而不是相鄰推力器CEX離子電流的遷移造成的。

調整一個離子推力器的出口平面方向，使得兩束束流形成大約10°的夾角，發現對束流特性影響很小。表明當兩束呈小角度的等離子體之間的耦合對整體束流特性影響較小。

盡管實驗驗證了離子推力器組工作的可靠性，但是實驗中每個推力器工作在相同的功率下，并沒有對離子推力器組在不同的工作功率情況下的束流特性進行研究。

圖4 “3+1”NEXT離子推力器組圖[25]Fig.4 “3+1” configuration NEXT multi-thruster array

3.2 霍爾推力器組束流耦合

霍爾推力器的組合可以顯著的提高推力系統的功率水平[27]，但是方案是否可靠需要對推力器束流之間的相互作用進行分析研究。Lobbia等[28]對霍爾推力器組的性能進行了實驗，如圖5所示，實驗測試了推力器工作在200 V、300 V、400 V情況下的表現，研究發現霍爾推力器組中每個推力器產生的推力都大于推力器單獨工作時的推力水平。并且隨著推力器數量的增多，也會對推力器的推力帶來較小的增幅，實驗也表明推力器組同時工作也會使推力器的效率和比沖產生微小的增幅。

圖5 四霍爾推力器組實驗圖[28]Fig.5 Acluster of four Hall thrusters[28]

美國密歇根大學Beal等[29]使用阻滯勢分析儀（RPA）對霍爾推力器組束流等離子體耦合情況進行了研究，試驗表明，相比單一霍爾推力器，霍爾推力器組束流中低能離子數量顯著增大，尤其在推力器中心線低角度附近。研究者認為這一現象主要是兩束束流中初始離子的彈性散射造成的。

Saverdi等[30]在對HT-100霍爾推力器組進行試驗時發現，當兩臺并列放置的霍爾推力器同時工作時，其束流發散角相比單一推力器時增大9%，同時在偏離推力器中心線的低角度區域，低能離子質量濃度有明顯的增加趨勢。

3.3 混合推力器束流耦合關鍵問題

離子霍爾混合推力器獨特的設計，使得束流特性更加復雜，主要有以下幾個問題。

3.3.1 束流電勢差異

前文提到離子霍爾混合推力器束流不同于傳統單一電推力器產品，在推力器混合工作模式下，將有兩束能量、密度、電勢不同的束流從兩個同軸通道口噴出，電推力器束流服從高斯分布：

其中ρ0是初始離子質量濃度。

式中：M為離子質量；Ib為束流電流；vb為束電壓，R為束流初始寬度；h(t)為描述束流發散的量。由于5 kW級環型離子推力器屏柵極電壓為1 500 V，加速柵電壓為-250 V，霍爾推力器陽極電壓為300 V，由此可以知道離子霍爾混合推力器中離子推力器出口離子質量濃度低于霍爾推力器出口離子質量濃度，如圖6所示。這種情況導致最直接的問題是在推力器出口位置徑向形成電勢梯度，為了防止離子霍爾兩個部分產生電勢差而對推力器產生損傷，應該首先啟動中和陰極，保證推力器發射的離子能夠快速中和形成準中性等離子體。

圖6 5 kW級環型離子推力器及1.5 kW級霍爾推力器束流質量濃度分布圖Fig.6 The distribution of 5 kW ion thruster and 1.5 kW Hall thruster's beam mass density

3.3.2 束流離子速度差異

另一個特殊的問題是數值上環型離子通道口噴出的離子速度是霍爾通道噴出的2倍多，因此束流耦合模型中需要處理快離子和慢離子的耦合過程。一價慢離子和快離子之間碰撞如果產生電離，屬于基態重離子之間的碰撞電離，這類電離一般比電子碰撞電離更難發生，需要的膨脹能量在數百電子伏以上。假設考慮離子間的彈性碰撞，能量損失系數為：

由式（3）可知，質量為M，動能為ε的離子即使與另一離子最外層質量為m的電子碰撞，也只能把很小的一部分動能傳遞給電子。所以快離子和慢離子之間碰撞產生的電離可以忽略，只需要分析之間的彈性碰撞造成的能量變化對束流特性的影響，結合研究分析，當兩束速度不同的束流相互作用時，低能離子的密度會在束流邊界區域呈現增加的趨勢。

離子霍爾混合推力器有多種工作模式，離子工作模式、霍爾工作模式、離子霍爾工作模式。在推力器工作模式轉換的過程中，束流耦合的瞬態變化以及如何自適應反應是相當復雜的問題。

結合離子推力器組和霍爾推力器組實驗結果，認為當兩束束流作用時，首先電荷交換離子密度會明顯上升，并會在橫向快速下降。當霍爾部分工作過程中，啟動環型離子部分，推力器下游的背景氣體壓力的瞬變，會使霍爾束流的離子能量發生瞬變，這個變化可能引起電子的振蕩會造成霍爾部分的放電不穩定；反之，當環型離子工作時，突然啟動霍爾部分，背景氣體壓力的變化對離子束流的影響較小。這對推力器共用陰極的設計提出了更高的要求。

4 結論

在調研國內外離子霍爾混合推力器研究現狀的基礎上，著重關注由新型推力器帶來的離子和霍爾推力器束流耦合問題，介紹了離子推力器組和霍爾推力器束流組耦合研究進展；分析了束流中和研究方法的演變及成果，初步梳理了離子霍爾混合推力器束流耦合的關鍵技術難點及可能的解決辦法，為推動我國離子霍爾混合推力器的研發提供參考。

等離子體推力器束流耦合問題是電推力器技術經典問題之一，但是目前仿真研究手段及研究成果都不成熟。對離子束流耦合的深入研究不僅有利于離子霍爾混合推力器的技術發展，也將促進電推力器羽流理論的不斷成熟。