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

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

非均匀炸药冲击起爆和起爆后的行为

 介绍并分析了Campbell 等人及其他作者研究非均匀炸药冲击起爆和起爆后行为所获得的实验结果,但不涉及其冲击起爆条件。足够强的冲击波进入非均匀炸药后,爆轰将瞬时（指不经过感应时间）且直接（指不经过其他过程,如爆燃）被引发;非均匀炸药起爆后,其中传播的自始至终是一个不断增长的爆轰波,直至发展为正常爆轰,整个过程都是爆轰的增长（新定义）过程。不存在由反应冲击波不断增长并转变为爆轰波的所谓向爆轰的增长。所谓向爆轰的增长,实际上是爆轰的增长（按新定义）的初期;Craig原定义的爆轰的增长,实际上是爆轰的增长（按新定义）的后期;而所谓反应冲击波,实际上是增长中的初期爆轰波。爆轰的增长（按新定义）是所有猛炸药的特性,炸药反应不充分并逐渐趋于充分是爆轰的增长的化学机制。     

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

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

Detonation simulations in supersonic flow under circumstances of injection and mixing

The unsteady, reactive Navier-Stokes equations with a detailed chemical mechanism of 11 species and 27 steps were employed to simulate the mixing, flame acceleration and deflagration-to-detonation transition (DDT) triggered by transverse jet obstacles. Results show that multiple transverse jet obstacles ejecting into the chamber can be used to activate DDT. But the occurrence of DDT is tremendously difficult in a non-uniform supersonic mixture so that it required several groups of transverse jets with increasing stagnation pressure. The jets introduce flow turbulence and produce oblique and bow shock waves even in an inhomogeneous supersonic mixture. The DDT is enhanced by multiple explosion points that are generated by the intense shock wave focusing of the leading flame front. It is found that the partial detonation front decouples into shock and flame, which is mainly caused by the fuel deficiency, nevertheless the decoupled shock wave is strong enough to reignite the mixture to detonation conditions. The resulting transverse wave leads to further mixing and burning of the downstream non-equilibrium chemical reaction, resulting in a high combustion temperature and intense flow instabilities. Additionally, the longitudinal and transverse gradients of the non-uniform supersonic mixture induce highly dynamic behaviors with sudden propagation speed increase and detonation front instabilities.     

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

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

气相爆轰波在分叉管中传播现象的数值研究

数值研究气相爆轰波在分叉管中的传播现象.用二阶附加半隐龙格-库塔法和5阶WENO格式求解二维欧拉方程,用基元反应描述爆轰化学反应过程,得到了密度、压力、温度、典型组元质量分数场及数值胞格结构和爆轰波平均速度.结果表明:气相爆轰波在分叉管中传播,分叉口左尖点的稀疏波导致诱导激波后压力、温度急剧下降,诱导激波和化学反应区分离,爆轰波衰减为爆燃波(即爆轰熄灭).分离后的诱导激波在垂直支管右壁面反射,并导致二次起爆.畸变的诱导激波在水平和垂直支管中均发生马赫反射.分叉口上游均匀胞格区和分叉口附近大胞格区的边界不是直线,其起点通常位于分叉口左尖点上游或恰在左尖点.水平支管中马赫反射三波点迹线始于右尖点下游.分叉口左尖点附近的流场中出现了复杂的旋涡结构、未反应区及激波与旋涡作用.旋涡加速了未反应区的化学反应速率.反射激波与旋涡作用并使旋涡破碎.反射激波与未反应区作用,加速其反应消耗,并形成一个内嵌的射流.数值计算得到的波系演变和胞格结构与实验定性一致.     

Review of the Forest fire model for high explosives

Forest fire is a burn model for reactive hydro simulations that is widely used to describe both initiation and propagation of detonation waves in solid high explosives. Here we thoroughly review the assumptions of the model. In analogy with Whitham–Chisnell shock dynamics, we derive the Forest fire rate based on a characteristic analysis of a shock-to-detonation transition. Code issues with implementing the Forest fire model are also discussed. The analysis leading to the Forest fire rate is more general than the simplifying assumptions used by the model. The Forest fire analysis provides important constraints needed for any burn model to fit shock initiation data. This data is fundamental to ignition sensitivity of high explosives.     

Aluminum oxidation in shock and detonation waves

Experimental data on the detonation velocity of aluminized explosives and the temperature of the material behind the shock wave front in condensed media, including aluminum-oxidizer mixtures were examined. It was demonstrated that the oxidation of aluminum to the highest oxide behind the front of shock and detonation waves is limited by the dissociation of aluminum oxides at temperatures above 3.5 kK.     

