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
-
Poly- and Perfluorinated alkyl substances (PFAS) pose environmental and public health concerns. While incineration remains the most common PFAS remediation method, the complete combustion and pyrolysis mechanism of PFAS is unknown. This study aims to expand our understanding of the kinetics of gas-phase PFAS incineration by measuring the effect of difluoromethane (CHF) on propane ignition delay times (IDTs). The ignition delay times were measured by OH* emission and end-wall pressure time histories behind the reflected shock wave. Different concentrations of CH2F2 were mixed with fuel-lean propane-oxygen mixtures diluted in argon. Experiments were conducted at a nominal reflected shock pressure of P5 = 1 atm and reflected shock temperatures of 1200 < T5 < 1800 K. A new detailed chemical kinetic mechanism is presented. 135 new rate constants were computed using RRKM/ME theory, based upon stationary points computed using ANL0. The new mechanism is in excellent agreement with the measured ignition delay time. A novel sensitivity analysis helps to explain the elementary steps by which CH2F2 increases the ignition delay time.more » « less
-
Viscous shocks are a particular type of extreme event in nonlinear multiscale systems, and their representation requires small scales. Model reduction can thus play an essential role in reducing the computational cost for the prediction of shocks. Yet, reduced models typically aim to approximate large-scale dominating dynamics, which do not resolve the small scales by design. To resolve this representation barrier, we introduce a new qualitative characterization of the space–time locations of shocks, named the “shock trace,” via a space–time indicator function based on an empirical resolution-adaptive threshold. Unlike exact shocks, the shock traces can be captured within the representation capacity of the large scales, thus facilitating the forecast of the timing and locations of the shocks utilizing reduced models. Within the context of a viscous stochastic Burgers equation, we show that a data-driven reduced model, in the form of nonlinear autoregression (NAR) time series models, can accurately predict the random shock traces, with relatively low rates of false predictions. Furthermore, the NAR model, which includes nonlinear closure terms to approximate the feedback from the small scales, significantly outperforms the corresponding Galerkin truncated model in the scenario of either noiseless or noisy observations. The results illustrate the importance of the data-driven closure terms in the NAR model, which account for the effects of the unresolved dynamics brought by nonlinear interactions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0187735
