The influence of multi-layer confinement on detonation propagation in condensed-phase explosives

We examine, via multi-material simulation in a two-dimensional planar geometry, the effects on steady detonation propagation of the presence of a low-density intermediate layer between a condensed-phase high explosive (HE) and a high-density metallic confiner of finite thickness. Such elastomer intermediate layers are often added to eliminate air-gaps and the associated jetting effects that can arise due to machining imprecisions, or to prevent HE cracking due to environmental changes. Without an intermediate layer, the flow structure of a steady detonation/metal confiner interaction is well understood. In particular, there is no reflected wave passed into the HE due to the metal confinement. With the elastomer layer present, we find that, as the intermediate layer width increases, a complex wave interaction and communication path develops between the HE, intermediate, and metal layers. For thin intermediate layers, a shock-driven subsonic flow develops in the intermediate layer, passing information from the metal layer to the HE, with the detonation speed decreasing as the intermediate layer width increases. For wider intermediate layers, a Mach stem configuration develops in the intermediate layer, forcing a shock to be reflected into the HE. Simultaneously, a localized Prandtl-Meyer fan emerges from the intersection of the detonation shock with the HE-intermediate layer material interface. These HE structures are shown to have a substantial effect on the structure of the detonation driving zone. The Prandtl-Meyer fan becomes the dominant structure for critically large intermediate layer widths, wherein the presence of the metal layer does not affect the detonation propagation. We examine the detonation propagation speed and reaction and driving zone structure as a function of varying intermediate layer width. Two confinement metals are examined, along with two high explosive and three metal layer widths.     

Experimental measurement of energy release from an initiating layer in an insensitive explosive

When subjected to a shock of insufficient strength to trigger prompt reaction, heterogenous condensed phase explosives can form regions where significant amounts of the explosive remain unreacted for times much greater than the reaction time of the detonating explosive. This phenomena is observed for the explosive PBX 9502 (95 wt% TATB) both for planar and oblique input shocks. In this work, we build on previous results by performing cylinder expansion (CYLEX) tests where the explosive charge is comprised of a faster core of PBX 9501 (95 wt% HMX) inside a slower annulus of PBX 9502. The detonation in the faster PBX 9501 drives an oblique shock into the adjacent PBX 9502, and an annular transverse initiating layer (IL) results. In the test geometry, the IL travels steadily down the length of the test after a short run distance. At radial positions beyond the IL, an annular region of detonating PBX 9502 is observed. Using standard CYLEX test diagnostics, we infer the total energy release of this experiment. By making the assumptions that (1) the combined energy release is comprised of contributions from detonating PBX 9501, detonating PBX 9502, and the IL in the PBX 9502 and (2) mass-specific energy release for the detonating explosives is approximately the same as typically observed for each explosive, the IL energy release and reaction efficiency can be computed. Results are compared to prior results for a similar geometry, and indicate that while shock deadened PBX 9502 does not detonate promptly, it does eventually release a significant portion of its chemical potential energy over longer timescales on the order of 10 µs.     

Peel-off case failure in thermal explosions observed by the deflagration cylinder test

We have discovered a previously unidentified thermal explosion mode using the Los Alamos deflagration cylinder test (DFCT). The DFCT is a “pipe bomb”-style test similar to the detonation cylinder test (DTCT), which has been used for many years to calibrate detonation product equations of state. The shot is heated in an oven to a uniform test temperature. The pre-heated high explosive (HE) is triggered by a hot wire initiator on one end. The tube is back-illuminated by a bright light source, and its combustion-driven deformation and subsequent break-up are observed by a high-speed framing camera. Like the DTCT, the DFCT tube wall motion provides the primary diagnostic. A variety of reactive responses are possible, including quasi-steady deflagration and deflagration-to-detonation transition. This paper focuses on the behavior of the HMX-based explosive PBX 9501 at 155 °C. Under this condition burning appeared to occur only at the HE/tube interface, causing the tube to peel away from the HE core. Peel-off propagated as a wave that traveled along the tube at 500 m/s. This failure mode resulted in vigorous case venting, but the response was otherwise benign. We derive a steady peel-off-wave model that reproduces the essential observed features for realistic PBX 9501 parameter values.     

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.     

