Cell-like structure of unstable oblique detonation wave from high-resolution numerical simulation

A comprehensive numerical study was carried out to investigate the unsteady cell-like structures of oblique detonation waves (ODWs) for a fixed Mach 7 inlet flow over a wedge of 30° turning angle. The effects of grid resolution and activation energy were examined systematically at a dimensionless heat addition of 10. The ODW front remains stable for a low activation energy regardless of grid resolution, but becomes unstable for a high activation energy featuring a cell-like wave front structure. Similar to the situation with an ordinary normal detonation wave (NDW), a continuous increase in the activation energy eventually causes the wave-front oscillation to transit from a regular to an irregular pattern. The wave structure of an unstable ODW, however, differs considerably from that of a NDW. Under the present flow condition, triple points and transverse waves propagate downstream, and the numerical smoke-foil record exhibits traces of triple points that rarely intersect with each other. Several instability-driving mechanisms were conjectured from the highly refined results. Since the reaction front behind a shock wave can be easily destabilized by disturbance inherent in the flowfield, the ODW front becomes unstable and displays cell-like structures due to the local pressure oscillations and/or the reflected shock waves originating from the triple points. The combined effects of various instability sources give rise to a highly unstable and complex flow structure behind an unstable ODW front.     

Combustion wave propagation and detonation initiation in the vicinity of closed-tube end walls

There are not many studies on DDT with no obstacles and the initiation of DDT near the end of a closed tube. Therefore in the present study we experimentally investigate the mechanism of the combustion wave transition to a detonation wave when there are no obstacles. In particular, we show that a local explosion near the tube wall is necessary for the initiation of a detonation. Parameters that we varied are the wall configuration, distance between the ignition point and the wall, and initial filling pressure. The combustion waves and the compression waves are visualized using the Schlieren optical system. From the results, we found it is necessary for the combustion wave to reach four walls so that the detonation could be initiated by the local explosion. In the conditions of the present experiment, we exhibited that the local explosion did not occur in the vicinity of a single wall and four orthogonal walls; instead, the local explosion occurred in a situation with five orthogonal walls. The time of the local explosion and the detonation initiation is 2.6 ± 1.1 and 2.0 ± 0.1 times the characteristic time for the combustion wave to propagate hemispherically from an ignitor and reach the four walls.     

Unstable combustion induced by oblique shock waves at the non-attaching condition of the oblique detonation wave

The instability of oblique shock wave (OSW) induced combustion is examined for a wedge with a flow turning angle greater than the maximum attach angle of the oblique detonation wave (ODW), where archival results rarely exist for this case in previous literatures. Numerical simulations were carried out for wedges of different length scales to account for the ratio of the chemical and fluid dynamic time scales. The results reveal three different regimes of combustion. (1) No ignition or decoupled combustion was observed if a fluid dynamic time is shorter than a chemical time behind an OSW. (2) Oscillatory combustion was observed behind an OSW if a fluid dynamic time is longer than a chemical time behind an OSW and the fluid dynamic time is shorter than the chemical time behind a normal shock wave (NSW) at the same Mach number. (3) Detached bow shock-induced combustion (or detached overdriven detonation wave) was observed if a fluid dynamic time is longer than a chemical time behind a NSW. Since no ignition or decoupled combustion occurs as a very slow reaction and the detached wave occurs as an infinitely fast reaction, the finite rate chemistry is considered to be the key for the oscillating combustion induced by an OSW over a wedge of a finite length with a flow turning angle greater than the maximum attach angle for an ODW. Since this case has not been previously reported, grid independency was tested intensively to account for the interaction between the shock and reaction waves and to determine the critical time scale where the oscillating combustion can be observed.     

Shock wave and detonation propagation through U-bend tubes

