翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

低速增压风洞支架干扰的数值研究

对气动院FL-9低速增压风洞尾撑支架干扰进行数值计算.通过合理的拓扑结构,生成了YF-16模型的分区对接网格.计算对模型的典型安装状态进行研究,获得了全机气动特性和支架干扰量,探讨了此安装状态下弯刀对飞机气动力的干扰特性,并重点研究了Re数对干扰量的影响规律.     

螺旋桨瞬态流场数值模拟

建立完善的螺旋桨绕流理论异常困难,这一方面是由于螺旋桨流场整体旋转性、非定常非线性、三维效应、气流粘性、涡桨干扰等原因;另一方面是螺旋桨流场的流场结构有待于进一步认知。文章采用CFD方法对螺旋桨的瞬态流场结构进行了重点研究,详细绘制了自启动瞬间直至充分发展的螺旋桨瞬态流场结构图及所受气动力的变化图,分析了流场变化特性及其对升力的影响。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

“嫦娥三号”探测器月尘环境分析及试验要求

"嫦娥三号"是我国首个成功实施地外天体表面就位探测的航天器,在其研制过程中充分考虑了月面环境对探测器的影响,开展了针对性设计及地面验证工作。文章综述了探测器针对月尘环境方面的设计因素,总结了月面月尘激扬的因素及其相关数值分析工作,对月尘验证试验的设计及条件制定进行了说明。相关结果或工作思路可为后续月面探测器及火星表面探测器的研制提供参考。     

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

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

压气机稳定边界的数值预测研究

为研究不同压气机稳定性预测方法的准确性,采用最大静压升系数法与商用CFD软件,对某4级低速压气机的稳定边界进行了数值模拟,并与试验结果对比。计算结果表明:最大静压升系数法在预测压气机稳定边界上有一定的准确性。使用商用CFD软件NUMECA对某4级压气机进行模拟计算,考虑了不同的网格划分方案、使用不同的湍流模型和差分格式等因素,以出现数值发散作为压气机失稳的标准,结果表明:主流区的网格对计算结果影响较大。不同差分格式的结果显示,中心差分格式的结果最好。     

Numerical simulation of the three-dimensional flow-path in a ram-rotor

A kind of three-dimensional flow-path of the ram-rotor was designed in this paper according to the design method of the typical supersonic intake of aircraft and missiles,especially the design method of the compression ramp,the throat and the subsonic diffuser.The design Mach number of the flow-path was based on the averaged relative Mach number,regardless of the change of incoming Mach number along radial direction of the flow-path.Then,the flowfield of the designed flow-path was simulated numerically with computational fluid dynamics software at design point.The purpose of simulation was to obtain the details of the flowfield,including the structure of the shock wave system,position of flow separation zone,and so on.The performance of the ram-rotor was also evaluated.The numerical results show the structure of the shock wave system in the flow-path is similar to that in the supersonic intake.To improve the overall performance of the ram-rotor,the reasonable compromise between the total pressure ratio and the isentropic adiabatic efficiency should be reached.The structure optimization of the flow-path and the configuration of the shock wave system have great influence on the performance of the ram-rotor.      

适用于飞机结冰适航取证加速的分析验证/取证思路

飞机结冰问题研究目前主要有3大手段:数值模拟、风洞试验和飞行试验。这3大手段有着各自的局限性,其中,又以飞行试验最为耗费时间和资金。文中总结了机结冰相关研究方法现状,分析了各种方法的优势与局限性,介绍了模拟分析方法(Certification/Qualification by Analysis,CQbA)方法,给出了可用于CQbA相关工作的软件和设备简介。CQbA方法主要以数值模拟和风洞试验相结合的形式进行研究,强调各方法自身优点的协同关系,进行飞机结冰问题的研究和适航相关结果的认证,以实现高效、经济的研究过程。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。