High resolution numerical simulation of the structure of 2-D gaseous detonations

Two dimensional numerical simulation of the structure of gaseous detonation is investigated by utilizing the single step Arrhenius kinetic reaction mechanism in both high and low activation energy mixtures, characterized by their irregular and regular detonation structure, respectively. All the computations are performed on a small Beowulf cluster with six nodes. The dependency of the structure on the grid resolution is performed and it is found that, resolution of more than 300 cells per hrl is required to demonstrate the role of hydrodynamic instabilities, (KH and RM instabilities) in detonation propagation in irregular structures, while due to the absence of fine-scale structures, resolution of 50 cells per hrl, gives the physical structure of detonation with regular structures. Results show that the transverse waves in irregular structure are significantly stronger than the transverse wave in regular structure detonation, which can enhance the burning rate of the unburned pockets behind the shock front. Results for resolution of 600 cells per hrl illustrate that, in addition to the primary mode, the interaction of large vortices with the shock front provides secondary modes in the structure which leads to the irregularity of the structure in high activation energy mixture. In contrast with the results obtained for regular structure, which no unburned gas pockets and vortices observed behind the front, the results for irregular structure reveal that most portions of the gases, escape from shock compression and create large unburned gas pockets behind the both weak section of the Mach stem and the incident wave, which will burn eventually by the turbulent mixing due to the vortices associated with hydrodynamic instabilities. Therefore, the ignition mechanism in irregular structure is due to the both shock compression and by turbulent mixing associated with hydrodynamic instabilities, while the shock compression yields the ignition mechanism in regular structure detonation.     

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.     

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

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

Cellular Cell Bifurcation of Cylindrical Detonations

Cellular cell pattern evolution of cylindrically-diverging detonations is numerically simulated successfully by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. From thesimulation, three cell bifurcation modes are observed during the evolution and referred to as concave front focusing, kinked and wrinkled wave front instability, and self-merging of cellular cells. Numerical research demonstrates that the wave front expansion resulted from detonation front diverging plays a major role in the cellular cell bifurcation, which can disturb the nonlinearlyself-sustained mechanism of detonations and finally lead to cell bifurcations.     

无粘流中能捕捉较弱的接触间断和激波的反扩散差分格式

在研究弱入射激波遇到对称楔以后的马赫反射现象时,激波管实验不易测出很弱的接触间断,也不易捕捉到马赫反射与正规反射转换的条件。文章一方面研究了可压流体力学欧拉方程的数值方法,首先是用反扩散法改进接触间断的计算;另一方面根据格式粘性的特性和它引出的很微小的熵的变化规律来显示很弱的接触间断和反射激波。这样才易于将对三波点的分析推进一步。文[5,6]曾预言了一种反散波是连续的压缩波的新的激波反射类型。我们设想并根据计算初步确认这新类型反射实际应该是简单马赫反射,反射波虽弱仍是激波。     

Numerical study of three-dimensional detonation wave dynamics in a circular tube

The three-dimensional structures of a detonation wave propagating in a circular tube were investigated using a one-step irreversible Arrhenius kinetics model. A series of parametric studies were carried out to investigate the different modes of cell structure formation by changing the pre-exponential factor. Maximum pressure trace was recorded along the tube wall to investigate the detonation cell structures. The unsteady results obtained in three dimensions revealed the generation mechanism of the wave front structures of two-, three- and four-cell mode detonations. A six-cell mode detonation could be obtained using a finer grid. With the increase in pre-exponential factor, it was found that the number of detonation cells is increased while the cell size is reduced accordingly. In all the multi-cell modes, the detonation wave structures and smoked-foil records on the wall are formed by the propagation of transverse waves along the wall in clockwise and counter-clockwise directions, while the slapping wave moves in the radial direction. The presence of the slapping wave further strengthens the wave interactions in three-dimensional simulation. Comparison with two-dimensional simulation confirms the effect of the slapping wave in the radial direction. As a result, the detonation wave front structures changes from the polygonal shape to the multi-bladed fan shape, periodically.     

Critical diffraction of irregular structure detonations and their predictability from experimentally obtained D−κ data

