脉冲等离子体推力器羽流的一体化粒子模拟

将DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法与一维MHD双温放电模型相结合,一体化模拟NASA Glenn PPT羽流。对不同出口偏转角的羽流场进行模拟,给出了不同偏转角度下离子、电势、温度的变化情况,并与实验结果进行了比较。结果显示,该模型具有一体化(从工作过程到羽流)预测脉冲等离子体推力器羽流的能力。人工偏转角的加入,增强了羽流的扩散,使得羽流径向扩散、轴向收缩、回流减弱,同时对放电模型带来的不足有一定改善作用。     

脉冲等离子体推力器羽流回流影响分析

将三维MHD双温入口模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法,模拟实验室PPT样机羽流。验证计算显示,该模型具有模拟脉冲等离子体推力器羽流的能力。对不同初始电压和电容下的羽流场进行了模拟,给出了出口平面返流质量流率的变化情况和20μs羽流中离子和CEX离子的分布情况。计算结果显示,高能量状态对应高质量流率,高质量流率对应高动量,高动量离子和中性粒子对航天器撞击会造成更强影响。     

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

以国内研制的稳态等离子体推力器羽流为研究对象,建立了基于粒子单元法和直接蒙特卡罗模拟法的数值仿真模型,简化了电子动量方程,运用Boltzmann关系式求解了等离子体电势,研究了羽流电子数密度、离子电流密度及流场温度分布。通过与实验数据和已有数值结果对比,以验证模型和编制的程序。研究表明,羽流出口下游0.3~1.0 m处的羽流中心电子数密度约为2.5×1016~5×1015m-3,电流密度约为65~1 A/m2。文中工作对分析等离子体羽流污染具有参考意义。     

脉冲等离子体推力器的设计评估仿真

运用磁流体动力学(Magnetohydrodynamic,MHD)的方法,对脉冲等离子体推力器推力室工作过程开展三维双温MHD数值仿真,并对推力室的工作机理进行了分析研究。通过该模型计算得到的推进剂烧蚀质量和原冲量与实验结果相符,同时对不同初始电压、电容、电极长度及推进剂高宽比下的情况进行了评估,为推力器设计提供指导依据。结果显示,高放电能量、长电极长度及大高宽比对推力器的性能都有一定提升,但这种提升是有限的,在进行推力器设计时,这些因素应进行综合考虑。     

脉冲等离子体推力器研究综述

简要介绍了脉冲等离子体推力器(PPT)的基本工作原理;回顾了国内外脉冲等离子推力器的发展历史;阐述了它的优势与面临的挑战;分析了需要研究的关键技术与发展方向。     

推力器真空羽流热效应计算模型修正及误差分析

开展羽流热效应计算模型精确评估是合理进行热控设计的前提和条件。文章采用数值求解N-S方程和基于工程经验的点源法相结合的方法得到羽流内外流场;而后采用理论经验公式对羽流热效应计算模型中的关键参数(壁面热适应系数)进行分析,并讨论它们的取值范围;最后利用MBB 10 N推力器的地面试验数据,对羽流热效应关键参数的取值范围进行了精确评估,使得羽流热效应计算误差减小到工程允许的范围内,确认了修正后计算模型的正确性和可信度。     

脉冲等离子体推力器羽流三维粒子模拟

以脉冲等离子体推力器(PPT)羽流为研究对象,基于单元粒子法和直接模拟蒙特卡罗法,建立三维物理模型,模拟羽流瞬态流动.假定电子等温流动将其动量方程简化并离散,运用交替隐式法求解等离子体电势.通过与实验数据和轴对称仿真结果对比验证仿真模型,分析流场中电势、浓度分布随时问变化情况.研究表明,本文所建模型能在平行和垂直于电极的两个甲面内有效预测等离子体羽流参数,且仿真结果优于二维模拟.PPT羽流中推进剂返流质量份额极小.研究方法、结果对PPT羽流更深入研究具有参考价值.     

霍尔推力器羽流离子能量实验研究

霍尔推力器羽流中的离子能量分布情况对于评估推力器羽流影响,优化推力器在航天器上的布局具有重要意义。本文使用阻滞势分析器对霍尔推力器羽流的离子能量分布进行了实验研究,获得了推力器在不同工况下羽流场中关注位置的离子能量分布状况。实验结果表明:霍尔推力器羽流离子主要由电荷交换碰撞产生的低能量离子和高能量源离子组成;高能量源离子的分布在推力器轴线达到最大值,低能量离子的分布随着与推力器轴线夹角的减小呈先增后减态势;随着霍尔推力器放电电压的提高,羽流源离子能量分布会相应向高能量方向偏移。     

真空羽流的数值模拟及真空羽流污染的数值模拟分析研究

本文用ＤＳＭＣ方法对真空羽流场进行了数值模拟，所得数值模拟结果符合流动规律。建立了真空羽流污染的基本理论模型，并进行了真空羽流污染的数值计算与分析研究，所得结论与理论分析一致。本文的研究工作为进一步深入开展羽流污染规律研究打下了良好的基础，对工程设计具有重要的参考价值。     

脉冲等离子体推力器内部流场的数值分析

