航空发动机排气系统红外辐射特征数值计算研究

基于N—S方程,建立了某型发动机排气系统的数值计算模型;利用辐射传输方程（RTE）积分法所编制的红外辐射特征计算程序,计算了该发动机在地面工作状态下的红外辐射特征分布;在加力状态下,考虑了soot粒子的光谱吸收与发射。分析了水蒸气、二氧化碳、一氧化氮和一氧化碳气体介质的红外光谱吸收与发射,并考虑了红外辐射在大气中的传输,最后给出了在3.0～5．0μm波段内排气系统的红外辐射强度分布。     

Energetic neutral atom (ENA) and charged particle periodicities in Saturn’s magnetosphere

The Magnetospheric Imaging Instrument (MIMI) on the Cassini spacecraft has observed energetic neutral atoms (ENA) and charged particles at Saturn from mid-2004 to the present. The particles often but not always reveal striking periodic behavior that seems to depend on the type of particle and spacecraft location. When subjected to a Lomb periodogram analysis, energetic electrons (>150 keV) exhibited strong frequency peaks near 10.80 h (the nominal or “base” period of Saturn kilometric radiation) during 2006–2008, but essentially no periodicity during 2005. The electron periodograms also show pronounced “double” frequency peaks in 2007 and 2008. Energetic protons (3–26 keV) show strong peaks near the same period for 2005–2007, but none for 2008. Oxygen ions at the same energies display strong peaks for 2005 and 2006, but not for 2007 and 2008. By projecting the ENA images onto Saturn’s equatorial plane or onto a plane perpendicular to the equatorial plane and then summing the data in the appropriate dimension, “strip” images can be constructed from which a time history can be derived. These time histories of ENA emissions are also subjected to a Lomb periodogram analyses. The energetic hydrogen neutrals (20–50 keV) exhibited periodic behavior only during 2007, while energetic oxygen neutrals (64–144 keV) displayed a strong SKR-like period in 2005 and 2006 but not for 2007 or 2008. Some of this behavior may be due to changing spacecraft aspect relative to the ENA emissions, and some of it may be real. This periodic behavior may be consistent with a rotating anomaly that “flashes” brightly in the midnight-to-dawn sector once per 10.8 h, with the flash parameters depending on particle species and energy.     

航天器发动机羽流对敏感器热效应仿真研究

采用差分求解N-S方程与DSMC方法相结合的方法,研究了航天器单台发动机连续工作情况下真空羽流对航天器敏感器的热效应。首先通过求解N-S方程,获得发动机喷管的内流场,然后应用DSMC方法对喷管出口外轴对称羽流场进行计算,最后将轴对称羽流场计算结果作为模拟粒子入口边界条件,在并行计算机平台上进行三维羽流场和热效应计算得到航天器单发动机连续工作情况下羽流场对敏感器的热效应。以两个敏感器为例,对仿真结果进行了分析和比较,并得出了相应结论。     

微波等离子体推力器流动模拟(英文)

采用N-S方程求解了100 W微波等离子体推力器(MPT)选用不同推进工质时的性能参数;并采用直接蒙特卡洛模拟方法(DsMC)对MPT羽流进行了数值模拟.结果表明,几种工质的推力变化不大,氮气为23.6 mN,氮气为24.8mN,氩气为24.8 nuN;但比冲区别较大,氮气为565.2 s,氮气为243.7 8,氢气为180.2 s.羽流场中,密度、压强及温度沿轴向和径向均逐渐减小;轴向速度在轴线附近变化不大,采用氩气工质时,约1 700 m/s,在远离轴线区域,沿流动方向逐渐增大,沿径向逐渐减小;径向速度沿轴向变化不大,沿径向逐渐增大,并在接近流动区域边界时迅速减小.     

霍尔推力器羽流对S波段电磁波的衰减研究

霍尔推力器羽流对S波段电磁波信号有衰减作用,这对星上推力器与天线的布置策略有重要影响。为研究S波段电磁波的衰减程度,建立几何光路算法,对电磁波在等离子体羽流媒质中的传输衰减进行数值描述,并以试验所测的信噪比数据来验证算法精度。验证试验采用空间透射波法,在相同工况的试验与计算结果对比中,最大计算误差为26.1%,平均误差为7.3%;在趋势上,试验结果与计算结果保持一致。在此基础上,针对84个在轨典型工况进行S波段电磁波信噪比衰减计算,结果显示：羽流对电磁波的吸收作用是影响衰减的主要因素,衰减最大区域出现在天线位置0.3～0.6m、方向角15°～30°的区域,衰减值在-7~-4dB之间。该结果可为霍尔推力器与天线在卫星的布置提供参考依据。     

用Langmuir探针诊断MPT羽流

为了探索微波等离子推力器(MPT)羽流的特性,建立了一套Langmuir探针诊断系统并对MPT的羽流进行了诊断,得出了MPT喷管羽流的电子温度和电子数密度等重要特性参数。实验结果表明,所测参数反映了等离子体的特性。与理论计算结果基本相符,并且诊断手段简单易行,为MPT羽流诊断提供了一种实验手段。     

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

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

离子火箭发动机羽流特性分析

分析了离子火箭发动机羽流组成和各种粒子的产生机理,建立了束离子、中性推进剂原子、交换电荷离子的物理模型,并以20cm氙离子火箭发动机为例计算了羽流中各种粒子的空间分布。通过对羽流组分分布特性的分析,提出了羽流污染防护的有效措施和方法。     

对称负重合型气动伺服阀零位流动状态分析

针对负重合型气动伺服阀(PSV)零位特性设计缺少理论依据问题,根据单个节流口的气体质量流量公式建立起气动伺服阀滑阀级的数学模型,采用假设求证法分析了零位时气体流经上、下游节流口的流动状态与影响因素之间的关系。结果表明在不同供、排气压力比下,负重合量不均等系数小于0.5283的气动伺服阀有3种可能的零位流动状态,即上、下游节流口均为亚声速流动,或者均为声速流动,或者上游节流口为声速时下游节流口为亚声速流动;负重合量不均等系数不小于0.5283的气动伺服阀有两种可能的零位流动状态,即上、下游节流口均为亚声速流动,以及上游节流口为亚声速时下游节流口为声速流动。可知气动伺服阀的零位流动状态由负重合量不均等系数与供、排气压力比共同决定。以负重合量不均等系数分别为0.5、1、2的气动伺服阀为例进行了计算和实验验证,实验结果与理论分析相吻合。     

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

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