爆轰波与激波对撞的实验研究

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波,以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.      

气体-燃料液滴两相系统爆轰的数值模拟

用两相流体力学模型对气体 燃料液滴系统进行了研究。数值模拟了点火后两相系统爆轰波的发展过程，得到爆轰波的结构和参数。数值模拟结果表明气体 燃料液滴系统爆轰波有较宽的反应区，因而两相爆轰波的曲率对爆速的影响效应十分明显。进行了燃料液滴尺寸对爆轰波的结构和参数的影响的数值模拟。除了很小的液滴外，燃料液滴在爆轰波前导激波面和ＣＪ面间不能完全气化。随着液滴尺寸的增加，燃料液滴在爆轰波前导激波面和ＣＪ面间释放出的能量随之减少，爆轰参数也随之下降。     

Interaction of a high-speed combustion front with a closely packed bed of spheres

The head-on collision of a combustion front with a closely packed bed of ceramic-oxide spheres was investigated in a vertical 76.2 mm diameter tube containing a nitrogen diluted stoichiometric ethylene–oxygen mixture. A layer of spherical beads in the diameter range of 3–12.7 mm was placed at the bottom of the tube and a flame was ignited at the top endplate. Four orifice plates spaced at one tube diameter were placed at the ignition end of the tube in order to accelerate the flame to either a “fast-flame” or a detonation wave before the bead layer face. The mixture reactivity was adjusted by varying the initial mixture pressure between 10 and 100 kPa absolute. The pressure before and within the bead layer was measured by flush wall-mounted pressure transducers. For initial pressures where a fast-flame interacts with the bead layer peak pressures recorded at the bead layer face were as high as five times the reflected Chapman–Jouget detonation pressure. The explosion resulting from the interaction developed by two distinct mechanisms; one due to the shock reflection off the bead layer face, and the other due to shock transmission and mixing of burned and unburned gas inside the bead layer. The measured explosion delay time (time after shock reflection from the bead layer face) was found to be independent of the incident shock velocity. As a result, the explosion initiation is not the direct result of the shock reflection process but instead is more likely due to the interaction of the reflected shock wave and the trailing flame. The bead layer was found to be very effective in attenuating the explosion front transmitted through the bead layer and thus isolating the tube endplate. This paper is based on work that was presented at the 21th International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

Modelling detonation waves in condensed energetic materials: multiphase CJ conditions and multidimensional computations

A hyperbolic multiphase flow model with a single pressure and a single velocity but several temperatures is proposed to deal with the detonation dynamics of condensed energetic materials. Temperature non-equilibrium effects are mandatory in order to deal with wave propagation (shocks, detonations) in heterogeneous mixtures. The model is obtained as the asymptotic limit of a total non-equilibrium multiphase flow model in the limit of stiff mechanical relaxation only (Kapila et al. in Phys Fluids 13:3002–3024, 2001). Special attention is given to mass transfer modelling, that is obtained on the basis of entropy production analysis in each phase and in the system (Saurel et al. in J Fluid Mech 607:313–350, 2008). With the help of the shock relations given in Saurel et al. (Shock Waves 16:209–232, 2007) the model is closed and provides a generalized ZND formulation for condensed energetic materials. In particular, generalized CJ conditions are obtained. They are based on a balance between the chemical reaction energy release and internal heat exchanges among phases. Moreover, the sound speed that appears at sonic surface corresponds to the one of Wood (A textbook of sound, G. Bell and Sons LTD, London, 1930) that presents a non-monotonic behaviour versus volume fraction. Therefore, non-conventional reaction zone structure is observed. When heat exchanges are absent, the conventional ZND model with conventional CJ conditions is recovered. When heat exchanges are involved interesting features are observed. The flow behaviour presents similarities with non ideal detonations (Wood and Kirkwood in J Chem Phys 22:1920–1924, 1950) and pathological detonations (Von Neuman in Theory of detonation waves, 1942; Guenoche et al. in AIAA Prog Astron Aeronaut 75: 387–407, 1981). It also present non-conventional behaviour with detonation velocity eventually greater than the CJ one. Multidimensional resolution of the corresponding model is then addressed. This poses serious difficulties related to the presence of material interfaces and shock propagation in multiphase mixtures. The first issue is solved by an extension of the method derived in Saurel et al. (J Comput Phys 228(5):1678–1712, 2009) in the presence of heat and mass transfers. The second issue poses the difficult mathematical question of numerical approximation of non-conservative systems in the presence of shocks associated to the physical question of energy partition among phases for a multiphase shock. A novel approach is used, based on extra evolution equations used to retain the information of the material initial state. This method insures convergence in the post-shock state. Thanks to these various theoretical and numerical ingredients, one-dimensional and multidimensional unsteady detonation waves computations are done, eventually in the presence of material interfaces. Convergence of the numerical hyperbolic solver against ZND multiphase solution is reached. Material interfaces, shocks, detonations are solved with a unified formulation where the same equations are solved everywhere with the same numerical scheme.     

