Propagation of shock and detonation waves in dust-laden gases

This article examines the flows of a two-phase mixture of a gas with solid particles arising as a result of the propagation of shock waves or detonation waves through a homogeneous medium at rest. It is assumed that the basic assumptions of the mechanics of mutually penetrating continua hold [1], whereby it is possible to describe the flow of each phase of the mixture within the framework of the mechanics of a continuous medium. We assume that the solid phase consists of identical, incompressible, and nondeformable particles of spherical shape. It is assumed that the temperature inside the particles is homogeneous. Collisions between particles and their Brownian motion are ignored. It is assumed that the carrier phase is an ideal gas (the viscosity is only allowed for in the interaction forces between phases). The contribution of the volume of the particles is not considered. On the basis of these assumptions, the following problems are considered: the propagation of a detonation wave in a mixture of a detonating gas and chemically inert particles and the motion of a dust-gas mixture in a shock tube in the presence of combustion of the particles.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 93–99, November–December, 1984.     

Ionization in gaseous detonation waves

Using numerical simulations the spatial and temporal evolution of the electrical conductivity in ionized gas mixtures produced by detonation is investigated in mixtures with and without potassium seeding. Without flow seeding, the magnitude of the electrical conductivity is observed to be too small for significant magnetohydrodynamic forces to affect the conducting flow. This is consistent with past observations. With potassium seeding the electrical conductivity can be increased by 4–5 orders of magnitude, however, a critical percentage of seed particles is observed for which detonation is no longer sustainable. Therefore, there are limits to the electrical conductivity that can be achieved in burned gas mixtures. To further understand the dynamics of the physics involved a parametric study is conducted by varying the ambient pressure, the nitrogen dilution and the potassium seeding percentage. In order to reach these conclusions, a detailed kinetic mechanism of 26 species and 65 reactions has been compiled from the available literature and validated for applications of detonation and ionization chemistry. From the computed mass fractions the mixture-averaged electrical conductivity is then computed.     

A chemical shock tube driven by detonation

A chemical shock tube driven by a detonation driver is described in the present paper. This shock tube can produce a single controlled high-temperature pulse for studies of gas-phase reaction kinetics, but the difficulty associated with the timing for the rupture of diaphragms in the conventional chemical shock tube is overcome, because the detonation wave in the driver section can be predicted correctly and shows a good repeatability. In addition, this shock tube is capable of providing higher temperature conditions for the test gas than the conventional high-pressure shock tube, owing to the inherently high-driving capability of the detonation driver. The feasibility of this shock tube is examined by numerical simulations and preliminary experiments.     

Re-initiation phenomenon of gaseous detonation induced by shock reflection

The two-dimensional, time-dependent, reactive Navier–Stokes equations including the effects of viscosity, thermal conduction and molecular diffusion were solved to reveal the wave evolution and chemical dynamics involved in the re-initiation process. The computation was performed for hydrogen–oxygen–argon mixtures at the low initial pressure (8.00 kPa), using detailed chemical reaction model. The results show that, the decoupled leading shock reflects on the right wall of the vertical branch. High temperature and pressure behind the reflected shock induce the generation of hot spots and local explosion. Therefore, the re-initiation of gaseous detonation occurs. In the re-initiation area, there exist very high OH concentration and no H ₂ concentration. However, in front of reflected shock, there exist relatively high H ₂ concentration and no OH radicals. Additionally, the shock–flame interaction induces RM instability. This results in the fast mixing between hot reacted gas mixture and the relatively cold unreacted gas mixture and accelerates the chemical reactions. However, the shock–flame interaction contributes much less to the re-initiation, in contrast with shock reflection. The transition of leading shock from regular reflection to Mach reflection happens during the re-initiation. The computed evolution of wave structures involved in the re-initiation is qualitatively agreeable with that from the experimental schlieren images.      

气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明：气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。

     

Performance of a detonation driven shock tunnel

