缸内激波对锥顶型燃烧室的活塞破坏机理

以二维数值模拟为基础,研究了锥顶型燃烧室内的冲击波发展的震荡过程,得到作用于活塞不同位置处的超压分布。模拟结果表明:由于燃烧室结构的独特性,导致冲击波能在特定区域进行汇聚,致使该区域超压明显高于其他区域。将该模拟结果与实际破坏失效的活塞进行对比,发现冲击波汇聚区域往往就是活塞被破坏的地方。数值模拟结果与实际破坏结果吻合很好。这为设计燃烧室形状以避免冲击波对活塞造成破坏提供了理论基础。     

Formation of blast waves as a result of a low pressure domain decomposition

In this article the formation and propagation of blast wave as a result of a focusing of shock wave in a domain with low pressure and density are examined in the frame of an ideal gas model. We consider the decomposition of a pressure-density discontinuity on the boundary of the spherical (or cylindrical) domain which is filled by a gas whose pressure and density are lower than the pressure and density of gas that filling the external space. At initial moments of this decomposition the rarefaction wave propagates in the external space and the converging shock wave is formed in the domain of low pressure and density. The intensity of the converging shock wave gradually increases, and the wave transitions to the self-similar regime. After implosion, a diverging shock wave is formed which propagates through the disturbed gas. The values of wave intensity, wave impulse and other parameters for some magnitudes of initial parameters have been determined by means of numerical calculations. Received 10 August 1997 / Accepted 13 July 1998     

Numerical investigation of shock wave reflections near the head ends of rotating detonation engines

The influence of various chamber geometries on shock wave reflections near the head end of rotating detonation engines was investigated. A hydrogen/air one-step chemical reaction model was used. The results demonstrated that the variation in flow field along the radial direction was not obvious when the chamber width was small, but became progressively more obvious as the chamber width increased. The thrust increased linearly, and the detonation height and the fuel-based gross specific impulse were almost constant as the chamber width increased. Near the head end, shock waves reflected repeatedly between the inner and outer walls. Both regular and Mach reflections were found near the head end. The length of the Mach stem increased as the chamber length increased. When the chamber width, chamber length and injection parameters were the same, the larger inner radius resulted in more shock wave reflections between the inner and outer walls. The greater the ratio of the chamber width to the inner radius, the weaker the shock wave reflection near the head end. The detonation height on the outer wall and the thrust, both increased correspondingly, while the specific impulse was almost constant as the inner radius of the chamber increased. The numerical shock wave reflection phenomena coincided qualitatively with the experimental results.     

The detonation burning of a fuel mixture resulting from a double explosion

