Publications

We evaluate the use of reference modules for monitoring effective irradiance in PV power plants, as compared with traditional plane-of-array (POA) irradiance sensors, for PV monitoring and capacity tests. Common POA sensors such as pyranometers and reference cells are unable to capture module-level irradiance nonuniformity and require several correction factors to accurately represent the conditions for fielded modules. These problems are compounded for bifacial systems, where the power loss due to rear side shading and rear-side plane-of-array (RPOA) irradiance gradients are greater and more difficult to quantify. The resulting inaccuracy can have costly real-world consequences, particularly when the data are used to perform power ratings and capacity tests. Here we analyze data from a bifacial single-axis tracking PV power plant, (175.6 MWdc) using 5 meteorological (MET) stations, located on corresponding inverter blocks with capacities over 4 MWdc. Each MET station consists of bifacial reference modules as well pyranometers mounted in traditional POA and RPOA installations across the PV power plant. Short circuit current measurements of the reference modules are converted to effective irradiance with temperature correction and scaling based on flash test or nameplate short circuit values. Our work shows that bifacial effective irradiance measured by pyranometers averages 3.6% higher than the effective irradiance measured by bifacial reference modules, even when accounting for spectral, angle of incidence, and irradiance nonuniformity. We also performed capacity tests using effective irradiance measured by pyranometers and reference modules for each of the 5 bifacial single-axis tracking inverter blocks mentioned above. These capacity tests evaluated bifacial plant performance at ∼3.9% lower when using bifacial effective irradiance from pyranometers as compared to the same calculation performed with reference modules.

Abstract not provided.

Abstract not provided.

At Sandia National Laboratories, QSCOUT (the Quantum Scientific Computing Open User Testbed) is an ion-trap based quantum computer built for the purpose of allowing users low-level access to quantum hardware. Commands are executed on the hardware using Jaqal (Just Another Quantum Assembly Language), a programming language designed in-house to support the unique capabilities of QSCOUT. In this work, we describe a batching implementation of our custom software that speeds the experimental run-time through the reduction of communication and upload times. Reducing the code upload time during experimental runs improves system performance by mitigating the effects of drift. We demonstrate this implementation through a set of quantum chemistry experiments using a variational quantum eigensolver (VQE). While developed specifically for this testbed, this idea finds application across many similar experimental platforms that seek greater hardware control or reduced overhead.

The III-nitride semiconductors are attractive for on-chip, solid-state vacuum nanoelectronics, having high thermal and chemical stability, low electron affinity, and high breakdown fields. Here we report top-down fabricated, lateral gallium nitride (GaN)-based nanoscale vacuum electron diodes operable in air, with ultra-low turn-on voltages down to ~0.24 V, and stable high field emission currents, tested up to several microamps for single-emitter devices. We present gap-size and pressure dependent studies which provide insights into the design of future nanogap vacuum electron devices. The vacuum nanodiodes also show high resistance to damage from 2.5 MeV proton exposure. Preliminary results on the fabrication and characteristics of lateral GaN nano vacuum transistors will also be presented. The results show promise for a new class of robust, integrated, III-nitride based vacuum nanoelectronics.

This paper applies sensitivity and uncertainty analysis to compare two model alternatives for fuel matrix degradation for performance assessment of a generic crystalline repository. The results show that this model choice has little effect on uncertainty in the peak 129I concentration. The small impact of this choice is likely due to the higher importance of uncertainty in the instantaneous release fraction and differences in epistemic uncertainty between the alternatives.

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 violated, the existing system security margins will be largely unknown. As an alternative to the use of complex optimizations over dynamic conditions, this work employs the use of Reinforcement-based 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. The approach requires the training of a Machine Learning Agent to accomplish this task using modeled data and employs it as a decision support tool under severe, near-blackout conditions.

High - temperature particle receivers are being pursued to enable next - generation concentrating solar thermal power (CSP) systems that can achieve higher temperatures (> 700 C) to enable more efficient power cycles, lower overall system costs, and emerging CSP - based process - heat applications. The objective of this work was to develop characterization methods to quantify the particle and heat losses from the open aperture of the particle receiver. Novel camera - based imaging methods were developed and applied to both laboratory - scale and larger 1 MW t on - sun tests at the National Solar Thermal Test Facility in Albuquerque, New Mexico. Validation of the imaging methods was performed using gravimetric and calorimetric methods. In addition, conventional particle - sampling methods using volumetric particle - air samplers were applied to the on - sun tests to compare particle emission rates with regulatory standards for worker safety and pollution. Novel particle sampling methods using 3 - D printed tipping buckets and tethered balloons were also developed and applied to the on - sun particle - receiver tests. Finally, models were developed to simulate the impact of particle size and wind on particle emissions and concentrations as a function of location. Results showed that particle emissions and concentrations were well below regulatory standards for worker safety and pollution. In addition, estimated particle temperatures and advective heat losses from the camera - based imaging methods correlated well with measured values during the on - sun tests.

Type 2 high-pressure hydrogen vessels for storage at hydrogen refueling stations are designed assuming a predefined operational pressure cycle and targeted autofrettage conditions. However, the resulting finite life depends significantly on variables associated with the autofrettage process and the pressure cycles actually realized during service, which many times are not to the full range of the design. Clear guidance for cycle counting is lacking, therefore industry often defaults to counting every repressurization as a full range pressure cycle, which is an overly conservative approach. In-service pressure cycles used to predict the growth of cracks in operational pressure vessels results in significantly longer life, since most in-service pressure cycles are only a fraction of the full design pressure range. Fatigue crack growth rates can vary widely for a given pressure range depending on the details of the residual strains imparted during the autofrettage process because of their influence on crack driving forces. Small changes in variables associated with the autofrettage process, e.g., the target autofrettage overburden pressure, can result in large changes in the residual stress profile leading to possibly degraded fatigue life. In this paper, computational simulation was used for sensitivity studies to evaluate the effect of both operating conditions and autofrettage conditions on fatigue life for Type 2 highpressure hydrogen vessels. The analysis in this paper explores these sensitivities, and the results are used to provide guidance on cycle counting. In particular, we identify the pressure cycle ranges that can be ignored over the life of the vessel as having negligible effect on fatigue life. This study also examines the sensitivity of design life to the autofrettage process and the impact on life if the targeted residual strain is not achieved during manufacturing.

Monitoring cavern leaching after each calendar year of oil sales is necessary to support cavern stability efforts and long-term availability for oil drawdowns in the U.S. Strategic Petroleum Reserve. Modeling results from the SANSMIC code and recent sonars are compared to show projected changes in the cavern’s geometry due to leaching from raw-water injections. This report aims to give background on the importance of monitoring cavern leaching and provide a detailed explanation of the process used to create the leaching plots used to monitor cavern leaching. In the past, generating leaching plots for each cavern in a given leaching year was done manually, and every cavern had to be processed individually. A Python script, compatible with Earth Volumetric Studio, was created to automate most of the process. The script makes a total of 26 plots per cavern to show leaching history, axisymmetric representation of leaching, and SANSMIC modeling of future leaching. The current run time for the script is one hour, replacing 40-50 hours of the monitoring cavern leaching process.

The Jet Propulsion Laboratory has a keen interest in exploring icy moons in the solar system, particularly Jupiter's Europa. Successful exploration of the moon's surface includes planetary protection initiatives to prevent the introduction of viable organisms from Earth to Europa. To that end, the Europa lander requires a Terminal Sterilization Subsystem (TSS) to rid the lander of viable organisms that would potentially contaminate the moon's environment. Sandia National Laboratories has been developing a TSS architecture, relying heavily on computational models to support TSS development. Sandia's TSS design approach involves using energetic material to thermally sterilize lander components at the end of the mission. A hierarchical modeling approach was used for system development and analysis, where simplified systems were constructed to perform empirical tests for evaluating energetic material formulation development and assist in developing computational models with multiple tiers of physics fidelity. Computational models have been developed using multiple Sandia-native computational tools. Three experimental systems and corresponding computational models have been developed: Tube, Sub-Box Small, and Sub-Box Large systems. This paper presents an explanation of the application context of the TSS along with an overview description of a small portion of the TSS development from a modeling and simulation perspective, specifically highlighting verification, validation, and uncertainty quantification (VVUQ) aspects of the modeling and simulation work. Multiple VVUQ approaches were implemented during TSS development, including solution verification, calibration, uncertainty quantification, global sensitivity analysis, and validation. This paper is not intended to express the design results or parameter values used to model the TSS but to communicate the approaches used and how the results of the VVUQ efforts were used and interpreted to assist system development.