The present work reports new experiments of detonation diffraction in a 2D channel configuration in stoichiometric mixtures of ethylene, ethane, and methane with oxygen as oxidizer. The flow field details are obtained using high-speed schlieren near the critical conditions of diffraction. The critical initial pressure for successful diffraction is reported for the ethylene, ethane and methane mixtures. The flow field details revealed that the lateral portion of the wave results in a zone of quenched ignition. The dynamics of the laterally diffracting shock front are found in good agreement with the recent model developed by Radulescu et al. (Physics of Fluids 2021). The model provides noticeable improvement over the local models using Whitham’s characteristic rule and Wescott, Bdzil and Stewart’s model for weakly curved reactive shocks. These models provide a link between the critical channel height and the critical wave curvature. The critical channel heights and global curvatures are found in very good agreement with the critical curvatures measured independently by Xiao and Radulescu (Combust. Flame 2020) in quasi-steady experiments in exponential horns for three mixtures tested. Furthermore, critical curvature data obtained by others in the literature was found to provide a good prediction of critical diffraction in 2D. These findings suggest that the critical diffraction of unstable detonations may be well predicted by a model based on the maximum curvature of the detonation front, where the latter is to be measured experimentally and account for the role of the cellular structure in the burning mechanism. This finding provides support to the view that models for unstable detonations at a meso-scale larger than the cell size, i.e., hydrodynamic average models, are meaningful.     

Stability and ambiguous representation of shock wave discontinuity in thermodynamically nonideal media

The nonlinear analysis of the behavior of a shock wave on a Hugoniot curve fragment that allows for the ambiguous representation of shock wave discontinuity has been performed. The fragment under consideration includes a section where the condition L > 1 + 2M is satisfied, which is a linear criterion of the instability of the shock wave in media with an arbitrary equation of state. The calculations in the model of a viscous heat-conductive gas show that solutions with an instable shock wave are not implemented. In the one-dimensional model, the shock wave decays into two shock waves or a shock wave and a rarefaction wave, which propagate in opposite directions, or can remain in the initial state. The choice of the solution depends on the parameters of the shock wave (position on the Hugoniot curve), as well as on the form and intensity of its perturbation. In the two-dimensional and three-dimensional calculations with a periodic perturbation of the shock wave, a “cellular” structure is formed on the shock front with a finite amplitude of perturbations that does not decrease and increase in time. Such behavior of the shock wave is attributed to the appearance of the triple configurations in the inclined sections of the perturbed shock wave, which interact with each other in the process of propagation along its front.     

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.     

激光驱动的冲击波自生磁场以及外加磁场的冲击波放大研究

冲击波是天体物理观测中常见的现象, 其对粒子的加速被认为是高能宇宙射线的来源. 宇宙中冲击波周围往往存在很强的磁场, 但人们对于此类强磁场的产生放大过程的理解并不充分. 本文利用二维粒子模拟程序研究了激光与磁化或者非磁化等离子体相互作用产生的冲击波现象, 给出了冲击波波前处磁场的产生放大特性. 研究发现, 作用过程中的自生磁场可以储存能量, 从而进一步加速电子; 当存在外加磁场时, 由冲击波加速的电子和离子的能量都比同条件下非磁化等离子体的能量高; 而且外加磁场藉由冲击波放大倍数则与其值有极大关系. 与天文观测中推断的磁场与背景磁场相比放大千倍这一研究结果的比较可以看出, 天体冲击波周围磁场放大主要是由局域内生磁场导致的.     

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

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

多孔脆性介质冲击波压缩破坏的细观机理和图像

本文采用一种具有良好定量性质的离散元模型研究了带孔洞的各向同性脆性介质在细观尺度上的压缩破坏特征. 通过对孤立孔洞、三种简单的孔洞排布方式和大量孔洞随机排布等几种情况的模拟, 认识到了剪切破坏和局域拉伸破坏是冲击波压缩下多孔介质的基本破坏模式; 孔洞之间的损伤贯通会促进孔洞在较低应力下发生塌缩, 但损伤区的应力松弛过程却会对一定范围内的介质起到损伤屏蔽作用; 不同区域中损伤促进和损伤屏蔽的综合效果是在多孔脆性介质中形成一种高损伤区与低损伤区间错排布的奇特损伤分布. 本文的研究结果为深入理解脆性材料冲击波压缩破坏的演化过程和机理提供了细观尺度上的初步物理图像.     

Shock waves and formation of carbon dioxide hydrate at an increased pressure in the gas-liquid medium

The processes of breaking, solution, and formation of hydrates behind a shock wave of moderate amplitude were studied experimentally in water with carbon dioxide bubbles under different initial static pressures. It is shown that an increase in the static pressure in a gas-liquid medium leads to reduction of critical relative amplitude of the shock wave, corresponding to starting development of Kelvin — Helmholtz instability and bubble splitting into small gas inclusions behind the shock wave front. It is shown that the rates of carbon dioxide solution and hydrate formation behind the shock wave front are close by the value; their dependences on medium and wave parameters are determined. Calculations by the model of gas hydration behind the shock wave are presented. The work was financially supported by the Russian Foundation for Basic Research (grants Nos. 06-01-00142 and 06-08-00657).     

Shock wave chemistry

Abstract

