An official website of the United States government
Abstract Predictive simulations of the shock‐to‐detonation transition (SDT) in heterogeneous energetic materials (EM) are vital to the design and control of their energy release and sensitivity. Due to the complexity of the thermo‐mechanics of EM during the SDT, both macro‐scale response and sub‐grid mesoscale energy localization must be captured accurately. This work proposes an efficient and accurate multiscale framework for SDT simulations of EM. We introduce a new approach for SDT simulation by using deep learning to model the mesoscale energy localization of shock‐initiated EM microstructures. The proposed multiscale modeling framework is divided into two stages. First, a physics‐aware recurrent convolutional neural network (PARC) is used to model the mesoscale energy localization of shock‐initiated heterogeneous EM microstructures. PARC is trained using direct numerical simulations (DNS) of hotspot ignition and growth within microstructures of pressed HMX material subjected to different input shock strengths. After training, PARC is employed to supply hotspot ignition and growth rates for macroscale SDT simulations. We show that PARC can play the role of a surrogate model in a multiscale simulation framework, while drastically reducing the computation cost and providing improved representations of the sub‐grid physics. The proposed multiscale modeling approach will provide a new tool for material scientists in designing high‐performance and safer energetic materials.
more »
« less
Multi-scale modeling of shock initiation of a pressed energetic material III: Effect of Arrhenius chemical kinetic rates on macro-scale shock sensitivity
Multi-scale predictive models for the shock sensitivity of energetic materials connect energy localization (“hotspots”) in the microstructure to macro-scale detonation phenomena. Calculations of hotspot ignition and growth rely on models for chemical reaction rates expressed in Arrhenius forms; these chemical kinetic models, therefore, are foundational to the construction of physics-based, simulation-derived meso-informed closure (reactive burn) models. However, even for commonly used energetic materials (e.g., HMX in this paper) there are a wide variety of reaction rate models available. These available reaction rate models produce reaction time scales that vary by several orders of magnitude. From a multi-scale modeling standpoint, it is important to determine which model best represents the reactive response of the material. In this paper, we examine three global Arrhenius-form rate models that span the range of reaction time scales, namely, the Tarver 3-equation, the Henson 1-equation, and the Menikoff 1-equation models. They are employed in a meso-informed ignition and growth model which allows for connecting meso-scale hotspot dynamics to macro-scale shock-to-detonation transition. The ability of the three reaction models to reproduce experimentally observed sensitivity is assessed by comparing the predicted criticality envelope (Walker–Wasley curve) with experimental data for pressed HMX Class V microstructures. The results provide a guideline for model developers on the plausible range of time-to-ignition that are produced by physically correct Arrhenius rate models for HMX.
more »
« less
- Award ID(s):
- 2203580
- PAR ID:
- 10545274
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Journal of Applied Physics
- Volume:
- 135
- Issue:
- 8
- ISSN:
- 0021-8979
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract It has been suggested before that small-scale magnetic flux rope (SMFR) structures in the solar wind can temporarily trap energetic charged particles. We present the derivation of a new fractional Parker equation for energetic-particle interaction with SMFRs from our pitch-angle-dependent fractional diffusion-advection equation that can account for such trapping effects. The latter was derived previously in le Roux & Zank from the first principles starting with the standard focused transport equation. The new equation features anomalous advection and diffusion terms. It suggests that energetic-particle parallel transport occurs with a decaying efficiency of advection effects as parallel superdiffusion becomes more dominant at late times. Parallel superdiffusion can be linked back to underlying anomalous pitch-angle transport, which might be subdiffusive during interaction with quasi-helical coherent SMFRs. We apply the new equation to time-dependent superdiffusive shock acceleration at a parallel shock. The results show that the superdiffusive-shock-acceleration timescale is fractional, the net fractional differential particle flux is conserved across the shock ignoring particle injection at the shock, and the accelerated particle spectrum at the shock converges to the familiar power-law spectrum predicted by standard steady-state diffusive-shock-acceleration theory at late times. Upstream, as parallel superdiffusion progressively dominates the advection of energetic particles, their spatial distributions decay on spatial scales that grow with time. Furthermore, superdiffusive parallel shock acceleration is found to be less efficient if parallel anomalous diffusion is more superdiffusive, while perpendicular particle escape from the shock, thought to be subdiffusive during SMFR interaction, is reduced when increasingly subdiffusive.more » « less
-
We study numerical approximations of the reactive Euler equations of gas dynamics. In addition to shock, contact and rarefaction waves, these equations admit detonation waves appearing at the interface between different fractions of the reacting species. It is well-known that in order to resolve the reaction zone numerically, one has to take both space and time stepsizes to be proportional to the reaction time, which may cause the numerical method to become very computationally expensive or even impractical when the reaction is fast. Therefore, it is necessary to develop underresolved numerical methods, which are capable of accurately predicting locations of the detonation waves without resolving their detailed structure. One can distinguish between two different degrees of stiffness. In the stiff case, the reaction time is very small, while in the extremely stiff case, the reaction is assumed to occur instantaneously. In [A. Kurganov, in Hyperbolic problems: theory, numerics, applications, Springer, Berlin, 2003], we proposed a simple underresolved method—an accurate deterministic projection (ADP) method—for one-dimensional hyperbolic systems with stiff source terms including the reactive Euler equations in the extremely stiff regime. In this paper, we extend the ADP method to the (non-extremely) stiff case, multispecies detonation models, and the two-dimensional reactive Euler equations in all of the aforementioned regimes. We also investigate ways to distinguish between different regimes in practice as well as study the limitations of the proposed ADP methods with respect to the ignition temperature. We demonstrate the accuracy and robustness of the ADP methods in a number of numerical experiments with both relatively low and large ignition temperature, and illustrate the difficulties one may face when the ignition temperature is low.more » « less
-
null (Ed.)Evaporation and breakup of liquid droplets are common in many applications of the shock-driven multiphase instability (SDMI), such as in liquid-fuelled detonation engines, multiphase ejector pumps and turbines and explosive dispersal of liquid particles (i.e. chemical or biological agents). In this paper, the effects of evaporation and breakup of droplets on the mixing induced by the SDMI are considered through simulations and compared with experimental results. The evaporation model is validated against previous experimental data. The capabilities of the simulations and particle models are then demonstrated through a qualitative comparison with experimental results where breakup effects are negligible (i.e. small droplets). The simulation results are explored further to quantify the effects of evaporation (i.e. mixing enhancement) in the SDMI, providing further insight into the experimental results. A new breakup model, derived from previous works, is then presented for low Reynolds number (below 500), low Weber number (below 100) droplets in a shock-driven multiphase instability. The breakup model capabilities are then demonstrated through a comparison with experimental results where breakup effects are significant (larger droplet sizes). Finally, the simulation results are used to highlight the importance of breakup parameters on the evaporation rate and large-scale mixing in the SDMI. Overall, it is shown that evaporation is enhanced by the large-scale hydrodynamics instability, the SDMI, and that breakup of the droplets significantly increases the strength of the instability, and rate of droplet evaporation.more » « less
-
Nitrogen-rich heterocycles are essential for designing novel energetic green materials with the combination of high explosive performance and acceptable mechanical sensitivities. In this work, two sets of high nitrogen-azoles, derived from tetrazoles and triazole assemblies with N -trinitromethane, 5,5′-(2-(trinitromethyl)-2 H -1,2,3-triazole-4,5-diyl)bis(1 H -tetrazole) (TBTN) and N -methylene tetrazole, 5,5′-(2-((1 H -tetrazol-5-yl)methyl)-2 H -1,2,3-triazole-4,5-diyl)bis(1 H -tetrazole) (TBTT) are described. Their molecular structures were confirmed using multinuclear ( 1 H, 13 C, and 15 N) NMR spectra and single-crystal X-ray diffraction analysis. These molecules are attention attracting results emanating from methodologies utilized to access a unique class of tri-ionic salts in reaction with nitrogen-rich bases. The thermostabilities, mechanical sensitivities, and detonation properties of all new compounds were determined. Surprisingly, the nitro-based tri-cationic salts, 5b (Dv = 9376 m s −1 ) and 5c (Dv = 9418 m s −1 ), have excellent detonation velocities relative to HMX (Dv = 9144 m s −1 ), while those of the nitro-free tri-cationic salts, 8b·H2O (Dv = 8998 m s −1 ) and 8c·0.5H2O (Dv = 9058 m s −1 ), are superior to that of RDX (Dv = 8795 m s −1 ) and approach HMX values. Additionally, nearly all new compounds are insensitive to mechanical stimuli because of the high percentage of hydrogen bond interactions (HBs) between the anions and cations, which are evaluated using two-dimensional (2D) fingerprint and Hirshfeld surface analyses. It is believed that the work presented here is the first example of high-performing and insensitive tri-cationic energetic salts, which may establish a discovery platform for the “green” synthesis of future energetic materials.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0187735
