The Energetic Neutrons campaign led by Sandia National Laboratories (SNL) had a successful year testing electronic devices and printed circuit boards (PCBs) under 14 MeV neutron irradiation at OMEGA. During FY21 Sandia’s Neutron Effects Diagnostics (NEDs) and data acquisition systems were upgraded to test novel commercial off-the-shelf and Sandia-fabricated electronic components that support SNL’s National Security mission. The upgrades to the Sandia platform consisted of new cable chains, sample mount fixtures and a new fiber optics platform for testing optoelectronic devices.
Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.
Bar-Nes, Gabriela B.; Klein-BenDavid, Ofra K.; Kosson, David K.; Gruber, Chen G.; Taylor, Autumn T.; Brown, Kevin B.; Delapp, Rossane D.; Brown, Lesa B.; Ayers, John A.; Meeusen, Hans M.; Matteo, Edward N.; Jove Colon, Carlos F.; Mitchell, Chven A.; Pyrak-Nolte, Laura J.
The Computer Science Research Institute (CSRI) brings university faculty and students to Sandia National Laboratories for focused collaborative research on Department of Energy (DOE) computer and computational science problems. The institute provides an opportunity for university researches to learn about problems in computer and computational science at DOE laboratories, and help transfer results of their research to programs at the labs. Some specific CSRI research interest areas are: scalable solvers, optimization, algebraic preconditioners, graph-based, discrete, and combinatorial algorithms, uncertainty estimation, validation and verification methods, mesh generation, dynamic load-balancing, virus and other malicious-code defense, visualization, scalable cluster computers, beyond Moore’s Law computing, exascale computing tools and application design, reduced order and multiscale modeling, parallel input/output, and theoretical computer science. The CSRI Summer Program is organized by CSRI and includes a weekly seminar series and the publication of a summer proceedings.
This report describes an assessment of flamelet based soot models in a laminar ethylene coflow flame with a good selection of measurements suitable for model validation. Overall flow field and temperature predictions were in good agreement with available measurements. Soot profiles were in good agreement within the flame except for near the centerline where imperfections with the acetylene-based soot-production model are expected to be greatest. The model was challenged to predict the transition between non-sooting and sooting conditions with non-negligible soot emissions predicted even down to small flow rates or flame sizes. This suggests some possible deficiency in the soot oxidation models that might alter the amount of smoke emissions from flames, though this study cannot quantify the magnitude of the effect for large fires.
The TChem open-source software is a toolkit for computing thermodynamic properties, source term, and source term’s Jacobian matrix for chemical kinetic models that involve gas and surface reactions.
This report presents an analysis of the performance of deployable energy systems comprised of wind energy systems integrated with diesel generators, photovoltaic systems, and battery storage to meet the load requirements of a representative U.S. Army forward operating base. The analysis is conducted using HOMER, a microgrid analysis software that can search through a wide range of parameters to design and optimize microgrid power systems. The search parameters include the system architecture, the wind and solar resources, and the availability of diesel fuel. The results of the analysis measure the relative performance of the different systems and environments in terms of the overall transportation cost to deploy the system and the ability to provide resilience in terms of meeting mission critical loads.
At the Z Facility at Sandia National Laboratories, the magnetized liner inertial fusion (MagLIF) program aims to study the inertial confinement fusion in deuterium-filled gas cells by implementing a three-step process on the fuel: premagnetization, laser preheat, and Z-pinch compression. In the laser preheat stage, the Z-Beamlet laser focuses through a thin polyimide window to enter the gas cell and heat the fusion fuel. However, it is known that the presence of the few μm thick window reduces the amount of laser energy that enters the gas and causes window material to mix into the fuel. These effects are detrimental to achieving fusion; therefore, a windowless target is desired. The Lasergate concept is designed to accomplish this by "cutting"the window and allowing the interior gas pressure to push the window material out of the beam path just before the heating laser arrives. In this work, we present the proof-of-principle experiments to evaluate a laser-cutting approach to Lasergate and explore the subsequent window and gas dynamics. Further, an experimental comparison of gas preheat with and without Lasergate gives clear indications of an energy deposition advantage using the Lasergate concept, as well as other observed and hypothesized benefits. While Lasergate was conceived with MagLIF in mind, the method is applicable to any laser or diagnostic application requiring direct line of sight to the interior of gas cell targets.
The myriad detectors and instruments now available for materials characterization provide researchers with an ever-growing suite of tools to probe material behavior. Progress in the development of instrumentation and workflows that enable the collection, and leverage the potential, of various data modalities have provided novel insights into material behavior. Using data across multiple length scales, or performing complementary analyses of in situ and ex situ data, can help reveal a more complete picture of dynamic processes or material structure. However, the accurate combination, or fusion, of these disparate data modalities presents new challenges. Differences in resolution, as well as the varying length scales at which physical phenomena are exploited to generate these data, necessitate novel approaches to accurately interpret and combine these data. Furthermore, the papers within this special topic focus on the collection and fusion of multimodal data to better understand structural materials. From new frameworks and workflows for data segmentation and analysis, process monitoring, enhancing simulations, or interrogating mechanical response, these papers reveal the potential benefits of utilizing multimodal data.
Liquefied petroleum gas (LPG) is a viable, cleaner alternative to traditional diesel fuel used in busses and other heavy-duty vehicles and could play a role in helping the US meet its lower emission goals. While the LPG industry has focused efforts on developing vehicles and fueling infrastructure, we must also establish safe parameters for maintenance facilities which are servicing LPG fueled vehicles. Current safety standards aid in the design of maintenance facilities, but additional quantitative analysis is needed to prove safeguards are adequate and suggest improvements where needed. In this report we aim to quantify the amount of flammable mass associated with propane releases from vehicle mounted fuel vessels within enclosed garages. Furthermore, we seek to qualify harm mitigation with variable ventilations and facility layout. To accomplish this we leverage validated computational resources at Sandia National Laboratories to simulate various release scenarios representative of real world vehicles and maintenance facilities. Flow solvers are used to predict the dynamics of fuel systems as well as the evolution of propane during release events. From our simulated results we observe that both inflow and outflow ventilation locations play a critical role in reducing flammable cloud size and potential overpressure values during a possible combustion event.
This white paper describes the work performed by Sandia National Laboratories in the New Mexico Small Business Agreement with BayoTech. BayoTech is a hydrogen generation and distribution company that is located in Albuquerque, NM. Their goal is to distribute hydrogen via their hydrogen systems which utilize the core design that was developed by Sandia. However, because the hydrogen economy is in its nascency, the safety and operation of the generating systems require independent validation. Additionally, in their pursuit of permitting at various locations around the nation, they require fire protection engineering support in discussions with local fire marshals and neighboring industrial entities. Sandia National Laboratories has subject matter expertise in hydrogen risk modeling of consequence (overpressure and dispersion) as well as fire protection engineering. Throughout this project, Sandia has worked with BayoTech to provide our expertise in these subject areas to facilitate the market entry of their hydrogen generation project to address the dire need for decarbonization due to climate change. The general approach of the support by Sandia is outlined in the main body, while the location specific evaluation for the Port of Stockton is contained in Appendix A.
Simultaneous data of the quasi-static compaction and electrical conductivity of porous, binary powder mixtures have been collected as a function of bulk density. The powder mixtures consist of a metal conductor, either titanium or iron, an insulator, and pores filled with ambient air. The data show a dependency of the conductivity in terms of relative bulk density and metal volume fraction on conductor type and conductor particle characteristics of size and shape. Finite element models using particle domains generated by discrete element method are used to simulate the bulk conductivity near its threshold while the general effective media equation is used to model the conductivity across the compression regime.
The U.S. Department of Energy Solar Energy Technologies Office initiated the Generation 3 Concentrating Solar Power (CSP) program to achieve higher operating temperatures (>700 °C) to enable next-generation CSP high-temperature power cycles such as the supercritical CO2 (sCO2) Brayton Cycle. Three teams were selected to pursue high-temperature gas, liquid, and solid pathways for the heat-transfer media. Phases 1 and 2, which lasted from 2018 – 2020, consisted of design, modeling, and testing activities to further de-risk each of the technologies and develop a design for construction, commissioning, and operation of a pilot-scale facility in Phase 3 (2021 – 2024). This report summarizes the activities in Phases 1 and 2 for the solid-particle pathway led by Sandia National Laboratories. In Phases 1 and 2, Sandia successfully de-risked key elements of the proposed Gen 3 Particle Pilot Plant (G3P3) by improving the design, operation, and performance of key particle component technologies including the receiver, storage bins, particle-to-sCO2 heat exchanger, particle lift, and data acquisition and controls. Modeling and testing of critical components have led to optimized designs that meet desired performance metrics. Detailed drawings, piping and instrumentation diagrams, and process flow diagrams were generated for the integrated system, and structural analyses of the assembled tower structure were performed to demonstrate compliance with relevant codes and standards. Instrumentation and control systems of key subsystems were also demonstrated. Together with Bridgers & Paxton, Bohannan Huston, and Sandia Facilities, we have completed a 100% G3P3 tower design package with stamped engineering drawings suitable for construction bid in Phase 3.
The following report contains data and data summaries collected for the SkySun LLC elevated Ganged PV arrays. These arrays were fabricated as a series of PV panels in various orientations, suspended by cables, at the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories (SNL). Starting in February of 2021, Sandia personnel have collected power and accelerometer data for these arrays to assess design and operational efficacy of varying ganged- PV configurations. The purpose of this power data collection was to see how the various array orientations compare in power collection capability depending on the time of day, year, and the specific daily solar direct normal irradiance (DNI). The power data was collected as a measurement of the power output from the various series strings. The project team measured direct current (DC) voltage and current from the respective arrays. The accelerometer data was collected with the purpose of demonstrating potential destructive mode shapes that could take place with each of the arrays when exposed to high winds. This allowed the team to evaluate whether impacts with respect to specific array orientations using suspended cables is a safe design. All data collection was performed during calendar year 2021.
Naughton, Brian T.; Jimenez, Tony; Preus, Robert; Summerville, Brent; Whipple, Bradley; Reen, Dylan; Gentle, Jake; Lang, Eric
This document aims to provide guidance on the design and operation of deployable wind systems that provide maximum value to missions in defense and disaster relief. Common characteristics of these missions are shorter planning and execution time horizons and a global scope of potential locations. Compared to conventional wind turbine applications, defense and disaster response applications place a premium on rapid shipping and installation, short-duration operation (days to months), and quick teardown upon mission completion. Furthermore, defense and disaster response applications are less concerned with cost of energy than conventional wind turbine applications. These factors impart design drivers that depart from the features found in conventional distributed wind turbines, thus necessitating unique design guidance. The supporting information for this guidance comes from available relevant references, technical analyses, and input from industry and military stakeholders. This document is not intended to be a comprehensive, prescriptive design specification. This document is intended to serve as a written record of an ongoing discussion of stakeholders about the best currently available design guidance for deployable wind turbines to help facilitate the effective development and acquisition of technology solutions to support mission success. The document is generally organized to provide high-level, focused guidance in the main body, with more extensive supporting details available in the referenced appendices. Section 2 begins with a brief qualitative description of the design guidelines being considered for the deployable wind turbines. Section 3 provides an overview of the characteristics of the mobile power systems commonly used in U.S. military missions. Section 4 covers current military and industry standards and specifications that are relevant to a deployable wind turbine design. Section 5 presents the deployable turbine design guidelines for the application cases.
Adams, Brian H.; Bohnhoff, William J.; Dalbey, Keith R.; Ebeida, Mohamed S.; Eddy, John P.; Eldred, Michael S.; Hooper, Russell W.; Hough, Patricia D.; Hu, Kenneth T.; Jakeman, John D.; Khalil, Mohammad; Maupin, Kathryn A.; Monschke, Jason A.; Ridgway, Elliott M.; Rushdi, Ahmad A.; Seidl, Daniel T.; Stephens, John A.; Swiler, Laura P.; Laros, James H.; Winokur, Justin G.
The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a user's manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.
Theristis, Marios; Livera, Andreas; Micheli, Leonardo; Ascencio-Vasquez, Julian; Makrides, George; Georghiou, George E.; Stein, Joshua S.
A linear performance drop is generally assumed during the photovoltaic (PV) lifetime. However, operational data demonstrate that the PV module degradation rate (Rd) is often nonlinear, which, if neglected, may increase the financial uncertainty. Although nonlinear behavior has been the subject of numerous publications, it was only recently that statistical models able to detect change-points and extract multiple Rd values from PV performance time-series were introduced. A comparative analysis of six open-source libraries, which can detect change-points and calculate nonlinear Rd, is presented in this article. Since the real Rd and change-point locations are unknown in field data, 960 synthetic datasets from six locations and two PV module technologies have been generated using different aggregation and normalization decisions and nonlinear degradation rate patterns. The results demonstrated that coarser temporal aggregation (i.e., monthly vs. weekly), temperature correction, and both PV module technologies and climates with lower seasonality can benefit the change-point detection and Rd extraction. This also raises a concern that statistical models typically deployed for Rd analysis may be highly climatic-and technology-dependent. The comparative analysis of the six approaches demonstrated median mean absolute errors (MAE) ranging from 0.06 to 0.26%/year, given a maximum absolute Rd of 2.9%/year. The median MAE in change-point position detection varied from 3.5 months to 6 years.
Single case comparisons between severe accident simulations can provide detailed insights into severe accident model behavior, however, they cannot offer insights into model uncertainty, sensitivity to uncertain parameters, or underlying model biases. In this analysis, the single case benchmark comparison of the MELCOR material interaction models for a station blackout (SBO) scenario of a boiling water reactor (BWR) using representative Fukushima Daiichi Unit 1 boundary conditions is expanded to include an uncertainty analysis. As part of this uncertainty analysis, 1200 simulations are performed for each material interaction model (2400 total), with random sampling of 14 uncertain MELCOR input parameters. Input parameters are selected for their impact on models representing core degradation processes. These include candling, fuel rod failure, debris quenching and dryout. The analysis performed here is not a traditional “best-estimate” uncertainty analysis that uses best-estimate parameters or identifies best-estimate figure of merit distributions. Instead, it is an exploratory uncertainty analysis that identifies and interrogates underlying model form biases of the two material interaction models (eutectics and interactive materials models). Uniform distributions are applied to all uncertain parameters to ensure coverage of the model parameter uncertainty space. Key findings from this study include underlying model form biases exhibited by material interaction models, and notable differences in accident progression outcomes between the material interaction models. This uncertainty study extends and confirms the conclusions from the first part of this study, which compared the impact of material interaction modeling on simulation of a short-term station blackout scenario with representative Fukushima Daiichi Unit I boundary conditions. In particular, this study confirms that the eutectics model generally exhibits accelerated degradation and failure of fuel components, the core plate, and the lower head. The eutectics model also has a tendency to exhibit a greater degree of core degradation, greater debris mass formation, and larger debris mass ejection. Finally, the eutectics model exhibits higher maximum temperatures for fuel, cladding, particulate debris, oxidic molten pool, and metallic molten pool components than the interactive materials model; interactive materials model simulations exhibit a soft “limitation” on maximum temperatures that is related to the temperature at which material relocation occurs.
The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane infrastructure and transportation systems. HyRAM+ is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ includes generic probabilities for equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of release and flame physics. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impact on people and structures. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards. HyRAM+ is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM+ version 4.0. The most significant change for HyRAM+ version 4.0 from HyRAM version 3.1 is the incorporation of other alternative fuels, namely methane (as a proxy for natural gas) and propane into the toolkit. This change necessitated significant changes to the installable graphical user interface as well as changes to the back-end Python models. A second major change is the inclusion of physics models for the overpressure associated with the delayed ignition of an unconfined jet/plume of flammable gas.
Acceptance risk may be estimated using cumulative probability distribution functions applied to populations of measurands and test measurement values. “Simple acceptance” is a decision rule that sets the acceptance range of a test result equal to the tolerance range specification. While acceptance risk is comprised of both consumer risk and producer risk, this paper compares the effects of simple acceptance decision rules and guardbanding decision rules on consumer risk. Consumer risk is also known as the probability of false acceptance. The terms describe the risk of accepting test results as passing when the actual values exceed specification limits. False acceptance is only possible when the true value of a measurand is out of tolerance and the test result indicates that the measurand is within tolerance. Metrologists generally have some information about the test uncertainty regarding a specific acceptance test result. Along with a general lack of knowledge regarding the measurand population parameters, the complicated interplay between risk, dispersion, and central values generally prevents calibration laboratories from fully characterizing acceptance probabilities. Organizations that model the measurand populations can reduce consumer risk by avoiding certification of items whose measurand populations are not well centered. The models in this paper present a recurring trend: guardbanding reduces nonnegligible risks of false acceptance when compared to simple acceptance. However, guardbanding is not the most effective means of acceptance risk mitigation. Systematic characterization of measurand populations can provide the information a calibration laboratory needs to reliably control acceptance risk.
This work, building on previous efforts, develops a suite of new graph neural network machine learning architectures that generate data-driven prolongators for use in Algebraic Multigrid (AMG). Algebraic Multigrid is a powerful and common technique for solving large, sparse linear systems. Its effectiveness is problem dependent and heavily depends on the choice of the prolongation operator, which interpolates the coarse mesh results onto a finer mesh. Previous work has used recent developments in graph neural networks to learn a prolongation operator from a given coefficient matrix. In this paper, we expand on previous work by exploring architectural enhancements of graph neural networks. A new method for generating a training set is developed which more closely aligns to the test set. Asymptotic error reduction factors are compared on a test suite of 3-dimensional Poisson problems with varying degrees of element stretching. Results show modest improvements in asymptotic error factor over both commonly chosen baselines and learning methods from previous work.
It is hypothesized that dark matter is composed of particles called quark nuggets, and further that these particles have a permanent magnetic dipole moment. If this hypothesis is true, calculations predict that a magnetized quark nugget (MQN) will oscillate when encountering the Earth's magnetosphere, and emit RF radiation between 30kHz and 30MHz. To support testing this hypothesis, a loop antenna sensor was designed and developed, which is described in this report. This sensor operates between 300kHz and 3MHz and achieves about -11dBfT/vHz sensitivity at 1.5MHz.