化学反应和球面聚心爆轰波的相互作用

考察球面爆轰波聚心传播过程中,波阵面附近压力和温度不断升高引起化学反应进程的改变;对比氢氧可燃气体与氮气的数值模拟结果,分析化学反应对波面温度和压力的影响,从而考察Zeldovich理论预测聚心爆轰波后参数的精确性.数值结果表明,爆轰波聚心传播初期,放热的燃烧反应对波后热力学参数起主导作用;传播后期,波阵面趋近于对称中心时,吸热的气体解离反应变得非常活跃,解离反应对后期的汇聚压力影响不大,但会在很大程度上限制汇聚温度的升高.     

一维爆轰波在扰动来流中传播的数值研究

爆轰波在静止气体或定常来流中的传播得到了广泛研究, 然而在扰动来流中的传播研究较少。这方面的研究不仅是爆轰传播机制的重要组成部分, 还可为爆轰发动机的应用提供参考。文章基于两步诱导-放热总包反应模型, 开展了一维爆轰波在正弦密度扰动来流中的传播数值模拟。通过对数值结果分析, 获得了放热反应控制参数与爆轰波内在不稳定性的关系, 并在此基础上研究了扰动波长和幅值对一维爆轰波动力学过程的影响。研究发现, 在波前施加连续扰动会诱导爆轰波表现出更复杂的动力学行为, 且影响过程与爆轰波的内在不稳定性相关。对于稳定爆轰波, 扰动只在特定波长范围内引起前导激波后的压力振荡。对于不稳定爆轰波, 扰动会进一步强化其内在不稳定性。扰动幅值越大, 对爆轰波动力学过程的影响越显著。      

Numerical investigation on detonation cell evolution in a channel with area-changing cross section

The two-dimensional cellular detonation propagating in a channel with area- changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.     

气相爆轰波沿胞格的动力学机理研究

 建立了一种以子单元分析为基础,研究气相爆轰波沿胞格运动时的动力学机理的新方法。根据该子单元的性质和斜冲击波关系,首先推导了对撞前后前导冲击波沿胞格对称轴的马赫数之比和入射冲击波入射角及胞格几何性质的关系,求解了胞格结构中的三波点对撞问题。然后,采用爆炸波模拟前导冲击波的自持运动过程,求解气相爆轰波沿胞格的动力学过程,理论分析表明,气相爆轰波在胞格起点首先经历一个增长过程,然后才出现衰减。理论分析结果与实验和数值计算结果的比较表明符合得较好。     

The re-initiation of cellular detonations downstream of an inert layer

The current work aims to examine how the nature of cellular instabilities controls the re-initiation capability and dynamics of a gaseous detonation transmitting across a layer of inert (or non-detonable) gases. This canonical problem is tackled via computational analysis based on the two-dimensional, reactive Euler equations. Two different chemical kinetic models were used, a simplified two-step induction-reaction model and a detailed model for hydrogen-air. For the two-step model, cases with relatively high and low activation energies, representing highly and weakly unstable cellular detonations, respectively, are considered. For the weakly unstable case, two distinct types of re-initiation mechanisms were observed. (1) For thin inert layers, at the exit of the layer the detonation wave front has not fully decayed and thus the transverse waves are still relatively strong. Detonation re-initiation in the reactive gas downstream of the inert layer occurs at the gas compressed by the collision of the transverse waves, and thus is referred to as a cellular-instability-controlled re-initiation. (2) If an inert layer is sufficiently thick, the detonation wave front has fully decayed to a planar shock when it exits the inert layer, and re-initiation still occurs downstream as a result of planar shock compression only, which is thus referred to as a planar-shock-induced re-initiation. Between these two regimes there is a transition region where the wave front is not yet fully planar, and thus perturbations by the transverse waves still play a role in the re-initiation. For the highly unstable case, re-initiation only occurs via the cellular-instability-controlled mechanisms below a critical thickness of the inert layer. Additional simulations considering detailed chemical kinetics demonstrate that the critical re-initiation behaviors of an unstable stoichiometric mixture of hydrogen-air at 1 atm and 295 K are consistent with the finding from the two-step kinetic model for a highly unstable reactive mixture.     

