基于CFD的10 cm氙離子推力器陽極推進劑供給方式優(yōu)化

胡竟，楊福全，郭德洲，高俊，鄭茂繁

（1.蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點實驗室，蘭州730000；2.北京航空航天大學 材料科學與工程學院，北京100083）

離子推力器以其壽命長［1-2］、比沖高［3］、工作模式多［4］、推力精確可調(diào)［5］等優(yōu)點在航天器軌控任務中得到廣泛應用。以滿足重力梯度測量衛(wèi)星無拖曳飛行任務及超低軌高分辨率地球觀測衛(wèi)星軌道維持任務為目標，中國開展了500 W 級10 cm氙離子推力器研制［6］。

處于離子推力器放電室結(jié)構(gòu)上游的陽極環(huán)，其陽極推進劑氙氣的供給方式是影響放電室內(nèi)陽極推進劑濃度分布及其變化梯度的核心，并最終影響到放電室內(nèi)氙氣的電離效率及放電穩(wěn)定性，因此，以陽極推進劑供給方式設計為核心的放電室內(nèi)陽極推進劑濃度分布及其變化梯度的設計是放電室放電模式可靠性設計的關鍵技術(shù)之一。而10 cm氙離子推力器之所以能夠在1～20 mN范圍內(nèi)以12μN的推力分辨率實現(xiàn)連續(xù)可調(diào)，就是因為推力器陽極環(huán)拓撲結(jié)構(gòu)設計綜合考慮了放電室磁場變化、放電電流變化及陰極工作點變化對離子推力器工作效率、工作穩(wěn)定性及氙氣電離效率的影響。

美國NASA［7-8］研究了陽極推進劑供給方式對放電室內(nèi)推進劑電離效率的影響，通過優(yōu)化NSTAR-30 cm離子推力器放電室的陽極環(huán)結(jié)構(gòu)，有效改善了引出束流的平直度。Nakayama和Narisawa［9］開展了推進劑自陽極中心單出口供給和多出口環(huán)管供給對其在放電室內(nèi)分布影響的研究，發(fā)現(xiàn)推進劑的多出口環(huán)管供給可有效解決放電室內(nèi)推進劑分布不均勻且平均壓強低的問題，極大提高了推進劑在放電室的駐留時間，提升了整體電離效率，改善了推進劑利用率。Nakayama等［10］通過優(yōu)化陽極環(huán)結(jié)構(gòu)提高了14 cm離子推力器的工作性能。Reid［11］采取數(shù)值模擬與試驗驗證相結(jié)合的方式開展了陽極環(huán)結(jié)構(gòu)對推力器工作性能的影響分析研究。Janes和Lowder［12］采取等離子體診斷的手段開展了放電室內(nèi)推進劑分布均勻性的影響因素研究。

現(xiàn)有陽極推進劑供給方式的影響主要從宏觀角度對陽極環(huán)結(jié)構(gòu)、安裝位置與推力器工作性能影響作用關系開展研究，而對于供給方式中較為微觀的推進劑輸入方式、推進劑分配方式研究較少，有待進一步深入研究。目前，中國的離子推力器放電室陽極環(huán)結(jié)構(gòu)設計主要依靠工程經(jīng)驗，借助試驗測試，缺少理論指導。有鑒于此，以陽極環(huán)拓撲結(jié)構(gòu)優(yōu)化為基礎，研究分析陽極推進劑供給方式對離子推力器放電室內(nèi)陽極推進劑分布特性及電離效率的影響具有重要的意義。針對上述問題，以滿足超低軌高分辨率地球觀測衛(wèi)星在軌多目標、差異化軌控任務對10 cm氙離子推力器的應用需求為目標，通過開展以提高放電室內(nèi)陽極推進劑供給均勻性為目標的陽極推進劑供給方式優(yōu)化研究，實現(xiàn)陽極推進劑利用率的提升，進而實現(xiàn)整機效能的提升。為此，本文從離子推力器工程應用角度出發(fā)，運用計算流體動力學（Computational Fluid Dynamics，CFD）理論［13］，研究了陽極推進劑供給方式對陽極環(huán)和放電室空腔內(nèi)陽極推進劑分布特性的影響作用關系，并將優(yōu)化分析結(jié)果在10 cm氙離子推力器中進行對比驗證。

1 離子推力器原理與陽極推進劑供給

