基于IFE-PIC方法的Kaufman型离子推进器数值模拟研究.pdfVIP

摘要

随着探索太空的范围愈加广阔深远，Kaufman型离子推进器的设计研发也

朝着更高功率和更大推力的方向发展，将推进器送上太空工作之前需要进行各

方面的功能性验证，相关实验检测耗时较长，试错成本较高。为延长推进器在

深空探测的工作年限，需要重点关注最容易腐蚀失效的栅极器件，而放电室内

部的物理特征对栅极的腐蚀具有重要的影响。为此，采用数值模拟的辅助手段

研究放电室内的微观变化，并采用粒子融合算法降低了三维模型的计算量。具

体研究内容如下：

使用三维各向异性与各向同性IFE算法求解等离子体控制方程，并考虑二

次电子温度对放电室中等离子体扩散与运输的影响，将二次电子温度方程与等

离子体控制方程进行耦合迭代求解，获得放电稳定状态下各种物理量的分布。

并针对LIPS-300推进器在4kW与5kW两种工况下进行了仿真模拟，与实验

结果进行对比验证，得到的束流仿真值均与实验值相接近。

考虑栅极上游等离子体密度的分布对栅极寿命的影响较为关键，为减小栅

极上游等离子体密度分布的方差，重新设计了栅极结构，将圆形栅极改为了环

形栅极。为消除偏置阴极管环形推进器引出羽流所带来的力矩，提高推进器工

作的稳定性，将阴极管对称地布置在内外环磁极连线处。同时，分别对比了有

无二次电子温度梯度情况下，环形通道放电室内等离子体扩散密度分布与扩散

电势分布的差异。

针对放电室内部的阳极壁面容易腐蚀失效的问题，使用浸入式有限元基函

数插值方法获取了离子推进器阳极壁面的离子通量等物理参数，对比了圆台部

分与圆柱部分的离子通量大小，其中在非磁极区域内，圆台部分的离子能量要

高于圆柱部分，而在磁极区域内，圆柱部分的离子能量要高于圆台部分。据此

分析得出，放电室内阳极壁面的腐蚀最容易发生在圆柱部分的磁极位置处，此

区域需要通过特定的结构与材料设计以延长放电室的使用寿命。

关键词：Kaufman型离子推进器；二次电子温度；IFE基函数插值；能量沉积

Keywords:Kaufman-typeionthruster,secondaryelectrontemperature,IFEbasis

functioninterpolation,energydeposition

目录

摘要I

ABSTRACTII

第1章绪论1

1.1课题背景及研究目的1

1.2国内外研究现状2

1.2.1离子推进器研究现状2

1.2.2放电室数值模拟研究3

1.2.3离子推进器放电室内壁面等离子体研究5

1.3主要研究内容6

第2章基于PIC-MCC算法的三维放电室解耦迭代模型的研究8

2.1引言8

2.2放电室三维数值模拟控制方程的导出与基本假设8

2.2.1等离子体控制方程的导出与基本假设8

2.2.2中性原子控制方程的导出与基本假设10

2.3二次电子能量守恒方程的推导与化简12

2.4变权重粒子的PIC-MCC算法14

2.5本章小结18

第3章离子推进器放电室解耦迭代模型的求解19

3.1引言19

3.2放电室内控制方程的求解19

3.2.1三维各向异性有限元算法推导19

3.2.2三维各向同性有限元算法推导21

3.3基于二次电子温度均匀分布下的电场求解22

3.3.1LIPS-300推进器的建模22

3.3.2二次电子温度均匀分布下的求解结果24

3.4基于二次电子温度不均匀分布下的场求解27

3.4.1耦合能量方程的迭代过程27

3.4.2耦合能量方程后的仿真变化29

3.5本章小结33

第4章环形离子推进器放电室的模拟仿真34

4.1引言34

4.2偏置环形离子推进器的模型设置34

对称环形离子推进器的模型设置及其初步计算38

4.4耦合能量方程后对环形离子推进器放电室的模拟仿真39

4.