凝聚炸药爆轰波在高声速材料界面上的折射现象分析

凝聚炸药爆轰在边界高声速材料约束下传播时,爆轰波会在约束材料界面上产生复杂的折射现象.本文针对凝聚炸药爆轰波在高声速材料界面上的折射现象展开理论和数值模拟分析.首先通过建立在爆轰ZND模型上的改进爆轰波极曲线理论给出爆轰波折射类型,然后发展一种求解爆轰反应流动方程的基于特征理论的二阶单元中心型Lagrange计算方法来数值模拟典型的爆轰波折射过程.从改进爆轰波极曲线理论和二阶Lagrange方法数值模拟给出的结果看出,凝聚炸药爆轰波在高声速材料界面上的折射类型有四种:反射冲击波的正规折射、带束缚前驱波的非正规折射、带双Mach反射的非正规折射、带λ波结构的非正规折射.     

Nonlinear dynamics of self-sustained supersonic reaction waves: Fickett's detonation analogue

The present study investigates the spatiotemporal variability in the dynamics of self-sustained supersonic reaction waves propagating through an excitable medium. The model is an extension of Fickett's detonation model with a state-dependent energy addition term. Stable and pulsating supersonic waves are predicted. With increasing sensitivity of the reaction rate, the reaction wave transits from steady propagation to stable limit cycles and eventually to chaos through the classical Feigenbaum route. The physical pulsation mechanism is explained by the coherence between internal wave motion and energy release. The results obtained clarify the physical origin of detonation wave instability in chemical detonations previously observed experimentally.     

Detonation re-initiation mechanism following the Mach reflection of a quenched detonation

This experimental study addresses the re-initiation mechanism of detonation waves following the Mach reflection of a shock–flame complex. The detonation diffraction around a cylinder is used to reproducibly generate the shock–flame complex of interest. The experiments are performed in methane–oxygen. We use a novel experimental technique of coupling a two-in-line-spark flash system with a double-frame camera in order to obtain microsecond time resolution permitting accurate schlieren velocimetry. The first series of experiments compares the non-reactive sequence of shock reflections with the reflection over a rough wall under identical conditions. It was found that the hot reaction products generated along the rough wall are entrained by the wall jet into a large vortex structure behind the Mach stem. The second series of experiments performed in more sensitive mixtures addressed the sequence of events leading to the detonation establishment along the Mach and transverse waves. Following ignition and jet entrainment, a detonation first appears along the Mach stem while the transverse wave remains non-reactive. The structure of the unburned tongue however indicates local instabilities and hot spot formation, leading to the rapid reaction of this gas. Numerical simulations are also reported, confirming the sequence of ignition events obtained experimentally.     

Nonequilibrium Zeldovich-von Neumann-Doring theory and reactive flow modeling of detonation

This paper discusses the Nonequilibrium Zeldovich-von Neumann-Doring (NEZND) theory of self-sustaining detonation waves and the Ignition and Growth reactive flow model of shock initiation and detonation wave propagation in solid explosives. The NEZND theory identified the nonequilibrium excitation processes that precede and follow the exothermic decomposition of a large high explosive molecule into several small reaction product molecules. The thermal energy deposited by the leading shock wave must be distributed to the vibrational modes of the explosive molecule before chemical reactions can occur. The induction time for the onset of the initial endothermic reactions can be calculated using high pressure-high temperature transition state theory. Since the chemical energy is released well behind the leading shock front of a detonation wave, a physical mechanism is required for this chemical energy to reinforce the leading shock front and maintain its overall constant velocity. This mechanism is the amplification of pressure wavelets in the reaction zone by the process of de-excitation of the initially highly vibrationally excited reaction product molecules. This process leads to the development of the three-dimensional structure of detonation waves observed for all explosives. For practical predictions of shock initiation and detonation in hydrodynamic codes, phenomenological reactive flow models have been developed. The Ignition and Growth reactive flow model of shock initiation and detonation in solid explosives has been very successful in describing the overall flow measured by embedded gauges and laser interferometry. This reactive flow model uses pressure and compression dependent reaction rates, because time-resolved experimental temperature data is not yet available. Since all chemical reaction rates are ultimately controlled by temperature, the next generation of reactive flow models will use temperature dependent reaction rates. Progress on a statistical hot spot ignition and growth reactive flow model with multistep Arrhenius chemical reaction pathways is discussed. The text was submitted by the authors in English.     

