一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

前缘机动襟翼气动力设计计算研究

本文应用欧拉方程对翼身组合体加装前缘机动襟翼构型进行数值计算，通过跨音速气动特性的计算和分析，得出前缘机动襟翼对翼身组合体气动特性的影响规律，为某型机加装机动襟翼的改进设计提供数值计算依据。     

强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

轴流压气机转静二维非定常薄层N-S数值模拟

采用非定常薄层Navier-Stokes紊流运动方程来描述一跨声轴流压气机中径处的二维转-静非定常流场，速度脉动的相关分析表明转子尾迹对静子的热流流场几乎没有什么影响，而粘性的影响取决于叶片排之间的轴向距离。     

三维极坐标视线动力学

阐述了极坐标视线动力学方程的研究现状。建立了三维极坐标视线动力学方程，讨论了该动力学方程的简化问题，得到了特殊平面（攻击平面）内的简化形式，该简化形式在特定的条件下与三维极坐标视线动力学方程是等价的。     

涡轮叶栅冷气掺混计算的数学模型与方法

讨论了目前涡轮叶栅冷气掺混的几种主要计算方法,提出采用多组分守恒方法求解油气比,从而进行变比热容和变气体常数的冷气掺混数值计算。     

跟踪问题最优控制律精细积分

构成有限时间最优跟踪系统的控制律需要求解 Riccati微分方程及外部控制输入向量满足的微分方程 ,前者是非线性矩阵微分方程 ,后者是变系数线性微分方程。在结构力学与最优控制的模拟理论基础上所发展的精细积分方法借鉴了计算结构力学中的算法 ,可以精确有效地求解这些微分方程。这种方法的特点之一在于步长幅度变化较大时 ,Riccati微分方程的数值解仍可以保持很高的精度 ,并且变系数线性微分方程的求解亦可纳入其体系而不必用通常的差分方法。本文介绍了用精细积分方法求解这些方程的过程 ,并给出了数值算例。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

马赫数对侧压式高超音速进气道及等直隔离段三维内流场的影响的数值分析

用 N-S方程和 RNG k-ε紊流模型计算了 RBCC用侧压式高超音速进气道三维内流场 ,重点分析了马赫数对流场的影响。从出口截面上气流参数的均匀程度来看 ,Ma3 时比 Ma6时的更均匀 ;从压缩气流的主要部位来看 ,Ma6情况下在支板前肩点之后到喉部截面之前的这个等宽度收敛通道内 ,而 Ma3 时主要是在前肩点之前的收敛通道内 ;从对气流的压缩程度来看 ,Ma6比 Ma3 更大 ;从唇口激波的强度来看 ,Ma6比 Ma3 更强 ;从进气道的流量捕获率来看 ,Ma6比 Ma3 更高 ;高超音速进气道内靠近侧壁以及侧壁和顶板的角区附近是最有可能出现亚音速的区域。