Open Access
Issue
JNWPU
Volume 42, Number 5, October 2024
Page(s) 793 - 800
DOI https://doi.org/10.1051/jnwpu/20244250793
Published online 06 December 2024

© 2024 Journal of Northwestern Polytechnical University. All rights reserved.

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

卫星运行于超低地球轨道(very low Earth orbit, VLEO, < 250 km)具有发射成本低、通信损耗小和对地观测精度高等优势[13]。然而, VLEO面临着较大的大气阻力, 维持轨道需要不断地消耗推进剂。这限制了卫星在该轨道的长期驻留[4]。为解决这一问题, 使用轨道稀薄大气作推进剂的吸气式电推进(atmosphere breathing electric propulsion, ABEP)技术被提出。VLEO大气主要由氮分子和氧原子组成, 氧对阴极的腐蚀是ABEP研究的一个关键挑战[5]。电子回旋共振(electron cyclotron resonance, ECR)离子源采用无碰撞的ECR机制产生等离子体, 无需热阴极, 有应用于ABEP的潜力[6]。在前期研究中, 为降低推进剂供给压强, 本课题组提出了一种栅格进气道ECR离子源[7]。研究表明, 该进气道的推进剂供给压强仅为传统小孔进气道的0.26。

流量对ECR离子源的离子束流影响显著。国内外ECR离子源的研究经验表明存在临界流量[810], 当流量低于临界流量(简称“临界流量前”)时, 离子束流随流量的增加而快速增加; 当流量高于临界流量时(简称“临界流量后”), 随流量增加, 离子束流会趋于饱和或快速下降。Tani等[9]通过磁路优化提高了临界流量, 进而将离子源的束流提高了16%。Tsukizaki等[1112]通过光学探针诊断了不同流量下离子源轴线处的电场和激发态氙密度的分布。研究表明, 临界流量后, 圆波导内电场和激发态氙密度高于临界流量前; 等离子体对微波的衰减可能是限制离子束流增长的主要因素。

上述研究仅限于氙工质ECR离子源, 且激发态氙的密度分布与等离子体的密度分布并不完全等同, 但其为揭示限制离子束流增长的内在机理指明了方向——等离子体诊断。Langmuir探针可应用于低气压、非平衡等离子体的诊断, 能够在宽泛的实验条件下获得局部等离子体的参数分布。Funaki等[13]使用Langmuir探针诊断了氙工质ECR离子源内的等离子体。但是, 受限于磁场对诊断结果的影响, 仅得到了离子饱和电流和悬浮/空间电势。保持探针轴线与磁力线垂直, 可将磁场对诊断结果的影响降到最低[1416]。基于这一理论, 文献[1718]采用弯曲的Langmuir探针研究了磁路对2 cm和10 cm ECR离子源内等离子体分布的影响。

针对VLEO大气的主要成分氮气, 基于前述栅格进气道ECR离子源, 本文开展其在临界流量前后的弯曲Langmuir探针诊断实验, 并进一步探究了临界流量前后限制ECR离子源离子束流增长的主要原因。

1 ECR离子源的结构和Langmuir探针诊断实验

1.1 ECR离子源的结构

图 1所示, 栅格进气道ECR离子源主要由进气道、进气栅格、圆波导、球形天线、放电室、内外磁环和栅极系统组成。进气栅格安装于进气道和圆波导之间, 兼顾进气和防止微波能量泄漏的功能[7]。球形天线安装于圆波导内部, 微波从此处注入离子源。内、外磁环嵌入放电室内壁, 它们在离子源内产生一定的磁场, 形成ECR区。实验使用的微波频率为4.2 GHz, 对应ECR区的磁通密度为0.15 T。栅极系统由屏栅和加速栅组成, 具体可见文献[7]。

thumbnail 图1

栅格进气道ECR离子源

1.2 Langmuir探针诊断实验

1.2.1 诊断实验系统

图 2a)所示, 诊断实验系统主要由固态微波源、环形器、衰减器、功率计、ECR离子源、Langmuir探针、2D坐标架和扫描电源组成。实验前对离子源的入射微波功率进行校准, 反射功率通过功率计测量[8]。将探针插入带有高电压的栅极系统来诊断离子源内的等离子体是困难的[13]。因此, 实验在不引出离子束流(栅极接地)的条件下进行, 并采用如图 2b)所示的开槽栅极代替栅极系统以便探针通过和移动。开槽栅极的结构参数如表 1所示。如图 2b)所示, Langmuir探针由陶瓷管和钨丝组成, 电流收集段钨丝的直径和长度分别为0.2和7 mm; 2D坐标架控制探针在x, z方向逐步移动。诊断点位于如图 2c)所示网格的格点处, x, z方向相邻诊断点的间距均为5 mm。图中的红色线和带箭头的黑色线分别代表ECR区、磁力线, 其分布通过COMSOL Multiphysics软件计算获得。

