Publications

2021-2022 Storm Water Annual Report And Supporting Documents.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Crossed-field diodes (CFDs) are used in multiple high-power applications and are characterized by an applied magnetic field orthogonal to the electric field, induced by the applied voltage across the anode-cathode gap. In vacuum, the Hull cutoff magnetic field (HCMF) represents the maximum applied magnetic field for which an electron from the cathode can reach the anode. This study investigates the effects of non-vacuum conditions on electron trajectories by introducing electron mobility, which represents particle collisions. We used numerical solutions of the electron force law and particle-in-cell simulations (XPDP1) to assess electron motion for various electron mobilities. For magnetic fields above the HCMF in vacuum, reducing the electron mobility increases the time for an electron emitted from the cathode to reach the anode. Reducing mobility below 22 C s/kg eliminates the HCMF for any magnetic field, meaning that an emitted electron will always cross the gap. We derived the magnetic field, mobility, and electron transit time corresponding to this condition by solving for the condition when the electron velocity in the direction across the anode-cathode gap going to zero at the anode. A parametric study of these conditions using theory and XPDP1 is performed under different gap distances, voltages, and magnetic fields.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. The article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Lead–acid batteries are important to modern society because of their wide usage and low cost. The primary source for production of new lead–acid batteries is from recycling spent lead–acid batteries. In spent lead–acid batteries, lead is primarily present as lead pastes. In lead pastes, the dominant component is lead sulfate (PbSO4, mineral name anglesite) and lead oxide sulfate (PbO•PbSO4, mineral name lanarkite), which accounts for more than 60% of lead pastes. In the recycling process for lead–acid batteries, the desulphurization of lead sulfate is the key part to the overall process. In this work, the thermodynamic constraints for desulphurization via the hydrometallurgical route for recycling lead pastes are presented. The thermodynamic constraints are established according to the thermodynamic model that is applicable and important to recycling of lead pastes via hydrometallurgical routes in high ionic strength solutions that are expected to be in industrial processes. The thermodynamic database is based on the Pitzer equations for calculations of activity coefficients of aqueous species. The desulphurization of lead sulfates represented by PbSO4 can be achieved through the following routes. (1) conversion to lead oxalate in oxalate-bearing solutions; (2) conversion to lead monoxide in alkaline solutions; and (3) conversion to lead carbonate in carbonate solutions. Among the above three routes, the conversion to lead oxalate is environmentally friendly and has a strong thermodynamic driving force. Oxalate-bearing solutions such as oxalic acid and potassium oxalate solutions will provide high activities of oxalate that are many orders of magnitude higher than those required for conversion of anglesite or lanarkite to lead oxalate, in accordance with the thermodynamic model established for the oxalate system. An additional advantage of the oxalate conversion route is that no additional reductant is needed to reduce lead dioxide to lead oxide or lead sulfate, as there is a strong thermodynamic force to convert lead dioxide directly to lead oxalate. As lanarkite is an important sulfate-bearing phase in lead pastes, this study evaluates the solubility constant for lanarkite regarding the following reaction, based on the solubility data, PbO•PbSO4 + 2H+ ⇌ 2Pb2+ + SO42− + H2O(l).

Abstract not provided.

Thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio-posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens-is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300-430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Write-optimized data structures (WODS), offer the potential to keep up with cyberstream event rates and give sub-second query response for key items like IP addresses. These data structures organize logs as the events are observed. To work in a real-world environment and not fill up the disk, WODS must efficiently expire older events. As the basis for our research into organizing security monitoring data, we implemented a tool, called Diventi, to index IP addresses in connection logs using RocksDB (a write-optimized LSM tree). We extended Diventi to automatically expire data as part of the data structures’ normal operations. We guarantee that Diventi always tracks the N most recent events and tracks no more than N+ k events for a parameter k< N, while ensuring the index is opportunistically pruned. To test Diventi at scale in a controlled environment, we used anonymized traces of IP communications collected at SuperComputing 2019. We synthetically extended the 2.4 billion connection events to 100 billion events. We tested Diventi vs. Elasticsearch, a common log indexing tool. In our test environment, Elasticsearch saw an ingestion rate of at best 37,000 events/s while Diventi sustained ingestion rates greater than 171,000 events/s. Our query response times were as much as 100 times faster, typically answering queries in under 80 ms. Furthermore, we saw no noticeable degradation in Diventi from expiration. We have deployed Diventi for many months where it has performed well and supported new security analysis capabilities.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Solid state nuclear magnetic resonance (NMR) spectroscopy and small-to wide-angle X-ray scattering (SWAXS) methods were used to characterize the heterogeneous dynamics and polymer domain structure in rubber modified thermoset materials containing the diglycidyl ether of bisphenol A (DGEBA) epoxy resin and a mixture of Jeffamine reactive rubber and 4,4-diaminodicyclohexylmethane (PACM) amine curing agent. The polymer chain dynamics and morphologies as a function of the PACM/Jeffamine ratio were determined. Using dipolar-filtered NMR experiments, the resulting networks are shown to be composed of mobile and rigid regions that are separated on nanometer length scales, along with a dynamically immobilized interface region. Proton NMR spin diffusion experiments measured the dimensions of the mobile phase to range between 9 and 66 nm and varied with the relative PACM concentration. Solid state 13C magic angle spinning NMR experiments show that the highly mobile phase is composed entirely of the dynamically flexible polyether chains of the Jeffamine rubber, the immobilized interface region is a mixture of DGEBA, PACM, and the Jeffamine rubber, with the PACM cross-linked to DGEBA predominantly residing in the rigid phase. The SWAXS results showed compositional nanophase separation spanning the 11–77 nm range. These measurements of the nanoscale compositional and dynamic heterogeneity provide molecular level insight into the very broad and controllable glass transition temperature distributions observed for these highly cross-linked polymer networks.

Abstract not provided.

Abstract not provided.

The chloride-induced stress corrosion cracking (CISCC) mechanism of cold spray (CS) coating on a galvanically similar substrate is investigated. Arc welded 304L stainless steel (SS) specimens are loaded into four-point bend fixtures, cold sprayed with 304L SS, then immersed in boiling MgCl2. Interconnected porosity forms through crevice corrosion along CS splat boundaries, allowing corrosive species to penetrate through the CS layer. Nevertheless, the substrate is resistant to CISCC likely because of compressive stress introduced by peening during CS particle impacts. These findings underscore the importance of residual stress in the environmental degradation of CS coatings or repairs of engineering structures.

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.

Many engineering design problems can be formulated as decisions between two possible options. This is the case, for example, when a quantity of interest must be maintained below or above some threshold. The threshold thereby determines which input parameters lead to which option, and creates a boundary between the two options known as the decision boundary. This report details a machine learning approach for estimating decision boundaries, based on support vector machines (SVMs), that is amenable to large scale computational simulations. Because it is computationally expensive to evaluate each training sample, the approach iteratively estimates the decision boundary in a manner that requires relatively few training samples to glean useful estimates. The approach is then demonstrated on three example problems from structural mechanics and heat transport.

The Multi-scenario extreme weather simulator (MEWS) is a stochastic weather generation tool. The MEWS algorithm uses 50 or more years of National Oceanic and Atmospheric Association (NOAA) daily summaries [1] for maximum and minimum temperature and NOAA climate norms [2] to calculate historical heat wave and cold snap statistics. The algorithm takes these statistics and shifts them according to multiplication factors provided in the Intergovernmental Panel on Climate Change (IPCC) physical basis technical summary [3] for heat waves.

Abstract not provided.

Abstract not provided.

The COVID-19 pandemic has forced many organizations—from national laboratories to private companies—to change their workforce model to incorporate remote work. This study and the summarized results sought to understand the experiences of remote workers and the ways that remote work can impact recruitment and retention, employee engagement, and career development. Sandia, like many companies, has committed to establishing a hybrid work model that will persist postpandemic, and more Sandia employees than ever before have initiated remote work agreements. This parallels the nationwide increase in remote employment and motivates this study on remote work as an enduring part of workforce models.

Abstract not provided.

