预混气爆轰马赫反射实验研究

为研究斜爆轰马赫反射的实验现象以及不稳定性对马赫反射的影响,在矩形爆轰管道中,对3种预混气的马赫反射现象进行了实验研究,取得烟膜记录。烟膜结果显示,三波点高度在楔角附近与无反应三激波理论类似,随后逐渐偏向反应三激波理论,三波点实际高度由反应实际放热量决定,三种气体偏转始末距离与结束偏转距离之比平均值分别为0.42,0.49,0.31,与三者活化能大小关系相符合。胞格结构随初始压力增大而渐趋规律,胞格尺寸也逐渐减小,由CJ区域到过驱区域胞格结构存在明显的转变。横波是影响斜爆轰稳定的重要因素,压力增大可抑制分子的无规则运动,使横波出现的位置更加集中,胞格更为规律。预混气爆轰不稳定度大的气体横波出现的时空位置更不易确定,胞格结构更不规律,胞格尺寸较小。     

振动激发对激波反射的影响

定常激波反射分为规则反射和马赫反射,在不同条件下2种反射结构之间会相互转变。高超声速流动中的激波反射问题常面临高温气体效应,随着温度逐渐升高,最先出现的是空气分子振动激发。通过理论分析和定量计算,研究了振动激发对激波反射及其转变规律的影响。首先给出考虑振动激发的空气热力学模型,并分析其与量热完全气体的差异以及对激波关系的影响;接着分析在规则反射和马赫反射中,振动激发对激波反射流场的影响规律;最后讨论振动激发对激波反射转变2个准则点的影响。研究结果表明,振动激发使激波极线的整体轮廓变大,且这种差异在经过一次激波反射之后被明显放大,会对激波反射的流场产生重大影响;对于激波反射的转变,振动激发使转变的2个准则点都变大,且对规则反射向马赫反射转变的脱体准则影响更大。     

气相爆震波马赫反射转变机理的实验研究

在矩形爆震管道中,选用具有规则胞格结构的高浓度氩气稀释稳态气体C2H2-2.5O2-8.17Ar及不规则胞格结构的非稳态气体C2H2-5N2O,利用纹影和烟膜实验对CJ爆震状态下的爆震波发生马赫反射的转变机理进行了研究。在楔角为30°,5k Pa初始条件下,C2H2-2.5O2-8.17Ar马赫反射转变时马赫杆临界高度约为1.25cm,过驱胞格结构出现的最小距离略小于胞格长度的1/3;C2H2-5N2O马赫反射转变时马赫杆临界高度约为1.5cm,由于其胞格结构不规则,故难以衡量过驱胞格结构出现的位置。实验结果表明:爆震波马赫反射三波点运动轨迹特征为在初始阶段遵循无反应冲击波理论,随后马赫反射发生转变使得运动轨迹逐渐平行于反应冲击波理论,从而验证了爆震波开始形成前导激波马赫反射,由于CJ区域内三波点的介入引起扰动,最终在马赫反射区形成过驱爆震的过程。     

拉瓦尔喷管外发生激波反射工况详细分析

平面拉瓦尔喷管的环境压力与入口总压的比值从设计工况到第三极限工况(管口处为正激波)逐渐增大时,喷管外的激波现象不同.根据已知的激波反射理论,将设计工况到第三极限工况间的工作状况进一步细分为强斜激波工况、马赫反射工况、双解区和正规反射工况、并将由气流偏转角所表示的发生正规反射和马赫反射的条件用环境压力与总压的比值来表示,找出了强、弱斜激波的环境压力与总压比值的临界值,得到这些临界值随着喷管面积比的变化趋势.最后针对面积比为5的平面拉瓦尔喷管,由理论分析得到其管外发生正规反射、马赫反射、强斜激波反射时对应的压比分别为:0.0841,0.0959和0.2005,并用数值模拟来验证理论分析的正确性.      

气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.      

用于爆震引燃的激波聚焦无反应流场数值模拟

