A Reactive Burn Model for Shock Initiation in a PBX


Book Description

In the formulation of a reactive burn model for shock initiation, we endeavor to incorporate a number of effects based on the underlying physical concept of hot spot ignition followed by the growth of reaction due to diverging deflagration fronts. The passage of a shock front sets the initial condition for reaction, leading to a fraction of the hot spots that completely burn while others will quench. The form of the rate model is chosen to incorporate approximations based on the physical picture. In particular, the approximations imply scaling relations that are then used to mathematically separate various contributions. That is, the model is modular and refinements can be applied separately without changing the other contributions. For example, the effect of initial temperature, porosity, etc. predominantly enter the characterization of the non-quenching hot spot distribution. A large collection of velocity gauge data is shown to be well represented by the model with a very small number of parameters.




Modeling the Shock Initiation of PBX 9501 in ALE3D.


Book Description

The SMIS (Specific Munitions Impact Scenario) experimental series performed at Los Alamos National Laboratory has determined the 3-dimensional shock initiation behavior of the HMX based heterogeneous high explosive, PBX9501, which has a PMMA case and a steel impact cover. The SMIS real-world shot scenario creates a unique test-bed because many of the fragments arrive at the impact plate off-center and at an angle of impact. The goal of this model validation experiments is to demonstrate the predictive capability of the Tarver-Lee Ignition and Growth (I & G) reactive flow model in this fully 3-dimensional regime of Shock to Detonation Transition (SDT).




Deflagration Wave Profiles


Book Description

Shock initiation in a plastic-bonded explosives (PBX) is due to hot spots. Current reactive burn models are based, at least heuristically, on the ignition and growth concept. The ignition phase occurs when a small localized region of high temperature (or hot spot) burns on a fast time scale. This is followed by a growth phase in which a reactive front spreads out from the hot spot. Propagating reactive fronts are deflagration waves. A key question is the deflagration speed in a PBX compressed and heated by a shock wave that generated the hot spot. Here, the ODEs for a steady deflagration wave profile in a compressible fluid are derived, along with the needed thermodynamic quantities of realistic equations of state corresponding to the reactants and products of a PBX. The properties of the wave profile equations are analyzed and an algorithm is derived for computing the deflagration speed. As an illustrative example, the algorithm is applied to compute the deflagration speed in shock compressed PBX 9501 as a function of shock pressure. The calculated deflagration speed, even at the CJ pressure, is low compared to the detonation speed. The implication of this are briefly discussed.




Modeling Three-Dimensional Shock Initiation of PBX 9501 in ALE3D.


Book Description

A recent SMIS (Specific Munitions Impact Scenario) experimental series performed at Los Alamos National Laboratory has provided 3-dimensional shock initiation behavior of the HMX-based heterogeneous high explosive, PBX 9501. A series of finite element impact calculations have been performed in the ALE3D [1] hydrodynamic code and compared to the SMIS results to validate and study code predictions. These SMIS tests used a powder gun to shoot scaled NATO standard fragments into a cylinder of PBX 9501, which has a PMMA case and a steel impact cover. This SMIS real-world shot scenario creates a unique test-bed because (1) SMIS tests facilitate the investigation of 3D Shock to Detonation Transition (SDT) within the context of a considerable suite of diagnostics, and (2) many of the fragments arrive at the impact plate off-center and at an angle of impact. A particular goal of these model validation experiments is to demonstrate the predictive capability of the ALE3D implementation of the Tarver-Lee Ignition and Growth reactive flow model [2] within a fully 3-dimensional regime of SDT. The 3-dimensional Arbitrary Lagrange Eulerian (ALE) hydrodynamic model in ALE3D applies the Ignition and Growth (I & G) reactive flow model with PBX 9501 parameters derived from historical 1-dimensional experimental data. The model includes the off-center and angle of impact variations seen in the experiments. Qualitatively, the ALE3D I & G calculations reproduce observed 'Go/No-Go' 3D Shock to Detonation Transition (SDT) reaction in the explosive, as well as the case expansion recorded by a high-speed optical camera. Quantitatively, the calculations show good agreement with the shock time of arrival at internal and external diagnostic pins. This exercise demonstrates the utility of the Ignition and Growth model applied for the response of heterogeneous high explosives in the SDT regime.




Numerical Calculation of Shock-induced Initiation of Detonations. [PBX 9501].


Book Description

The results of some numerical calculations of the impact of steel cylinders and spheres on the plastic-bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two-dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor-acceptor-cell method to compute mixed cell fluxes. The mechanism of shock initiation to detonation in heterogeneous explosives is best described as local decomposition at hot spots that are formed by shock interactions with density discontinuities. The liberated energy strengthens the shock so that as it interacts with additional inhomogeneities, hotter hot spots are formed, and more of the explosive is decomposed. The shock wave grows stronger until a detonation begins. This mechanism of initiation has been described numerically by the Forest Fire burn model, which gives the rate of explosive decomposition as a function of local pressure. The parameters in the Forest Fire burn model have been developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results have been compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying sizes, shapes, and velocities. We find that there is good agreement between the calculated and experimental data.







Testing and Modeling of PBX-9591 Shock Initiation


Book Description

