Purdue University Graduate School
Sakano_Thesis_Final.pdf (18.06 MB)


Download (18.06 MB)
posted on 2021-07-23, 01:51 authored by Michael N SakanoMichael N Sakano
For decades, the response of high-energy (HE) density materials at extreme conditions of pressure and temperature from strong insults like burning or impact have been studied in depth by the shock community. Shock physicists aim to develop a fundamental understanding for coupled chemical and physical processes across orders of magnitude spatial and temporal regimes. In order to succeed, this requires extensive collaboration between experiments and simulations, ranging from the electronic to the engineering scales. The end goals would be to develop predictive multiscale models capable of explaining ignition and initiation of HE systems and composites. The collected works in this thesis detail my contributions to the field of HE materials, specifically addressing the chemical reactivity at the atomistic level using reactive molecular dynamics (MD) simulations.
Through this endeavor, we aim to develop a critical understanding for the decomposition processes of HE materials. We begin with a validation the reactive force field, ReaxFF, by addressing the very strong anisotropic shock sensitivity in 2,2-Bis[(nitrooxy)methyl]propane-1,3-diyl dinitrate (PETN) through direct comparison of time-evolved spectra between experiments and simulations. Such strong orientation dependence is thought to relate to the initial decomposition events. Therefore we compare spectra at three different shock pressures, where we observe similar timescales for the disappearance of the NO2 symmetric and antisymmetric stretch modes. A more detailed chemical species analysis indicates that the NO2 molecular species could be considered the primary intermediate which initiates the decomposition process. Furthermore, these results suggest that the combination of explicit MD simulations and ultrafast spectroscopy will be key to the development of a detailed understanding of chemistry at extreme conditions.
Following the validation study, we further our understanding of reactivity in HE systems by investigating the differences in kinetics between an ordered and disordered system. It has been shown that shocked material is often severely strained, causing a loss in crystalline order. This in turn results in the disordered materials, such as amorphous solids, having
faster reactivity due to their higher internal energy and/or lower thermal conductivity. Our results indicate that extra energy is required to break the long-range order in bulk crystalline systems, thus resulting in slower decomposition rates. Further analyses of thermal hotspots point towards slightly faster chemical propagation in the amorphous samples due to lower thermal conductivity. These results provide an understanding for how molecular disorder can be attributed to increased reactivity.
After developing an understanding for the initial decomposition processes of HE materials, we turn our attention to a growing interest in the community which is the developing reduced order chemistry models for use in multiscale efforts. Many schemes report mechanisms that are obtained from experiments, which can have large error bars depending on the apparatus and/or extraction technique, or from gas phase simulations, which may not be relevant at shock conditions. To circumvent these issues, we develop a coarse-grained chemical kinetics model from all-atom reactive MD simulations by taking advantage of an unsupervised dimensionality reduction machine learning technique called non-negative matrix factorization. Doing so allows us to represent the overall decomposition chemistry as latent concentrations akin to reactants, intermediates, and products, which we then use to extract kinetics parameters and heats of reaction. These values are implemented into a continuum model, where we could simulate the criticality of thermal hotspots at regimes beyond the reach of MD, as well as verify how uncertainties in the parameters vary as a function of hotspot sizes.
Finally, we close with significant progress made towards on-going and future work, where we address two of the most challenging ideas in the field of HE materials: 1) developing definitive chemistry models at extreme conditions, and 2) improving coarse-grained descriptions for multiscale modeling.


US Office of Naval Research grant number N00014-16-1-2557


Degree Type

  • Doctor of Philosophy


  • Materials Engineering

Campus location

  • West Lafayette

Advisor/Supervisor/Committee Chair

Alejandro Strachan

Additional Committee Member 2

Marisol Koslowski

Additional Committee Member 3

Steven F. Son

Additional Committee Member 4

David F. Bahr