SEARCH · Engineering Papers

Results for “Variational”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

58 GEOSCIENCES↗

DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗

DOE OSTI2024

Long-Term Infrasound Sensor Calibration and Characterization

Previous testing has shown that infrasound sensors deployed in the field can exhibit notable deviations from their nominal, lab-based calibrations. These variations may be due to changes in environmental conditions, long-term sensor drift, or other unresolved features. In early 2018, we installed two identical infrasound elements with five infrasound sensors at each element (Chaparral M50A, Chaparral M64LN, CEA/Martec MB2005, CEA/Seismowave MB3a, and Hyperion IFS-5113A). These sensors were accepted or under consideration for use in the International Monitoring System network of the Comprehensive Nuclear-Test-Ban Treaty. Each element had all sensors connected to a single digitizer and port to the atmosphere. We also recorded instrument enclosure air temperature and humidity and external air temperature. Using the MB2005 as the reference, we examine the relative sensor response (both magnitude and phase) as a function of time and frequency and compare it with quarterly laboratory calibrations and environmental conditions. Here, we find that the magnitude response for all sensors exhibits some variability in both the lab and field, with the amplitude variations often >5%. The field-based variations are more severe and occur on both long-term (months) and short-term (hours) timescales. Short-term variability correlates with changes in environmental conditions and is considerable (up to 25%) for the Chaparral M50A and noticeable (∼5%) for the French Alternative Energies and Atomic Energy Commission (CEA) MB3a. Long-term magnitude variability for the Chaparral M50A was up to 50% and does not clearly correlate with environmental conditions. The other sensors show some long-term magnitude offsets, but they have relatively stable responses in the conditions we examined. The MB3a also displays some frequency-dependent magnitude variability and shows a minor dependence on temperature. Phase estimates are stable and near zero for all sensors tested. These results strongly suggest sensor response and variability due to environmental conditions should be considered in future infrasound data interpretation and sensor selection and development.

Fee, David↗

DOE OSTI2021

Learning the ground state of a non-stoquastic quantum Hamiltonian in a rugged neural network landscape

Strongly interacting quantum systems described by non-stoquastic Hamiltonians exhibit rich low-temperature physics. Yet, their study poses a formidable challenge, even for state-of-the-art numerical techniques. Here, we investigate systematically the performance of a class of universal variational wave-functions based on artificial neural networks, by considering the frustrated spin- 1/2 1 / 2 J_1-J_2 J 1 − J 2 Heisenberg model on the square lattice. Focusing on neural network architectures without physics-informed input, we argue in favor of using an ansatz consisting of two decoupled real-valued networks, one for the amplitude and the other for the phase of the variational wavefunction. By introducing concrete mitigation strategies against inherent numerical instabilities in the stochastic reconfiguration algorithm we obtain a variational energy comparable to that reported recently with neural networks that incorporate knowledge about the physical system. Through a detailed analysis of the individual components of the algorithm, we conclude that the rugged nature of the energy landscape constitutes the major obstacle in finding a satisfactory approximation to the ground state wavefunction, and prevents learning the correct sign structure. In particular, we show that in the present setup the neural network expressivity and Monte Carlo sampling are not primary limiting factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Utility Sector Impacts of Reduced Electricity Demand: Updates to Methodology and Results

This report presents revisions to the methodology for estimating the utility sector impacts associated with electricity demand reductions, first published in our 2014 study. The earlier report described how to calculate utility impact factors based on data from the Energy Information Agency's Annual Energy Outlook (AEO) projections of electricity supply and demand. The 2014 methodology relied on data from the AEO Reference case, and side cases implementing demand-side policies with a substantial impact on total demand. Beginning in 2015, the AEO publication schedule was revised to limit the number of side cases published with each AEO, and no suitable demand-focused side cases have been released since 2014. Hence, this paper presents a generalization of the methodology that allows a uniform approach irrespective of the nature of the side cases published with the AEO. The impact factors we calculate represent coeffcients of variation between a unit reduction to site electricity demand and corresponding supply-side changes, including reductions to electric power generation and primary fuel consumption by fuel type, emissions of criteria pollutants (Hg, NO x and SO 2 ) and greenhouse gases (CO 2 , N 2 O and CH 4 ), and total installed capacity by plant type. The impact factors are time-dependent, and also depend on the demand sector and end use. In this report, we apply the methodology to all AEO publications from 2014-2019, and all scenarios. We quantify the variation in factors across different side cases within each AEO, and across different AEO publications, and compare this to the variation due to changes in methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Technical Gaps in Hanford High-Level Waste Solids Settling Behavior and Settling Time Evaluation for Direct Feed High-Level Waste (DFHLW) Operations

