Search

Meso-scale calculations of energy localization and initiation in energetic material microstructures must capture the deformation and collapse of pores and high-temperature shear bands, which lead to hotspots. Because chemical reaction rates depend sensitively on temperature, predictive continuum models need to get the pore-collapse dynamics and resulting hotspot temperatures right; this imposes stringent demands on the fidelity of thermophysical model forms and parameters and on the numerical methods employed to perform high-resolution meso-scale calculations. Here, continuum material models for β-HMX are examined in the context of shock-induced pore collapse, treating predictions from all-atom molecular dynamics (MD) simulations as ground truth. Using atomistics-consistent material properties, we show that the currently available strength models for HMX fail to correctly capture pore collapse and hotspot temperatures. Insights from MD are then employed to advance a Modified Johnson–Cook (M-JC) strength model, which is shown to capture key aspects of the physics of shock-induced localization in HMX. The study culminates in a MD-guided strength model for β-HMX that produces continuum pore-collapse results in better alignment on several aspects with those predicted by MD, including pore-collapse mechanism and rate, shear-band formation in the collapse zone, and temperature, strain, and stress fields in the hotspot zone and the surrounding material. The resulting MD-informed/MD-determined M-JC model should improve the fidelity of meso-scale simulations to predict the detonation initiation of HMX-based energetic materials in microstructure-aware multi-scale frameworks. [ABSTRACT FROM AUTHOR]

High-fidelity constitutive models are critical for accurate meso-scale continuum modeling and prediction of shock initiation of crystalline energetic materials (EMs). While empirically calibrated or atomistic-guided anisotropic elastoplastic models of EM such as cyclotetramethylene-tetranitramine (HMX) and cyclotrimethylene-trinitramine (RDX) capture important micromechanical phenomena (such as dislocation evolution, slip-resistance, and anisotropic elasticity), the computational cost of using anisotropic single-crystal plasticity models can become prohibitive for meso-scale computations of void-collapse and hotspot formation in microstructures. Thermo-mechanically representative, isotropic, pressure, temperature, and rate-dependent material constitutive models are practical alternatives for meso-scale simulations of the shock response of microstructures. To this end, this work constructs physically consistent isotropic plasticity from anisotropic single-crystal plasticity models for HMX and RDX. State-of-the-art crystal plasticity models for HMX and RDX are used to compute the stress states in single crystals oriented in three different directions relative to shocks generated by impact at velocities ranging from 100 to 1000 m/s. Post-shock von Mises stress fields for the three orientations are then used to calibrate the strain-rate hardening coefficient and the reference strain rate for a rate-dependent Johnson–Cook (JC) yield surface model. We compare the pressures and the post-shock von Mises stresses between the JC and the anisotropic models to show that the isotropic computations closely approximate the averaged deformation response of the three different crystal orientations. We then model the interaction of a shock generated by a 500 m/s impact with a 0.5 μm void and show that the pressures and the deviatoric stresses obtained using the isotropic model closely match those computed from anisotropic models for both HMX and RDX. The resulting isotropic J2 plastic flow model for HMX and RDX can be employed to perform meso-scale simulations for energy localization due to shear bands and void collapse in the two materials. [ABSTRACT FROM AUTHOR]

Material models for single-crystal β-HMX are systematically examined in the context of continuum pore-collapse simulations. Continuum predictions using five different isotropic material models are compared head-to-head with molecular dynamics (MD) predictions for a 50 nm cylindrical pore in β-HMX subject to a range of shock strengths. Shock waves were generated using a reverse-ballistic configuration, propagating along [010] in the MD simulations. The continuum models are improved hierarchically, drawing on temperature- and pressure-dependent MD-derived material parameters. This procedure reveals the sensitivity of the continuum predictions of pore collapse to the underlying thermophysical models. The study culminates in an MD-calibrated isotropic rate- and temperature-dependent strength model, which includes appropriate submodels for the temperature-dependent melting point of β-HMX [M. P. Kroonblawd and R. A. Austin, Mech. Mater. 152, 103644 (2021)], pressure-dependent shear modulus [A. Pereverzev and T. Sewell, Crystals 10, 1123 (2020)], and temperature-dependent specific heat, that produces continuum pore-collapse results similar to those predicted by MD. The resulting MD-informed model should improve the fidelity of simulations to predict the detonation initiation of HMX-based energetic materials containing micrometer-scale pores. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) and correlated molecular orbital electronic structure calculations were used to study the Al + CO 2 → AlO + CO reaction on the electronic ground-state potential-energy surface (PES). Geometries were optimized using DFT (M11/jun-cc-pV(Q+d)Z) and more accurate energies were obtained using the composite Weizmann-1 theory with Brueckner doubles (W1BD). The results comprise the most complete, most systematic characterization of the Al + CO 2 reaction surface to date and are based on consistent application of high-level methods for all stationary points identified. The pathways from Al + CO 2 to AlO + CO on the electronic ground-state PES all involve formation of one or more stable AlCO 2 complexes denoted η-AlCO 2 , trans-AlCO 2 , and C 2v -AlCO 2 , among which η-AlCO 2 and C 2v -AlCO 2 are the least and most stable, respectively. We report a new minimum-energy pathway for the overall reaction, namely formation of η-AlCO 2 from reactants and dissociation of that same complex to products via a bond-insertion reaction that passes through a fourth (weakly metastable) AlCO 2 complex denoted cis-OAlCO. Natural Bond Orbital analysis was applied to study trends in charge distribution and the degree of charge transfer in key structures along the minimum-energy pathway. The process of aluminum insertion into CO 2 is discussed in the context of analogous processes for boron and first-row transition metals.