Supersonic H2-air combustions behind oblique shock waves

In order to study the mechanisms of initiation and stabilization of H₂-Air combustions (stoechiometric mixture initially atT ₀=293 K andp ₀=0.5 bar) in supersonic flow conditions behind an oblique shock wave (OSW), an original technique is used where OSW is generated in this mixture by the lateral expansion of the burnt gas behind a normal CJ gaseous detonation propagating into a bounding reactive mixture. Four Mach numberM of propagation of OSW are considered in the study, namelyM=7.7-6.1-4.4 and 3. Depending on the Mach numberM and inclinaison angle of OSW different regimes of combustion may occur in the driven mixture. For high values ofM (6.1 and 7.7) delayed steady overdriven oblique detonation waves (SODW) were obtained with a near CJ detonation wave as the critical regime. It was found that SODW obtained correspond quite well to prediction of the polar method. When thermal conditions behind the OSW are lower, either for high Mach number 6.1 and 7.7 for smaller angle than the previous case, or for lower Mach number, 4.4 and 3, the flame initiated at the apex is stabilized as a turbulent oblique flame behind the OSW. With much lower conditions, no combustion appears in the H₂-Air mixture.     

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.     

Influence of mixture sensitivity and pore size on detonation velocities in porous media

Experiments have been carried out to determine the dependence of the detonation velocity in porous media, on mixture sensitivity and pore size. A detonation is established at the top end of a vertical tube and allowed to propagate to the bottom section housing the porous bed, comprised of alumina spheres of equal diameter (1–32 mm). Several of the common detonable fuels were tested at atmospheric initial pressure. Results indicate the existence of a continuous range of velocities with change in Φ, spanning the lean and the rich propagation limits. For all fuels in a given porous bed, the velocity decreases from a maximum value at the most sensitive mixture near Φ≈1 (minimum induction length), toV/V _CJ≈0.3 at the limits. A decrease in pore size brings about a reduction inV/V _CJ and a narrowing of the detonability range for each fuel. For porous media comprised of spherical particles, it was possible to correlate the velocity data corresponding to a variety of different mixtures and for a broad range of particle sizes, using the following empirical expression:V/V _CJ=[1–0.35 log(d _c /d _p)]±0.1. The critical tube diameterd _c is used as a measure of mixture sensitivity andd _p denotes the pore diameter. An examination of the phenomenon at the composition limits, suggests that wave failure is controlled by a turbulent quenching mechanism.     

Behavior of the vorticity vector in supersonic axisymmetric swirl flows behind a detonation wave

The behavior of the vorticity vector on a discontinuity surface arising in a supersonic nonuniform combustible gas flow with the formation of a shock or detonation wave is studied. In the general case, it is a vortex flow with prescribed distributions of parameters. It is demonstrated that the ratio of the tangential component of vorticity to density remains continuous in passing through the discontinuity surface, while the quantities proper become discontinuous. Results calculated for flow vorticity behind a steady-state detonation wave in an axisymmetric supersonic flow of a combustible mixture of gases are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 15–21, November–December, 2007     

Periodic oscillation and fine structure of wedge-induced oblique detonation waves

An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.     

Interaction of a gasdynamic shock with small disturbances

The stability of shock wave based on the definition of Landau and Lifschitz^[1] is treated in this paper. This is tantamount to solving the problem of interaction of small disturbances with a shock wave. Small disturbances are introduced on both sides of a steady, non-dissipative, plane shock wave. Landau et al.^[1] obtained the stability criterionM ₁>1,M ₂<1 for small disturbances which are travelling in the direction perpendicular to the shock wave. In the present paper, we assume that the small disturbances may be two dimensional, i.e. they may be propagating in the direction inclined to the shock wave. The conclusions obtained are: regardless of whether the incident wave and diverging wave are defined according to the direction of the phase velocity or the group velocity, the shock wave is unstable for some frequencies and longitudinal wave lengths of the disturbances, even if the conditionsM ₁>1,M ₂<1 are fulfilled. Then several experiments are proposed, and the problem of ways to define the incident wave and diverging wave is discussed. The meaning of this problem is illustrated. The same results can be obtained for the steady shock wave in a tube.     

Temperature measurement of reflected shock wave by using chemical indicator

This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction of small amount of CF₄ premixed in the test gas. The final product C₂F₄ is used as the temperature indicator, which is sampled and detected by a gas chromatography in the experiment. The detected concentration of C₂F₄ is correlated to the temperature of the reflected shock wave with the initial pressureP ₁ and test time τ as parameters in the temperature range 3 300 K<T<5 600 K, pressure range 5 kPa<P ₁<12 kPa and τ≅0.4 ms. The project supported by the China Aerodynamics Project for Basic Researches (J13.5.2 ZK04)     