Abstract not provided.

The development of the High-Resolution Wavelet Transform (HRWT) is driven by the need of increasing the high-frequency resolution of widely used discrete Wavelet Transforms (WTs). Based on the Stationary Wavelet Transform (SWT), which is a modification of the Discrete Wavelet Transform (DWT), a novel WT that increases the number of decomposition levels (therefore increasing the previously mentioned frequency resolution) is proposed. In order to show the validity of the HRWT, this paper encompasses a theoretical comparison with other discrete WT methods. First, a summary of the DWT and the SWT, along with a brief explanation of the WT theory, is provided. Then, the concept of the HRWT is presented, followed by a discussion of the adherence of this new method to the WT's common properties. Finally, an example of the application is performed on a transient waveform analysis from a power system fault event, outlining the benefits that can be obtained from its usage compared to the SWT.

Wind turbine wakes are characterized by helical trailing tip vortices that are highly stable initially and act as a shield against mixing with the ambient flow and thereby delay wake recovery until destructive mutual interference of the vortices begins. Delayed wake recovery in turn reduces the power production of downstream turbines that are positioned in the wakes of upstream turbines. The long natural decay length forces wind farms to have large distances between turbines to yield sufficient wake recovery. Herein, we tested a new concept aimed at accelerating the breakdown of wind turbine tip vortices by causing the vortices to interact with one another almost immediately behind the rotor. By adding a spire behind the rotor, essentially a blockage to perturb the paths of the tip vortices, we hypothesized that the altered paths of the tip vortices would cause their destructive interference process to begin sooner. The concept of a nacelle-mounted spire was tested in high-fidelity large-eddy simulations using Nalu-Wind. Four different spires were modeled with wall-resolved meshes behind the rotor of a wind turbine with another turbine five diameters downstream. We compared power and wake data against baseline results to determine whether the spires accelerated wake recovery of the upstream turbine and thereby increased the power of the downstream turbine. The results showed no change in the total power of the two turbines for any spire compared to its respective baseline. These results were further explored by testing at higher spatial resolution and without turbulence in the inflow. The increased spatial resolution increased the apparent stability of the tip vortices while the lack of turbulence did not. We conclude that the spires’ geometry and size were inadequate to alter the helical paths of the trailing tip vortices and that modeling of the formation and decay of tip vortices may be highly sensitive to model parameters.

Pre-chamber ignition has demonstrated capability to increase internal combustion engine in-cylinder burn rates and enable the use of low engine-out pollutant emission combustion strategies. In the present study, newly designed passive pre-chambers with different nozzle-hole patterns - that featured combinations of radial and axial nozzles - were experimentally investigated in an optically accessible, single-cylinder research engine. The pre-chambers analyzed had a narrow throat geometry to increase the velocity of the ejected jets. In addition to a conventional inductive spark igniter, a nanosecond spark ignition system that promotes faster early burn rates was also investigated. Time-resolved visualization of ignition and combustion processes was accomplished through high-speed hydroxyl radical (OH*) chemiluminescence imaging. Pressure was measured during the engine cycle in both the main chamber and pre-chamber to monitor respective combustion progress. Experimental heat release rates (HRR) calculated from the measured pressure profiles were used as inputs for two different GT-Power 1D simulations to evaluate the pre-chamber jet-exit momentum and penetration distance. The first simulation used both the calculated main-chamber and pre-chamber HRR, while the second used only the main chamber HRR with the pre-chamber HRR modeled. Results show discrepancies between the models mainly in the pressurization of the pre-chamber which in turn affected jet penetration rate and highlights the sensitivity of the simulation results to proper input selection. Experimental results further show increased pressurization, with an associated acceleration of jet penetration, when operating with nanosecond spark ignition systems regardless of the pre-chamber tip geometry used.

This paper presents a visualization technique for incorporating eigenvector estimates with geospatial data to create inter-area mode shape maps. For each point of measurement, the method specifies the radius, color, and angular orientation of a circular map marker. These characteristics are determined by the elements of the right eigenvector corresponding to the mode of interest. The markers are then overlaid on a map of the system to create a physically intuitive visualization of the mode shape. This technique serves as a valuable tool for differentiating oscillatory modes that have similar frequencies but different shapes. This work was conducted within the Western Interconnection Modes Review Group (WIMRG) in the Western Electric Coordinating Council (WECC). For testing, we employ the WECC 2021 Heavy Summer base case, which features a high-fidelity, industry standard dynamic model of the North American Western Interconnection. Mode estimates are produced via eigen-decomposition of a reduced-order state matrix identified from simulated ringdown data. The results provide improved physical intuition about the spatial characteristics of the inter-area modes. In addition to offline applications, this visualization technique could also enhance situational awareness for system operators when paired with online mode shape estimates.

A high-speed, two-color pyrometer was developed and employed to characterize the temperature of the ejecta from pyrotechnic igniters. The pyrometer used a single objective lens, beamsplitter, and two high-speed cameras to maximize the spatial and temporal resolutions. The pyrometer used the integrated intensity of under-resolved particles to maintain a large region of interest to capture more particles. The spectral response of the pyrometer was determined based on the response of each optical component and the total system was calibrated using a black body source to ensure accurate intensity ratios over the range of interest.

Abstract not provided.

Soft error rates (SER) are characterized for the 5-nm bulk FinFET D flip-flops for alpha particles, thermal neutrons, and high-energy neutrons as a function of supply voltage. At nominal operating voltage, the 5-nm node has higher SER than the 7-nm node for all three particle types, with increases of 148%, 168%, and 26%, respectively. The overall SER for the 5-nm node was ~2X greater than that of the 7-nm node, because the reduction in critical charge was higher than that in collected charge. For alpha particle exposures, temperature effects on SER were more prominent for the 5-nm node than both the 7-nm and 16-nm node. Relative contribution of alpha particle SER increases with scaling, and it accounts for 13% of the overall SER at the 5-nm node.

For the protection engineer, it is often the case, that full coverage and thus perfect selectivity of the system is not an option for protection devices. This is because perfect selectivity requires protection devices on every line section of the network. Due to cost limitation, relays may not be placed on each branch of a network. Therefore, a method is needed to allow for optimal coordination of relays with sparse relay placement. In this paper, methods for optimal coordination of networks with sparse relay placement introduced in prior work are applied to a system where both overcurrent and distance relays are present. Additionally, a method for defining primary (Zone 1) and secondary (Zone 2) protection zones for the distance relays in such a sparse system is proposed. The proposed method is applied to the IEEE 123-bus test case. The proposed method is found to successfully coordinate the system while also limiting the maximum relay operating time to 1.78s which approaches the theoretical lower bound of 1.75s.

Analog in-memory computing is a method to improve the efficiency of deep neural network inference by orders of magnitude, by utilizing analog properties of a nonvolatile memory. This places new requirements on the memory device, which physically represent neural net weights as analog states. By carefully considering the algorithm implications when mapping weights to physical states, it is possible to achieve precision very close to that of a digital accelerator using a 40nm embedded SONOS.

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged Navier–Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal model validation with uncertainty considerations that leverages experimental data from the HIFiRE-1 wind tunnel tests. The geometry is a multi-conic shape that produces complex flow phenomena under hypersonic conditions. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

Snow and ice accumulation on photovoltaic (PV) panels is a recognized-but poorly quantified-contributor to PV performance, not only in geographic areas that see persistent snow in winter but also at lower latitudes, where frozen precipitation and 'snowmageddon' events can wreak havoc with the solar infrastructure. In addition, research on the impact of snow and cold on PV systems has not kept pace with the proliferation of new technologies, the rapid deployment of PV in northern latitudes, and experiences with long-term field performance. This paper describes the value of a dedicated outdoor research facility for longitudinal performance and reliability studies of emerging technologies in cold climates.

