Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic load model”

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.

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

At least 991 records · Page 55

Influence of Pore Water on the Fracture of Silica Sand at Particle Scale

The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.

cavitation↗

Modeling of ATR Fuel in DOE Standard Canisters with Helium Backfill

One pathway for road-ready and final disposition packaging configurations for the aluminum-clad spent nuclear fuel (ANSF) fuel is storage within helium backfilled sealed Department of Energy (DOE) standard canisters. The typical packaging configuration for the 15-foot DOE standard canisters places 10 advance test reactor (ATR) elements a basket, and three baskets are loaded within each DOE canister. During in-reactor operations and cooling pond storage conditions, oxyhydroxide layers form on the surface of the aluminum clad fuel. These layers produce hydrogen gas over time due to the fuel’s radiation field. As part of the packing procedure, the ATR fuel should be dried to remove any residual physio-/chemi- sorbed water bound to the surface. A 50-year CFD model of the DOE canister packaged with fuel was developed to provide a temperature profile for coupled chemical modeling of the conditions within the canister. The results of this modeling include results at fully saturated and fully dried fuel cladding conditions. In the associated experimental work, radiolysis experiments tests were completed in a helium environment, and G-values for the radiolytic production of hydrogen from the oxyhydroxide layers were provided. That reaction was coupled with the thermal profiles and gas-phase reactions to develop a 50-year model of the conditions within a sealed DOE canister with ATR fuel. For a nominal scenario of stored ATR fuel, after 50 years the model results give a 1.36 atm total pressure, 7% mole percent hydrogen, for the upper decay heat, 1.51 atm total pressure, 16% mole percent hydrogen, and for upper decay heat with undried fuel 2.6 atm total pressure, 15% mole percent hydrogen. No case modeled yields significant oxygen, and for the lower decay heat case that is modeled, hydrogen concentrations are under the 4% flammability limit after 50 years of storage. The modeled pressures for all cases modeled are below the pressure limit for the DOE standard sealed canister.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Design and optimization of flexible decoupled high-temperature gas-cooled reactor plants with thermal energy storage

Advanced nuclear power plants are well-positioned for future zero-carbon grids, however, the need for flexible power generation will be required over the traditional emphasis on baseload generation for meeting historical demands. To achieve such flexibility, this work examines viable configurations for coupling nuclear energy production with thermal energy storage. Previous designs on nuclear-thermal energy storage configurations for advanced reactor designs, which utilized reactor steam as the heat source for charging the thermal energy storage, are restricted by the heat diversion ratio and efficiency losses, thus their impacts can be limited. In this context, this study proposes configurations for fully decoupling the nuclear reactor from the power cycle and positioning the storage as an intermediate loop, thereby achieving an unconstrained heat diversion ratio and improved efficiency. Compared with a standard high-temperature gas-cooled reactor’s power cycle, steady-state thermodynamic modeling and dispatch optimizations quantify the benefits of a steam reheat cycle within the fully decoupled thermal energy system to separate the plant cycle from the high-pressure primary side. These benefits are further detailed, compatible with required high-temperature and high-pressure conditions, through (1) open-source dynamic transient models that examine the impact of off-design operation on the systems, (2) the investigation of components design and costing and finally (3) sizing and dispatch optimization. The fully-decoupled design achieves a cycle efficiency of 43.1%, an enhancement over the vendor’s standard efficiency of 42.2% (Xe-100 design). Here, the proposed design offers strengthened physical barriers from the nuclear island as well as superior operational flexibility and power boosting. Dispatch optimization and market analysis reveal that thermal energy storage size is highly dependent on the peak patterns of electricity prices and the minimum generation level constraint imposed on the balance of plant. Evaluation of off-design operation demonstrates that the full decoupling design with the suggested fail-safe control mechanisms ensures a minimal impact on reactor parameters, even during rapid power ramping.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Wildfire Risk Evaluation Framework for Grid Operations and Planning

