甲烷-空气爆轰波传播的数值模拟

研究了甲烷-空气混合物爆轰波在二维突扩通道中传播的各种复杂行为。结果表明：爆轰波在进入突扩通道的初始阶段有局部向爆燃转变的现象；爆炸波在壁面发生马赫反射形成的高温高压区域将直接驱动自持爆轰波的重新形成。     

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。     

凝聚炸药中爆轰波马赫反射的实验研究

在一平板装药中,使两平面爆轰波以不同的角度相互碰撞。在一定条件下,产生非正常反射,形成超压马赫爆轰波。本文介绍了可以直接观测马赫爆轰波的形成、生长规律及其传播过程的实验方法,并对马赫爆轰波参数进行了计算。     

气相爆轰波传播过程中的自点火效应

基于基元反应模型和单步反应模型,对直管道中H₂-air混合气体中爆轰波的传播过程进行了数值模拟,揭示了气相爆轰波传播过程中的自点火效应。利用数值模拟方法计算了不同爆轰模型的点火延迟时间,并得到了爆轰波三波点的传播过程以及所形成胞格结构的尺寸。结果表明,胞格宽度与点火延迟时间成正比;爆轰波诱导区内气体的点火延迟时间与三波点的运动周期基本一致。进一步对结果分析可知,爆轰波的自维持传播取决于点火延迟时间(表征化学反应的特征时间)和三波点的运动周期(表征流动的特征时间)的匹配;当二者相匹配时,经过前导激波压缩后形成的高温高压爆轰气体,在短时间内实现了自点火,同时释放出大量的能量推动了爆轰波的前进,即爆轰波的稳定自维持传播依靠其自点火机制。     

圆柱形爆轰波的二维数值模拟

基于带化学反应的二维Euler方程,对圆柱形爆轰波的直接起爆和传播过程进行了二维数值模拟研究,拟分析起爆条件和初始压强对圆柱形爆轰波形成和传播的影响。研究发现,圆柱形爆轰波起爆成功向外传播的过程中,新的三波结构的生成标志着爆轰波进入稳定传播阶段。在起爆能量足够的情况下,起爆半径(曲率)的大小决定着三波结构初始形成时的数目和传播半径,起爆压强对其基本不产生影响;起爆半径大(曲率小)时,三波结构初始形成时的传播半径大、数目多,圆柱形爆轰波进入稳定传播阶段的传播距离长;数值模拟中,初始压强的提高,有助于圆柱形爆轰在较短的传播距离内进入稳定传播阶段。     

气相爆轰波在障碍物上Mach反射的实验验证

本文公布了气相爆轰波沿收缩管道传播时发生Mach反射的实验证据。在爆轰波通过的管道中安装不同楔角的楔块,形成管道的收缩。爆轰波在通过楔块时会发生Mach反射。利用烟熏玻璃片记录到了爆轰波Mach反射时形成的三波点迹线及其两侧胞格尺寸和密度的变化。据我们掌握的资料,这是首次用胞格结构变化的记录证实,气相爆轰波与无化学反应的空气中的冲击波一样,在一定的入射条件下会发生Mach反射。这一实验结果可使我们更深入了解爆轰波的本质,也为数值模拟气相爆轰波在障碍物上Mach反射现象提供了可对比的依据。     

气体爆轰波在声学吸收壁下游的再加强过程

实验研究了爆轰波在声学吸收段下游发生的再恢复过程。利用声学吸收段首先将一个稳定的爆轰波衰减,于是产生一个非稳定爆轰波,速度降为 45% ~60?J。该非稳定爆轰波传出声学吸收段后再次加强,并在距声学吸收段出口端一定距离处突然发展成超驱动爆轰。实验利用烟迹技术捕捉了爆轰波从发生衰减到再次加强的全过程,获得了爆轰波再加强距离随初压变化的实验规律。由此,本文建立了再加强段内加速传播火焰的理论模型,称为激波 -爆燃波复合态,并解释了火焰加速和超驱动爆轰形成的机理。     

爆轰波平掠惰性气体界面及其解耦现象的数值研究