In this work we introduce Bootstrapped Paired Neural Networks (BPNN), a semi-supervised, low-shot model with uncertainty quantification (UQ). BPNN can be used for classification and target detection problems commonly encountered when working with aerospace imagery data, such as hyperspectral imagery (HSI) data. When collecting aerospace imaging data, there is often large amounts of data which can be costly to label, so we would like to supplement labeled data with the vast unlabeled data (often > 90% of the data) available, we do this using semi-supervised techniques (Exponential Average Adversarial Training). Often, it is difficult and costly to obtain the sample size necessary to train a deep learning model to a new class or target; using paired neural networks (PNN), our model is generalized to low-and no-shot learning by learning an embedding space for which the underlying data population lives, this way additional labeled data may not be necessary to detect for targets or classes which weren't originally trained on. Finally, by bootstrapping the PNN, the BPNN model gives an uncertainty score on predicted classifications with minimal statistical distributional assumptions. Uncertainty is necessary in the high consequence problems that many applications in aerospace endure. The model's ability to provide uncertainty for its own predictions can be used to reduce false alarms rates, provide explainability to black box models, and help design efficient future data collection campaigns. Although models exist to contain two of these three qualities, to our knowledge no model contains all three: semi-supervised, low-shot, and uncertainty quantification. We generate a HSI scene using a high fidelity data simulator that gives us ground truth radiance spectra, allowing us to fully assess the quality of our model and compare to other common models. When applying PBNN to our HSI scene, it outperforms in target detection against classic methods, such as Adaptive Cosine Estimator (ACE), simple deep learning models without low-shot or semi-supervised capabilities, and models using only low-shot learning techniques such as regular PNN. When extending to targets not originally trained on, the model again outperforms the regular PNN. Using the UQ of predictions, we create 'high confidence sets' which contain predictions which are reliably correct and can help suppress false alarms. This is shown by the higher performance of the 'high confidence set' at particular constant false alarm rates. They also provide an avenue for automation while other predictions in high consequence situations might need to be analyzed further. BPNN is a powerful new predictive model that could be used to maximize the data collected by aerial assets while instilling confidence in model predictions for high consequence situations and being flexible enough to find previously unobserved targets.

The application of traveling wave principles for fault detection in distribution systems is challenging because of multiple reflections from the laterals and other lumped elements, particularly when we consider communication-free applications. We propose and explore the use of Shapelets to characterize fault signatures and a data-driven machine learning model to accurately classify the faults based on their distance. Studies of a simple 5-bus system suggest that the use of Shapelets for detecting faults is promising. The application to practical three-phase distribution feeders is the subject of continuing research.

Abstract not provided.

Calibrating a finite element model to test data is often required to accurately characterize a joint, predict its dynamic behavior, and determine fastener fatigue life. In this work, modal testing, model calibration, and fatigue analysis are performed for a bolted structure, and various joint modeling techniques are compared. The structure is designed to test a single bolt to fatigue failure by utilizing an electrodynamic modal shaker to axially force the bolted joint at resonance. Modal testing is done to obtain the dynamic properties, evaluate finite element joint modeling techniques, and assess the effectiveness of a vibration approach to fatigue testing of bolts. Results show that common joint models can be inaccurate in predicting bolt loads, and even when updated using modal test data, linear structural models alone may be insufficient in evaluating fastener fatigue.

Incorrect modeling of control characteristics for inverter-based resources (IBRs) can affect the accuracy of electric power system studies. In many distribution system contexts, the control settings for behind-the-meter (BTM) IBRs are unknown. This paper presents an efficient method for selecting a small number of time series samples from net load meter data that can be used for reconstructing or classifying the control settings of BTM IBRs. Sparse approximation techniques are used to select the time series samples that cause the inversion of a matrix of candidate responses to be as well-conditioned as possible. We verify these methods on 451 actual advanced metering infrastructure (AMI) datasets from loads with BTM IBRs. Selecting 60 15-minute granularity time series samples, we recover BTM control characteristics with a mean error less than 0.2 kVAR.

Propagating thermal runaway events are a significant threat to utility-scale storage installations. A propagating thermal runaway event is a cascading series of failures in which energy released from a failed cell triggers subsequent failures in nearby cells. Without intervention, propagation can turn an otherwise manageable single cell failure into a full system conflagration. This study presents a method of mitigating the severity of propagating thermal runaway events in utility-scale storage systems by leveraging the capabilities of a module-interfaced power conversion architecture. The method involves strategic depletion of storage modules to delay or arrest propagation, reducing the total thermal energy released in the failure event. The feasibility of the method is assessed through simulations of propagating thermal runaway events in a 160 kW/80 kWh energy storage system.

Abstract not provided.

Variable energy resources (VERs) like wind and solar are the future of electricity generation as we gradually phase out fossil fuel due to environmental concerns. Nations across the globe are also making significant strides in integrating VERs into their power grids as we strive toward a greener future. However, integration of VERs leads to several challenges due to their variable nature and low inertia characteristics. In this paper, we discuss the hurdles faced by the power grid due to high penetration of wind power generation and how energy storage system (ESSs) can be used at the grid-level to overcome these hurdles. We propose a new planning strategy using which ESSs can be sized appropriately to provide inertial support as well as aid in variability mitigation, thus minimizing load curtailment. A probabilistic framework is developed for this purpose, which takes into consideration the outage of generators and the replacement of conventional units with wind farms. Wind speed is modeled using an autoregressive moving average technique. The efficacy of the proposed methodology is demonstrated on the WSCC 9-bus test system.

Abstract not provided.

Abstract not provided.

This work investigates both avalanche behavior and failure mechanism of 3 kV GaN-on-GaN vertical P-N diodes, that were fabricated and later tested under unclamped inductive switching (UIS) stress. The goal of this study is to use the particular avalanche characteristics and the failure mechanism to identify issues with the field termination and then provide feedback to improve the device design. DC breakdown is measured at the different temperatures to confirm the avalanche breakdown. Diode's avalanche robustness is measured on-wafer using a UIS test set-up which was integrated with a wafer chuck and CCD camera. Post failure analysis of the diode is done using SEM and optical microscopy to gain insight into the device failure physics.

Injecting CO2 into a deep geological formation (i.e., geological carbon storage, GCS) can induce earthquakes along preexisting faults in the earth's upper crust. Seismic survey and regional geo-structure analysis are typically employed to map the faults prone to earthquakes prior to injection. However, earthquakes induced by fluid injection from other subsurface energy storage and recovery activities show that systematic evaluation of the potential of induced seismicity associated with GCS is necessary. This study mechanistically investigates how multiphysical interaction among injected CO2, preexisting pore fluids and rock matrix alters stress states on faults and which physical mechanisms can nucleate earthquakes along the faults. Increased injection pressure is needed to overcome capillary entry pressure of the fault zone, driven by the contrast of fluids' wetting characteristics. Accumulated CO2 within the reservoir delays post shut-in reduction in pressure and stress fields along the fault that may enhance the potential for earthquake nucleation after terminating injection operations. Elastic energy generated by coupled processes transfers to low-permeability or hydraulically isolated basement faults, which can initiate slip of the faults. Our findings from generic studies suggest that geomechanical simulations integrated with multiphase flow system are essential to detect deformation-driven signals and mitigate potential seismic hazards associated with CO2 injection.

We investigate the space complexity of two graph streaming problems: MAX-CUT and its quantum analogue, QUANTUM MAX-CUT. Previous work by Kapralov and Krachun [STOC 19] resolved the classical complexity of the classical problem, showing that any (2 - ?)-approximation requires O(n) space (a 2-approximation is trivial with O(log n) space). We generalize both of these qualifiers, demonstrating O(n) space lower bounds for (2 - ?)-approximating MAX-CUT and QUANTUM MAX-CUT, even if the algorithm is allowed to maintain a quantum state. As the trivial approximation algorithm for QUANTUM MAX-CUT only gives a 4-approximation, we show tightness with an algorithm that returns a (2 + ?)-approximation to the QUANTUM MAX-CUT value of a graph in O(log n) space. Our work resolves the quantum and classical approximability of quantum and classical Max-Cut using o(n) space.We prove our lower bounds through the techniques of Boolean Fourier analysis. We give the first application of these methods to sequential one-way quantum communication, in which each player receives a quantum message from the previous player, and can then perform arbitrary quantum operations on it before sending it to the next. To this end, we show how Fourier-analytic techniques may be used to understand the application of a quantum channel.

