膨胀波对斜劈稳定斜爆震波特性影响的数值研究

为探究燃烧室内可能存在的膨胀波对斜爆震稳定燃烧的影响，基于块结构自适应网格加密的AMROC程序，求解多组分可压缩化学反应流Euler方程，研究了膨胀波对斜劈诱导斜爆震流场的影响。发现无限长斜劈模型中，斜爆震波面角度是缓慢增加的，膨胀波影响下，斜爆震波面角度明显下降。从爆震波面角度变化中可以明显看出膨胀波的影响范围和爆震波的衰减程度。来流的低静压不会改变膨胀波影响范围，但容易导致爆震波衰减直至发生解耦，因此高马赫数低静压来流条件下斜爆震燃烧的稳定性和燃烧效率需要着重考虑。此外，使用普朗特-梅耶膨胀波基本原理对膨胀区进行分析，发现爆燃区流场参数与理论值吻合较好，且壁面附近前马赫线角度与近似膨胀波前沿较为接近；基于此，发展了一种定性评估斜劈末端膨胀波影响范围的手段，在膨胀波前马赫线角度基础上适当增加4°～10°，可以近似得到膨胀波的前沿位置。     

爆震波在开缝式弯管内的传播特性研究

为了解旋转爆震波在单侧部分受限情况下产生的波系结构及环缝内压力波动特点,以恰化学当量比的氢气/氧气混合物作为反应物,采用9组分19步的基元反应模型,忽略粘性、热传导和扩散等输运效应,对爆震波在开缝式弯管内的传播情况进行了数值模拟。研究结果表明,爆震波在开缝式弯管内传播时,爆震波及爆震区内的反射激波后的高压区都会在环缝内产生透射激波,但环缝内激波不单纯依赖弯管内高压区的侧向膨胀产生,且环缝内的一些反射激波产生的波后高压区还会反过来在爆震区内产生透射激波;环缝内第二道斜激波会被接触面分为两部分,在接触面前方的惰性气体一侧具有更大的倾斜角;环缝位于外侧时爆震波的倾斜程度比位于内侧时更大;当环缝处于外侧和内侧时,会由于激波相交、反射产生比第一道斜激波更大的压力峰值,内侧开缝时环缝内最大压力峰值约为第一道斜激波后的两倍,且内侧开缝时环缝内压力大于4倍填充压力的高压区域范围比外侧开缝时大36.4%。     

环形射流初始压力对激波聚焦起爆的影响分析

为了研究射流入射压力对环形激波聚焦爆震起爆的影响,以氢气和空气混合物为例,对不同射流入射压力下环形激波聚焦爆震起爆过程进行了数值模拟,结果表明,当环形射流入射强度足够大时,其形成的激波聚焦能够形成高温、高压区域,从而直接起爆爆震波;爆震管推力壁对激波具有反射加强作用,有助于爆震波的形成;环形射流入射压力存在一个临界值,低于此值时,则不能起爆爆震波.      

激波绕射触发爆震波的试验研究

在单爆震试验平台上对激波绕射触发爆震波的特性进行了试验研究.针对不同的激波反射器结构尺寸,测量了爆震管内气体压力变化历程及火焰传播速度,获得了不同工况下激波绕射触发爆震波的特性.结果表明:激波反射器对爆震波的形成有促进作用,并取决于反射器中的反射锥及反射孔板的结构与尺寸;在反射锥堵塞比一定的情况下,存在最佳反射锥锥角使得爆震波的形成距离最小;而反射锥与反射孔板的间距增大,爆震波形成距离增大,一旦形成稳定的爆震波后,反射器对爆震波压力的影响不大.      

当量比和间隙尺寸对爆震波传播过程的影响

为探索微型脉冲爆震发动机推进系统的可行性,进行了微空间内爆震极限特性的研究.通过对预混氢气/氧气爆震火焰在平板狭缝中的传播过程进行测量,分析了爆震波压力、速度随当量比和间隙尺寸的变化趋势,发现在极限范围外,来自稳流段的爆震波都能在间隙内通过压力和速度的自我调整后达到稳定状态,调整过程所需的距离随间隙尺寸减小而增长.根据爆震波速度衰减的定义,给出了激波和火焰锋面速度沿间隙通道方向的变化特点,观察到了微尺寸下爆震波的4种传播模式:稳定爆震波、准稳定爆震波、低速爆震波和非爆震波.      

脉冲爆震发动机快起爆的二维数值模拟

为研究脉冲爆震发动机的点火过程，用二维欧拉方程和Korobeinikov两步爆震模型对爆震波快起爆过程进行了数值模拟研究。用高温高压未燃气体点火方式来模拟试验中的电火花塞点火。研究表明：电火花塞的能量不足以直接点燃爆震波，爆震波是在经历了一系列反射激波的相互作用后最终建立的。激波．爆震波转捩是快点火过程，是爆震波建立的一个非常重要的机理。在设计点火装置时，火花塞要尽可能靠近壁面和拐角，以便产生高强度的反射激波。同时还要控制反射激波的方向，使之向未燃气体中顺利传播，才能起到很好的点火效果。     

