XHVRB reactive flow modeling of the PBX 9501 gap-stick experiment
Explore the source record for details and available documents.
Damm, David Lee [Sandia National Laboratories (SNL↗
SEARCH · Engineering Papers
Results for “PBX 9501”
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.
Detonation Front Interaction With Meso-Scale 2 Void Systems in PBX 9501
Explore the source record for details and available documents.
36 MATERIALS SCIENCE↗
Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids
Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.
42 ENGINEERING↗
Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state
Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.
36 MATERIALS SCIENCE↗
Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization
Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.
97 MATHEMATICS AND COMPUTING↗
Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models
Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.
97 MATHEMATICS AND COMPUTING↗
Water Sorption Study of Various Polymers and Polymeric Composites
In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.
36 MATERIALS SCIENCE↗
Numerical calculation of shock-induced initiation of detonations
Results of 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 parameters in the Forest Fire burn model are 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 are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.
Cort, G. E.↗
A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots
While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.
97 MATHEMATICS AND COMPUTING↗