1.1 離子推力器工作原理

如圖1所示，給中和器、主陰極和放電室陽極供給不同流量的推進劑氙氣后，主陰極產(chǎn)生的電子受電場力作用進入放電室。在放電室里，主陰極發(fā)射的原初電子、氙氣被碰撞并電離后產(chǎn)生的二次電子在電場和磁場的共同作用下，繞磁力線以螺旋振蕩運動的形式向放電室陽極運動。磁場的存在，通過延長電子的運動路徑，進而提高了電子與氙氣的碰撞電離幾率，最終增加了放電室內(nèi)氙離子的濃度。至此，放電室內(nèi)就形成了由氙氣、電子和氙離子組成的放電等離子體云。等離子體云中的電子基本被陽極吸收，而離子則受推力器結(jié)構(gòu)下游的柵極組件電場作用，被聚焦、加速、引出并產(chǎn)生推力。被引出的離子束流與中和器發(fā)射的電子會進行中和，形成近中性的粒子束流，以保證航天器呈電中性。至此，以推進劑氙氣作為介質(zhì)完成了電能到動能的轉(zhuǎn)換并產(chǎn)生推力的過程。

由此發(fā)現(xiàn)，延長電子運動路徑、提高放電室推進劑氙氣駐留時間可有效提升電子與氙氣的碰撞電離幾率，并最終改善放電室陽極推進劑利用率。

1.2 放電室陽極推進劑分布特性影響分析

對離子推力器而言，放電室內(nèi)陽極推進劑分布特性是決定陽極推進劑電離位置及等離子體云宏觀形態(tài)的關鍵因素。離子推力器的放電室陽極推進劑利用率（η）及電離損耗（ζ）與推力器結(jié)構(gòu)參數(shù)、工作參數(shù)密切相關，可分別表示為

式中：Vd、Vc和Vp分別為放電電壓、主陰極觸持電壓和等離子體電壓；Ib、Ip和I*分別為屏柵電流、放電室中總的離子電流和放電室中激發(fā)態(tài)離子電流；Ts、TeV分別為屏柵有效透明度、電子溫度；Ap、As和Aas分別為原初電子碰撞陽極的面積、屏柵面積和陽極面積；U＋、U*分別為推進劑氙氣的電離電勢、激發(fā)電勢；np、ne分別為放電室內(nèi)原初電子濃度、電子總濃度；vp、va分別為放電室內(nèi)原初電子速度、離子速度；fc為離子運動速度因子，與玻姆速度有關；φ為陽極鞘層電勢。

由式（1）、式（2）可知，放電室陽極推進劑利用率及電離損耗是關于離子推力器結(jié)構(gòu)參數(shù)和放電室工作參數(shù)的函數(shù)。對整體結(jié)構(gòu)參數(shù)已確定的10 cm氙離子推力器，其只與放電室工作參數(shù)有關。

通過優(yōu)化陽極推進劑在放電室空腔內(nèi)的分布特性，從而改變陽極推進劑電離位置，實現(xiàn)陽極推進劑在有效駐留時間內(nèi)被碰撞電離幾率的提高，構(gòu)建高效、穩(wěn)定電離條件下的等離子體放電電壓關系，進而實現(xiàn)放電室陽極推進劑利用率的提升，而伴隨著等離子體放電電壓下降的陽極推進劑利用率提升也將聯(lián)動引起離子推力器其他電參數(shù)的變化，使得陽極推進劑電離損耗出現(xiàn)下降。與此同時，實現(xiàn)陽極推進劑在放電室空腔內(nèi)良好的濃度分布與梯度變化也是提升等離子體云宏觀對稱性及其內(nèi)部離子密度均勻性的關鍵，從而為提高柵極組件引出離子束流的平直度奠定基礎，最終為減緩柵極組件的離子刻蝕速率、提高離子推力器在軌工作壽命提供保障。

1.3 陽極推進劑供給與流動特性分析

10 cm氙離子推力器工作過程中，陽極推進劑供給單元通過減壓閥將推進劑氙氣壓強調(diào)節(jié)至0.3 MPa左右，并以該壓強作為推進劑流量調(diào)控的前級供給壓強。此后，流量控制裝置通過三路供氣管路分別實現(xiàn)對中和器、主陰極和放電室陽極的推進劑供給和調(diào)控。10 cm氙離子推力器的放電室陽極推進劑供給流量為0.5 mg/s，經(jīng)φ1 mm內(nèi)徑的供氣管路進入陽極氣體分配環(huán)，并由陽極氣體分配環(huán)將其均分在放電室內(nèi)部。