Structural alloys may experience corrosion when exposed to molten chloride salts due to selective dissolution of active alloying elements. One way to prevent this is to make the molten salt reducing. For the KCl + MgCl2 eutectic salt mixture, pure Mg can be added to achieve this. However, Mg can form intermetallic compounds with nickel at high temperatures, which may cause alloy embrittlement. This study shows that an optimum level of excess Mg could be added to the molten salt which will prevent corrosion of alloys like 316 H, while not forming any detectable Ni-Mg intermetallic phases on Ni-rich alloy surfaces.

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.4 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

Abstract not provided.

Neural networks (NN)s have been increasingly proposed as surrogates for approximation of systems with computationally expensive physics for rapid online evaluation or exploration. As these surrogate models are integrated into larger optimization problems used for decision making, there is a need to verify their behavior to ensure adequate performance over the desired parameter space. We extend the ideas of optimization-based neural network verification to provide guarantees of surrogate performance over the feasible optimization space. In doing so, we present formulations to represent neural networks within decision-making problems, and we develop verification approaches that use model constraints to provide increasingly tight error estimates. We demonstrate the capabilities on a simple steady-state reactor design problem.

Abstract not provided.

For many industries addressing varied customer needs means producing a family of products that satisfy a range of design requirements. Manufacturers seek to design this family of products while exploiting opportunities for shared components to reduce manufacturing cost and complexity. We present a mixed-integer programming formulation that determines the optimal design for each product, the number and design of shared components, and the allocation of those shared components across the products in the family. This formulation and workflow for product family design has created significant business impact on the industrial design of product families for large-scale commercial HVAC chillers in Carrier Global Corporation. We demonstrate the approach on an open case study based on a transcritical CO2 refrigeration cycle. This case study and our industrial experience show that the formulation is computationally tractable and can significantly reduce engineering time by replacing the manual design process with an automated approach.

Quantifying gas-surface interactions for hypersonic reentry applications remains a challenging and complex problem where credible models are needed to design and analyze thermal protection systems. A flexible sensitivity analysis approach is demonstrated to analyze finite-rate ablation models to identify reaction parameters and mechanisms of influence on predicted quantities of interest. Simulations of hypersonic flow over a sphere-cone are presented using parameterized Park, Zhluktov and Abe (ZA), and MURI finite-rate models that describe the oxidation and sublimation of carbon. The results presented in this study emphasize the importance of characterizing model inputs that are shown to have a high impact on predicted quantities and build evidence to assess credibility of these models.

Abstract not provided.

In high temperature (HT) environments often encountered in geothermal wells, data rate transfers for downhole instrumentation are relatively limited due to transmission line bandwidth and insertion loss and the processing speed of HT microcontrollers. In previous research, Sandia National Laboratory Geothermal Department obtained 3.8 Mbps data rates over 1524 m (5000 ft) for single conductor wireline cable with less than a 1x10-8 bit error rate utilizing low temperature NITM hardware (formerly National InstrumentsTM). Our protocol technique was a combination of orthogonal frequency-division multiplexing and quadrature amplitude modulation across the bandwidth of the single conductor wireline. This showed it is possible to obtain high data rates in low bandwidth wirelines. This paper focuses on commercial HT microcontrollers (µC), rather than low temperature NITM modules, to enable high-speed communication in an HT environment. As part of this effort, four devices were evaluated, and an optimal device (SM320F28335-HT) was selected for its high clock rates, floating-point unit, and on-board analog-to-digital converter. A printed circuit board was assembled with the HT µC, an HT resistor digital-to-analog converter, and an HT line driver. The board was tested at the microcontroller's rated maximum temperature (210°C) for a week while transmitting through a 1524 m (5000 ft) wireline. A final test was conducted to the point of failure at elevated temperatures. This paper will discuss communication methods, achieved data rates, and hardware selection. This effort contributes to the enhancement of HT instrumentation by enabling greater sensor counts and improving data accuracy and transfer rates.

The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. The test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.

Calibrating a finite element model to test data is often required to accurately characterize a joint, predict its dynamic behavior, and determine fastener fatigue life. In this work, modal testing, model calibration, and fatigue analysis are performed for a bolted structure, and various joint modeling techniques are compared. The structure is designed to test a single bolt to fatigue failure by utilizing an electrodynamic modal shaker to axially force the bolted joint at resonance. Modal testing is done to obtain the dynamic properties, evaluate finite element joint modeling techniques, and assess the effectiveness of a vibration approach to fatigue testing of bolts. Results show that common joint models can be inaccurate in predicting bolt loads, and even when updated using modal test data, linear structural models alone may be insufficient in evaluating fastener fatigue.

In accident scenarios involving release of tritium during handling and storage, the level of risk to human health is dominated by the extent to which radioactive tritium is oxidized to the water form (T2O or THO). At some facilities, tritium inventories consist of very small quantities stored at sub-atmospheric pressure, which means that tritium release accident scenarios will likely produce concentrations in air that are well below the lower flammability limit. It is known that isotope effects on reaction rates should result in slower oxidation rates for heavier isotopes of hydrogen, but this effect has not previously been quantified for oxidation at concentrations well below the lower flammability limit for hydrogen. This work describes hydrogen isotope oxidation measurements in an atmospheric tube furnace reactor. These measurements consist of five concentration levels between 0.01% and 1% protium or deuterium and two residence times. Oxidation is observed to occur between about 550°C and 800°C, with higher levels of conversion achieved at lower temperatures for protium with respect to deuterium at the same volumetric inlet concentration and residence time. Computational fluid dynamics simulations of the experiments were used to customize reaction orders and Arrhenius parameters in a 1-step oxidation mechanism. The trends in the rates for protium and deuterium are extrapolated based on guidance from literature to produce kinetic rate parameters appropriate for tritium oxidation at low concentrations.

Variables estimated by Battery Management Systems (BMSs) such as the State of Charge (SoC) may be vulnerable to False Data Injection Attacks (FDIAs). Bad actors could use FDIAs to manipulate sensor readings, which could degrade Battery Energy Storage Systems (BESSs) or result in poor system performance. In this paper we propose a method for accurate SoC estimation for series-connected stacks of batteries and detection of FDIA in cell and stack voltage sensors using physics-based models, an Extended Kalman Filter (EKF), and a Cumulative Sum (CUSUM) algorithm. Utilizing additional sensors in the battery stack allowed the system to remain observable in the event of a single sensor failure, allowing the system to continue to accurately estimate states when one sensor at a time was offline. A priori residual data for each voltage sensor was used in the CUSUM algorithm to find the minimum detectable attack (500 μV) with no false positives.

Wind turbine wakes are characterized by helical trailing tip vortices that are highly stable initially and act as a shield against mixing with the ambient flow and thereby delay wake recovery until destructive mutual interference of the vortices begins. Delayed wake recovery in turn reduces the power production of downstream turbines that are positioned in the wakes of upstream turbines. The long natural decay length forces wind farms to have large distances between turbines to yield sufficient wake recovery. Herein, we tested a new concept aimed at accelerating the breakdown of wind turbine tip vortices by causing the vortices to interact with one another almost immediately behind the rotor. By adding a spire behind the rotor, essentially a blockage to perturb the paths of the tip vortices, we hypothesized that the altered paths of the tip vortices would cause their destructive interference process to begin sooner. The concept of a nacelle-mounted spire was tested in high-fidelity large-eddy simulations using Nalu-Wind. Four different spires were modeled with wall-resolved meshes behind the rotor of a wind turbine with another turbine five diameters downstream. We compared power and wake data against baseline results to determine whether the spires accelerated wake recovery of the upstream turbine and thereby increased the power of the downstream turbine. The results showed no change in the total power of the two turbines for any spire compared to its respective baseline. These results were further explored by testing at higher spatial resolution and without turbulence in the inflow. The increased spatial resolution increased the apparent stability of the tip vortices while the lack of turbulence did not. We conclude that the spires’ geometry and size were inadequate to alter the helical paths of the trailing tip vortices and that modeling of the formation and decay of tip vortices may be highly sensitive to model parameters.