Shock wave chemistry, a new scientific trend, deals with investigations of chemical aspects of the substance state under this new type of effect. Indeed, shock wave effect is not a greater imposition than pressure and temperature actions. Characteristic features of the effect are the tremendous rates of substance loading and subsequent unloading. The effects result in a substance in a strongly non- equilibrium state. The lifetime of the state is governed by the relaxation process of those phenomena which are provoked by shock waves in the substance. For instance, in the case of substance consisting of complex molecules with a large number of internal degrees of freedom, differing strongly in excitation times, all kinetic parts of the shock energy are at first absorbed by the translational degrees of freedom inside the shock wave front. Then, the energy is redistributed to the vibrational degrees of freedom. The non-equilibrium state time is not longer than the excitation time of the most slowly excited vibrational degrees of freedom (10¹⁰-10^?9 s). The same order of magnitude is the relaxation time of liquid substance polarization caused by dipolar molecules mechanically turning under the shock discontinuity zone effect. In polymers the zone turns some separate groups of polymer molecule atoms. In such a case the relaxation period, on the contrary, may last as long as it can. As far as “hot are concerned, their lifetime is determined by thermal relaxation regularities and it depends on their size. The hot spots in solids appear during the shock compression process at the sites of an imperfect substance structure. In liquids the hot spots can orighate when a shock wave front passes through negative density fluctuations. It transforms the fluctuations of very small size and of high probability into some positive temperature regions of large size and extremely low probability at equilibrium state behind the wave front. The hot spots in perfect solids (possibly in liquids too) appear due to the effect of shear stresses in shock front. Pointed and lengthy defects of solid structure occur under the effect. The lengthy defects appear in the shock wave front due to the transition from one-dimensional to volume compression. The transition takes place if the wave intensity is larger than the dynamic elastic limit of the solid under investigation. In brittle materials the transition results in their grinding into fragments and in the relative displacement of the fragments. Some liquid melted layers of substance appear between the fragments in the process of displacement. Their lifetime is also determined by the thermal relaxation regularities and probably is small. Nevertheless, the layers obviously govern the spall strength of brittle solids and promote solid-phase shock reactions. The defects created in solids by the shock effect can exist for a very long time if the solid substance residual temperature is lower than its recrystallization temperature. Therefore, solid substance treatment by shocks of proper intensity can increase their chemical reactivity.     

Observation of indirect parallel pumping of magneto-elastic modes in layered YIG/GGG structures

In this work we report new experimental results on nonlinear excitation of magnetoelastic (ME) modes in layered YIG/GGG waveguide structures at GHz frequencies, obtained by a guided-wave light scattering technique. It is shown that the fundamental spin-dipole wave (SDW) mode induces a secondary microwave field at double frequency that can efficiently excite shear elastic modes of the structure. The distinctive feature of a mechanism for ME coupling proposed is that it is free from the selection rules and provides continuous excitation of elastic modes within a wide frequency range by means of the standard microstripe line excitation system.     

Rayleigh mechanism of local elevation of energy in propagation of a shock wave in an active medium

A mechanism of action of a shock wave on an active medium, which leads to an additional energy release source, is considered. This source moves together with the shock wave front and depends on the magnitude and direction of the electric field applied to the plasma and on the current density in the plasma. The study is a continuation of an earlier publication devoted to the propagation of weak shock waves. Here, we consider shock waves of an arbitrary intensity with an arbitrary mechanism of formation of an additional energy contribution due to variation of the parameters of the medium as a result of passage of the shock wave. Special cases of this effect are the propagation of a shock wave in a plasma and detonation burning.     

On the mechanism of the deflagration-to-detonation transition in a hydrogen-oxygen mixture

The flame acceleration and the physical mechanism underlying the deflagration-to-detonation transition (DDT) have been studied experimentally, theoretically, and using a two-dimensional gasdynamic model for a hydrogen-oxygen gas mixture by taking into account the chain chemical reaction kinetics for eight components. A flame accelerating in a tube is shown to generate shock waves that are formed directly at the flame front just before DDT occurred, producing a layer of compressed gas adjacent to the flame front. A mixture with a density higher than that of the initial gas enters the flame front, is heated, and enters into reaction. As a result, a high-amplitude pressure peak is formed at the flame front. An increase in pressure and density at the leading edge of the flame front accelerates the chemical reaction, causing amplification of the compression wave and an exponentially rapid growth of the pressure peak, which “drags” the flame behind. A high-amplitude compression wave produces a strong shock immediately ahead of the reaction zone, generating a detonation wave. The theory and numerical simulations of the flame acceleration and the new physical mechanism of DDT are in complete agreement with the experimentally observed flame acceleration, shock formation, and DDT in a hydrogen-oxygen gas mixture.