Publications

Disordered iron germanium (FeGe) has recently garnered interest as a testbed for a variety of magnetic phenomena as well as for use in magnetic memory and logic applications. This is partially owing to its ability to host skyrmions and antiskyrmions—nanoscale whirlpools of magnetic moments that could serve as information carriers in spintronic devices. In particular, a tunable skyrmion-antiskyrmion system may be created through precise control of the defect landscape in B20-phase FeGe, motivating the development of methods to systematically tune disorder in this material and understand the ensuing structural properties. To this end, we investigate a route for modifying magnetic properties in FeGe. In particular, we irradiate epitaxial B20-phase FeGe films with 2.8 MeV Au4+ ions, which creates a dispersion of amorphized regions that may preferentially host antiskyrmions at densities controlled by the irradiation fluence. To further tune the disorder landscape, we conduct a systematic electron diffraction study with in situ annealing, demonstrating the ability to recrystallize controllable fractions of the material at temperatures ranging from ∼150 to 250 °C. Finally, we describe the crystallization kinetics using the Johnson-Mehl-Avrami-Kolmogorov model, finding that the growth of crystalline grains is consistent with diffusion-controlled one-to-two dimensional growth with a decreasing nucleation rate.

The present article summarizes the submissions to the Tribomechadynamics Research Challenge announced in 2021. The task was a blind prediction of the vibration behavior of a system comprising a thin plate clamped on two sides via bolted joints. Both geometric and frictional contact nonlinearities are expected to be relevant. Provided were the CAD models and technical drawings of all parts as well as assembly instructions. The main objective was to predict the frequency and damping ratio of the lowest-frequency mode as function of the amplitude. Many different prediction approaches were pursued, ranging from well-known methods to very recently developed ones. After the submission deadline, the system has been fabricated and tested. The aim of this article is to evaluate the current state of the art in modeling and vibration prediction, and to provide directions for future methodological advancements.

Water security and climate change are important priorities for communities and regions worldwide. The intersections between water and climate change extend across many environmental and human activities. This Primer is intended as an introduction, grounded in examples, for students and others considering the interactions between climate, water, and society. In this Primer, we summarize key intersections between water and climate across four sectors: environment; drinking water, sanitation, and hygiene; food and agriculture; and energy. We begin with an overview of the fundamental water dynamics within each of these four sectors, and then discuss how climate change is impacting water and society within and across these sectors. Emphasizing the relationships and interconnectedness between water and climate change can encourage systems thinking, which can show how activities in one sector may influence activities or outcomes in other sectors. We argue that to achieve a resilient and sustainable water future under climate change, proposed solutions must consider the water–climate nexus to ensure the interconnected roles of water across sectors are not overlooked. Toward that end, we offer an initial set of guiding questions that can be used to inform the development of more holistic climate solutions. This article is categorized under: Science of Water > Water and Environmental Change Engineering Water > Water, Health, and Sanitation Human Water > Value of Water.

Sandia will construct a cold flow receiver test platform in order to characterize falling particle curtain receivers at commercially relevant scales. While Sandia has extensive experience in R&D of falling particle curtain receivers, most have been at pilot scale and smaller - on the order of 1 MW_th with characteristic dimensions of nominally 1-2 m to adequately collect solar energy from the heliostat field at the National Solar Thermal Test Facility (NSTTF). However, scaling up receivers to commercially relevant scales (25 MW_th and above) will require a thorough understanding of particle curtain dynamics at larger scales, especially longer drop heights, for design certainty. The goal of this project will be to construct a cold falling particle curtain test rig capable of simulating particle characteristics that are expected in a commercial scale CSP plant, namely the drop height, curtain thickness, and particle mass flow rate (normalized by length of curtain). This will enable data collection on curtain opacity and spread, both of which are correlated to receiver efficiency and reliable construction, for commercially relevant scales. It will also permit validation of numerical models that will enable detailed receiver characterization and design past currently validated scales.

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

