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

在矩形爆震管道中,选用具有规则胞格结构的高浓度氩气稀释稳态气体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区域内三波点的介入引起扰动,最终在马赫反射区形成过驱爆震的过程。     

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

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

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

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

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

气相爆轰波传播特性的数值模拟及实验对照

本文应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波进行了二维数值模拟.氢氧混合物的化学反应模型考虑了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.本文首先对爆轰波数值结果和实验结果进行了对照验证,然后对爆轰波在楔面反射由马赫反射向规则反射转变的过程进行了数值分析,得到了反射转变临界角,并和实验结果及理论分析进行了比较,结果是令人满意的;本文还对爆轰波的多波结构进行了初步的数值分析.     

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征,基于带化学反应的Navier-Stokes方程,对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式,时间项采用二阶Runge-kutta方法。结果表明,对于弹头：（1）在亚爆轰条件下,能够模拟氢气-空气预混诱导爆轰流场的精细结构;（2）在超爆轰条件下,通过精细调整网格,能够很好地分辨强烈耦合的激波和燃烧波,且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构,预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好;通过对流场波系结构变化过程的研究,获得了流场三波点随时间的演化过程。     

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

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

高温区碰并诱导气相爆轰二次起爆的数值研究

针对心形管内气相爆轰二次起爆过程进行了数值研究。采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维欧拉方程,采用基元反应来描述爆轰化学反应过程。结果表明:气相爆轰在心形管中传播,由于壁面几何结构的作用,爆轰波熄灭后会发生二次起爆。该二次起爆由高温区碰并所诱导,且在管道不同位置诱导了两处二次起爆。爆轰流场中存在未反应气囊、射流、爆轰波马赫反射、激波与火焰作用及其所诱导的不稳定性等现象。计算得到的二次起爆过程胞格结构演变与实验基本一致。     

烟板法研究马赫反射三波点轨迹

本文介绍激波管中用烟板法测量马赫反射三波点轨迹的一种实验方法。在ST130×130毫米激波管中的实验研究证实了该方法是成功的。这种方法简单而且易行。     

狭缝内乙烯/氧气预混气体爆轰几何极限的实验

为研究狭缝内爆轰波传播极限,实验得到了不同初始压力(0.004~0.04MPa)下化学恰当比的乙烯/氧气预混气体在狭缝高度为1.0~4.0mm狭缝内的爆轰性能.采用烟膜板记录爆轰波运行轨迹,高速摄影捕捉火焰面.结果表明:狭缝高度越小,爆轰极限对应的临界初始压力越高.对于近极限条件内不稳定爆轰传播模式,包括"结巴"爆轰和驰振爆轰,其对应的初始压力范围随着狭缝高度的降低而变宽.考虑初始和边界条件,将水力直径与胞格宽度之比作为合理的爆轰敏感性参数描述预混气体爆轰特性,得出在不同狭缝通道内爆轰极限范围为,即该比值在0.326~0.403之间.      

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

尖劈诱导的斜爆轰波的精细结构及其影响因素

为了揭示斜爆轰的精细结构及其影响因素,基于带化学反应的Euler方程,利用五阶WENO格式,对尖劈诱导的斜爆轰进行了数值模拟。结果表明,在劈面的压缩下,斜爆轰波阵面被分为三种结构,依次是类ZND结构、单三波点阵面结构和双三波点阵面结构。同时,根据三波点的轨迹绘制了胞格结构,其中单三波点的轨迹为平行直线,而双三波点的轨迹为倾斜的蜂窝状结构。讨论了来流马赫数和尖劈角度对劈面压缩性的影响,随着马赫数和劈角的增加,尖劈的压缩更厉害,爆轰波阵面也就越光滑。     

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

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

A theoretical method for solving shock relations coupled with chemical equilibrium and its applications

In this study, a theoretical method is proposed to solve shock relations coupled with chemical equilibrium. Not only shock waves in dissociated flows but also detonation waves in combustive mixtures can be solved. The global iterative solving process is specially designed to mimic the physical and chemical process in reactive shock waves to ensure good stability and fast convergence in the proposed method. Within each global step, the single-variable equations of normal and oblique shock relatio...     

Viscosity effects on weak irregular reflection of shock waves in steady flow

The viscosity effects on weak shock wave reflection are investigated with the Navier–Stokes and DSMC flow solvers. It is shown that the viscosity plays a crucial role in the vicinity of three-shock intersection. Instead of a singular triple point, in viscous flow there is a smooth three shock transition zone, where one-dimensional shock jump relations cannot be applied. At the flow parameters corresponding to the von Neumann reflection, when no inviscid three-shock solution exists, the computations predict an irregular shock-wave configuration similar to that observed previously in experiments. The existence of a viscous zone in the region of shock-wave interaction allows a continuous transition from the parameters behind the Mach stem to the parameters behind the reflected shock, which is impossible in the inviscid three-shock theory. Results of numerical simulations agree qualitatively with the conclusions of the approximate viscous theory, which describes the flow in the vicinity of the three-shock intersection.     

汽油/空气两相脉冲爆震发动机触发爆震的研究

有效的触发爆震波是脉冲爆震发动机正常工作的关键所在。为了在多循环两相脉冲爆震发动机的缓燃到爆震转变过程中,有效的控制激波反射,成功的触发爆震。通过在爆震室内合理的排放障碍物来促进激波与火焰的相互作用,组织激波反射,适时的触发爆震,从而有效地缩短促发爆震的距离,并最终获得充分发展的爆震波。试验发现爆震管内同样的障碍物在不同的位置起到的作用不同。在直径58 mm,长度1275 mm的汽油/空气PDE爆震管内成功触发频率为50 Hz的爆震。      

爆轰波穿越惰性气团时的透射激波参数分析

爆轰波穿越惰性气团再进入可燃混合气体时,透射激波的参数主要受到三个基本过程的影响:头激波两次穿越惰性气团与可燃气体之间接触面导致的激波强度下降、放热反应中止、爆轰波后稀疏波对激波的衰减作用。放热反应中止是透射激波衰减的主要原因,不同种类的惰性气体对透射激波参数影响不大。     

火焰加速和爆燃向爆震转变过程的数值模拟及试验验证

为了研究火焰加速现象、爆燃向爆震转变和不稳定爆震向稳定爆震的转变过程,对带环形孔板的爆震室进行了数值模拟.研究发现用较低的点火能量对爆震室中的燃料和氧化剂点火产生层流火焰,在孔板的阻碍作用和火焰诱导激波以及反射波的加速作用下,经过几个孔板的阻碍加强作用,在火焰和强激波之间的未燃物中形成爆炸中心,最终引爆未燃混气.同时对爆炸波向稳定爆震转变过程中遇到孔板产生三波点,以及马赫波向入射波转变的全过程进行了分析;对不同燃料在不同当量比下的起爆距离进行了研究,并与试验结果进行了对比.通过分析,对火焰加速和爆燃向爆震转变的过程有了更加全面的认识,为进一步试验提供了参考.