Generation of magnetic field by detonation waves

The generation of a magnetic field by a system of detonation waves in a condensed explosive is reported. The convergence of the detonation waves, which exhibit a high conductivity in the chemical reaction zone, increases the magnetic field at the axis of the system. The fact of magnetic field generation is demonstrated experimentally. Features of the new method of magnetic cumulation are discussed. A simple compression model that qualitatively agrees with experimental data is proposed.     

Features of the dynamics of postdetonation waves

We present the results of numerical investigations of the parameters of postdetonation waves forming at a passage from the zone occupied with a bubbly liquid formed by the detonation wave to a zone filled with a liquid without bubbles. The dependence of the pressure amplitude of detonation waves and postdetonation waves on the gas volumetric content of bubbles has been studied. A possibility of the detonation transfer through the layer of a bubble-free liquid separating the regions of the bubbly liquid has been shown, the map of possible situations at the detonation transfer through the layer of this liquid has been presented.     

Detonation development in hydrogen/air mixtures inside a closed chamber: role of a cold wall

Detonation development from a hot spot has been extensively studied, where ignition occurs earlier than that in the surrounding mixtures. It has also been reported that a cool spot can induce detonation for large hydrocarbon fuels with Negative Temperature Coefficient (NTC) behavior, since ignition could happen earlier at lower temperatures. In this work we find that even for hydrogen/air mixtures without NTC behaviors, a cold wall can still initiate and promote detonation. End-wall reflection of the pressure wave and wall heat loss introduce an exothermic center outside the boundary layer, and then autoignitive reaction fronts on both sides may evolve into detonation waves. The right branch can be further strengthened by appropriate temperature gradient near the cold wall, and exhibits different dynamics at various initial conditions. The small excitation time and the large diffusivity of hydrogen provide the possibility for detonation development within the limited space between the autoignition kernel and the cold wall. Moreover, detonation may also develop near the flame front, which may or may not co-exist with detonation waves from the cold wall. Correspondingly, wall heat flux evolution exhibits different responses to detailed dynamic structures. Finally, we propose a regime diagram describing different combustion modes including normal flame, autoignition, and detonation from the wall and/or the reaction front. The boundary of normal flame regime qualitatively agrees with the prediction by the Livengood-Wu Integral method, while the detonation development from both the end wall and the reaction front observes Zel'dovich mechanism. Compared to hydrocarbons, hydrogen is resistant to knock onset but it is more prone to superknock development. The latter mode becomes more destructive in the presence of wall heat loss. This study isolates and identifies the role of wall heat loss on a potential mechanism for superknock development in hydrogen-fueled spark-ignition engines.     

Quantum effects in the kinetics of the initiation of detonation condensation waves

The features of the kinetics of the initiation of detonation condensation waves in carbon suboxide and acetylene have been experimentally studied at high pressures near the low-temperature limits. The role of quantum effects in the expansion of detonation limits has been analyzed. Quantum corrections to the endothermic reaction rates, which are caused by an increase in the high-energy tails of the momentum distribution functions at high pressures due to the manifestation of the uncertainty principle for the energy of colliding particles at a high collision frequency, have been quantitatively estimated. It has been shown that experimentally observed deviations in Arrhenius dependences of the induction periods of the initiation of detonation condensation waves are well described by the proposed quantum corrections.     

Peculiarities of the variation of the detonation front surface area for DINA-Acetonitrile charges with a near-critical diameter

The previously revealed similarity between the passage of detonation in a weak liquid HE (a 68.5: 31.5 wt % DINA-A solution, close to the critical one) from a tube into a unconfined volume and in powerful liquid HE, such as nitroglycerin is discussed. It was concluded that the similarity in the behavior of liquid HEs so different in power is not associated with the kinetic stability of the detonation front and is determined by the phenomenon of the termination of the chemical reaction in the front of the corresponding waves by rarefaction waves.     