气相爆震波衍射现象的数值研究

对不同条件下气相爆震波在突扩管道和带扩散段管道中的衍射过程进行了数值研究.结果表明:起爆管的直径和临界爆震管直径对爆震波能否成功衍射传入主爆管起到决定性的作用.壁面反射激波与爆震波相互作用在邻近壁面处产生高温、高压区域,有助于减少爆震的诱导时间以及横波间距.扩散段扩散角越小时,反射激波强度越大.高强度的反射激波会将爆震波分为3段,在分界点处产生高温、高压点,对爆震波的衍射传播起到支持作用.      

混气活性对爆震波在突扩段传播特性的影响研究

为研究爆震波经过突扩段的传播特性,实验采用低活性混气C2H2+2.5O2+70%Ar和高活性混气C2H2+2.5O2,测量了爆震波在扩张比为1.34的圆管间传播的速度变化,并运用化学反应动力学分析爆震波动力学参数。结果表明:爆震波在小管出口的速度至多降至0.6倍的Chapman-Jouguet爆震波速。高活性混气触发的爆震波的成功传播是由于在大管产生局部爆炸,随后形成了过驱爆震波。相应地,低活性混气触发的爆震波不会形成过驱爆震波,且波速的波动范围较小。当预混气体压力小于爆震波从小管向大管成功传播的最低压力时,低活性混气触发的爆震波更快地衰减为爆燃波。通过化学反应动力学分析可知,高活性混气的爆震波稳定性参数较高,更容易形成局部爆炸。在最低压力处两种混气的爆震波诱导区长度接近0.8mm,管径与胞格宽度之比接近于1,因此爆震波在圆管中传播的准则可以用于预估爆震波在扩张比为1.34的圆管内能否成功传播。     

鼓包诱导斜爆震波的数值研究

为实现低马赫数条件下斜爆震波的快速起爆,以及提高斜爆震波对来流扰动的抗干扰能力,提出了通过在诱导区内设置一个小鼓包来控制斜爆震波起爆的方案,并通过数值求解二维多组分欧拉方程研究了楔面和鼓包诱导的斜爆震波。结果表明：通过在诱导区内设置鼓包可实现斜爆震波的提前起爆;且在较宽马赫数范围内,其起爆位置相对固定。鼓包诱导的斜爆震波主要有强耦合和弱耦合两种起爆形式,起爆形式与鼓包和斜激波的相对位置有关。此外,来流扰动会使楔面诱导的斜爆震波剧烈振荡,不利于斜爆震发动机的稳定工作;鼓包的存在对来流扰动所引发的斜爆震波振荡具有一定的抑制作用。     

高温超声速流中爆震波特性数值研究

通过对高温超声速流中爆震波性质的研究,评估其在高超声速冲压发动机燃烧室的燃烧组织中应用的可行性,并通过数值模拟对分析结论进行了验证。提出了一种新的爆震波起爆机制,注入高温超声速流中的燃料混气可通过自身缓慢的释热使流动进入局域热壅塞状态,进而借助局域热壅塞产生的激波实现爆震波的起爆。计算结果表明在适宜的温度与马赫数条件下,注入高温超声速流中的燃料可通过新的起爆机制在超声速流中形成一道稳定的驻定爆震波。表明在高超声速冲压发动机燃烧室中存在着通过驻定爆震波实现火焰稳定的可能性。     

Oblique detonation wave triggered by a double wedge in hypersonic flow

Pressure-gain combustion has gained attention for airbreathing ramjet engine applications owing to its better thermodynamic efficiency and fuel consumption rate. In contrast with traditional detonation induced by a single wedge, the present study considers oblique shock interactions attached to double wedges in a hypersonic combustible flow. The temperature/pressure increases sharply across the interaction zone that initiates an exothermic reaction, finally resulting in an Oblique Detonation Wav...     

双楔面诱导的斜爆轰波阵面的磁流体控制

基于高超声速飞行器爆轰推进发展需要,数值研究了由双楔面诱导的斜爆轰波类型及其磁流体(MHD)控制的可行性,探讨了后楔倾角变化对双楔斜爆轰流场与MHD控制的影响.结果表明:对于前、后楔倾角分别为15°与20°的双楔面诱导的斜爆轰,Lorentz力可以将不同来流马赫数条件下的稳定斜爆轰波阵面恢复到设计位置,并可以使不稳定斜爆轰波趋于稳定,但是,无法将不稳定斜爆轰波恢复到设计位置.而当后楔倾角稍微增加时,若爆轰流场稳定,同样可对斜爆轰波阵面进行有效控制.但当后楔倾角大于25°时,稳定斜爆轰波流场失稳,此时MHD无法控制流场稳定性.但若增加气体反应速率,虽然斜爆轰波不稳定,MHD控制却可使其趋于稳定,但无法将爆轰波阵面恢复到设计位置.      

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

