Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Photovoltaic (PV) performance is affected by reversible and irreversible losses. These can typically be mitigated through responsive and proactive operations and maintenance (O&M) activities. However, to generate profit, the cost of O&M must be lower than the value of the recovered electricity. This value depends both on the amount of recovered energy and on the electricity prices, which can vary significantly over time in spot markets. The present work investigates the impact of the electricity price variability on the PV profitability and on the related O&M activities in Italy, Portugal, and Spain. It is found that the PV revenues varied by 1.6 × to 1.8 × within the investigated countries in the last 5 years. Moreover, forecasts predict higher average prices in the current decade compared to the previous one. These will increase the future PV revenues by up to 60% by 2030 compared to their 2015–2020 mean values. These higher revenues will make more funds available for better maintenance and for higher quality components, potentially leading to even higher energy yield and profits. Linearly growing or constant price assumptions cannot fully reproduce these expected price trends. Furthermore, significant price fluctuations can lead to unexpected scenarios and alter the predictions.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The use of simple models for response prediction of building structures is preferred in earthquake engineering for risk evaluations at regional scales, as they make computational studies more feasible. The primary impediment in their gainful use presently is the lack of viable methods for quantifying (and reducing upon) the modeling errors/uncertainties they bear. This study presents a Bayesian calibration method wherein the modeling error is embedded into the parameters of the model. The method is specifically described for coupled shear-flexural beam models here, but it can be applied to any parametric surrogate model. The major benefit the method offers is the ability to consider the modeling uncertainty in the forward prediction of any degree-of-freedom or composite response regardless of the data used in calibration. The method is extensively verified using two synthetic examples. In the first example, the beam model is calibrated to represent a similar beam model but with enforced modeling errors. In the second example, the beam model is used to represent the detailed finite element model of a 52-story building. Both examples show the capability of the proposed solution to provide realistic uncertainty estimation around the mean prediction.

Battery storage systems are increasingly being installed at photovoltaic (PV) sites to address supply-demand balancing needs. Although there is some understanding of costs associated with PV operations and maintenance (O&M), costs associated with emerging technologies such as PV plus storage lack details about the specific systems and/or activities that contribute to the cost values. This study aims to address this gap by exploring the specific factors and drivers contributing to utility-scale PV plus storage systems (UPVS) O&M activities costs, including how technology selection, data collection, and related and ongoing challenges. Specifically, we used semi-structured interviews and questionnaires to collect information and insights from utility-scale owners and operators. Data was collected from 14 semi-structured interviews and questionnaires representing 51.1 MW with 64.1 MWh of installed battery storage capacity within the United States (U.S.). Differences in degradation rate, expected life cycle, and capital costs are observed across different storage technologies. Most O&M activities at UPVS related to correcting under-performance. Fires and venting issues are leading safety concerns, and owner operators have installed additional systems to mitigate these issues. There are ongoing O&M challenges due the lack of storage-specific performance metrics as well as poor vendor reliability and parts availability. Insights from this work will improve our understanding of O&M consideration at PV plus storage sites.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A mathematical framework is provided for a substructuring-based domain decomposition (DD) approach for nonlocal problems that features interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation (PDE) problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal DD methods.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The recent discovery of bright, room-temperature, single photon emitters in GaN leads to an appealing alternative to diamond best single photon emitters given the widespread use and technological maturity of III-nitrides for optoelectronics (e.g. blue LEDs, lasers) and high-speed, high-power electronics. This discovery opens the door to on-chip and on-demand single photon sources integrated with detectors and electronics. Currently, little is known about the underlying defect structure nor is there a sense of how such an emitter might be controllably created. A detailed understanding of the origin of the SPEs in GaN and a path to deterministically introduce them is required. In this project, we develop new experimental capabilities to then investigate single photon emission from GaN nanowires and both GAN and AlN wafers. We ion implant our wafers with the ion implanted with our focused ion beam nanoimplantation capabilities at Sandia, to go beyond typical broad beam implantation and create single photon emitting defects with nanometer precision. We've created light emitting sources using Li⁺ and He⁺, but single photon emission has yet to be demonstrated. In parallel, we calculate the energy levels of defects and transition metal substitutions in GaN to gain a better understanding of the sources of single photon emission in GaN and AlN. The combined experimental and theoretical capabilities developed throughout this project will enable further investigation into the origins of single photon emission from defects in GaN, AlN, and other wide bandgap semiconductors.