Optimization of automated data collection is gaining increased interest for the purposes of enabling closed-loop self-correcting systems that inherently maximize operational efficiencies and reduce waste. Many data collection systems have several variables which influence data accuracy or consistency and which can require frequent user interaction to be monitored and maintained. Operating upon a Robo-MET.3D™ automated mechanical serial-sectioning system, a run-to-run control algorithm has been developed to accelerate data collection and reduce data inconsistency. Using historical data amassed over a decade of experiments, a linear regression model of the deterministic system dynamics is created and used to employ a run-to-run control algorithm that optimizes selected system inputs to reduce operator intervention and increase efficacy while reducing variance of system output.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this work we employ an encoder–decoder convolutional neural network to predict the failure locations of porous metal tension specimens based only on their initial porosities. The process we model is complex, with a progression from initial void nucleation, to saturation, and ultimately failure. The objective of predicting failure locations presents an extreme case of class imbalance since most of the material in the specimens does not fail. In response to this challenge, we develop and demonstrate the effectiveness of data- and loss-based regularization methods. Since there is considerable sensitivity of the failure location to the particular configuration of voids, we also use variational inference to provide uncertainties for the neural network predictions. We connect the deterministic and Bayesian convolutional neural network formulations to explain how variational inference regularizes the training and predictions. We demonstrate that the resulting predicted variances are effective in ranking the locations that are most likely to fail in any given specimen.

Abstract not provided.

Abstract not provided.

Stoichiometric spark-ignition engines suffer efficiency penalties due to throttling losses at low loads, a low specific-heat ratio of the stoichiometric working fluid, and limits on compression ratio due to end-gas autoignition leading to undesirable knocking. Mixed-Mode Combustion (MMC) mitigates these shortcomings by using a lean working fluid where a spark-initiated pilot-stabilized deflagrative flame front is followed by controlled end-gas autoignition. This MMC study investigates the effects of initial conditions (intake air temperature, intake pressure, equivalence ratio, and intake oxygen fraction) on autoignition tendency of four gasoline-range fuels with varying properties and composition. The use of fuels with varying octane sensitivity (S) allowed exploring the importance of low-temperature heat release in triggering autoignition. Fuels with high S were less reactive for conditions that promote low-temperature chemistry (operation at high intake air pressure or without N2 dilution). Conversely, an Alkylate fuel with low S showed a greater autoignition resistance at operating conditions that were unfavorable for low-temperature chemistry. Next, the effect of residual gas composition on autoignition tendency of fuels was examined with a chemical-kinetics model. Among the various molecules in the residual gas, nitric oxide (NO) enhanced the low-temperature chemistry and increased the autoignition tendency most significantly. The fuels’ autoignition response to increasing NO amount corroborates the experimental observations. Next, the sequential autoignition of the end-gas was assessed to be less impacted by thermal stratification because of lean mixtures showing relatively less low-temperature chemistry, when compared to stoichiometric mixtures. Next, the effect of changing equivalence ratio on the autoignition was found to be similar for all fuels, regardless of their S. With changing intake air temperature, the response of fuels’ autoignition tendency depended on the dilution level used. At high dilution (i.e. low intake [O2]), fuels’ reactivity increased with increasing intake air temperature. In contrast, for operation without dilution, the autoignition tendency of the low-S Alkylate fuel decreased with increasing intake air temperature, while that of high-S High Cycloalkane fuel still increased with increasing intake air temperature. In conclusion, conventional octane metrics (RON and MON) have utility in assessing the autoignition tendency under lean MMC operation. Moreover, the fuel requirements for MMC align with that of stoichiometric operation: i.e., high RON and high S fuels are desirable for stable non-knocking operation.

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 Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2022 accomplishments by the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2022, the PFLOTRAN development team made several advancements to our software infrastructure, code performance, and process modeling capabilities.

Abstract not provided.

The Material Protection, Accounting, and Control Technologies (MPACT) program utilizes modeling and simulation to assess Material Control and Accountability (MC&A) concerns for a variety of nuclear facilities. Single analyst tools allow for rapid design and evaluation of advanced approaches for new and existing nuclear facilities. A low enriched uranium (LEU) fuel conversion and fabrication facility simulator has been developed to assist with MC&A for existing LEU fuel fabrication for light water reactors. Simulated measurement blocks were added to the model (consistent with current best practices). Material balance calculations and statistical tests have also been added to the model.