Settling of high-level waste (HLW) solids in process vessels is a key conceptual process step in providing HLW feed directly to the Hanford Waste Treatment and Immobilization Plant (WTP) HLW Vitrification Facility. Direct Feed High-Level Waste (DFHLW) is a potential approach to initiating HLW vitrification prior to completing of the WTP Pretreatment Facility. Settling would be used with subsequent supernatant decant to concentrate HLW feed. To support planning for DFHLW, Washington River Protection Solutions (WRPS) requested support from the Pacific Northwest National Laboratory to evaluate the current data set available to predict the time needed for HLW solids to settle, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling time. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers including: Gap 1: In-Tank Settling Rates Faster than Settling of Laboratory Samples, Gap 2: Effect of Sludge Leaching/Washing on Predicted Settling Times, Gap 3: Predicting Waste Settling from Waste Chemistry (Waste Type), Gap 4: Predicting Waste Settling from Particle Size and Density Distributions (PSDDs), Gap 5: Insufficient Laboratory and In-Tank Settling Data to Represent Hanford Waste, Gap 6: Methods for Real-Time, In-Tank Tracking of Settling, Gap 7: Prediction of Sediment Erosion Resistance as a Function of Settling Time, and Gap 8: Prediction of Sediment Solids Content as a Function of Settling Time. In addition to the data gaps, an overarching observation of the settling rate and settled layer data is the significant variation in behavior. At similar solids concentrations, settling rates can vary by as much as 3 orders of magnitude depending on the source waste tank, and significantly different settling rates are noted between laboratory and in situ tests for the same waste tank. The range of average solids concentration in existing HLW sediment, which may have been quiescent for decades, can vary from less than 7 wt% to greater than 74 wt% solids. The shear strengths (or yield stresses) measured on laboratory samples range from less than 27 Pa to greater than 6400 Pa. These variations can challenge process planning for the application of a settle/decant process for DFHLW. This report describes the significance of the gaps to the settle/decant process and presents uncertainties by way of examples. Potential technical approaches for resolving these gaps are described and the estimated difficulty in resolving these gaps is evaluated. Based on the significance of the gap and the difficulty of resolution, recommendations are made to address specific gaps. Scoping estimates of the potential settling times for DFHLW solids have been made based on the existing data set with its associated gaps. Depending on the process vessel depth and final sediment concentration, substantial fractions of the scoping estimate results for settling times for characterized HLW exceed the 2-week period that has been previously assumed for process planning. There is also significant disparity, potentially greater than a factor of 5000 difference, in the estimated settling times depending on process vessel depth and final sediment solid concentration. This variation in results underscores the significance of the identified gaps and uncertainties with respect to process planning for utilizing settle/decant operations for DFHLW.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Paleo-Megadroughts and Abrupt Climate Changes in the Speleothem Records. Final report

