Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pop plot”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Numerical computation of Pop plot

The Pop plot — distance-of-run to detonation versus initial shock pressure — is a key characterization of shock initiation in a heterogeneous explosive. Reactive burn models for high explosives (HE) must reproduce the experimental Pop plot to have any chance of accurately predicting shock initiation phenomena. This report describes a methodology for automating the computation of a Pop plot for a specific explosive with a given HE model. Illustrative examples of the computation are shown for PBX 9502 with three burn models (SURF, WSD and Forest Fire) utilizing the xRage code, which is the Eulerian ASC hydrocode at LANL. Comparison of the numerical and experimental Pop plot can be the basis for a validation test or as an aid in calibrating the burn rate of an HE model. Issues with calibration are discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2018

Review of PBX 9502 Pop Plot data

Pop plot data (distance of run and time to detonation as a function of shock pressure) is a simple characterization of shock initiation of a high explosive. It is used, along with other data, to calibrate reactive burn models. For the insensitive TATB based explosive PBX 9502, there are data at 5 initial temperatures (77 K, -55 C, 23 C, 76 C, 130 C) and for several lots. Here we review the dependence of the Pop plot on temperature and the variation with lot. There is a significant increase in shock sensitivity with temperature.Since the Pop plot is important for modeling, fit parameters are interpolated to cover the temperature range –55 < T < 76 C. The scatter in the data and variation with lot is important for assessing the accuracy to which a model can be calibrated.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗
DOE OSTI2019

TATB Sensitivity to Shocks from Electrical Arcs

Use of insensitive high explosives (IHEs) has significantly improved ammunition safety because of their remarkable insensitivity to violent cook–off, shock and impact. Triamino–trinitrobenzene (TATB) is the IHE used in many modern munitions. Previously, lightning simulations in different test configurations have shown that the required detonation threshold for standard density TATB at ambient and elevated temperatures (250 °C) has a sufficient margin over the shock caused by an arc from the most severe lightning. In this paper, the Braginskii model with Lee–More channel conductivity prescription is used to demonstrate how electrical arcs from lightning could cause detonation in TATB. The steep rise and slow decay in typical lightning pulse are used in demonstrating that the shock pressure from an electrical arc, after reaching the peak, falls off faster than the inverse of the arc radius. For detonation to occur, two necessary detonation conditions must be met: the Pop–Plot criterion and minimum spot size requirement. The relevant Pop–Plot for TATB at 250 °C was converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. The arc cross–section was required to meet the minimum detonation spot size reported in the literature. Furthermore, one caveat is that when the shock pressure exceeds the detonation pressure the Pop–Plot may not be applicable, and the minimum spot size requirement may be smaller.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

SURF Model Calibration Strategy

SURF and SURFplus are high explosive reactive burn models for shock initiation and propagation of detonation waves. They are engineering models motivated by the ignition & growth concept of high spots and for SURFplus a second slow reaction for the energy release from carbon clustering. A key feature of the SURF model is that there is a partial decoupling between model parameters and detonation properties. This enables reduced sets of independent parameters to be calibrated sequentially for the initiation and propagation regimes. Here we focus on a methodology for tting the initiation parameters to Pop plot data based on 1-D simulations to compute a numerical Pop plot. In addition, the strategy for tting the remaining parameters for the propagation regime and failure diameter is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

Relating microstructure, temperature, and chemistry to explosive ignition and shock sensitivity

