应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

针栓式喷注器锥形液膜破碎特性试验

采用高速摄影获得了针栓式喷注器在不同喷注压降和结构参数下的表面波破碎图像,测量了锥形液膜的破碎长度和破碎时间,研究了变工况时液膜破碎长度和破碎时间的变化规律.试验结果表明:在喷注压降不变的条件下,针栓式喷注器能够实现流量的线性调节.针栓式喷注器设计时,在合理的推进剂动量比范围内,狭缝宽度应尽量取小.液膜在低工况时破碎得更快.液膜破碎长度和破碎时间均随喷注压降的增加而减小.      

固体推进剂二维绝热微尺度燃烧模型的数值研究

为实现固体微推力器工作过程的一体化模拟,基于Fluent计算软件的二次开发功能(UDF)和简化化学动力学模型,实现了固体推进剂的二维气-凝相绝热微尺度燃烧模型的建立,该模型针对固体微推力器所用双基推进剂,包含两步凝相反应和五步气相反应,燃速、推进剂表面温度和组分质量分数基于燃面物理特性计算得到,并考虑了粘性作用对气相和凝相反应的影响。针对0.5MPa,1.0MPa,2.0MPa和5.1MPa四种工况进行了计算,结果表明,高压工作环境下出现发光火焰区,且随表面压力增大而逐渐靠近壁面,凝相反应区厚度和嘶嘶区、暗区主要反应物在燃面的质量分数随推进剂表面压力增大而减小。对称面处推进剂燃速,推进剂表面温度和气相火焰结构与实验结果基本一致。由于壁面附近较高的粘性作用,气相火焰在壁面位置更加靠近推进剂燃面,并导致壁面位置推进剂燃速高于对称面位置。该模型实现了二维环境下考虑分步凝相反应的推进剂绝热燃烧模型的一体化计算,较好地拓展了原模型的应用范围。     

基于模型的民机液压告警系统验证方法研究

利用Simulink/Labview建立民用客机液压系统告警逻辑模型,通过注入顶层定义的根源故障获得可视化告警信息指示,符合设计分析结果,且具有故障叠加告警显示等优点,为验证复杂系统集成中面临的根源派生故障告警信息显示抑制和排序功能提供一种有意义的方法。     

航空燃油类型对催化惰化系统性能的影响

在设计了一种催化惰化系统流程并描述其工作原理的基础上,以从油箱中抽吸气体的摩尔流量为基准,推导了流经催化反应器后各气体组分的流量关系,通过质量守恒方程及气体平衡溶解关系,建立了油箱气相空间气体浓度变化的数学模型。选择了RP-3、RP-5和RP-6燃油作为对象,用所建立的数学模型计算了不同载油率和催化反应器效率下的气相空间氧浓度变化关系。研究显示,由于3种燃油的蒸汽压不同,造成从外界环境补气及进入油箱的混合惰气流量不同,从而导致气相空间氧浓度的变化规律差异远大于采用中空纤维膜的机载惰化系统。因此,在设计催化惰化系统时要充分考虑燃油类型对惰化系统性能的影响。     

基于全局灵敏度分析的侧向气动导数不确定性对侧向飞行载荷的影响

考虑飞行载荷计算中使用的气动导数存在不确定性,利用基于方差的全局灵敏度分析(GSA)方法,结合偏航机动的动力学模型,分析了侧向气动导数不确定性对侧向飞行载荷的影响。以某型飞机为例,运用该方法得到侧向气动导数的全局灵敏度排序。结果表明:侧滑角、阻尼贡献的垂尾侧向载荷及垂尾侧向总载荷受全机侧力系数对侧滑角导数的影响最大,受航向静稳定导数及方向舵操纵效能的影响次之;方向舵偏度贡献的垂尾侧向载荷只受全机侧力系数对方向舵偏度导数的影响;无尾飞机侧向载荷主要受航向静稳定导数、方向舵操纵效能及无尾飞机侧力系数对侧滑角导数的影响;偏航阻尼导数基本不影响各侧向飞行载荷。同时也验证了方法的有效性,对提高飞行载荷的计算精度有一定的指导意义。     

Powered-descent trajectory optimization scheme for Mars landing

This paper presents a trajectory optimization scheme for powered-descent phase of Mars landing with considerations of disturbance. Firstly, θ–D$\theta ''–''D$ method is applied to design a suboptimal control law with descent model in the absence of disturbance. Secondly, disturbance is estimated by disturbance observer, and the disturbance estimation is as feedforward compensation. Then, semi-global stability analysis of the composite controller consisting of the nonlinear suboptimal controller and the disturbance feedforward compensation is proposed. Finally, to verify the effectiveness of proposed control scheme, an application including relevant simulations on a Mars landing mission is demonstrated.     

基于知识推理的航天器自主故障诊断方法

面向深空探测这类型号任务,提出一种基于知识推理的航天器自主故障诊断的智能方法,采取直接线性存储的知识库结构、主动监测事实更新并触发推理以及基于散列算法的事实管理三项关键技术,适用于航天器硬件处理能力和存储容量受限的应用环境.采用C语言在基于TSC695处理器的嵌入式系统上实现了该故障诊断方法,实验的结果表明该方法有效.     

航天器嵌入式操作系统研究与设计

随着航天器结构和功能的复杂化,对作为系统资源管理者的操作系统的功能和可靠性等要求也日益增强.分析嵌入式实时操作系统在国内外航天领域的研究与应用现状,提出面向航天器应用特点的具有多任务管理、内存管理、文件系统、容错和故障管理等功能的高可靠、高性能的嵌入式实时操作系统的研究与设计方案,并展望其在未来网络化及分布式安全航天器中的发展前景.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.