This project is motivated by the speleothem isotope records in Asia, which show regional responses in the hydrologic cycle to different climate forcings. Speleothem isotopic records are typically interpreted in terms of local precipitation variations or monsoon intensity. Our study demonstrates that non-local processes also play an important role. We started this project to understand the regional difference in speleothem isotopic composition between the Last Glacial Maximum (LGM) and the present-day. The record in Southwest China showed greater depletion during the LGM compared to those in East China. Our modeling and analysis showed that speleothems record, in addition, large scale changes in atmospheric circulation and moisture transport and their subsequent impact on precipitation. We developed an algorithm to partition total precipitation according to their formation dynamics, namely into frontal and non-frontal precipitation, and showed that the two have different trends and hence different causal mechanism. We then focused our subsequent attention on circulation impacts on precipitation changes. We applied a machine learning algorithm to detect rainbands in the ERA-Interim reanalysis product, and showed that the seasonal migrations of the rainbands are tied to the seasonal migrations of the jet stream, in particular the northerlies of the jet meanders. These northerlies, in turn, are partly topographic Rossby waves excited as the upstream westerlies impinge on the Tibetan Plateau. The seasonal variations of these upstream westerlies thus contribute to the seasonal movements of the rainbands and regional precipitation changes. Our analysis of the modern precipitation isotope record further confirms the importance of jet stream changes in the isotopic variations and shows that isotope-enriched years have reduced summer seasonality, with less pronounced northward migration of the jet.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Noisy Intermediate-Scale Quantum Applications on a Pathfinder System

Work performed under this one-year LDRD was concerned with estimating resource requirements for small quantum test beds that are expected to be available in the near future. This work represents a preliminary demonstration of our ability to leverage quantum hardware for solving small quantum simulation problems in areas of interest to the DOE. The algorithms enabling such studies are hybrid quantum-classical variational algorithms, in particular the widely-used variational quantum eigensolver (VQE). Employing this hybrid algorithm, in which the quantum computer complements the classical one, we implemented an end-to-end application-level toolchain that allows the user to specify a molecule of interest and compute the ground state energy using the VQE approach. We found significant limitations attributable to the classical portion of the hybrid system, including a greater than greater-than-quartic power scaling of the classical memory requirements with the system size. Current VQE approaches would require an exascale machine for solving any molecule with size greater than 150 nuclei. Our findings include several improvements that we implemented into the VQE toolchain, including a new classical optimizer that is decades old but hadn't been considered before in the VQE ecosystem. Our findings suggest limitations to variational hybrid approaches to simulation that further motivate the need for a gate-based fault-tolerant quantum processor that can implement larger problems using the fully digital quantum phase estimation algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Sensitivity and Uncertainty of the IFR-1 BISON Benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This report summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment. For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. The resulting variations in benchmark predictions were normalized to calculate sensitivities. The strongest sensitivities were identified and used to inform input selections for the UQ. The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because of the wide confidence intervals associated with the benchmark predictions. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark development and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Comprehensive Review of Multi-arm Caliper Data for the Big Hill SPR Site

The Big Hill SPR site has a rich data set consisting of multi-arm caliper (MAC) logs collected from the cavern wells. This data set provides insight into the on-going casing deformation at the Big Hill site. This report summarizes the MAC surveys for each well and presents well longevity estimates where possible. Included in the report is an examination of the well twins for each cavern and a discussion on what may or may not be responsible for the different levels of deformation between some of the well twins. The report also takes a systematic view of the MAC data presenting spatial patterns of casing deformation and deformation orientation in an effort to better understand the underlying causes. The conclusions present a hypothesis suggesting the small-scale variations in casing deformation are attributable to similar scale variations in the character of the salt-caprock interface. These variations do not appear directly related to shear zones or faults.

58 GEOSCIENCES↗

DOE OSTI2023

Synthetic Calcium Carbonate Production by Carbon Dioxide (CO 2 ) Mineralization of Industrial Waste Brines