A study is made of the propagation of a multifront detonation burning in a fuel mixture consisting of a gaseous fuel and an oxidant with additions of combustible solid or liquid particles arising as a result of a double point explosion. In such combustible media it is possible for there to be propagation of several detonation or burning fronts following one after the other. The easily igniting gaseous fuel burns in the first detonation wave, which propagates in the gaseous mixture with particles which are heated by the products of the explosion, ignite and burn in the second detonation wave or in the flame front. Self-similar regimes of propagation of such waves in an idealized formulation were studied in [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–131, March–April, 1985.     

Numerical study on three-dimensional flow field of continuously rotating detonation in a toroidal chamber

Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonationshock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interesting properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.     

舱内爆炸角隅汇聚反射冲击波超压特性研究

为研究密闭舱室内爆角隅汇聚反射冲击波超压特性,利用缩比模型进行了某典型舱室内爆试验,得到远离角隅、两面角隅和三面角隅处的冲击波载荷,结合数值模拟研究了3种特征位置处冲击波传播规律及载荷特征。研究结果表明:远离角隅处壁面反射冲击波超压曲线呈现单峰结构,反射冲击波以球面波传播;距两面角隅一定范围内冲击波超压曲线呈现双峰结构,两面角隅冲击波超压曲线呈现单峰结构,角隅汇聚反射冲击波以椭球状传播;距三面角隅一定范围内冲击波超压曲线呈现多峰结构,三面角隅冲击波超压曲线呈现单峰结构,角隅汇聚反射冲击波以球面波传播;在合理假设条件下,根据量纲分析及数值模拟结果,得到首次冲击时角隅汇聚反射冲击波载荷经验计算公式。     

Coupling Relation Between Air Shockwave and High-Temperature Flow from Explosion of Methane in Air

After a methane-in-air explosion in a coal mine tunnel, a secondary explosion of coal dust is prone to happen. The shockwave in the gas explosion produces a coal dust suspension, and the peak temperature band may detonate that suspension. This secondary detonation depends on the space-time relation between the shockwave and the peak temperature band. This paper presents a methodology to estimate the coupling relation between the air shockwave and high-temperature flow from the explosion of methane in air. The commercial software package AutoReaGas was used to carry out the numerical simulation for the explosion processes of methane in air in the tunnel. Based on the numerical simulation and its analysis, the coupling relation between the leading shock wave and high-temperature flow was demonstrated for a methane-in-air explosion in a tunnel. In the near field of the ignition point, the deflagration wave transmits energy by heat, and the temperature load is in the front of the pressure wave. With development of deflagration and deflagration-to-detonation transition, the corresponding mechanism of energy transmission is changed from heat conduction to shock compression, and a precursor pressure wave is formed gradually. The time interval between the precursor pressure wave and high-temperature flow behind the wave increases with distance. Attenuation of the precursor shock wave and high-temperature flow depends on the length of the methane-in-air space in a tunnel. Beyond the methane-in-air space, the quantitative relation of the time interval between the precursor shock wave and high-temperature flow with axial distance from ignition and the length of methane-in-air space was proposed.     

增强型喷射器对爆轰波DDT过程的影响

在激波、火焰及射流同时存在的流场中，组织燃烧转爆轰过程是脉冲爆震发动机实现点火、起爆的关键问题。设计一类喷射器，采用C₂H₂/O₂/Ar反应，数值验证了该喷射器能增强爆震室燃料燃烧转爆轰的可行性，并讨论了流场中热点的点火机制。结果显示:该装置在流场中可激发不稳定性，产生漩涡，加速能量、质量的交换。流场产生热点，促进火焰速度加快，追赶前导激波。喷射器位置影响前导激波的运动速度。在一定范围内，前导激波速度越大，碰撞产生的热点越容易激发燃烧转爆轰过程。     

Behavior of detonation propagation in mixtures with concentration gradients

Behavior of detonation waves in mixtures with concentration gradients normal to the propagation direction was studied experimentally. Mixtures with various concentration gradients were formed by sliding the separation plate which divides a detonation chamber from a diffusion chamber in which a diffusion gas was initially introduced. A stoichiometric hydrogen–oxygen mixture was charged in the detonation chamber, while oxygen or nitrogen was filled in the diffusion gas chamber. Temporal concentration measurement was conducted by the infrared absorption method using ethane as alternate of oxygen. Smoked foil records show a deformation of regular diamond cells to parallelogram ones, which well corresponds to local mixture concentration. Schlieren photographs reveal the tilted wave front whose angle is consistent with the deflection angle of the detonation front obtained from trajectories of the triple point. The local deflection angle increases with increase in local concentration gradient. Calculation of wave trajectory based on the ray tracing theory predicts formation of the tilted wave front from an initial planar front.      

爆震室壁面条件对脉冲爆震发动机爆震性能的影响

在爆震室内快速形成稳定传播的爆轰波是脉冲爆震发动机的关键．本文利用有限速率化学反应模型,考虑粘性、热对流,基于N-S方程对氢气与空气/氧气为反应混合物的爆震发动机爆震室内流场进行计算．从流场压力、速度、涡量、湍流动能等方面研究爆震室壁面条件对燃烧爆轰性能的影响,分析流场爆轰波压力与流场湍动能的关系,讨论可燃气体燃烧转爆轰的机理．结果表明：爆震室内燃烧爆轰机理受到化学反应能量释放、壁面摩擦效应、壁面与外界热交换的影响．在文中讨论的范围内,相比于半圆形和三角形的爆震室装置,矩形的爆震室增强装置能在更短的时间内得到较高的爆轰波压力和湍动能峰值．壁面粗糙层高度（粗糙度）影响爆震室的燃烧爆轰性质．当壁面粗糙度为0.15mm时,粗糙度对爆轰的激励作用大于抑制作用,能较快形成稳定的爆轰波,且推力为35.5N;随着壁面对流换热系数的增大,爆震室壁面的散热加剧．当壁面对流换热系数大于临界值2.6W/(m2·K)时,爆震室内不能形成稳定的爆震波．     

调控燃烧室燃料初始分布建立稳定旋转爆轰波的方法

旋转爆轰发动机具有比传统航空航天发动机更高的燃烧效率，近年来引起人们的关注。其中，点火启动过程尤为重要。为达到一次点火就能在燃烧室内建立稳定旋转爆轰波的目的，本文提出通过控制点火前燃料初始分布来建立稳定旋转爆轰波的方法，并基于纳维-斯托克斯方程与10组分27可逆反应基元化学反应模型的数值模拟验证了该方法的可行性。对旋转爆轰波传播特性的研究表明，燃料在发动机燃烧室中的分布是影响旋转爆轰波建立的关键。在燃料喷注压力较低时此影响尤为明显，它决定了爆轰波发展第一周期内波前燃料层厚度。而波前燃料层与波的稳定传播密切相关。基于该方法，本文对燃烧室初始流速为360 m/s，喷注总压0.4 MPa的旋转爆轰发动机实现了点火至稳定爆轰，得到的爆轰波传播平均速度为1 604 m/s，频率为5 347.6 Hz。此外，燃料初始填充率作为燃料初始分布的量化指标，文中给出了它建立稳定旋转爆轰时的临界范围。     

一维与二维爆轰传播的时空关联特性数值研究

一维爆轰传播的理论完备、计算准确, 二维斜爆轰传播由于壁面与黏性效应, 大尺度、高精度预测还有一定难度. 利用Euler方程和H₂-Air基元反应模型, 对二维有限长楔面诱导的斜爆轰和活塞驱动一维非定常正爆轰进行计算比较研究, 从时空两个维度方面, 分析了两者在起爆过程、稀疏波传播、爆轰波面演化中的关联特性. 研究结果表明: 在过驱动度相同的条件下, 经过时空变换的活塞驱动一维爆轰传播与二维驻定斜爆轰在起爆区波系结构、波面演化特征和主要参数分布规律方面无论定性或者定量对比均符合较好, 所以, 一维非定常爆轰和二维驻定斜爆轰具有时空相关性. 两者的差异主要体现在过驱动斜爆轰受稀疏波影响过渡到近Chapman-Jouguet (C-J)爆轰状态所需的弛豫时间不同, 原因可能是起源于活塞和壁面稀疏波强度的差异. 本文提出的一维与二维爆轰传播的时空关联方法不仅有助于认知斜爆轰起爆、过驱爆轰产生、胞格爆轰演化的三阶段规律, 还可以对比揭示壁面、边界层和黏性效应的影响, 应用在斜爆轰发动机燃烧室设计中能够有效节约计算时间和成本, 并降低复杂度.      

On the shock-induced explosion resulting from an impact by a planar detonation

We analyse the processes which occur when a planar detonation propagating from the fixed end of a donor explosive charge impacts on an acceptor homogeneous explosive. We propose a model for estimating the minimal length of the donor charge for which an explosion can be generated in the acceptor. We show that the self-similarity of the donor flow imposes a minimal length on the donor charge so its piston effect be capable of keeping the volumetric-expansion rate of the shocked acceptor to small-enough values and, thereby, of triggering explosion in a finite time. The donor detonation is represented as a Chapman-Jouguet discontinuity; the chemical decomposition in the acceptor is described by the Arrhenius global rate law. The model reproduces the experimental trend according to which the smaller minimal lengths are obtained with donor explosives that have larger heats of reaction and initial pressures. The minimal lengths predicted by the model agree well with those obtained by means of one-dimensional numerical simulations. Additional simulations show that the minimal length for generating an explosion is smaller than, but perhaps of the same order as, the minimal length for generating a transition to detonation. Further work is necessary to (i) analyse the case of donor explosives with finite reaction rates, and to (ii) account for the detonation cellular structure in the simulations of shock-to-detonation transitions. Received 21 December 2001 / Accepted 15 July 2002 Published online 4 November 2002 Correspondence to: Pierre Vidal An abridged version of this paper was presented at the 18th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Seattle, USA, from July 29 to August 3, 2001     

孔板扰动对爆轰波胞格结构特性影响的实验研究

实验采用稳定预混气2H₂+O₂+3Ar及不稳定预混气C₂H₂+5N₂O和CH₄+2O₂，在圆形爆轰管内通过烟膜手段记录了爆轰波的胞格结构，得到了胞格尺寸与初始压力之间的关系式；研究了胞格结构在扰动上下游的变化过程，分析了胞格不稳定性对胞格结构特征的影响，获得了爆轰波经过扰动后重新恢复至平衡状态的特征尺度。结果表明:爆轰波经过扰动后，对于稳定预混气，在扰动下游主胞格结构变得不规则，没有出现次生胞格；对于不稳定预混气，扰动下游伊始爆轰波的次生模态被抑制，由于爆轰波自身的不稳定性，随后出现了局部爆炸点及精细胞格结构；爆轰波在扰动下游传播了一段距离后恢复至平衡状态，该长度在8～15倍之间的胞格尺寸范围内变化，并且随初始压力的变化趋势并不明显。研究结果反映出爆轰波经过孔板扰动后恢复至平衡态所需的长度与爆轰波流体动力学厚度相当。     

Numerical and experimental study of a micro-blast wave generated by pulsed-laser beam focusing

This paper describes a numerical and experimental study of a micro-blast wave which is produced from the source of several tens microns in dia. and propagates in the length scale of a few centimeter in diameter. The micro-blast wave was generated by focusing a Nd:Glass pulsed-laser beam in ambient air. Its propagation and reflection were visualized by using double exposure holographic interferometry and simulated numerically using the dispersion-controlled scheme to solve the Euler and Navier-Stokes equations with initial conditions of a point-source explosion specified with the Taylor similarity law. Good agreement was obtained between numerical solutions and experimental results, and this spherical micro-blast wave was shown to be a handy model of blast waves created in large scale explosions. Received 28 October 1997 / Accepted 30 April 1998     

Behavior of detonation waves at low pressures

With respect to stability of gaseous detonations, unsteady behavior of galloping detonations and re-initiation process of hydrogen-oxygen mixtures are studied using a detonation tube of 14 m in length and 45 mm i.d. The arrival of the shock wave and the reaction front is detected individually by a double probe combining of a pressure and an ion probe. The experimental results show that there are two different types of the re-initiation mechanism. One is essentially the same as that of deflagration to detonation transition in the sense that a shock wave generated by flame acceleration causes a local explosion. From calculated values of ignition delay behind the shock wave decoupled from the reaction front, the other is found to be closely related with spontaneous ignition. In this case, the fundamental propagation mode shows a spinning detonation. Received 10 March 1997 / Accepted 8 June 1997     

Investigation of the detonation regimes in gaseous mixtures with suspended starch particles

The existence of a secondary discontinuity at the rear of a detonation front shown in experiments by Peraldi and Veyssiere (1986) in stoichiometric hydrogen-oxygen mixtures with suspended 20-m starch particles has not been explained satisfactorily. Recently Veyssiere et al. (1997) analyzed these results using a one-dimensional (1-D) numerical model, and concluded that the heat release rate provided by the burning of starch particles in gaseous detonation products is too weak to support a double-front detonation (DFD), in contrast to the case of hybrid mixtures of hydrogen-air with suspended aluminium particles in which a double-front detonation structure was observed by Veyssiere (1986). A two-dimensional (2-D) numerical model was used in the present work to investigate abovementioned experimental results for hybrid mixtures with starch particles. The formation and propagation of the detonation has been examined in the geometry similar to the experimental tube of Peraldi and Veyssiere (1986), which has an area change after 2 m of propagation from the ignition point from a 69 mm dia. section to a 53 mm 53 mm square cross section corresponding to a 33% area contraction. It is shown that the detonation propagation regime in these experiments has a different nature from the double-front detonation observed in hybrid mixtures with aluminium particles. The detonation propagates as a pseudo-gas detonation (PGD) because starch particles release their heat downstream of the CJ plane giving rise to a non-stationary compression wave. The discontinuity wave at the rear of the detonation front is due to the interaction of the leading detonation front with the tube contraction, and is detected at the farthest pressure gauge location because the tube length is insufficient for the perturbation generated by the tube contraction to decay. Thus, numerical simulations explain experimental observations made by Peraldi and Veyssiere (1986). Received 5 July 1997 / Accepted 13 July 1998     

Relaxation of overdriven detonation waves with finite reaction rate

A numerical study is made of the interaction of a detonation wave having finite reaction velocity with a rarefaction wave of different intensity which approaches it from the rear, for the Zeldovich-Neumann-Doring (ZND) model with a single irreversible reaction A B. It is found that, for a fixed value of the parameter characterizing the initial supercompression (depending on the activation energy and the heating value of the mixture), the considered interaction leads either to a gradual relaxation of the detonation wave and its transition to the Chapman-Jouguet (CJ) regime, or to the development of undamped oscillations.Interest in the problems of detonation and supersonic combustion has increased in recent years. This is associated with the appearance and development of new experimental and theoretical techniques; it is also associated with the further development of air-breathing reaction engines, and other practical requirements. The present state of detonation theory is reflected in the survey [1].It has been established [2] that the detonation wave in gases nearly always has a complex nonuniform structure. Transverse disturbances are observed under a wide range of conditions and differ both in amplitude and wavelength. At the same time, behind the detonation leading front there is a region of uncompletely burned gas corresponding to the effective ignition induction period [3]. In spinning detonation the induction period is significantly longer than the heat release period and transverse detonation waves traveling in the induction zone of the head wave appear [3, 4]. Such a secondary detonation wave is free of transverse disturbances. The same is true of the detonation waves observed in the wake behind a body moving at high speed in a combustible medium [5] or in a gas which has been preheated by a shock wave [6].Although it is possible, under favorable conditions, to study in detail the system of discontinuities accompanying detonation, information on the extensive zones in which heat release takes place is scarce, the mechanism of detonation wave autonomy (in particular, the role of the rarefaction zone behind the wave) is not entirely clear, and the fact that, in spite of the complex structure, an autonomous detonation propagates with the CJ velocity calculated on the basis of one-dimensional theory has not yet been explained.In studying the nonlinear phenomena associated with the finite reaction rate it is quite acceptable to investigate only the simple one-dimensional detonation model, with which it is convenient to restrict ourselves to a single effective chemical reaction. This model is particularly reasonable since, in certain cases, the real detonation is virtually one-dimensional.The question of the stability of the one-dimensional detonation wave to disturbances of its structure has been examined by several authors [7–13]. The use of computers makes possible the direct computation of flows with heat release and the study of their properties. This method has been used in [11–13] to study the stability problem for a detonation wave with respect to finite disturbances.In the present paper we present a numerical study of the interaction of a detonation wave having finite chemical reaction rate with a rarefaction wave of different intensity approaching it from the rear for the ZND model with a single irreversible reaction A B. It is found that for a fixed value of the parameter characterizing the difference between detonation and the CJ waves, depending on the activation energy E and the mixture heating value Q_m, the interaction in question leads either to a gradual relaxation of the detonation wave and its transition to the CJ regime (this relaxation may be accompanied by decaying oscillations) or to the appearance of undamped oscillations (the unstable regime). The parameters E and Q_m affect the wave stability differently: with increase of Q_m, the wave is stabilized; with increase of E, it is destabilized. The boundary between the stable and unstable detonation wave propagation regimes is found. This boundary has a weak dependence on the rarefaction wave intensity. Estimates and calculated examples show that the amplitude of the unstable wave oscillations is finite and that the average detonation propagation velocity is close to the CJ velocity computed for the given heating value Qm.The author wishes to thank G. G. Chernyi for his guidance and L. A. Chudov for advice on computational questions.     

驻定斜爆轰波并行数值模拟

采用多组分化学反应Euler方程组对驻定在高速飞行弹丸上的斜爆轰波流场进行了数值模拟。计算中分别采用TVD格式和基元反应模型,并基于并行编程模型MPI（message passing interface）实现了非结构网格上的并行计算,对流项和化学反应项用时间分裂法进行处理。计算结果表明并行计算能有效地提高计算速度,扩展计算规模,为进一步研究超驱爆轰推进技术奠定基础。     

Numerical prediction of rock mass damage due to accidental explosions in an underground ammunition storage chamber

Accidental detonations in an underground ammunition storage chamber inside a rock mass may cause severe damage to the rock mass around the chamber, adjacent tunnels and chambers, ground surface, and in the worst case cause sympathetic detonation of explosives in adjacent storage chambers. To prevent such damage, underground ammunition storage chambers are often situated at minimum depth below the ground surface, and spaced at minimum distance from each other, so that damage, should it occur, is limited to the accidental chamber. Different codes and regulations for ammunition storage chambers specify minimum embedment depth and separation distance for underground ammunition storage chambers. They are usually given in terms of the rock mass properties and the weight of explosive stored in chambers. Some empirical formulae, usually based on the peak particle velocity of the stress wave or the maximum strain of the rock mass, are also available to estimate the damage zones in the rock mass from an explosion. All these empirical methods do not include the effects of explosion details, such as the loading density, chamber geometry and explosive distribution. In this paper, a previously calibrated numerical model is used to estimate the damage zones in a granite mass resulting from an accidental explosion in an underground ammunition storage chamber. Effects of various explosion conditions on rock mass damage are investigated. On the basis of the numerical results, some empirical formulae are derived to predict damage zones around the explosion chamber, as well as safe embedment depth of the storage chamber and safe separation distance between adjacent chambers. The numerical results are also compared with available empirical formulae and code specifications. It should be noted that the characteristics of stress wave propagation around an ammunition storage chamber has been published in a preceding paper (Int. J. Blast. Fragm. 5:57–90, 2001.