应用所建立的广义拉格朗日乘子形式的磁流体力学模型对平行电极型脉冲等离子体推力器的内部流场进行了数值研究.对等离子体流动过程的分析表明,放电初期等离子体在很强的电磁力作用下高速喷出推力器,而到放电后期等离子体受到的气动力比电磁力大得多,等离子体的运动非常缓慢.电磁冲量占推力器元冲量的大部分,增强电磁加速作用可以进一步提高推力器的性能.     

用于微小卫星推进装置的脉冲等离子体推力器

论述了脉冲等离子体推力器(PPT)的工作原理和研究状况。在分析其主要优缺点的基础上，讨论了PPT技术的发展方向，并就PPT的应用前景作了展望。分析认为PPT是一种先进的小推力动力装置，特别适合于微小卫星的轨道转移、姿态控制、位置保持以及星座编队飞行。     

霍尔推力器等离子体羽流粒子模拟

建立了霍尔推力器羽流仿真模型,用单元粒子-直接模拟蒙特卡罗（PIC-DSMC）混合方法对SPT-7推力器的流场进行数值模拟,分析背压、扩张角和电子温度对流场的影响。结果表明：背压粒子增加了回流区内离子和高速粒子,加重羽流污染。SPT-70推力器羽流出口处扩张角约为30°。实验数据验证了仿真模型的正确和方法的可行,对电推力器及其羽流污染等的研究有一定的参考价值。     

离子推力器羽流特性及其污染分析

介绍了与电推进系统有关的空间环境效应的形成原因及其对航天器性能、寿命等的影响。阐述了离子火箭发动机羽流内束离子、中性推进剂原子、交换电荷（CEX）离子和电子等主要成分与航天器相互作用的过程及机理。分析表明，离子推力器出口处的中性推进剂原子与高速束离子流碰撞后产生的CEX离子Xe^＋，以及带电离子轰击推力器组件特别是加速极所产生的金属CEX离子，是造成离子火箭发动机羽流污染的主要成分。在此基础上提出了若干防污染措施。     

轨控发动机羽流红外辐射的数值仿真

针对工作在高真空环境下的轨控发动机,数值模拟了其羽流红外特性。首先计算了考虑化学反应的轨控发动机喷管的内流场和外流场,得到了温度、压力、组分浓度等参数的分布。基于HITRAN 2008和HITEMP 2010数据库采用逐线积分法编程计算了气体光谱吸收系数。最后在此基础上用有限体积法建立的求解辐射传输方程的模型,计算得到羽流红外辐射强度在2~10μm范围内随波长变化的曲线。分析了羽流气体组分、波长、探测角度对羽流光谱辐射的影响,与同类文献中的计算结果进行了比较,结果表明：本文的计算模型和方法能较好地模拟轨控发动机羽流的红外辐射特性。     

氢氧发动机模型真空羽流场试验和仿真研究

研制了一个用于模拟中国长征火箭二级的60 N推力氢氧发动机的缩比模型,并在北京航空航天大学真空羽流效应实验系统进行了试验。使用皮托管阵列测量了羽流压力场,结果显示当距发动机喷管出口的距离从140 mm增加到600 mm时,羽流场的最大压力从12 400 Pa降到了400 Pa。为验证CFD-DSMC混合的数值仿真方法,将试验结果与仿真结果进行了对比分析,二者一致性非常好。对比结果显示数值仿真方法在羽流效应分析方面的强大功能。研究获得了模型发动机羽流场的压力分布特性,可用于原型发动机的羽流效应分析。     

电弧放电等离子体对超声速边界层影响的数值模拟

基于电弧放电等离子体热阻塞机理,对等离子体超声速流动控制过程进行了数值模拟,研究了等离子体对边界层的影响,分析了放电区大小、温度等对其作用效果的影响。结果显示：在高温等离子体放电区的上下游近区发生了边界层分离及漩涡运动;在放电区内有两种边界层,即高温等离子体与外界低温气流之间的温度边界层和气流与壁面之间的粘性边界层;放电区内形成漩涡运动的原因有两个,即内外压差和边界层分离;上游边界层的分离点y轴坐标随温度的增大而减小;增大来流速度,放电区上游分离点y轴坐标呈先增大后减小的趋势、放电区内漩涡运动加剧、下游近区边界层分离点y轴坐标减小。     

脉冲等离子体源控制航天器表面充电电位的研究

在空间环境运行的航天器存在表面充电现象,而航天器表面充电引发的静电放电是导致航天器异常及故障的重要原因之一。因此,在航天器设计和应用中,必须对航天器表面电位采取必要的控制和防护措施。文章介绍了用脉冲等离子体源进行航天器表面充电电位主动控制的研究。通过模拟实验和实验数据分析,证实了用脉冲等离子体源能有效地控制航天器表面充...     

等离子体电弧激励器特性及减阻的仿真研究

设计了一种用于高超声速流动控制的等离子体直流电弧激励器,运用多物理场有限元方法对激励器的工作特性进行了数值模拟,得到了激励器准确热源分布和不同工况下电势、温度等的分布情况。其中,表面电弧激励器在电流7 A、阴阳极距离3.5 mm时在阴极表面最高能得到2.690 5×10~(10) W/m~3的热源值。运用有限体积方法对激励器的斜坡减阻效果进行了数值模拟,得到马赫数为7时总阻力系数最多减小32%的结论。