作为一种面向未来高马赫数吸气推进技术的新概念发动机,斜爆轰发动机的运行依赖其稳定的宏观波系结构。本文针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析爆轰波系结构转变的临界条件。研究结果发现,随着楔面角度的增加,依次出现四种结构：激波诱导燃烧,斜爆轰波双规则反射,回流区马赫反射,楔面燃烧。对于稳定的波系结构,楔面压缩角同时存在下临界、亚临界以及超临界三种极限条件,在波系演变过程还伴随激波规则反射和马赫反射的转变。     

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

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

Investigation of hot jet on active control of oblique detonation waves

The hot jet injection is utilized to actively control the oblique detonation wave, such as initiating and stabilizing an oblique detonation wave at a desired position that is shorter than the length of induction zone, and adjust the height of the oblique detonation wave at the exit of combustor when the oblique detonation wave engine is working on off-design flight conditions. The fifth order Weighted Essentially Non-Oscillatory (WENO) scheme and a two-step reversible reaction mechanism of the stoichiometric H₂/Air are adopted in the simulations. With the help of hot jet injection, the transition from inert oblique shock wave to the oblique detonation wave immediately occurs near the position of hot jet injection, and consequently the length of combustor can be reduced. The angle of oblique detonation wave also decreases as the hot jet injection approaches the nose of wedge. Additionally, the height of the oblique detonation wave at the exit of combustor can be flexibly adjusted, and also depends on the injection position and the strength of the hot jet. If the velocity of the hot jet is too weak to directly trigger the overall oblique detonation wave at the position of injection, increasing the injection pressure will improve the strength of the hot jet and results in a successful transition.     

斜爆轰波的波角和法向速度-曲率关系初探

为研究斜劈诱导斜爆轰波的波阵面弯曲效应，以期为斜爆轰的不稳定性及其演化规律提供新的见解，基于加权本质无振荡（WENO）格式空间离散和附加Runge-Kutta方法时间离散的求解器，针对不同的化学反应参数（释热量、放热速率和化学反应区参考长度）条件，开展斜爆轰波的数值计算研究。结果表明斜爆轰波沿波阵面的波角变化可分为3个区域：区域I，波角平滑减小；区域II，波角跃升后衰减；区域III，波角有规律振荡。波阵面法向速度-曲率关系在区域I呈现准垂直直线变化趋势，并伴随着爆轰波强度的不断衰减；在区域III则呈现出"D"形曲线，即由极曲线段、光滑水平变化段和拟线性变化段组成，为类胞格结构的周期性演变；区域II可认为是以上两个区域特征的耦合。不同的化学反应参数对斜爆轰波波阵面的弯曲效应影响存在较大差别。     

Initiation characteristics of oblique detonation waves from a finite wedge under argon dilution

In this paper, the flow field characteristics of Oblique Detonation Waves (ODWs) induced by a finite wedge under argon dilution are studied by solving the Euler equations with a detailed chemical model of hydrogen and air. First, the effects of the expansion waves, argon concentration, geometric parameters, and Mach number on the ODW are discussed. The results show that the changes of these parameters may make the oblique detonation not be initiated. Then, the ODW initiation criterion of the finite wedge is summarized, as the characteristic length of the induction zone L_C and the characteristic length of the oblique wedge L_W meet the condition L_C/L_W < 1, the initiation of the ODW occurs; otherwise, it does not occur. What’s more, the Constant Volume Combustion (CVC) theory is applied to study the characteristic length of induction zone. It is found that CVC theory is more suitable for the “smooth transition” type of ODW flow field, the theoretical and numerical characteristic length in induction regions are in good agreement. This work is of great significance for the design of oblique detonation engines for hypersonic vehicles.     

环形爆震室中火焰加速数值模拟及验证

对带不同数量孔板的环形爆震室进行数值模拟,并通过试验对数值计算进行验证,来研究火焰加速现象、爆燃向爆震转变过程和不同当量比下起爆距离.数值计算采用二维轴对称非定常Navier-Stokes方程来模拟流体动力学过程.研究发现用较低的点火能量对火焰混合区的可燃气点火产生低速火焰,低速火焰向环形爆震室射流并改变方向向出口传播,火焰在孔板的阻碍作用以及火焰诱导激波和反射波的加速作用下,由层流变为湍流,湍流火焰与其诱导激波相互加强,最终引爆未燃混气;还对爆震波在孔板区的传播过程进行了分析,对不同当量比下的火焰速度和起爆距离进行了模拟研究.      

连续旋转爆震波细致结构及自持机理

基于化学非平衡流解耦方法和详细化学反应机理,在0.1mm的网格尺度下对连续旋转爆震波内流场和爆震波平掠惰性气体界面的传播过程进行了数值模拟。详细分析了连续旋转爆震波的空间结构,由于在出口侧无几何约束,受一系列稀疏波的影响,旋转爆震波强度受到不同程度的削弱,形状也发生了相应的变化。爆震波平掠惰性气体界面时的流场分布与连续旋转爆震波结构极为类似,对比研究了波前可燃气体层高度对爆震波结构和传播过程的影响,当爆震波高度较小时,其抵御惰性气体侧稀疏波衰减的能力弱,波阵面严重变形,爆震波解耦。若爆震波要稳定自持传播,其必须要达到一定的临界高度来抵御稀疏波的影响。