等离子体密度对栅极系统束流影响分析

针对现有离子推力器栅极系统束流引出分析中未考虑放电室等离子体密度均匀性影响的不足,基于数值方法对该问题进行了研究。采用单元内粒子 （Particle In Cell,PIC） 方法对栅极系统束流引出进行了数值模拟。PIC方法中电场求解采用有限差分方法,离子加速方法采用跳蛙格式。以国产20cm Xe离子推力器作为算例,计算得到了栅极中心和边缘孔引出束流离子和电场的分布。结果显示栅极中心孔相对边缘孔束流聚焦性好,束流发散角小,但中心孔轴线的鞍点电势较边缘孔高24.8V。因此,放电室等离子体密度分布对束流引出具有重要影响,电子返流现象最先出现在栅极中心位置。     

背景压强对电推进羽流场影响的数值模拟

真空舱内背景压强是电推力器地面试验过程中影响工作性能评估和羽流场参数诊断的重要参数。针对LIPS-200型离子推力器羽流场参数的数值仿真中采用的背景压强建立方法进行了仿真分析。仿真中采用混合粒子网格(PIC)方法和直接模拟蒙特卡罗(DSMC)方法处理羽流场中等离子体运动和粒子间碰撞,分别采用虚拟粒子和计算粒子建立压强的方式,对电推进羽流场进行了数值模拟,并与绝对真空环境进行对比分析。结果表明:背景压强的存在导致中性粒子和电荷交换离子数密度较绝对真空环境高1个量级以上。虚拟粒子可大幅提高计算效率,获得的流场中电荷交换离子分布与计算粒子结果相近,但中性粒子分布相差较大,虚拟粒子无法表征壁面及真空泵的影响。      

稳态等离子体推进器羽流场数值模拟

采用二维轴对称模型,使用粒子网格法(PIC)和直接模拟蒙特卡洛法(DSMC)相结合的方法,对稳态等离子体推进器(SPT)羽流场进行了数值模拟.采用DSMC方法中的随机取样频率法(RSF)求解粒子碰撞过程,并对比了不同的分配电荷方式、电子运动模型及SPT出口条件时的羽流场.将不同条件计算得到的羽流场中距SPT出口0.1?m,0.5?m及1.0?m处的离子电流密度和电荷密度与实验结果进行了对比,得出在采用面积权重法分配电荷、等熵模型描述电子运动和用实验值设定发动机出口参数时对SPT羽流场数值模拟的电流密度和轴向附近的电荷密度结果与实验结果符合程度较好的结论.     

碘工质射频离子推力器栅极系统束流特性仿真

为分析碘工质射频离子推力器束流特性,应用粒子云网格算法（PIC）,对碘工质射频离子推力器栅极系统进行三维数值仿真。碘工质电离产物中除了大量存在的I⁺,还包含少量I₂⁺、I²⁺、I³⁺几种多价离子。对添加多价离子前后束流特性的仿真结果进行对比,得出该栅极系统离子空间分布、电势分布、离子相空间分布以及束流、束流发散角,并对变密度工况下束流大小进行统计。计算结果表明,程序能较好地模拟离子在栅极系统中的运动情况：添加多价离子后,等离子体悬浮电势有所上升,鞍点电势有所下降,但整体电势分布的变化幅度较小;添加多价离子后束流大小略有增加,束流发散角略有减小,通过理论分析可知仿真所得理论推力及理论比冲均有小幅度增加;放电室等离子体数密度增加到约2.4×10¹⁷m^-3时,该栅极系统达到了束流引出的极限,后续增大等离子体密度引出束流不增反降。模拟结果可以为碘工质射频离子推力器栅极系统的设计提供参考。     

30cm氙离子推力器磁场特性分析与优化设计