Calculation of shock wave propagation in water containing reactive gas bubbles

The entry of a shock wave from air into water containing reactive gas (stoichiometric acetylene–oxygen mixture) bubbles uniformly distributed over the volume of the liquid has been numerically investigated using equations describing two-phase compressible viscous reactive flow. It has been demonstrated that a steady-state supersonic self-sustaining reaction front with rapid and complete fuel burnout in the leading shock wave can propagate in this bubbly medium. This reaction front can be treated as a detonation-like front or “bubble detonation.” The calculated and measured velocities of the bubble detonation wave have been compared at initial gas volume fraction of 2 to 6%. The observed and calculated data are in satisfactory qualitative and quantitative agreement. The structure of the bubble detonation wave has been numerically studied. In this wave, the gas volume fraction behind the leading front is approximately 3–4 times higher than in the pressure wave that propagates in water with air bubbles when the other initial conditions are the same. The bubble detonation wave can form after the penetration of the shock wave to a small depth (~300 mm) into the column of the bubbly medium. The model suggested here can be used to find optimum conditions for maximizing the efficiency of momentum transfer from the pressure wave to the bubbly medium in promising hydrojet pulse detonation engines.     

The SURF model and the curvature effect for PBX 9502

SURF is a high explosive burn model based on the ignition & growth concept of hot-spot reaction. For the TATB based explosive PBX 9502, the model has been calibrated to shock-to-detonation transition experiments. To apply the SURF model for propagating detonation waves, the rate has to be extended to a higher pressure regime than is sampled by shock initiation experiments. The experimentally measured curvature effect – detonation speed as a function of front curvature or D _n(κ) – provides the appropriate data for calibrating the propagation regime. The calibration to the curvature effect is based on the ODEs for the reaction zone profile of a detonation wave in conjunction with a shooting algorithm to determine the rate model parameters, for a given κ, needed to obtain a specified detonation speed. A complication for calibrating PBX 9502 rate models arises from the kink in the experimentally measured D _n(κ) curve. This results from the fast and slow reactions that TATB exhibits. To account for this, we use an extension of the SURF model that utilises a sequence of two reactions. The first, with a fast rate, is due to molecular decomposition and is described by the original SURF formulation. The second, with a slow rate, is due to carbon clustering and is used to contribute additional energy from the formation of carbon bonds. The wave profile equations are generalised to the SURF-plus model. Model parameters are then determined for the propagation regime to fit the curvature effect data. The extended model is applicable to both the shock initiation regime and the propagating detonation wave regime.     

悬浮RDX炸药粉尘爆轰的数值模拟

用两相流模型对悬浮RDX炸药粉尘爆轰波进行了数值模拟。RDX炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,颗粒表面发生熔化。参考液滴在高速气流作用下剥离的效应,假设炸药熔化部分在高速气流的作用下发生剥离,破碎成极小的颗粒,瞬时发生分解反应,释放出能量支持爆轰波传播。数值模拟了在不同粒径和浓度的悬浮RDX炸药粉尘中爆轰波的发展与传播过程,得到了爆轰波流场中气-固两相的物理量分布,并确定了爆轰波参数。在较低的RDX粉尘浓度条件下,爆轰波阵面压力的峰值曲线出现振荡。当RDX粉尘浓度在80~150 g/m3时,数值模拟得到的爆轰波阵面压力峰值曲线的振荡是规则的;当RDX粉尘浓度为70 g/m3时,爆轰波阵面压力峰值曲线出现不规则振荡。     

CE/SE方法数值模拟炸药粉尘爆轰

采用CE/SE方法数值模拟悬浮在空气中的RDX炸药粉尘的两相爆轰过程.炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,释放能量支持爆轰波传播.数值模拟爆轰波管中的粉尘爆轰,得到爆轰波流场中的物理量分布,确定爆轰参数,数值结果与文献符合较好.数值模拟复杂通道中的炸药粉尘爆轰,预测了爆轰波的发展和传播过程以及爆轰波后的流场演化.数值结果表明CE/SE方法能成功模拟气体-固体两相爆轰,为粉尘爆轰的研究提供了新的数值预测手段.     

基于Mie-Grüneisen状态方程的爆轰波运动数值模拟