Self-assembled quantum dots (QDs) embedded in InGaAs quantum wells (QWs) are used as active regions for photonic-crystal surface-emitting lasers (PCSELs). An epitaxial regrowth method is developed to fabricate the dot-in-well (DWELL) PCSELs. The epitaxial regrowth starts with the growth of a partial laser structure containing bottom cladding, waveguide, active region, and the photonic crystal (PC) layer. The PC layer is patterned to realize the cavity. Subsequently a top cladding layer is regrown to complete the laser structure. During the regrowth of the top cladding layer, the partial laser structure is subjected to high growth temperatures in excess of 600 °C resulting in an unintentional annealing of the active region. This annealing of the active region can alter the QDs by changing their size resulting in a blue shift in photoluminescence (PL) and narrowing PL emission. This effect results in the misaligning of the gain peak and the cavity resonance, resulting in sub-optimal lasing performance. DWELL active regions are known to have better thermal stability compared to both QDs and QWs and could be an ideal candidate for regrown PCSELs. We successfully demonstrate an optically-pumped epitaxially-regrown DWELL PCSEL with an emission wavelength of 1230 nm operating at room temperature. Furthermore, the DWELL active region shows excellent emission wavelength stability and intensity despite the high temperature regrowth process.

Traditional Monte Carlo methods for particle transport utilize source iteration to express the solution, the flux density, of the transport equation as a Neumann series. Our contribution is to show that the particle paths simulated within source iteration are associated with the adjoint flux density and the adjoint particle paths are associated with the flux density. We make our assertion rigorous through the use of stochastic calculus by representing the particle path used in source iteration as a solution to a stochastic differential equation (SDE). The solution to the adjoint Boltzmann equation is then expressed in terms of the same SDE, and the solution to the Boltzmann equation is expressed in terms of the SDE associated with the adjoint particle process. An important consequence is that the particle paths used within source iteration simultaneously provide Monte Carlo samples of the flux density and adjoint flux density in the detector and source regions, respectively. The significant practical implication is that particle trajectories can be reused to obtain both forward and adjoint quantities of interest. To the best our knowledge, the reuse of entire particles paths has not appeared in the literature. Monte Carlo simulations are presented to support the reuse of the particle paths.

Electrochemical random access memory (ECRAM) is an emerging three-terminal nonvolatile memory (NVM) with highly controllable channel conductance which is promising for use as an analog memory (or synapse) in analog in-memory computing (IMC) systems. Energy-efficient analog IMC computing is particularly desirable for power-constrained, high-radiation environments such as satellites. However, little is known about the suitability of ECRAM for use in a total ionizing dose (TID) environment. This work investigates the effect of Co-60 gamma radiation on the channel conductance and noise—two properties critical for analog IMC systems—of a TaOx-based ECRAM up to 17.3 Mrad(SiO2) for both low- and high-channel-conductance state devices. A transient increase in conductance is observed in response to radiation which consists of two elements: an immediate increase in conductivity due to photocurrent and a secondary increase in conductivity, which has a slower rise and saturation and can persist for hours after exposure. This secondary, persistent photoconductivity is attributed to charging caused by hole trapping. These transient effects would not likely occur in a space environment due to the low dose rate compared with this experiment. No permanent change is found in the low conductance state (LCS) following exposure and the minor shift in the high conductance change would be less significant than the regular retention decay in this state. A permanent increase in the random telegraph noise is observed, possibly due to increased traps created in the channel. This work demonstrates that TaOx-based ECRAM is suitable for use in spaceborne analog IMC systems that are subject to significant TID.

In this paper, we propose a learning-based method utilizing the Soft Actor-Critic (SAC) algorithm to train a binary Support Vector Machine (SVM) classifier. This classifier is designed to identify valid input spaces in high-dimensional, highly constrained systems while minimizing the total runtime of offline simulations. The simulations adapt their runtime based on the likelihood that a given training input will be informative to the classifier. Furthermore, we introduce a method for using the trained SAC model to predict whether a desired system input is likely to violate constraints, along with a technique to adjust the input as necessary. Additionally, we explore the potential of this model to detect faults or adversarial attacks within the system. The effectiveness of our approach is demonstrated through various simulations of challenging classification problems and a constrained quadrotor model.

