-
In view of the application requirements of electric propulsion system for China’s asteroid deep space exploration mission, the overall scheme is designed, an ion thruster prototype model is established by using a four-ring-cusp field discharge chamber, 30-cm beam current extraction diameter three-grid ion optics system. Reasonableness and compatibility of discharge chamber and grid design are verified experimentally and theoretically . The test results are shown below. The ion thruster can operate steadily over an input power envelope of 277–3120 W, thrust increases linearly from 9.9 to 117.2 mN, specific impulse rises from 1269 to 3492 s, the beam divergence angle drops from 30.7° to 26.8° and stabilizes above a certain power value, the thrust vector angle is less than 1.5° and beam flatness parameter is greater than 0.75 at different operating points. The maximum percentage reduction in grid gap aberration is 90% with the strain relief molybdenum mounting ring thermal design. This research provides a reference for multi-mode ion thruster design and in-orbit engineering applications.
1. 引 言
离子电推进因其高比冲、长寿命、推力宽范围精确可调等优势成为先进卫星平台的标志, 为航天器提供动力, 在星际深空探测领域, 多模式离子推力器更是受到各国的普遍重视并得到工程应用[1-3]. 美国DS-1和DAWN号深空探测器配备NSTAR-30推力器, 口径为30 cm, 放电室设计采用三极环型会切场, 离子光学系统采用凸面钼双栅, 推力在20—92 mN宽范围可调[4-6], 随后研发36 cm NEXT推力器, 放电室采用多极环切场, 并采取多余物抑制措施, 离子光学系统为钼双栅结构, 推进剂采用放电室反向供给, 推力在26—236 mN宽范围可调, 比冲最高4200 s[7]. 同时针对深空探测美国也开发了口径为57 cm的NEXIS推力器, 放电室采用六极环切场, 离子光学系统采用碳碳复合材料, 功率覆盖范围13—28 kW[8]. 22 cm T6推力器采用发散场磁场构型放电室, 钡钨空心阴极, 变孔径钼屏栅+碳加速栅凹面双栅离子光学系统, 在2.43—4.50 kW功率下, 推力75—145 mN, 比冲最高达到4120 s[9].
国内外学者就多模式离子推力器等离子体运动、碰撞、激发、电离、输送、约束等开展了大量理论研究[10-14], 放电室仿真模拟方面, Brophy和Wilbur [15]建立了描述放电室性能的零维模型, Arakawa和Wilbur [16]与Mahalingam和Menart [17]基于PIC-MCC算法通过二维模型探究了磁场结构对放电室等离子体的约束. Wang等[18]、陈茂林等[19]、龙建飞等[20]以及Wirz 和Goebel[21]通过仿真研究了多模式离子推力器离子光学系统设计及其对推力器的影响, 多模式离子推力器在国外已实现在轨应用, 公开研究成果聚集于在轨表现, 地面磨损测试等方面, 但对离子推力器多模式化研究、详细工程设计报道较少.
我国研究起步较晚, 兰州空间技术物理研究所研制的单工况点LIPS-200、双工况点LIPS-300推力器适用于地球同步轨道卫星, 分别于2017年、2020年实现在轨工程应用[22-24], 针对我国首次小天体探测任务探测器与太阳距离1—3.5 AU大范围变化需求, 设计了300—3000 W多工作点电推进系统, 为小天体探测器从地球向主带彗星转移段提供必要的推力 [25].
本文基于我国首次小天体探测任务对推力器的实际需求, 提出总体方案设计, 开展放电室宽范围适应性、离子光学系统屏栅变孔径及热稳定性设计, 搭建性能测试、热稳定性和束流分布特性测试平台, 验证推力器放电室、离子光学系统设计的合理性及匹配性, 本文的研究对离子推力器多模式化设计及工程化样机研制有重要意义.
2. 多模式离子推力器设计
离子推力器多模式设计要求提供原初电子的主阴极具有宽范围电子发射能力, 放电室兼顾宽范围适应性和小功率放电稳定性, 离子光学系统要解决大功率模式下栅极中心束流欠聚焦和小功率模式下栅极边缘束流过聚焦问题, 最终通过放电室与离子光学系统匹配设计实现较高的束流平直度, 离子推力器在结构、几何参数一定的条件下, 通过电气参数的控制实现多模式调节和稳定输出.
2.1 总体方案设计
推力器总体方案设计: 为了提升放电室放电效率、改善放电室下游束流分布特性, 放电室采用四极环形会切场构型, 永磁铁为耐高温低磁损钐钴材料, 空心阴极采用20 A全石墨触持极六硼化镧发射体方案, 离子光学系统采用高引出、低溅射率三栅极技术方案以提高工作寿命, 并采用应力释放结构及钼材料安装环增加了热稳定性, 供电采用高低压分组及大电流供电内引线三线并联技术方案, 总体承力结构采用主支撑环集成式框架设计.
2.2 放电室
放电室是产生和约束等离子体的场所, 小天体探测多模式离子推力器放电室的设计需要重点考虑功率宽范围适应性和低功率下放电稳定性.
束电流口径设计是推力器设计的初始条件, 由Child-Langmuir方程可得栅极组件的最大束电流密度为[26]
Jmax=49ε0√2emiV3/2tl2e, (1) 式中,
ε0 为真空介电常数;e 为元电荷量;mi 为工质气体离子质量;le 为栅极有效栅间距,le=√(lg+ts)2+d2s4 ,lg 为栅极间距,ts 为屏栅厚度,ds 为屏栅孔直径;Vt 为屏栅与加速栅之间的电压.考虑推力器可靠性设计裕度, 取极限束电流密度的80%为推力器工作最大电流密度, 则束电流口径:
Db=2√Ibmax(Jmax×80%)⋅π. (2) 推力器最大工况下
Ibmax=2.1 A , 结合下文栅极设计参数, 束电流口径需大于24.6 cm, 考虑设计裕度及栅极抗力学特性, 栅极口径设计为30 cm.其中放电室口径
Dd 取值一般较束电流口径略大[27], 取:Dd=(1—1.2)Db. (3) 放电室长度设计需要兼顾工质利用率和电离效率, 长径比
Ld/Dd 一般在0.4—1范围内分布, 放电室几何构型采用直段阳极筒+锥段阳极筒的方式, 这主要是为了增加原初电子的约束长度, 提高推力器效率和束流均匀性, 减少对栅极的溅射刻蚀, 提升推力器寿命.多模式离子推力器整体结构设计参数来源于最大功率约束条件下束流口径, 极限束电流密度等. 最小功率点下, 离子推力器放电室内工质气体流率较小, 为了兼顾最大功率点设计要求, 放电室几何参数较大, 小流率下放电室中性气体浓度低, 当阴极发射电流一定, 放电室等离子体密度低容易引发放电不稳定.
因此主阴极发射电流和工质流量大小也是影响推力器小功率下放电稳定性的重要因素, 通过增加放电室气体工质流量, 牺牲推力器工质利用率, 增大主阴极触持电流牺牲效率可提高多模式离子推力器小功率下放电稳定性. 但为了提升推力器性能, 拓宽推力器多模式调节能力, 需要开展磁场的优化设计.
多模式离子推力器放电室磁场由四组磁环组成, 如图1所示, 分别为上磁极、中磁极、锥段磁极和下磁极, 通过偶数的四极场可以使靠近栅极的区域趋向为无场区, 使得喷口处等离子体密度均匀性得到改善, 满足功率宽范围多工作点输出要求.
环形会切场放电室将摆放永磁铁的四周打开, 南北极方向交替的永磁体组成的磁场结构如图2所示, 从图可以看出, 永磁铁会形成尖头形磁力线结构, 靠近永磁铁的磁感应强度有一个最大值, 并沿着腔内方向变小, 这种结构使得放电室等离子体空间的大部分为无磁场区.
其中放电室磁场强度的空间分布为[26]
By(x,y)=πwBo2dcos(πxd)e−πy/d. (4) By 是磁极表面的磁场强度, d是磁极间的距离,w 是磁极的宽度,x 是沿着阳极壁面距离磁极的距离,y 是与阳极壁垂直的方向.磁铁间中心线总的磁场强度为
B(y)=2(2qδ)yd(d24+y2)2=2Myd(d24+y2)2. (5) q 是磁偶极子数量,δ 是磁铁高度的一半,M 是每个磁铁单位长度的磁偶极子磁场强度, 等于磁偶极子的数量乘以磁铁长度:M=BrVm4πw, (6) 其中
Br 是磁铁剩磁,Vm 是磁铁体积.由上式可知, 当
y=d2√3 时, 磁场强度最大, 为B(ymax)=5.2Md2. (7) 放电室磁场设计要兼顾放电高效性、均匀性、小功率下稳定性, 磁场强度高可提升推力器效率, 但太高会导致阳极吸收面积较小, 阳极面积小到一定值时等离子体电位将变为负值, 放电电压一定, 负等离子体电位降低了原初电子能量, 电子、原子的碰撞频率降低, 电离率下降, 放电室在小功率下放电震荡变大, 甚至导致推力器熄弧.
离子推力器放电电流:
Ia=14(8kBTπm)1/122eneAaexp(−eϕ/kBTe), (8) 其中kB为玻尔兹曼常数, 1.3807 × 10–23J/K ;
Te 为电子温度, eV ; m为电子质量, kg; e为电子电量, 1.602 × 10–19 C;ne 为等离子体电子密度, m–3; Aa为阳极表面能吸收电子和离子的面积大小, m2;ϕ 为阳极的等离子电势, V.Aa与电子和离子的拉莫尔半径有关:
Aa = 4rhLc = 4√reriLc, (9) 其中rh为混合拉莫尔半径, m; Lc为磁极的总长度, m; re为电子拉莫尔半径, m; ri为离子拉莫尔半径, m.
从(9)式可以看出, 随着磁场强度减小, 阳极表面能吸收电子和离子的面积Aa增大, 放电电流增大, 即放电电流与磁场强度成反比. 为了避免小功率下推力器出现放电不稳定甚至熄弧现象, 通过减弱放电室磁场强度可增加磁空区, 进而提高放电均匀性, 同时由于磁场对原初电子运动约束能力的降低, 放电电流增加, 通过减弱磁场与放电电流的设计匹配, 可满足小天体多模式离子推力器低功率下稳定放电需求.
基于离子推力器口径和放电室几何构型, 尺寸, 磁场设计中磁极数量和强度的选择要考虑将磁场控制在阳极表面区域, 增大放电室内部的磁空区, 迭代优化后放电室内部磁感应强度分布如图3所示, 兼顾宽范围放电室均匀性和小功率下放电稳定性, 多模式离子推力器放电室设计参数如表1所列.
表 1 离子推力器放电室设计参数Table 1. Design parameters of the ion thruster discharge chamber.几何构型 直段 + 锥段阳极筒 放电室口径和栅极直径比 1—1.2 长径比 0.4—1 磁极数 4 无量纲后磁极间距离 0.6—0.85 闭合磁等势线/Gs 50—60 磁体体积宽度比/cm2 0.4—0.8 永磁体剩磁/Gs ≥9500 永磁体矫顽力/(kA·m–1) ≥700 2.3 离子光学系统
离子光学系统其引出性能取决于本身的设计特性和与上游放电室等离子体分布的匹配性. 离子光学系统引出束流因空间电荷效应限制而遵守Child-Langmuir方程, 一定总加速电压下, 用单孔导流系数p表征离子光学系统的引出能力
p=ibV3/2t=πε09√2emid2sl2e. (10) ib 为单孔束流强度, 因此, 要使离子光学系统获得较大的导流系数就需要栅极的孔径大于栅间距, 理想的离子光学系统要求是束流密度高、发散角小、均匀度高, 选取合适的栅极孔径和栅间距离可以获得适宜的导流系数, 从而获得较好的离子引出, 减小对栅极的轰击.凹面栅极会挤占放电室空间, 沿栅极表面想要设计成无场区较困难, 束流平直度较差, 进而引起栅极温度梯度增加, 溅射加剧, 同时选择球面凸栅可保证栅极在宽范围引出热循环工况下形变方向的一致性, 因此选择面向小天体探测的多模式离子推力器采用凸面三栅结构.
束流均匀性直接关系着推力器的寿命和可靠性, 为了提升推力器束流均匀性, 采用变孔径屏栅孔径
ds(r) 与等离子体密度ni(r) 和电子温度Te(r) 分布相匹配, 即屏栅上不同径向位置的屏栅孔径按该位置的ni 和Te 来确定, 单孔束流ib 与等离子体密度ni 、电子温度Te 、损失因子exp(−2ts/2tsdsds) 和屏栅孔径ds 各个参数的关系为[27]ib∝ni√Teexp(−2tsds)d2s. (11) 将屏栅进行分区, 然后计算出各区域
ni 和Te 的相对平均值, 若将(11)式指数项展开为级数, 可得exp(−2tsds)=1−2tsds+2t2sd2s−⋯. (12) 一般
tsds 取0.15—0.2,2t2sd2s≈0.06 , 一级近似可忽略二次方以上项, 根据(11)式和(12)式可得ib,kib,j=ni,k√Te,k(1−2tsds,k)d2s,kni,j√Te,j(1−2tsds,j)d2s,j, (13) 式中, 下标k代表k区域, 下标j代表j区域. 根据各区域内单孔引出束流相同, 即
ib,k=ib,j , 确定不同区域的屏栅孔径ds(r) , 参数如表2所列, 通过不同的孔径引出相同的束流来改善束流密度均匀性, 其示意图如图4所示.表 2 离子光学系统屏栅变孔径分区及归一化后参数Table 2. Normalized parameters for screen grids with variable aperture zones.分区 分区范围 孔径 0—OA 0—3 1 OA—OB 3—5 1.05 OB—OC 5—12 1 OC—OD 12—15 0.94 推力器热态点火时, 栅极材料受热膨胀造成屏栅、加速栅和减速栅之间的间距沿着半径方向发生变化, 由于各电极间厚度、几何透明度、能量沉积的差异性, 三栅极工作时各电极的变形速率和变形量不同, 栅间距分布不一致使得导流系数发生改变, 栅孔对中性差使得离子运动轨迹发生偏移, 影响束流均匀性和推力器寿命. 栅极组件初始设计采用无应力释放钛合金支撑环方案, 点火中推力器双栅间束流闪烁较多, 在高功率点火后甚至出现栅片边缘翘曲现象, 如图5所示.
为了满足多模式离子推力器栅极组件宽范围束流引出需求, 提高束流均匀性, 延长寿命. 采取以下措施提高栅极组件热稳定性: 1) 将栅极组件安装环设计在推力器放电室外部, 以降低安装环处温度, 减小放电室对栅极组件边缘的热影响; 2) 栅极安装环采用钼材料代替之前的钛合金材料, 使工作时安装环与栅片间的热变形协调一致; 3) 栅极安装环设计应力释放槽, 提高栅极在高温下的形变一致性, 图5给出了优化前采用无应力释放钛安装环和优化后采用应力释放钼安装环栅极边缘变形示意图, 可以看出优化后推力器高温工况下栅片边缘无翘曲, 热稳定性好.
基于多模式离子推力器离子光学系统束流宽范围引出特性, 设计参数如表3所列.
表 3 离子光学系统设计参数Table 3. Design parameters of the ion optical system for ion thrusters.参数名称 参数指标 栅极材料 Mo 栅极直径/cm 30 束流直径/cm 28.6 屏栅加速栅厚度比 1∶1.1 加速栅减速栅厚度比 1∶1.1 屏栅加速栅孔径比 3∶2 屏栅减速栅孔径比 9∶8 栅间距 1 屏栅加速栅透明度比 1.45 屏栅减速栅透明度比 2.52 空心阴极设计发射电流20 A, 为了满足长寿命需要, 阴极触持极为全石墨材料, 空心阴极发射体为六硼化镧材料. 图6为推力器原理样机照片, 图7为推力器现场工作图.
3. 试验系统和方法
离子推力器性能测试实验是在TS-7 B电推进真空实验系统上开展的, 设备参数参考文献[3], 试验组成如图8所示, 采用激光干涉推力测量系统对推力进行实测, 对比冲进行计算.
束流发散角和推力偏角是束流特性的重要表征量, 其测试方法是: 通过探针采集所在测点的电流密度来确定束流面积和束流中心, 根据束流面积计算束发散角γ, 根据束流中心计算推力偏心角λ, 束流发散角通过测量离子推力器出口下游1000 mm处束流密度分布得到, 束流均匀性分布测量采用单法拉第探针在束流下游通过线扫描的方式进行测量.
束发散角:
γ=2arctgΦT−ΦB2L, (14) 式中, γ为推力器束流发散全角, 一般指束流覆盖90%的区域; ФB为推力器束径(300 mm); ФT为束发散角测试直径; L为法拉第筒到推力器出口的距离.
推力偏角:
λ=tan−1dL0, (15) 其中d为法拉第探测面处束流中心到推力器中心轴线的距离; L0为束流发散角所在直线与推力器中心轴线的交点到法拉第探针探测面的距离.
性能测试试验: 在额定供电、供气参数下进行23个工作点下的测试试验[25], 采用激光干涉推力测量系统进行推力实测和比冲计算.
栅极热稳定性设计验证: 在推力器不引束流的冷态情况下开展栅极间距测试, 测试时, 在栅片上设置有13个间距测量点, 其中1—12测量区域为栅极半径OF (200 mm)处所在直径15 mm的圆心区域, 如上图3所示, 表征栅极四周区域间距的变化, 13测量区域为圆心O所在直径15 mm区域, 表征栅极中心区域的间距变化.
4. 结果和讨论
4.1 功率与束电流、推力、比冲
离子推力器输入功率为工作时所输入的所有电功率之和:
Pin=VbIb+VdId+|Va|Ia+VkIk+VnIn, (16) 式中, Vb, Vd, Va, Vk, Vn分别为屏栅电压、放电电压、加速电压、阴极和中和器触持电压, 单位为V; Ib, Id, Ia, Ik, In分别为束电流、放电电流、加速电流、阴极和中和器触持电流, 单位为A.
离子推力器工作参数与性能的关系如下[23]:
推力:
F=αβ√2MeUbIb, (17) 其中
α 为发散角修正系数;β 为双荷离子修正系数; M是氙离子的质量, 2.18 × 10–25 kg; g为重力加速度, 9.8 m/s2; e是电子的电量, 1.6 × 10–19 C;Ub 为束电压;θ 为90%束流区域内的发散角半角.采用激光干涉推力测量系统对推力进行实测并与(17)式对比修正, α取0.972, β取0.977. 比冲按(18)式计算:
Isp = F˙mg. (18) 从图9和图10可以看出, 多模式离子推力器设计合理, 实现了功率宽范围稳定放电和束流可靠引出, 能够满足多工作点变推力需求, 束电流随功率增大, 从0.3 A不断上升到2.1 A, 推力随功率近线性增加, 在277—3120 W功率下能够实现9.9—117.2 mN变推力宽范围工作, 推力实测值与计算值最大误差小于2.7%, 比冲随功率增加开始快速增大, 台阶上升, 到达一定值后趋于稳定.
功率宽范围调节策略分为2步, 首先是保持屏栅电压一定, 精确调节阳极电流和放电室流率, 在其屏栅电压下引出不同束流, 实现功率小范围调节, 再通过大步长调节屏栅电压实现功率的宽范围调节, 离子推力器工质流量供给分为3路, 即中和器、主阴极和阳极, 一定栅极电压下, 离子推力器的比冲取决于流率, 推力器宽范围调节流率设置原则是: 工作点一定, 主阴极、中和器工作在点状模式且比冲最大化, 即中和器、阴极工作在点状模式的流率最小值为该工作点流率设计值, 阳极流率设置要求是放电室工质利用率不高于放电损耗突增转折点, 且最高不超过90%, 通过电气参数大小步长迭代调节能够实现比冲在1269—3492 s宽范围工作.
4.2 栅极热稳定性设计验证
栅间距和栅间距分布一致性是离子推力器多模式调节和长寿命实现的重要影响因素, 这是由于栅极几何参数不变, 一定束电流下, 受离子光学系统交叉限制和导流系数限制, 栅极可正常聚焦的总加速电压范围一定, 即栅极束流正常聚焦引出的限制其实为束电流和总加速电压之间的限制关系[26]. 如果栅间距不一致, 则交叉限制屏栅电压受限于栅极的最小热态栅间距, 同时导流系数限制屏栅电压又受限于栅极的最大热态栅间距, 栅极可正常聚焦多模式调节范围变小, 当工作在最优导流系数附近, 栅极通过静电场的作用对离子聚焦、定向和加速, 能够避免栅极对离子的拦截. 当工作电压或者是电流超过极限范围后, 由于栅极间的空间电荷饱和效应会造成离子发散角急剧增加, 并导致高能离子对加速栅和减速栅直接轰击, 降低推力器工作寿命.
离子推力器栅极为薄壁多孔弱刚度结构, 由于受到等离子轰击存在热变形, 在工程应用中为了保证栅极组件在多次热冲击下的稳定性, 通常在栅极初始装配后开展栅极稳定化热处理, 方法是将初装的栅极组件装配与离子推力器上按照低功率-高功率, 高功率-低功率顺序开展束流引出, 通过高低温循环释放栅极加工、制造、装配产生的应力, 因此在热稳定性验证实验中分别测量初始装配和热稳定化后栅极间距用以评估热设计.
图11给出了栅极组件采用无应力释放钛安装环和应力释放钼安装环屏栅-加速栅各测量点间距测量值. 优化前, 栅极组件初始装配后屏栅-加速栅二者间距为0.912 mm, 稳定化后释放了应力, 二者间距增大到1.140 mm, 初始装配与稳定化后间距对比, 四周变化百分比25%, 中心变化百分比25.8%. 优化后, 栅极组件初始装配后屏栅-加速栅二者间距为0.892 mm, 稳定化后变化为1.023 mm, 四周前后变化百分比14.6%, 中心变化百分比18.5%.
图12给出了优化前后加速-减速栅各测量点间距测量值, 采用无应力释放钛安装环, 栅极初始装配和稳定化处理后加速-减速栅二者间距分别为0.798 mm和1.004 mm, 采用应力释放钼安装环栅极初始装配后加速-减速栅二者间距为0.804 mm, 稳定化后变为0.825 mm, 四周前后变化百分比从优化前的25.8%降低到2.6%, 中心前后变化百分比从优化前的20.8%降低到7.4%. 表4给出了优化前后栅极稳定化处理后栅间距变化对比数据. 从表4可以看出, 采用应力释放结构优化和钼安装环材料优化热设计提升了栅极热稳定性, 屏栅-加速四周极差从0.1 mm降低到0.01 mm, 中心极差从0.330 mm降低到0.180 mm, 加速-减速四周极差从0.290 mm减低到0.075 mm, 中心极差从0.290 mm降低到0.115 mm, 中心极差较四周大, 主要是因为离子推力器栅片连接采用四周固定, 在栅极组件经历高温工况后远离固定位置的中心发生较大的偏移.
表 4 优化前后栅极稳定化处理后栅间距变化量Table 4. Variation in grid gap with stabilization before and after optimization.四周变化量/mm 中心变化量/mm 四周极差/mm 中心极差/mm 屏栅-加速 优化前 0.228 0.235 0.100 0.330 优化后 0.131 0.165 0.010 0.180 加速-减速 优化前 0.206 0.166 0.290 0.290 优化后 0.021 0.060 0.075 0.115 4.3 束流分布特性
图13给出了推力器不同工作点下的束流发散角、推力偏角测试结果, 可以看出最小功率工作点277 W时, 束流发散角为30.7°, 其余工作点束流发散角均小于30°, 矢量偏角均在1.5°以内, 符合设计预期, 有利于长寿命实现, 且随着功率的增加, 束流发散角呈现变小趋势并趋于稳定, 分析认为该现象与栅极组件引束流的聚焦性有关, 离子鞘层形状相对栅极位置并非固定不变, 而是受等离子体密度, 栅间电压, 栅极材料特性、栅孔孔径和间距等各种因素的影响, 理想工作条件下, 鞘层弯曲成球冠状, 凹面对着屏栅孔, 使得鞘面发射的离子全部进入屏栅孔内, 离子有理想聚焦的运动轨迹, 小功率工作点下等离子体密度低, 鞘层与离子加速区域较远, 部分离子束过聚焦导致束流发散角大[23]. 此外, 高功率下束流发散角较小也与导流系数p有关, 随着功率的增大, 等离子体密度增大, 沉积在栅极组件的能量变大, 栅极有效间距
le 减小促使栅极导流能力增强, 束流发散角变小, 后续采用边缘变孔径设计, 以解决小功率下边缘区域束流欠聚焦问题.离子束流的均匀性是离子推力器设计的关键指标, 是影响离子推力器宽范围调节和长寿命实现的关键因素. 推力器束流密度径向分布测试是准确评估束流平直度的有效方法, 束流密度径向分布曲线是束流均匀性的表征量, 可为推力器性能改善和寿命优化提供数据, 图14给出了推力器功率最大工况点1.2倍拉偏下栅极热稳定性设计优化前后束流密度对比曲线, 可以看出, 改进前离子推力器束流密度分布呈双山峰形, 改进后密度分布曲线基本为平顶形, 束流均匀性较好, 束流平直度提升了10%. 图15给出了优化后7个典型工作点下束流密度径向分布特性测试结果, 可以看出多模式离子推力器随着功率增大, 束电流密度增大, 在不同的工作点下束流密度分布曲线无尖峰, 较平缓, 束流均匀性较好, 束流平直度达到0.75, 说明多模式离子推力器放电室磁场和离子光学系统设计合理且匹配良好, 能够满足电推进系统多模式, 长寿命设计要求.
5. 结 论
1) 基于小天体探测任务对电推进系统的设计需求, 提出四极环形会切场放电室, 凸面三栅变孔径栅极, 全石墨触持极六硼化镧发射体空心阴极、集成式主支撑环承力结构的总体设计方案, 开展了功率宽范围适应性、低功率放电稳定性的放电室参数详细设计;
2) 通过变孔径屏栅设计方案实现了引出束流与放电室电子温度分布相匹配, 束流均匀性好, 各工况下束流平直度大于0.75, 提出应力释放钼安装环热设计并开展试验验证, 结果表明: 栅极稳定化处理是控制栅间距一致性的重要方法, 优化后栅极变形得到有效控制, 栅极间距极差最大下降百分比为90%, 优化前后屏栅-加速, 加速-减速四周极差分别从0.10 mm, 0.290 mm下降到0.01 mm和0.075 mm, 中心极差分别从0.330 mm, 0.290 mm下降到0.180 mm和0.115 mm;
3) 搭建了离子推力器性能测试和束流特性测试平台, 试验结果表明: 多模式离子推力器推力随功率近线性增加, 比冲随功率台阶上升, 束流发散角随功率呈现变小趋势并趋于稳定, 在277—3120 W功率下能够实现推力在9.9—117.2 mN, 比冲在1269—3492 s宽范围稳定工作, 各工作点束流发散角小于30°, 矢量偏角小于1.5°.
[1] Burak K K, Deborah A L 2017 J. Propul. Power 33 264
Google Scholar
[2] Li J X, Wang Z H, Zhang Y B, Fu H M, Liu C R, Krishnaswamy S 2016 J. Propul. Power 32 948
Google Scholar
[3] 李建鹏, 靳伍银, 赵以德 2022 物理学报 71 015202
Li J P, Jin W Y, Zhao Y D 2022 Acta Phys. Sin. 71 015202
[4] Brophy J R, Mareucei M G, Ganapathi C B, Garner C E, Henry M D, Nakazono B, Noon D 2003 Presented at the 39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit Huntsville, USA, July 20–23, 2003 p2003-4542-1
[5] Rayman M D, Varghese P, Lehman D H, Livesay L 2000 Acta Astronaut 47 475
Google Scholar
[6] Garner C E, Rayman M D, Brophy J R, Mikes S C 2011 Presented at the 47th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit San Diego, USA, July 31−August 03, 2011 p2011-5661-1
[7] Malone S P, Soulas G C 2004 Presented at the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit Fort Lauderdale, USA, July 11–14, 2004 p2004-3784-1
[8] Herman D A , Soulas G C , Patterson M J 2007 Presented at the 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit Cincinnati, USA, July 08–11, 2007 p2007-5212-1
[9] Goebel D M, Martinez-Lavin M, Bond T A, King M 2002 Presented at the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Joint Propulsion Conferences Indianapolis, USA, July 7–10, 2002 p2002-4348-1
[10] Snyder J S, Goebel D M, Hofer R R, Polk J E 2012 J. Propul. Power 28 371
Google Scholar
[11] Jahn R G, Von J W 2006 Physics of Electric Propulsion (New York: Dover Pubns) p68
[12] Farnell C C, Williams J D 2011 Plasma Sources Sci. Technol 20 025006
Google Scholar
[13] Bittencourt J A 1980 Fundamentals of Plasma Physics (New York: Springer) p95
[14] Piel A, Brown M 2011 Phys. Today 64 55
[15] Brophy J R, Wilbur P J 1985 AIAA J 23 1731
Google Scholar
[16] Arakawa Y, Wilbur P J 1991 J. Propul. Power 7 125
Google Scholar
[17] Mahalingam S, Menart J A 2002 Presented at the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Joint Propulsion Conferences Indianapolis, USA, July 7–10, 2002 p2002-4262-1
[18] Wang J, Polk J, Brophy J, Katz I 2003 J. Propul. Power 19 1192
Google Scholar
[19] 陈茂林, 夏广庆, 毛根旺 2014 物理学报 63 182901
Chen M L, Xia G Q, Mao G W 2014 Acta Phys. Sin. 63 182901
[20] 龙建飞, 张天平, 李娟, 贾艳辉 2017 物理学报 66 162901
Long J F, Zhang T P, Li J, Jia Y H 2017 Acta Phys. Sin. 66 162901
[21] Wirz R, Goebel D M 2008 Plasma Sources Sci. Technol 17 035010
Google Scholar
[22] Zhang T P, Wang X Y, Jiang H C 2013 Presented at the 33th International Electric Propulsion Conference Washington, USA, October 6–10, 2013 p2013-48-1
[23] 赵以德, 李娟, 吴宗海, 黄永杰, 李建鹏, 张天平 2020 物理学报 69 115203
Zhao Y D, Li J, Wu Z H, Huang Y J, Li J P, Zhang T P 2020 Acta Phys. Sin. 69 115203
[24] 李建鹏, 张天平, 赵以德, 李娟, 郭德洲, 胡竟 2021 推进技术 42 1435
Li J P, Zhang T P, ZhaoY D, Li J, Guo D Z, Hu J 2021 J. Propul. Technol. 42 1435
[25] 李建鹏, 靳伍银, 赵以德 2022 物理学报 71 075203
Li J P, Jin W Y, Zhao Y D 2022 Acta Phys. Sin. 71 075203
[26] Goebel D M, Katz I 2008 Fundamentals of Electric Propulsion: Ion and Hall Thruster (Hoboken: John Wiley and Sons) p115
[27] 张天平, 杨福全, 李娟 2020 离子电推进技术 (上海: 科学出版社) 第91页
Zhang T P, Yang F Q, Li J 2020 Technology of ion electric propulsion (Shanghai: Science Press) p91 (in Chinese)
-
表 1 离子推力器放电室设计参数
Table 1. Design parameters of the ion thruster discharge chamber.
几何构型 直段 + 锥段阳极筒 放电室口径和栅极直径比 1—1.2 长径比 0.4—1 磁极数 4 无量纲后磁极间距离 0.6—0.85 闭合磁等势线/Gs 50—60 磁体体积宽度比/cm2 0.4—0.8 永磁体剩磁/Gs ≥9500 永磁体矫顽力/(kA·m–1) ≥700 表 2 离子光学系统屏栅变孔径分区及归一化后参数
Table 2. Normalized parameters for screen grids with variable aperture zones.
分区 分区范围 孔径 0—OA 0—3 1 OA—OB 3—5 1.05 OB—OC 5—12 1 OC—OD 12—15 0.94 表 3 离子光学系统设计参数
Table 3. Design parameters of the ion optical system for ion thrusters.
参数名称 参数指标 栅极材料 Mo 栅极直径/cm 30 束流直径/cm 28.6 屏栅加速栅厚度比 1∶1.1 加速栅减速栅厚度比 1∶1.1 屏栅加速栅孔径比 3∶2 屏栅减速栅孔径比 9∶8 栅间距 1 屏栅加速栅透明度比 1.45 屏栅减速栅透明度比 2.52 表 4 优化前后栅极稳定化处理后栅间距变化量
Table 4. Variation in grid gap with stabilization before and after optimization.
四周变化量/mm 中心变化量/mm 四周极差/mm 中心极差/mm 屏栅-加速 优化前 0.228 0.235 0.100 0.330 优化后 0.131 0.165 0.010 0.180 加速-减速 优化前 0.206 0.166 0.290 0.290 优化后 0.021 0.060 0.075 0.115 -
[1] Burak K K, Deborah A L 2017 J. Propul. Power 33 264
Google Scholar
[2] Li J X, Wang Z H, Zhang Y B, Fu H M, Liu C R, Krishnaswamy S 2016 J. Propul. Power 32 948
Google Scholar
[3] 李建鹏, 靳伍银, 赵以德 2022 物理学报 71 015202
Li J P, Jin W Y, Zhao Y D 2022 Acta Phys. Sin. 71 015202
[4] Brophy J R, Mareucei M G, Ganapathi C B, Garner C E, Henry M D, Nakazono B, Noon D 2003 Presented at the 39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit Huntsville, USA, July 20–23, 2003 p2003-4542-1
[5] Rayman M D, Varghese P, Lehman D H, Livesay L 2000 Acta Astronaut 47 475
Google Scholar
[6] Garner C E, Rayman M D, Brophy J R, Mikes S C 2011 Presented at the 47th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit San Diego, USA, July 31−August 03, 2011 p2011-5661-1
[7] Malone S P, Soulas G C 2004 Presented at the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit Fort Lauderdale, USA, July 11–14, 2004 p2004-3784-1
[8] Herman D A , Soulas G C , Patterson M J 2007 Presented at the 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit Cincinnati, USA, July 08–11, 2007 p2007-5212-1
[9] Goebel D M, Martinez-Lavin M, Bond T A, King M 2002 Presented at the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Joint Propulsion Conferences Indianapolis, USA, July 7–10, 2002 p2002-4348-1
[10] Snyder J S, Goebel D M, Hofer R R, Polk J E 2012 J. Propul. Power 28 371
Google Scholar
[11] Jahn R G, Von J W 2006 Physics of Electric Propulsion (New York: Dover Pubns) p68
[12] Farnell C C, Williams J D 2011 Plasma Sources Sci. Technol 20 025006
Google Scholar
[13] Bittencourt J A 1980 Fundamentals of Plasma Physics (New York: Springer) p95
[14] Piel A, Brown M 2011 Phys. Today 64 55
[15] Brophy J R, Wilbur P J 1985 AIAA J 23 1731
Google Scholar
[16] Arakawa Y, Wilbur P J 1991 J. Propul. Power 7 125
Google Scholar
[17] Mahalingam S, Menart J A 2002 Presented at the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Joint Propulsion Conferences Indianapolis, USA, July 7–10, 2002 p2002-4262-1
[18] Wang J, Polk J, Brophy J, Katz I 2003 J. Propul. Power 19 1192
Google Scholar
[19] 陈茂林, 夏广庆, 毛根旺 2014 物理学报 63 182901
Chen M L, Xia G Q, Mao G W 2014 Acta Phys. Sin. 63 182901
[20] 龙建飞, 张天平, 李娟, 贾艳辉 2017 物理学报 66 162901
Long J F, Zhang T P, Li J, Jia Y H 2017 Acta Phys. Sin. 66 162901
[21] Wirz R, Goebel D M 2008 Plasma Sources Sci. Technol 17 035010
Google Scholar
[22] Zhang T P, Wang X Y, Jiang H C 2013 Presented at the 33th International Electric Propulsion Conference Washington, USA, October 6–10, 2013 p2013-48-1
[23] 赵以德, 李娟, 吴宗海, 黄永杰, 李建鹏, 张天平 2020 物理学报 69 115203
Zhao Y D, Li J, Wu Z H, Huang Y J, Li J P, Zhang T P 2020 Acta Phys. Sin. 69 115203
[24] 李建鹏, 张天平, 赵以德, 李娟, 郭德洲, 胡竟 2021 推进技术 42 1435
Li J P, Zhang T P, ZhaoY D, Li J, Guo D Z, Hu J 2021 J. Propul. Technol. 42 1435
[25] 李建鹏, 靳伍银, 赵以德 2022 物理学报 71 075203
Li J P, Jin W Y, Zhao Y D 2022 Acta Phys. Sin. 71 075203
[26] Goebel D M, Katz I 2008 Fundamentals of Electric Propulsion: Ion and Hall Thruster (Hoboken: John Wiley and Sons) p115
[27] 张天平, 杨福全, 李娟 2020 离子电推进技术 (上海: 科学出版社) 第91页
Zhang T P, Yang F Q, Li J 2020 Technology of ion electric propulsion (Shanghai: Science Press) p91 (in Chinese)
期刊类型引用(1)
1. 高斌,李娟,陈娟娟,李如月,王彦龙. 离子推力器C-C复合栅多工况热态微位移研究. 真空. 2024(05): 90-96 . 百度学术
其他类型引用(0)
Catalog
Metrics
- Abstract views: 4675
- PDF Downloads: 69
- Cited By: 1