A detonation driven shock tunnel is useful as a ground test facility for hypersonic flow research. By attaching a convergent section ahead of the primary diaphragm in the driver section, the downstream operation mode became available to generate a high-enthalpy test flow. A 100 mm diameter shock tunnel was for the first time installed in the Laboratory of High-Temperature-Gas Dynamics (LHD), Institute of Mechanics, Chinese Academy of Sciences, and after its continuous refitments, a high performance detonation driven shock tunnel was achieved to generate high-enthalpy and high-Reynolds number test flows. A new method to burst a metal diaphragm with the downstream operation mode is discussed.Received: 13 December 2003, Accepted: 26 August 2004, Published online: 26 November 2004[/PUBLISHED]W. Zhao: Correspondence to     

爆轰波对碰区产物驱动金属圆管的研究

利用闪光X光照相、光学高速分幅照相实验,测得了金属圆管在两个爆轰波对碰后产物驱动外壳的膨胀、变形过程,该过程表明爆轰波对碰作用使得对碰处金属圆管相对贯穿断裂时间明显提前,对碰处金属圆管径向膨胀速度较非对碰处明显增加。利用拉格朗日二维流体动力学程序TTD2C,数值模拟了对碰驱动过程,计算结果与实验结果符合很好,且用Taylor断裂判据得到的金属圆管相对贯穿断裂时间也明显提前。     

Oxyhydrogen combustion and detonation driven shock tube

The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incidence shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation. The project supported by State Science and Technology Committee, National Natural Foundation of Science of China (19082012), Chinese Academy of Sciences and Project of National High Technology of China. In memory of academician Kuo Yonghuai's 90th anniversary.     

Propagation of weak shock waves

An asymptotic method is developed in order to treat the evolution of weak shock waves. One obtains a geometrical theory according to which weak shock waves propagate along rays and satisfy a transport law.     

Comparative characteristics of strong shock and detonation waves in bubble media by an electrical wire explosion

Strong shock and detonation waves in inert and chemically active bubble media, which are generated by a wire explosion initiated by a capacitor with a stored energy $W_0 =12.3$ –1,600 J, is experimentally studied. The measurements are performed near the wire and far from the wire in a vertical shock tube 4.5 m long with a volume fraction of the gas in the medium $\beta _0 =1$ –4 %. It is shown that in inert bubble medium, a short intensely decaying shock wave (SW) with intense pressure oscillations is formed in the vicinity of wire explosion point; near the explosion point at $\beta _0 \le 2$  % the SW propagates with the velocity of sound in a liquid. In chemically active bubble medium, an unsteady detonation wave generated by a wire explosion is formed. The pressure amplitude and the velocity of this wave are greater and the length is smaller than those of SW in an inert bubble medium in the same range of explosion energy. It is found that in the interval of low energy explosion from ${\sim }12$ to 64 J, the formation of the bubble detonation wave occurs faster than that at high energies ( $3\times 10^{2}$ – $10^{3}$  J).     

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

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

Propagation of detonation waves in plane channels with obstacles

The detonation wave propagation in plane channels filled with a stoichiometric hydrogen-air mixture at rest under standard conditions is numerically modeled with account for the actual kinetics of the chemical interaction. The calculations show that the stable cellular structure of the detonation wave formed in a plane channel with parallel walls is not always uniquely determined by its width. The effect of transverse walls and sharp expansion of the channel on the propagation of the cellular detonation wave is studied. The conditions of conservation and restoration of detonation are determined.     

激波绕射碰撞加速诱导爆轰的数值模拟

基于单步化学反应的Euler方程和对激波(爆轰波)、接触间断具有良好捕捉效果的Roe/HLL混合格式以及自适应网格技术,模拟了激波在方形管中与方块障碍物相互作用,并发生绕射碰撞来诱导爆轰的过程.结果表明,弱激波在绕经方块时,形成上、下绕射激波并在方块尾部发生碰撞,生成局部高温高压点,可加快爆轰的形成;而当管内阻塞比超过...     

Dynamic characteristics of spherically converging detonation waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure. Zonglin Jiang is presently a visiting professor at McGill University, Canada.     

Advances in the analytical study of the dynamics of gaseous detonation waves