The United States (US) environmental protection agency's (EPA's) climate change indicators for wildfires show a long-term trend of increased annual wildfire activity, larger wildfire size, and more variable dynamics in wildfire behavior. This has caused more frequent preemptive public safety power shutoff (PSPS) events in the regions with recognized high wildfire risk. These preemptive power shutoffs attempt to prevent the ignition of a wildfire but it nonetheless renders the transmission line non-operational, which often sheds load of downstream communities. As wildfires and PSPS events become more frequent, it is crucial to find the communities most at risk of load shedding. To that end, this paper proposes the wildfire risk evaluation of the system (WiRES) framework, which is a performance-based framework that translates extreme weather-related and PSPS event probabilities into a cumulative probability of a non-operational transmission line. This study also provides a geospatial visualization tool that breaks down the entire western electrical coordinating council (WECC) region into 50 km grid cells which can be used to 1) filter out transmission lines with higher than a threshold outage probability, and 2) graphically discover the affected regions and their biophysical and socioeconomic metrics. %such as the social vulnerability index, population density, gross domestic product, etc. Lastly, an impact assessment study is conducted which connects the results of the proposed framework to python and powerworld-based contingency analysis to highlight the applicability of the framework.

Wildfire, Resilience, Extreme Events, Grid Reslien↗

Grid Resiliency with a 100% Renewable Microgrid

San Diego Gas & Electric Company (SDG&E) installed America’s first and largest utility-scale microgrid in Borrego Springs in 2013. The first generation Borrego Springs Microgrid utilized diesel generators to form and stabilize the microgrid island, with support from grid-scale batteries and local solar photovoltaic (PV) generation. In this project, SDG&E in partnership with National Renewable Energy Laboratory (NREL) demonstrated through modeling, simulation and utility field testing that blackstart and islanding of the microgrid can be led with 100% renewable, inverter based resources (IBRs), to help reduce community reliance on conventional generation resources. Through equipment upgrades, grid-forming island leader capability was transitioned to a battery IBR instead of the Borrego Springs Microgrid diesel generators. A new microgrid controller was integrated to the microgrid and programmed to control and manage multiple energy storage systems. Synchrophasor and other power quality data verified autonomous, high-speed response of the IBRs through blackstart, islanding, and load step testing. Results of project field evaluations provide distribution systems operators (DSO) with increased confidence that renewable, IBR can replace traditional generators to blackstart and island microgrids and rapidly establish stable island frequency with rapid changes in peak power demand. Importantly, the project validated the integration feasibility of a distributed energy resource management system (DERMS) controller that manages multiple grid-forming and grid-following IBRs, establishing a standard design interface to reduce the complexity of integrating new DERs in the future and supporting replication by the industry. As a result of learnings in this project, SDG&E has implemented the microgrid controller strategy at multiple other microgrid sites, thereby validating the replicability of the solution. Hardware-in-the-loop (HIL) simulations including power and controller HIL hardware — along with electromagnetic transient (EMT) simulations of Borrego Springs Microgrid —informed adjustments to inverter parameters and were important to characterize the performance of the IBRs in relevant operating conditions before deployment. The EMT and HIL simulations of islanding the entire community are important contributions in providing confidence in IBR performance prior to future islanding of the community in the field. High-fidelity EMT and/or HIL simulation of IBRs can de-risk field operations, and its relevance and importance as a tool is increasing as distribution grids and microgrids become more complex and dynamic with an increasing proportion of renewable generation, distributed energy storage, and two-way power and energy flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Understanding Transient Combustion Phenomena in Low-NO x Gas Turbines

The goal of this work was to understand, quantify, and eventually predict combustion instability during operational transients in low-NO x gas turbine combustors. Transient operation is defined as a mode where one or more key operating parameters, including equivalence ratio, fuel composition, or fuel splitting between nozzles in a combustor, are varying in time. The timescale over which this variation occurs will be referred to as the transient timescale, the magnitude of the change in the transient parameter will be referred to as the transient amplitude, and the direction of the variation will be referred to as the transient direction. Transients occur regularly in power generation gas turbine operation. For example, load-following gas turbines vary load to meet power demands while accounting for variation in unsteady grid contributors like wind and solar power. Operators achieve this by varying equivalence ratio and fuel splitting inside the combustors to produce more or less power. Another example of transient operation is during fuel-switching operations, a particular issue for turbines operating with syngas fuels in a co-generation facility. In this case, the composition of the fuel stream may vary in time either as a result of fluctuations in the pipeline, or in the case of a co-generation plant, variations in feedstock or gasification processes. This fluctuation in fuel composition can have significant implications for combustion processes, particularly with high-hydrogen fuels. A small number of studies have investigated transient behavior in gas turbine combustors, particularly in the context of flame blow-off; here, investigators have typically focused on issues of transient amplitude, though without the depth and physics-based experimental design that is proposed in this program. To our knowledge, there have been no studies that have quantified the impact of the transient timescale on combustion. The work helped us understand the impact of several types of transients on combustion stability. Both static and dynamic combustion stability will be considered. The transients will be quantified using their amplitude, direction, and timescale, and we will consider transients in equivalence ratio, fuel composition, and fuel splitting. There were two main objectives of this work: 1. Provide a fundamental understanding of combustion stability during transient operation using data obtained from a number of high-speed diagnostic techniques on experimental configurations that capture important gas turbine combustor features and operating conditions. This goal was met in conjunction with our industrial partner, GE Global Research, with whom we worked closely throughout the program. 2. Develop a stability prediction framework that will allow for the prediction, and possible control, of instabilities as a function of the transient amplitude, direction, and timescale. This input/output model should be physics-based and allow for screening of operating conditions and prediction of unstable behavior at the end of a transient by monitoring combustor health and precursor signals during operation.

