Jeffery Allen Leiding

Reliable prediction of detonation behavior in condensed-phase explosives requires models that offer quantitative robustness while maintaining computational efficiency for practically sized (centimeter-scale) systems. Advanced Detonation Shock Dynamics (ADSD) builds on the classical Detonation Shock Dynamics (DSD) framework by incorporating shock acceleration. In this work, we introduce methodological advances to ADSD to improve its ability to describe shock-to-detonation transition (SDT), diameter effect, and detonation failure with reactive-burn level accuracy at a substantially reduced computational cost. We then perform a comprehensive uncertainty quantification (UQ) analysis of the ADSD model parameters for PBX 9502 using experimental data from three experimental benchmarks: SDT, the diameter effect, and the shape effect. Statistical characterization is performed using UQ, yielding posterior parameter distributions that capture uncertainty and correlations; notably, many parameters showed low standard deviations, suggesting high confidence in their inferred values. The UQ-informed parameter sets are subsequently used to simulate corner-turning experiments, yielding predicted dead zones that show strong agreement with experimental observations. This combined optimization and Bayesian UQ framework provides the first systematic quantification of uncertainty in ADSD, enhancing confidence in its predictive capability and laying the groundwork for broader application of reduced-order detonation models in safety-critical engineering contexts.

Additive manufacturing, and particularly the vat photopolymerization process, enables the fabrication of complex geometries at high resolution and small length scales, making it well-suited for fabricating cellular structures (e.g., foams and lattices). Among these, elastomeric cellular structures are of growing interest due to their tunable compliance and energy dissipation. However, comprehensive data on the compressive behavior of these structures remains limited, especially for investigating the structure-property effects from changing the density and distribution of material within the cellular structure. This study explores how the mechanical response of polyurethane-based simple cubic structures changes when varying volume fraction, unit cell length, and unit cell patterning, which have not been systematically investigated previously in additively manufactured elastomers. Increasing volume fraction from 10% to 50% yielded significant changes in compressive stress–strain performance (decreasing strain at 0.5 MPa by 41.6% and increasing energy absorption density by 3962.5%). Although changing the unit cell length between 2.5 and 7 mm in ~30 mm parts did not result in statistically different stress–strain responses, modifying the configuration of struts of different thicknesses across designs with 30% volume fraction altered the stress–strain behavior (differences of 12.5% in strain at 0.5 MPa and 109.4% for energy absorption density). Power law relationships were developed to understand the interactions between volume fraction, unit cell length, and elastic modulus, and experimental data showed strong fits (R² > 0.91). These findings enhance the understanding of how multiple structural design aspects influence the performance of elastomeric cellular materials, providing a foundation for informing strategic design of tailorable materials for diverse mechanical applications.

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow model for high explosive PBX 9404 is developed. We specifically calibrate an AWSD model for PBX 9404 by fitting equations of state for reactants and detonation products to the results of thermochemical calculations and to experimental data from multiple sources. The calibrated equations of state are then coupled with an Arrhenius rate law based on shock temperature that describes the reaction progress during PBX 9404 detonation. The parameter values in the rate law are calibrated to experimental gas-gun data and diameter effect results. The results of the calibrated AWSD model are in strong agreement with available experimental data for PBX 9404. A similar level of agreement between predicted and experimental results is observed when the calibrated model is validated on data that were not used in the model parameterization procedure. Our results illustrate that the AWSD model is capable of accurately describing the many important properties and observables in the reactive burn of PBX 9404. Because of the historical significance of PBX 9404 in high explosives research and its current use in aging studies, this work provides an important model of a legacy material, which can be used to make comparisons to new high explosive formulations.

Due to its ability to achieve geometric complexity at high resolution and low length scales, additive manufacturing (AM) has increasingly been used for fabricating cellular structures (e.g., foams and lattices) for a variety of applications. Specifically, elastomeric cellular structures offer tunability of compliance as well as energy absorption and dissipation characteristics. However, there are limited data available on compression properties for printed elastomeric cellular structures of different designs and testing parameters. In this work, the authors evaluate how unit cell topology, part size, the rate of compression, and aging affect the compressive response of polyurethane-based simple cubic, body-centered, and gyroid structures formed by vat photopolymerization AM. Finite element simulations incorporating hyperelastic and viscoelastic models were used to describe the data, and the simulated results compared well with the experimental data. Of the designs tested, only the parts with the body-centered unit cell exhibited differences in stress–strain responses at different part sizes. Of the compression rates tested, the highest displacement rate (1000 mm/min) often caused stiffer compressive behavior, indicating deviation from the quasi-static assumption and approaching the intermediate rate response. The cellular structures did not change in compression properties across five weeks of aging time, which is desirable for cushioning applications. This work advances knowledge on the structure–property relationships of printed elastomeric cellular materials, which will enable more predictable compressive properties that can be traced to specific unit cell designs.

FeH is one of the most challenging diatomic molecules to study under electronic structure theory. Here, we have successfully studied 22 electronic states of FeH using ab initio multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] levels of theory. We report their potential energy curves (PECs), excitation energies, dissociation energies, equilibrium electronic configurations, and a series of spectroscopic constants with the use of augmented triple-ζ, quadruple-ζ, and quintuple-ζ quality correlation consistent basis sets. The scalar relativistic effects and active space and core electron correlation contribution on the properties of FeH are also explored. The use of a large CASSCF active space that includes 4s, 4p, 3d, and 4d orbitals of Fe and the 1s of H is critical for producing accurate full PECs with proper dissociations and predicting the exact order of the electronic states. Our findings are in harmony with the experimental results available in the literature and will serve as reference values for future studies of FeH. Furthermore, with the use of PECs, the total internal partition function sum (TIPS) of FeH was calculated across a range of temperatures. Finally, we exploited the single-reference nature of the a⁶Δ of FeH and its ionized product FeH⁺ (X⁵Δ) to evaluate the associated density functional theory (DFT) errors on their dissociation energies and spectroscopic parameters.

The U.S. Department of Energy's Office of Scientific and Technical Information

The U.S. Department of Energy's Office of Scientific and Technical Information

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.