The objective of the research outlined in this paper is to provide experimental and computational data on initiation, propagation, and stability of gaseous fuel–air detonations in tubes with U-bends implying their use for design optimization of pulse detonation engines (PDEs). The experimental results with the U-bends of two curvatures indicate that, on the one hand, the U-bend of the tube promotes the shock-induced detonation initiation. On the other hand, the detonation wave propagating through the U-bend is subjected to complete decay or temporary attenuation followed by the complete recovery in the straight tube section downstream from the U-bend. Numerical simulation of the process reveals some salient features of transient phenomena in U-tubes.     

Initiation of detonation by conical projectiles

Initiation of detonation by a hypersonic conical projectile launched into a combustible gas mixture is investigated. From analytic considerations of the flowfield, energetic and kinetic limits are proposed to predict the conditions required to initiate an oblique detonation wave in the mixture. To experimentally investigate these limits, projectiles with cone half angles varying from 15° to 60° were launched into a stoichiometric mixture of hydrogen/oxygen with 70% argon dilution at initial pressure between 10 and 200 kPa. The projectiles were launched from a combustion-driven gas gun at velocities as great as 2.5 km/s (corresponding to 150% of the Chapman Jouguet velocity). Pictures of the flowfields generated by the projectiles were taken via schlieren photography. Five combustion regimes could be observed about the projectile ranging from a prompt and delayed oblique detonation wave formation, combustion instabilities, a wave splitting, and an inert shock wave. The two theoretical limits provide a means to interpret the observed flowfield regimes and are in satisfactory agreement with the experimental results.     

Calibration of the Pseudo-Reaction-Zone model for detonation wave propagation

An approach for the calibration of an advanced programmed burn (PB) model for detonation performance calculations in high explosive systems is detailed. Programmed burn methods split the detonation performance calculation into two components: timing and energy release. For the timing, the PB model uses a Detonation Shock Dynamics (DSD) surface propagation model, where the normal surface speed is a function of local surface curvature. For the energy release calculation and subsequent hydrodynamic flow evolution, a Pseudo-Reaction-Zone (PRZ) model is used. The PRZ model is similar to a reactive burn model in that it converts reactants into products at a finite rate, but it has a reaction rate dependent on the normal surface speed derived from the DSD calculation. The PRZ reaction rate parameters must be calibrated in such a way that the rate of energy release due to reaction in multi-dimensional geometries is consistent with the timing calculation provided by the DSD model. Our strategy for achieving this is to run the PRZ model in a detonation shock-attached frame in a compliant 2D planar slab geometry in an equivalent way to a reactive burn model, from which we can generate detonation front shapes and detonation phase speed variations with slab thickness. In this case, the D _n field used by the PRZ model is then simply the normal detonation shock speed rather than the DSD surface normal speed. The PRZ rate parameters are then iterated on to match the equivalent surface front shapes and surface phase speed variations with slab thickness derived from the target DSD model. For the purposes of this paper, the target DSD model is fitted to the performance properties of an idealised condensed-phase reactive burn model, which allows us to compare the detonation structure of the calibrated PRZ model to that of the originating idealised-condensed phase model.     

多管PDE爆轰噪声传播过程物理特性实验研究

为研究多管脉冲爆轰发动机爆轰噪声传播过程及声波物理特性,对单管至四管脉冲爆轰发动机爆轰噪声物理特性开展实验研究,获得了多管脉冲爆轰发动机爆轰噪声的波形、声压衰减规律、辐射特性、持续时间和频谱特性等物理参数。结果表明:在管口区域,爆轰噪声峰值衰减较快,在远离管口区域,衰减速率逐渐放缓。随着爆轰管数量的增加,爆轰峰值噪声在各方位角上的声压级逐渐增大。在远场区域内爆轰噪声指向性图呈“M”形。随着爆轰管数量的增加,爆轰噪声的A/B持续时间均增加。爆轰噪声是能量主要集中在中低频段的宽频噪声,其频谱峰值和峰值所在频点随爆轰管数量的增加而改变。     

Effect of debris fragments on direct initiation of spherical detonation waves in stoichiometric oxygen/hydrogen mixtures