Self-determination has been an on-going effort for Native American people and gained much traction with the passing of The Energy Policy Act of 2005, which included the Indian Tribal Energy Development and Self-Determination Act. Congress passed this act to assist Native American tribes and Alaska Native villages with planning, development, and assistance to achieve their energy goals. The Ute Mountain Ute Tribe (UMUT) has relied on oil and natural gas for economic support the last 70 years. Burning fossil fuels, along with oil and gas development, decreases the quality of air and leads to increased greenhouse gas emissions. Subsequently, the burning of fossil fuels to produce energy is now more costly than many renewable energy sources, including solar photovoltaic (PV) systems. Environmental stewardship, along with the need to maintain revenue generation, has led UMUT’s efforts to achieve energy self-determinism employing PV and exploring other technology. In the past, the tribe completed a 1 megawatt PV project near Towaoc, Colorado, which serves as a case study on the tribe’s energy goals: a future where renewables will dominate their energy landscape. This paper explores UMUT’s past and on-going efforts toward energy independence and how it relates to the broader landscape of Native American energy sovereignty.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Barium titanate (BTO) is a ferroelectric material used in capacitors because of its high bulk dielectric constant. However, the impact of the size of BTO on its dielectric constant is not yet fully understood and is highly contested. Here, we present an investigation into the dielectric constant of BTO nanoparticles with diameters ranging between 50 and 500 nm. BTO nanoparticles were incorporated into acrylonitrile butadiene styrene and injection molded into parallel plate capacitors, which were used to determine nanocomposite dielectric constants. The dielectric constants of BTO nanoparticles were obtained by combining experimental measurements with computational results from COMSOL simulations of ABS-matrix nanocomposites containing BTO. The dielectric constant of BTO was observed to be relatively constant at nanoparticle diameters as small as 200 nm but sharply declined at smaller nanoparticle sizes. These results will be useful in the development of improved energy storage and power conditioning systems utilizing BTO nanoparticles.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Grid operating security studies are typically employed to establish operating boundaries, ensuring secure and stable operation for a range of operation under NERC guidelines. However, if these boundaries are severely violated, existing system security margins will be largely unknown, as would be a secure incremental dispatch path to higher security margins while continuing to serve load. As an alternative to the use of complex optimizations over dynamic conditions, this work employs the use of machine learning to identify a sequence of secure state transitions which place the grid in a higher degree of operating security with greater static and dynamic stability margins. Several reinforcement learning solution methods were developed using deep learning neural networks, including Deep Q-learning, Mu-Zero, and the continuous algorithms Proximal Reinforcement Learning, and Advantage Actor Critic Learning. The work is demonstrated on a power grid with three control dimensions but can be scaled in size and dimensionality, which is the subject of ongoing research.

The wide variety of inverter control settings for solar photovoltaics (PV) causes the accurate knowledge of these settings to be difficult to obtain in practice. This paper addresses the problem of determining inverter reactive power control settings from net load advanced metering infrastructure (AMI) data. The estimation is first cast as fitting parameterized control curves. We argue for an intuitive and practical approach to preprocess the AMI data, which exposes the setting to be extracted. We then develop a more general approach with a data-driven reactive power disaggregation algorithm, reframing the problem as a maximum likelihood estimation for the native load reactive power. These methods form the first approach for reconstructing reactive power control settings of solar PV inverters from net load data. The constrained curve fitting algorithm is tested on 701 loads with behind-the-meter (BTM) PV systems with identical control settings. The settings are accurately reconstructed with mean absolute percentage errors between 0.425% and 2.870%. The disaggregation-based approach is then tested on 451 loads with variable BTM PV control settings. Different configurations of this algorithm reconstruct the PV inverter reactive power timeseries with root mean squared errors between 0.173 and 0.198 kVAR.

Abstract not provided.

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, natural gas, and autogas systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying leak sizes and rates, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 5.0. The most significant change for HyRAM+ version 5.0 from HyRAM+ version 4.1 is the ability to model blends of different fuels. HyRAM+ was previously only suitable for use with hydrogen, methane, or propane, with users having the ability to use methane as a proxy for natural gas and propane as a proxy for autogas/liquefied petroleum gas. In version 5.0, real natural gas or autogas compositions can be modeled as the fuel, or even blends of natural gas with hydrogen. These blends can be used in the standalone physics models, but not yet in the quantitative risk assessment mode of HyRAM+.