推進劑氙氣在上述過程的流動特性可根據(jù)Reynold數(shù)（Re）［14］和Knudsen數(shù)（Kn）［15］進行判別。

式中：ρ為推進劑氙氣密度；u為推進劑氙氣流速；d為供氣管路直徑；μ為推進劑氙氣的黏滯系數(shù)［16］，常溫下其黏滯系數(shù)為2.017×10-4g·cm-1·s-1；Q為推進劑氙氣體積流量；P為供氣管路內(nèi)推進劑氙氣壓強；M為推進劑氙氣原子質(zhì)量；A為供氣管路截面面積；k為玻耳茲曼常數(shù)；T為推進劑氙氣溫度；λ為溫度和壓強確定條件下的推進劑氙氣平均自由程；n為單位體積內(nèi)的推進劑氙氣分子數(shù)；σ為推進劑氙氣碰撞截面。當供氣管路內(nèi)推進劑氙氣壓強低于500 Pa時，Kn<0.1，推進劑流動處于連續(xù)介質(zhì)區(qū)及速度滑移區(qū)；而Re<3 000，推進劑流動處于層流或低雷諾數(shù)流動區(qū)域。

以式（3）、式（4）為基礎，結(jié)合10 cm氙離子推力器陽極推進劑供給管路結(jié)構(gòu)參數(shù)，分析發(fā)現(xiàn)：當陽極推進劑流量為0.5 mg/s時，流量控制裝置出口壓強約為400 Pa，此時包括陽極氣體分配環(huán)在內(nèi)的陽極供氣管路中推進劑的流動狀態(tài)為層流，整體處于連續(xù)介質(zhì)區(qū)及速度滑移區(qū)。該分析結(jié)論與Ketsdever［17］及Lilly［18-19］等開展的空間微推力測量中類似尺度的微型管路中不同來流壓強條件下He工質(zhì)氣體流動分析的研究結(jié)果相一致。結(jié)合上述分析，本文采用基于CFD方法的FLUENT［20］開展相似來流條件和幾何尺度下的陽極推進劑動態(tài)力學行為研究。

2 陽極推進劑供給優(yōu)化與分析

2.1 陽極推進劑供給方式設計過程

陽極推進劑供給方式的差異直接體現(xiàn)在陽極環(huán)的基本結(jié)構(gòu)上，為有效分析陽極推進劑供給方式對陽極推進劑利用率的影響，在研究陽極環(huán)結(jié)構(gòu)特征的基礎上，將陽極環(huán)分為進氣管和分配環(huán)2部分（見圖2）。受流量控制裝置調(diào)控的陽極推進劑通過進氣管流入分配環(huán)后，自分配環(huán)上均布的微孔進入放電室。

圖2 陽極環(huán)結(jié)構(gòu)分解示意圖Fig.2 Schematic diagram of structural decomposition of anode-ring

通過調(diào)整進氣管數(shù)量、分配環(huán)數(shù)量及推進劑噴射方向?qū)崿F(xiàn)陽極推進劑供給方式的優(yōu)化?；趥鹘y(tǒng)Kaufman型離子推力器陽極環(huán)結(jié)構(gòu)，從分析方案正確、結(jié)構(gòu)設計合理、制造工藝可行的角度，通過優(yōu)化分配環(huán)和推進劑噴口設計形成了8種可實現(xiàn)推進劑差異化供給的陽極環(huán)：①一個進氣管、一個分配環(huán)且推進劑沿陽極環(huán)軸向直接噴射的陽極環(huán)結(jié)構(gòu)（簡稱單管單環(huán)直噴結(jié)構(gòu)）；②兩個進氣管、兩個分配環(huán)且推進劑向陽極環(huán)內(nèi)側(cè)與外側(cè)雙向側(cè)噴的陽極環(huán)結(jié)構(gòu)（簡稱雙管雙環(huán)側(cè)噴結(jié)構(gòu)）；③單管單環(huán)側(cè)噴結(jié)構(gòu)；④單管雙環(huán)側(cè)噴結(jié)構(gòu)；⑤單管雙環(huán)直噴結(jié)構(gòu)；⑥雙管單環(huán)直噴結(jié)構(gòu)；⑦雙管單環(huán)側(cè)噴結(jié)構(gòu)；⑧雙管雙環(huán)直噴結(jié)構(gòu)（見圖3）。