Identifying the location of faults in a fast and accurate manner is critical for effective protection and restoration of distribution networks. This paper describes an efficient method for detecting, localizing, and classifying faults using advanced signal processing and machine learning tools. The method uses an Isolation Forest technique to detect the fault. Then Continuous Wavelet Transform (CWT) is used to analyze the traveling waves produced by the faults. The CWT coefficients of the current signals at the time of arrival of the traveling wave present unique characteristics for different fault types and locations. These CWT coefficients are fed into a Convolutional Neural Network (CNN) to train and classify fault events. The results show that for multiple fault scenarios and solar PV conditions, the method is able to determine the fault type and location with high accuracy.

Abstract not provided.

Turbine generator power from simulations using Actuator Line Models and Actuator Disk Models with a Filtered Lifting Line Correction are compared to field data of a V27 turbine. Preliminary results of the wake characteristics are also presented. Turbine quantities of interest from traditional ALM and ADM with the Gaussian kernel (ϵ) set at the optimum value for matching power production and that resolve the kernel at all mesh sizes are also presented. The atmospheric boundary layer is simulated using Nalu-Wind, a Large Eddy Simulation code which is part of the ExaWind code suite. The effect of mesh resolution on quantities of interest is also examined.

Metasurface lenses are fabricated using membrane projection lithography following a CMOS-compatible process flow. The lenses are 10-mm in diameter and employ 3-dimensional unit cells designed to function in the mid-infrared spectral range.

Downtown low-voltage (LV) distribution networks are generally protected with network protectors that detect faults by restricting reverse power flow out of the network. This creates protection challenges for protecting the system as new smart grid technologies and distributed generation are installed. This report summarizes well-established methods for the control and protection of LV secondary network systems and spot networks, including operating features of network relays. Some current challenges and findings are presented from interviews with three utilities, PHI PEPCO, Oncor Energy Delivery, and Consolidated Edison Company of New York. Opportunities for technical exploration are presented with an assessment of the importance or value and the difficulty or cost. Finally, this leads to some recommendations for research to improve protection in secondary networks.

Abstract not provided.

We present a procedure for randomly generating realistic steady-state contingency scenarios based on the historical outage data from a particular event. First, we divide generation into classes and fit a probability distribution of outage magnitude for each class. Second, we provide a method for randomly synthesizing generator resilience levels in a way that preserves the data-driven probability distributions of outage magnitude. Finally, we devise a simple method of scaling the storm effects based on a single global parameter. We apply our methods using data from historical Winter Storm Uri to simulate contingency events for the ACTIVSg2000 synthetic grid on the footprint of Texas.

Abstract not provided.

This paper presents the formulation, implementation, and demonstration of a new, largely phenomenological, model for the damage-free (micro-crack-free) thermomechanical behavior of rock salt. Unlike most salt constitutive models, the new model includes both drag stress (isotropic) and back stress (kinematic) hardening. The implementation utilizes a semi-implicit scheme and a fall-back fully-implicit scheme to numerically integrate the model's differential equations. Particular attention was paid to the initial guesses for the fully-implicit scheme. Of the four guesses investigated, an initial guess that interpolated between the previous converged state and the fully saturated hardening state had the best performance. The numerical implementation was then used in simulations that highlighted the difference between drag stress hardening versus combined drag and back stress hardening. Simulations of multi-stage constant stress tests showed that only combined hardening could qualitatively represent reverse (inverse transient) creep, as well as the large transient strains experimentally observed upon switching from axisymmetric compression to axisymmetric extension. Simulations of a gas storage cavern subjected to high and low gas pressure cycles showed that combined hardening led to substantially greater volume loss over time than drag stress hardening alone.

Abstract not provided.

Neural networks (NN) have become almost ubiquitous with image classification, but in their standard form produce point estimates, with no measure of confidence. Bayesian neural networks (BNN) provide uncertainty quantification (UQ) for NN predictions and estimates through the posterior distribution. As NN are applied in more high-consequence applications, UQ is becoming a requirement. Automating systems can save time and money, but only if the operator can trust what the system outputs. BNN provide a solution to this problem by not only giving accurate predictions and estimates, but also an interval that includes reasonable values within a desired probability. Despite their positive attributes, BNN are notoriously difficult and time consuming to train. Traditional Bayesian methods use Markov Chain Monte Carlo (MCMC), but this is often brushed aside as being too slow. The most common method is variational inference (VI) due to its fast computation, but there are multiple concerns with its efficacy. MCMC is the gold standard and given enough time, will produce the correct result. VI, alternatively, is an approximation that converges asymptotically. Unfortunately (or fortunately), high consequence problems often do not live in the land of asymtopia so solutions like MCMC are preferable to approximations. We apply and compare MCMC-and VI-trained BNN in the context of target detection in hyperspectral imagery (HSI), where materials of interest can be identified by their unique spectral signature. This is a challenging field, due to the numerous permuting effects practical collection of HSI has on measured spectra. Both models are trained using out-of-the-box tools on a high fidelity HSI target detection scene. Both MCMC-and VI-trained BNN perform well overall at target detection on a simulated HSI scene. Splitting the test set predictions into two classes, high confidence and low confidence predictions, presents a path to automation. For the MCMC-trained BNN, the high confidence predictions have a 0.95 probability of detection with a false alarm rate of 0.05 when considering pixels with target abundance of 0.2. VI-trained BNN have a 0.25 probability of detection for the same, but its performance on high confidence sets matched MCMC for abundances >0.4. However, the VI-trained BNN on this scene required significant expert tuning to get these results while MCMC worked immediately. On neither scene was MCMC prohibitively time consuming, as is often assumed, but the networks we used were relatively small. This paper provides an example of how to utilize the benefits of UQ, but also to increase awareness that different training methods can give different results for the same model. If sufficient computational resources are available, the best approach rather than the fastest or most efficient should be used, especially for high consequence problems.

Remote assessment of physiological parameters has enabled patient diagnostics without the need for a medical professional to become exposed to potential communicable diseases. In particular, early detection of oxygen saturation, abnormal body temperature, heart rate, and/or blood pressure could affect treatment protocols. The modeling effort in this work uses an adding-doubling radiative transfer model of a seven-layer human skin structure to describe absorption and reflection of incident light within each layer. The model was validated using both abiotic and biotic systems to understand light interactions associated with surfaces consisting of complex topography as well as multiple illumination sources. Using literature-based property values for human skin thickness, absorption, and scattering, an average deviation of 7.7% between model prediction and experimental reflectivity was observed in the wavelength range of 500-1000 nm.

Here, we describe the isolation of 18 unique anti SARS-CoV-2 human single-chain antibodies from an antibody library derived from healthy donors. The selection used a combination of phage and yeast display technologies and included counter-selection strategies meant to direct the selection of the receptor-binding motif (RBM) of SARS-CoV-2 spike protein’s receptor binding domain (RBD2). Selected antibodies were characterized in various formats including IgG, using flow cytometry, ELISA, high throughput SPR, and fluorescence microscopy. We report antibodies’ RBD2 recognition specificity, binding affinity, and epitope diversity, as well as ability to block RBD2 binding to the human receptor angiotensin-converting enzyme 2 (ACE2) and to neutralize authentic SARS-CoV-2 virus infection in vitro. We present evidence supporting that: 1) most of our antibodies (16 out of 18) selectively recognize RBD2; 2) the best performing 8 antibodies target eight different epitopes of RBD2; 3) one of the pairs tested in sandwich assays detects RBD2 with sub-picomolar sensitivity; and 4) two antibody pairs inhibit SARS-CoV-2 infection at low nanomolar half neutralization titers. Based on these results, we conclude that our antibodies have high potential for therapeutic and diagnostic applications. Importantly, our results indicate that readily available non immune (naïve) antibody libraries obtained from healthy donors can be used to select high-quality monoclonal antibodies, bypassing the need for blood of infected patients, and offering a widely accessible and low-cost alternative to more sophisticated and expensive antibody selection approaches (e.g. single B cell analysis and natural evolution in humanized mice).