This paper describes an ongoing effort to develop a detonation sensitivity test for PBX-9501 that is suitable for studying pristine and damaged HE. The approach involves testing and comparing the sensitivities of HE pressed to various densities and those of pre-damaged samples with similar porosities. The ultimate objectives are to understand the response of pre-damaged HE to shock impacts and to develop practical computational models for use in system analysis codes for HE safety studies. Computer simulation with the CTH shock physics code is used to aid the experimental design and analyze the test results. In the calculations, initiation and growth or failure of detonation are modeled with the empirical HVRB model. The historical LANL SSGT and LSGT were reviewed and it was determined that a new, modified gap test be developed to satisfy the current requirements. In the new test, the donor/spacer/acceptor assembly is placed in a holder that is designed to work with fixtures for pre-damaging the acceptor sample. CTH simulations were made of the gap test with PBX-9501 samples pressed to three different densities. The calculated sensitivities were validated by test observations. The agreement between the computed and experimental critical gap thicknesses, ranging from 9 to 21 mm under various test conditions, is well within 1 mm. These results show that the numerical modeling is a valuable complement to the experimental efforts in studying and understanding shock initiation of PBX-9501.




THREE-DIMENSIONAL IGNITION AND GROWTH REACTIVE FLOW MODELING OF PRISM FAILURE TESTS ON PBX 9502


Book Description

The Ignition and Growth reactive flow model for shock initiation and detonation of solid explosives based on triaminotirnitrobenzene (TATB) is applied to three-dimensional detonation wave propagation. The most comprehensive set of three-dimensional detonation wave propagation data is that measured using the trapezoidal prism test. In this test, a PBX 9501 (95% HMX, 2.5% Estane, and 2.5% BDNPA/F) line detonator initiates a detonation wave along the trapezoidal face of a PBX 9502 (95% TATB and 5% Kel-F binder) prism. The failure thickness, which has been shown experimentally to be roughly half of the failure diameter of a long cylindrical charge, is measured after 50 mm of detonation wave propagation by impact with an aluminum witness plate. The effects of confinement impedance on the PBX 9502 failure thickness have been measured using air (unconfined), water, PMMA, magnesium, aluminum, lead, and copper placed in contact with the rectangular faces of the prism parallel to the direction of detonation propagation. These prism test results are modeled using the two-dimensional PBX 9502 Ignition and Growth model parameters determined by calculating failure diameter and tested on recent corner turning experiments. Good agreement between experimentally measured and calculated prism failure thicknesses for unconfined and confined PBX 9502 is reported.




Reactive Burn Models and Ignition & Growth Concept


Book Description

Plastic-bonded explosives are heterogeneous materials. Experimentally, shock initiation is sensitive to small amounts of porosity, due to the formation of hot spots (small localized regions of high temperature). This leads to the Ignition and Growth concept, introduced by Lee and Tarver in 1980, as the basis for reactive burn models. A homogeneized burn rate needs to account for three mesoscale physical effects (i) the density of burnt hot spots, which depends on the lead shock strength; (ii) the growth of the burn fronts triggered by hot spots, which depends on the local deflagration speed; (iii) a geometric factor that accounts for the overlap of deflagration wavelets from adjacent hot spots. These effects can be combined and the burn model defined by specifying the reaction progress variable [lambda](t) as a function of a dimensionless reaction length [tau]{sub hs}(t)/l{sub hs}, rather than by xpecifying an explicit burn rate. The length scale l{sub hs} is the average distance between hot spots, which is proportional to [N{sub hs}(P{sub s})]−13, where N{sub hs} is the number density of hot spots activated by the lead shock. The reaction length [tau]{sub hs}(t) = {line_integral}0{sup t} D(P(t'))dt' is the distance the burn front propagates from a single hot spot, where D is the deflagration speed and t is the time since the shock arrival. A key implementation issue is how to determine the lead shock strength in conjunction with a shock capturing scheme. They have developed a robust algorithm for this purpose based on the Hugoniot jump condition for the energy. The algorithm utilizes the time dependence of density, pressure and energy within each cell. The method is independent of the numerical dissipation used for shock capturing. It is local and can be used in one or more space dimensions. The burn model has a small number of parameters which can be calibrated to fit velocity gauge data from shock initiation experiments.




SHOCK INITIATION EXPERIMENTS ON PBX 9501 EXPLOSIVE AT PRESSURES BELOW 3 GPa WITH ASSOCIATED IGNITION AND GROWTH MODELING.


Book Description

Shock initiation experiments on the explosive PBX 9501 (95% HMX, 2.5% estane, and 2.5% nitroplasticizer by weight) were performed at pressures below 3 GPa to obtain in-situ pressure gauge data, run-distance-to-detonation thresholds, and Ignition and Growth modeling parameters. Propellant driven gas guns (101 mm and 155 mm) were utilized to initiate the PBX 9501 explosive with manganin piezoresistive pressure gauge packages placed between sample slices. The run-distance-to-detonation points on the Pop-plot for these experiments showed agreement with previously published data and Ignition and Growth modeling parameters were obtained with a good fit to the experimental data. This parameter set will allow accurate code predictions to be calculated for safety scenarios in the low-pressure regime (below 3 GPa) involving PBX 9501 explosive.