利用上述不同規(guī)格類型的陽極環(huán)，運用CFD理論，建立了包括陽極推進劑、進氣管和分配環(huán)在內(nèi)的CFD陽極環(huán)模型，研究了陽極推進劑未發(fā)生放電情況下，當保持流道當量面積相同時，陽極推進劑供給方式對其流速與壓強及其在放電室內(nèi)分布的影響作用關系。

2.2 陽極環(huán)內(nèi)陽極推進劑壓強與流速變化

針對10 cm氙離子推力器陽極推進劑供給流量及其入口壓強條件，以陽極環(huán)內(nèi)陽極推進劑流場為基礎，采用CFD陽極環(huán)模型，研究了單管單環(huán)直噴結(jié)構(gòu)、雙管雙環(huán)側(cè)噴結(jié)構(gòu)等8種規(guī)格陽極環(huán)內(nèi)陽極推進劑流速與壓強變化規(guī)律。圖4給出了部分規(guī)格的陽極環(huán)取進氣管所在方位θ=45°、135°、225°及315°共計4個截面處的陽極推進劑壓強與流速分布云圖。

圖3 陽極推進劑供給方案Fig.3 Propellant allocation scheme of anode

圖4 不同規(guī)格陽極環(huán)出口處陽極推進劑流速與壓強分布Fig.4 Velocity and pressure distribution of anode propellant at exit for different types of anode-ring

如圖4所示，不同規(guī)格陽極環(huán)的陽極推進劑壓強與流速變化規(guī)律不同，不同方位角間陽極推進劑壓強與流速分布差異明顯，并均呈現(xiàn)強烈的非線性特征分布。對比圖4（a）、圖4（c）及圖4（b）、圖4（d）可以發(fā)現(xiàn)，雙管結(jié)構(gòu)陽極環(huán)不同方位角間的陽極推進劑壓強分布的均勻性明顯優(yōu)于單管結(jié)構(gòu)陽極環(huán)，這種優(yōu)勢對于單管結(jié)構(gòu)陽極環(huán)尤為明顯。受分配環(huán)管道直徑小、流導低的影響，單管單環(huán)結(jié)構(gòu)陽極環(huán)進氣管附近位置方位角的陽極推進劑壓強明顯高于其他位置；而雙管單環(huán)結(jié)構(gòu)陽極環(huán)在采取兩端進氣后，分配環(huán)內(nèi)部不同方位角的陽極推進劑壓強分布則較為均勻。對于增加一層分配環(huán)的雙環(huán)結(jié)構(gòu)陽極環(huán)，受分配環(huán)內(nèi)陽極推進劑二次分配的影響，不同方位角間的陽極推進劑壓強分布較為均勻，陽極推進劑在分配環(huán)內(nèi)的分布較少受管道流導的影響。與此同時，受壓強分布均勻性變化的影響，管道內(nèi)及出口位置的陽極推進劑流速同樣存在相應的變化。

此外，鑒于雙環(huán)分壓后的陽極環(huán)內(nèi)部陽極推進劑壓強低于單環(huán)結(jié)構(gòu)，其出口處的陽極推進劑流速低于單環(huán)結(jié)構(gòu)。可見，通過增加分配環(huán)，可有效改善陽極推進劑在陽極環(huán)內(nèi)部分布的均勻性，并降低流場不同區(qū)域間的壓強變化梯度，進而影響到不同區(qū)域的陽極推進劑噴射速度。通過采取雙層分配環(huán)，可有效控制陽極推進劑噴射速度，降低了高速氣流進入放電室后尚未電離便快速透過柵極組件所造成陽極推進劑損耗，提高了陽極推進劑在放電室內(nèi)的駐留時間，為實現(xiàn)電子與陽極推進劑的有效碰撞電離奠定基礎。

