两种前体压缩方式对斜爆震燃烧影响的数值研究

为推动工程应用,探究斜爆震发动机前体压缩程度对斜爆震燃烧的影响,本文建立了两道等强激波和斜激波-等熵两种前体压缩简化模型,通过数值模拟对比了飞行马赫数8~10条件下,两种压缩方式对斜爆震波结构以及斜爆震波总压损失的影响。结果表明,前体压缩方式的差异会引起斜爆震波起爆区结构以及起爆位置的变化,且随着飞行马赫数降低,压缩方式对起爆区结构影响减小,对起爆位置影响增大。两道等强激波的前体压缩方式对应的斜爆震燃烧过程总压损失更小,同时可减小点火起爆距离,有利于缩短燃烧室长度。但综合考虑进气压缩与燃烧过程,斜激波-等熵的前体压缩方式对应燃烧室出口气流总压更大。斜爆震发动机的设计需要综合考虑前体压缩与斜爆震燃烧损失以实现发动机总体性能最优。     

流向涡与斜激波相互作用的试验研究和数值模拟

为研究流向涡与斜激波相互作用在超声速燃烧中的应用,进行了由翼产生的流向涡与楔块产生的二维斜激波相互作用的燃烧室冷流试验研究.在不同激波强度下,纹影仪捕捉到了强、中、弱不同的涡/波作用现象.仿真与试验结果符合得较好.试验结果表明:对马赫数为2.3流场,在翼攻角10°时,能产生强流向涡,此工况下,锲角越大,涡/波作用越强.仿真结果表明:马赫数对涡/波作用影响较大,总压影响不明显,总温可影响亚声速回流区的尺寸.      

高超声速二元进气道前体曲线激波逆向设计

为实现高超声速二元进气道前体激波的压缩能力、压缩效率、流量捕获特性和结构长度可控,提出了平面曲线激波逆向设计方法。采用B-Spline曲线控制激波形状。利用有旋特征线法,求解激波的影响域及决定激波的壁面。设计了一凹激波,并对波后流场进行CFD无黏数值模拟,CFD结果和设计结果一致,验证了设计方法的可行性。此外,还设计了直激波、凹激波和凸激波3种激波,并对其在设计点和非设计点处的性能开展了数值研究。在设计点处,分析了激波的压缩比、激波的总压恢复系数、激波压缩区的压升、激波压缩区出口马赫数和流动角随激波控制角的变化规律。在非设计点处,分析了激波压缩区的流量系数和总压恢复系数随攻角和马赫数的变化规律。     

跨声速叶栅抽吸流、激波以及分离流相干效应

以某高负荷、跨声速压气机叶栅为研究对象,应用数值模拟手段探讨通过抽吸控制激波从而控制附面层发展的可行方法。研究结果表明：随着抽吸量的增加吸力面马赫数峰值提高,激波损失增加,同时使得吸力面马赫数峰值点位置后移,附面层分离减弱,分离的减弱所导致的总压恢复系数增加量要远大于激波强度增加所导致的总压恢复系数减小量;抽吸对叶栅性能改善存在一个最佳抽吸量1.2%;在保证叶栅静压压升不变的前提下,相对于未抽吸条件1.2%抽吸使得叶栅总压恢复系数提高10%,扩散因子降低18%,落后角减小5°;通道激波后实施附面层小流量抽吸不能有效改善附面层内部流动参数,当实现前缘入射斜激波投射点位于通道激波上游时,叶表附面层流动得到较大改善。     

二维高超声速压缩面斜激波系的优化配置

选择进气道外压缩段的气流转折角为设计变量,以进气道的压升比为约束条件,分别以总压恢复系数和阻力系数为优化目标,基于一维气体动力学,在飞行马赫数6,飞行高度25 km,进气道出口气流方向与自由来流夹角等于0,°2,°4°和6°的条件下,对二维混压式高超声速进气道进行了优化设计,得到了优化的斜激波配置。优化结果表明,在所讨论的各种情况下,以阻力系数最小为目标得到的外压缩段斜激波系与以总压恢复系数最大为目标得到的外压缩段斜激波系基本相同,都近似为等强度斜激波系。运用数值模拟手段验证了等强度斜激波系配置原则,本文提供的方法可以用于二维混压式高超声速进气道的初步设计。     

基于膨胀波效应的高超声速进气道肩部流动分离控制研究

为改善高超声速进气道唇口激波/附面层干扰诱导的肩部流动分离，从膨胀波及激波理论出发，推导出了膨胀波效应影响下的斜激波附面层干扰理论公式，获得了影响斜激波诱导分离的主要因素：膨胀角梯度、激波角及波前马赫数。在此基础上，开展了膨胀波效应影响下的流动分离控制研究，给出了膨胀波效应影响下斜激波诱导分离的判别及预测方法。结果表明：增大激波入射点处膨胀角梯度，可以显著减小甚至消除肩部流动分离；随着激波角增大，激波强度及逆压力梯度增加，分离区尺寸显著增大。而波前马赫数对分离区尺寸的影响不显著；在进口马赫数3.57～5.18，唇罩角度6°～10°范围内，当激波入射点处逆压比梯度小于250 m-1时，斜激波诱导的流动分离消失，可为改善超声速/高超声速进气道内流道流动分离提供技术支撑。     

双弯曲入射激波的可控中心体内收缩基准流场设计