图 3为典型的Langmuir探针I-V曲线。该曲线可被探针悬浮电势Vf、等离子体电势Vp分为三段。其中, Vf为收集电流为0的点; Vp为电子电流段曲线的拐点, 可通过计算Ⅰ-Ⅴ曲线的二阶导数确定[19]。根据轨道运动限制方法, 离子密度和电子温度可分别通过(1)~(2)式获得[2022]。

式中, Ap, e, k, M, IiIt分别为探针表面积、电子电荷量、玻尔兹曼常数、离子质量、离子饱和区的收集电流、过渡区域的收集电流。

thumbnail 图2

诊断实验系统

表1

开槽栅极的结构参数

thumbnail 图3

Langmuir探针的I-V曲线

1.2.2 诊断实验工作参数设置

屏栅电压1 350 V、加速栅电压-350 V时, 不同微波功率下, 离子束流随流量的变化如图 4所示。

Pi, q分别为微波功率和推进剂流量。当Pi=25 W, q=3 mL/min时, 离子源可以正常启动, 但无法维持设定电压下的离子引出, 因而离子束流和推进剂利用率记为0。不同微波功率下, 离子束流和推进剂利用率随流量的变化趋势相近。当q≤4 mL/min时, 离子束流和推进剂利用率随着流量的增加而增加; 当q>4 mL/min时, 随着流量增加, 离子束流缓慢增加、趋于饱和, 而推进剂利用率则快速下降。综上所述, 氮工质栅格进气道ECR离子源的临界流量在4 mL/min附近。

诊断实验工作参数如表 2所示。

实验在不引出离子束流的条件下进行, 并使用开槽栅极代替栅极系统。因此, 相同流量下诊断实验和引出实验所对应的离子源内氮气分子数密度nn并不相等。为了验证诊断实验设置的流量是否落在预期区间, 即3, 6和9 mL/min分列临界流量工况两侧, (3)~(4)式被联立以求得不同流量对应的nn[23], 计算结果如图 5所示。

式中, Nin, vn, Ib, Ag, C, L, H, A分别为单位时间进入离子源的分子数、分子热运动速度、离子束流、栅格开孔面积、栅格流导、栅孔周长、栅孔厚度、栅孔面积。

图 5中2条蓝色虚线从上到下依次为图 4Pi=25 W, Pi=35 W时临界流量4 mL/min所对应的nn。从图中可以看出, 诊断实验使用3 mL/min流量时, 其nn始终处于2条蓝色虚线的下方, 离子源处于临界流量前的状态; 使用6, 9 mL/min流量时, 其nn始终处于2条蓝色虚线的上方, 离子源处于临界流量后的状态。综上所述, 诊断实验的流量设置与表 2中的预期相符。

thumbnail 图4

离子束流和推进剂利用率随推进剂流量的变化

表2

诊断实验工作参数

thumbnail 图5

离子源内氮气分子数密度随推进剂流量的变化

2 结果与讨论

2.1 诊断实验结果

Pi=35 W时, 离子源的诊断结果如图 6所示。图中酒红色线所包围的区域为寻常(ordinary, O)波截止区, 其截止条件为 。其中, ω, ε0, me分别为微波频率、真空介电常数和电子质量。受磁场影响, 从Langmuir探针的I-V曲线获得准确的电子数密度是困难的[14]。因此, 在计算中使用准中性假设, 以离子数密度代替电子数密度。由于诊断区域远离壁面, 这一假设是合理的, 并不会给截止区的计算带来显著差异。对于右手螺旋偏振波和异常波的截止区, 计算发现这两者随流量变化并不显著, 因此图中并未给出。