Abstract not provided.

Sedimentary-hosted geothermal energy systems are permeable structural, structural-stratigraphic, and/or stratigraphic horizons with sufficient temperature for direct use and/or electricity generation. Sedimentary-hosted (i.e., stratigraphic) geothermal reservoirs may be present in multiple locations across the central and eastern Great Basin of the USA, thereby constituting a potentially large base of untapped, economically accessible energy resources. Sandia National Laboratories has partnered with a multi-disciplinary group of collaborators to evaluate a stratigraphic system in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The goal of this study is to inform an optimized strategy for subsequent exploration and development of this and analogous resources. Building from prior Nevada Play Fairway Analysis (PFA), this team is primarily 1) collecting additional geophysical data, 2) employing novel joint geophysical inversion/modeling techniques to update existing 3D geologic models, and 3) integrating the geophysical results to produce a working, geologically constrained thermo-hydrological reservoir model. Prior PFA work highlights Steptoe Valley as a favorable resource basin that likely has both sedimentary and hydrothermal characteristics. However, there remains significant uncertainty on the nature and architecture of the resource(s) at depth, which increases the risk in exploratory drilling. Newly acquired gravity, magnetic, magnetotelluric, and controlled-source electromagnetic data, in conjunction with new and preceding geoscientific measurements and observations, are being integrated and evaluated in this study for efficacy in understanding stratigraphic geothermal resources and mitigating exploration risk. Furthermore, the influence of hydrothermal activity on sedimentary-hosted reservoirs in favorable structural settings (i.e., whether fault-controlled systems may locally enhance temperature and permeability in some deep stratigraphic reservoirs) will also be evaluated. This paper provides details and current updates on the course of this study in-progress.

Community, corporate, and government organizations are being targeted by disinformation attacks at an unprecedented rate. These attacks interrupt the ability of organizations to make high-consequence decisions and can lower their confidence in datasets and analytics. New interdisciplinary research approaches are being actively developed to expand resilience theory applications to organizations, and to determine the metrics and mitigations needed to increase resilience against disinformation. This paper presents initial ideas on adapting resilience methodologies for organizations and disinformation, highlighting key areas that require further exploration in this emerging field of research.

This paper presents Energy Storage-based Packetized Delivery of Electricity (ES-PDE) that is radically different from the operation of today's grid. Under ES-PDE, loads are powered by energy storage systems (ESS) most of the time and only receive packets of electricity periodically to power themselves and charge their ESSs. Therefore, grid operators can schedule the delivery of electricity in a manner that utilizes existing grid infrastructure. Since customers are powered by the co-located ESSs, when grid outages occur, they can be self-powered for some time before the grid is fully restored.In this paper, two operating schemes for ES-PDE are proposed. A Mixed-Integer-Linear-Programming (MILP) optimization is developed to find the optimal packet delivery schedule for each operating scheme. A case study is conducted to demonstrate the operation of ES-PDE.

The state of charge (SoC) estimated by Battery Management Systems (BMSs) could be vulnerable to False Data Injection Attacks (FDIAs), which aim to disturb state estimation. Inaccurate SoC estimation, due to attacks or suboptimal estimators, could lead to thermal runaway, accelerated degradation of batteries, and other undesirable events. In this paper, an ambient temperature-dependent model is adopted to represent the physics of a stack of three series-connected battery cells, and an Unscented Kalman Filter (UKF) is utilized to estimate the SoC for each cell. A Cumulative Sum (CUSUM) algorithm is used to detect FDIAs targeting the voltage sensors in the battery stack. The UKF was more accurate in state and measurement estimation than the Extended Kalman Filter (EKF) for Maximum Absolute Error (MAE) and Root Mean Squared Error (RMSE). The CUSUM algorithm described in this paper was able to detect attacks as low as ±1 mV when one or more voltage sensor was attacked under various ambient temperatures and attack injection times.

This manuscript presents the recent advances in Mixed-Integer Nonlinear Programming (MINLP) and Generalized Disjunctive Programming (GDP) with a particular scope for superstructure optimization within Process Systems Engineering (PSE). We present an environment of open-source software packages written in Python and based on the algebraic modeling language Pyomo. These packages include MindtPy, a solver for MINLP that implements decomposition algorithms for such problems, CORAMIN, a toolset for MINLP algorithms providing relaxation generators for nonlinear constraints, Pyomo.GDP, a modeling extension for Generalized Disjunctive Programming that allows users to represent their problem as a GDP natively, and GDPOpt, a collection of algorithms explicitly tailored for GDP problems. Combining these tools has allowed us to solve several problems relevant to PSE, which we have gathered in an easily installable and accessible library, GDPLib. We show two examples of these models and how the flexibility of modeling given by Pyomo.GDP allows for efficient solutions to these complex optimization problems. Finally, we show an example of integrating these tools with the framework IDAES PSE, leading to optimal process synthesis and conceptual design with advanced multi-scale PSE modeling systems.

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 violated, the existing system security margins will be largely unknown. As an alternative to the use of complex optimizations over dynamic conditions, this work employs the use of Reinforcement-based 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. The approach requires the training of a Machine Learning Agent to accomplish this task using modeled data and employs it as a decision support tool under severe, near-blackout conditions.

This paper presents a set of tests on a bolted benchmark structure called the S4 beam with a focus on evaluating coupling between the first two modes due to nonlinearity. Bolted joints are of interest in dynamically loaded structures because frictional slipping at the contact interface can introduce amplitude-dependent nonlinearities into the system, where the frequency of the structure decreases, and the damping increases. The challenge to model this phenomenon is even more difficult if the modes of the structure become coupled, violating a common assumption of mode orthogonality. This work presents a detailed set of measurements in which the nonlinearities of a bolted structure are highly coupled for the first two modes. Two nominally identical bolted structures are excited using an impact hammer test. The nonlinear damping curves for each beam are calculated using the Hilbert transform. Although the two structures have different frequency and damping characteristics, the mode coupling relationship between the first two modes of the structures is shown to be consistent and significant. The data is intended as a challenge problem for interested researchers; all data from these tests are available upon request.

The possibility of estimating the effective resistance at contact points along a seam in a cylindrical vessel is investigated. The vessel is formed from two top-hat structures bolted together at a flange. Aluminum shims at the bolt locations ensure a nearly constant 5-mil gap or slot between the flanges. Cavity modes are excited with a short monopole antenna inside the structure, and external near fields 5 mm away from the slot are probed around the vessel circumference. Comparison of CST and FDTD simulations with measurements reveals that the shape of the field-vs-angle curve is strongly dependent on the contact resistance, indicating that meaningful estimates can be extracted.

This paper uses co-located wind and photovoltaic generation, along with battery energy storage, as a single plant and introduces a method to provide a flexible synthetic inertia (SI) response based on plant-wide settings. The proposed controller accounts for variable resources and correctly adjusts device responses when an inverter-based resource (IBR) may become unavailable to provide a consistent plant level SI response. The flexible SI response is shown to adequately replace the lost synchronous inertial response from equivalent conventional generation when IBR penetration is approximately 25% in a small power system. Furthermore, it is shown that a high gain SI response provided by the combined IBR plant can reduce the rate of change of frequency magnitude over 50% from the equivalently rated conventional generation response.

This chapter deals with experimental dynamic substructures which are reduced order models that can be coupled with each other or with finite element derived substructures to estimate the system response of the coupled substructures. A unifying theoretical framework in the physical, modal or frequency domain is reviewed with examples. The major issues that have hindered experimental based substructures are addressed. An example is demonstrated with the transmission simulator method that overcomes the major historical difficulties. Guidelines for the transmission simulator design are presented.