180°圆弧形钝感炸药的爆轰波传播

 采用贴体坐标下与Level Set方法相结合的爆轰冲击波动力学（DSD）计算方法,研究了180°圆弧形钝感炸药中非理想爆轰波的传播过程。通过数值模拟计算和实验测量的对比分析,得到了180°圆弧形炸药中爆轰波传播的一些规律：圆弧形钝感炸药可以实现定常爆轰,即在极坐标中整个爆轰波以固定角速度转动。这种定常阵面的形状和角速度与圆弧的外半径无关,定常体系依赖于圆弧形炸药的内半径和覆盖圆弧的外壳物质。对描述圆弧形炸药中爆轰波传播规律的经验公式进行了研究,结果表明这些经验公式能够准确描述爆轰波速度的变化,在实验测量和预估方面具有一定的参考价值。     

固体炸药的磁驱动准等熵压缩实验研究

较宽压力范围内未反应炸药的本构关系和状态方程对于深入和精确认识压缩波作用下炸药组分间相互作用的力学过程和起爆热点的形成机制具有重要意义。较之冲击压缩,磁驱动准等熵压缩加载(无冲击压缩)是获取较宽压力范围内未反应炸药的动态压缩力学特性更有效的手段。基于大电流产生的电磁力作用原理,在国内率先实现了炸药的磁驱动无冲击压缩实验技术。通过对负载电极、炸药样品参数的优化设计和安装工艺的控制,实现了5 GPa载荷内JO-9159炸药的磁驱动准等熵压缩加载。基于激光位移干涉测量技术和Lagrange数据处理方法,获得了JO-9159炸药的速度响应历史和准等熵压缩线。所得结果与文献数据进行了比较,结果表明,在实验压力范围内,JO-9159炸药的等熵压缩线与PBX9501炸药的等熵压缩线一致。     

弹粘塑性双球壳塌缩热点反应模型

 基于Kim的弹粘塑性单球壳塌缩模型,考虑PBX炸药中的粘结剂效应,假设炸药和粘结剂均为弹粘塑性材料,建立了弹粘塑性双球壳塌缩热点反应模型,给出了炸药球壳在冲击压力作用下的速度、应变、温度和化学反应速率的时空分布,以及新的热点反应速率理论表达式。把新的热点反应项与Kim的低压下慢反应项和张震宇提出的高压反应速率方程相结合,得到了新的冲击起爆三项式细观反应速率模型。把该模型加入DYNA2D中,模拟了PBX-9501炸药的一维冲击起爆过程,结果表明：该模型除了可以解释炸药颗粒度和孔隙度的影响外,还可以较好地描述粘结剂强度和含量对PBX炸药冲击起爆感度的影响。     

Reaction propagation modes in millimeter-scale tubes for ethylene/oxygen mixtures

Effects of tube diameter and equivalence ratio on reaction front propagations of ethylene/oxygen mixtures in capillary tubes were experimentally analyzed using high speed cinematography. The inner diameters of the tubes investigated were 0.5, 1, 2 and 3 mm. The flame was ignited at the center of the 1.5 m long smooth tube under ambient pressure and temperature before propagated towards the exits in the opposite directions. A total of five reaction propagation scenarios, including deflagration-to-detonation transition followed by steady detonation wave transmission (DDT/C–J detonation), oscillating flame, steady deflagration, galloping detonation and quenching flame, were identified. DDT/C–J detonation mode was observed for all tubes for equivalence ratios in the vicinity of stoichiometry. The velocity for the steady detonation wave propagation was approximately Chapman–Jouguet velocity for 1, 2, and 3 mm I.D. tubes; however, a velocity deficit of 5% was found for the case in 0.5 mm I.D. tube. For leaner mixtures, an oscillating flame mode was found for tubes with diameters of 1 to 3 mm, and the reaction front travelled in a steady deflagrative flame mode with velocities around 2–3 m/s when the mixture equivalence ratio becomes even leaner. Galloping detonation wave propagation was the dominant mode for the fuel lean regime in the 0.5 mm I.D. tube. For rich mixtures beyond the detonation limits, a fast flame followed by flame quenching was observed.     