图 6可以发现, 不同流量下, 离子和热电子主要分布于与ECR区相交且磁力线闭合的磁场区域, 这个区域可定义为“主等离子体产生区(main plasma production aera, MPPA)”。在MPPA中, 电子受磁镜约束, 可以多次穿越ECR区获得能量[18, 24]。同时, ECR离子源中的电离主要由热电子引发[25], 所以离子和热电子的分布合理。从图 6中还可以发现, 电子温度随流量增加呈现逐渐下降的趋势, 其中MPPA电子温度的下降速度远高于其他区域。O波截止区以及离子密度随流量增加的变化趋势是非单调的。当流量由3 mL/min增加至6 mL/min, O波截止区由栅极一侧几乎扩展至整个MPPA; 离子密度呈现爆发性增长, 密度峰值由3.38×1017 m-3增长至5.35×1017 m-3。当流量进一步增加至9 mL/min, O波截止区向内磁环一侧收缩; 离子密度小幅下降, 峰值降至4.70×1017 m-3

Pi=25 W时, 离子源的诊断结果如图 7所示。当q=3 mL/min时, 离子源内离子密度峰值仅为1.62×1017 m-3, 不足以截止O波; 蓝色线包围区域电子温度超出氮气电离能(15.58 eV), 这表明当前流量下氮气密度是限制等离子体产生的主要因素。随流量增加, 离子源内等离子体参数分布和O波截止区的变化规律与Pi=35 W的工况基本一致。结合图 6图 7, 可以发现, 对于MPPA以外的区域, 高、低功率离子源电子温度无明显差异; 在MPPA中, 低功率离子源电子温度明显更高, 且更高的区域主要分布在两者O波截止区覆盖存在差异的区域。

thumbnail 图6

Pi=35 W时离子源诊断实验结果

thumbnail 图7

Pi=25 W时离子源诊断实验结果

2.2 讨论

离子束流取决于离子源的引出能力和等离子体产生能力。其中, 引出能力由推力器头部尺寸和栅极电压共同决定, 流量对其影响可以忽略不计。因此, 流量对离子束流的影响等同于其对等离子体产生能力的影响。等离子产生能力可通过电离反应速率来评估。电离反应速率定义为

式中: Ne, σion, cref分别为电子数目、电子和氮气分子的电离碰撞截面、电子和氮气分子的相对运动速度; σion, cref仅为电子温度Te的函数。流量对等离子产生能力的影响可分解为其对Ne, nn, Te 3个参数的影响。

达到临界流量前, 限制离子源等离子体产生的主要因素为nn。如图 4所示, 临界流量前, 推进剂利用率随流量增加而增加。如图 6~7所示, 临界流量前, Pi=35 W和Pi=25 W的工况中, MPPA内大部分区域的电子温度在12~14 eV、12~16 eV范围内。假设电子温度服从麦克斯韦分布, 则这些区域能量超出氮分子电离能的电子占比分别为0.27~0.34, 0.27~0.40。这两者表明, 临界流量前的等离子体产生中, 氮气分子密度过低。

达到临界流量后, 限制离子源等离子体产生的主要因素为Te。如图 4所示, 临界流量后, 推进剂利用率随流量增加而下降。如图 6~7所示, 临界流量后, MPPA中离子密度骤升, 电子温度骤降。这两者表明, 临界流量后的等离子体产生中, 氮气分子和电子充足, 但电子温度过低。

流量对电子温度的影响主要体现在电子能量的损失与获得。一方面, 随着流量增加, 电子和氮气分子的各类碰撞加剧, 电子能量损失频率升高。另一方面, 随着流量增加, 达到临界流量前等离子体产生的限制被解除, 离子密度快速增加, 并在离子源内形成截止区阻碍电子能量的获得。这也解释了MPPA(O波截止区主要存在区域)中更大的电子温度下降幅度, 以及低功率离子源在高低功率离子源O波截止区覆盖存在差异区域更高的电子温度。

令人困惑的是, O波是一种微波电场与静磁场方向平行的波, 并不直接参与电子的加热。为验证O波截止区对电子温度的影响是否有效, 有必要计算不同工况下的电子获能。电子单次穿越ECR区获得的平均能量可通过 计算, E, v, 分别为电场垂直磁力线方向的分量、电子速度平行磁力线方向的分量、ECR区磁场关于磁力线的梯度。仅与离子源结构和工作参数有关, 常用作衡量电子获能的指标。其中, E可基于冷等离子体假设采用文献[18]的计算方法得到。

电子获能指标计算结果如图 8所示。临界流量前, 低功率工况, 不存在O波截止区, 内外磁环中间区域电子获能指标最高; 高功率工况, 不存在O波截止区的内磁环区域电子获能指标最高。当流量增加至6 mL/min时, 高、低功率离子源中, O波截止区大幅扩张, 电子获能指标平均值大幅下降, 分别由1.53×108, 1.92×108降至2.75×107, 3.32×107。当流量进一步增加至9 mL/min时, 高、低功率离子源中, O波截止区小幅缩小, 电子获能指标平均值小幅上升, 分别升至4.06×107, 8.97×107。上述结果表明, O波截止区对电子获能指标的影响显著, 区域的扩张和收缩分别对应电子获能指标平均值的下降和升高。电子获能指标决定了电子的能量获取效率, 对Te的影响重大。因此O波截止区对Te的影响是切实有效的。