03 NATURAL GAS↗

High-Performance Integrated Window and Façade Solutions for California

The researchers developed a new generation of high-performance façade systems and supporting design and management tools to support industry in meeting California’s greenhouse gas reduction targets, reduce energy consumption, and enable an adaptable response to minimize real-time demands on the electricity grid. The project resulted in five outcomes: (1) The research team developed an R-5, 1-inch thick, triplepane, insulating glass unit with a novel low-conductance aluminum frame. This technology can help significantly reduce residential cooling and heating loads, particularly during the evening. (2) The team developed a prototype of a windowintegrated local ventilation and energy recovery device that provides clean, dry fresh air through the façade with minimal energy requirements. (3) A daylight-redirecting louver system was prototyped to redirect sunlight 15–40 feet from the window. Simulations estimated that lighting energy use could be reduced by 35–54 percent without glare. (4) A control system incorporating physics-based equations and a mathematical solver was prototyped and field tested to demonstrate feasibility. Simulations estimated that total electricity costs could be reduced by 9-28 percent on sunny summer days through adaptive control of operable shading and daylighting components and the thermostat compared to state-of-the-art automatic façade controls in commercial building perimeter zones. (5) Supporting models and tools needed by industry for technology R&D and market transformation activities were validated. Attaining California’s clean energy goals require making a fundamental shift from today’s ad-hoc assemblages of static components to turnkey, intelligent, responsive, integrated building façade systems. These systems offered significant reductions in energy use, peak demand, and operating cost in California.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

Grid-Connected Modular Soft-Switching Solid State Transformers (M-S4T)

The objective of this project is to develop and verify the concept of a flexible and modular soft-switching solid-state transformer (M-S4T) for direct grid-connected applications. The ability to directly connect power electronics converters to the medium voltage grid (4 kV – 13 kV), and to potentially replace the passive and bulky, but ubiquitous 60 hertz service transformer in the 25 kVA to 100 kVA range, with a more flexible and controllable device, has been regarded as the ‘holy grail’ in grid control. However, this has proven to be extremely difficult. This project has developed the solutions to several key challenges of the direct grid-connected power electronics and realized a 7.2 kV M-S4T prototype. First, a protection method to protect the M-S4T from the high voltages (110 kV for the 13 kV system) that occur on the grid due to transients and lightning strikes have been developed and experimentally verified. Second, the realization and the operation of the M-S4T based on high-voltage SiC devices (>3.3 kV) and a medium-frequency medium-voltage low-leakage transformer in a single-stage solid-state transformer with zero-voltage switching, low dv/dt, and low electromagnetic interference has been successfully demonstrated up to 7.5 kV peak. Third, an oil-cooling system and stable communication and distributed control system for converter module voltage sharing have been developed and experimentally verified. The developed M-S4T has realized a modular universal high-performance power conversion system. This conversion system is scalable to different voltage and power levels and adaptable to four-quadrant bidirectional operation. Moreover, the use of passive cooling techniques meets the equipment life requirements, and the lightning protection scheme fulfills the basic insulation level specifications for direct grid connection. Such power conversion system opens up near-term opportunities, including energy storage, solar PV, or electric vehicle charging with significant cost and footprint savings. In the longer term, the possibility of replacing the utility distribution transformer with an M-S4T will be transformative for future distribution grids with a compact footprint and full controllability to enable high renewable energy and storage penetration. In addition to the main project, this report expands on the Plus-Up projected including as part of the main award. This project developed and demonstrated the technology for autonomous collaborative inverters that can be connected in an ad hoc manner to the grid. The aim of the project was to: (1) evaluate the existing techniques for grid-connected inverters and find their limitations; (2) develop detailed requirements for grid-connected inverters in the modern grid with millions of active nodes; (3) design a unified control strategy that brings more autonomy and intelligence to grid-connected inverters, and addresses parts of the issues with the existing techniques. The proposed technique, called UniCon, enables inverters to 1) connect/disconnect to/from the grid in an ad hoc manner; (2) work based on local sensing. Slow communication could be used for a more optimized behavior; (3) work automatically in both grid-forming/grid-following mode; (4) handle large disturbances, e.g., big load step and fault, in an oscillation-free manner; (5) work collaboratively with other inverters in steady-state and during transients. UniCon can be implemented in the middle-level control; hence it is agnostic to the vendor and to the implementation of the inner voltage/current and protection loops. Furthermore, a new synchronization scheme, based on deep learning, was developed that can extract the grid voltage phase and amplitude in a stable manner. The method is cheap to implement can improve the dynamic performance of the grid-connected inverters during fast transients, e.g., fault. The proposed control scheme was validated by (1) MATLAB/Simulink; (2) hardware-in-the-loop results, and; (3) experimental results using three inverters that form a microgrid in a down-scaled feeder. Lastly, both the M-S4T and UniCon have achieved promising tangible paths to markets. In the case of the M-S4T, the underlying technology — the Soft Switching Solid State Transformer (S4T) developed at the Georgia Tech Center for Distributed Energy (GT-CDE) has been licensed by GridBlock from the Georgia Tech Research Corporation, and GridBlock has been working with manufacturing partner Jabil (one of the largest US-based contract manufacturers) and system integrator Power Secure (largest deployer of microgrids in the US with 4.7 GW under management), to meet the strong initial demand. Similarly, GridBlock has an exclusive license to the UniCon technology, developed under this award by GT-CDE. The UniCon provides an intermediate control layer that enables the implementation of the higher-level ‘transactive’ control commands for the system. The architecture of the system - slow communications with the cloud for system optimization and setpoints, and the use of locally measured quantities for real-time control, provide a very robust and secure way of implementing a real-time must-run grid that is also secure and stable. This is a brand-new functionality that is critical for the future grid and key to GridBlock’s business model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗

TRANSENSOR: Transformer Real-time Assessment INtelligent System with Embedded Network of Sensors and Optical Readout. Final Report

Utilities across the world are wrestling with evolving market dynamics, population growth and climate change. Distributed Energy Resources (DERs) such as solar photovoltaics, distributed generators and energy storage systems are becoming important parts of the U.S. energy mix. These factors are driving an increasing need for low-cost grid asset monitoring. Aside from being costly, traditional utility monitoring systems are not sufficiently robust and do not provide real-time visibility into the condition of grid assets such as transformers. Lack of accurate real-time measurements on performance has resulted in the use of lagging indicators, such as oil sample analysis. To address this need, an innovative embedded optical sensing technology, Transformer Real-time Assessment Intelligent System (TRANSENSOR) was developed, validated and demonstrated in this project. To date, TRANSENSOR has focused on transformers but is extendable to other grid assets. In addition to the technology being embedded into new transformers during manufacturing, a retrofit configuration that can be installed on existing transformers in the field was also developed. It is anticipated that the ability to retrofit existing transformers will help accelerate adoption of the underlying TRANSENSOR technology. Phase 1 of the project focused on laboratory development and qualification of the technology. Following iterations and exploration of relevant optical sensing modalities and multiplexed configurations of interest for the transformer environment, an effective candidate configuration was agreed upon, down-selected and custom-designed for embedding into General Electric (GE) network transformers. Two new GE 500 kilovolt ampere (kVA), 27 kilovolt (kV) distribution network transformers were built with embedded fiber-optic (FO) sensors and successfully qualified per industry standards at GE’s Shreveport, Louisiana facility during Phase 1 of the project. Following the successful completion of Phase 1, the team proceeded to Phase 2, which focused on a field demonstration of the technology. Over Phase 2, the first new GE transformer built with embedded fiber-optic sensors was installed in an above-ground cage and connected to the grid at Con Edison’s Astoria facility. A second transformer equipped with TRANSENSOR was installed in an underground vault. Additionally, an older (1982 year model) GE transformer in an above-ground cage was retrofitted with fiber- optic sensors and reconnected to the grid at ConEd’s facility. Analysis of data acquired from the sensors showed interesting correlations with transformer loading. Additionally, key events such as transformer low-voltage network connection, primary-side energizing, and a pressure loss event from an oil sampling were detected by the TRANSENSOR system. Online data processing/feature extraction algorithms for the second transformer in the underground vault detected key features and event alerts that were transmitted through a wireless 4G connection. The remote deployment concept showed promising results with data collection running for a total of 8 months across the three (3) GE transformers instrumented for the Phase 2 field demonstration. This sets the stage well for further development and commercialization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on the Banshee Distribution Network: Preprint

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on Banshee Distribution Network

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling of ATR fuel in DOE Standard Canisters with Helium Backfilled Condition

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel dictates storage within helium backfilled sealed DOE standard canisters. These sealed canisters are intended for extened (>50 year) dry storage). The typical packaging configuration for the 15-foot DOE canisters places 10 ATR elements within a Type 1a basket, and three baskets are loaded within each DOE canister. During in-reactor operations and cooling pond storage conditions, oxyhydroxide layers form on the surface of the aluminum clad fuel. These layers produce hydrogen gas over time due to the fuel’s radiation field. As part of the packing procedure, the ATR fuel should be dried to remove any residual physio-/chemi- sorbed water from the surface. As testing to the effectiveness of the drying procedure is still underway, this modeling will include results at fully saturated and fully dried conditions. In previous modeling efforts, the G-value for the production of hydrogen from the oxyhydroxide layers was assumed to be in argon environments as measured by Task 2 - Oxyhydroxide Layer Radiolytic Gas Generation Resolution. Previous experimental testing showed differences in the hydrogen generated based on the gaseous environment. In the associated experimental work, Task 2 - Oxyhydroxide Layer Radiolytic Gas Generation Resolution, additional tests were completed in a helium environment, and updated G-values for the radiolytic production of hydrogen from the oxyhydroxide layers were provided. These values are 28% and 58% higher than values for argon. In addition, a change to the modeling of the oxyhydroxide radiolysis has been made from previous reports. This change assumes the dependency of the dose rate on the overall reaction rate is applied to the total weight of the sample, rather than to just the weight of the oxyhydroxide layer. This decreases the dependency of the radiolytic reaction on the thickness of the oxyhydroxide layer. For a nominal scenario of stored ATR fuel, assuming the chemi-/physio- sorbed water have been fully removed, the internal canister pressure increases to 1.61 atm over a 50 year period, with a hydrogen mole percentage of 21%. As in previous modeling, any oxygen present is in negligible amounts (<1 ppt). If a small amount residual air is present, nitric acid can form up to 1300 ppm. For a scenario with high fuel decay heat, the model shows internal pressure increasing to 2.1 atm, with 39.3 mole percentage of hydrogen. In a scenario where significant chemi-/physio- sorbed water is present within the corrosion layer, the nominal scenario shows a pressure increase to 2.54 atm, with 21.1 mole percent hydrogen. The high decay heat case shows a pressure increase to 3.18 atm with 39.9 mole percent hydrogen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ARPA-E Grid Optimization (GO) Competition Challenge 2