The global scale of CO 2 emission has crossed 36 B tons, and the United State represents 14% of the total emissions. In light of the high cost associated with current CO 2 capture processes ($60-to-150 per ton of CO 2 ), carbon capture and utilization (CCU), wherein CO 2 is converted to beneficial products, provides a pragmatic path to overcome the economic barrier for CO 2 emissions control. In particular, CO 2 mineralization offers an attractive route as it produces high-value mineral carbonates that sequester CO 2 in a stable form. For instance, fine precipitated calcium carbonate, a valuable product with unit price in the range of $230-280/t and a global market projected to reach 99 M tons in 2020, can be produced by capturing CO 2 within aqueous Ca 2+ solution. However, two critical challenges need to be overcome—(i) the need for costly processes such as electrolysis or addition of alkali hydroxides to maintain alkalinity during mineralization, and (ii) the geographic availability of Ca-rich solutions (or brines) which often renders efficient integration with power plants impractical. To simultaneously address these challenges, this project developed two novel routes to enable carbonate production based on efficient CO 2 mineralization. In one variation, the project developed a new alkaline carbonation process to capture CO 2 and produce calcite precipitates from coal ashes. Herein, coal ash is carbonated first in a sodium carbonate solution. The carbonation reaction produces sodium hydroxide and raises the solution pH. This high-pH hydroxide solution is then used for CO 2 capture, which converts the sodium hydroxide solution back to sodium carbonate to repeat the carbonation cycle. The carbonated ash residue is refined with a CO 2 pressure swing step to produce high-purity precipitated calcium carbonate. In another variation, Ca-rich produced water brine serves as the Ca-source to mineralize CO 2 . An H + /Na + ion-exchange cycle was designed to provision alkalinity during mineralization and regenerate the ion-exchange reagent in brines with high salinity. Ca-depleted brines can then be treated within a centralized water treatment facility. Both variations beneficially utilize reject streams—such as coal ashes and brines from oil and gas extraction or CO 2 storage operations—that are available at substantial quantities in the vicinity of coal power plants within the U.S. In addition to technology development, this project included techno-economic and life-cycle analysis to identify technically and economically appropriate solutions for power plants in different geographic locations. Taken together, the project provides a unique route to integrate CO 2 emissions control and waste handling/treatment for coal power plants, while producing a high-value product to offset the economic burden associated with waste management. The developed processes utilize post-desulfurization flue gas from coal-fired power plants as is. The CO 2 conversion reactions are performed in alkaline brines at ambient pressure from the flue gas, thereby minimizing the energy burden. By beneficiation of industrial waste streams with low energy input, the processes offers significant technical advantages in energy and CO 2 footprint over the current paradigm of precipitated calcium carbonate production, which involves calcination of limestone at processing temperatures in excess of 800 °C. Overall, this project developed viable CO 2 mineralization processes while strategically maximizing the economic benefit.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Investigation of oxygen ion mobility through Zr x Ta 2 O 2x+5 (ZTOx) and Hf 6 Ta 2 O 17 (HTO) at elevated temperatures

Thermal Protection Systems (TPS) allow for technology to withstand high heat, which is especially critical for aerospace applications. As technology advances, better TPS are required to keep up with more extreme operating temperatures. Ultra-High Temperature Ceramics (UHTCs) are of interest because of their high melting points above 3000°C, but they suffer from oxidation when exposed to the atmosphere. A standard remedy to prevent oxidation is the use an oxide coating to shield from oxygen in the atmosphere while still allowing high operating temperatures and thermal cycling. However, standard oxide coatings such as yttria stabilized zirconia (YSZ) have limited operating temperatures before they allow oxygen to diffuse into the substrate material. A recently proposed alternative material has the formula A 6 B 2 O 17 (A: Zr or Hf) (B: Ta or Nb). In this study, determination of the oxygen protection ability of Hf 6 Ta 2 O 17 (HTO), Zr 6 Ta 2 O 17 (ZTO6), and various Zr x Ta 2 O 2x+5 (ZTOx) was determined using Electrochemical Impedance Spectroscopy (EIS). Optimal synthesis and sintering conditions were determined for ZTO6 and several of its ZTO variations. Activation energy for oxygen mobility was also determined. The importance of ball milling was found to be crucial to achieve high sintering density with these ceramics. The sintering study revealed that the ideal temperature of the ZTO variations increased with increased Zr/Ta ratio in the formulation. It was determined that both HTO and the ZTO variations conduct oxygen ions less than current ceramic coatings such as YSZ, with HTO conducting the least, with results confirming the validity of both HTO and ZTO materials as candidates for TPS coatings based on their low oxygen ion conductivity.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗

DOE OSTI2023

The Energy and Operational Impacts of Using 0-10V Control for LED Streetlights

