壁面抽吸条件下的爆震波传播特性

为研究静止气中壁面抽吸对爆震波传播特性的影响,采用数值模拟方法研究了多孔抽吸壁面条件下爆震波的流场结构、传播速度等特性的变化规律。结果表明多孔抽吸壁面对爆震波传播特性有两方面的影响。首先,流场与多孔壁面发生碰撞产生弧形激波,弧形激波对爆震波的横波结构造成直接破坏,导致靠近多孔壁面区域的爆震波产生明显的速度亏损,爆震波强度削弱甚至熄爆;其次,多孔壁面作用导致了流场不稳定性增强,流场与孔板碰撞产生高温高压点,对爆震波的传播有促进作用,尤其在临界条件下,可能导致爆震波熄爆后重新起爆。在保持抽吸压力和边界条件不变情况下,对不同活性及不同抽吸距离时的爆震波传播特性进行研究,发现预混气活性降低、抽吸距离增长时,爆震波结构和传播速度受壁面抽吸影响增强。在加长抽吸距离条件下,随着预混气活性的降低,存在3种爆震波传播现象,即自持传播、熄爆后重新起爆以及完全熄爆。将相应工况下的实验研究与数值模拟结果进行对比,验证了数值模拟结果的正确性。     

来流条件对热流组分扩散项影响效应分析

使用数值模拟的方法研究了高超声速热化学非平衡流动中不同来流条件对热流组分扩散项的影响。以RAM-CII飞行试验外形为例,使用AUSM+-up格式耦合LU-SGS方法求解带化学反应源项的多组分NS方程。得到如下结论:在非平衡热环境数值模拟研究中,完全催化壁时组分扩散项热流所占比率表征了壁面催化效应的强弱;在同一高度下,随着马赫数升高热流组分扩散项比例越来越显著,而马赫数相同时,随着高度增高热流组分扩散项所占比例越来越小,主要原因在于当地化学反应进行程度不同。     

用于可压缩壁面湍流LES及DNS入口边界条件的"回收-调节"方法综述

为了减小湍流流动问题的大涡模拟和直接数值模拟的入口湍流脉动边界条件对计算结果的影响,叙述了用于可压缩壁面湍流入口脉动生成的一种方法--"回收-调节"方法(recycling-rescaling method),叙述了这种方法的理论基础及简要推导过程,介绍了几种常用的"回收-调节"方法并指出了实现中的几个关键问题,分析了...     

二维有限长度柔性壁面上T-S波演化的数值研究

受自然界鸟类羽毛的柔性特征启发,利用数值模拟的方式开展了柔性壁面对亚声速边界层中T-S波演化的影响研究。刚性壁面上的数值结果与线性理论吻合得很好,验证了数值方法的可靠性。在此基础上,将部分刚性壁面替换为柔性壁面,结果表明,柔性壁面可以抑制T-S波在空间上的增长,从而推迟边界层流动转捩。壁面的变形不只跟随T-S波的波形,还因为柔性段与刚性段相接的前缘和后缘引起与扰动源频率相同的更大尺度的壁面波动,壁面的实际变形由这几种波叠加而成。开展的参数研究结果表明,增大表面的质量密度对于柔性壁面衰减扰动的效果几乎没有影响;增大表面张力和增加底部支撑的弹性系数可以增加壁面的刚性,减小壁面变形的幅度;增加阻尼可以抑制柔性段前后缘产生的大尺度壁面波动的传播,而对跟随T-S波的变形影响不大。总体上,柔性壁面的变形程度越大,其扰动的抑制效果越强。      

带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

前体边界层状态对高超声速进气道流动结构及性能的影响

采用CFD方法研究了前体边界层状态对轴对称高超声速进气道的波系结构、激波/边界层干扰特性以及流量捕获能力、总压恢复系数等的影响,获得了不同马赫数下进气道性能参数随一级锥边界层强制转捩位置的变化规律,并对采用二维粗糙体作为强制转捩措施给流场所带来的干扰进行了分析,研究中还对比了具体进气道设计方案间的差别.      

基于欧拉方程的高超声速飞行器的壁面流线生成计算

提出了一种新的利用表面流函数法计算高超声速飞行器表面流线分布的方法.首先提出了表面流函数的概念,并通过理论推导,得到了表面流函数与表面流线的关系;然后运用结构化网格求解三维Eu-ler方程,计算得到高超声速钝头体的边界层外缘外部无粘流场气流参数;最后利用无粘流场气流参数和表面流函数的方法计算了高超声速飞行器的精确表面流线分布.结果表明,在无攻角和攻角小于20°的情况下均可以得到较好的壁面流线分布.高精度的表面流线的得到,为利用边界层内粘性主导区域的积分方程等方法进一步精确预测高超声速飞行器表面的气动加热奠定了基础.      

基于线性稳定性理论的转捩控制试验

 使用线化稳定性理论研究了某后掠机翼试验模型在风洞试验条件下可能的试验结果。在当前试验条件下,三维边界层中横流驻波的不同波长扰动发展趋势不同,根据理论分析和数值计算给出了最大放大率波长为4.0 mm。通过升华法试验验证了理论分析与数值计算结果。在试验中使用丝网印刷技术在后掠机翼前缘添加粗糙带,成功地引入了理论计算的最大放大率波长。根据预期影响了三维边界层的转捩。     

非对称来流下隔离段动态压力特性

为了探讨非对称来流下矩形隔离段内动态压力特性,用直联实验方式以及动态压力测量技术进行了试验,并用统计分析方法分析了数据。实验结果表明,入口来流的非对称性对上下壁面压力脉动大小以及传播平均速度有较大影响,而对频率无多大影响。激波串区域内壁面压力脉动向下游传播比向上游传播衰减得快。在隔离段出口超声速条件下,观测到的压力脉动频率主要在70 Hz以下,而在出口亚声速条件时,压力脉动的频率不仅有70Hz以下的部分,而且还有100~200 Hz之间的部分。     

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.