An analysis is put forth that relates explosive properties including microstructure, temperature, and chemistry to explosive ignition and sensitivity. Unlike approaches that focus only on elucidating mesoscopic mechanisms important to ignition, the present analysis seeks a methodology directly applicable to continuum explosive burn models. Because the Scaled Uniform Reactive Front (SURF) burn model has a microstructural basis, it is chosen as the starting point of the present analysis. To build upon the SURF framework, a literal translation is performed of the high-level conceptual notions for which SURF is based, to concrete ignition–combustion parameters, a statistical description of the microstructure, and other thermo-chemical data acquired by non-shock experiments. The analysis requires a void volume fraction distribution acquired from ultra-small angle neutron scattering (USANS). The shock response of PBX 9502 is used to illustrate the theory. While the analysis requires calibration to a shock experiment (the pop plot, for example), the results are a self-consistent set of realistic physical quantities that contribute to the shock initiation process including, initial temperature, chemical kinetics, a statistical description of the microstructure, and the hot spot size, spacing, and temperature. The utility of having a mesoscale based theory is that once the critical hot spot temperature, as a function of shock pressure, is known for a given explosive and initial porosity distribution, the entire set of meso- and macro-scale results, including the pop plot can be calculated for any other porosity. Thus, one can understand changes in sensitivity at different densities (void size distributions), relying only on the assumption that the hot spot temperature curve would not differ significantly if the morphology of the hots spots were similar. Another important utility of the present analysis is to address the question of the role of initial temperature on the observed shock sensitivity of PBX 9502 and other explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

PETN spark-gap detonators

The well-developed theory of Lorentz plasma that is dominated by electron–ion interactions is used to calculate the PETN arc characteristics. The spark-gap discharge current is a ramp with 10 to 25 ns rise time to peak and remaining constant subsequently. The approximate formulas for the arc channel conductivity, arc temperature, arc radius, and shock pressure from the arc are obtained from a system of nonlinear ordinary differential equations, which is the similarity solution of hydrodynamic equations similar to the Braginskii approximation. These arc parameters are given for the peak current ranging from 100 A to 1000 A and with different rise times. Representative cases are compared to the nonlinear ordinary differential equation code results. The shock pressures at the peak current are comparable to those from a typical commercial EBW bridgewire burst reported in the literature; the arc radius at the peak current is comparable to a typical bridgewire diameter of 0.0375 mm (e.g., RISI detonators, RP-1, and RP-80). The relevant Pop-Plot for low-density PETN is converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. Finally, this criterion is then used to determine the detonation thresholds, which are comparable with test data obtained by Tucker, et al.

42 ENGINEERING↗
DOE OSTI2017

Comparison Between Surf and Multi-Shock Forest Fire High Explosive Burn Models

PAGOSA1 has several different burn models used to model high explosive detonation. Two of these, Multi-Shock Forest Fire and Surf, are capable of modeling shock initiation. Accurately calculating shock initiation of a high explosive is important because it is a mechanism for detonation in many accident scenarios (i.e. fragment impact). Comparing the models to pop-plot data give confidence that the models are accurately calculating detonation or lack thereof. To compare the performance of these models, pop-plots2 were created from simulations where one two cm block of PBX 9502 collides with another block of PBX 9502.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Prediction of Probabilistic Detonation Threshold via Millimeter-Scale Microstructure-Explicit and Void-Explicit Simulations

We present an approach and relevant models for predicting the probabilistic shoc-to-detonation transition (SDT) behavior and Pop plot (PP) of heterogeneous energetic materials (HEM) via mesoscopic microstructure-explicit (ME) and void explicit (VE) simulations at the millimeter (mm) sample size scale. Although the framework here is general, the particular material considered in this paper is pressed Octahydro-1,3,5,7-tetranitro-1,2,3,5-tetrazocine (HMX). To systematically delineate the effects of material heterogeneities, four material cases are considered. These cases are homogeneous material, material with granular microstructure but no voids, homogeneous material with voids, and material with both granular microstructure and voids. Statistically equivalent microstructure sample sets (SEMSS) are generated and used. Eulerian hydrocode simulations explicitly resolve the material heterogeneities, voids, and the coupled mechanical-thermal-chemical processes. In particular, it is found that both microstructure and voids strongly influence the SDT behavior and PP. The effects of different combinations of microstructure heterogeneity and voids on the SDT process and PP are quantified and rank-ordered. The overall framework uses the Mie–Grüneisen equation of state and a history variable reactive burn model (HVRB). A novel probabilistic representation for quantifying the PP is developed, allowing the calculation of (1) the probability of observing SDT at a given combination of shock pressure and run distance, (2) the run-distance to detonation under a given combination of shock pressure and prescribed probability, and (3) the shock pressure required for achieving SDT at a given run distance with a prescribed probability. The results are in agreement with general trends in experimental data in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2018

