电推力器在卫星上作为主推力器在国外已经在轨应用^[1-2]，国内已有化学推力器与电推力器混合应用的卫星，全电推进卫星也已在研^[3]。总之电推力器应用于卫星，对卫星在轨工作环境的影响情况仍需研究，这是因为传统的化学推力器燃烧产物在理论上是电中性的气体和水^[4]，它们没有在卫星表面易于吸附的理论基础，但电推力器的羽流产物是等离子体，是电离态的带电粒子，理论上有在卫星表面吸附的倾向。尤其工程实际的离子电推力器正离子束与电子束是分立设置，工作时通过分立的电子束与离子束的汇聚形成宏观电中性的等离子体，汇聚的过程及结果有待研究^[5-6]。离子推力器的离子束是在束流系统的引导下射出的，如果离子束流不能完全聚焦，离子束流将侵蚀离子推力器的钼栅极组件而溅射出钼粒子^[7]，在羽流中携带有钼的中性粒子，这些中性粒子一部分通过离子交换电荷碰撞生成钼离子，一部分通过与电子直接碰撞离化形成钼离子^[8]。在各种粒子碰撞和推力器羽流局部电势的作用下返流到卫星表面形成污染而影响卫星热控性能^[9]。羽流的等离子体不仅可以在卫星表面吸附，它的离子在卫星表面也可能激发出二次溅射，更为有害。文献^[8]中通过数值模拟得到了钼原子的空间分布，文献^[13]中研究了均匀厚度钼膜的光学特性，文献^[7]结合试验并根据经验模型得到了栅极材料的原子溅射率，但是没有实现从钼原子的出射到钼原子的沉积直到对卫星OSR（Optical Solar Reflector）片热控性能影响的完整分析。为了研究离子电推力器羽流中各种成分的相互作用和各种粒子的空间分布情况、电场分布情况，本文以某实践卫星的LIPS-200氙离子推力器应用为研究对象，结合离子推力器在卫星中的布局参数和工作参数建立了离子束与电子束分立设置的工程模型，研究包含电子束与离子束的汇聚并形成宏观电中性的等离子体的过程，得到离子、电子、中性粒子的空间分布和电场分布，以及钼粒子在卫星表面的分布及沉积厚度，预示卫星表面热性能的变化，从而使模型的建立更加贴近工程实际，得到客观的分析结果。

卫星的热性能在地面要通过热真空实验加以确定，电推力器羽流中钼原子沉积对卫星热控性能的影响也可以通过地面热真空实验加以研究确定，但是这种地面热真空实验很难测得准确的结果。原因是，要想实验模拟太空中的真空环境需要至少10-3 Pa的真空度，热真空实验要在密闭的真空罐内进行，空罐内体积有限，一旦电推力器点火喷出羽流，真空罐内的真空度难以维持，而电推羽流打到罐壁上产生的溅射对类似OSR片这种光学表面性能的影响无法回避^[10]，因此，通过地面实验准确确定电推羽流对OSR片光学表面性能的定量影响难于实现。所以，本文拟采用数值建模模拟的方法，辅以实验结果在一定程度上考察模型的正确性，来研究电推羽流对卫星热控性能的影响。

为了建立卫星内部仪器的工作环境，需要对卫星进行热控。一般卫星热控是以被动热控为主，主动热控为辅，即通过调控卫星表面对阳光的吸收率和其红外发射率来建立起卫星内部温度与卫星在轨的空间低温环境的平衡^[11]，这就是被动热控的基本思想。主动热控是在被动热控的基础上，提高或降低局部温度。可见，卫星表面的发射率和吸收率的改变，将是对卫星热控基础的改变，其影响是长期的和单向的，需要给予足够重视。

现有卫星的热控表面主要有两种，一是热控多层材料，二是OSR片。OSR片（光学太阳反射镜）是一种以玻璃为透明薄膜表面层的第二表面镜型热控涂层，其薄膜表面层可用石英玻璃或铈玻璃，金属底层由玻璃背面真空蒸发-沉积得到，通常为银或铝。热控多层主要是绝热作用，OSR片是一种发射率和吸收率可设计的光学表面，是通过发射率和吸收率的设计来改变卫星内部热的流出和外部热的流入，来调节卫星的内部温度。如果热控涂层性能退化，卫星入轨后温度升高，用于卫星温度场数值计算的参数设计初值与参数的实际值之间就会存在一定程度的偏差，以致卫星不能正常工作^[9]。发射率和吸收率是在OSR片制造中即确定了的，卫星热控则通过改变OSR片的铺设面积调节卫星的内部温度。可见，OSR片是卫星内部温度调节的热通道，而电推力器的安装位置一般也恰好在OSR片附近，这也在一定程度上加重了电推力器羽流对卫星热控性能的影响。

随束流飞出的钼粒子，部分沉积在卫星表面及太阳翼表面，成为电推进羽流沉积污染的主要因素。从对电推力器羽流的初步分析可知^[12]，电推力器羽流的钼粒子沉积厚度为几纳米到几百纳米，恰和光学膜系厚度的量值相当，可能会改变光学膜系的性质。因此，需要深入研究热控涂层性能在电推力器羽流钼原子沉积的影响下的退化规律。

由于钼粒子空间分布较准确的实验结果不易得到，这里将电推力器羽流的离子分布模拟结果与实验结果进行比较来检验模型的正确与否。