Appearing as one of the key-components of lithium-ion batteries (LIBs), this work specifically focuses on the additive manufacturing (AM) of custom-shape separators, facilitated by the filament material extrusion process, also called fused deposition modeling (FDM). The development and optimization of composite thermoplastic filament feedstocks combining polypropylene and paraffin wax, followed by the 3D printing of the separator membranes is shown. A post-processing step, based on thermal induced phase separation (TIPS), is introduced to promote porosity formation through removal of the paraffin wax sacrificial phase within the 3D printed items. Separators with different polypropylene/paraffin wax ratios are developed and the impact on printability, mechanical strength, porosity, and electrochemical performances, is thoroughly discussed. X-ray micro-computed tomography is employed to assess the geometric fidelity and to detect printing defects in a complex 3D lattice structure. The performance of the 3D printed porous separators is also compared to a commercial separator. This pioneering research establishes a foundation for the creation of porous separators that can adapt to and conform into 3D printed battery architectures with novel form factors, and also creates opportunities for the use of FDM and TIPS for a wide range of applications that employ porous structures beyond the energy storage field.

Tank farm workers involved in nuclear cleanup activities perform physically demanding tasks, typically while wearing heavy personal protective equipment (PPE). Exoskeleton devices have the potential to bring considerable benefit to this industry but have not been thoroughly studied in the context of nuclear cleanup. In this paper, we examine the performance of exoskeletons during a series of tasks emulating jobs performed on tank farms while participants wore PPE commonly deployed by tank farm workers. The goal of this study was to evaluate the effects of commercially available lower-body exoskeletons on a user’s gait kinematics and user perceptions. Three participants each tested three lower-body exoskeletons in a 70-min protocol consisting of level treadmill walking, incline treadmill walking, weighted treadmill walking, a weight lifting session, and a hand tool dexterity task. Results were compared to a no exoskeleton baseline condition and evaluated as individual case studies. The three participants showed a wide spectrum of user preferences and adaptations toward the devices. Individual case studies revealed that some users quickly adapted to select devices for certain tasks while others remained hesitant to use the devices. Temporal effects on gait change and perception were also observed for select participants in device usage over the course of the device session. Device benefit varied between tasks, but no conclusive aggregate trends were observed across devices for all tasks. Evidence suggests that device benefits observed for specific tasks may have been overshadowed by the wide array of tasks used in the protocol.

This paper is concerned with inserting three-dimensional computer-aided design (CAD) geometries into meshes composed of hexahedral elements using a volume fraction representation. An adaptive procedure for doing so is presented. The procedure consists of two steps. The first step performs spatial acceleration using a k-d tree. The second step involves subdividing individual hexahedra in an adaptive mesh refinement (AMR)-like fashion and approximating the CAD geometry linearly (as a plane) at the finest subdivision. The procedure requires only two geometric queries from a CAD kernel: determining whether or not a queried spatial coordinate is inside or outside the CAD geometry and determining the closest point on the CAD geometry’s surface from a given spatial coordinate. We prove that the procedure is second-order accurate for sufficiently smooth geometries and sufficiently refined background meshes. We demonstrate the expected order of accuracy is achieved with several verification tests and illustrate the procedure’s effectiveness for several exemplar CAD geometries.

The Na+ super ion conductor (NaSICON, Na1+xZr2SixP3-xO12) is a solid electrolyte well-known for fast, selective Na+ transport at low temperatures, uniquely enabling sodium-based batteries. Producing high-quality NaSICON from solid-state methods, especially when cost-effective, potentially hygroscopic precursors are used, is not trivial. To understand and eliminate the influence of humidity during processing, a scheme was developed to reproducibly yield a high Na+ conductivity (3.75 mS/cm at 25°C, 81.7 mS/cm at 150°C), high density (97%), and machinable NaSICON without the use of binders, sintering aids, or dopants. Controlled humidity studies over 20%–50% RH coupled with thermal, structural, and electrical analysis reveal that calcination temperatures < 1000°C leave NaSICON processing susceptible to water absorption at > 20% RH due to the presence of hygroscopic Na3PO4 and Na2CO3 during shaping, pressing, and sintering. Water absorption results in NaSICON with lower densities, machinability, and Na+ conductivity, due to impaired intergranular Na+ transport. At the other extreme, fully converting precursor to the NaSICON phase at 1230°C before pressing and sintering leads to poor conductivity and density. By calcining at 1000°C, excellent quality NaSICON may be produced under a range of laboratory environments, enabling low-cost production of high-conductivity, machinable NaSICON necessary the ever-growing energy storage market.