Lighting systems have historically utilized either a proprietary control method, or one of a handful of standardized methods (e.g., 0-10V, DALI, DMX512) aimed at facilitating vendor interoperability. The utility and market success of the standardized methods has been limited for a variety of reasons. Although 0-10V is presently the most commonly available control interface for indoor and outdoor products in North America – even after the emergence of LED products and “connected lighting systems” with more modern network interfaces and luminaire-level sensors and intelligence – the use of 0-10V methods has significant tradeoffs. For instance, it is difficult to predict relative luminaire light output and input power at any particular control voltage, and performance across luminaires is inconsistent. This study characterizes 23 LED streetlights that claimed dimmability via a 0-10V interface, quantifies the performance variation found in market-available LED drivers, and explores the potential impact of the most recent 0-10V voluntary standard (ANSI C137.1-2022). Notably, the tested products were all manufactured prior to the release of this standard, and thus do not make compliance claims. Variation in response to 0-10V control voltages was expected to be more significant across different make/model luminaires than different units of the same make/model; the test population was structured based on this expectation. The 23 tested luminaires comprised 21 unique make/model streetlights from 14 different manufacturers. Three samples of one make/model streetlight were included in the test population to facilitate a limited exploration of unit-to-unit variation. The 21 unique make/model luminaires contained 20 unique make/model LED drivers from eight different driver manufacturers. Luminaire input power and current were measured at 11 control voltages (0.5 V and 1–10 V in 1 V increments).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Rapid AI-based Dissection of Ashes using Raman and XRF Spectroscopy (RADAR-X)