放电室磁场设计直接影响放电室的放电稳定性及推力器在轨工作寿命,针对多种工作模式下30cm氙离子推力器磁场设计的复杂性问题,对推力器电磁体磁场向永磁体磁场转换中放电室的磁场特性进行了研究,并对永磁体磁场的关键参数进行优化设计。建立30cm氙离子推力器放电室磁场转换的磁路模型,运用有限元分析理论,利用实际工程数据验证磁路模型计算结果的正确性与方法的可行性。在此基础上,分析获得给定磁路构型下产生要求磁感应强度的永磁体关键尺寸。以放电室工作阳极震荡电压、减速栅极电流、加速栅极电流和磁路系统质量为目标,采用多目标粒子群优化算法,对永磁体的关键参数进行优化,得到30cm氙离子推力器设计性能目标下的磁路构型最优结果。本研究可为高效、稳定工作的离子推力器磁路设计及优化提供方法。     

离子发动机羽流二维轴对称数值模型与验证

为了解离子发动机羽流特别是交换电荷离子(CEX,Charge-Exchange)的分布和流动特性,建立了离子发动机羽流的物理模型,采用粒子网格(PCI,Particle in Cell)方法对2种型号的离子发动机羽流场进行数值模拟计算,与其地面实验数据进行对比分析.结果表明:在CEX离子密度大小及分布、电势的大小及最大电势梯度的位置、CEX离子流动角方面,模拟结果同实验结果符合得相当好.在电势结构方面,由于舱壁电势的影响,模拟结果同实验结果相比有一些差别.羽流模型和计算结果为相关羽流实验和数值模拟研究提供参考.      

脉冲等离子体推力器放电电离二维PIC建模与仿真

摘要： 针对脉冲等离子体推力器（PPT）的放电过程,利用粒子网格 蒙特卡洛（PIC MCC）方法建立仿真计算模型.以LES 6 PPT为例,加入电离碰撞进行电离仿真.通过粒子运动碰撞与电磁场耦合仿真计算得到电流与电路总电阻的变化规律,揭示了PPT放电过程中等离子体密度分布情况.通过对比不加入粒子预分布与加入粒子预分布的两种条件下的计算结果,得到了加入粒子预分布使带电粒子密度计算结果更接近实验结果的结论.根据PPT的工作过程,在放电之前推力器内存在等离子体,所以在仿真研究中应进行粒子的预分布.文中的研究方法对PPT的粒子方法模拟具有一定的参考意义.     

离子推力器羽流热效应仿真分析

离子推力器工作时向外喷出的羽流与航天器表面碰撞,会引起敏感材料热变形等热效应,严重时会导致航天任务失败。针对兰州空间技术物理研究所研制的LIPS-200型离子推力器羽流热效应进行了仿真分析。仿真中,使用粒子网格（PIC）方法处理等离子体运动,使用直接模拟蒙特卡罗（DSMC）方法处理粒子间碰撞,使用Maxwell模型处理粒子与壁面的能量交换,对电推进羽流热效应测量中的部分测点进行了数值模拟。结果表明,仿真结果与实验数据符合较好,离子推力器出口轴线上滞止热流仿真值与实验测量值误差小于17.0%。此外,热流计对流场的影响主要集中在热流计附近0.1 m范围内,对整体流场影响较小。      

多级会切场推力器振荡特性模拟研究

目前多级会切场推力器数值模拟研究采用的主要手段是全粒子模拟方法,但是这种全粒子模拟对计算机的要求较高,因此文章通过建立HEMPT的一维流体模型,再现HEMPT的振荡现象,对放电电流和等离子体参数之间的关联性及等离子体参数分布规律进行研究。沿着两条典型的电子传导路径,利用MATLAB中Simulink模块建立一维流体模型。结果分析表明,推力器内部等离子体参数的振荡周期与放电电流振荡周期一致。推力器的主电离区在两个磁尖端之间,电势降主要集中在出口附近。     

磁拓扑结构对环型离子推力器放电性能的影响