The Open-Source Offshore (OSO) airfoils have been developed for research purposes for offshore wind turbines, offering a set of airfoils that align with modern turbine design requirements and industry design practices without proprietary constraints on research use. The eventual airfoil family will target the IEA 22 MW reference wind turbine, which was originally developed with the FFA airfoils. The two airfoils summarized in Table 1 (OSO-21-WT1 and OSO-30-WT1) started development as part of a family of airfoils being designed to target the IEA 22 MW wind turbine. The criteria used to design these airfoils are summarized in Table 1, which aim to encapsulate requirements of modern airfoils for offshore wind turbine applications, and were developed with feedback from industry and research experts. The airfoils were designed using XFOIL and candidate airfoils were then analyzed in RFOIL, which is considered more accurate than XFoil for high lift predictions of thicker airfoils. The design process for a preliminary family of airfoils is available, including a more detailed explanation of the design requirements and metrics similar to those used for these airfoils. Most of the design criteria are met for these two airfoils, with two exceptions. For both airfoils, the L/D Roughness Loss metric is exceeded (42% > 40% goal) and the desired lift coefficient margin over the design value (“CL_Margin”) was moderately exceeded (0.43 > 0.3) while smooth-stall characteristics (computed) were achieved. Note that all of the metrics were computed using RFOIL, and like other new airfoils, these will need to be experimentally validated at a range of Reynolds numbers. The airfoil coordinates will be shared publicly on Sadia National Laboratories’ public Github repository:

Abstract not provided.

Code verification is an essential part of credibility analysis for computational models. It assesses whether the mathematical model is implemented correctly into the code and whether the numerical methods behave consistently, and is done before solution verification and validation. Robust guidance for code verification exists in the literature. However, there is no known, concise guide for selecting the approach based on the model form that also presents an overview of the common elements. This document was written to address this gap as an accessible reference for beginning a code-verification effort.

Hardware-in-the-Loop (HIL) methodologies for test system development often require simulations capable of running at MHz speeds on FPGAs. The stringent memory and speed constraints necessitate compromises between model fidelity and execution speed. Numerically solving the underlying governing equations represents the highest accuracy, but slowest responding approach. By storing pre-determined results in look-up-tables (LUT), one can balance speed and accuracy.

In recent years, the procurement of essential materials such as Stainless Steels, Ni-Alloys, and Kovar has become increasingly challenging. Lead times for these materials have extended to over 24 months, placing significant stress on the reliability of attaining these alloys. This situation is exacerbated by the supplier’s minimum purchase requirement of excessively large quantities to meet unique purchase order (PO) requirements which is usually above the maximum requirements of Sandia National Laboratories (SNL). Additionally, the reluctance of foundries to engage in small lot development for SNL specific needs further complicates the procurement process. This white paper explores the reconstitution of SNL’ Melt Laboratory as a strategic solution to mitigate these challenges.

Hydrogen Extremely Low Probability of Rupture (HELPR) is a modular probabilistic fracture mechanics modeling platform developed to assess structural integrity of pipelines for transmission and distribution of hydrogen. HELPR couples fatigue and fracture engineering models with probabilistic methods to generate fast predictions and enables quantification of prediction uncertainty and sensitivity. This user manual serves as a guide through the various analysis features HELPR contains.