thumbnail 图8

不同流量下的电子获得指标对比

3 结论

本文采用弯曲的Langmuir探针, 测量了临界流量前后氮气工质栅格进气道ECR离子源离子密度和电子温度的分布, 并分析了限制离子束流增长的主要因素。结论如下:

1) 离子和热电子主要分布于MPPA, 电子温度随着流量的增加而逐渐下降, MPPA中电子温度的下降幅度远高于其他区域。达到临界流量前, 随着流量增加, 离子密度爆发性增长, O波截止区大幅扩张; 临界流量后, 随着流量增加, 离子密度下降, O波截止区小幅收缩。

2) 临界流量前后限制离子束流增长的主要因素分别为中性气体密度和电子温度。限制因素的不同是造成离子束流随流量变化趋势改变的根本原因。O波截止区通过阻碍电子能量获得影响电子温度, 这一机制主要作用于临界流量后。此时, 对应O波截止区的扩张, 代表电子能量获取效率的电子获能指标相比临界流量前下降了一个数量级。

References

  1. BERTOLUCCI G, BARATO F, TOSON E, et al.. Impact of propulsion system characteristics on the potential for reduction of Earth observation missions at very low altitudes cost[J]. Acta Astronautica, 2020, 176: 173–191 [Article] [Google Scholar]
  2. CRISP N H, ROBERTS P C E, LIVADIOTTI S, et al. The benefits of very low Earth orbit for Earth observation missions[J]. Progress in Aerospace Sciences, 2020, 117: 100619 [Article] [CrossRef] [Google Scholar]
  3. JOSEP V L, ROBERTS P C E, ZHOU H, et al. Very low Earth orbit mission concepts for Earth observation benefits and challenges[C]//Proceedings of the 12th Reinventing Space Conference, London, 2014 [Google Scholar]
  4. DRINKWATER M R, FLOBERGHAGEN R, HAAGMANS Ret al. GOCE: ESA’s first Earth explorer core mission[J]. Space Science Reviews, 2003, 17: 419–432 [CrossRef] [Google Scholar]
  5. ZHENG P, WU J J, ZHANG Y, et al. A Comprehensive review of atmosphere-breathing electric propulsion systems[J]. International Journal of Aerospace Engineering, 2020, 2020: 8811847 [Google Scholar]
  6. NISHIYAMA K. Air breathing ion engine concept[C]//Proceedings of the 54th International Astronautical Congress of the International Astronautical Federation, Bremen, 2003 [Google Scholar]
  7. YANG Juan, MOU Hao, GENG Hai, et al.. Calculation study of electromagnetic fields and neutral particles of a wide inlet electron cyclotron resonance ion source[J]. Journal of Propulsion Technology, 2023, 44(7): 2208097 (in Chinese) [Google Scholar]
  8. JIN Yizhou, YANG Juan, FENG Bingbinget al. Ion extraction experiment for electron cyclotron resonance ion source with different magnetic topology[J]. Acta Physica Sinica, 2016, 65(4): 045201 (in Chinese) [CrossRef] [Google Scholar]
  9. TANI Y, TSUKIZAKI R, KODA D, et al. Performance improvement of the μ10 microwave discharge ion thruster by expansion of the plasma production volume[J]. Acta Astronautica, 2019, 157: 425–434 [Article] [Google Scholar]
  10. BARQUERO S, TABATA K, TSUKIZAKI R, et al.. Performance characterization of the μ10 electron-cyclotron-resonance ion thruster using alternative propellants: Krypton vs Xenon[J]. Acta Astronautica, 2023, 211: 750–754 [Article] [Google Scholar]
  11. TSUKIZAKI R, KOIZUMI H, NISHIYAMA K, et al. Measurement of axial neutral density profiles in a microwave discharge ion thruster by laser absorption spectroscopy with optical fiber probes[J]. Review of Scientific Instruments, 2011, 82(12): 123103 [Article] [CrossRef] [Google Scholar]
  12. TSUKIZAKI R, KOIZUMI H, NISHIYAMA K, et al. Plasma diagnostics inside the ECR ion thruster μ10 by laser absorption spectroscopy with optical fiber probes[C]//Proceedings of the 48th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Atlanta, 2012 [Google Scholar]
  13. FUNAKI I, KUNINAKA H, TOKI K. Plasma characterization of a 10 cm diameter microwave discharge ion thruster[J]. Journal of Propulsion and Power, 2004, 20(4): 718–727 [Article] [CrossRef] [Google Scholar]
  14. YANG J, FU Y L, LIU X C, et al. Bended probe diagnostics of the plasma characteristics within an ECR ion source with a rectangular waveguide[J]. Plasma Science and Technology, 2018, 20(8): 085402 [Article] [CrossRef] [Google Scholar]
  15. DEMIDOV V I, RATYNSKAIA S V, RYPDAL K. An electrostatic probe with end plugs for measurements in a magnetized plasma[J]. Review of Scientific Instruments, 1999, 70(11): 4266–4269 [Article] [NASA ADS] [CrossRef] [Google Scholar]
  16. SHIKAMA T, KITAOKA H, HASUO M. Detection of electron energy distribution function anisotropy in a magnetized electron cyclotron resonance plasma by using a directional Langmuir probe[J]. Physics of Plasmas, 2014, 21(7): 073510 [Article] [NASA ADS] [CrossRef] [Google Scholar]
  17. XIA X, YANG J, JIN Y Z, et al. The influence of magnetic circuit and operating parameters on the plasma property of 2 cm ECRIT ion source[J]. Vacuum, 2020, 179: 109517 [Article] [CrossRef] [Google Scholar]
  18. MOU H, JIN Y Z, YANG J, et al. Interaction between plasma and electromagnetic field in ion source of 10 cm ECR ion thruster[J]. Chinese Physics B, 2022, 31(7): 075202 [Article] [CrossRef] [Google Scholar]
  19. LIEBERMAN M A, LICHTENBERG A J. Principles of plasma discharges and materials processing[M]. New Jersey: John Wiley & Sons, 2006 [Google Scholar]
  20. DEMIDOV V I, RATYNSKAIA S V, RYPDAL K. Electric probes for plasmas: the link between theory and instrument[J]. Review of Scientific Instruments, 2002, 73(10): 3409–3439 [Article] [NASA ADS] [CrossRef] [Google Scholar]
  21. CHEN F F. Langmuir probe analysis for high density plasmas[J]. Physics of Plasmas, 8(6): 3029–3041 [Article] [Google Scholar]
  22. JIN Yizhou, YANG Juan, LUO Litao, et al. Diagnosis of electron cyclotron resonance plasma inside ion source[J]. High Voltage Engineering, 2015, 41(9): 2950–2597 (in Chinese) [Google Scholar]
  23. GOEBEL D M, KATZ I. Fundamentals of electric propulsion: ion and hall thrusters[J]. New Jersey: John Wiley & Sons, 2008 [CrossRef] [Google Scholar]
  24. CORAL G, TSUKIZAKI R, NISHIYAMA Ket al. Microwave power absorption to high energy electrons in the ECR ion thruster[J]. Plasma Sources Science and Technology, 2018, 27(9): 095015 [Article] [CrossRef] [Google Scholar]
  25. YAMASHITA Y, TANI Y, TSUKIZAKI Ret al. Numerical investigation of plasma properties for the microwave discharge ion thruster μ10 using PIC-MCC simulation[J]. Physics of Plasmas, 2019, 26: 073510 [Article] [NASA ADS] [CrossRef] [Google Scholar]
  26. XIA Xu, YANG Juan, JIN Yizhou, et al. Experimental study of magnetic circuit and antenna position influence on performance of 2 cm electron cyclotron resonance ion thruster[J]. Acta Physica Sinica, 68(23): 235202 (in Chinese) [Google Scholar]

All Tables

表1

开槽栅极的结构参数

表2

诊断实验工作参数

All Figures

thumbnail 图1

栅格进气道ECR离子源

In the text
thumbnail 图2

诊断实验系统

In the text
thumbnail 图3

Langmuir探针的I-V曲线

In the text
thumbnail 图4

离子束流和推进剂利用率随推进剂流量的变化

In the text
thumbnail 图5

离子源内氮气分子数密度随推进剂流量的变化

In the text
thumbnail 图6

Pi=35 W时离子源诊断实验结果

In the text
thumbnail 图7

Pi=25 W时离子源诊断实验结果

In the text
thumbnail 图8

不同流量下的电子获得指标对比

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.