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

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

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

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

超声速可调进气道内流双解现象及其节流特性

为研究超声速可调进气道喉道调节过程中的内流结构及节流特性，设计了工作马赫数范围为0～4的超声速可调进气道，在来流马赫数为2.9的风洞中借助高速纹影观测系统和动态压力测量系统开展了试验研究。结果表明当内收缩比（ICR）为1.79时，进气道通流流场存在设计和非设计流态的内流双解现象。其中，在设计流态下进气道内通道中的波系结构正常建立；而在非设计流态下进气道内收缩段中存在局部分离诱导的复杂波系结构，导致其整体压升高于设计流态，并存在宽频、低频的小幅振荡。此外设计和非设计流态下进气道的下游节流性能相当，其临界堵塞度分别为42.4%和41.7%，临界压比分别为15.8和16.0。在两类流态的下游节流过程中扰动均以结尾斜激波串的形式向上游传播，且临界状态下激波串头波刚好位于喉道附近，但两类流态的结尾激波串在空间分布特征和振荡特性上均存在明显区别。     

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

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

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

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

连续式风洞二喉道调节马赫数控制策略

为了降低连续式跨超声速风洞压力波动,提高马赫数稳定性,需要对二喉道调节马赫数控制方式进行研究,针对现有文献鲜少对该控制策略描述等问题,以0.6m连续式跨声速风洞为例,对二喉道控制马赫数的原理进行分析,基于运动控制器加伺服驱动器双PID(proportion-integral-derivative)控制模式实现二喉道位置精确控制,提出了二喉道和压缩机转速的组合控制流程,并采用分段变参数模糊PID加串级控制的算法实现马赫数精确控制,最后进行了试验验证。结果表明马赫数控制精度优于0.001,且每个马赫数极曲线(9个攻角阶梯)的时间可控制在4min以内,证明所提出的控制策略是有效的,可为连续式跨超声速风洞的设计调试提供参考。      

连续旋转爆轰发动机外流场三维数值模拟

为研究连续旋转爆轰发动机(CRDE)内外流场的变化特性,采用氢气-空气单步有限速率化学反应模型,对内径为40 mm、外径为60 mm、长度为50 mm的连续旋转爆轰发动机进行三维数值模拟,获得了CRDE内外流场结构特征和旋转爆轰波相关参数的变化特性,分析了不同进气总压条件对流场结构和发动机性能的影响。结果表明:爆轰产物在燃烧室出口附近膨胀加速,压力和温度大幅降低,在流场下游产生激波使压力回升,且随进气总压的升高,激波距燃烧室出口距离增加;出口附近羽流中心形成低压高温区域,中心平面上的平均压力低于环境压力,给发动机推力带来了副作用;羽流外围的空气受出口处斜激波的扰动,压力呈现出周期性变化;发动机推力随进气总压的升高而呈线性增加,进气总压为0.55MPa时,发动机推力达到了1160 N。计算仿真结果对掌握连续旋转爆轰发动机外流场特性具有一定的参考价值。     

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

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

烧蚀后C-SiC复合材料隔离段内激波串长度预测

为了研究畸变气流影响下烧蚀后的C-SiC复合材料隔离段内激波串长度预测问题,设计了能够模拟进气道喉道流场畸变的隔离段直连台,并开展了马赫数2.5和3.5来流的两种不同烧蚀程度的C-SiC复合材料隔离段直连实验,获得了有效的纹影和压力数据,结合实验和理论分析,采用无量纲不平度hs/Dhy来表征烧蚀后C-SiC隔离段内壁面严重不规则粗糙特征,尝试在Waltrup公式中增加包含无量纲不平度的因子(1+4hs/Dhy),以此形成新的公式来预测烧蚀后C-SiC隔离段内激波串长度,预测精度较原来公式大幅提高,按照10%的误差带分析,修正后的Waltrup公式预测精度提高了33.3%。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。