同時，由圖4（a）、圖4（b）及圖4（c）、圖4（d）發(fā)現(xiàn)，單/雙環(huán)結(jié)構(gòu)所導致的陽極推進劑流速存在差異，不同方位角間陽極推進劑在雙向側(cè)噴結(jié)構(gòu)下的流速一致性優(yōu)于單向直噴結(jié)構(gòu)下的流速一致性。因此，采取雙向側(cè)噴方式可顯著提升陽極推進劑在放電室中心區(qū)域、邊緣區(qū)域及周向整體濃度分布的一致性，進而提高放電室內(nèi)等離子體分布的均勻性，提升陽極推進劑電離效率與放電穩(wěn)定性，降低高能雙荷離子產(chǎn)生幾率，提高離子推力器工作壽命。

基于上述分析，雙環(huán)、側(cè)噴結(jié)構(gòu)陽極環(huán)對改善放電室內(nèi)陽極推進劑分布的均勻性，提升放電室陽極推進劑利用率有著顯著作用。

2.3 供給方式對放電室內(nèi)陽極推進劑分布影響

綜上所述，單/雙管結(jié)構(gòu)的陽極推進劑供給方式對有著陽極推進劑二次分配、調(diào)壓功能的雙環(huán)結(jié)構(gòu)陽極推進劑分布特性影響較小，但雙管結(jié)構(gòu)的陽極環(huán)較單管結(jié)構(gòu)陽極環(huán)將大幅增加離子推力器放電室陽極推進劑供給管路走向設計與結(jié)構(gòu)布局的復雜性，并提高離子推力器整機質(zhì)量。因此，雙環(huán)、側(cè)噴狀態(tài)下的陽極環(huán)單管結(jié)構(gòu)比雙管結(jié)構(gòu)具有更強的結(jié)構(gòu)可實施性和工藝可操作性。

有鑒于此，對分析后優(yōu)選的單管雙環(huán)側(cè)噴結(jié)構(gòu)陽極環(huán)和現(xiàn)階段10 cm氙離子推力器使用的單環(huán)直噴結(jié)構(gòu)陽極環(huán)2種不同供給方式下陽極推進劑在放電室空腔內(nèi)壓強分布變化規(guī)律開展了研究。圖5給出了2種規(guī)格狀態(tài)下的放電室陽極推進劑壓強分布。

如圖5所示，在陽極推進劑供給總量相同的情況下，受供給方式的影響，陽極推進劑在放電室內(nèi)的壓強場分布規(guī)律不同。放電室內(nèi)單管單環(huán)直噴結(jié)構(gòu)陽極環(huán)的進氣管附近位置陽極推進劑壓強明顯高于其他區(qū)域，陽極推進劑在放電室內(nèi)的周向分布存在顯著的不對稱性。區(qū)別與此，雙環(huán)側(cè)噴陽極環(huán)在放電室內(nèi)產(chǎn)生的陽極推進劑壓強周向分布較為均勻，其徑向壓強變化梯度一致性較好。

圖5 兩種不同規(guī)格陽極環(huán)下放電室陽極推進劑壓強分布Fig.5 Pressure distribution of anode propellant in chamber for two types of anode-ring

3 優(yōu)化結(jié)果驗證

針對上述分析，為驗證CFD優(yōu)化分析結(jié)果的正確性和方法的可行性，以圖6所示10 cm氙離子推力器為基礎，在除陽極環(huán)外的離子推力器結(jié)構(gòu)狀態(tài)、設備試驗條件及過程質(zhì)量控制均保持一致的情況下，對離子推力器在單管單環(huán)直噴結(jié)構(gòu)陽極環(huán)和單管雙環(huán)側(cè)噴結(jié)構(gòu)陽極環(huán)2種狀態(tài)下的放電室性能開展了工作性能測試。圖7為2種規(guī)格陽極環(huán)的分配環(huán)實物圖。

圖6 10 cm氙離子推力器產(chǎn)品外觀Fig.6 Product appearance of 10 cm xenon ion thruster

圖7 兩種不同規(guī)格陽極環(huán)的分配環(huán)實物圖Fig.7 Picture of distribution ring for two types of anode-ring

試驗測試在蘭州空間技術(shù)物理研究所0.5 k W級離子推力器性能測試設備中開展，設備真空室為4 m（直段）×φ1.5 m（內(nèi)徑），真空抽氣系統(tǒng)配備2臺對氮氣抽速為1×104L/s的外置式低溫泵和2臺對氙氣抽速為1.49×104L/s內(nèi)置式低溫泵作為設備主泵，粗抽采用螺桿泵、羅茨泵及分子泵聯(lián)動的復合抽氣機組。設備空載抽空24 h后的極限真空度優(yōu)于5×10-5Pa，加載0.983 mg/s氙氣時，工作真空度優(yōu)于1×10-3Pa。