Detonation Initiation by Shock Reflection from an Orifice Plate

Results from an experimental investigation of the interaction of a “non-ideal” shock wave and a single obstacle are reported. The shock wave is produced ahead of an accelerated flame in a 14 cm inner-diameter tube partially filled with orifice plates. The shock wave interacts with a single larger blockage orifice plate placed 15–45 cm after the last orifice plate in the flame acceleration section of the tube. Experiments were performed with stoichiometric ethylene–oxygen mixtures with varying amounts of nitrogen dilution at atmospheric pressure and temperature. The critical nitrogen dilution was found for detonation initiation. It is shown that detonation initiation occurs if the chemical induction time based on the reflected shock state is shorter than the time required for an acoustic wave to traverse the orifice plate upstream surface, from the inner to the outer diameter. The similarity between the present results and those obtained from previous investigators looking at detonation initiation by ideal shock reflection produced in a shock tube indicates that the phenomenon is not sensitive to the detailed structure of the shock front but only on the average shock strength.This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems, Montreal, Canada, July 31–August 5, 2005.     

On the propagation mechanism of a detonation wave in a round tube with orifice plates

This study deals with the investigation of the detonation propagation mechanism in a circular tube with orifice plates. Experiments were performed with hydrogen air in a 10-cm-inner-diameter tube with the second half of the tube filled with equally spaced orifice plates. A self-sustained Chapman–Jouguet (CJ) detonation wave was initiated in the smooth first half of the tube and transmitted into the orifice-plate-laden second half of the tube. The details of the propagation were obtained using the soot-foil technique. Two types of foils were used between obstacles, a wall-foil placed on the tube wall, and a flat-foil (sooted on both sides) placed horizontally across the diameter of the tube. When placed after the first orifice plate, the flat foil shows symmetric detonation wave diffraction and failure, while the wall foil shows re-initiation via multiple local hot spots created when the decoupled shock wave interacts with the tube wall. At the end of the tube, where the detonation propagated at an average velocity much lower than the theoretical CJ value, the detonation propagation is much more asymmetric with only a few hot spots on the tube wall leading to local detonation initiation. Consecutive foils also show that the detonation structure changes after each obstacle interaction. For a mixture near the detonation propagation limit, detonation re-initiation occurs at a single wall hot spot producing a patch of small detonation cells. The local overdriven detonation wave is short lived, but is sufficient to keep the global explosion front propagating. Results associated with the effect of orifice plate blockage and spacing on the detonation propagation mechanism are also presented.     

爆轰波透射孔栅形成的高速爆燃波的结构和行为