Paper reports a result of experiments of spherical shock waves generated by explosions of micro-explosives weighing from 1 to 10 mg ignited by the irradiation of Q-switched laser beam and direct initiation to a spherical detonation wave in stoichiometric oxygen/hydrogen mixtures at 10–200 kPa. We visualized the interaction of debris particles ejected micro-explosives’ surface with shock waves by using double exposure holographic interferometry and high-speed video recording. Upon explosion, minute inert debris launched supersonically from micro-charge surface precursory to shock waves initiated spherical detonation waves. To examine this effect we attached 0.5–2.0 μm diameter SiO₂ particles densely on micro-explosive surfaces and observed that the supersonic particles, significantly promoted the direct initiation of spherical detonation waves. The domain and boundary of detonation wave initiations were experimentally obtained at various initial pressures and the amount of micro-charges.     

The quantum effects of the spin and the Bohm potential in the oblique propagation of magnetosonic waves

We study the quantum corrections to the oblique propagation of the magnetosonic waves in a warm quantum magnetoplasma composed by mobile ions and electrons. We use a fluid formalism to include quantum corrections due to the Bohm potential and to the spin magnetization energy of electrons. The effects of both quantum corrections are shown in the dispersion relation for perpendicular, parallel and oblique propagation. We find that the quantum contributions to the low frequency depend on the type in the oblique propagation with respect to the background magnetic field. The relevance in astrophysical scenarios is exemplified.     

Experimental study of the influence of annular nozzle on acoustic characteristics of detonation sound wave generated by pulse detonation engine

Acoustic characteristics of the detonation sound wave generated by a pulse detonation engine with an annular nozzle, including peak sound pressure, directivity, and A duration, are experimentally investigated while utilizing gasoline as fuel and oxygen-enriched air as oxidizer. Three annular nozzle geometries are evaluated by varying the ratio of inner cone diameter to detonation tube exit diameter from 0.36 to 0.68. The experimental results show that the annular nozzles have a significant effect on the acoustic characteristics of the detonation sound wave. The annular nozzles can amplify the peak sound pressure of the detonation sound wave at 90° while reducing it at 0° and 30°. The directivity angle of the detonation sound wave is changed by annular nozzles from 30° to 90°. The A duration of the detonation sound wave at 90° is also increased by the annular nozzles. These changes indicate that the annular nozzles have an important influence on the acoustic energy distribution of the detonation sound wave, which amplify the acoustic energy in a direction perpendicular to the tube axis and weaken it along the direction of the tube axis.     

爆轰波在炸药-金属界面上的折射分析

针对爆轰波在炸药-金属界面上折射时由实验获得的金属折射冲击波压力与经典爆轰波极曲线理论预测的压力存在显著差异这一问题, 本文展开了进一步的理论和数值模拟分析研究. 首先通过分析指出经典爆轰波极曲线理论的缺陷, 并对爆轰波极曲线理论进行了改进, 改进爆轰波极曲线理论给出了炸药爆轰波折射类型以及折射冲击作用点处的压力值. 然后发展了一个基于次特征理论来数值求解爆轰反应流动控制方程的二阶中心型Lagrange方法, 并数值模拟了一个典型的炸药爆轰波折射实验. 改进爆轰波极曲线理论和数值模拟分析结果表明, 爆轰波折射类型有三种:反射冲击波的正规折射、带Mach反射的非正规折射、无反射波的正规折射, 并且金属折射冲击波压力值随入射角增大而单调减小.     

An experimental investigation of the propagation mechanism of critical deflagration waves that lead to the onset of detonation