采用基元化学反应模型和迎风TVD格式,数值研究了爆轰波平掠惰性气体界面时的物理现象及其作用机制,并用点隐方法克服化学反应源项引起的计算刚性。数值结果显示,当爆轰波平掠过惰性气体界面时,形成了爆轰波、界面、透射激波以及稀疏波相互作用的现象。在高N2比例稀释的可燃混合气体情况下,当爆轰波平掠过特定惰性气体界面时,它与惰性气体界面相互作用产生的稀疏波可以导致爆轰波的解耦。     

多障碍物通道中激波诱导气相爆轰的数值研究

应用数值模拟方法，研究了直通道中激波经过多块矩形障碍物时诱导 H2/O2混合气体起爆的物理机制. 研究表明：在前导激波强度不足以诱导波后气 体直接起爆的条件下，经过激波压缩的可燃气体也可能在远离激波的障碍物之间的凹槽部位 起爆；障碍物表面产生的压缩波、膨胀波和气流滑移面对可燃气体的起爆、爆轰波的形成和 传播过程有重要的影响；添加不同稀释比的氮气可以影响爆轰波后流场的温度分布；增加障 碍物的间距可以改变可燃气起爆位置.     

基于CE/SE方法模拟空气中RDX-Al悬浮粉尘的两相爆轰

通过CE/SE方法模拟了空气中炸药-铝粉尘的两相爆轰过程, 研究了双粉尘爆轰过程中粒子不同密度对爆轰波速度、压力的影响, 得到密度与波速、爆压间的线性关系。模拟得到悬浮粉尘在复杂通道中的爆轰波传播过程, 研究了双粉尘爆轰的流场演化过程, 选取流场中的一些点对该处流场的压力及温度随时间的变化进行重点研究, 对比了单铝粉尘在同种条件下的爆轰过程, 发现双粉尘爆轰明显提高了爆轰波波速和流场的压力及温度。模拟结果表明CE/SE方法可以成功模拟双粉尘的爆轰过程, 可为多粉尘爆轰的研究提供参考。     

带隔板装药爆轰波马赫反射理论研究和数值模拟

基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。     

Numerical simulation of Mach reflection of cellular detonations

The Mach reflection of cellular detonation waves on a wedge is investigated numerically in an attempt to elucidate the effect of cellular instabilities on Mach reflection, the dependence of self-similarity on the thickness of a detonation wave, and the initial development of the Mach stem near the wedge apex. A two-step chain-branching reaction model is used to give a thermally neutral induction zone followed by a chemical reaction zone for the detonation wave. A sufficiently large distance of travel of the Mach stem is computed to observe the asymptotic behavior in the far field. Depending on the scale at which the Mach reflection process occurs, it is found that the Mach reflection of a cellular detonation behaves essentially in the same way as a planar ZND detonation wave. The cellular instabilities, however, cause the triple-point trajectory to fluctuate. The fluctuations are due to interactions of the triple point of the Mach stem with the transverse waves of cellular instabilities. In the vicinity of the wedge apex, the Mach reflection is found to be self-similar and corresponds to that of a shock wave of the same strength, since the Mach stem is highly overdriven initially. In the far field, the triple-point trajectory approaches a straight line, indicating that the Mach reflection becomes self-similar asymptotically. The distance of the approach to self-similarity is found to decrease rapidly with decreasing thickness of the detonation front.     

Wave dynamic processes in cellular detonation reflection from wedges

When the cell width of the incident detonation wave (IDW) is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H ₂/O ₂ diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes. The English text was polished by Yumming Chen.     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体，在初始压强20kPa的条件下， 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动 力学过程.     

A numerical simulation of reflection processes of a detonation wave on a wedge