以实验为主,研究光滑直管中乙炔氧气爆轰波透射孔栅形成的高速爆燃波结构和行为。实验以基于纹影平台的高速转鼓摄影记录孔栅下游近场内初始爆燃波的结构和发展,并以压力传感器跟踪其后继走向。研究发现,高速爆燃波是前驱激波和火焰的组合结构,随着初始压力的提升,它分别表现为紧随于前驱激波的层流和湍流燃烧火焰。入射爆轰胞格尺度小于或与孔栅扰动尺度相当时,湍流燃烧在下游近场迅速形成;层流结构的爆燃波通常无法抵制背景稀疏波而走向衰弱,而湍流燃烧结构可发生加速和向爆轰的转捩;他们之间存在一个不稳定的临界状态,高速爆燃波得以以0.5~0.6倍CJ爆轰速度传播较长距离,这一状态对应于双间断Rankine-Hugoniot关系的等容燃烧解。     

Shock wave attenuation by grids and orifice plates

The interaction of weak shock waves with porous barriers of different geometries and porosities is examined. Installing a barrier inside the shock tube test section will cause the development of the following wave pattern upon a head-on collision between the incident shock wave and the barrier: a reflected shock from the barrier and a transmitted shock propagating towards the shock tube end wall. Once the transmitted shock wave reaches the end wall it is reflected back towards the barrier. This is the beginning of multiple reflections between the barrier and the end wall. This full cycle of shock reflections/interactions resulting from the incident shock wave collision with the barrier can be studied in a single shock tube test. A one-dimensional (1D), inviscid flow model was proposed for simulating the flow resulting from the initial collision of the incident shock wave with the barrier. Fairly good agreement is found between experimental findings and simulations based on a 1D flow model. Based on obtained numerical and experimental findings an optimal design procedure for shock wave attenuator is suggested. The suggested attenuator may ensure the safety of the shelter’s ventilation systems.     

Head on collision of a planar shock wave with a dusty gas layer

The head-on collision of a planar shock wave with a dust-air suspension is studied numerically. In this study the suspension is placed inside a conduit adjacent to its rigid end-wall. It is shown that as a result of this collision two different types of transmitted shock waves are possible, depending on the strength of the incident shock wave and the dust loading ratio in the suspension. One possibility is a partially dispersed shock wave, the other is a compression wave. The flow fields resulting in these two options are investigated. It is shown that in both cases, at late times after the head-on reflection of the transmitted shock wave from the conduit end-wall a negative flow (away from the end-wall) is evident. The observed flow behavior may suggest a kind of dust particle lifting mechanism that could shed new light on the complex phenomenon of dust entrainment behind sliding shock waves.      

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.     

Experimental study and kinetic modeling of the thermal decomposition of gaseous monomethylhydrazine. Application to detonation sensitivity

The thermal decomposition of gaseous monomethylhydrazine has been studied in a 38.4 mm i.d. shock tube behind a reflected shock wave at 1040–1370 K, 140–455 kPa and in mixtures containing 97 to 99 mol% argon, by using MMH absorption at 220 nm. A chemical kinetic model based on MMH decomposition profiles has been developed. This model has been used, with some assumptions, to evaluate the detonation sensitivity of pure gaseous MMH. This compound is found to be much less sensitive to detonation than hydrazine. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

Stability of Detonation Profiles in the ZND Limit

Confirming a conjecture of Lyng–Raoofi–Texier–Zumbrun, we show that stability of strong detonation waves in the ZND, or small-viscosity, limit is equivalent to stability of the limiting ZND detonation together with stability of the viscous profile associated with the component Neumann shock. Moreover, on bounded frequencies the nonstable eigenvalues of the viscous detonation wave converge to those of the limiting ZND detonation, while on frequencies of order one over viscosity, they converge to one over viscosity times those of the associated viscous Neumann shock. This yields immediately a number of examples of instability and Hopf bifurcation of reacting Navier–Stokes detonations through the extensive numerical studies of ZND stability in the detonation literature.     

The criterion of the existence or inexistence of transverse shock wave at wedge supported oblique detonation wave