有限谱ENO格式在爆轰波数值模拟中的应用

使用有限谱ENO(EssentiallyNon Oscillatory)格式 ,采用V .P .Korobeinikov二步化学反应模型 ,对稀释的化学当量的氢气和氧气混合物中非定常自维持爆轰波进行了二维数值模拟 ,研究了波的产生和演化机理 ,分析了爆轰波的三波结构和传播过程。计算得到的爆轰波参数和结构与以前的计算和实验结果一致。研究表明 :有限谱ENO格式可以成功地模拟非定常自维持爆轰波。     

Mach reflection of a plane shock wave passing a mountain of 45° inclination

The transition from regular reflection (RR) to Mach reflection (MR) as a plane shock wave diffracts around a triangular mountain of 45° inclination is analysed in this paper, both by optical measurement in a shock tube and by numerical simulation the numerical method developed by Li Yingfan^[1] is of the FLIC type with triangular mesh. The dependence of the critical transition point Lk ofRR→MR on shock Mach numberM _i is analyzed and the variations of the incidence angle ω _i of the impinging shock and the reflection angle ω _r with the distanceL ^* are investigated. Our experimental and numerical results agree well with the theoretical results of Iton and Italya.     

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters he, v, w, L, wb. The φ800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2-35)GHz (ω=2πf, wave length A =15cm-8mm). The electron density in the plasma is ne = (3&#215;10^10-1&#215;10^14) cm^-3. The collision frequency v = (1&#215;10^8-6&#215;10^10) Hz. The thickness of the plasma layer L = (2-80)cm. The electron circular frequency ωb=eBo/me, magnetic flux density B0 = (0-0.84)T. The experimental results show that when the plasma layer is thick (such as L/λ≥10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters ne,v,ω, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and A are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters ne,v,ω, L. In fact, if ω&lt;ωp, v^2&lt;&lt;ω^2, the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if ω&gt;ωp, v^2&lt;&lt;ω^2 (just v≈f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.     

Wave angle for oblique detonation waves

The flow field associated with a steady, planar, oblique detonation wave is discussed. A revision is provided for- diagrams, where is the wave angle and is the ramp angle. A new solution is proposed for weak underdriven detonation waves that does not violate the second law. A Taylor wave, encountered in unsteady detonation waves, is required. Uniqueness and hysteresis effects are also discussed.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Shock wave triple-point morphology

The existence and characteristics of shock wave triple points are examined. The analysis utilizes a single flow plane for the three shocks and is local to the triple point. It applies when the flow is unsteady, three-dimensional, and the upstream flow is nonuniform. Under more restrictive conditions, a relation is also derived for the ratio of the curvature of the Mach stem to that of the reflected shock. For given values of the ratio of specific heats, γ, and the upstream Mach number, M ₁, a solution window is established. A parametric set of solutions is generated within the window for γ = 1, 1.4, and 5/3 and for 16 values of M ₁ ranging from solution onset to M ₁ = 6.A solution can be one of three types, these stem from the velocity tangency condition along the slip stream. Topics are addressed such as solution multiplicity, shock wave and slip stream orientation, the nature of the reflected wave (weak, strong, inverted, normal), the nature of the Mach stem (weak, strong, normal), and differences due to changes in γ and M ₁.     

Detonation wave initiated by explosive condensation of supersaturated carbon vapor

An experimental study of the influence of condensation of supersaturated carbon vapor formed behind reflected shock waves on the process of propagation of a shock wave and formation of a detonation wave of condensation is carried out. Highly supersaturated carbon vapor was formed from thermal decay of unstable carbon suboxide C₃O₂ → C + 2CO behind a shock wave in mixtures containing 10–30% C₃O₂ in Ar. This reaction was followed by fast growth of condensed carbon particles, accompanied by heat release. Experiments have shown a considerable temperature and pressure increase in the narrow zone behind the wave front, resulting in shock wave amplification and transition to a detonation-like regime. An analysis of the kinetics and heat release in the given conditions and calculations based upon one-dimensional detonation theory have shown that in a mixture of 10% C₃O₂ + Ar, insufficient heat release resulted in a regime of “overdriven detonation”. In a mixture of 20% C₃O₂ + Ar a very good coincidence of measured values of pressure and wave velocity with calculated Chapman–Jouguet parameters is observed. In a 30% C₃O₂ + Ar mixture, an excess heat release caused a slow down of the effective condensation rate and a regime of “underdriven detonation” is observed.     

Detonability of THDCPD-exo–air mixtures

In the frame of industrial risk and propulsive application, the detonability study of JP10–air mixtures was performed. The simulation and measurements of detonation parameters were performed for THDCPD-exo/air mixtures at various initial pressure (1 bar < P ₀ < 3 bar) and equivalence ratio (0.8 < Φ < 1.6) in a heated tube (T ₀ ~ 375 K). Numerical simulations of the detonation were performed with the STANJAN code and a detailed kinetic scheme of the combustion of THDCPD. The experimental study deals with the measurements of detonation velocity and cell size λ. The measured velocity is in a good agreement with the calculated theoretical values. The cell size measurements show a minimum value for Φ ~ 1.2 at every level of initial pressure studied and the calculated induction length L _i corresponds to cell size value with a coefficient k = λ/L _i = 24 at P ₀ = 1 bar. Based on the comparison between the results obtained during this study and those available in the literature on the critical initiation energy E _c, critical tube diameter d _c and deflagration to detonation transition length L _DDT, we can conclude that the detonability of THDCPD–air mixtures corresponds to that of hydrocarbon–air mixtures.

---

This paper is based on the work presented at the 33rd International Pyrotechnics Seminar, IPS 2006, Fort Collins, July 16–21, 2006.     

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.     

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟。研究结果表明，标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波，爆轰波与激波和接触间断相互作用产生了复杂的波系结构；爆轰波爆点位置在对称轴上并不是固定的点，而是随着初始激波马赫数的变化而发生移动；可燃气体初始温度和压力对起爆临界马赫数都有影响，但是初始温度的影响大得多。     

爆炸去磁脉冲功率发生器的实验和理论计算

设计了两种结构的爆炸去磁脉冲功率发生器:炸药透镜平面波冲击加载圆柱形磁体和中心装药冲击加载圆环形磁体，磁体材料都采用Nd2Fe14B永磁体。进行了爆炸冲击加载实验，测量了爆炸去磁脉冲发生器的输出电流、电压以及Nd2Fe14B永磁体中的冲击波速度。通过对爆炸去磁脉冲功率发生器基本原理的分析，建立了感生电动势的理论计算方法。用Maxwell 3D电磁场分析软件，对Nd2Fe14B永磁体进行了静磁场分析，给出了Nd2Fe14B永磁体的磁场强度分布规律。实验结果表明:爆炸去磁脉冲功率发生器每匝线圈能够产生约67 V的感生电动势,对于0.56 H的电感负载能够输出脉宽78.3 s、峰值1 212 A的脉冲电流。     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Precursor vacuum-ultraviolet radiation ahead of strong shock waves in a shock tube

The profile and excitation mechanism of vacuum-ultraviolet radiation emitted from shock wave is investigated in a shock tube. For shock wave in argon, the rdiation is due to resonant transition excited by argon-argon collision in the shock front with excitation cross section coefficientS ^*=1.0×10⁻¹⁷ cm²·ev⁻¹ and activation energyE ^*=11.4 ev. For shock wave in air the radition is emitted from a very thin shock layer in which the mechanism ofX ¹∑→b ¹∑ of N₂ is excited with excitation cross sectionQ=2×10⁻¹⁶cm² and activation energyE ^*=12.1 ev. Institute of Mechanics, Academia Sinica     

Vorticity and its rate of change just downstream of a curved shock

A theory is developed for the vorticity and its substantial derivative just downstream of a curved shock wave, the resulting formulas are exact, algebraic, and explicit. Analysis is for a cylinder-wedge or sphere-cone body, at zero incidence, whose downstream half-angle is θ_b. Derived formulas directly depend only on the ratio of specific heats, γ, the freestream Mach number, M ₁, the local slope and curvature of the shock, and the dimensionality parameter, σ, which is zero for a two-dimensional shock and unity for an axisymmetric shock. In turn, the slope and curvature depend on γ, M ₁, and θ_b. Numerical results are provided for a bow shock in which θ_b is 5°, 10°, or 15°, M ₁ is 2, 4, or 6, and γ = 1.4. There is little dependence on the half angle but a strong dependence on the freestream Mach number and on dimensionality. For vorticity and its substantial derivative, the dimensionality dependence gradually decreases with increasing Mach number. In comparison to the two-dimensional case, an axisymmetric shock generates considerable vorticity in a region relatively close to the symmetry axis. Moreover, the magnitude of the vorticity, in this region, is further enhanced in the flow downstream of the shock. This dimensionality difference in vorticity and its substantial derivative is attributed to the three-dimensional relief effect in an axisymmetric flow.

---