固体火箭冲压发动机含硼推进剂燃烧性能实验研究

利用地面直连试验系统对含硼推进剂在某固体火箭冲压发动燃烧过程性能进行试验研究。通过对含硼推进剂燃烧后凝聚相样品的SEM,EDS和XPS分析,探讨推进剂燃烧过程。在凝聚相EDS分析中,硼元素含量随着远离推进剂轴向方向显著降低,氧元素含量显著升高。在补燃室中,由于补充富氧空气,一次燃烧产物进一步反应,导致环境温度上升。由于高温,硼颗粒发生燃烧,产生大量气态硼化物,从而导致硼元素含量下降。二次燃烧凝聚相产物中,硼的非完全氧化物占比在40%以上,氮化硼占比在20%以上,硼颗粒占比7%以下。研究结果表明,随着推进剂在燃气发生器和补燃室内的一、二次燃烧,硼颗粒逐渐减少,并分别与环境中的C和O等元素发生化学反应,在凝聚相中的含量逐渐降低,氧元素在补燃室之后显著增加,氮化硼为凝聚相主要成分之一,存在于各特征位置。推进剂中的硼颗粒没有被完全燃烧,燃烧效率有待于进一步提高。     

基于惯性空间稳定的双轴惯导旋转调制方法

陀螺标度因数误差是影响长航时船用旋转调制惯导系统的关键误差源,其与地球自转和载体运动的耦合误差,可导致惯导系统误差发散。针对此问题,结合船用惯导使用特点,采用外航向、内俯仰的双轴旋转框架结构。在此基础上,提出了一种基于惯性系的双轴旋转惯导系统多位置转停调制方案,通过补偿地球自转和载体运动在双轴旋转惯导内外框架旋转轴上的投影分量,可显著降低陀螺标度因数误差对长航时导航精度的影响。数学仿真和船载试验结果表明,在载体航向角运动的场景下,该方法与传统的双轴旋转调制方案相比可有效抑制地球周期项振幅的增大,系统导航位置误差的发散也降低50%以上。     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

低压工况下10N姿控推力器燃烧室工作过程仿真

描述了发动机燃烧室工作过程仿真模型的主要参数,应用此模型对10N推力器在不同供应管路压强条件下燃烧室的脉冲工作过程进行了仿真.对比分析了长脉冲和短脉冲两种工作模式下管路压强逐渐降低时燃烧室压强、推力等参数的变化规律.通过长脉冲工况的仿真理论推力与实验测量推力反推出管路压强为1.48,1.0,0.6MPa时推力效率分别为94%,91%和90%.短脉冲工况仿真结果表明:管路压强为0.4MPa时推力器已经无法正常点火,应避免推力器工作于此工况.d:\PDF\.pdf      

高超声速气动力及激波位置快速计算方法研究

采用快速计算方法进行高超声速气动力计算时,影响计算精度的关键问题主要在于模型面网格的划分和计算方法的选取。采用一种灵活实用的结构化面网格划分策略,使得模型的各个部件能分别选择合适的计算方法;发展一种基于近似流线的二阶激波膨胀波方法,该方法可以用于多种具有三维流场特性的部件,不仅降低对使用者的经验依赖,还能提高计算精度;配合激波位置计算方法,可以较为准确地计算模型的激波位置,保证边界层外缘参数的计算精度;粘性力计算使用基于起始面元修正的Spalding-Chi方法和参考温度方法。通过对四个典型算例的计算与分析,表明本文发展的高超声速气动力计算方法具有较高的计算精度,能够作为高超声速飞行器初步设计阶段的气动力快速分析工具。     

基于真空羽流试验的洁净真空系统设计

研制的洁净真空系统可用于真空羽流试验及热真空试验研究.详细介绍了真空系统的组成、方案设计及工作模式;计算了不同工作模式下的真空抽气时间;分别对罗茨泵机组、分子泵及低温泵进行了容器极限抽速测试,测试结果表明,满足技术指标;最后对真空容器压力维持能力进行了测试,测得容器平均压升率为0.111Pa/h,对应的平均漏气率为7.24Pa·L/s,真空容器达到了较高的工艺水平.真空系统有多种工作模式供选择,可根据试验需求及故障模式选取,增加了“真空羽流效应实验系统”运行的实用性和经济性.      

辐射板强化喷管换热与红外抑制效果试验研究

试验研究了在发动机喷管中加装金属辐射板前后,喷管壁面温度、热喷流温度与喷管红外辐射特征的变化.结果表明,加装金属辐射板后,热喷流与喷管壁面之间的热量传递显著增强,热喷流中心温度降低,壁面温度明显升高,在90°方向上,热喷流3 ～ 5μm波段的红外辐射强度降低了38.5％.文中从热喷流、喷管壁面以及金属辐射板等相关部件的温度变化情况对红外辐射强度的变化原因进行了解释.     

双模态超燃冲压发动机流量匹配的临界面积法

建立基于控制体法的0-D守恒方程(流量连续方程、动量守恒方程和能量守恒方程)的热力喉道(Ma=1.0)求解集总参数模型,该模型包括了壁面摩擦模型、燃油喷射模型、燃油雾化掺混模型和化学动力学模型.并针对集总参数模型求解热力喉道的特殊性,提出了基于流量平衡的临界面积法,应用于隔离段和燃烧室的1-D流场计算,实现了双模态超燃冲压发动机具有热力喉道时的各种模态的隔离段和燃烧室的流量平衡计算,精确捕捉到热力学喉道,确定了隔离段流动状态和燃烧室的工作模态及相关的流路沿程参数.计算结果表明:采用了临界面积法的集总参数计算模型能解决热力学喉道求解问题,其计算精度达10-4,单点工况计算时间小于0.1s.      

静止可燃气热射流爆震起爆块结构自适应网格加密数值模拟

为了深入认识理解静止可燃气热射流爆震起爆的机理过程,采用块结构自适应网格加密方法的开源程序AMROC,进行二次开发,在LINUX系统中建立小型集群并行计算系统,应用于静止可燃气热射流爆震起爆过程高效精细数值模拟,并根据流场特性划分三个阶段进行分析。结果表明:对于一定条件下的热射流,静止可燃气中的热射流爆震起爆对壁面的依赖性很强,正是通过壁面的反射作用形成局部爆震燃烧;通过左壁面和上壁面的有效二次壁面反射,上壁面处的局部爆震燃烧进一步增强,强化的局部爆震燃烧逐渐向整个等直管道中传播;上壁面和下壁面的局部爆震燃烧同时形成之后,双三波点结构在管壁中来回碰撞碰,形成周期性的稳定爆震燃烧。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.