A simplified theoretic method and numerical simulations were carried out to investigate the characterization of propagation of transverse shock wave at wedge supported oblique detonation wave.After solution validation,a criterion which is associated with the ratio Φ (u 2 /u CJ) of existence or inexistence of the transverse shock wave at the region of the primary triple was deduced systematically by 38 cases.It is observed that for abrupt oblique shock wave (OSW)/oblique detonation wave (ODW) transition,a transverse shock wave is generated at the region of the primary triple when Φ < 1,however,such a transverse shock wave does not take place for the smooth OSW/ODW transition when Φ > 1.The parameter Φ can be expressed as the Mach number behind the ODW front for stable CJ detonation.When 0.9 < Φ < 1.0,the reflected shock wave can pass across the contact discontinuity and interact with transverse waves which are originating from the ODW front.When 0.8 < Φ < 0.9,the reflected shock wave can not pass across the contact discontinuity and only reflects at the contact discontinuity.The condition (0.8 < Φ < 0.9) agrees well with the ratio (D ave /D CJ) in the critical detonation.     

A numerical investigation of transient detonation in granulated material

A two-phase model based upon principles of continuum mixture theory is numerically solved to predict the evolution of detonation in a granulated reactive material. Shock to detonation transition (SDT) is considered whereby combustion is initiated due to compression of the material by a moving piston. In particular, this study demonstrates the existence of a SDT event which gives rise to a steady two-phase Chapman-Jouguet (CJ) detonation structure consisting of a single lead shock in the gas and an unshocked solid; this structure has previously been independently predicted by a steady-state theory. The unsteady model equations, which constitute a non-strictly hyperbolic system, are numerically solved using a modern high-resolution method. The numerical method is based on Godunov's method, and utilizes an approximate solution for the two-phase Riemann problem. Comparisons are made between numerical predictions and known theoretical results for 1) an inert two-phase shock tube problem, 2) an inert compaction wave structure, and 3) a reactive two-phase detonation structure; in all cases, good agreement exists. Received 4 August 1995 / Accepted 17 February 1996     

高速弹丸诱导斜爆轰激波结构实验研究

斜爆轰推进系统在高超声速推进领域具有广阔的应用前景,其释热迅速、比冲高、燃烧室结构简单的优点吸引研究人员的持续关注.然而,斜爆轰的地面试验同时涉及到高速试验环境模拟、燃料与氧化剂混合、高温燃烧流场结构测量等技术难点,当前国内外系统的试验研究仍然十分有限,难以支撑斜爆轰发动机的研制.为了研究自持传播的斜爆轰激波结构与波面流动特性,基于爆轰驱动二级轻气炮开展了高速弹丸诱导斜爆轰实验研究,使用直径30 mm球头圆柱形弹丸发射进入充满氢/氧可燃混合气体的实验舱中以起爆斜爆轰波,并采用两种阴影技术对实验流动结构进行测量.实验中在不同速度、不同充气压力下观察到三种弹丸诱导激波结构,即激波诱导燃烧、弹丸起爆爆轰波和相对弹丸驻定的斜爆轰波,实验舱充气压力下降则会造成爆轰横波尺度增加与波面流动失稳.实验中,斜爆轰激波角与理论分析结果吻合较好,弹丸气动不稳定带来较大的弹丸攻角会对激波角测量带来一定偏差.通过对斜爆轰波波面法向传播速度的测量发现,随着远离弹丸,斜爆轰传播速度由弹丸飞行速度衰减至接近实验气体CJ速度,弹丸速度的降低会加速斜爆轰波传播速度的衰减.     

喷管对水下爆轰气泡形态和压力特征影响的实验研究