Thermal runaway and its propagation are major safety issues in containerized lithium-ion battery energy storage systems. While conduction-driven propagation has received much attention, the thermal hazards associated with propagation via hot gases vented from failing cells are still not fully understood. Vented gases can lead to global safety issues in containerized systems, via heat transfer to other parts of the system and potential combustion hazards. In this work, we validate the characteristics of vented gases from cells undergoing thermal runaway in the thermal propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway). In particular, we assess the evolution of vented gases, venting time, and temperature profiles of single cell and multi-cell arrays based on experiments performed in Archibald et al (Fire Technology, 2020). While several metrics for estimating the venting time are assessed, a metric based on the CO2 generation results in consistent predictions. Vented gas evolution, and venting times predicted by the simulations are consistent with those estimated during the experiments. The simulation resolution and other model parameters, especially the use of an intra-particle diffusion limiter, have a large role in prediction of venting time.

In order to evaluate the time evolution of avalanche breakdown in wide and ultra-wide bandgap devices, we have developed a cable pulser experimental setup that can evaluate the time-evolution of the terminating impedance for a semiconductor device with a time resolution of 130 ps. We have utilized this pulser setup to evaluate the time-to-breakdown of vertical Gallium Nitride and Silicon Carbide diodes for possible use as protection elements in the electrical grid against fast transient voltage pulses (such as those induced by an electromagnetic pulse event). We have found that the Gallium Nitride device demonstrated faster dynamics compared to the Silicon Carbide device, achieving 90% conduction within 1.37 ns compared to the SiC device response time of 2.98 ns. While the Gallium Nitride device did not demonstrate significant dependence of breakdown time with applied voltage, the Silicon Carbide device breakdown time was strongly dependent on applied voltage, ranging from a value of 2.97 ns at 1.33 kV to 0.78 ns at 2.6 kV. The fast response time (< 5 ns) of both the Gallium Nitride and Silicon Carbide devices indicate that both materials systems could meet the stringent response time requirements and may be appropriate for implementation as protection elements against electromagnetic pulse transients.

We propose a new benchmark for high-performance (HP) computers. Similar to High Performance Conjugate Gradient (HPCG), the new benchmark is designed to rank computers based on how fast they can solve a sparse linear system of equations, exhibiting computational and communication requirements typical in many scientific applications. The main novelty of the new benchmark is that it is now based on Generalized Minimum Residual method (GMRES) (combined with Geometric Multi-Grid preconditioner and Gauss-Seidel smoother) and provides the flexibility to utilize lower precision arithmetic. This is motivated by new hardware architectures that deliver lower-precision arithmetic at higher performance. There are other machines that do not follow this trend. However, using a lower-precision arithmetic reduces the required amount of data transfer, which alone could improve solver performance. Considering these trends, an HP benchmark that allows the use of different precisions for solving important scientific problems will be valuable for many different disciplines, and we also hope to promote the design of future HP computers that can utilize mixed-precision arithmetic for achieving high application performance. We present our initial design of the new benchmark, its reference implementation, and the performance of the reference mixed (double and single) precision Geometric Multi-Grid solvers on current top-ranked architectures. We also discuss challenges of designing such a benchmark, along with our preliminary numerical results using 16-bit numerical values (half and bfloat precisions) for solving a sparse linear system of equations.

There are unique challenges associated with protection and self-healing of microgrids energized by multiple inverterbased distributed energy resources. In this study, prioritized undervoltage load shedding and undervoltage-supervised overcurrent (UVOC) for fault isolation are demonstrated using PSCAD. The PSCAD implementations of these relays are described in detail, and their operation in a self-healing microgrid is demonstrated.

Area efficient self-correcting flip-flops for use with triple modular redundant (TMR) soft-error hardened logic are implemented in a 12-nm finFET process technology. The TMR flip-flop slave latches self-correct in the clock low phase using Muller C-elements in the latch feedback. These C-elements are driven by the two redundant stored values and not by the slave latch itself, saving area over a similar implementation using majority gate feedback. These flip-flops are implemented as large shift-register arrays on a test chip and have been experimentally tested for their soft-error mitigation in static and dynamic modes of operation using heavy ions and protons. We show how high clock skew can result in susceptibility to soft-errors in the dynamic mode, and explain the potential failure mechanism.

As presented above, because similar existing DOE-managed SNF (DSNF) from previous reactors have been evaluated for disposal pathways, we use this knowledge/experience as a broad reference point for initial technical bases for preliminary dispositioning of potential AR SNF. The strategy for developing fully-formed gap analyses for AR SNF entails the primary step of first obtaining all the defining characteristics of the AR SNF waste stream from the AR developers. Utilizing specific and accurate information/data for developing the potential disposal inventory to be evaluated is a key principle start for success. Once the AR SNF waste streams are defined, the initial assessments would be based on comparison to appropriate existing SNF/waste forms previously analyzed (prior experience) to make a determination on feasibility of direct disposal, or the need to further evaluate due to differences specific to the AR SNF. Assessments of criticality potential and controls would also be performed to assess any R&D gaps to be addressed in that regard as well. Although some AR SNF may need additional treatment for waste form development, these aspects may also be constrained and evaluated within the context of disposal options, including detailed gap analysis to identify further R&D activities to close the gaps.

Abstract not provided.

High-performing teams learn intelligent and efficient communication and coordination strategies to maximize their joint utility. These teams implicitly understand the different roles of heterogeneous team members and adapt their communication protocols accordingly. Multi-Agent Reinforcement Learning (MARL) seeks to develop computational methods for synthesizing such coordination strategies, but formulating models for heterogeneous teams with different state, action, and observation spaces has remained an open problem. Without properly modeling agent heterogeneity, as in prior MARL work that leverages homogeneous graph networks, communication becomes less helpful and can even deteriorate the cooperativity and team performance. We propose Heterogeneous Policy Networks (HetNet) to learn efficient and diverse communication models for coordinating cooperative heterogeneous teams. Building on heterogeneous graph-attention networks, we show that HetNet not only facilitates learning heterogeneous collaborative policies per existing agent-class but also enables end-to-end training for learning highly efficient binarized messaging. Our empirical evaluation shows that HetNet sets a new state of the art in learning coordination and communication strategies for heterogeneous multi-agent teams by achieving an 8.1% to 434.7% performance improvement over the next-best baseline across multiple domains while simultaneously achieving a 200× reduction in the required communication bandwidth.

This paper provides a study of the potential impacts of climate change on intermittent renewable energy resources, battery storage, and resource adequacy in Public Service Company of New Mexico's Integrated Resource Plan for 2020 - 2040. Climate change models and available data were first evaluated to determine uncertainty and potential changes in solar irradiance, temperature, and wind speed in NM in the coming decades. These changes were then implemented in solar and wind energy models to determine impacts on renewable energy resources in NM. Results for the extreme climate-change scenario show that the projected wind power may decrease by ~13% due to projected decreases in wind speed. Projected solar power may decrease by ~4% due to decreases in irradiance and increases in temperature in NM. Uncertainty in these climateinduced changes in wind and solar resources was accommodated in probabilistic models assuming uniform distributions in the annual reductions in solar and wind resources. Uncertainty in battery storage performance was also evaluated based on increased temperature, capacity fade, and degradation in roundtrip efficiency. The hourly energy balance was determined throughout the year given uncertainties in the renewable energy resources and energy storage. The loss of load expectation (LOLE) was evaluated for the 2040 No New Combustion portfolio and found to increase from 0 days/year to a median value of ~2 days/year due to potential reductions in renewable energy resources and battery storage performance and capacity. A rank-regression analyses revealed that battery round-trip efficiency was the most significant parameter that impacted LOLE, followed by solar resource, wind resource, and battery fade. An increase in battery storage capacity to ~30,000 MWh from a baseline value of ~14,000 MWh was required to reduce the median value of LOLE to ~0.2 days/year with consideration of potential climate impacts and battery degradation.

Abstract not provided.

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged Navier–Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code-and solution-verification activities. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities.