During the development of high-consequence items, test systems should be capable of differenti ating between test failures resulting from narrowly missing requirements versus those indicating potentially catastrophic faults. In many instances, classifying the data corresponds to simply identifying whether measured waveforms have approximately the anticipated shape. Cast in this light, the problem reduces to converting raw data into a form optimal for use with neural network classifiers. This manuscript investigates different means of representing raw data for image classification. Raw data plots and Short Time Fourier Transform (STFT) spectrograms are classified by both custom built, small-scale, Convolution Neural Networks (CNN) and open-source, multi-million parameter, pre-trained deep CNNs. In the case of time varying frequency content, the STFTs provide images with greater detail and can be accurately classified with simpler networks. This requires less mem ory and runs faster than classifying the raw data using the more sophisticated options—making STFTs optimal for applications with memory constraints. STFTs are not a panacea. In some cases the time-domain signal contains useful information that should not be discarded. Rather than using raw data or STFTs, the images can be constructed from both by using red and green channels of an RGB image to visualize the real and imaginary components of the transform, with the raw data occupying the blue channel.

Over the past six years, an informal working group has developed to investigate existing sensitivity analysis methods, examine new methods, and identify best practices. The focus is on the use of sensitivity analysis in case studies involving geologic disposal of spent nuclear fuel or nuclear waste. Three additional case studies are presented in this Volume 2 report, including more nonlinear behavior, outputs which exhibit bifurcation, regime changes, and nested sampling.

Continued development of the additive conductivity material model [1], used to simulate changes in heat transfer that occurs in void generating foam decomposition, has resulted in an improved model and new features. The previous version of the model was calibrated against the Aria Bulk Fluid Element (BFE) solution and proposed a third-order polynomial correction term best captured the increased heat transfer due to voids in the foam. An investigation of the Fuego Conjugate Heat Transfer (CHT) and Aria BFE solutions at several geometries revealed the CHT solution and BFE solution had differing behavior across length scales, especially at smaller scales. Five calibration studies, using the Fuego CHT as the calibration data, were carried out with polynomial functions of 4-th, 3-rd, 2-nd, 1-st and 0-th orders to determine the best correction function that generalized well across length scales. Each polynomial function was calibrated/trained on six different sized geometries and then tested on three uniquely sized geometries. This study revealed that the 1-st order additive conductivity model performed the best. A new feature of void formation scaling was implemented to more realistically capture the heat transfer as voids are created. A scaling term was added to the model to activate the conductivity correction as decomposition progresses.

Code verification is an essential part of credibility analysis for computational models. It assesses whether the mathematical model is implemented correctly into the code and whether the numerical methods behave consistently, and is done before solution verification and validation. Robust guidance for code verification exists in the literature. However, there is no known, concise guide for selecting the approach based on the model form that also presents an overview of the common elements. This document was written to address this gap as an accessible reference for beginning a code-verification effort.

Vertical-axis wind turbines (VAWTs) have been the subject of research and development for nearly a century. However, this turbine architecture has fallen in and out of favor on multiple occasions. Beginning in the late 1970s, the U.S. Department of Energy sponsored an extensive experimental program through Sandia National Laboratories which produced a mass of experimental data from several highly instrumented turbines. Turbines designed, built, and tested include the 2 meter, 5 meter, 17 meter, and 34 meter and their respective configurations. This program kicked off a commercial collaboration and resulted in the FloWind turbines. The FloWind turbines had several notable design changes from the experimental turbines that, in conjunction with a general lack of understanding regarding predicting fatigue at the time, led to the majority of the turbines failing prematurely during the late 80s.

This project is intended to support the development of new traction drive systems that meet the targets of 100 kW/L for power electronics and 50 kW/L for electric machines with reliable operation to 300,000 miles. To meet these goals, new designs must be identified that make use of state-of-the-art and next-generation electronic materials and design methods. Designs must exploit synergies between components, for example converters designed for high-frequency switching using wide band gap (WBG) devices and ceramic capacitors. This project included: (1) a survey of available technologies; (2) investigating new technologies, that for example, reduce volume of thermal management or magnetic components; (3) the development of computer aided design tools that consider the converter volume, reliability, and electrical performance; (4) exercising the design software to evaluate performance gaps and predict the impact of certain technologies and design approaches, i.e. GaN semiconductors, ceramic capacitors, ceramic thermal management components, and select topologies; (5) building and testing hardware prototypes to validate models and concepts. The design tools enable co-optimization of the power module and passive elements and provide some design guidance. At the end of the project, new advanced computing methods, such as machine learning approaches, were considered.