The reflection of a CJ detonation from a perforated plate is used to generate high speed deflagrations downstream in order to investigate the critical conditions that lead to the onset of detonation. Different perforated plates were used to control the turbulence in the downstream deflagration waves. Streak Schlieren photography, ionization probes and pressure transducers are used to monitor the flow field and the transition to detonation. Stoichiometric mixtures of acetylene–oxygen and propane–oxygen were tested at low initial pressures. In some cases, acetylene–oxygen was diluted with 80% argon in order to render the mixture more “stable” (i.e., more regular detonation cell structure). The results show that prior to successful detonation initiation, a deflagration is formed that propagates at about half the CJ detonation velocity of the mixture. This “critical” deflagration (which propagates at a relatively constant velocity for a certain duration prior to the onset of detonation) is comprised of a leading shock wave followed by an extended turbulent reaction zone. The critical deflagration speed is not dependent on the turbulence characteristics of the perforated plate but rather on the energetics of the mixture like a CJ detonation (i.e., the deflagration front is driven by the expansion of the combustion products). Hence, the critical deflagration is identified as a CJ deflagration. The high intensity turbulence that is required to sustain its propagation is maintained via chemical instabilities in the reaction zone due to the coupling of pressure fluctuations with the energy release. Therefore, in “unstable” mixtures, critical deflagrations can be supported for long durations, whereas in “stable” mixtures, deflagrations decay as the initial plate generated turbulence decays. The eventual onset of detonation is postulated to be a result of the amplification of pressure waves (i.e., turbulence) that leads to the formation of local explosion centers via the SWACER mechanism during the pre-detonation period.     

冲击波物理与爆轰物理研究进展

在中国工程物理研究院(以下简称中物院)建院50周年之际,文章作者对近年来在中物院开展的冲击波物理与爆轰物理研究工作进展情况作了介绍.首先,对该项研究工作的意义作了说明.其次,简单回顾了过去完成的几个代表性研究工作.然后,对近期主要的研究方向,其中包括冲击波温度测量和高压熔化规律研究、高密度气体物态方程研究、固体物质的宽区物态方程实验和理论研究、高压本构关系研究、材料的损伤演化及动态破坏研究、爆轰波传播规律与驱动飞层研究,以及高能炸药起爆性能与热点形成机理研究在实验技术、规律性探索到理论建模方面取得的近期主要进展作了重点介绍.文章最后对未来研究作了展望.     

Acoustic backscattering enhancements resulting from the interaction of an obliquely incident plane wave with an infinite cylinder

Background and objective

The analysis of the acoustic backscattering enhancements from tilted cylinders is of particular importance in determining some of the (visco)elastic properties of the cylinder, and/or its surrounding fluid in ultrasonic non-destructive evaluation (NDE) and imaging (NDI) applications. Previous related investigations on an aluminum cylinder limited to incidence angles varying from 0° to 40°, revealed the existence of an anomalous “pseudo-Rayleigh” mode (above the critical Rayleigh angle) identified as the rigid-body translational dipole (n = 1) mode. The objective here is to provide a complete investigation on the backscattering enhancements for incidence angles larger than 40° for various elastic and viscoelastic cylinder materials.

Method

Using the partial-wave series solution for the linear scattering by an infinite circular cylinder, the acoustic backscattering from isotropic elastic and viscoelastic (polymer-type) cylinders excited by an obliquely incident plane acoustic wave is investigated. Total and resonance backscattering form functions are calculated for several elastic and viscoelastic cylinder materials immersed in water versus the angle of incidence 0° ? α < 90°. The “pure” resonance peaks are isolated by subtracting a rigid background from the total form function, so the associated resonance modes are properly identified.

Results and conclusion

The plots of the partial-wave series reveal acoustic backscattering enhancements (not shown in previous investigations) generally occurring at ka? 0.1 at a critical angle α_c bounded by the longitudinal and shear waves coupling angles θ_L=sin^-1(c/c_L) and θ_S=sin^-1(c/c_S) such that θ_L<α_c<θ_S (where c_L and c_S are the phase velocities of the longitudinal and shear waves inside the elastic cylinder, and c is the speed of sound in the surrounding medium). It is shown here that the backscattering enhancements with a critical angle θ_L<α_c<θ_S result from the excitation of the monopole (n = 0) resonance mode. Moreover, additional acoustic backscattering enhancements still occur in the range 1 ?ka? 6 even though the angle of tilt is greater than the Rayleigh wave coupling angle θ_R=sin^-1(c/c_R) (where c_R is the Rayleigh wave velocity in an elastic half-space). The resonance scattering theory shows that such additional enhancements are associated with the excitation of a dipole (n = 1) resonance mode which may result from the interference of meridional and/or helical waves propagating along the cylinder’s surface. It is therefore essential to consider tilt angles ranging from normal to end-on incidence for a complete analysis of the backscattering by elastic and viscoelastic cylinders.     

