基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

左行运动激波主导的管内流动特性研究

为了分析左行运动激波主导的管内流动特征，本文采用非定常数值仿真方法，对亚声速进口条件下等直管道内左行运动激波传播与演化特性、左行运动激波/边界层干扰特征开展研究。研究结果表明：在出口周期性强压力脉动干扰下管内存在连续的左行运动激波，该左行运动激波传播特征具有相似性，激波强度、传播速度按幂函数规律衰减。气流经过左行运动激波后总压、总温、静压阶跃式升高，随后受膨胀波影响气流总压、总温、静压下降；左行运动激波/边界层干扰诱发形成翼型回流区，该回流区随运动激波强度衰减逐渐减小。理论与数值分析表明存在左行运动激波后速度为零和运动激波两侧总压相等的两个临界状态。波前马赫数低于临界值或左行运动激波强度高于临界值时，左行运动激波后为倒流、波后总压高于波前。     

不同直径圆球诱导燃烧的振荡机制与频率特性

为了分析不同直径圆球诱导振荡燃烧的规律，并揭示圆球大小在振荡燃烧现象中所发挥的深层次作用，本文采用二维轴对称Euler方程和基元反应模型，对不同直径的圆球在H2/air预混气体中诱导振荡燃烧的现象开展数值模拟研究。研究发现，振荡频率并不是简单地随直径增大而逐渐从高频向低频连续过渡，而是存在两次突变，形成了超高频、高频以及低频三种振荡燃烧模态。在两种模态间过渡时，振荡达到稳定状态前，会存在一段双频耦合的振荡阶段。三种不同振荡燃烧模态的产生是受到了不同振荡机制的作用，而两种模态间过渡时的双频耦合现象则是两种机制相互竞争的结果。     

背景波系下的隔离段激波串运动特性及其流动机理研究进展

对高超声速进气道-隔离段激波串在复杂背景波系下的突跳运动特性及其流动机理的最新研究进展进行了综述，涵盖了背景波系作用下的激波串运动特性、突跳机理和突跳运动特性的数学描述方法，以期对高超声速进气道相关研究工作提供一定的参考。首先，对固定背景波系和变化背景波系下的激波串运动特性和突跳机制进行阐述，指出隔离段壁面压力顺压力梯度和逆压力梯度的交替变化是激波串突跳特性产生的内在物理机制。其次，对背景流场下隔离段激波串突跳运动的触发机理和触发条件进行了讨论。最后，基于对运动特性和突跳机制的认识，尝试给出了背景波系作用下的隔离段激波串运动特性的数学模型，为激波串前缘位置控制提供参考。     

圆弧激波缓聚点火及后续燃烧传播

采用激波管实验和准一维数值模拟的方法，对预混可燃气体中圆弧汇聚激波的自点火现象及后续燃烧波的传播特性进行研究。其中圆弧汇聚激波由平面运动激波通过精确设计的弧形过渡管段转变得到。研究表明:收缩段中圆弧汇聚激波波后的非均匀梯度环境由激波在平直段、弧形过渡段和扇形收缩段中传播所分别诱导的3个梯度区共同构成。随着圆弧汇聚激波的不断增强，圆弧激波后某处首先形成一个无激波的温和反应区。该反应区逆流锋面的初期运动速度远超Chapman-Jouguet（CJ）爆轰波速，而反应产物区流动则呈现出一定弱爆轰波特征。进一步分析发现，该反应锋面本质上是一种"自发反应波"（spontaneous reaction wave），而非常规意义上的动力学波，其速度与汇聚激波波后气流点火时间梯度的倒数吻合。而后，反应区的扩张速度很快降至CJ爆轰波速以下，伴随反应锋面附近激波的产生以及激波-火焰复合结构的形成。激波-火焰结构最终加速演变为反向传播的爆轰波。在一定的条件下，由于入射激波转变过程和汇聚所构造的特定点火环境，自发反应波可再次赶超爆轰波，成为新的燃烧波前；而当自发反应波速度再次低于CJ爆轰波速时，它将再次转变为爆轰波；在此过程中，原先的爆轰波阵面蜕变为反应产物中传播的激波。     

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

扩压叶栅叶尖非定常流动的涡动力学机理

为了揭示压气机叶尖区旋涡结构与流动非定常性之间的关联，采用URANS对一亚声速平面扩压叶栅在不同攻角下的流场进行了求解，并借助Q判据提取了叶尖瞬态涡系结构。结果表明：泄漏涡的破碎现象能够通过诱导新的涡结构间接作用于相邻通道的叶尖流动，是导致叶尖流场失稳的关键因素。在-0.3°和+0.7°攻角下，叶尖泄漏涡发生了螺旋破碎，并伴随有非定常诱导涡的出现，诱导涡对相邻叶片载荷的影响使得叶尖泄漏涡发生周期性摆动；在+1.7°攻角下，泄漏涡破碎会导致反流涡的形成，反流涡的输运会给叶片载荷和来流攻角带来非定常扰动，反过来又会作用于泄漏涡的破碎和反流涡的生成，最终表现为一种自维持的非定常流动现象。     

基于气体放电的高超声速激波结构显示技术

基于气体放电辐射强度与气体密度的相关性，在高超声速脉冲风洞FD-20中搭建了气体放电流场显示系统，并分别以平板模型、平板-方块模型和简化进气道模型为试验模型，在来流马赫数Ma=12.16、来流静压p≈106Pa的流场条件下开展气体放电流场显示技术研究。在平板实验中，气体放电方法较准确地观测到了电极之间的平板前缘激波结构，与纹影技术测得激波角相差仅为0.21°。在平板-方块实验中，气体放电方法观测到了2个截面（对称面和远离对称面截面）的激波结构，对称面波系结构与纹影和数值计算所得结果基本一致，远离对称面截面的波系结构与数值计算结果基本一致。在简化进气道实验中，气体放电方法观测到了内流道激波交叉形成的菱形结构，且尺寸与数值计算结果相差较小，约为7.9%。这些实验结果表明，在高超声速脉冲风洞中，采用气体放电方法可以获得清晰准确的激波结构，不仅可进行分截面激波结构观测，还可对被模型遮挡的内部区域激波结构进行显示，而且特别适合用于局部复杂流动波系结构的观测。