为得到环型离子推力器最佳磁拓扑结构以有效提高推力器放电效率、降低放电损耗,对不同类型磁拓扑结构下的放电通道气体放电过程进行研究。采用PIC MCC数值计算方法对等离子体产生及运输过程进行数值模拟,分析磁场分布对等离子体密度分布、电子损耗率及放电稳定性的影响,结合统计结果,得到推力器放电性能曲线,最后进行试验验证。研究结果表明,相较多极场结构,环尖场磁构型能更好地约束电子运动,大幅降低其在阳极壁面损耗率,增加电子与中性原子碰撞概率,显著提高放电效率、降低放电损耗。     

面向深空探测的30cm离子推力器寿命预估

为评价30cm离子推力器的寿命,提出了有限寿命考核结合栅极仿真模型的推力器寿命预估方法,利用粒子-蒙特卡洛（PIC-MCC）方法建立了栅孔溅射腐蚀模型,开展了30cm离子推力器寿命预估研究,分析了栅孔刻蚀速率、单孔电场及离子引出特性,给出了30cm推力器寿命预估值。结果表明,栅孔直径仿真与实测值一致性较好,误差在20%以内,基于每个寿命小节栅孔实测值对模型的修正是有效的;根据仿真结果,10000h寿命考核后,减速栅最先失效,30cm离子推力器在最大功率（3kW）工况下预估寿命为37540h,能够满足小天体探测任务对电推进系统的长寿命需求。     

同轴微阴极电弧推力器的粒子网格法数值模拟

采用粒子网格方法对同轴微阴极电弧推力器μCAT工作过程进行了模拟研究, 并应用自相似方法对模型进行简化,获得了推力器羽流区的电子数密度分布、离子数密度分布、电势分布及离子轴向平均速度,通过改变磁感应强度和位形分析磁场对推力器内等离子体运动特性及推力器性能的影响。计算结果表明,电子被外加磁场捕获约束在磁力线附近,低速离子与高速电子形成的双极扩散电场加速离子喷出;在相同流量情况下,磁感应强度002T时,离子返流严重,磁感应强度005~030T时,磁感应强度变化对速度影响较小;磁场位形对离子运动和推力器性能有较大影响,磁力线与轴线夹角较小时离子速度下降明显,夹角较大时离子返流严重。     

霍尔推力器分割高偏压电极等离子体放电特性

霍尔推力器通道等离子体与壁面有很强的相互作用,为了降低壁面腐蚀,提高在轨寿命, 针对推力器全通道放电过程建立二维物理模型,采用粒子模拟方法（Particle In Cell）,数值研究了电离区壁面分割高于阳极偏压的低发射石墨电极对推力器放电特性的影响,讨论了放电通道电势、离子数密度、电子温度、电离速率及比冲的变化规律。结果表明：在电离区不同位置分割高偏压电极对等离子体放电特性影响明显,电极位置在电离区前端时,电极偏压高于阳极电压60V时通道内放电等离子体参数几乎不变。而电极位置在电离区末端,电极偏压高于阳极电压18V时就会导致加速区轴向扩张,离子聚焦效果强,电子温度显著升高,电子与壁面相互作用减弱,羽流发散角减小。由此推力器比冲提升约12%,寿命延长,性能提高。     

放电电压和屏栅电压对离子推力器性能的影响

利用试验和数值模拟相结合的方法研究6 cm Kaufman离子推力器放电电压和屏栅电压的变化对其工作性能的影响。试验中,离子推力器使用氩气作为推进剂,测量了多组不同工况下的性能参数。此外,基于Goebel的理论模型模拟了放电电压对束流电流和推进剂利用率的影响;采用单元内粒子 蒙特卡罗碰撞（PIC-MCC方法模拟屏栅电压对束流电流、推进剂利用率和加速栅极电流的影响。试验和数值模拟结果一致,发现当放电电压逐渐增大时,引出的束流电流和推进剂利用率先增加然后趋于稳定;当屏栅电压逐渐增大时,引出的束流电流和推进剂利用率先增加然后趋于稳定,加速栅极电流先减小后趋于稳定。研究可以为提高多模式离子推力器的性能提供参考。