One dimensional shock initiation of hot isostatically pressed composition B

Twelve one dimensional shock experiments have been carried out on the UK isostatically pressed Composition B (59.5% RDX, 39.5% TNT, 1% wax) comprising of seven sustained pulse experiments and three short shock experiments using the embedded gauge technique at the Los Alamos National Laboratory gas gun facility. The evolution of the reactive growth at and behind the shock front has been measured along with the run to detonation distance. Two low pressure experiments have also been carried out to obtain non-reactive equation of state data. These data have been used to create the Pop plot and Hugoniot states for the UK Comp B. The shock initiation behavior of the UK Comp B has been compared to that of the equivalent US composition. The reactive growth shows a feature that was observed in the US composition in which the wave profiles display a high level of pre-detonation noise. This was hypothesized to be due to a piezoelectric effect in the RDX crystals. The results of these experiments have shown that this effect may be localized in the ε phase at shock pressures in the region of 4.15 GPa and above.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Three-dimensional microstructure-explicit and void-explicit mesoscale simulations of detonation of HMX at millimeter sample size scale

Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.

Miller, Christopher (ORCID:0000000206844381)↗
DOE OSTI2020

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shock initiation of cyclotol (75/25) at both ambient temperature and 70 °C

Here, six shock initiation experiments have been carried out to study the shock sensitivity and develop the unreacted Hugoniot for the Cyclotol composition 75/25 %wt RDX/TNT. The experiments were carried out both at ambient temperature, and elevated to 70 °C, close to the TNT melt temperature of 80 °C. Two sets of three experiments were fired at nominal pressures of 4.1, 6.0 and 7.7 GPa. Comparison has been made to the previously studied (80/20 %wt RDX/TNT) composition and analysis shows that the five percent decrease in RDX content has resulted in a significant increase in shock sensitivity. The increase in initial temperature has yielded an unexpectedly small increase in shock sensitivity, particularly relative to the effects of the compositional change. Wave profiles will be presented, along with a discussion of the unreacted Hugoniot and relative sensitivities presented in the Pop plot.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A study of shock initiation experiments for the explosive PBX 9502 using three reactive burn models

Shock to detonation transition (SDT) experiments are essential in calibrating and validating reactive burn models for explosives. This work investigates the large collection of SDT test data for the explosive PBX 9502 at ambient temperature that was presented by Gustavsen, Sheffield, and Alcon [Journal of Applied Physics, 99, 114907 (2006)]. We first analyze the experimental data and compare two different methods of determining the shock transition time/distance (namely, the bilinear method and the single-curve method). This reveals some of the uncertainty in estimating shock transition points, which contributes to scatter in Pop-plot data. Next, we compare the WSD, AWSD, and SURFplus reactive burn models for a collection of approximately 20 experimental shots using the FLAG hydrocode. Error estimates are used to quantify how well each reactive burn model (and their respective parameter calibrations) performs at predicting the SDT process for a range of loading conditions. Additionally, the importance of mesh resolution and numerical dissipation in SDT simulations will be assessed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Shock initiation of low density polymer bonded explosive LX-14: A study of two morphologies