All Eulerian method of computing elastic response of explosively pressurised metal tube

We present an all Eulerian approach to simulate the elastic response of a metal tube loaded explosively by a gaseous detonation. The high strain rate deformation of the metal tube subjected to high explosive detonation is mathematically described by hyperbolic processes where the characteristics of existing wave motions were correlated with the local particle velocities through the speed of sound in the metal. This is a favourable case for the hydrocode which is based on a compressible gas dynamics solver and for simulating a high strain rate and dominantly plastic response of a material subject to an explosive loading. The hydrocodes fall substantially short of predicting elastic motion without the plastic flow of the confining material, for relatively minor pressure loadings due to a gaseous explosion as opposed to a high explosive detonation of a charged tube. The corresponding loading pressure due to gaseous explosion is a few orders of magnitude lower than those resulting from high explosive loadings. Utilising a hydrocode designed to handle the reactive process leading to a plastic flow of the confining materials is of great interest and a significant challenge. The new technique, based on the Eulerian framework, preserves the feature of a Lagrangian code while utilising all the benefits of an Eulerian solver that uses fixed grids with the level-sets for defining the multi-material interfaces. The hybrid particle level-set algorithm is combined with a hydrodynamic solver that adds an elasticity correction when handling the structural response while the overall scheme remained hyperbolic during the entire reactive flow. Several unseen dynamics of detonation flow associated with the elastically loaded tube of finite thickness are reported by using the present method for analysing the highly pressurised vessel.     

激光作用铝靶爆轰波流场观测与理论分析

进行了强激光作用铝靶实验,采用纹影照相法,观察爆轰波流场演化过程,分析了爆轰波衰减规律。根据冲击波运动的自相似性,采用点爆炸模型描述了激光作用下爆轰波流场的演化,计算了波阵面速度、压力和温度。结果显示:爆轰波阵面沿迎着激光光源方向较快传播,波阵面形状由最初的半椭球形逐渐向半球形转变,在演化过程中扰动区结构复杂,存在多个密度间断层。在爆轰波开始传播阶段,波阵面的压力和温度较高但下降很快。     

Flame spread across surfaces of PBX 9501

There is little flame spread data for homogeneous energetic materials and no data for nitramines. We report the results of flame spread experiments of PBX 9501 (HMX (cyclotetramethylenetetranitramine) based explosive). The horizontal flame spread rate, S_f, is of the same order of magnitude as normal deflagration and varies nearly as the square root of pressure, as our scaling analysis presented here predicts. In the vertical orientation, the flame propagation downward was observed to be slightly faster than horizontal flame spread, presumably because of the melt layer flowing downward on the sample. In an accident scenario, a charge may be fractured or the surface roughened. Consequently, we also examined the effect of roughness. Minor roughness created by explosives machining was found to have a negligible affect on flame spread. However, more significant roughness can increase the rate between two and three times over normal flame spread for the conditions considered here. In addition we examine the effect of sample edges and configuration. Corners result in more favorable heat loss and therefore affect flame spread rate. We argue that the increased spread rate on edges and rough surfaces is because of favorable heat transfer convergence.     

On the structure and accuracy of programmed burn