通过对用于脉冲爆震引燃的内楔表面激波聚焦冷态流场的数值模拟，分析了不同角度内楔表面在激波聚焦时流场结构的不同特点。数值计算的结果表明：对于大内楔角的情况，激波进入内楔角之后产生规则反射，温度的大幅度升高发生在激波于内楔顶点反射之后；对于小内楔角的情况，激波进入内楔角之后产生马赫反射，温度的升高在入射激波到达内楔顶点之前就已经发生，而在入射激波于内楔顶点反射之后得到进一步的提升。同时，还对不同角度内楔的反射激波后气流温度进行了比较。     

斜爆轰波系在受限空间内的演变及其临界条件的数值研究

针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析了斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析了爆轰波系结构转变的临界条件.研究结果发现,随着楔面角度的增加,依次出现四种结构:激波诱导燃烧、斜爆轰波双规则反射、回流区马赫反射、楔面燃烧.对于稳定的波系结构,楔面压缩角同时存在下临界...     

展向宽度对激波入射平板反射类型的影响

激波干扰是高超声速飞行器气动布局和超燃冲压发动机设计中需重点考虑的局部干扰现象,当该现象发生时会产生复杂的波系结构,对流场行为特征产生影响,进而对飞行器物面载荷及发动机性能产生显著影响。采用数值计算方法针对斜激波入射平板问题,在固定来流马赫数5、单位雷诺数7.12×10⁶ m^-1不变的条件下,通过改变上下平板展向宽度研究了6个不同状态下的激波反射类型。结果表明在受限空间内必须考虑侧向溢流对激波反射类型的影响。随平板展向宽度增加,激波反射类型从正则反射逐渐过渡至马赫反射,且马赫杆长度变长并逐渐前移,直至导致内流道堵塞形成脱体激波。采用激波极曲线方法在激波入射角度固定的条件下对两种激波干扰类型的产生机制进行了分析,发现随激波强度增加两束透射激波极曲线上移缩小,进而造成波后流场参数匹配需求,激波反射类型逐渐从正则反射向马赫反射转变,得到了与数值计算结果一致的结论。     

球面激波对平面的马赫反射的简化处理

本文利用Whitham绕射理论,给出了球面激波对平面绕射时马赫杆的近似解,同时求得了马赫反射三波点迹线及壁面超压。另外还导出了新的面积函数表达式。在壁面超压的计算上,这些理论结果与用Whitham方程求得的数值解相符合。     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

斜激波总压损失率极小值理论解与物理意义

以法向马赫数作为激波强度表征量,对斜激波关系式进行重新推导,得到了穿过激波总压损失率极小值的理论解。控制方程表达式为激波角对物面角的线性函数。依据斜激波总压损失率极小值解析公式,首先,绘制了针对超声速流总压损失率应用的楔形角-激波角-马赫数的斜激波效率图。其次,通过生成斜激波三维总压损失率等值线图,呈现了总压损失率在楔形角-激波角-特征马赫数空间上的分布规律。此外,利用斜激波效率图,揭示了等总压损失率条件下马赫数与激波角的对称双解现象。     

高超声速内收缩进气道轴对称基准流场改进

针对轴对称基准流场中前缘曲激波靠近中心体的部分激波强度过大现象,基于马赫数分布可控反设计方法,将这道前缘曲激波分解为一道较弱弯曲激波和部分等熵压缩波,改进的基准流场存在“四波四区”结构且压缩效率明显提高.基于该改进的基准流场和常规“两波三区”基准流场分别设计了圆形进口的内收缩进气道并对其流场特点和性能进行数值研究.结果表明:改进的进气道的流场能较好保持基准流场的特点;在来流马赫数为4.0~7.0范围内具有较高的压缩效率和良好的流量捕获能力,设计点的出口压比和总压恢复系数分别为17.56和0.540;改进的进气道采用来流马赫数从高到低前缘弯曲激波和汇集的等熵压缩波依次封口的设计概念,在提高流量捕获能力的同时减小了总压损失,总体性能优于常规进气道,来流马赫数为7.0时总压恢复系数相对提高了23.6%,来流马赫数为4.0时流量系数相对提高了5.7%.      

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

Influence of shock wave impinging region on supersonic film cooling