We evaluate the use of reference modules for monitoring effective irradiance in PV power plants, as compared with traditional plane-of-array (POA) irradiance sensors, for PV monitoring and capacity tests. Common POA sensors such as pyranometers and reference cells are unable to capture module-level irradiance nonuniformity and require several correction factors to accurately represent the conditions for fielded modules. These problems are compounded for bifacial systems, where the power loss due to rear side shading and rear-side plane-of-array (RPOA) irradiance gradients are greater and more difficult to quantify. The resulting inaccuracy can have costly real-world consequences, particularly when the data are used to perform power ratings and capacity tests. Here we analyze data from a bifacial single-axis tracking PV power plant, (175.6 MWdc) using 5 meteorological (MET) stations, located on corresponding inverter blocks with capacities over 4 MWdc. Each MET station consists of bifacial reference modules as well pyranometers mounted in traditional POA and RPOA installations across the PV power plant. Short circuit current measurements of the reference modules are converted to effective irradiance with temperature correction and scaling based on flash test or nameplate short circuit values. Our work shows that bifacial effective irradiance measured by pyranometers averages 3.6% higher than the effective irradiance measured by bifacial reference modules, even when accounting for spectral, angle of incidence, and irradiance nonuniformity. We also performed capacity tests using effective irradiance measured by pyranometers and reference modules for each of the 5 bifacial single-axis tracking inverter blocks mentioned above. These capacity tests evaluated bifacial plant performance at ∼3.9% lower when using bifacial effective irradiance from pyranometers as compared to the same calculation performed with reference modules.

Abstract not provided.

Abstract not provided.

We develop methods that could be used to qualify a training dataset and a data-driven turbulence closure trained on it. By qualify, we mean identify the kind of turbulent physics that could be simulated by the data-driven closure. We limit ourselves to closures for the Reynolds-Averaged Navier Stokes (RANS) equations. We build on our previous work on assembling feature-spaces, clustering and characterizing Direct Numerical Simulation datasets that are typically pooled to constitute training datasets. In this paper, we develop an alternative way to assemble feature-spaces and thus check the correctness and completeness of our previous method. We then use the characterization of our training dataset to identify if a data-driven turbulence closure learned on it would generalize to an unseen flow configuration – an impinging jet in our case. Finally, we train a RANS closure architected as a neural network, and develop an explanation i.e., an interpretable approximation, using generalized linear mixed-effects models and check whether the explanation resembles a contemporary closure from turbulence modeling.

We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poisson’s and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints, and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both Poisson’s problem and the linear peridynamic solid model.

A crucial component of field testing is the utilization of numerical models to better understand the system and the experimental data being collected. Meshing and modeling field tests is a complex and computationally demanding problem. Hexahedral elements cannot always reproduce experimental dimensions leading to grid orientation or geometric errors. Voronoi meshes can match complex geometries without sacrificing orthogonality. As a result, here we present a high-resolution 3D numerical study for the BATS heater test at the WIPP that compares both a standard non-deformed cartesian mesh along with a Voronoi mesh to match field data collected during a salt heater experiment.

The DOE GeoVision study identified that Enhanced Geothermal Systems (EGS) resources have the potential to provide a significant contribution toward achieving the goal of converting the U.S. electricity system to 100% clean energy over the next few decades. To further the implementation of commercial EGS development, DOE's Geothermal Technologies Office (GTO) initiated the Wells of Opportunity (WOO) Amplify program, where unproductive wells in selected geothermal fields are to be stimulated using EGS technologies, resulting in increased power production from these resources. As part of the WOO-Amplify project, GTO assembled the Amplify Monitoring Team (AMT), whose role is to provide in-field and near-field seismic monitoring design, deployment and data analysis for stimulations under the WOO-Amplify initiative. This team, consisting of scientists and engineers from Lawrence Berkeley National Laboratory (LBNL), Sandia National Laboratories (SNL), and the US Geological Survey (USGS), is working with WOO-Amplify EGS Operators in Nevada to develop and deploy optimized seismic monitoring systems at four geothermal fields where WOO-Amplify well stimulations are planned: Don A. Campbell, Tungsten Mountain and Jersey Valley operated by Ormat Technologies, and Patua operated by Cyrq Patua Acquisition Company LLC. Using geologic and geophysical field data provided by the WOO-Amplify teams, the focus of the AMT is to develop advanced simulation and modeling techniques, design targeted seismic monitoring arrays, develop innovative and cost-effective methodologies for drilling seismic monitoring boreholes, deploy effective seismic instrumentation, and facilitate the use of microseismic data to monitor well stimulation and flow within the geothermal reservoir. Realtime seismic data from the four WOO-Amplify sites will be streamed to a publicly accessible Amplify Monitoring website. AMT's advanced simulations and template matching techniques applied during pre-stimulation phases can help improve understanding of potential seismic hazard and inform the Operator's Induced Seismicity Mitigation Protocol (ISMP). Over the next two years, AMT will be drilling, instrumenting, and recording seismic data at the WOO-Amplify field sites, telemetering the seismic waveform data to AMT's central processing system and providing the processed location data to the WOO Amplify Operator teams. These data and monitoring systems will be critical for effective monitoring of the effects of planned well stimulation and extended flow tests during the next stage of the WOO-Amplify project.

Abstract not provided.

We measured the Hugoniot, Hugoniot elastic limit (HEL), and spallation strength of laser powder bed fusion (LPBF) AlSi10Mg via uniaxial plate-impact experiments to stresses greater than 13 GPa. Despite its complex anisotropic microstructure, the LPBF AlSi10Mg did not exhibit significant orientation dependence or sample-to-sample variability in these measured quantities. We found that the Hugoniot response of the LPBF AlSi10Mg is similar to that of other Al-based alloys and is well approximated by a linear relationship: u_s = 5.49 + 1.39u_p. Additionally, the measured HELs ranged from 0.25 to 0.30 GPa and spallation strengths ranged from 1.16 to 1.45 GPa, consistent with values reported in other studies of LPBF AlSi10Mg and Al-based alloys. Furthermore, strain-rate and stress dependence of the spallation strength were also observed.

Sandia National Laboratories has been tasked to operate and maintain the National Solar Thermal Test Facility (NSTTF) located in Kirtland Airforce Base near Albuquerque, NM. The NSTTF provides established test platforms and experienced researchers and technologists in the field of Concentrating Solar Technologies (CST) and Concentrating Solar Power (CSP). This three-year project seeks to maintain the NSTTF for development, testing, and application of new CSP technologies that are instrumental in advancing the state-of-the-art in support of SunShot and Generation 3 CSP technology goals. In turn, these technologies will form the foundation of the global CSP industry and continue to advance the technology to new levels of efficiency, higher temperatures, lower costs, lower risk, and higher reliability. The NSTTF provides established test platforms and highly experienced researchers and technologists in the CSP field.

Lithium/fluorinated graphite (Li/CF _x ) primary batteries show great promise for applications in a wide range of energy storage systems due to their high energy density (>2100 Wh kg ^–1 ) and low self‐discharge rate (<0.5% per year at 25 °C). While the electrochemical performance of the CF _x cathode is indeed promising, the discharge reaction mechanism is not thoroughly understood to date. In this article, a multiscale investigation of the CF _x discharge mechanism is performed using a novel cathode structure to minimize the carbon and fluorine additives for precise cathode characterizations. Titration gas chromatography, X‐ray diffraction, Raman spectroscopy, X‐ray photoelectron spectroscopy, scanning electron microscopy, cross‐sectional focused ion beam, high‐resolution transmission electron microscopy, and scanning transmission electron microscopy with electron energy loss spectroscopy are utilized to investigate this system. Results show no metallic lithium deposition or intercalation during the discharge reaction. Crystalline lithium fluoride particles uniformly distributed with <10 nm sizes into the CF _x layers, and carbon with lower sp ² content similar to the hard‐carbon structure are the products during discharge. This work deepens the understanding of CF _x as a high energy density cathode material and highlights the need for future investigations on primary battery materials to advance performance.

Abstract not provided.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User's Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user's notes and of course the material in the open literature.

The How To Manual supplements the User’s Manual and the Theory Manual. The goal of the How To Manual is to reduce learning time for complex end to end analyses. These documents are intended to be used together. See the User’s Manual for a complete list of the options for a solution case. All the examples are part of the Sierra/SD test suite. Each runs as is. The organization is similar to the other documents: How to run, Commands, Solution cases, Materials, Elements, Boundary conditions, and then Contact. The table of contents and index are indispensable. The Geometric Rigid Body Modes section is shared with the Users Manual.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user's guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.