This manual describes the installation and use of the Xyce™ XDM Netlist Translator. XDM simplifies the translation of netlists generated by commercial circuit simulator tools into Xyce-compatible netlists. XDM currently supports translation from PSpice, HSPICE, and Spectre netlists into Xyce™ netlists.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document presents the facility-recommended characterization of the neutron, prompt gamma ray, and delayed gamma ray radiation fields in the University of Texas at Austin Nuclear Engineering Teaching Laboratory (NETL) TRIGA reactor for the beam port 1/5 free-field environment at the 128-inch location adjacent to the core centerline. The designation for this environment is NETL-FF-BP1/5-128-cca. The neutron, prompt gamma ray, and delayed gamma ray energy spectra, uncertainties, and covariance matrices are presented as well as radial and axial neutron and gamma ray fluence profiles within the experiment area of the cavity. Recommended constants are given to facilitate the conversion of various dosimetry readings into radiation metrics desired by experimenters. Representative pulse operations are presented with conversion examples.

Abstract not provided.

This report examines the localization of high frequency electromagnetic fields in general three-dimensional cavities along periodic paths between opposing sides of the cavity. The focus is on the case where the mirrors at the ends of the orbit are concave and have two different radii of curvature. The cases where these orbits lead to unstable localized modes are known as scars. The ellipsoidal coordinate system is utilized in the construction of the scarred modes. The field at the interior foci is examined as well as trigonometric projections along the periodic scarred ray path.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Niobium doped lead-tin-zirconate-titanate ceramics near the PZT 95/5 orthorhombic AFE – rhombohedral FE morphotropic phase boundary Pb1-0.5y(Zr0.865-xTixSn0.135)1-yNbyO3 were prepared according to a 22+1 factorial design with x = 0.05, 0.07 and y = 0.0155, 0.0195. The ceramics were prepared by a traditional solid-state synthesis route and sintered to near full density at 1250°C for 6 h. All compositions were ∼98% dense with no detectable secondary phases by X-ray diffraction. The ceramics exhibited equiaxed grains with intergranular porosity, and grain size was ∼5 µm, decreasing with niobium substitution. Compositions exhibited remnant polarization values of ∼32 µC/cm2, increasing with Ti substitution. Depolarization by the hydrostatic pressure induced FE-AFE phase transition was drastically affected by variation of the Ti and Nb substitution, increasing at a rate of 113 MPa /1% Ti and 21 MPa/1% Nb. Total depolarization output was insensitive to the change in Ti and Nb substitution, ∼32.8 µC/cm2 for the PSZT ceramics. The R3c-R3m and R3m-Pm3m phase transition temperatures on heating ranged from 90 to 105°C and 183 to 191°C, respectively. Ti substitution stabilized the R3c and R3m phases to higher temperatures, while Nb substitution stabilized the Pm3m phase to lower temperatures. Thermal hysteresis of the phase transitions was also observed in the ceramics, with transition temperature on cooling being as much as 10°C lower.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present a proof-of-concept demonstration of a narrow linewidth $^{87}$Rb magneto-optical trap (MOT) operating on the narrow linewidth $5S_{1/2}$ → $6P_{3/2}$ transition at 420 nm. We stabilized the absolute frequency of the 420 nm laser to an atomic transition in $^{87}$Rb and demonstrate a MOT using 420 nm light driving the $5S_{1/2}$, $F = 2$ → $6P_{3/2}, F' = 3$ transition. We then use tome-of-flight measurements to characterize the 420 nm MOT temperature, observing a minimum temperature of about $T^{(420)}_{horizontal}$ = 150μK and $T^{(420)}_{vertical}$ = 250μK before the opportunity to perform significant characterization and optimization. Although this temperature is significantly higher then the expected 420 nm Doppler cooling limit ($T_D^{(420)}$ ≈ 34 μK), these are already approaching the Doppler limit of a standard 780 nm MOT ($T_D^{(780)}$ ≈ 146 μK). We believe that with further optimization the Doppler cooling limit of ≈ 34 μK can be achieved. This initial result answers our key research question and demonstrates the viability of applying narrow linewidth laser cooling as a robust technique for future fieldable quantum sensors.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

On October 1, 2022, sound level measurements were taken at various locations throughout Kirtland Air Force Base (KAFB) and Southeastern Albuquerque. The purpose was to support sound propagation modeling predictions and sound regulations for public exposure during the detonation of an approximately 300-pound energetic experiment. Ground Zero was located on Range 7 of Sandia Thunder Range (06647). A total of 8 measurement locations were identified (e.g., 5 on KAFB and 3 in the Southeastern Albuquerque neighborhoods).

