Worldwide growth in electric vehicle use is prompting new installations of private and public electric vehicle supply equipment (EVSE). EVSE devices support the electrification of the transportation industry but also represent a linchpin for power systems and transportation infras-tructures. Cybersecurity researchers have recently identified several vulnerabilities that exist in EVSE devices, communications to electric vehicles (EVs), and upstream services, such as EVSE vendor cloud services, third party systems, and grid operators. The potential impact of attacks on these systems stretches from localized, relatively minor effects to long-term national disruptions. Fortunately, there is a strong and expanding collection of information technology (IT) and operational technology (OT) cybersecurity best practices that may be applied to the EVSE environment to secure this equipment. In this paper, we survey publicly disclosed EVSE vulnerabilities, the impact of EV charger cyberattacks, and proposed security protections for EV charging technologies.
Rezaul Karim, Mohammad; Narasimhachary, Santosh; Radaelli, Francesco; Amann, Christian; Dayal, Kaushik; Silling, Stewart A.; Germann, Timothy C.
We present a computational study and framework that allows us to study and understand the crack nucleation process from forging flaws. Forging flaws may be present in large steel rotor components commonly used for rotating power generation equipment including gas turbines, electrical generators, and steam turbines. The service life of these components is often limited by crack nucleation and subsequent growth from such forging flaws, which frequently exhibit themselves as non-metallic oxide inclusions. The fatigue crack growth process can be described by established engineering fracture mechanics methods. However, the initial crack nucleation process from a forging flaw is challenging for traditional engineering methods to quantify as it depends on the details of the flaw, including flaw morphology. We adopt the peridynamics method to describe and study this crack nucleation process. For a specific industrial gas turbine rotor steel, we present how we integrate and fit commonly known base material property data such as elastic properties, yield strength, and S-N curves, as well as fatigue crack growth data into a peridynamic model. The obtained model is then utilized in a series of high-performance two-dimensional peridynamic simulations to study the crack nucleation process from forging flaws for ambient and elevated temperatures in a rectangular simulation cell specimen. The simulations reveal an initial local nucleation at multiple small oxide inclusions followed by micro-crack propagation, arrest, coalescence, and eventual emergence of a dominant micro-crack that governs the crack nucleation process. The dependence on temperature and density of oxide inclusions of both the details of the microscopic processes and cycles to crack nucleation is also observed. The results are compared with fatigue experiments performed with specimens containing forging flaws of the same rotor steel.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user's guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.
For decades, it has been observed that the commonly used Borgnakke-Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke-Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke-Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke-Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). The consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.
To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. The novel SM proposed here is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.
Electrochemical characteristics and semiconducting behavior of additively manufactured electron beam melted (EBM) and wrought (WR) Ti–6Al–4V (Ti-G5) are compared in Ringer’s physiological solution. X-ray diffraction (XRD) and field emission scanning electron microscopy (FE-SEM) confirmed the α + β structure of the tested materials, with two different microstructure types of “bimodal” and “basket-weave” for WR and EBM, respectively. Potentiodynamic polarization (PDP) revealed that the corrosion current density for EBM (icorr = 0.27 ± 0.06 μA cm−2) is less than the WR (icorr = 0.70 ± 0.05 μA cm−2). Moreover, potentiostatic polarization (PS) that was employed to form the passive layers at three different potentials of 300, 500, and 700 mVAg/AgCl, showed that the passive films on the EBM sample are thinner. This finding was confirmed by electrochemical impedance spectroscopy (EIS). Furthermore, through Mott–Schottky (M–S) analysis, donor densities on WR passive films were found to be ~ 1.5 times larger than EBM. Although PS and EIS confirmed that the passive layer on EBM is thinner, it provides higher corrosion resistance than WR. The passive layer on both samples were found to have n-type characteristics with a duplex structure. Graphical abstract: [Figure not available: see fulltext.]
A projection-based reduced order model (pROM) methodology has been developed for transient heat transfer problems involving coupled conduction and enclosure radiation. The approach was demonstrated on two test problems of varying complexity. The reduced order models demonstrated substantial speedups (up to 185×) relative to the full order model with good accuracy (less than 3% L∞ error). An attractive feature of pROMs is that there is a natural error indicator for the ROM solution: the final residual norm at each time-step of the converged ROM solution. Using example test cases, we discuss how to interpret this error indicator to assess the accuracy of the ROM solution. The approach shows promise for many-query applications, such as uncertainty quantification and optimization. The reduced computational cost of the ROM relative to the full-order model (FOM) can enable the analysis of larger and more complex systems as well as the exploration of larger parameter spaces.
Garces, Milton A.; Bowman, Daniel B.; Zeiler, Cleat; Christe, Anthony; Yoshiyama, Tyler; Williams, Brian; Colet, Meritxell; Takazawa, Samuel; Popenhagen, Sarah
A smartphone plummeted from a stratospheric height of 36 km, providing a near-real-time record of its rapid descent and ground impact. An app recorded and streamed useful internal multi-sensor data at high sample rates. Signal fusion with external and internal sensor systems permitted a more detailed reconstruction of the Skyfall chronology, including its descent speed, rotation rate, and impact deceleration. Our results reinforce the potential of smartphones as an agile and versatile geophysical data collection system for environmental and disaster monitoring IoT applications. We discuss mobile environmental sensing capabilities and present a flexible data model to record and stream signals of interest. The Skyfall case study can be used as a guide to smartphone signal processing methods that are transportable to other hardware platforms and operating systems.
Response to elongational flow is fundamental to soft matter and directly impacts new developments in a broad range of technologies form polymer processing and microfluidics to controlled flow in biosystems. Of particular significance are the effects of elongational flow on self-assembled systems where the interactions between the fundamental building blocks control their adaptation. Here we probe the effects of associating groups on the structure and dynamics of linear polymer melts in uniaxial elongation using molecular dynamics simulations. We study model polymers with randomly incorporated backbone associations with interaction strengths varying from 1kBT to 10kBT. These associating groups drive the formation of clusters in equilibrium with an average size that increases with interaction strength. Flow drives these clusters to continuously break and reform as chains stretch. These flow-driven cluster dynamics drive a qualitative transition in polymer elongation dynamics from homogeneous to nanoscale localized yield and cavitation as the association strength increases.
This article describes a calculation of the spontaneous emission limited linewidth of a semiconductor laser consisting of hybrid or heterogeneously integrated, silicon and III–V intracavity components. Central to the approach are a) description of the multi-element laser cavity in terms of composite laser/free-space eigenmodes, b) use of multimode laser theory to treat mode competition and multiwave mixing, and c) incorporation of quantum-optical contributions to account for spontaneous emission effects. Application of the model is illustrated for the case of linewidth narrowing in an InAs quantum-dot laser coupled to a high- (Formula presented.) SiN cavity.
This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.8 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.
Errors in quantum logic gates are usually modeled by quantum process matrices (CPTP maps). But process matrices can be opaque and unwieldy. We show how to transform the process matrix of a gate into an error generator that represents the same information more usefully. We construct a basis of simple and physically intuitive elementary error generators, classify them, and show how to represent the error generator of any gate as a mixture of elementary error generators with various rates. Finally, we show how to build a large variety of reduced models for gate errors by combining elementary error generators and/or entire subsectors of generator space. We conclude with a few examples of reduced models, including one with just 9N2 parameters that describes almost all commonly predicted errors on an N-qubit processor.
Babiniec, Sean M.; Reinholz, Emilee L.; Coker, Eric N.; Larsen, Marin E.
Intumescent materials are in wide use as protective coatings in fire protection or thermal management applications. These materials undergo chemical reactions occurring from approximately 300°C to 900°C, which outgas and expand the material, providing an appreciable increase in insulative performance. However, the complicated chemical mechanisms and large changes in materials properties complicate the incorporation of these materials into predictive thermal models. This document serves to outline the thermochemical characterization of select intumescent materials, the extraction of relevant parameters, and the incorporation of these parameters into the ChemEQ reaction model implemented in Aria. This work was performed in 2016 and documented in a draft SAND report in March 2017. In 2022, the draft SAND report was discovered and put through R&A.
Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User's Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user's notes and of course the material in the open literature. This page intentionally left blank.
This report documents the progress and current results of the MELCOR spent fuel cask input model. The MELCOR model is being developed to investigate aerosol transport and deposition given the aerosol physical phenomena models within MELCOR. To perform the analyses, a general portrayal of the MAGNASTOR® cask system has been employed; however, this system was selected based on available information to provide a reasonable representation of a spent fuel cask. The analytical results are not intended to characterize the performance of the MAGNASTOR® cask. Instead, the provided results are intended to enhance our general understanding of the aerosol behavior within casks and the validity of current models. The current model efforts are being performed to investigate hypothetical UO2 release from failed fuel pins within a spent fuel cask. The existing MELCOR model of the MAGNASTOR® cask system has been adapted to permit future comparative analyses with the GOTHIC representation of the MAGNASTOR® cask. To support this comparison, the PNNL model characteristics that are unrelated to the aerosol modeling were applied to the MELCOR model. These characteristics included improved comparability of the axial fidelity, total spent fuel power, fuel pin axial power profile, and heat losses from cannister. The thermal-hydraulic solutions are improved within the capability of the MELCOR code and will permit better overall agreement with the GOTHIC results. Detailed results are presented on the thermal-hydraulic analysis of the MELCOR cask as well as characterization of UO2 aerosol dispersion and deposition within the cask.
The “Decel” platform at Sandia National Laboratories investigates the Richtmyer–Meshkov instability (RMI) in converging geometry under high energy density conditions [Knapp et al., Phys. Plasmas 27, 092707 (2020)]. In Decel, the Z machine magnetically implodes a cylindrical beryllium liner filled with liquid deuterium, launching a converging shock toward an on-axis beryllium rod machined with sinusoidal perturbations. The passage of the shock deposits vorticity along the Be/D2 interface, causing the perturbations to grow. In this paper, we present platform improvements along with recent experimental results. To improve the stability of the imploding liner to the magneto Rayleigh–Taylor instability, we modified its acceleration history by shortening the Z electrical current pulse. Next, we introduce a “split rod” configuration that allows two axial modes to be fielded simultaneously in different axial locations along the rod, doubling our data per experiment. We then demonstrate that asymmetric slots in the return current structure modify the magnetic drive pressure on the surface of the liner, advancing the evolution on one side of the rod by multiple ns compared to its 180° counterpart. This effectively enables two snapshots of the instability at different stages of evolution per radiograph with small deviations of the cross-sectional profile of the rod from the circular. Using this platform, we acquired RMI data at 272 and 157 μm wavelengths during the single shock stage. Finally, we demonstrate the utility of these data for benchmarking simulations by comparing calculations using ALEGRA MHD and RageRunner.
On March 30th and 31st, 2022, the University of Texas at Austin (UT) Office of the Vice President for Research (OVPR) hosted Sandia National Laboratories (Sandia) for “Sandia Day at UT Austin” to understand the status of the strategic partnership and explore opportunities for partnership growth. The event brought together more than 115 UT and Sandia participants including executive leadership, researchers, faculty, staff, and students. Sandia Day primarily consisted of a half-day leadership meeting, a research poster session and networking event, and three break-out sessions focused on strategic priority areas: Microelectronics, Energy and Climate Security, and High-Performance and Edge Computing. Appendix A contains the full Sandia Day agenda. Additional meetings and workshops (adjunct meetings) were held in conjunction with Sandia Day to maximize partnership exploration. Adjunct meetings were Hypersonics, Decarbonization, Disinformation, and Battery Workshops. A summary of Sandia Day events, sessions, and meetings follows.
Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Results showed trends & errors consistent with related studies. Area validation metric provides way to quantify model form uncertainty. AVM shows that more work could be done to understand how model form uncertainty varies with mesh resolution. There is a possible atypical use of MAVM on time-series data. AVM shows mismatch between predicted flame height and experimental value less sensitive to variations in mixture fraction than temperature. Mismatch about experimental value also more symmetric for mixture fraction. Our analysis showed that mixture fraction is preferable for this application.
Elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene−propylene−diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.
Various laboratory-focused tools and methodologies for completing a safety risk assessment have been published, yet few similar resources to address chemical security exist.Herein, we describe a chemical security risk assessment case study at a university in a developing country.In this case study, we demonstrate a chemical security risk assessment for a university chemistry department, using an original inventory of 645 entries which was condensed to 295 chemicals after removing duplicates and erroneous entries.We then prioritized to highlight 83 chemicals of interest based on hazardous or dual-use properties that could lead to unacceptable consequences.We further refined to a list of 34 high-risk chemicals that required action, 48 chemicals that may need further justification and consideration for additional protection, and 1 chemical that did not need further consideration for additional protection.
Jackson, Stuart L.; Hinshelwood, David D.; Kaiser, Eric R.; Swanekamp, Stephen B.; Richardson, Andrew S.; Schumer, Joseph W.; Johnson, Michael J.; Foster, John E.; Durot, Christopher J.
Time-dependent visible and near-ultraviolet emission spectroscopy was used to track the presence of excited and ionized states induced as an electron beam transited a cavity filled with air at low pressures. An electron beam was produced in vacuum using a Febetron pulsed-power generator modified to produce a peak voltage of 100 kV, a peak current of 4.5 kA, and a pulse width of 100 ns. The beam electrons then passed through thin anode and pressure foils and transited a cavity filled with dry air at either 0.1 Torr or 1 Torr, exciting and ionizing the air along the way. The spectral measurements were combined with information from electrical, interferometric, and imaging diagnostics to reveal different ionization dynamics for the two cases under study.
A well-posed physics-based compact model for a three-terminal silicon–oxide–nitride–oxide–silicon (SONOS) synaptic circuit element is presented for use by neuromorphic circuit/system engineers. Based on technology computer aided design (TCAD) simulations of a SONOS device, the model contains a nonvolatile memristor with the state variable QM representing the memristor charge under the gate of the three-terminal element. By incorporating the exponential dependence of the memristance on QM and the applied bias V for the gate, the compact model agrees quantitatively with the results from TCAD simulations as well as experimental measurements for the drain current. The compact model is implemented through VerilogA in the circuit simulation package Cadence Spectre and reproduces the experimental training behavior for the source–drain conductance of a SONOS device after applying writing pulses ranging from –12 V to +11 V, with an accuracy higher than 90%.
We report pure-rotational N2-N2, N2-air, and O2-air S-branch linewidths for temperatures of 80-200 K by measuring the time-dependent decay of rotational Raman coherences in an isentropic free-jet expansion from a sonic nozzle. We recorded pure-rotational hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) spectra along the axial centerline of the underexpanded jet, within the barrel shock region upstream of the Mach disk. The dephasing of the pure-rotational Raman coherence was monitored using probe-time-delay scans at different axial positions in the jet, corresponding to varying local temperatures and pressures. The local temperature was obtained by fitting CARS spectra acquired at zero probe time delay, where the impact of collisions was minimal. The measured decay of each available Raman transition was fit to a dephasing constant and corrected for the local pressure, which was obtained from the CARS-measured static temperature and thermodynamic relationships for isentropic expansion from the known stagnation state. Nitrogen self-broadened transitions decayed more rapidly than those broadened in air for all temperatures, corresponding to higher Raman linewidths. In general, the measured S-branch linewidths deviated significantly in absolute and relative magnitudes from those predicted by extrapolating the modified exponential gap model to low temperatures. The temperature dependence of the Raman linewidth for each measured rotational state of nitrogen (J ≤ 10) and oxygen (N ≤ 11) was fit to a temperature-dependent power law over the measurable temperature domain (80-200 K) and extrapolated to both higher rotational states and room temperature. The measured and modeled low-temperature linewidth data provided here will aid low temperature gas-phase pressure measurements with fs/ps CARS.
Inverse prediction models have commonly been developed to handle scalar data from physical experiments. However, it is not uncommon for data to be collected in functional form. When data are collected in functional form, it must be aggregated to fit the form of traditional methods, which often results in a loss of information. For expensive experiments, this loss of information can be costly. In this study, we introduce the functional inverse prediction (FIP) framework, a general approach which uses the full information in functional response data to provide inverse predictions with probabilistic prediction uncertainties obtained with the bootstrap. The FIP framework is a general methodology that can be modified by practitioners to accommodate many different applications and types of data. We demonstrate the framework, highlighting points of flexibility, with a simulation example and applications to weather data and to nuclear forensics. Results show how functional models can improve the accuracy and precision of predictions.
Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.
Michelsen, Hope A.; Campbell, Matthew F.; Tran, Ich C.; Johansson, K.O.; Schrader, Paul E.; Bambha, Ray B.; Hammons, Joshua A.; Schaible, Eric; Zhu, Chenhui; Van Buuren, Anthony
We have developed a strategy for distinguishing between small-angle X-ray scattering (SAXS) from gas-phase species and newly formed nanoparticles in mixed gas- and particle-phase reacting flows. This methodology explicitly accounts for temperature-dependent scattering from gases. We measured SAXS in situ in a sooting linear laminar partially premixed co-flow ethylene/air diffusion flame. The scattering signal demonstrates a downward curvature as a function of the momentum transfer (q) at q values of 0.2-0.57 Å-1. The q-dependent curvature is consistent with the Debye equation and the independent-atom model for gas-phase scattering. This behavior can also be modeled using the Guinier approximation and could be characterized as a Guinier knee for gas-phase scattering. The Guinier functional form can be fit to the scattering signal in this q range without a priori knowledge of the gas-phase composition, enabling estimation of the gas-phase contribution to the scattering signal while accounting for changes in the gas-phase composition and temperature. We coupled the SAXS measurements with in situ temperature measurements using coherent anti-Stokes Raman spectroscopy. This approach to characterizing the gas-phase SAXS signal provides a physical basis for distinguishing among the contributions to the scattering signal from the instrument function, flame gases, and nanoparticles. The results are particularly important for the analysis of the SAXS signal in the q range associated with particles in the size range of 1-6 nm.
Li, Chang; Shyamsunder, Abhinandan; Hoane, Alexis G.; Long, Daniel M.; Kwok, Chun Y.; Kotula, Paul G.; Zavadil, Kevin R.; Gewirth, Andrew A.; Nazar, Linda F.
Aqueous zinc-metal batteries are plagued by poor Zn reversibility owing to zinc dendrite and layered double hydroxide (LDH) formation. Here, we introduce a novel additive—N,N-dimethylformamidium trifluoromethanesulfonate (DOTf)—in a low-cost aqueous electrolyte that can very effectively address these issues. The initial water-assisted dissociation of DOTf into triflic superacid creates a robust nanostructured solid-electrolyte interface (SEI)—revealed by operando spectroscopy and cryomicroscopy—which excludes water and enables dense Zn deposition. We demonstrate excellent Zn plating/stripping in a Zn||Cu asymmetric cell for more than 3,500 cycles. Furthermore, near 100% CE is realized at a combined high current density of 4 mA cm−2 and an areal capacity of 4 mAh cm−2 over long-term cycling. Zn||Zn0.25V2O5·nH2O full cells retain ∼83% of their capacity after 1,000 cycles with mass-limited Zn anodes. By restricting the depth of discharge, the cathodes exhibit less proton intercalation and LDH formation with an extended lifetime of 2,000 cycles.
We present an in-depth study of metal-insulator interfaces within granular metal (GM) films and correlate their interfacial interactions with structural and electrical transport properties. Nominally 100 nm thick GM films of Co and Mo dispersed within yttria-stabilized zirconia (YSZ), with volumetric metal fractions (φ) from 0.2-0.8, were grown by radio frequency co-sputtering from individual metal and YSZ targets. Scanning transmission electron microscopy and DC transport measurements find that the resulting metal islands are well-defined with 1.7-2.6 nm average diameters and percolation thresholds between φ = 0.4-0.5. The room temperature conductivities for the φ = 0.2 samples are several orders of magnitude larger than previously-reported for GMs. X-ray photoemission spectroscopy indicates both oxygen vacancy formation within the YSZ and band-bending at metal-insulator interfaces. The higher-than-predicted conductivity is largely attributed to these interface interactions. In agreement with recent theory, interactions that reduce the change in conductivity across the metal-insulator interface are seen to prevent sharp conductivity drops when the metal concentration decreases below the percolation threshold. These interface interactions help interpret the broad range of conductivities reported throughout the literature and can be used to tune the conductivities of future GMs.
Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.
We present a new analysis methodology that allows for the self-consistent integration of multiple diagnostics including nuclear measurements, x-ray imaging, and x-ray power detectors to determine the primary stagnation parameters, such as temperature, pressure, stagnation volume, and mix fraction in magnetized liner inertial fusion (MagLIF) experiments. The analysis uses a simplified model of the stagnation plasma in conjunction with a Bayesian inference framework to determine the most probable configuration that describes the experimental observations while simultaneously revealing the principal uncertainties in the analysis. We validate the approach by using a range of tests including analytic and three-dimensional MHD models. An ensemble of MagLIF experiments is analyzed, and the generalized Lawson criterion χ is estimated for all experiments.
Dereka, Bogdan; Lewis, Nicholas H.C.; Zhang, Yong; Hahn, Nathan H.; Keim, Jonathan H.; Snyder, Scott A.; Maginn, Edward J.; Tokmakoff, Andrei
Understanding the mechanisms of charge transport in batteries is important for the rational design of new electrolyte formulations. Persistent questions about ion transport mechanisms in battery electrolytes are often framed in terms of vehicular diffusion by persistent ion-solvent complexes versus structural diffusion through the breaking and reformation of ion-solvent contacts, i.e., solvent exchange events. Ultrafast two-dimensional (2D) IR spectroscopy can probe exchange processes directly via the evolution of the cross-peaks on picosecond time scales. However, vibrational energy transfer in the absence of solvent exchange gives rise to the same spectral signatures, hiding the desired processes. We employ 2D IR on solvent resonances of a mixture of acetonitrile isotopologues to differentiate chemical exchange and energy-transfer dynamics in a comprehensive series of Li+, Mg2+, Zn2+, Ca2+, and Ba2+ bis(trifluoromethylsulfonyl)imide electrolytes from the dilute to the superconcentrated regime. No exchange phenomena occur within at least 100 ps, regardless of the ion identity, salt concentration, and presence of water. All of the observed spectral dynamics originate from the intermolecular energy transfer. These results place the lower experimental boundary on the ion-solvent residence times to several hundred picoseconds, much slower than previously suggested. With the help of MD simulations and conductivity measurements on the Li+ and Zn2+ systems, we discuss these results as a continuum of vehicular and structural modalities that vary with concentration and emphasize the importance of collective electrolyte motions to ion transport. These results hold broadly applicable to many battery-relevant ions and solvents.
Using the thesis of W.R. Nolan (cite) as a guide, Cobalt Iron (CoFe) powders were reacted with 0.1 wt.% and 0.2 wt.% phosphoric acid in a 20:1 ratio of acetone to phosphoric acid. The powders were then dried at room temperature. The resulting phosphate coated CoFe was then mixed with 0.75 wt.% of the lubricant N,N' ethylene bis-stearamide (trade name: Acrawax C) and hot pressed to form a consolidated soft magnetic material referred to as CoFeP. With an avenue of synthesis for CoFeP determined, a proper amount of stock was synthesized for continuous “brick” production. While under current optimization, these 1x1 mm magnetic bricks will ultimately be placed and secured along the inside wall of each MK Magnetics transformer core by an appropriate CoFeP dispersed epoxy. As of now each brick has been produced though a pressing and annealing process via square 2x2 cm die. Before a brick is made a pressure calculation is run to ensure the dies maximum operating pressure is not exceeded. Figure 1. ensures the user’s safety by showing that the tons-on-ram required for a 2x2 cm square die to reach 760 MPa is below the point of die failure.
Over the past 50 years, the Renewable Energy Program at Sandia has advanced research in the field with a focus on three key goals; 1) reduce the cost, 2) improve resilience and reliability and, 3) decrease the regulatory burden of renewable energy. Sandia’s expertise, coupled with the Village of Questa’s expanding renewable energy portfolio, presents the opportunity to deploy the Labs’ deep science and engineering capabilities towards the energy goals of KCEC and the Village of Questa. Preliminary research efforts by Sandia technical staff has broadly identified early opportunities for further research, development, and demonstration in the emerging renewable energy segment of agrivoltaics. Agrivoltaics is an emerging and promising area of photovoltaics which entails land use considerations as well as concerns regarding landscape transformation, biodiversity, and ecosystem well-being. In recent years, agrivoltaics systems have been the subject of numerous studies due to their potential in the food-energy (and water) nexus. This document is a preliminary evaluation of the projects performance opportunities of agrivoltaics as a renewable energy technology strategy in the region of Questa, NM.
As concerns about the illicit movement of radioactive materials across international borders increase, so too has the need for increased protection of those borders both foreign and domestic. The challenge is not only to detect hidden radioactive materials, but also to distinguish them from legitimate radionuclides such as radio-pharmaceuticals that are often transported across borders and shipped throughout a country. With more than 600 U.S. border sites to protect, screening imported radioactive material requires a careful balance of high throughput and high search efficiency. However, these requirements are difficult to meet as rapid screening operations leave less time for radiation detectors to efficiently evaluate materials. In support of border security, Sandia developed organic glass scintillators.
The nuclear spins of low-density implanted Ga atoms in Ge are interesting candidates for solid state-based qubits. To date, activation studies of implanted Ga in Ge have focused on high densities. Here, we extend activation studies into the low-density regime. We use spreading resistance profiling and secondary ion mass spectrometry to derive electrical activation of Ga ions implanted into Ge as a function of the rapid thermal anneal temperature and implant density. We show that for our implant conditions, the activation is best for anneal temperatures between 400 and 650 °C with a maximum activation of 69% at the highest fluence. Below 400 °C, remaining implant damage results in defects that act as superfluous carriers, and above 650 °C, surface roughening and loss of Ga ions are observed. The activation increased monotonically from 10% to 69% as the implant fluence increased from 6 × 10 10 to 6 × 10 12 cm-2. The results provide thermal anneal conditions to be used for initial studies of using low-density Ga atoms in Ge as nuclear spin qubits.
A fast matrix-free formulation of the control volume finite element method is presented, requiring much less memory and computational work than previous efforts. The method is implemented and evaluated as a solver for low-Mach flow, including the evaluation of a preconditioning strategy for the pressure Poisson equation. The efficiency and scaling with polynomial order is evaluated on simple turbulent flows of interest, with appropriate solution quality metrics, and compared with a reference node-centered finite volume discretization. For a turbulent channel flow test, we show improvement in computational work for a given accuracy with the high-order scheme. The performance on a GPU accelerated platform is also investigated, with benefit shown for the matrix-free discretization.
This report is the revised (Revision 8) Task F specification for DECOVALEX-2023. Task F is a comparison of the models and methods used in deep geologic repository performance assessment. The task proposes to develop a reference case for a mined repository in a fractured crystalline host rock and a reference case for a mined repository in a salt formation. Teams may choose to participate in the comparison for either or both of the reference cases. For each reference case, a common set of conceptual models and parameters describing features, events, and processes that impact performance will be given, and teams will be responsible for determining how best to implement and couple the models. The comparison will be conducted in stages, beginning with a comparison of key outputs of individual process models, followed by a comparison of a single deterministic simulation of the full reference case, and moving on to uncertainty propagation and uncertainty and sensitivity analysis. This report provides background information, a summary of the proposed reference cases, and a staged plan for the analysis.
Particle fragmentation influences thermochemical energy conversion processes in different ways and is of significance in energy generation technologies. Different reactive material formulations trigger varied thermal response in extreme environments such as high velocity impact. This study investigated optical thermal response of powder gun launched intermetallic (Al:Zr) and thermite (Al:MoO3) projectiles using pyrometry and thermography. Projectiles were launched at 1250 m/s into an air-filled chamber and impacted a steel witness plate to create a dust explosion. The pyrometer was configured to measure temperatures directly at the point of impact, while the thermographic system measured temperatures throughout the explosion chamber. Results show that impact temperatures ranged between 3500 and 4000 K, but that the dynamics of energy conversion were different for the intermetallic and thermite projectiles. The intermetallic exhibited secondary reactions due to fragmented debris impacting the walls of the chamber. The thermite exhibited greater gas generation, propelling the debris field, and producing a more stochastic response with faster spreading and dissipation of thermal energy. Unique features such as microexplosions within fragmented particles were also analyzed. While both reactive materials produce similar temperatures, their mechanisms of energy conversion and release are different, indicating the potential of these materials for different ballistic applications.
Easily measured metrics that could assign quantifiable values to coating batches for quality control have started to be developed. High-density is an attribute of quality films. Increased density results in harder, more wear resistant coatings in inert and humid environments. Denser films are more resistant to oxidation from aging, limiting the severity and depth of oxide into the coating. Future work includes using metrics for quality. The next step is to develop in house deposition capabilities to develop process-structure relationships.
Precise temperature determination is a significant challenge in extreme environments of dynamic compression studies. How can radiance measurements taken in high-pressure shock experiments constrain temperature in a meaningful and physically consistent way? Experiments maintaining sample compression against a transparent window can be tailored to present a uniform measurement area with uncertain spectral emissivity. We compare several methods to analyze radiance collected at multiple wavelengths, applying statistical methods and physical principles to improve temperature inference. With proper radiance collection and analysis, dynamic temperature uncertainties become comparable to thermomechanical ambiguities of the emitting surface.
Mannion, Owen M.; Crilly, Aiden J.; Forrest, Chad J.; Appelbe, Brian D.; Betti, Riccardo; Glebov, V.Y.; Gopalaswamy, Varchas; Knauer, James P.; Mohamed, Zaarah L.; Stoeckl, Christian; Chittenden, Jerry P.; Regan, Sean P.
The apparent ion temperature and mean velocity of the dense deuterium tritium fuel layer of an inertial confinement fusion target near peak compression have been measured using backscatter neutron spectroscopy. The average isotropic residual kinetic energy of the dense deuterium tritium fuel is estimated using the mean velocity measurement to be ~103 J across an ensemble of experiments. The apparent ion-temperature measurements from high-implosion velocity experiments are larger than expected from radiation-hydrodynamic simulations and are consistent with enhanced levels of shell decompression. Furthermore, these results suggest that high-mode instabilities may saturate the scaling of implosion performance with the implosion velocity for laser-direct-drive implosions.
Design and construction of a particle-to-sCO2 heat exchanger test station is described in this paper. The purpose of this test station is to make steady-state measurements of thermal performance with sCO2 as the working fluid. While the test station was initially constructed to test a 20 kW heat exchanger developed under the Gen3 Particle pilot plant (G3P3) project, it also was designed to accommodate testing of other heat exchanger configurations. Improvements for this test station design is based on lessons learned from prior heat exchanger testing. Maximum pressure and temperature ratings are based on the desire to use primarily stainless steel in the construction to reduce cost and lead time of components. Construction of the test station was completed and commissioning and initial testing took place during the October to November 2020 timeframe.
Particle-based heat transfer materials used in concentrating solar power systems benefit from gravity-fed arrangements such as vertically integrated components inside the receiver tower which can eliminate the need for conveyance machinery. However, the amount of particles required for commercial scale systems near 100 MWe can require towers with very thick walls that must be built with high-strength concrete. Cost models for particle-based receiver towers with internal particle storage are being developed in this work and compared to well-established cost models that have been used to estimate tower costs for molten salt systems with external storage tanks. New cost models were developed to accommodate the high-temperature applications required for CSP. Further research is needed to directly compare costs between tower-integrated and external storage. For now, a method is proposed to superimpose increased storage costs with existing molten salt CSP towers. For instances where suitable materials are unavailable or do not meet the structural requirements, ground based storage bins must be used in concert with mechanical conveyance systems. Ground based storage vessels have been shown to be consistent with low thermal energy storage cost and heat loss goals. Ground based storage vessels are well-established in industry.
This investigation explores thermal-fluid flow phenomena in a proportional flow control valve (FCV) within a 2 in. ID high-temperature piping transport system. The FCVs are critical components to ensure flexible nominal operation of a 2 MWth concentrating solar power (CSP) pilot-scale system in present development at Sandia National Laboratories (SNL). A computational fluid dynamics (CFD) / finite element analysis (FEA) model was developed in ANSYS that investigates multifluid phase-change transport within various sections of an FCV to explore plating and subsequent thermal-mechanical stress challenges that can exist with operations as high as 730°C. Results from the thermal-fluid model in development suggest salt vapor phase change in the N2 gas purge lines as low as approximately 476°C, which can have a negative impact on valve reliability.
This project extends Sandia's experience in Light Detection And Ranging (LiDAR) to gain an understanding of the abilities and limits of using 3D laser scanning to capture the relative canting angles between heliostat mirror surfaces in 3D space to an accuracy sufficient to measure canting errors. To the authors' knowledge, this approach has never been developed or implemented for this purpose. The goal is to be able to automatically perform a 3D scan, retrieve the data, and use computational geometry and a-priori mechanical knowledge of the heliostats (facet arrangement and size) to filter and isolate the facets, and fit planar models to the facet surfaces. FARO FocusS70 laser range scanners are used, which provide a dense data coverage of the scan area in the form of a 3D point-cloud. Each point has the 3D coordinates of the surface position illuminated by the device as it scans the laser beam over an area, both in azimuth and elevation. These scans can contain millions of points in total. The initial plan was to primarily use the back side of the heliostat to capture the mirror (the back side being opaque). It was not expected to capture high-quality data from the reflective front side. The discovery that the front side did, indeed, yield surface data was surprising. This is a function of the soiling, or collected dust, on the mirror surface. Typical point counts on the mirror facets are seen to be between 10k - 100k points per facet, depending on the facet area and the scan point density. By collecting facet surface points, the data can be used to calculate an individual planar fit per facet, the normals of which correlate directly with the facet pointing angle. Comparisons with neighboring facets yield the canting angles. The process includes software which automatically: 1) controls the LiDAR scanner and downloads the resultant scan data, 2) isolates the heliostat data from the full scan, 3) filters the points associated with each individual facet, and 4) calculates the planar fit and relative canting angles for each facet. The goal of this work has been to develop this system to measure heliostat canting errors to less than 0.25 mrad accuracy, with processing time under 5 minutes per heliostat. A future goal is to place this or a comparable sensor on an autonomous platform, along with the software system, to collect and analyze heliostats in the field for tracking and canting errors in real time. This work complements Sandia's strategic thrust in autonomy for CSP collector systems.
This paper summarizes the evolution of the Gen 3 Particle Pilot Plant (G3P3) receiver design with the goal of reducing heat losses and increasing thermal efficiencies. New features that were investigated included aperture covers and shrouds, active airflow, multistage catch-and-release devices (stairs), and optimization of receiver cavity geometry. Simulations and ground-based testing showed that a reduced receiver volume and aperture shroud could reduce advective heat losses by ∼40 - 50%, and stairs could increase opacity and reduce backwall temperatures. The reduced volume receiver and stairs were selected for on-sun testing, and receiver efficiencies up to 80 - 90% were achieved in the current test campaign. The receiver thermal efficiency generally increased as a function of incident power and particle mass flow rates. In addition, particle outlet temperatures were maintained to within ±10 °C of a prescribed setpoint temperature up to ∼700 °C using a PID controller that adjusted the particle mass flow rate into the receiver in response to the measured particle outlet temperatures.
Padawer-Blatt, Aviv; Ducatel, Jordan; Bogan, Alex; Gaudreau, Louis; Zawadzki, Piotr; Austing, D.G.; Sachrajda, Andrew S.; Studenikin, Sergei; Tracy, Lisa A.; Reno, John; Hargett, Terry H.
Difference in g factors in multidot structures can form the basis of dot-selective spin manipulation under global microwave irradiation. Here employing electric dipole spin resonance facilitated by strong spin-orbit interaction (SOI), we observe differences in the extracted values of the single-hole effective g factors of the constituent quantum dots of a GaAs/AlGaAs double quantum dot device at the level of ~ 5 %–10%. We examine the continuous change in the hole g factor with electrical detuning over a wide range of interdot tunnel couplings and for different out-of-plane magnetic fields. The observed tendency of the quantum dot effective g factors to steadily increase on decreasing the interdot coupling or on increasing the magnetic field is attributed to the impact on the SOI of changing the dot confinement potential and heavy-hole light-hole mixing.
In the California Industrial General Permit (IGP) 2014-0057-DWQ for storm water monitoring, effective July 1, 2015, there are 21 contaminants that have been assigned NAL (Numeric Action Level) values, both annual and instantaneous. For annual NALs, an exceedance occurs when the average of all analytical results from all samples taken at a facility during a reporting year for a given parameter exceeds an annual NAL value listed in Table 2 of the General Permit. For instantaneous maximum NALs, an exceedance occurs when two or more analytical results from samples taken for any parameter within a reporting year exceed the instantaneous maximum NAL value (for TSS and O&G), or are outside of the instantaneous maximum NAL range (for pH) listed in Table 2.
For systems that require complete metallic enclosures (e.g., containment buildings for nuclear reactors), it is impossible to access interior sensors and equipment using standard electromagnetic techniques. A viable way to communicate and supply power through metallic barriers is the use of elastic waves and ultrasonic transducers, introducing several design challenges that must be addressed. Specifically, the use of multiple communication channels on the same enclosure introduces an additional mechanism for signal crosstalk between channels: guided waves propagating in the barrier between channels. This work numerically and experimentally investigates a machined phononic crystal to block MHz Lamb wave propagation between ultrasonic communication channels, greatly reducing wave propagation and the resulting crosstalk voltage. Blind grooves are machined into one or both sides of a metallic barrier to introduce a periodic unit cell, greatly altering the guided wave dispersion in the barrier. Numerical simulations are used to determine a set of groove geometries for testing, and experiments were performed to characterize the wave-blocking performance of each design. The best-performing design was tested using piezoelectric transducers bonded to the barrier, showing a 14.4 dB reduction in crosstalk voltage. The proposed periodic grooving method is a promising technique for completely isolating ultrasonic power/data transfer systems operating in a narrow frequency range.
Previous efforts determined a set of calibrated, optimal model parameter values for Reynolds-averaged Navier–Stokes (RANS) simulations of a compressible jet in crossflow (JIC) using a $k–ε$ turbulence model. These parameters were derived by comparing simulation results to particle image velocimetry (PIV) data of a complementary JIC experiment under a limited set of flow conditions. Here, a $k–ε$ model using both nominal and calibrated parameters is validated against PIV data acquired from a much wider variety of JIC cases, including a realistic flight vehicle. The results from the simulations using the calibrated model parameters showed considerable improvements over those using the nominal values, even for cases that were not used in the calibration procedure that defined the optimal parameters. This improvement is demonstrated using a number of quality metrics that test the spatial alignment of the jet core, the magnitudes of multiple flow variables, and the location and strengths of vortices in the counter-rotating vortex cores on the PIV planes. These results suggest that the calibrated parameters have applicability well outside the specific flow case used in defining them and that with the right model parameters, RANS solutions for the JIC can be improved significantly over those obtained from the nominal model.
The study of charge carrier transport at nanoscale electrical contacts is crucial for the development of next-generation electronics. In this study, we review recent modeling efforts on quantum tunneling, current crowding, and contact resistance across electrical interfaces with nanometer scale dimensions. A generalized self-consistent model for quantum tunneling induced electron transport in metal–insulator–metal (MIM) junctions is summarized. Rectification of a dissimilar MIM junction is reviewed. A modified two-dimensional (2D) transmission line model is used to investigate the effects of spatially varying specific contact resistivity along the contact length. The model is applied to various types of electrical contacts, including ohmic contacts, MIM junction based tunneling contacts, and 2D-material-based Schottky contacts. Roughness engineering is recently proposed to offer a possible paradigm for reducing the contact resistance of 2D-material-based electrical contacts. Contact interface engineering, which can mitigate current crowding near electrical contacts by spatially designing the interface layer thickness or properties, without requiring an additional material or component, is briefly reviewed. Tunneling engineering is suggested to eliminate severe current crowding in highly conductive ohmic contacts by introducing a thin tunneling layer or gap between the contact members. Furthermore, unsolved problems and challenges are also discussed.
Epoxy underfills can be implemented in electronic packaging to enhance solder joint reliability of surface mounted components. However, it is important for an engineer to have a failure criterion that can be used for failure predictions and redesign of electronic assemblies. For this study, data from epoxy bond failure in mock electronic part assemblies were correlated to finite element analyses to predict adhesive failure initiation. Experiments were performed to determine failure loads for various loading locations and nonlinear viscoelastic analyses were performed for the same loading locations to determine a maximum principal strain failure parameter. Predictions showed that a maximum principal strain failure parameter defined from one test could be used as an indicator of adhesive failure of an epoxy bond undergoing other modes of loading. Failure initiation predictions matched experimental data using a maximum principal strain failure parameter for an epoxy bond undergoing mixed modes of loading for both unfilled and alumina oxide filled 828DEA epoxy. Such experimental setup is deemed appropriate for future epoxy testing.
This project will provide scientific understanding needed to design, optimize, and calibrate the next generations of off-road diesel engines that comply with increasingly stringent pollutant emission regulations while achieving thermal efficiencies exceeding 50%.
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. (2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. (3) Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.
Accelerated aging studies of β CL-20 thin films deposited on glass surfaces in different environments (N2, air, air/water) were conducted. Samples were analyzed with attenuated total reflectance infrared (ATR-IR) spectroscopy. Spectral features of molecular lattice inclusions in CL-20 films aged in an air/water environment were observed. The features occurred after β CL-20 polymorph transformation to α CL-20 and were accompanied by the appearance of lattice water peaks. To assist ATR-IR peak assignment, density functional theory studies were performed, and IR spectra of α CL-20 lattice inclusions of small molecules that were identified as degradation products of CL-20 were computed. Simulated spectra of NO2, HNCO, N2O, and CO2 incorporated into partially hydrated α CL-20 matched the experimental spectrum of β CL-20 thin films aged in air/water.