Accurate computation of the evolution of a (typically curved) detonation front in a complex geometry, and of the state behind it, is a practical problem in the design of devices that use high-energy explosives. Direct numerical simulations are infeasible: accuracy demands high resolution of the smallest scale (the reaction zone), which is typically several orders of magnitude smaller than the device scale. Programmed burn is an engineering alternative that has been shown to produce acceptable results at a fraction of the cost. The underlying algorithm prescribes the trajectory of the detonation front a priori and replaces the actual reaction zone by a mock up that is a few computational cells thick and in which the reaction rate is taken to be a constant. The state of the explosive at the end of the reaction zone is thereby computed at a relatively modest cost, and the bulk of the computational effort is reserved for the region behind the reaction zone wherein the products of detonation perform useful work. The reasons for the remarkable fidelity to which the physical situation is captured by the programmed burn are not well-understood. This investigation, aimed at achieving such an understanding, considers a model problem for a steady, curved detonation propagating down a rate stick. It examines the pseudo-reaction-zone structure of the programmed burn, studies the sensitivity of the state of the reaction products to the choice of the reaction zone length, and compares the results to those for the actual, physical reaction zone. Conclusions are drawn as to the causes behind the success of the programmed-burn algorithm. The analysis is based on the asymptotic limits of small front curvature and small departures from the Chapman–Jouguet speed. Results are presented for ideal as well as nonideal explosives.     

高能炸药爆速曲率关系的一种计算方法

研究平面爆轰波爆轰化学能和状态方程之间的关系,得到平面爆轰波关系式和CJ条件.基于柱坐标系下的一维流体力学方程组,采用坐标变换方法得到含曲率的爆轰波关系式及广义声速条件.在反应区内采用压力平衡和温度平衡的混合规则,在给定的复杂状态方程和复杂反应率下求解,得到PBX9404和PBX9502炸药的爆速和曲率的关系.已知炸药...     

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

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

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射

 采用高速转镜分幅相机和电探针技术研究了猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射图像和不同位置上的爆轰波传播时间。从研究得出，两种炸药都在拐角顶点附近绕射，爆轰波传播时间增长，爆速变小。但是两种炸药绕射爆轰波的状态不一样，钝感炸药IHE-2中爆轰波绕过直角时，在拐角顶点附近约10 mm范围内炸药未完全反应，猛炸药RHT-901中爆轰波绕过直角时未出现类似现象。两者相比，钝感炸药中绕射爆轰波速度变化大，波阵面曲率半径小，而猛炸药的绕射爆轰波速度变化小，波阵面曲率半径大。这说明炸药的爆轰波绕射与炸药的冲击感度、反应区宽度有关。     

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

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

拉氏程序使用减敏模型模拟炸药绕爆问题

利用二维拉氏非结构网格爆轰流体程序对高能炸药的爆轰波绕爆问题进行数值模拟.拉氏程序对绕爆问题进行模拟时,会引起网格严重扭曲,导致网格大变形.使用JWL状态方程和三项式反应率,拉氏程序得到的绕爆数值结果与实验定性相同,较好地模拟了爆轰波绕爆过程中的漩涡现象.为了模拟一些钝感炸药在绕爆过程中的死区现象,在炸药反应率中加入减敏因子,模拟结果捕捉到死区,并与国外文章中的数值模拟结果定性相同.     

散心爆轰的数值模拟

本文在瞬时稳态爆轰的假定下,简要叙述了散心爆轰的不定常性质。在散心爆轰的数值模拟方法中,我们改进了人造粘性并用Cochran反应率和JWL状态方程来描述炸药的引爆过程。计了几个典型的散心爆轰问题,得到了比较满意的结果。     

PBX-9502炸药爆轰反应速率模型研究

利用高能炸药爆轰反应的拉格朗日分析方法,得到了一个形式简单的PBX-9502炸药的爆轰反应速率方程。运用该速率方程,计算了PBX-9502炸药爆轰波反应区的纵向结构,并与三项式点火增长模型的计算结果及实验结果进行比较。结果表明:该爆轰反应速率方程能很好地模拟PBX-9502炸药爆轰波反应区的纵向结构。