标准线性固体材料中球面应力波传播特征研究

基于标准线性固体模型, 结合球面波波动方程, 给出了球面应力波的粒子速度v、粒子位移u、径向应力σ_r、切向应力σ_θ、径向应变ε_r、切向应变ε_θ、折合速度势、折合位移势在Laplace域的理论解. 采用基于Crump算法的Laplace数值逆变换方法分析了上述物理量的传播特征. Laplace数值反演结果表明, 线黏弹性材料对强间断球面应力波的初始响应为纯弹性响应, 强间断在传播过程中包含了几何衰减和本构黏性衰减, 应力、应变、粒子速度的衰减特性和粒子位移、应力、应变、折合位移势等物理量的稳态值同黏弹性球面波的理论预测一致. 折合速度势和折合位移势的峰值随波传播距离的增加逐渐衰减, 这与理想弹性理论给出的折合速度势和折合位移势不随传播距离变化的结论不同. 折合位移势的稳态值与介质的静态剪切模量成反比, 与稳态空腔压力成正比, 与空腔半径的三次方成正比.     

基于衍射条纹非对称性的表面波衰减系数的测量

通过采用将入射光掠入射到频率为几百赫兹的衰减低频液体表面波上,观察到了清晰的间距分布具有明显非对称性的衍射条纹。当入射光掠入射到低频液体表面上时,衍射图样非对称分布具有普遍性。理论上得到了考虑表面波衰减影响后观察屏处衍射光强分布的解析表达式。理论结果表明:衍射图样的这种非对称分布与衍射级次、入射光波波长、液体表面波波长及表面波衰减系数有关;正负级次衍射亮条纹距离中央零级亮条纹的间距体现表面波衰减信息;正衍射级次之间以及负衍射级次之间的条纹间距体现表面波波长信息。利用所提出的方法实现了对几百赫兹液体表面波的衰减系数的实时便捷测量。     

Spin wave propagation in the domain state of a random field magnet

Inelastic neutron scattering with high wave-vector resolution has characterized the propagation of transverse spin wave modes near the antiferromagnetic zone center in the metastable domain state of a random field Ising magnet. A well-defined, long wavelength excitation is observed despite the absence of long-range magnetic order. Direct comparisons with the spin wave dispersion in the long-range ordered antiferromagnetic state reveal no measurable effects from the domain structure. This result recalls analogous behavior in thermally disordered anisotropic spin chains but contrasts sharply with that of the phonon modes in relaxor ferroelectrics. Received 2 November 2002 / Received in final form 4 February 2003 Published online 11 April 2003 RID="a" ID="a"leheny@pha.jhu.edu     

Fano line shape and anti-crossing of Raman active E2g peaks in the charge density wave state of NbSe2

Low energy Raman scattering from the ab-plane of the 2H polytype single crystal NbSe₂ has been investigated in the normal (N), incommensurate charge density wave (ICDW) and superconducting (SC) phases. The temperature dependence of the polarization resolved Raman response has been obtained for the excitation wavelength of 647 nm and fitted to phenomenological models for the E_2g and A_1g symmetry channels. The A_1g response can be fitted by a simple damped oscillator peak superimposed on continuous background. The E_2g response displays an anti-resonance interference pattern between the inter-layer phonon and the CDW-induced mode such that a hybridized configuration (Fano line shape [1]) is required for modelling. The polarization specific peak maxima positions and line widths as a function of temperature, deduced in this manner, are presented. Partial suppression of the electronic continuum scattering in the Raman shift range up to 110 cm⁻¹ in the A_1g symmetry channel and beyond 300 cm⁻¹ in the E_2g symmetry channel is indicative of high energy electronic states far away from the Fermi surface participating in the ICDW formation.