Waste-to-energy (WTE) facilities incinerate ~35 million tons of municipal solid waste annually in the United States. The incineration process reduces the mass and the volume of the waste fraction by over 75 and 95%, respectively. The fraction left after incineration remains as ash residues and is referred to as WTE ash, compromising of bottom and fly ash. In the United States, ~10 million tons of WTE ashes are generated annually and predominantly landfilled because of the lack of secondary end-use pathways. This incurs a significant financial burden (landfilling costs) on U.S. WTE facilities and also results in the loss of materials to the landfill. The primary objective of this research is to understand better the elemental and mineralogical composition of WTE ashes from diverse origins and find composition dependent upcycling pathways for diverting these ashes from landfills. This primary objective was addressed through three research tasks –(Task I) An AI-led Multi-Modal Approach for Compositional Analysis, (Task II) Developing a Dissolution-Based Test for Real Time Analysis, and (Task III) Establishing composition-dependent end uses. The chemical composition of WTE ash is dependent on two factors, i.e., the input waste composition and the operational parameters of a WTE facility (combustion conditions). Amongst these two factors, the input waste composition will likely show spatial and temporal variation. As a result, the chemical composition of WTE ash will also fluctuate. To understand the spatial and temporal variation in WTE ash composition, in Task I, we collected 128 ash samples (62 bottom ash and 66 fly ash samples) from 11 WTE facilities located in 11 U.S. states and characterized them via X-ray Fluorescence, powder X-ray Diffraction, and Raman Spectroscopy. The findings from this extensive characterization work indicated that the key elements in WTE fly ashes are Ca, Cl (greater than 10 wt. %), Si, S, K, Zn ( between 1 and 10 wt. %), Mg, Al, P, Ti, Fe, Cu, Br, and Pb (between 0.1 and 1 wt. %). Similarly, the key elements in WTE bottom ash fraction finer than 45μm are Ca (greater than 10 wt. %), Mg, Al, Si, S, Cl, K, Ti, Fe, Zn (between 1 and 10 wt. %), P, V, Cr, Mn, Cu, Br, and Pb (between 0.1 and 1 wt. %). Here, we note that the dominant fraction of WTE bottom ash is the coarse fraction. The coarse WTE bottom ash fraction (rich in silicon) was not characterized in this study because of excessive grinding requirements and their unsuitability as a supplementary cementitious material due to their coarse nature. The elements in WTE bottom ashes are present as calcite, anhydrite, vaterite, hydroxyapatite, quartz, bassanite, gehlenite, akermanite, hydrocalumite, and portlandite. Similarly, the mineralogical species present in WTE fly ashes are calcium chloride hydroxide, halite, calcite, anhydrite, sylvite, hydrocalumite, vaterite, hannebachite, quartz, and bassanite. Temporal variation in ash composition may also result in significant fluctuations in chemical compositions. Therefore, a WTE facility may need to monitor the ash composition (elemental and mineralogical composition) in real time. In Task I, we evaluated the possibility of using a portable X-ray fluorescence (XRF) spectrometer to monitor the elemental composition in real-time. Specifically, we collected XRF data on identical specimens via a portable XRF spectrometer (low-end) and a lab-based benchtop XRF spectrometer (high-end). The collected data was used to train an A.I. algorithm (portable XRF data as an input and benchtop XRF data as an output) to predict accurate elemental composition using portable XRF data. Finally, we developed a 2-minute photobleaching protocol to monitor the mineralogical characteristics of WTE ashes via Raman spectroscopy. Overall, the activities in Task I improved our understanding of ash composition and developed techniques to monitor elemental and mineralogical composition in near real-time. Based on the findings of Task I, we find that WTE ashes exhibit wide variability in mineralogy. For ICP-based elemental analysis, all the mineralogical species in WTE ashes must be brought into solution. This is traditionally accomplished with acid digestion using a combination of multiple acids. However, acid digestion with multiple acids is time-consuming and often fails to ensure complete digestion of the ash matrix. To address this limitation, in Task II, we developed an alkali-fusion-based digestion protocol using a combination of lithium tetraborate, lithium metaborate, and their combinations as possible alkali fluxes for digesting WTE ashes entirely and rapidly. The validity of the developed method was evaluated on two standard ash specimens, i.e., SRM 1633c coal fly ash and BCR-176R incineration fly ash specimen. The findings suggest that the developed protocol can ensure complete digestion of elements such as Al, Ba, Ca, Cr, Cu, Mg, Mn, P, Sr, V, Zn, Be, K, and rare earth elements. The recent changes in the energy market towards renewables and increased metal recycling have resulted in reduced supplies of supplementary cementitious materials (coal fly ash and slag). Therefore, in Task III, we evaluated the possibility of employing WTE ashes as SCMs. As the chemical composition of WTE ashes varies temporally (on an hourly basis), there was also a need to develop tests that can evaluate the suitability of material to act as supplementary cementitious material rapidly, i.e., in a few minutes. Therefore, in Task II, we also developed a rapid test to assess the suitability of a material to act as an SCM in 5 minutes. This represents a significant advance over the state-of-the-art R 3 test, which takes ~144 hours. This test was initially validated on amorphous aluminosilicates, such as calcined clays, and could be extended to evaluate other industrial by-products, such as WTE ashes. In Task III, we evaluated the possibility of employing WTE ashes for two applications, i.e., as an SCM and a lime substitute for clay stabilization. The findings from Task I indicated that WTE ashes are enriched in chlorine and, therefore, cannot be used directly as an SCM due to corrosion-related risks and altered hydration kinetics. Accordingly, we developed an ash treatment protocol to reduce the solubility of chlorine-containing species in WTE ashes. The developed treatment protocol also immobilized lead in certain mineral forms. As a result of the treatment, WTE ashes can be used as SCMs without any corrosion or heavy metal leaching concerns. The second application examined in this study was clay stabilization. WTE ashes are calcium-rich and can be an adequate lime replacement for clay stabilization. Our findings reveal that the sum of the concentrations of Ca(OH) 2 and CaClOH controls the clay stabilization capability of WTE ashes. In summary, in this work, we evaluated the elemental and mineralogical characteristics of U.S. WTE ashes from diverse origins and developed tests to evaluate the chemical characteristics of these ashes in real time through a portable XRF and a benchtop Raman spectrometer. Based on the chemical characteristics of these ashes, we developed an ash treatment process to enable the use of WTE ashes as an SCM and also evaluated the possibility of employing these ashes for clay stabilization. Overall, the findings from this work enables the diversion of WTE ashes from landfills for multiple end-uses, i.e., as an SCM or a lime substitute for clay stabilization.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗

DOE OSTI2024

2024 Update of Comprehensive Review of Multi-arm Caliper Data for the Big Hill SPR Site

58 GEOSCIENCES↗