試驗供氣采用設備地面供氣系統(tǒng)，主要由氙氣瓶、調(diào)壓模塊及流量控制模塊組成，主陰極及中和器推進劑流量控制精度為±（0.8%讀數(shù)＋0.2%滿量程），放電室陽極推進劑流量控制精度為±（0.4%讀數(shù)＋0.2%滿量程），所有流量均可在控制器加電1s后實現(xiàn)精確的測量和控制。

試驗供電采用設備地面供電系統(tǒng)，系統(tǒng)屏柵電源信號調(diào)控區(qū)間為0～1.5 k V/0～400 mA，電壓輸出精度為滿量程±0.075%，穩(wěn)定度為±0.01%；陽極電源信號調(diào)控區(qū)間為0～500 V/0～3 A，電壓輸出精度為滿量程±0.075%，穩(wěn)定度為±0.01%；勵磁電源信號調(diào)控區(qū)間為0～100 V/0～2 A，電流輸出精度為滿量程±0.075%，穩(wěn)定度為±0.01%。圖8為0.5kW 級離子推力器性能測試設備。

圖9為采用單管雙環(huán)側(cè)噴結(jié)構(gòu)陽極環(huán)時的離子推力器工作狀態(tài)。表1為2種狀態(tài)下離子推力器工作性能對比。

圖8 0.5 k W 級離子推力器性能測試設備Fig.8 Performance testing equipment for 0.5 kW class ion thruster

圖9 單管雙環(huán)側(cè)噴結(jié)構(gòu)陽極環(huán)下的離子推力器工作狀態(tài)Fig.9 Operation state of propellant distributor of thruster with dual-stage side ejecting anode-ring

表1 兩種不同規(guī)格陽極環(huán)下離子推力器工作性能對比Table 1 Thrust operation perfor mance comparison of two types of anode-ring

如表1所示，對于10 cm氙離子推力器的20 mN設計工作點，在整機功率基本一致的情況下，通過優(yōu)化陽極推進劑供給方式，離子推力器效率和放電室陽極推進劑利用率分別由59.6%、91.7%提升至62.2%、98.4%，電離損耗由277.9 W/A降至241.2 W/A，各項性能指標得到全面提升，在性能滿足要求的前提下，較好地實現(xiàn)了整體優(yōu)化的目標。同時，也充分驗證了CFD計算優(yōu)化結(jié)果的正確性與方法的可行性。

4 結(jié) 論

本文利用CFD理論，研究了未發(fā)生氣體放電情況下，不同供給方式時陽極環(huán)內(nèi)陽極推進劑的壓強與流速變化情況，以及陽極推進劑供給方式對離子推力器放電室空腔內(nèi)陽極推進劑分布特性的影響作用關系，獲得了以陽極推進劑分配方式與噴射方向?qū)崿F(xiàn)最大限度優(yōu)化為基礎的陽極環(huán)結(jié)構(gòu)，有效提升了10 cm 氙離子推力器工作性能。得到如下結(jié)論：

1）離子推力器放電室的陽極供氣中，陽極推進劑在包括陽極氣體分配環(huán)在內(nèi)的陽極供氣管路中的流動狀態(tài)為層流，整體處于連續(xù)介質(zhì)區(qū)及速度滑移區(qū)。

2）采取雙環(huán)分壓、雙側(cè)開孔的陽極環(huán)，可有效降低陽極推進劑進入放電室的速度，陽極推進劑在放電室內(nèi)周向分布的均勻性得到有效改善，基于該結(jié)構(gòu)狀態(tài)的陽極環(huán)，10 cm氙離子推力器陽極推進劑電離損耗由277.9 W/A降至241.2 W/A，放電室陽極推進劑利用率由91.7% 提升至98.4%。

在后續(xù)研究中，將采取等離子在線診斷與磁流體（Magnetohydrodynamics，MHD）數(shù)值分析相結(jié)合的方式開展電磁場作用下的推進劑中性原子分布對電離過程的影響機理研究，進一步優(yōu)化放電室供氣及磁場，從而提高工作性能。