Recent theoretical results on the dynamics of gaseous detonations are presented. An asymptotic analysis is performed, retaining the physical mechanisms controlling the modifications to the inner structure of the detonation. As a result, the system of hyperbolic equations for the compressible fluid mechanics coupled with a detailed chemical kinetics of heat release is reduced to a single integral equation for the propagation velocity of the combustion wave versus time. Concerning the direct initiation of spherical detonations by a blast wave, curvature effects are shown to be responsible for a critical condition of initiation. Near criticality, the role of the unsteadiness of the inner structure is pointed out. The whole complexity of the critical dynamics is reproduced and explained by the integral equation. The necessary background knowledge in gaseous detonation is recalled in the two first sections of the article in order to facilitate the reading by non-specialists.     

Propagation of shock waves in free space

The problem of the exit of a shock wave from an axisymmetric channel and its propagation in a free space occupied by an ideal gas is examined. This problem has been studied earlier in [1], in which the shock wave front was considered planar, as well as in [2], in which the wave front was regarded as a surface of an ellipsoid of revolution. The solutions obtained in these studies assumed the presence of two regions in the wave-front surface: the region of the original shock wave and a region stemming from the decomposition of an infinitesimally thin annular discontinuity of the gas parameters, with the wave intensity over the front surface in each region being considered constant, i.e., the wave character of the process over the front was not considered. In this study a solution will be achieved by the method of characteristics [3–5] of the equations of motion of the shock-wave front, as obtained in [6, 7]. Flow fields are determined for the region immediately adjacent to the shock-wave front for a wide range of shock-wave Mach numbers M _a =1.6–20.0 for = 1.4. On the basis of the data obtained, by introduction of variables connected with the length of the undisturbed zone, as calculated from the channel cross-section along the x axis, together with the pressure transition at the wave front, relationships are proposed which approach self-similarity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 163–166, September–October, 1971.In conclusion, the author thanks S. S. Semenov for his valuable advice on this study.     

Propagation of shock waves in nonuniform plasma

Results of numerical simulation of the propagation of one-dimensional magnetohydrodynamic shock waves in a nonuniform plasma containing a magnetic field are discussed. Possible uses for the production of high velocities and temperatures and astrophysical applications are considered. The essential effect of the magnetic field is shown; acceleration of a shock wave is intensified in a medium with decreasing density.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 22–26, March–April, 1976.In conclusion, the authors are grateful to S. K. Godunov for a detailed discussion of the computational aspects and of the results, and to A. E. Voitenko for a discussion of experimental possibilities and of the results.     

爆轰波对碰驱动下金属圆管膨胀变形特性研究

采用高速分幅摄影和激光速度干涉仪（velocity interferometer system for any reflector, VISAR）联合测试技术,对HR-2钢管在爆轰波对碰驱动下的膨胀变形特性进行了研究。实验得到了圆管膨胀变形的过程图像以及圆管特定部位径向膨胀的位移-时间曲线和速度-时间曲线。采用LSDYNA-3D程序对HR-2钢管的膨胀变形过程进行了数值模拟,计算结果与实验结果吻合较好。     

爆轰波对碰驱动下加环紫铜圆管对碰区变形特性研究

针对爆轰波对碰驱动下金属圆管对碰区易出现过早断裂的问题,以加环紫铜（T2）圆管作为研究对象,利用高速分幅摄影、脉冲X光照相、VISAR三种测试技术,对其在爆轰波对碰驱动下的变形及破坏过程进行了观测。实验结果显示：在紫铜圆管对碰部位增加保护环结构,可以有效推迟该部位圆管的破裂时间。但同时也发现：加保护环后圆管对碰区会出现射流状超前突起现象。由高速摄影结果与X光照相结果对比分析证实：射流状超前突起物质为圆管表面微喷颗粒与空气冲击波混合形成的雾状结构。另外,采用DYNA3D程序对实验模型进行了模拟,计算结果与实验结果符合较好。最后,利用DYNA3D程序就保护环倒角变化对金属圆管对碰区壳体破裂的影响特征进行了初步分析。