图 10为离子数密度模拟结果与试验结果对比，可以看出两者相当一致。由于实验在先，没有建模分析的预示，实验中布置测点较少，所以这里也引用文献^[17]的结果给出一个较为细致的比较。

图 10 离子数密度模拟结果与试验结果对比
Fig.10 Comparison between simulation results and experiment results of ion density

在图 11中可见离子浓度在不同径向距离处的参数值都随轴向距离的增加而减小；在靠近推力器中心线位置处，参数值的变化幅度较大，在远离推力器中心线位置处，参数值的变化幅度较小。在径向140 mm曲线和径向180 mm曲线上均出现了明显偏离主束流特征的参数值，说明这些位置的点已经不在主束流区。模拟结果和试验结果的趋势和参数值范围基本一致，在此环节证明了建模分析的有效性。

图 11 离子浓度的轴向分布的模拟结果和实验结果对比
Fig.11 Comparison between simulation and experiment results of axial distributions of ion density

通过计算得到钼离子在电推力器附近卫星表面10 000 h的沉积厚度如图 12所示。可见，由于推力器两侧与推力器出口轴线的夹角相对较小，在推力器两侧会有较大量的沉积，沉积率最大值量级为2.2×10^-7 埃/s，10 000 h沉积量最大值为7.912 埃，而在推力器安装位置的背面（即图3中下部中心位置处）会有较少的沉积。

图 12 OSR表面钼离子10 000 h沉积厚度分布（已取对数，埃）
Fig.12 Deposition thickness of Mo ions on OSR surface in 10 000 h

参阅文献^[11]得到钼膜玻璃数据，见表2。

表2 钼膜玻璃的光学性质
Tab.2 Optical properties of molybdenum-film glass

结合上述分析，得到了基板的光学性质和钼膜的光学性质，复合膜的光学性质为复合模型每层膜的光学性质的组合。经过薄膜的光线会被薄膜反射和吸收，同样也会发生在玻璃-薄膜界面再次发射、透射、吸收。结合典型的OSR参数（吸收率为0.12，发射率为0.78）计算钼的复合膜特性，得到表3。

通过线性插值，得到OSR片复合膜吸收率分布如图 13所示。可以看到，散热板OSR片上部分区域的吸收率有所上升，在推力器附近上升最为明显，这将引起卫星散热板平衡温度的上升。

表3 钼的复合膜特性
Tab.3 Properties of molybdenum composite film

图 13 复合膜吸收率
Fig.13 Absorptance of composite film

结合已经得到的复合膜的光学性质对卫星的热平衡进行分析。OSR表面接收到的外热流有太阳辐射、地球红外辐射、地球反射辐射以及舱体辐射4种热流。忽略地球红外辐射、地球反射辐射以及舱体辐射，认为OSR表面接收到的外热流只有太阳辐射，OSR发出的热流主要是对外的热辐射。卫星南北板散热能力按照q=300 W/m²计算。根据能量守恒定律，可以列出卫星舱板温度计算表达式为

Q₁+Q₂=Q₃

式中：传入的外热流量Q₁=αS，其中α为膜的吸收率，太阳常数S为1 353 W/m₂；内热源Q₂=q；向深冷空间辐射散热量Q₃=εσT4，其中ε为膜的发射率，斯忒藩-波尔兹曼常数σ为5.67×10^-8W/（m²·K⁴），T为平衡温度。

将各热量的计算式代入上式，整理后得电池板能量守恒方程为

代入ε＝0.78，α＝0.12，得到卫星OSR片初始平衡温度为320 K。

由于真实卫星还有辐射换热等热导途径，所以实际结果会与这里的计算有差别，这里相当于绝热情况下的温度，目的是揭示羽流对热控的影响情况。

图 14为离子电推力器工作10 000 h后卫星舱板上钼粒子沉积之后的温度升高情况，可以看到最高温度变化包络为25 ℃，绝大部分区域温度升高不超过10 ℃。由于一般仪器设备的环境温度设计余量为±11 ℃^[11]，可见电推力器羽流将对卫星寿命末期的整星热控设计产生不利影响。因此，对温度变化较为敏感的设备应避免安装在推力器周围，另外通过采用局部加强导热等热控措施可以缓解局部小区域的温度升高。

图 14 OSR片温度变化分布
Fig.14 Distribution of OSR temperature variation

为了研究离子推力器在卫星上的应用可能，对卫星在轨热性能的影响进行了建模分析。建模中尽可能采用了工程化参数，并且针对离子电推力器的离子流和中和电子流分置的实际情况进行了二者的分立式建模。建模中模拟了正负粒子流汇流成羽流的过程，成功获得了离子推力器羽流的工质离子、电子、中性粒子、钼离子的空间分布，电场分布，得到了钼粒子在卫星表面的分布及沉积厚度。比较了模型计算的离子分布与实验获得的离子分布情况，说明了模型分析的正确性，给出了卫星表面热性能的变化及局部区域温升最大可达二十多度的结果。预示了离子推力器卫星在寿命末期应关注羽流产生的对热控的不利影响。由于推力器两侧与推力器出口轴线的夹角相对较小，在推力器两侧会有较大量的沉积。此区域的热平衡温度会有较大升高，最大温升包络可能在二十摄氏度量级，一些对温度变化较为敏感的设备应避免安装在推力器周围。通过采用局部加强导热等热控措施可以缓解局部小区域的温度升高。