The ARPA-E Grid Optimization (GO) Competition Challenge 2, from 2020 to 2021, expanded upon the problem posed in Challenge 1 by adding adjustable transformer tap ratios, phase shifting transformers, switchable shunts, price-responsive demand, ramp rate constrained generators and loads, and fast-start unit commitment. Furthermore, Challenge 2 was a maximization problem while Challenge 1 was a minimization problem. Specifically, the economic surplus, defined as the benefit of serving load minus the cost of generation, is being maximized. It was expected that the objective value of a given solution should be positive, representing economic gain, but negative objectives from poor solutions were possible. The two code submission feature of Challenge 1 was maintained. Additionally, Divisions 3 and 4 within the competition permitted on/off switching of transmission lines (Divisions 1 and 2 did not). After the initial release of the Problem Formulation on 7/20/2020, ARPA-E Director Lane Genatowski announced Challenge 2 on 9/12/2020. The final May 31, 2021, version of the Problem Formulation was 97 pages long with 299 equations. The Challenge proceeded with 2 non-prize Events and 2 prize Events. Teams receiving Challenge 1 FOA awards and prize money were required to use the prize money to fund their Challenge 2 efforts (Georgia Institute of Technology, Global Optimal Technology, Inc., Lawrence Livermore National Laboratory, Lehigh University, Northwestern University, Artelys, Columbia, Pearl Street Technologies, Pennsylvania State University, and University of Colorado Boulder). For more information on the competition and challenge 2 see the "GO Competition Challenge 2 Information" resource below. Challenge 1 and Challenge 3 information can be found in the resources linked below.

ACOPF↗