Abstract. A two dimensional numerical simulation has been performed to study reflection processes of detonation waves on a wedge. The numerical scheme adopted is the flux corrected transport scheme and a two-step chemical reaction is assumed for a stoichiometric oxyhydrogen mixture diluted with argon. Transverse wave structures of the detonation are produced by artificial disturbances situated in front of a one-dimensional Chapman-Jouguet detonation wave. Numerical grids are generated by solving a Laplace equation. Results show that in the case where Mach reflection occurs, the cells in the Mach stem are smaller than those in the incident wave and are distorted in shape. There is also an initiating stage during which the cells in the Mach stem are created. The critical angle beyond which Mach reflection cannot occur is discussed. Received 15 October 1999 / Accepted 27 March 2000     

Mechanisms of detonation transmission in layered H2-O2 mixtures

When a plane detonation propagating through an explosive comes into contact with a bounding explosive, different types of diffraction patterns, which may result in the transmission of a detonation into the bounding mixture, are observed. The nature of these diffraction patterns and the mode of detonation transmission depend on the properties of the primary and bounding explosives. An experimental and analytical study of such diffractions, which are fundamental to many explosive applications, has been conducted in a two channel shock tube, using H₂-O₂ mixtures of different equivalence ratios as the primary and bounding or secondary explosive. The combination of mixtures was varied from rich primary / lean secondary to lean primary / rich secondary since the nature of the diffraction was found to depend on whether the Chapman-Jouguet velocity of the primary mixture,D _p, was greater than or less than that of the secondary mixture,D _s. Schlieren framing photographs of the different diffraction patterns were obtained and used to measure shock and oblique detonation wave angles and velocities for the different diffraction patterns, and these were compared with the results of a steady-state shock-polar solution of the diffraction problem. Two basic types of diffraction and modes of detonation reinitiation were observed. WhenD _p>D _s, an oblique shock connecting the primary detonation to an oblique detonation in the secondary mixture was observed. WithD _p<D _s, two modes of reinitiation were observed. In some cases, ignition occurs behind the Mach reflection of the shock wave, which is transmitted into the secondary mixture when the primary detonation first comes into contact with it, from the walls of the shock tube. In other cases, a detonation is initiated in the secondary mixture when the reflected shock crosses the contact surface behind the incident detonation. These observed modes of Mach stem and contact surface ignition have also been observed in numerical simulations of layered detonation interactions, as has the combined oblique-shock oblique-detonation configuration whenD _p>D _s. WhenD _p>D _s, the primary wave acts like a wedge moving into the secondary mixture with velocityD _p after steady state has been reached, a configuration which also arises in oblique-detonation ramjets and hypervelocity drivers.     

激波在不同三角楔面内传播特性的数值研究

激波在收缩管内的反射与聚焦会形成高温高压区,点燃可燃混合气并诱导爆轰,因此对爆轰发动机的点火具有重要意义。本文基于二维N-S方程,结合五阶WENO格式,对马赫数为6的正激波在三角形楔面内的反射与聚焦现象进行了数值研究。结果表明,楔面顶角的变化对激波的反射类型以及聚焦均有明显的影响:随着顶角的增加,激波的反射类型从马赫反射向过渡马赫反射和双马赫反射转变,且壁面上的前向射流更加明显;三波点第一次碰撞产生的高温高压区足够满足可燃混合气体的点火条件,且其温度与压力值随顶角的增加而增大;当激波在楔面上发生临界双马赫反射时,温度与压力达到最大;当顶角增加到一定值时,激波在楔面反射转变为常规反射,不会产生激波对碰,因而没有高温高压区。     

气相爆轰波正反射激波加速研究

实验采用压力传感器测量了指定点压力时间曲线。数值模拟基于二维反应欧拉方程和基元反应模型,采用二阶附加半隐的龙格-库塔法和5阶WENO格式分别离散时间和空间导数项,获得了指定点数值压力时间曲线。理论分析基于爆轰理论和激波动力学,分析了气相爆轰波反射过程所涉及的复杂波系演变并获得了反射激波速度。结果表明:本文数值模拟和理论计算定性上重复并解释了实验现象。气相爆轰波在右壁面反射后,右行稀疏波加速反射激波。其加速原因是:尽管激波波前声速减小,但激波马赫数增大,波前气流速度减小。在低初压下,可能还由于爆轰波后未反应或部分反应气体的作用,导致反射激波加速幅度比高初压下大。