Radical recombination in a hydrocarbon-fueled scramjet nozzle

To reveal the radical recombination process in the scramjet nozzle flow and study the effects of various factors of the recombination, weighted essentially non-oscillatory(WENO)schemes are applied to solve the decoupled two-dimensional Euler equations with chemical reactions to simulate the hydrocarbon-fueled scramjet nozzle flow. The accuracy of the numerical method is verified with the measurements obtained by a shock tunnel experiment. The overall model length is nearly 0.5 m, with inlet static temperatures ranging from 2000 K to 3000 K, inlet static pressures ranging from 75 k Pa to 175 k Pa, and inlet Mach numbers of 2.0 ± 0.4 are involved.The fraction Damkohler number is defined as functions of static temperature and pressure to analyze the radical recombination progresses. Preliminary results indicate that the energy releasing process depends on different chemical reaction processes and species group contributions. In hydrocarbon-fueled scramjet nozzle flow, reactions with H have the greatest contribution during the chemical equilibrium shift. The contrast and analysis of the simulation results show that the radical recombination processes influenced by inflow conditions and nozzle scales are consistent with Damkohler numbers and potential dissociation energy release. The increase of inlet static temperature improves both of them, thus making the chemical non-equilibrium effects on the nozzle performance more significant. While the increase of inlet static pressure improves the former one and reduces the latter, it exerts little influence on the chemical non-equilibrium effects.     

非线性因素耦合下电动作动器的 动力学建模与仿真研究

针对飞行器半实物仿真领域对电动作动器高精度建模的需求,对电动作动器的各分系统分别进行动力学方程描述,建立了含间隙、刚度和摩擦等非线性因素的完整电动作动器模型。通过对比模型仿真与实测结果,表明所建分系统的动力学方程是合理正确的。建立的模型对于指导作动器的设计和优化均具有较大的参考价值。     

PR方程推算氮气/灭火剂蒸气等熵指数

氮气(N2)/灭火剂蒸气组成的混合气体的容积等熵指数和温度等熵指数随充填压力和温度的变化关系是影响机载灭火系统灭火剂释放过程的重要因素。基于PR(Peng Robinson)方程和范德瓦尔混合规则,编制了N2/灭火剂混合气体的容积等熵指数和温度等熵指数的计算程序,计算了三组混合气体N2/HFC227ea、N2/CF3I和N2/CF3Br充填压力分别为42MPa和25MPa,初始温度为293K时等熵指数随温度的变化曲线。结果表明:它们的容积等熵指数和温度等熵指数随着温度的升高而逐渐下降,呈近似线性关系。在相同温度下,充填压力为42MPa时的三组混合物的容积等熵指数和温度等熵指数均高于充填压力25MPa时的结果。在相同的充填压力和温度下,N2/CF3I的等熵指数最大,其次为N2/HFC227ea和N2/CF3Br的最小。      

喷射成形7055铝合金热变形行为模拟

为研究喷射成形7055铝合金的热变形行为,在应变速率为0．001～5 s －1、变形温度为300～450℃、工程应变量为50％条件下,在 Gleeble-3500热-力模拟试验机上进行热压缩实验。结果表明：喷射成形7055铝合金的流变应力随应变速率的增大而增大,随变形温度升高而减小。在应变速率为5s －1时由变形热引起的温升达25℃,经修正流变应力比实测值增高20 MPa。采用包含 Z 参数的 Arrhenius 双曲线正弦本构方程可准确描述喷射成形7055铝合金的热变形流变应力行为,变形激活能为146．91 kJ·mol －1。所建本构方程的平均相对误差（Er ）为2．89％,说明可准确预测喷射成形7055铝合金的热变形流变应力。     

氩气雾化FGH95合金的热模拟实验

采用Gleeble-1500D热模拟试验机进行氩气雾化FGH95合金的热压缩实验,在不同的温度和应变速率下,获得FGH95合金的变形应力应变曲线,根据变形数据,建立FGH95合金的变形本构方程,并基于动态材料模型,绘制合金的热加工图。计算得到氩气雾化FGH95合金的变形激活能Q=695.78 kJ/mol,通过建立的本构方程计算得出的峰值应力与实验值符合较好,平均误差范围约6%;根据热加工图,确定FGH95合金安全的热加工区域如下：1070~1100℃,0.01~0.001 s-1,当温度增加到1100℃以上后,应变速率可以增大到0.5 s-1。     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。     

Kepler方程的一种迭代加速算法

采用Steffensen迭代法求解椭圆轨道下的Kepler方程,克服了不动点迭代法收敛速度慢以及Newton迭代法需要重复计算导数值和迭代结构相对复杂等方面的不足。试验表明利用该方法求解Kepler方程既满足了高精度要求,而且收敛速度快。同时它还能满足不同的偏心率和不同的精度的要求,因而是求解Kepler方程的一种较理想算法。     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

双微射流作动器定向控制主流数值模拟

采用改进的Beam-Warming近似因式分解法求解二维、粘性、非定常、可压全N-S方程,对双微射流作动器定向控制主流的合成流场进行数值模拟。计算结果揭示了主流在双微射流作用下会比在单微射流作用下产生更大的偏转,并探讨了其产生偏转的原因。分析了主流的偏转程度与相邻作动器相位差之间的关系,即相邻作动器之间的相位差必须在一定的范围内,主流才能发生更大的偏转。     

多守恒性的差分格式的分析和比较

球面大气浅水波方程(在忽略摩擦和外源强迫的条件下)具有总能量、总质量、总位涡、总位涡拟能和总角动量守恒等五个重要的物理守恒性.然而,将大气浅水波方程差分离散化后,常常不能保持这些守恒性.本文对新发展的两种较好的差分格式进行比较;一种可称之为修正总能量守恒格式,它可以较好地保持上述五个守恒性中的四个;而另一种可称之为准辛格式,它可以使上述五种守恒性均近似保持.文中对这两类格式作了具体分析,并作相应的数值试验和比较,结果表明这两类格式均值得推广应用.