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

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

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

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

GPR格式推广用于多维复杂流动的数值计算

采用算子分裂技术，将原始求解一维流动问题的高精度、高分辨率ＧＲＰ格式推广用于包含各种不同间断解及其相互作用的多维复杂流动计算。以轴对称管端出口处激波绕射传播为例，作了数值模拟，得到了与实验结果吻合很好的流场精细结构。     

GRP格式推广用于多维复杂流动的数值计算

采用算子分裂技术，将原始求解一线流动问题的高精度、高分辨率GRP格式推广用于包含各种不同间断解及其相互作用的多维复杂流动计算。以轴对称管端出口处激波绕射传播为例，作了数值模拟，得到与实验结果吻合很好的流场精细结构。     

基于流固耦合方法模拟减压器动态特性

为了研究内部流场对减压器动态特性的影响,基于新离散几何守恒律和高精度界面流固耦合技术对其动态特性进行数值模拟,其中单自由度质量弹簧阻尼结构动力学模型采用Newmark法进行计算;三维非定常积分形式的ALE流体控制方程采用弹簧近似的动网格有限体积法进行求解。采用虚拟网格通气技术实现阀门内部运动部件接触/分离过程引起的计算区域拓扑变化。计算结果表明,在减压器正常工作状态下,出口压力存在振荡现象,其均值接近按照静态性能设计的理论值;进口压力增高,振幅减小;进口压力太低时,减压器无法实现出口压力的设计理论值。在正常工作状态的进口压力范围内,减压器存在颤振现象,出口压力远大于设计理论值。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

化学非平衡解耦算法精度分析及其改进

模拟化学非平衡流问题涉及到流动方程和化学反应方程两部分的求解,流动方程组中包含所有组元的偏微分方程,导致求解变量成数量级增加,点隐算法和全隐算法在求解源项刚性问题时,通常会涉及到矩阵求逆运算,这些因素带来的巨大计算量限制了隐式算法在复杂工程问题中的应用。1993年刘君提出了化学非平衡流解耦算法,将控制方程组分解为流动和化学反应两部分,流动方程的求解采用冻结流假设来描述流体微团沿流线的运动过程,化学反应方程的求解描述流体微团在随体坐标系下发生绝热、定容的爆炸过程。结合解耦算法和有限体积法的特点对这种解耦算法进行改进,提出的优化算法不需要求解组元变量所对应的偏微分方程组,只求解由5个基本变量构成的,形式上与量热完全气体近似的偏微分方程组,通过对比计算结果发现优化算法可以显著地提高计算效率。同时,将精细积分方法应用于化学非平衡流问题的求解中,通过与传统的VODE方法和α-QSS方法对比发现,精细积分方法的鲁棒性更优、精度对时间步长不敏感,适当的选取时间步长,可以充分发挥精细积分方法的优势。