Abstract not provided.

We present an elementary and self-contained proof that there are exactly four symmetry classes of the elasticity tensor in two dimensions: oblique, rectangular, square, and isotropic. In two dimensions, orthogonal transformations are either reflections or rotations. The proof is based on identification of constraints imposed by reflections and rotations on the elasticity tensor, and it simply employs elementary tools from trigonometry, making the proof accessible to a broad audience. For completeness, we identify the sets of transformations (rotations and reflections) for each symmetry class and report the corresponding equations of motions in classical linear elasticity.

Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. A computationally led, material design and synthesis methodology is presented here that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

ERMA is leveraging Sandia’s Microgrid Design Toolkit (MDT) [1] and adding significant new features to it. Development of the MDT was primarily funded by the Department of Energy, Office of Electricity Microgrid Program with some significant support coming from the U.S. Marine Corps. The MDT is a software program that runs on a Microsoft Windows PC. It is an amalgamation of several other software capabilities developed at Sandia and subsequently specialized for the purpose of microgrid design. The software capabilities include the Technology Management Optimization (TMO) application for optimal trade-space exploration, the Microgrid Performance and Reliability Model (PRM) for simulation of microgrid operations, and the Microgrid Sizing Capability (MSC) for preliminary sizing studies of distributed energy resources in a microgrid.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document is aimed at providing guidance to the National Nuclear Security Administration’s (NNSA) Office of International Nuclear Security’s (INS) country and regional teams for implementing effective physical protection systems (PPSs) for nuclear power plants (NPPs) to prevent the radiological consequences of sabotage. This recommendation document includes input from the Physical Protection Functional Team (PPFT), the Response Functional Team (RFT), and the Sabotage Functional Team (SFT) under INS. Specifically, this document provides insights into increasing and sustaining physical protection capabilities at INS partner countries’ NPP sites. Nuclear power plants should consider that the intent of this document is to provide a historical context as well as technologies and methodologies that may be applied to improve physical protection capabilities. It also refers to relevant guidance from the International Atomic Energy Agency (IAEA) and the U.S. Nuclear Regulatory Commission (NRC).

Abstract not provided.

Abstract not provided.

2021-2022 Annual storm water report for LAMP office C930.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Process parameter selection in laser powder bed fusion (LPBF) controls the as-printed dimensional tolerances, pore formation, surface quality and microstructure of printed metallic structures. Measuring the stochastic mechanical performance for a wide range of process parameters is cumbersome both in time and cost. In this study, we overcome these hurdles by using high-throughput tensile (HTT) testing of over 250 dogbone samples to examine process-driven performance of strut-like small features, ~1 mm2 in austenitic stainless steel (316 L). The output mechanical properties, porosity, surface roughness and dimensional accuracy were mapped across the printable range of laser powers and scan speeds using a continuous wave laser LPBF machine. Tradeoffs between ductility and strength are shown across the process space and their implications are discussed. While volumetric energy density deposited onto a substrate to create a melt-pool can be a useful metric for determining bulk properties, it was not found to directly correlate with output small feature performance.

Abstract not provided.

The D-value or dangerous quantity system was designed by the International Commission for Radiological Protection for the determination of source protection categories that can be used to reduce the likelihood of accidents, the consequences of which could result in harm to individuals or costly or expensive cleanup. The process includes multiple scenarios for exposure and two different approaches to the evaluation of detriment. This document provides an example calculation using ¹³⁷Cs to walk through the complex process of determining its D-value in the hopes of making the process easily understandable.

ReNCAT is a software application that suggests microgrid portfolios that reduce the impact of large-scale disruptions to power, as measured by the Social Burden Metric. ReNCAT examines a power distribution network to identify regions that can be isolated into microgrids that enable critical services to be provided even if the remainder of the study area is left without power. ReNCAT operates on a simplified representation of the power grid, one that aggregates and approximates loads and conductors. Microgrids are formed within the power network by setting switch states to split or join portions of the grid. ReNCAT identifies candidate microgrid portfolios with varying tradeoffs between cost and service availability.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.