Abstract not provided.

Abstract not provided.

Abstract not provided.

Metal hydride hydrogen compression utilizes a reversible heat-driven interaction of a hydride-forming metal alloy with hydrogen gas. This paper reports on the development of a laboratory scale two-stage Metal Hydride Compressor (MHC) system with a feed pressure of 150 bar delivering high purity H2 gas at outlet pressures up to 875 bar. Stage 1 and stage 2 AB₂ metal hydrides are identified based on experimental characterization of the pressure-composition-temperature (PCT) behavior of candidate materials. The selected metal hydrides are each combined with expanded natural graphite, increasing the thermal conductivity of the composites by an order of magnitude. These composites are integrated in two compressor beds with internal heat exchangers that alternate between hydrogenation and dehydrogenation cycles by thermally cycling between 20 °C and 150 °C. The prototype compressor achieved compression of hydrogen from 150 bar to 700 bar with an average flow rate of 33.6 g/hr.

A collection of x-ray computed tomography scans of candy.

Abstract not provided.

Abstract not provided.

The goal of the project was to protect US critical infrastructure and improve energy security through technical analysis of the risk landscape presented by the anticipated massive deployment of interoperable EV chargers.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Geologic Disposal Safety Assessment Framework is a state-of-the-art simulation software toolkit for probabilistic post-closure performance assessment of systems for deep geologic disposal of nuclear waste developed by the United States Department of Energy. This paper presents a generic reference case and shows how it is being used to develop and demonstrate performance assessment methods within the Geologic Disposal Safety Assessment Framework that mitigate some of the challenges posed by high uncertainty and limited computational resources. Variance-based global sensitivity analysis is applied to assess the effects of spatial heterogeneity using graph-based summary measures for scalar and time-varying quantities of interest. Behavior of the system with respect to spatial heterogeneity is further investigated using ratios of water fluxes. This analysis shows that spatial heterogeneity is a dominant uncertainty in predictions of repository performance which can be identified in global sensitivity analysis using proxy variables derived from graph descriptions of discrete fracture networks. New quantities of interest defined using water fluxes proved useful for better understanding overall system behavior.

Accurate representation of environmental controllers of soil organic carbon (SOC) stocks in Earth System Model (ESM) land models could reduce uncertainties in future carbon–climate feedback projections. Using empirical relationships between environmental factors and SOC stocks to evaluate land models can help modelers understand prediction biases beyond what can be achieved with the observed SOC stocks alone. In this study, we used 31 observed environmental factors, field SOC observations (n = 6,213) from the continental United States, and two machine learning approaches (random forest [RF] and generalized additive modeling [GAM]) to (a) select important environmental predictors of SOC stocks, (b) derive empirical relationships between environmental factors and SOC stocks, and (c) use the derived relationships to predict SOC stocks and compare the prediction accuracy of simpler model developed with the machine learning predictions. Out of the 31 environmental factors we investigated, 12 were identified as important predictors of SOC stocks by the RF approach. In contrast, the GAM approach identified six (of those 12) environmental factors as important controllers of SOC stocks: potential evapotranspiration, normalized difference vegetation index, soil drainage condition, precipitation, elevation, and net primary productivity. The GAM approach showed minimal SOC predictive importance of the remaining six environmental factors identified by the RF approach. Our derived empirical relations produced comparable prediction accuracy to the GAM and RF approach using only a subset of environmental factors. The empirical relationships we derived using the GAM approach can serve as important benchmarks to evaluate environmental control representations of SOC stocks in ESMs, which could reduce uncertainty in predicting future carbon–climate feedbacks.

Abstract not provided.

A collection of x-ray computed tomography scans of candy.

Abstract not provided.

Abstract not provided.

The Sandia National Laboratories site sustainability plan and its associated DOE Sustainability Dashboard data entries encompass Sandia National Laboratories contributions toward meeting the DOE sustainability goals. This site sustainability plan fulfills the contractual requirement for National Technology & Engineering Solutions of Sandia, LLC, the management and operating contractor for Sandia National Laboratories, to deliver an annual sustainability plan to the DOE National Nuclear Security Administration Sandia Field Office.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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 eficiency is achieved as the number of processors grows.

Abstract not provided.

Abstract not provided.

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 SiO₄ or PO₄ tetrahedra sharing common corners with ZrO₆ 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

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 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.

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 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.

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.

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.