针对具有不同类型喷管的爆轰管在水下爆轰中形成的燃气射流问题,搭建了爆轰实验平台,研究了单次爆轰过程中尾部喷管对水下气泡形态与压力特征的影响。采用数字粒子图像测速技术对高速摄影机拍摄得到的气泡脉动图片进行流场可视化分析,得到各喷管工况下的气泡速度场。为了确认爆轰管内是否形成稳定爆轰波,并观察爆轰波在气液两相界面上的透反射特性,爆轰管尾部安装有2个动态压力传感器,同时在距离喷管一定距离处设置一个水下爆炸传感器,以监测水中传播的压力波。结果表明：扩张喷管工况下的气泡脉动过程与直喷管工况基本一致,但扩张喷管提高了燃气射流速度,气泡膨胀体积更大;因为燃气射流的持续性,收敛喷管工况下的气泡脉动过程具有明显差异,气泡膨胀体积较小,但气泡二次脉动时长相较于一次脉动时长衰减更小;扩张喷管提高了气泡脉动强度,扩张喷管工况下的气泡脉动压力与透射冲击波压力远大于直喷管工况下的气泡脉动压力与透射冲击波压力;收敛喷管工况下的气泡脉动压力与透射冲击波压力都较小,但收敛喷管燃气射流的持续性减缓了气泡脉动压力的衰减速度。相比于直喷管,扩张喷管工况下的气泡脉动时间、气泡脉动压力与透射冲击波压力都更大。收敛喷管工况下的气泡...     

爆轰波与激波对撞爆轰状态演化机理研究

对当量比氢氧混合气体中爆轰波与激波的正面对撞过程进行了二维数值研究. 采用了二阶精度NND差分格式与改进的两阶段化学反应模型,并以数值x-t纹影图以及烟迹图记录了对撞过程. 数值研究表明,透射爆轰波受到膨胀影响首先会衰减,甚至产生局部解耦现象;然后由于三波点的碰撞又能再次重新耦合. 在爆轰波对撞过程中,由于燃烧不均匀性而产生的弱横波对爆轰胞格的形成起着重要作用.      

Re-initiation phenomenon of gaseous detonation induced by shock reflection

The two-dimensional, time-dependent, reactive Navier–Stokes equations including the effects of viscosity, thermal conduction and molecular diffusion were solved to reveal the wave evolution and chemical dynamics involved in the re-initiation process. The computation was performed for hydrogen–oxygen–argon mixtures at the low initial pressure (8.00 kPa), using detailed chemical reaction model. The results show that, the decoupled leading shock reflects on the right wall of the vertical branch. High temperature and pressure behind the reflected shock induce the generation of hot spots and local explosion. Therefore, the re-initiation of gaseous detonation occurs. In the re-initiation area, there exist very high OH concentration and no H ₂ concentration. However, in front of reflected shock, there exist relatively high H ₂ concentration and no OH radicals. Additionally, the shock–flame interaction induces RM instability. This results in the fast mixing between hot reacted gas mixture and the relatively cold unreacted gas mixture and accelerates the chemical reactions. However, the shock–flame interaction contributes much less to the re-initiation, in contrast with shock reflection. The transition of leading shock from regular reflection to Mach reflection happens during the re-initiation. The computed evolution of wave structures involved in the re-initiation is qualitatively agreeable with that from the experimental schlieren images.      

气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明:气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。     

含铝炸药能量释放的简化模型

为了在水下爆炸效应中反映出非理想爆轰特性的影响,建立了含铝炸药非理想爆轰能量释放的简化模型。该模型以CJ爆轰理论和二次反应理论为基础,把含铝炸药化学反应划分为快速反应和慢速反应,以释放的化学能和慢反应速率常数作为非理想特征参数,并应用于一维数值模拟。计算结果与基本实验结果一致,冲击波峰值的计算误差不大于10％,衰减时间常数的误差小于5％,冲击波能与实验值也具有良好的一致性。简化模型合理地描述了含铝炸药非理想爆轰的主要过程及非理想特性,可应用于含铝炸药的设计和爆炸效应的分析。