Engineering topics

Hamilton, Brenden William

Publications and source records attributed to Hamilton, Brenden William.

9 records

DOE OSTI2023

High Pressure Suppression of Plasticity in RDX [Slides]

RDX shear banding nucleates without other defects present and is important for both mechanical properties and for mechano-chemical effects. Above ~1.0 km/s particle velocity, the shear band mechanism ceases to occur. It is not replaced by a phase transformation or other plastic event. Shear band initial sites over-nucleate and lower the VM tress drastically. This results in a fully reversible, high pressure mechanism that suppresses the formation of shear bands.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Interplay of Mechanochemistry and Material Processes in the Graphite to Diamond Phase Transformation

The manifestation of intramolecular strains in covalent systems is widely known to accelerate chemical reactions and open alternative reaction paths. This process is moderately well understood for isolated molecules and unimolecular processes. However, in condensed matter processes such as phase transformations, material properties and structure may influence typical mechanochemical effects. Therefore, we utilize steered molecular dynamics to induce out of plane strains in graphite and compress the system under a constant strain rate to induce phase transformation. We show that the out of plane strain allows phase transformations to initiate at small amounts of compressive strain. Yet, in contrast to typical mechanochemical results, the sum of compressive and out of plane work needed to form a diamond has a local minimum due to altered defect formation processes during phase transformation. Additionally, these altered processes slow the kinetics of the phase transformation, taking longer from initiation to total material transformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Using limited neural networks to assess relative mechanistic influence on shock heating in granular solids

The rapid compaction of granular media results in localized heating that can induce chemical reactions, phase transformations, and melting. However, there are numerous mechanisms in play that can be dependent on a variety of microstructural features. Machine learning techniques such as neural networks offer a ubiquitous method to develop models for physical processes. Limiting what kind of microstructural information is used as an input and assessing normalized changes in network error, the relative importance of different mechanisms can be inferred. Here we utilize binned, initial density information as network inputs to predict local shock heating in a granular high explosive trained from large-scale molecular dynamics simulations. In this study, the spatial extent of the density field used in the network is altered to assess the importance and relevant length scales of the physical mechanisms in play, where different microstructural features result in different predictive capabilities.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Influence of Pore Surface Structure and Contents on Shock-Induced Collapse and Energy Localization

The majority of computational efforts in the shockwave driven collapse of porosity focus on idealized pore geometries and structures. These give significant insight into the physical mechanisms in play, but omit some potentially important effects of real energetic materials. Here, to address this, we simulate pore collapse in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for four types of pore surfaces and fillings: an idealized cylindrical pore, a pore surface coated with polymer, a melted/amorphous pore surface, and a pore filled with gas molecules. We find that these perturbations to the pore structure result in only minor changes to the highly complex collapse mechanisms. However, significant differences in hotspot temperature are observed, with peak temperatures ranging over several hundred Kelvin. The polymer-coated system has the most effect, lowering the hotspot temperature by 300 K, and also limiting the manifestation of intramolecular strain energy that is known to induce mechanochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗

DOE OSTI2023

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. Further, this new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