The long-standing problem of predicting the electronic structure of matter on ultra-large scales (beyond 100,000 atoms) is solved with machine learning.

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 CO₂ 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.

Abstract not provided.

Abstract not provided.

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 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.

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.

Surface flashover in vacuum imposes a substantial physical limit on modern, large-scale pulsed power. One of the ramifications is a minimum size requirement for new machines, which in itself becomes a hard barrier to the modernization and improvement of existing infrastructure. Pulsed power topologies require the physical mechanisms of both anode- and cathode-initiated flashover to be considered. Originally, the geometrical implications of field emission at the cathode triple junction (CTJ) motivated the usage of configurations that avoid electrons impinging on the insulating material. This will largely suppress the cathode-initiated flashover, which is best described by the secondary electron avalanche mechanism, gas desorption, and final breakdown in the desorbed gas. It depends on the cascade growth of a conducting plasma along the length of the insulator from the cathode. Mitigating the cathode-initiated flashover typically comes at the cost of a significant field enhancement at the anode triple junction (ATJ). In a typical implementation, the anode field may be three times higher than the cathode field for a given voltage, making the corresponding anode-initiated flashover much more common than otherwise. In the case of pulsed, anode-initiated flashover, experimental evidence suggests that charge is directly extracted from the insulator resulting in the insulator taking on a net positive charge advancing the anode potential. Along with accompanying gas desorption from the surface, the potential will then propagate from the anode toward the cathode until the effective length of the gap is sufficiently reduced to support flashover. The underlying physical mechanisms of cathode- and anode-directed flashover are discussed in light of previously gathered experimental data and recent experiments with pulsed, high-gradient, anode-initiated flashover.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Grey Zone Test Range (GZTR) social model operates as a piece of the overall GZTR modeling effort. It works in conjunction with supply models for resources, an electric grid model for power availability, and a traffic model for road congestion, as well as a general controller framework that allows external system effects. The social model functions as an aggregate model where the entire population of the city is divided into groups based on the Transportation Analysis Zones (TAZs), a common geospatial boundary present in all GZTR models. These groups will act as a singular community; each time step the state of the system around them will be assessed and then community will come up with a general plan of action that they will attempt to follow for the day. Additionally, they will track values for their general emotional state and memory about negative impacts in the near past.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Recent efforts at Sandia such as DataSEA are creating search engines that enable analysts to query the institution’s massive archive of simulation and experiment data. The benefit of this work is that analysts will be able to retrieve all historical information about a system component that the institution has amassed over the years and make better-informed decisions in current work. As DataSEA gains momentum, it faces multiple technical challenges relating to capacity storage. From a raw capacity perspective, data producers will rapidly overwhelm the system with massive amounts of data. From an accessibility perspective, analysts will expect to be able to retrieve any portion of the bulk data, from any system on the enterprise network. Sandia’s Institutional Computing is mitigating storage problems at the enterprise level by procuring new capacity storage systems that can be accessed from anywhere on the enterprise network. These systems use the simple storage service, or S3, API for data transfers. While S3 uses objects instead of files, users can access it from their desktops or Sandia’s high-performance computing (HPC) platforms. S3 is particularly well suited for bulk storage in DataSEA, as datasets can be decomposed into object that can be referenced and retrieved individually, as needed by an analyst. In this report we describe our experiences working with S3 storage and provide information about how developers can leverage Sandia’s current systems. We present performance results from two sets of experiments. First, we measure S3 throughput when exchanging data between four different HPC platforms and two different enterprise S3 storage systems on the Sandia Restricted Network (SRN). Second, we measure the performance of S3 when communicating with a custom-built Ceph storage system that was constructed from HPC components. Overall, while S3 storage is significantly slower than traditional HPC storage, it provides significant accessibility benefits that will be valuable for archiving and exploiting historical data. There are multiple opportunities that arise from this work, including enhancing DataSEA to leverage S3 for bulk storage and adding native S3 support to Sandia’s IOSS library.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This research effort examined the application of Nafion polymers in alcohol solvents as an anti-ice surface coating, as a mixture with hydrophilic polymers and freezing point depressant salt systems. Co-soluble systems of Nafion, polymer and salt were applied using dip coating methods to create smooth films for frost observation over a Peltier plate thermal system in ambient laboratory conditions. Cryo-DSC was applied to examine freezing events of the Nafion-surfactant mixtures, but the sensitivity of the measurement was insufficient to determine frost behavior. Collaborations with the Fog Chamber at Sandia-Albuquerque, and in environmental SAXS measurements with CINT-LANL were requested but were not able to be performed under the research duration. Since experimental characterization of these factors is difficult to achieve directly, computational modeling was used to guide the scientific basis for property improvement. Computational modeling was performed to improve understanding of the dynamic association between ionomer side groups and added molecules and deicing salts. The polyacrylic acid in water system was identified at the start of the project as a relevant system for exploring the effect of varying counterions on the properties of fully deprotonated polyacrylic acid (PAA) in the presence of water. Simulations were modeled with four different counterions, two monovalent counterions (K+ and Na+) and two divalent counterions (Ca2+ and Mg2+). The wt% of PAA in these systems was varied from ~10 to 80 wt% PAA for temperatures from 250K to 400K. In the second set of simulations, the interpenetration of water into a dry PAA film was studied for Na+ or Ca2+ counterions for temperatures between 300K and 400K. The result of this project is a sprayable Nafion film composite which resists ice nucleation at -20 °C for periods of greater than three hours. It is composed of Nafion polymer, hydrophilic polyethylene oxide polymer and CaCl2 anti-ice crosslinker. Durability and field performance properties remain to be determined.