Shock waves can significantly affect the film cooling for supersonic flow and shock waves may have different influence when impinging in different regions. The present study numerically compared the results of shock wave impinging in three different regions and analyzed the effect of impinging region. The shock wave generators were located at x/s = 5, 25, 45 with 4°, 7° and 10° shock wave incidence. The mainstream Mach number was 3.2 and the coolant Mach number was 1.2 or 1.5. The numerical results illustrated that the shock wave impinged in the further upstream region led to a larger high-pressure region and a larger vortex in the boundary layer. Moreover, placing the shock wave generator upstream resulted in the lower mass fraction of coolant in the downstream region. The velocity in the upstream part of the cooling layer was lower than the midstream and downstream part, which resulted in the less ability to resist the shock wave impingement. Therefore, the upstream impingement deteriorated the cooling performance to a greater extent. The study also manifested that the stronger shock wave had a larger effect on supersonic film cooling. Increasing the coolant inlet Mach number can increase the blowing ratio and reduce the mixing, which was of benefit to improve cooling effect.     

重型燃气轮机高雷诺数CDA叶型转捩特性数值计算

为研究重型燃气轮机的压气机叶片在高雷诺数工况下的气动性能,基于Gamma-Theta转捩模型的雷诺时均方程对某可 控扩散叶型进行了数值计算。通过对比不控制马赫数与控制马赫数,分析高雷诺数对可控扩散叶型气动性能及转捩特性的影响。 结果表明：在不控制马赫数条件下,在零攻角时,雷诺数从7×10 5 增大为9×10 5 ,总压损失增加了约391.95%;在高雷诺数工况下随 着雷诺数的增大,叶片流动损失不断增大,叶片可用攻角范围减小,同时在叶片吸力面出现激波,干扰转捩的产生。在控制马赫数 条件下,当Ma=0.6时,在零攻角工况下,雷诺数从8.2×10 5 增大为1×10 7 ,总压损失减小了约38.98%,吸力面转捩起始点从4.78%弦 长处前移至1.11%弦长处;在高雷诺数工况下,叶片流动损失随着雷诺数的增大不断减小,吸力面转捩位置前移。     

不同飞行状态下固体火箭发动机尾喷焰数值研究

为了研究飞行状态对固体火箭发动机尾喷焰的影响,建立了含化学反应项和组分输运项的N-S(Navier-Stokes)控制方程,并采用MUSCL(monotonic upstream-centered scheme for conversation laws) Roe格式的有限体积法进行求解.在此基础上,结合热力计算结果对某固体火箭发动机在不同飞行状态下的化学反应流和单一组分流进行仿真计算.结果表明:初始倾角和第一马赫波节长度随飞行高度的增加而增大,随来流马赫数的升高而减小.同单一组分流相比,化学反应流对第一马赫波节的影响较初始倾角显著.射流轴线温度因飞行状态的不同有较大涨落差异.初始倾角、第一马赫波节由来流总压决定,射流轴线温度由激波结构、复燃效应、来流条件等因素共同决定.      

Subsonic impulsively starting flow at a high angle of attack with shock wave and vortex interaction

Impulsively starting flow, by a sudden attainment of a large angle of attack, has been well studied for incompressible and supersonic flows, but less studied for subsonic flow. Recently, a preliminary numerical study for subsonic starting flow at a high angle of attack displays an advance of stall around a Mach number of 0.5, when compared to other Mach numbers. To see what happens in this special case, we conduct here in this paper a further study for this case, to display and analyze the full flow structures. We find that for a Mach number around 0.5, a local supersonic flow region repeatedly splits and merges, and a pair of left-going and right-going unsteady shock waves are embedded inside the leading edge vortex once it is sufficiently grown up and detached from the leading edge. The flow evolution during the formation of shock waves is displayed in detail. The reason for the formation of these shock waves is explained here using the Laval nozzle flow theory. The existence of this shock pair inside the vortex, for a Mach number only close to 0.5, may help the growing of the trailing edge vortex responsible for the advance of stall observed previously.     

变比热斜冲波解析解

首次推导了变比热斜冲波的基本公式及其衍生形式，分别就求解转角、波角、最大尖劈角、最小马赫数等典型问题，拟就解法及双精度计算机程序，一一作出例题，进入了变比热斜冲波园地。可以预见到的深入研究问题不少，尚待进一步探索。