依据C-J（Chapman-Jouguet）理论,对爆轰问题中的气态爆轰产物和未反应炸药分别考虑不同的参考状态,并根据参考状态选用特定的Mie-Grüneisen状态方程。忽略化学反应过程,爆轰产物厚度为零的前导激波面以界面的形式存在。数值模拟中,爆轰波的演化分为波面传播以及与未反应介质相互作用两个部分。传播过程中,爆轰波的传播速度即恒定的爆速,爆轰产物在传播过程中瞬间形成,而相互作用过程则是通过Mie-Grüneisen多介质混合模型来计算爆轰波的持续冲击作用。借助于Mie-Grüneisen状态方程以及Mie-Grüneisen多介质混合模型,可以很好地模拟爆轰波的运动过程。对比理论参数及文献的计算结果发现,模拟结果具备较好的准确度。     

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.     

氢化镁储氢型乳化炸药的爆炸特性研究

 通过理论计算和水下爆炸实验,初步研究了MgH₂敏化储氢型乳化炸药的爆炸特性和爆轰反应机理。结果表明：与玻璃微球敏化的乳化炸药相比,MgH₂敏化的乳化炸药水下爆炸的冲击波超压、比冲量、比冲击波能、比气泡能及水下爆炸比总能量显著增加,其中冲击波超压和水下爆炸总能量分别增加了20.5%和31.0%。MgH₂储氢型乳化炸药的爆轰机理与玻璃微球敏化乳化炸药不同,MgH₂在乳化炸药中起到了敏化剂和含能材料的双重作用,即MgH₂在乳化基质中水解产生均匀分布的氢气泡,起到了敏化作用,同时氢气参与爆炸反应,提高了炸药的爆炸能量和做功能力。     

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.     

隔板深度对激波聚焦起爆影响的数值模拟研究

考虑几何结构参数对激波聚焦触发爆轰波的复杂影响,对H2/Air预混气的环形射流激波聚焦起爆现象开展了数值模拟研究,详细分析了不同隔板深度条件下的激波聚焦过程、流场演化特征以及爆轰波参数变化规律。研究结果表明,凹腔内激波聚焦诱导的局部爆炸以及隔板前缘处射流形成"卷吸涡"是引起爆轰波触发的两个重要机制,而隔板深度是影响环形射流激波聚焦起爆性能的关键因素。随着隔板深度的增加,凹腔内激波聚焦的强度逐步增强,回传的能量损失有所减小,进而导致爆燃转爆轰的距离与时间显著缩短。此外,当隔板深度由1 mm逐渐增加至3 mm时,爆轰波自持传播稳定性呈现出先降低后升高的变化趋势,产生这一现象的主要原因是爆轰波强度与三波点运动的相互作用。     

一维过驱动爆轰波形成的数值研究

利用有限速率的基元反应模型对一维过驱动爆轰波的形成过程进行数值模拟.研究表明,在上游高温、高压、高速来流作用下首先会形成一道强激波,其波面方存在诱导区和放热区,然后诱导区和放热区界面会在来流扰动的作用下发生失稳,经过复杂的波系运动过程形成过驱动爆轰波.通过对不同初始条件下界面失稳过程的模拟和分析,研究了混合气体的组元、温度,来流的压力、温度、速度对过驱动爆轰波形成的影响.     

Effects of high activation energies on acoustic timescale detonation initiation

Acoustic timescale Deflagration-to-Detonation Transition (DDT) has been shown to occur through the generation of compression waves emitted by a hot spot or reaction centre where the pressure and temperature increase with little diminution of density. In order to compensate for the multi-scale nature of the physico-chemical processes, previous numerical simulations in this area have been limited to relatively small activation energies. In this work, a computational study investigates the effect of increased activation energy on the time required to form a detonation wave and the change in behaviour of each hot spot as the activation energy is increased. The simulations use a localised spatially distributed thermal power deposition of limited duration into a finite volume of reactive gas to facilitate DDT. The Adaptive Wavelet-Collocation Method is used to solve efficiently the 1-D reactive Euler equations with one-step Arrhenius kinetics. The DDT process as described in previous work is characterised by the formation of hot spots during an initial transient period, explosion of the hot spots and creation of an accelerating reaction front that reaches the lead shock and forms an overdriven detonation wave. Current results indicate that as the activation energy is raised the chemical heat release becomes more temporally distributed. Hot spots that produce an accelerating reaction front with low activation energies change behaviour with increased activation energy so that no accelerating reaction front is created. An acoustic timescale ratio is defined that characterises the change in behaviour of each hot spot.     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