As part of the development process, scaled testing of wave energy converter devices are necessary to prove a concept, study hydrodynamics, and validate control system approaches. Creating a low-cost, small, lightweight data acquisition system suitable for scaled testing is often a barrier for wave energy converter developers’ ability to test such devices. This paper outlines an open-source solution to these issues, which can be customized based on specific needs. This will help developers with limited resources along a path toward commercialization.

New approaches to preventing and treating infections, particularly of the respiratory tract, are needed. One promising strategy is to reconfigure microbial communities (microbiomes) within the host to improve defense against pathogens. Probiotics and prebiotics for gastrointestinal (GI) infections offer a template for success. We sought to develop comparable countermeasures for respiratory infections. First, we characterized interactions between the airway microbiome and a biodefense-related respiratory pathogen (Burkholderia thailandensis; Bt), using a mouse model of infection. Then, we recovered microbiome constituents from the airway and assessed their ability to re-colonize the airway and protect against respiratory Bt infection. We found that microbiome constituents belonging to Bacillus and related genuses frequently displayed colonization and anti-Bt activity. Comparative growth requirement profiling of these Bacillus strains vs Bt enabled identification of candidate prebiotics. This work serves as proof of concept for airway probiotics, as well as a strong foundation for development of airway prebiotics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The V31 containment vessel was procured by the US Army Recovered Chemical Material Directorate (RCMD) as a third - generation EDS containment vessel. It is the fifth EDS vessel to be fabricated under Code Case 2564 of the 2019 ASME Boiler and Pressure Vessel Code, which provides rules for the design of impulsively loaded vessels. The explosive rating for the vessel, based on the code case, is twenty-four (24) pounds TNT - equivalent for up to 1092 detonations. This report documents the results of explosive tests that were performed on the vessel at Sandia National Laboratories in Albuquerque, New Mexico to qualify the vessel for field operations use. There were three design basis configurations for qualification testing. Qualification test (1) consisted of a simulated M55 rocket motor and warhead assembly of 24lbs of Composition C-4 (30 lb TNT equivalent). This test was considered the maximum load case, based on modeling and simulation methods performed by Sandia prior to the vessel design phase. Qualification test (2) consisted of a regular, right circular cylinder, unitary charge, located central to the vessel interior of 19.2 lb of Composition C-4 (24 lb TNT equivalent). Qualification test (3) consisted of a 12-pack of regular, right circular cylinders of 2 lb each, distributed evenly inside the vessel (totaling 19.2 lb of C-4, or 24 lb TNT equivalent). All vessel acceptance criteria were met.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A collection of x-ray computed tomography scans of candy.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The purpose and scope of the viga span tables project for Rachel Wood Consulting (RWC) is focused on producing tabulated beam span tables for three species of wood vigas commonly used in New Mexico to allow producers, designers and builders to incorporate vigas into their building designs in a prescriptive manner similar to the span tables for sawn lumber incorporated into the International Residential Code (IRC) or the International Log Builders Association (ILBA) publication. The information provided in this report and the associated viga span tables also attempts to address and clarify questions raised by RWC during their review of the 2018 Los Alamos National Laboratory (LANL) New Mexico Small Business Assistance (NMSBA) program report by August Mosimann pertaining to span lengths, loading, deflection calculations, and log grading certification prior to submitting the span tables to the Construction Industries Division (CID) of New Mexico.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.