Uniform Acoustic Cavitation of Liquid in a Disk

A dynamical propagation model coupled to the oscillation of cavitation bubbles is applied to describe the imploding acoustic field in a cavitating liquid where the acoustic waves transmit from the outside to the inside of a circle disk. Numerical simulation shows that the imploding ability of a ring source can elevate the sound pressure or partly eliminate the decay due to both the bulk attenuation and the attenuation caused by cavitation. However,the imploding ability is limited and there exists a critical radius. When the radius of the disk is larger than the critical one, the imploding ability is not enough to eliminate the attenuation. Fortunately, the cavitation region can be effectively expanded if a hot plate is attached under the center of the disk because the cavitation threshold is related to the temperature of the liquid, which means that a region with good uniformity of cavitation can be enhanced by adjusting the temperature difference between the central and side liquid.     

Propagation of shock and detonation waves in channels with U-shaped bends of limiting curvature

Systematic experimental and theoretical studies of the propagation of shock and detonation waves in cylindrical tubes and planar channels with two U-shaped bends of limiting curvature were performed. It was demonstrated that U-shaped bends substantially facilitate detonation initiation in gases. The minimum shock wave velocity required to initiate the detonation of a stoichiometric propane-air mixture under normal conditions in a near-critical diameter tube with two U-shaped bends of limiting curvature was found to be ～800 m/s.     

Formation of transverse waves in oblique detonations

The structure of oblique detonation waves stabilized on a hypersonic wedge in mixtures characterized by a large activation energy is investigated via steady method of characteristics (MoC) calculations and unsteady computational flowfield simulations. The steady MoC solutions show that, after the transition from shock-induced combustion to an overdriven oblique detonation, the shock and reaction complex exhibit a spatial oscillation. The degree of overdrive required to suppress this oscillation was found to be nearly equal to the overdrive required to force a one-dimensional piston-driven detonation to be stable, demonstrating the equivalence of two-dimensional steady oblique detonations and one-dimensional unsteady detonations. Full unsteady computational simulations of the flowfield using an adaptive refinement scheme showed that these spatial oscillations are transient in nature, evolving in time into transverse waves on the leading shock front. The formation of left-running transverse waves (facing upstream) precedes the formation of right-running transverse waves (facing downstream). Both sets of waves are convected downstream away from the wedge in the supersonic flow behind the leading oblique front, such that the mechanism of instability must continuously generate new transverse waves from an initially uniform flow. Together, these waves define a cellular structure that is qualitatively similar to a normal propagating detonation.     

Transverse wave generation mechanism in rotating detonation

Detonation engines are expected to be included in a number of aerospace thrusters in the future. Several types of detonation engines are currently under examination, including the rotating detonation engine (RDE). Although the RDE has been explored experimentally, its rotating detonation propagation mechanism is not well understood. This paper clarifies the detonation mechanism and dynamics of the RDE by 2D and 3D simulation using compressible Euler equations with a full chemical reaction mechanism of H₂/O₂ and H₂/Air, especially from the triple-point and transverse detonation points of view. A total variation diminishing (TVD) scheme is used for the mixture of H₂/Air, and an advection upwind splitting method difference vector (AUSMDV) scheme is used for the mixture of H₂/O₂. The use of an AUSMDV scheme provides a much clearer detonation structure than does the TVD scheme. We focus on the complex interaction mechanism of the detonation front and burned mixture gases. We found out that at this interaction point, an unreacted gas pocket appears and ignites periodically to generate transverse waves at the detonation front and maintain detonation propagation.     

一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解：显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动.     

Experimental study of the relief of the front of steady detonation in liquid high explosives

The instability of detonation waves in a 81 : 19 nitromethane-acetone mixture is experimentally studied. The amplitude of the irregularities caused by detonation instability and the typical reaction time are estimated using a VISAR laser interferometer, which makes it possible to measure the particle velocity profiles at the explosive-window interface. The characteristic transverse dimensions of the irregularities are measured with the help of a CORDIN camera operating in the single shooting mode to record the intrinsic emission of the detonation front or the light of an outside source reflected from it.