Members of the Nuclear Criticality Safety (NCS) Program at Sandia National Laboratories (SNL) have updated the suite of benchmark problems developed to validate MCNP6 Version 2.0 for use in NCS applications. The updated NCS benchmark suite adds approximately 600 new benchmarks and includes peer review of all input files by two different NCS engineers (or one NCS engineer and one candidate NCS engineer). As with the originally released benchmark suite, the updated suite covers a broad range of fissile material types, material forms, moderators, reflectors, and neutron energy spectra. The benchmark suite provides a basis to establish a bias and bias uncertainty for use in NCS analyses at SNL.
Complex challenges across Sandia National Laboratories' (SNL) mission areas underscore the need for systems level thinking, resulting in a better understanding of the organizational work systems and environments in which our hardware and software will be used. SNL researchers have successfully used Activity Theory (AT) as a framework to clarify work systems, informing product design, delivery, acceptance, and use. To increase familiarity with AT, a working group assembled to select key resources on the topic and generate an annotated bibliography. The resources in this bibliography are arranged in six categories: 1) An introduction to AT; 2) Advanced readings in AT; 3) AT and human computer interaction (HCI); 4) Methodological resources for practitioners; 5) Case studies; and 6) Related frameworks that have been used to study work systems. This annotated bibliography is expected to improve the reader's understanding of AT and enable more efficient and effective application of it.
For 2D-temperature monitoring applications, a variant of EIT (Electrical Impedance Tomography) is evaluated computationally in this work. Literature examples of poor sensor performance in the center of the 2D domains away from the side electrodes motivated this study which seeks to overcome some of the previously noted shortcomings. In particular, the use of ‘sensing skins’ with novel tailored baseline conductivities were examined using the EIDORS package for EIT. It was found that the best approach for detecting a hot spot depends on several factors such as the current injection (stimulation) patterns, the measurement patterns, and the reconstruction algorithms. For a well-performing combination of these factors, tailored baseline conductivities were assessed and compared to the baseline uniform conductivity. It was discovered that for some EIT applications, a tailored distribution needs to be smooth and that sudden changes in the conductivity gradients should be avoided. Still, the benefits in terms of improved EIT performance were small for conditions for which the EIT measurements had been ‘optimized’ for the uniform baseline case. Within the limited scope of this study, only two specific cases showed benefits from tailored distributions. For one case, a smooth tailored distribution with increased baseline conductivity in the center provided a better separation of two centrally located hot spots. For another case, a smooth tailored distribution with reduced conductivity in the center provided better estimates of the magnitudes of two hot spots near the center of the sensing skin.
The DOE Office of Electricity views sodium batteries as a priority in pursuing a safe, resilient, and reliable grid. Improvements in solid-state electrolytes are key to realizing the potential of these large-scale batteries. NaSICON structure consists of SiO4 or PO4 tetrahedra sharing common corners with ZrO6 octahedra. Structure forms “tunnels” in three dimensions that can transport interstitial sodium ion. 3D structure provides higher ionic conductivity than other conductors (β’’-alumina), particularly at low temperature. Lower temperature (cheaper) processing compared to β’’-alumina. Our objective was to identify fundamental structure-processing-property relationships in NaSICON solid electrolytes to inform design for use in sodium batteries. In this work, the mechanical properties of NaSICON sodium ion conductors are affected by sodium conduction. Electrochemical cycling can alter modulus and hardness in NaSICON. Excessive cycling can lead to secondary phases and/or dendrite formation that change mechanical properties in NaSICON. Mechanical and electrochemical properties can be correlated with topographical features to further inform design decisions
Sanz-Matias, Ana; Roychoudhury, Subhayan; Feng, Xuefei; Yang, Feipeng; Cheng, Kao L.; Zavadil, Kevin R.; Guo, Jinghua; Prendergast, David
Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode-electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. To this end, we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.
The prevalence of COVID-19 is shaped by behavioral responses to recommendations and warnings. Available information on the disease determines the population’s perception of danger and thus its behavior; this information changes dynamically, and different sources may report conflicting information. We study the feedback between disease, information, and stay-at-home behavior using a hybrid agent-based-system dynamics model that incorporates evolving trust in sources of information. We use this model to investigate how divergent reporting and conflicting information can alter the trajectory of a public health crisis. The model shows that divergent reporting not only alters disease prevalence over time, but also increases polarization of the population’s behaviors and trust in different sources of information.
Grignard reagents of the general formula RMgX (X = Cl-, Br-, I-) have been utilized in various chemistries for over 100 years. We report that replacing the halide in a Grignard reagent with a reactive borohydride anion adds a new synthetic dimension for these influential compounds. We synthesized the series RMgBH4 (R = Et, n-Bu, Ph, Bn) and characterized the reactivity toward both organic and inorganic molecules. Using butylmagnesium borohydride (BuMgBH4) as an exemplar, we demonstrate that these compounds possess unique reactivity due to the presence of reducing borohydride groups, resulting in tandem reactivity with organic amides/esters to generate secondary and primary alcohols. Molecular dynamics simulations indicate the stability of BuMgBH4 is comparable to that of Mg(BH4)2 + MgBu2, validating the Schlenk equilibrium in borohydride Grignard compounds. Metadynamics simulations confirm that the equilibrium is kinetically accessible through solvent-mediated processes. BuMgBH4 also reacts with CO2 and NH3, revealing potential uses for CO2 utilization and as a mixed-anion metal borohydride/amide precursor.
The solution processability of ionogel solid electrolytes has recently garnered attention in the Li-ion battery community as a means to address the interface and fabrication issues commonly associated with most solid electrolytes. However, the trapped ionic liquid (ILE) component has hindered the electrochemical performance. Herein, we present a process to tune the properties by replacing the ILE in a silica-based ionogel after fabrication with a liquid component befitting the desired application. Electrochemical cycling under various conditions showcases gels containing different liquid components incorporated into LiFePO4 (LFP)/gel/Li cells: high power (455 W kg-1 at a 1 C discharge) systems using carbonates, low temperatures (-40 °C) using ethers, or high temperatures (100 °C) using ionic liquids. Fabrication of additive-manufactured cells utilizing the exchanged carbonate-based system is demonstrated in a planar LFP/Li4Ti5O12 (LTO) system, where a marked improvement over an ionogel is found in terms of rate capability, capacity, and cycle stability (118 vs 41 mA h g-1 at C/4). This process represents a promising route to create a separator-less cell, potentially in complex architectures, where the electrolyte properties can be facilely tuned to meet the required conditions for a wide range of battery chemistries while maintaining a uniform electrolyte access throughout cast electrodes.
The long-standing problem of predicting the electronic structure of matter on ultra-large scales (beyond 100,000 atoms) is solved with machine learning.
Artificial solid electrolyte interphases have provided a path to improved cycle life for high energy density, next-generation anodes like lithium metal. Although long cycle life is necessary for widespread implementation, understanding and mitigating the effects of aging and self-discharge are also required. Here, we investigate several coating materials and their role in calendar life aging of lithium. We find that the oxide coatings are electronically passivating whereas the LiF coating slows charge transfer kinetics. Furthermore, the Coulombic loss during self-discharge measurements improves with the oxide layers and worsens with the LiF layer. It is found that none of the coatings create a continuous conformal, electronically passivating layer on top of the deposited lithium nor are they likely to distribute evenly through a porous deposit, suggesting that none of the materials are acting as an artificial solid electrolyte interphase. Instead, they likely alter performance through modulating lithium nucleation and growth.
Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.
Understanding the selectivity of metal–organic frameworks (MOFs) to complex acid gas streams will enable their use in industrial applications. Herein, ab initio molecular dynamic simulations (AIMD) were used to simulate ternary gas mixtures (H2O-NO2-SO2) in rare earth 2,5-dihydroxyterephthalic acid (RE-DOBDC) MOFs. Stronger H2O gas-metal binding arose from thermal vibrations in the MOF sterically hindering access of SO2 and NO2 molecules to the metal sites. Gas-gas and gas-linker interactions within the MOF framework resulted in the formation of multiple secondary gas species including HONO, HNO2, NOSO, and HNO3−. Four gas adsorption sites were identified along with a new de-protonation reaction mechanism not observable through experiment. This study not only provides valuable information on competitive gas binding energies in the MOF, it also provides important chemical insights into transient chemical reactions and mechanisms.
Cryogenic plasma focused ion beam (PFIB) electron microscopy analysis is applied to visualizing ex situ (surface industrial) and in situ (subsurface geologic) carbonation products, to advance understanding of carbonation kinetics. Ex situ carbonation is investigated using NIST fly ash standard #2689 exposed to aqueous sodium bicarbonate solutions for brief periods of time. In situ carbonation pathways are investigated using volcanic flood basalt samples from Schaef et al. (2010) exposed to aqueous CO2 solutions by them. The fly ash reaction products at room temperature show small amounts of incipient carbonation, with calcite apparently forming via surface nucleation. Reaction products at 75° C show beginning stages of an iron carbonate phase, e.g., siderite or ankerite, common phases in subsurface carbon sequestration environments. This may suggest an alternative to calcite in carbonation low calcium-bearing fly ashes. Flood basalt carbonation reactions show distinct zonation with high calcium and calcium-magnesium bearing zones alternating with high iron-bearing zones. The calcium-magnesium zones are notable with occurrence of localized pore space. Oscillatory zoning in carbonate minerals is distinctly associated with far-from-equilibrium conditions where local chemical environments fluctuate via a coupling of reaction with transport. The high porosity zones may reflect a precursor phase (e.g., aragonite) with higher molar volume that then “ripens” to the high-Mg calcite phase-plus-porosity. These observations reveal that carbonation can proceed with evolving local chemical environments, formation and disappearance of metastable phases, and evolving reactive surface areas. Together this work shows that future application of cryo-PFIB in carbonation studies would provide advanced understanding of kinetic mechanisms for optimizing industrial-scale and commercial-scale applications.
Laser-induced photoemission of electrons offers opportunities to trigger and control plasmas and discharges. However, the underlying mechanisms are not sufficiently characterized to be fully utilized. Photoemission is highly nonlinear, achieved through multiphoton absorption, above threshold ionization, photo-assisted tunneling, etc., where the dominant process depends on the work function of the material, photon energy and associated fields, surface heating, background fields, etc. To characterize the effects of photoemission on breakdown, breakdown experiments were performed and interpreted using a 0D plasma discharge circuit model and quantum model of photoemission.
Monitoring of cooling tower performance in a nuclear reactor facility is necessary to ensure safe operation; however, instrumentation for measuring performance characteristics can be difficult to install and may malfunction or break down over long duration experiments. This paper describes employing a thermodynamic approach to quantify cooling tower performance, the Merkel model, which requires only five parameters, namely, inlet water temperature, outlet water temperature, liquid mass flowrate, gas mass flowrate, and wet bulb temperature. Using this model, a general method to determine cooling tower operation for a nuclear reactor was developed in situations when neither the outlet water temperature nor gas mass flowrate are available, the former being a critical piece of information to bound the Merkel integral. Furthermore, when multiple cooling tower cells are used in parallel (as would be in the case of large-scale cooling operations), only the average outlet temperature of the cooling system is used as feedback for fan speed control, increasing the difficulty of obtaining the outlet water temperature for each cell. To address these shortcomings, this paper describes a method to obtain individual cell outlet water temperatures for mechanical forced-air cooling towers via parametric analysis and optimization. In this method, the outlet water temperature for an individual cooling tower cell is acquired as a function of the liquid-to-gas ratio (L/G). Leveraging the tight tolerance on the average outlet water temperature, an error function is generated to describe the deviation of the parameterized L/G to the highly controlled average outlet temperature. The method was able to determine the gas flowrate at rated conditions to be within 3.9% from that obtained from the manufacturer’s specification, while the average error for the four individual cooling cell outlet water temperatures were 1.6 °C, -0.5 °C, -1.0 °C, and 0.3 °C.
Reliable climate predictions are important for making robust decisions in response to the changing climate. This project aims to reduce mis-modeling uncertainties arising from the representation of the land-atmosphere coupling in the Energy Exascale Earth System Model (E3SM) by using a machine learning approach. This approach will use an encoder-decoder architecture to represent the information that is developed in the land model and given to the atmosphere model. The simulated data will be taken from the E3SM simulation. However, the incorporation of observed data into the simulated dataset reduces mis-modeling uncertainties.
This memo serves as an initial deficiency study of current foam modeling approaches, to determine where model changes and/or improvements can be made to capture the phenomena of receding foam. We are looking for feedback on our approach and suggestions from interested internal customers.
Lithium-ion batteries (LIBs) have revolutionized our society in many respects, and we are expecting even more favorable changes in our lifestyles with newer battery technologies. In pursuing such eligible batteries, nanophase materials play some important roles in LIBs and beyond technologies. Stimulated by their beneficial effects of nanophase materials, we initiated this Focus. Excitingly, this Focus collects 13 excellent original research and review articles related to the applications of nanophase materials in various rechargeable batteries, ranging from nanostructured electrode materials, nanoscale interface tailoring, novel separators, computational calculations, and advanced characterizations.
The International Atomic Energy Agency (IAEA) applies safeguards to nuclear facilities that are not operating, including those undergoing decommissioning, and the IAEA’s effort in this area is both considerable and increasing. Specifically, the IAEA Department of Safeguards’ Division of Concepts and Planning (SGCP-003: Safeguards Approaches) identified an R&D need to “Develop safeguards implementation guidelines for facilities under decommissioning and safeguards concepts for post-accident facilities under decommissioning”. Nuclear facilities undergoing decommissioning are not exempt from safeguards agreements between the IAEA and Host State, and, accordingly, the requirement for verification of no diversion of nuclear material and detection of undeclared activities at decommissioned facilities remain even after facility shutdown. However, the effort required to meet safeguards objectives diminishes as nuclear material and essential equipment are removed during the decommissioning process which shifts the emphasis from verification of ever-diminishing fissile or source material inventories to verification of changes in facility design and equipment operability.
The New Mexico Small Business Assistance Program (NMSBA) has once again paired with Optical Radio Communications Technology (ORC Tech). A New Mexico startup Limited Liability Company (LLC), with Sandia National Laboratories (SNL) Engineers at the Sensors and Textiles Innovatively Tailored for Complex, High-Efficiency Detection (STITCHED) laboratory, to aid in the development of an ultra-passive, portable, deployable wireless signal booster technology.
This project matured a new understanding (a “modern synthesis”) of the structure and evolution of science and technology. It created an understanding and framework for how Sandia National Labs, the Department of Energy, and the nation, might improve their research productivity, with significant ramifications on national security and economic competitiveness.
This report documents our experience constructing a numerical method for the collisional Boltzmann equation that is capable of accurately capturing the collisionless through strongly collisional limits. We explore three different functional representations and present a detailed account of a numerical method based on a spatially dependent Gaussian mixture model (GMM). The Kullback-Leibler divergence is used as a closeness measure and various expectation maximization (EM) solution algorithms are implemented to find a compact representation in velocity space for distribution functions that exhibit significant non-Maxwellian character. We discuss issues that appear with this representation over a range of Knudsen numbers for a prototypical test problem and demonstrate that the strongly collisional limit recovers a solution to Euler's equations. Looking forward, this approach is broadly applicable to the non-relativistic and relativistic collisional Vlasov equations.