A series of six shock initiation experiments have been carried out on low density LX-14 powder in order to simulate a shock insult on the heavily damage polymer bonded explosive, LX-14. Two distinct morphologies were studied (tap density molding powder and machined swarf), both at the same density of 0.942 g/cm 3 , or 50.1% theoretical maximum density. The purpose of these experiments was to provide data to help make an assessment of the effects that damage has on the material sensitivity to a planar shock. This was achieved primarily by providing shock sensitivity data in the form of a Pop plot, and also reactants equation of state data, which aids in the determination of input conditions for both the experiments performed in this work and also future experiments of this material type. The experiments were of a cut-back format, consisting of four sample heights on each shot and diagnosed with optical velocimetry. The experiments were carried out at the Technical Area 40 Chamber 9 gas gun facility at the Los Alamos National Laboratory, where the LX-14 targets were subjected to Al 6061 and Oxygen Free High Conductivity copper impactors launched to velocities up to 2.14 km/s. Time corrected reactive growth wave profiles are presented in this paper along with the derivation of the following Hugoniot parameters for this explosive, where the molding powder and machine swarf are represented by the following linear equations, respectively, U s = 1.64 (±0.64)u p + 2.25 (±0.49) and U s = 0.87 (±.0.60)u p + 4.08 (±0.49). The results show that the steady increase in shock sensitivity with increasing void fraction reaches an inflection point beyond which the shock sensitivity begins to decease. This inflection point lies between 65% and 50% of the theoretical maximum density of the LX-14. In both cases, the damaged LX-14 was not as sensitive as expected, with a relative shock sensitivity of the molding powder being less than the pressed LX-14, and the machined swarf having a shock sensitivity that is comparable to pressed LX-14.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Failure Diameter of PBX 9502: Simulations with the SURFplus model

SURFplus is a reactive burn model for high explosives aimed at modelling shock initiation and propagation of detonation waves. It utilizes the SURF model for the fast hot-spot reaction plus a slow reaction for the energy released by carbon clustering. A feature of the SURF model is that there is a partially decoupling between burn rate parameters and detonation wave properties. Previously, parameters for PBX 9502 that control shock ini- tiation had been calibrated to Pop plot data (distance-of-run to detonation as a function of shock pressure initiating the detonation). Here burn rate parameters for the high pres- sure regime are adjusted to t the failure diameter and the limiting detonation speed just above the failure diameter. Simulated results are shown for an uncon ned rate stick when the 9502 diameter is slightly above and slightly below the failure diameter. Just above the failure diameter, in the rest frame of the detonation wave, the front is sonic at the PBX/air interface. As a consequence, the lead shock in the neighborhood of the interface is supported by the detonation pressure in the interior of the explosive rather than the reaction immediately behind the front. In the interior, the sonic point occurs near the end of the fast hot-spot reaction. Consequently, the slow carbon clustering reaction can not a ect the failure diameter. Below the failure diameter, the radial extent of the detonation front decreases starting from the PBX/air interface. That is, the failure starts at the PBX boundary and propagates inward to the axis of the rate stick.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

SURFplus Model Calibration for PBX 9502

The SURFplus reactive burn model is calibrated for the TATB based explosive PBX 9502 at three initial temperatures; hot (75 C), ambient (23 C) and cold (-55 C). The CJ state depends on the initial temperature due to the variation in the initial density and initial specific energy of the PBX reactants. For the reactants, a porosity model for full density TATB is used. This allows the initial PBX density to be set to its measured value even though the coeffcient of thermal expansion for the TATB and the PBX differ. The PBX products EOS is taken as independent of the initial PBX state. The initial temperature also affects the sensitivity to shock initiation. The model rate parameters are calibrated to Pop plot data, the failure diameter, the limiting detonation speed just above the failure diameters, and curvature effect data for small curvature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

Failure Diameter Resolution Study

Previously the SURFplus reactive burn model was calibrated for the TATB based explosive PBX 9502. The calibration was based on fitting Pop plot data, the failure diameter and the limiting detonation speed, and curvature effect data for small curvature. The model failure diameter is determined utilizing 2-D simulations of an unconfined rate stick to find the minimum diameter for which a detonation wave propagates. Here we examine the effect of mesh resolution on an unconfined rate stick with a diameter (10mm) slightly greater than the measured failure diameter (8 to 9 mm).

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