采用有旋特征线法设计了一种双弯曲入射激波的可控中心体内收缩基准流场,两道入射激波交于中心体起始点,入射激波和反射激波通过给定激波径向总压恢复系数分布进行反设计,壁面通过给定轴向马赫数分布规律进行反设计.该基准流场分为“三波四区”且压缩效率较高.基于该基准流场设计了圆形进口内收缩进气道并进行了黏性修正,数值计算结果表明:内收缩进气道设计点核心区的流场特征和激波形状与基准流场基本一致;在来流马赫数为4.0~7.0时进气道具有较高的压缩效率和良好的流量捕获能力,设计点喉道截面增压比和总压恢复系数分别为17.7和0.729;来流马赫数为5.0~7.0时内部总阻力系数变化平缓,从0.23下降为0.22.      

高超二元曲面压缩进气道前缘激波特性分析

不同的二元进气道前缘激波压缩方式直接影响进气道性能参数,尤其是流量系数。详细分析研究了两种类型的高超二元曲面压缩进气道前缘激波随攻角和来流马赫数变化的特性。计算结果显示,对于初始楔加等熵压缩的二元曲面进气道,在给定的初始楔角下,随着进气道攻角的增大前缘激波与压缩面之间的夹角先减小后增大,随着来流马赫数的增大,前缘激波与压缩面之间的夹角减小;而对于等压力梯度分布控制的弯曲激波压缩二元曲面进气道,随着进气道攻角或来流马赫数的增加,前缘激波与压缩面间的夹角都增大。通过理论和计算分析表明,是斜激波与马赫波不同的相交模式造成了两种曲面压缩进气道的前缘激波随进气道攻角和来流马赫数变化的特性不同。      

构型参数对悬臂斜坡喷注器混合性能的影响

为探索高马赫数下激波诱燃冲压发动机前体/进气道喷注器构型参数对燃料/空气混合的特性的影响,运用隐式RANS仿真,对不同构型的悬臂斜坡喷注器进行了三维数值模拟并进行了对比研究。结果表明,15°压缩角构型最终的混合效率比基准构型高9%且有较大的燃料穿透度,同时会付出30%的总压损失代价;15°压缩角构型的总压损失是10°压缩角构型的两倍;0°膨胀角构型的混合效率始终高出基准构型11%,但有较大的总压损失;后掠构型与基准构型相比混合效率,燃料穿透度基本不变,负的后掠角总压损失最小;三种不同喷注间距的喷注器总压损失大小基本相同,近流场的混合效率相差不大,混合终了阶段,20mm喷注间距喷注器的混合效率达到了较高水平;50～150mm是燃料/空气掺混剧烈发生的区域。     

大攻角范围超声压气机叶栅激波波系研究

对超声压气机叶栅的多攻角工况进行试验,利用纹影仪、油流试验及叶片表面等熵马赫数分布结果进行对比分析,观察到大攻角范围下叶栅激波波系结构发生了明显变化。为揭示激波结构变化原因,利用NUAA计算程序对叶栅进行仿真。研究发现,大攻角状态下叶栅通道中斜激波产生的原因,为前通道激波诱发附面层分离再附后,气流为沿叶片表面继续流动,从而形成斜激波;由于斜激波的增压降速,导致尾缘激波非常微弱甚至消失。     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

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

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

超声速气流中的斜激波诱导自点火

在来流马赫数为2.5的条件下对氢气和乙烯燃料的斜激波诱导自点火过程进行了试验研究.采用纹影和高速摄影获得了自点火初始火核的形成和发展过程.试验考察了不同的来流总温、喷注条件、凹腔长度、斜激波强度等因素对自点火过程和点火边界的影响.研究结果表明:入射斜激波能够提高气流的局部参数起到强化自点火的作用.凹腔内形成的低速回流区有助于初始火核形成后的火核逆流传播,并起到防止初始火核因气流振荡而吹脱的稳定源作用.增加喷注压力、增加凹腔长度、提高来流总温有助于自点火的发生.      

不同反压下椭圆形隔离段流场特征与气动性能

在反压条件下,对椭圆形等直隔离段进行了全三维数值仿真研究.结果表明:在对称的几何条件和边界条件以及均匀来流条件下,当来流马赫数较小时,流场与激波串结构基本对称;当来流马赫数较大(大于2)时,流场与激波串结构出现明显不对称现象,并且随着来流马赫数增大,这种不对称现象加剧.椭圆形隔离段抗反压能力随着长短半轴比减小逐渐增强,出口马赫数和温升比基本不受长短半轴比影响.不同来流马赫数下,当隔离段承受最大反压时,出口马赫数变化不大.      

变比热斜冲波解析解

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

设计马赫数为5的一级二元混压式进气道典型状态下气动特性

对一种设计马赫数为5的一级二元混压式进气道再入大气层过程典型状态进行了仿真和风洞试验,得到了该进气道典型状态下的气动特性.结果表明:当来流马赫数高于设计马赫数（为5）时,进气道外压斜激波系提前汇合,与唇罩入射斜激波相互作用,产生了波-波干扰;尽管发生了流动分离,但当来流马赫数为7和6时进气道出口上游气流紊流度分别不超过3.337和3.256,且流道内动态压力信号的功率谱密度呈现白噪声特征,不会对发动机造成结构损伤.因此,对于宽来流马赫数工作范围的进气道来讲,为了提供足够的流量,可以适当降低进气道的设计马赫数.      

Shock wave configurations and reflection hysteresis outside a planar Laval nozzle

When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection(RR) pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection(MR) pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited condition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these conditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RRMMR transition in the dual-solution domain.