Doucet, Mathieu; Browning, James F.; Doyle, Barney L.; Charlton, Timothy R.; Ambaye, Haile; Seo, Joohyun; Mazza, Alessandro R.; Wenzel, John F.; Burns, George B.; Wixom, Ryan R.; Veith, Gabriel M.
Haynes 230 nickel alloy is one of the main contenders for salt containment in the design of thermal energy storage systems based on molten salts. A key problem for these systems is understanding the corrosion phenomena at the alloy–salt interface, and, in particular, the role played by chromium in these processes. In this study, thin films of Haynes 230, which is also rich in chromium, were measured with polarized neutron reflectometry and Rutherford backscattering spectrometry as a function of annealing temperature. Migration of chromium to the surface was observed for films annealed at 400 and 600 °C. Combining the two techniques determined that more than 60% of chromium comprising the as-prepared Haynes 230 layer moves to the surface when annealed at 600 °C, where it forms an oxide layer.
The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet, the majority of computational studies aimed at predicting phenomena affected by these processes, such as the initiation and propagation of detonation waves in explosives or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed. The measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of the microstructure along with a fully dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on the dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide a clear insight into the nature of threshold behavior and are a way to understand complex physical phenomena.
Image processing and stereological techniques were used to characterize the heterogeneity of composite propellant and inform a predictive burn rate model. Composite propellant samples made up of ammonium perchlorate (AP), hydroxyl-terminated polybutadiene (HTPB), and aluminum (Al) were faced with an ion mill and imaged with a scanning electron microscope (SEM) and x-ray tomography (micro-CT). Properties of both the bulk and individual components of the composite propellant were determined from a variety of image processing tools. An algebraic model, based on the improved Beckstead-Derr-Price model developed by Cohen and Strand, was used to predict the steady-state burning of the aluminized composite propellant. In the presented model the presence of aluminum particles within the propellant was introduced. The thermal effects of aluminum particles are accounted for at the solid-gas propellant surface interface and aluminum combustion is considered in the gas phase using a single global reaction. Properties derived from image processing were used directly as model inputs, leading to a sample-specific predictive combustion model.
Hexanitrostilbene (HNS) films were deposited onto polycarbonate substrates using vacuum thermal sublimation. The deposition conditions were varied in order to alter porosity in the films, and the resulting microstructures were quantified by analyzing ion-polished cross-sections using scanning electron microscopy. The effects of these changes in microstructure on detonation velocity and the critical thickness needed to sustain detonation were determined. The polycarbonate substrates also acted as recording plates for detonation experiments, and films near the critical thickness displayed distinct patterns in the dent tracks that indicate instabilities in the detonation front when approaching failure conditions.
Physical vapor deposition is a technique that can be used to produce explosive films with controlled geometry and microstructure. Films of the high explosive hexanitroazobenzene (HNAB) were deposited by vacuum thermal evaporation. HNAB deposits in an amorphous state that crystallizes over time into a polycrystalline material with high density and a consistent porosity distribution. In previous work, we evaluated detonation critical thickness in HNAB films in an effectively infinite slab geometry with insignificant side losses. In this work, the effect of geometry on detonation failure was investigated by performing experiments on films with different thicknesses, while also changing lateral dimensions such that side losses became significant. The experimental failure thickness was determined to be 75.5 μm and 71.6 μm, for 400 μm and 1600 μm wide HNAB lines, respectively. It follows from this that the minimum width to achieve detonation behavior representing an infinite slab configuration is greater than 400 μm.
We report a series of time-resolved spectroscopic measurements that aim to characterize the reactions that occur during shock initiation of high explosives. The experiments employ time-and wavelength-resolved emission spectroscopy to analyze light emitted from detonating thin explosive films. This paper presents analysis of optical emission spectra from hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN) thin film samples. Both vibrationally resolved and broadband emission features are observed in the spectra and area as electronic transitions of intermediate species.
In both continuum hydrodynamics simulations and also multimillion atom reactive molecular dynamics simulations of shockwave propagation in single crystal pentaerythritol tetranitrate (PETN) containing a cylindrical void, we observed the formation of an initial radially symmetric hot spot. By extending the simulation time to the nanosecond scale, however, we observed the transformation of the small symmetric hot spot into a longitudinally asymmetric hot region extending over a much larger volume. Performing reactive molecular dynamics shock simulations using the reactive force field (ReaxFF) as implemented in the LAMMPS molecular dynamics package, we showed that the longitudinally asymmetric hot region was formed by coalescence of the primary radially symmetric hot spot with a secondary triangular hot zone. We showed that the triangular hot zone coincided with a double-shocked region where the primary planar shockwave was overtaken by a secondary cylindrical shockwave. The secondary cylindrical shockwave originated in void collapse after the primary planar shockwave had passed over the void. A similar phenomenon was observed in continuum hydrodynamics shock simulations using the CTH hydrodynamics package. The formation and growth of extended asymmetric hot regions on nanosecond timescales has important implications for shock initiation thresholds in energetic materials.
This report summarizes a project in which the authors sought to develop and deploy: (i) experimental techniques to elucidate the complex, multiscale nature of thermal transport in particle-based materials; and (ii) modeling approaches to address current challenges in predicting performance variability of materials (e.g., identifying and characterizing physical- chemical processes and their couplings across multiple length and time scales, modeling information transfer between scales, and statically and dynamically resolving material structure and its evolution during manufacturing and device performance). Experimentally, several capabilities were successfully advanced. As discussed in Chapter 2 a flash diffusivity capability for measuring homogeneous thermal conductivity of pyrotechnic powders (and beyond) was advanced; leading to enhanced characterization of pyrotechnic materials and properties impacting component development. Chapter 4 describes success for the first time, although preliminary, in resolving thermal fields at speeds and spatial scales relevant to energetic components. Chapter 7 summarizes the first ever (as far as the authors know) application of TDTR to actual pyrotechnic materials. This is the first attempt to actually characterize these materials at the interfacial scale. On the modeling side, new capabilities in image processing of experimental microstructures and direct numerical simulation on complicated structures were advanced (see Chapters 3 and 5). In addition, modeling work described in Chapter 8 led to improved prediction of interface thermal conductance from first principles calculations. Toward the second point, for a model system of packed particles, significant headway was made in implementing numerical algorithms and collecting data to justify the approach in terms of highlighting the phenomena at play and pointing the way forward in developing and informing the kind of modeling approach originally envisioned (see Chapter 6). In both cases much more remains to be accomplished.
Density Functional Theory (DFT) has emerged as an indispensable tool in materials research, since it can accurately predict properties of a wide variety of materials at both equilibrium and extreme conditions. However, for organic molecular crystal explosives, successful application of DFT has largely failed due to the inability of current exchange-correlation functionals to correctly describe intermolecular van der Waals (vdWs) forces. Despite this, we have discovered that even with no treatment of vdWs bonding, the AM05 functional and DFT based molecular dynamics (MD) could be used to study the properties of molecular crystals under compression. We have used DFT-MD to predict the unreacted Hugoniots for PETN and HNS and validated the results by comparison with crystalline and porous experimental data. Since we are also interested in applying DFT methods to study the equilibrium volume properties of explosives, we studied the nature of the vdWs bonding in pursuit of creating a new DFT functional capable of accurately describing equilibrium bonding of molecular crystals. In this report we discuss our results for computing shock Hugoniots of molecular crystals and also what was learned about the nature of bonding in these materials.
We have demonstrated the ability to control the microstructure of PETN films deposited using physical vapor deposition by altering the interface between the film and substrate. Evolution of surface morphology, average density, and surface roughness with film thickness were characterized using surface profilometry and scanning electron microscopy. While films on all of the substrates investigated showed a trend toward a lower average density with increasing film thickness, there were significant variations in density, pore size, and surface morphology in films deposited on different substrates.
In an effort to better understand the structural changes occurring during hydrogen loading of erbium target materials, we have performed in situ D{sub 2} loading of erbium metal (powder) at temperature (450 C) with simultaneous neutron diffraction analysis. This experiment tracked the conversion of Er metal to the {alpha} erbium deuteride (solid-solution) phase and then into the {beta} (fluorite) phase. Complete conversion to ErD{sub 2.0} was accomplished at 10 Torr D{sub 2} pressure with deuterium fully occupying the tetrahedral sites in the fluorite lattice.
Density Functional Theory (DFT) has over the last few years emerged as an indispensable tool for understanding the behavior of matter under extreme conditions. DFT based molecular dynamics simulations (MD) have for example confirmed experimental findings for shocked deuterium, enabled the first experimental evidence for a triple point in carbon above 850 GPa, and amended experimental data for constructing a global equation of state (EOS) for water, carrying implications for planetary physics. The ability to perform high-fidelity calculations is even more important for cases where experiments are impossible to perform, dangerous, and/or prohibitively expensive. For solid explosives, and other molecular crystals, similar success has been severely hampered by an inability of describing the materials at equilibrium. The binding mechanism of molecular crystals (van der Waals forces) is not well described within traditional DFT. Among widely used exchange-correlation functionals, neither LDA nor PBE balances the strong intra-molecular chemical bonding and the weak inter-molecular attraction, resulting in incorrect equilibrium density, negatively affecting the construction of EOS for undetonated high explosives. We are exploring a way of bypassing this problem by using the new Armiento-Mattsson 2005 (AM05) exchange-correlation functional. The AM05 functional is highly accurate for a wide range of solids, in particular in compression. In addition, AM05 does not include any van der Waals attraction, which can be advantageous compared to other functionals: Correcting for a fictitious van der Waals like attraction with unknown origin can be harder than correcting for a complete absence of all types of van der Waals attraction. We will show examples from other materials systems where van der Waals attraction plays a key role, where this scheme has worked well, and discuss preliminary results for molecular crystals and explosives.
The difficulty of calculating the ambient properties of molecular crystals, such as the explosive PETN, has long hampered much needed computational investigations of these materials. One reason for the shortcomings is that the exchange-correlation functionals available for Density Functional Theory (DFT) based calculations do not correctly describe the weak intermolecular van der Waals' forces present in molecular crystals. However, this weak interaction also poses other challenges for the computational schemes used. We will discuss these issues in the context of calculations of lattice constants and structure of PETN with a number of different functionals, and also discuss if these limitations can be circumvented for studies at non-ambient conditions.
In an effort to better understand the structural changes occurring during hydrogen loading of erbium target materials, we have performed D{sub 2} loading of erbium metal (powder) with simultaneous neutron diffraction analysis. This experiment tracked the conversion of Er metal to the {alpha} erbium deuteride (solid-solution) phase and then on to the {beta} (fluorite) phase. Complete conversion to ErD{sub 2.0} was accomplished at 10 Torr D{sub 2} pressure with deuterium fully occupying the tetrahedral sites in the fluorite lattice. Increased D{sub 2} pressure (up to 500 Torr at 450 C) revealed {approx}10 % deuterium occupation of the octahedral sites. Subsequent vacuum pumping of the sample at 450 C removed octahedral site occupancy while maintaining tetrahedral deuterium occupancy, thereby yielding stoichiometric ErD{sub 2.0} {beta} phase.
Sandia National Laboratories has cradle to grave responsibility for all neutron generators in the US nuclear weapons stockpile. As such, much research effort is exerted to develop a comprehensive understanding of all the major components of a neutron generator. One of the key components is the tritium containing target. The target is a thin metal tritide film. Sandia's research into metal tritides began in the early 1960's with a collaboration with the Denver Research Institute (DRI) and continues to this day with a major in house research effort. This document is an attempt to briefly summarize what is known about the aging of erbium tritide and to review the major publications conducted at Sandia in FY 07. First, a review of our knowledge of helium in erbium tritide will be presented. Second, executive summaries of the six major SAND reports regarding neutron tube targets published in FY07 by Department 2735, the Applied Science and Technology Maturation Department, and research partners are presented.
In this work, we examine the formation energies of interstitials in semiconductors obtained with four different pure functionals. Explicitely we investigate three silicon self-interstitials. All functionals give the same trend among those interstitials; the lowest energy being for formation of the <110>-split, somewhat higher energy for the formation of the hexagonal interstitial, while highest energy among the three is obtained for the meta-stable tetragonal configuration. However, the value for the formation energy for a specific interstitial differs substantially in calculations using different functionals. It is shown that the main contribution to these differences is stemming from the functionals different surface intrinsic errors. We also discuss the puzzle that the values obtained with the surface intrisic error free AM05 functional (Armiento and Mattsson, Phys. Rev. B 72, 085108 (2006)) gives values substantially lower than Quantum Monte Carlo results
Heterogeneous chemical reactions occurring at a gas/surface interface are fundamental in a variety of important applications, such as combustion, catalysis, chemical vapor deposition and plasma processing. Detailed simulation of these processes may involve complex, coupled fluid flow, heat transfer, gas-phase chemistry, in addition to heterogeneous reaction chemistry. This report documents the Surfkin program, which simulates the kinetics of heterogeneous chemical reactions. The program is designed for use with the Chemkin and Surface Chemkin (heterogeneous chemistry) programs. It calculates time-dependent or steady state surface site fractions and bulk-species production/destruction rates. The surface temperature may be specified as a function of time to simulate a temperature-programmed desorption experiment, for example. This report serves as a user's manual for the program, explaining the required input and format of the output. Two detailed example problems are included to further illustrate the use of this program.