Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Modelling the Fluctuations of Reactive Shock Waves in Heterogeneous Solid Explosives as Stochastic Processes

Journal Article · · Propellants, Explosives, Pyrotechnics
Although current phenomenological burn models are useful for describing the average or bulk reactive flow behaviour of heterogeneous explosives, one fundamental weakness inherent to these models is the loss of detailed microstructural information at the scale of the calculation. In order to include the effects of the microstructure, and in particular the underlying material heterogeneities that influence the build-up to detonation, a new paradigm is put forth for modelling sub-grid, reaction-induced fluctuations (i.e. “hot spots”) at the continuum level. This modelling approach assumes that the reaction rate is stochastic, rather than deterministic, and it uses Langevin-type equations with a mathematical framework built upon Itô calculus and Lambourn's CIM model for shocked heterogeneous explosives. This approach follows directly from our previous letter, Ref. [1], and is inspired by the probability density function (pdf) methods used in turbulent reactive flows. Here, the stochastic burn model is derived, implemented, and exercised far beyond what has been shown in previous work. New hydrocode simulation results demonstrate the role of stochastic fluctuations during shock initiation; these fluctuations are approximated by collections of discrete particles, that evolve with drift (i.e. deterministic) and diffusion (i.e. stochastic) coefficients. Furthermore, the particle values are propagated and averaged to calculate the heat release, yield strength, and material impedance in each computational cell. Hydrocode simulation results further show how the fluctuating hot spot energy may or may not be transmitted to the wave front, and result in a detonation wave-like structure. The fundamental stochastic nature of this model permits simulations to have varying outcomes with the same initial conditions; this allows for go/no-go estimation (e.g., marginal or failed detonations), which might possibly be calibrated using the statistical distributions from real materials. Mesoscale calculations of shocked heterogeneous explosives also show that these fluctuations are physically justified (i.e., Ref. [2]), and it is hypothesized that the pdf functions provide a link between the meso(grain) and continuum scales for practical engineering calculations. Hence, our approach is a paradigm shift for building efficient, reduced order continuum burn models, that represent the sub-grid stochastic behaviour of shocked heterogeneous solid explosives.
Research Organization:
Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
Grant/Contract Number:
AC04-94AL85000; NA0003525
OSTI ID:
1595423
Alternate ID(s):
OSTI ID: 1602904
Report Number(s):
SAND--2020-0010J; 681962
Journal Information:
Propellants, Explosives, Pyrotechnics, Journal Name: Propellants, Explosives, Pyrotechnics Journal Issue: 2 Vol. 45; ISSN 0721-3115
Publisher:
WileyCopyright Statement
Country of Publication:
United States
Language:
English

References (33)

A Consistent Hybrid Finite-Volume/Particle Method for the PDF Equations of Turbulent Reactive Flows journal September 1999
Stochastic Differential Equations book January 2003
On the effect of grain size on shock sensitivity of heterogeneous high explosives journal April 1997
PDF methods for turbulent reactive flows journal January 1985
Modeling heterogeneous energetic materials at the mesoscale journal February 2002
A measure of mixing from Lagrangian tracking and its application to granular and fluid flow systems journal December 2008
Relating microstructure, temperature, and chemistry to explosive ignition and shock sensitivity journal April 2018
Energy dissipation in polymer-bonded explosives with various levels of constituent plasticity and internal friction journal March 2019
A filtered mass density function approach for modeling separated two-phase flows for LES I: Mathematical formulation journal March 2006
A filtered mass density function approach for modeling separated two-phase flows for LES II: Simulation of a droplet laden temporally developing mixing layer journal August 2008
A review of predictive nonlinear theories for multiscale modeling of heterogeneous materials journal February 2017
On the use of reverse Brownian motion to accelerate hybrid simulations journal April 2017
Prediction of shock initiation thresholds and ignition probability of polymer-bonded explosives using mesoscale simulations journal May 2018
Coarse-Grain Model Simulations of Nonequilibrium Dynamics in Heterogeneous Materials journal June 2014
Shock Initiation of Solid Explosives journal January 1961
An Interpretation of Particle Velocity Histories during Growth to Detonation
  • Lambourn, B. D.
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2003: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings https://doi.org/10.1063/1.1780255
conference January 2004
Measurements of shock initiation in the tri-amino-tri-nitro-benzene based explosive PBX 9502: Wave forms from embedded gauges and comparison of four different material lots journal June 2006
On the systematics of particle velocity histories in the shock-to-detonation transition regime journal October 2006
The Hugoniot and shock sensitivity of a plastic‐bonded TATB explosive PBX 9502 journal May 1988
Mixture model for determination of shock equation of state journal April 2012
Mesoscale simulation of reactive pressed energetic materials under shock loading journal December 2015
Stochasticity in materials structure, properties, and processing—A review journal March 2018
Shock interactions with heterogeneous energetic materials journal March 2018
Letter: Modeling reactive shock waves in heterogeneous solids at the continuum level with stochastic differential equations journal May 2018
Quantification of probabilistic ignition thresholds of polymer-bonded explosives with microstructure defects journal October 2018
Multi-scale shock-to-detonation simulation of pressed energetic material: A meso-informed ignition and growth model journal August 2018
Phenomenological model of shock initiation in heterogeneous explosives journal January 1980
Highly scalable discrete-particle simulations with novel coarse-graining: accessing the microscale journal May 2018
Probabilistic models for reactive behaviour in heterogeneous condensed phase media journal February 2012
Multi-dimensional mesoscale simulations of detonation initiation in energetic materials with density-based kinetics journal December 2017
Initiation and Growth of Explosion in Liquids and Solids journal April 1952
Lagrangian PDF Methods for Turbulent Flows journal January 1994
A pseudo-random number generator for the System/360 journal January 1969

Similar Records

Letter: Modeling reactive shock waves in heterogeneous solids at the continuum level with stochastic differential equations
Journal Article · Tue May 08 20:00:00 EDT 2018 · Physics of Fluids · OSTI ID:1469625

A new approach to the reactive burn modeling of heterogeneous explosives
Conference · Wed Dec 31 23:00:00 EST 2003 · OSTI ID:977471