Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Microgrid design and multi-year dispatch optimization under climate-informed load and renewable resource uncertainty

Microgrids are an increasingly popular solution to provide energy resilience in response to increasing grid dependency and the growing impacts of climate change on grid operations. However, existing microgrid models do not currently consider the uncertain and long-term impacts of climate change when determining a set of design and operational decisions to minimize long-term costs or meet a resilience threshold. In this paper, we develop a novel scenario generation method that accounts for the uncertain effects of (i) climate change on variable renewable energy availability, (ii) extreme heat events on site load, and (iii) population and electrification trends on load growth. Additionally, we develop a two-stage stochastic programming extension of an existing microgrid design and dispatch optimization model to obtain uncertainty-informed and climate-resilient energy system decisions that minimizes long-term costs. Use of sample average approximation to validate our two case studies illustrates that the proposed methodology produces high-quality solutions that add resilience to systems with existing backup generation while reducing expected long-term costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Catalyzed oxidation of nuclear graphite by simulated fission products Sr, Eu, and I

The influence of three fission products Sr, Eu, and I on the oxidation of IG-110 nuclear graphite was studied in the temperature range of 400 to 1000 °C. Sr and Eu were introduced as chlorides, and I was introduced as NaI. The temperature dependence of both CO 2 and CO production during the graphite oxidation measured with mass spectroscopy and infrared spectrometry shows that the introduction of these three compounds to graphite significantly decreases the onset temperature for the oxidation of graphite. Among the three compounds, NaI is the most active towards the oxidation reaction, characterized by a significant decrease of the onset temperature from approximately 650 to 400 °C before and after its introduction to graphite. Separate measurements of CO 2 and CO concentration at varying temperatures enable the calculation of the activation energy for the formation of CO 2 and CO. The activation energies for the oxidation of pure and fission product-impregnated graphite samples decrease in the following order: standard IG-110 graphite, EuCl 3 -impregnated IG-110, SrCl 2 -impregnated IG-110, and NaI-impregnated IG-110. This trend indicates that the three compounds catalyze the oxidation of graphite at temperatures relevant to the operation of high-temperature gas-cooled reactors. Furthermore, it is found that the three compounds can also affect the molar ratio of reaction products CO 2 and CO, and the rates of the graphite oxidation. At temperatures higher than about 850 °C, the impregnated samples exhibit lower CO 2 : CO ratios than the pure graphite. Different from EuCl 3 and NaI, the introduction of SrCl 2 decreases the graphite oxidation rates at temperatures higher than about 770 °C. Their catalytic mechanism can be understood based on a redox cycle of the intermediate active species, promoting the dissociation of molecular oxygen and transfer to the carbon.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A Functional Reference Architecture and Assessment Thereof for the National Ignition Facility Industrial Controls Systems

The Industrial Control System (ICS) at the National Ignition Facility (NIF) has an effective, well-established architecture based off a conventional cyclical software paradigm focused on maintainability and the identification of issues. This architecture exhibits scalability in ensuring expansions of the system conform to the existing architecture, modularity enough to allow relatively easy integrations of such expansions and works as a successful tool to introduce control system engineers new to the NIF ICS to the structure of the system at each layer. This architecture, like most software architectures, is object-oriented, lending itself to ease of understanding by control systems engineers and software engineers familiar with an object-oriented perspective. There are occasions, however, where engineers of other disciplines require insight into the functionality and structure of the ICS for the purposes of understanding fundamentally how their own system is or will be governed by the ICS, without the need for the details of operation of the ICS or the object-oriented view. For this reason, a functional architecture of the ICS could be a potent tool for communicating this insight. Even more powerful, a generalization of this proposed functional ICS architecture in the form of a National Ignition Facility and Photon Science (NIF & PS) Industrial Controls Reference Architecture could communicate this insight not just to systems governed by the ICS in the NIF proper, but across entirety of the NIF & PS Principal Associate Directorate (PAD), anywhere an instance of the ICS architecture is present, such as the approximately 40 “small labs” distributed across the directorate. Such a tool will provide an alternative means of understanding the implementation of these control systems, conducive to a larger variety of engineering and scientific disciplines.

42 ENGINEERING↗
DOE OSTI2024

Creep Resistant Martensitic Steels for Operation at High-Temperatures in Power Generation Applications

Increasing the temperature capabilities of ferritic/martensitic 9-12% Cr steels can help in increasing the operating temperature of land-based turbines and minimize the use of expensive high-temperature alloys in the hot section. A creep resistant martensitic steel, JMP, was developed with the potential to operate at or above 650°C. The design of the alloys originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Cobalt was used for increased solid solution strengthening, Si for oxidation resistance and different W and Mo concentrations for matrix strength and stability. The JMP steels showed increases in creep life compared to MARBN/SAVE12AD at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN steel. Approximately 21 years of cumulative creep data are reported for the JMP steels which encompasses various compositions. The relationships between composition-microstructure-creep properties are discussed including characterization of microstructures after >20,000 hours in creep.

Detrois, Martin↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Computational investigation of the impact of metal–organic framework topology on hydrogen storage capacity

Metal–organic frameworks (MOFs) are promising, tunable materials for hydrogen storage. For application under cryogenic operating conditions, past work has run into a ceiling on performance due to a trade-off in the volumetric deliverable capacity (VDC) versus the gravimetric deliverable capacity (GDC). In this study, we computationally constructed and screened 105 230 MOF structures based on 529 nets to explore the effect of underlying topology on the hydrogen storage performance of the resulting materials. A machine learning model was developed based on simulated hydrogen uptake to facilitate screening of the entire dataset, and it successfully identified the top 10% of materials with a root-mean-square error of approximately 1 g L −1 as validated by subsequent grand canonical Monte Carlo simulations. We identified a promising structure based on the tsx topology that exhibits both VDC and GDC higher than the current benchmark material, MOF-5. Our data-driven analysis indicates that nets with higher net density yield MOFs with enhanced volumetric and gravimetric surface areas, thereby improving maximum VDC while shifting the capacity trade-off toward higher GDC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Universal control in bosonic systems with weak Kerr nonlinearities

Resonators with weak single-photon self-Kerr nonlinearities can theoretically be used to prepare Fock states in the presence of a loss much larger than their nonlinearities. Two necessary ingredients are large displacements and a two-photon (parametric) drive. Here, in this study, we find that these systems can be controlled to achieve any desired gate operation in a finite-dimensional subspace (whose dimensionality can be chosen at will). Moreover, we show that the two-photon driving requirement can be relaxed and that full controllability is achievable with only one-photon (linear) drives. We make use of both Trotter-Suzuki decompositions and gradient-based optimization to find control pulses for a desired gate, which reduces the computational overhead by using a small blockaded subspace. We also discuss the infidelity arising from input power limitations in realistic settings, as well as from corrections to the rotating-wave approximation. Our universal control protocol opens the possibility for quantum information processing using a wide range of lossy systems with weak nonlinearities.

Yuan, Ming [Univ. of Chicago, IL (United States)] ↗
DOE OSTI2022

Techno-Economic Analysis of a Cryogenic Flux Capacitor for Grid Storage

Abstract The Cryogenic Flux Capacitor (CFC) is a cold, dense fluid storage core with integrated design features that afford the designer flexibility and provide new possibilities for the storage and discharge of energy. The stored energy, in this case, is represented by hydrogen physically bonded within the nanoscale pores within the aerogel composite blanket material, and the process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (∼1,000 m2/g) allows for storage densities close to, or in some cases exceeding, that of normal boiling point liquids. Its performance easily exceeds what can be achieved via ambient temperature, high-pressure gas storage for an equivalent volume. CFC storage is predicted to be easily scalable, is constructed from readily available commercial materials, lends itself to a range of pressure applications, and is geometry insensitive. The high energy density of CFC-stored hydrogen allows long durations of storage, such as daily to monthly cycling, which corresponds to approximately 10 to 100 hours of duration, respectively. The techno-economic analysis examines different types of hydrogen storage for commercial viability. The assessment includes the use of hydrogen blended with natural gas in a combined cycle gas turbine. The study provides results for a range of blends from 100% natural gas to 100% hydrogen. The required fuel usages are estimated for the intended application in a grid scenario that features large amounts of renewable penetration, 10 hours, and as a baseload, 100 hours. The study presents preliminary estimates of the capital cost of storage and operating cost of storage. In addition, the cost of firing natural gas and capturing the carbon dioxide (CO2) in a carbon capture and sequestration (CCS) system is assessed. Based on the study, it is shown that CFC can competitively meet performance and cost needs for grid-based energy storage.

08 HYDROGEN↗
DOE OSTI2024

Predicting the Rate of Degradation Related to Oxygen Electrode Delamination in Solid Oxide-Ion Electrolyzers

One of the leading causes for the performance degradation in H 2 -producing solid oxide electrolytic cells (SOECs) is the gradual delamination of oxygen electrode (OE) from the electrolyte under a strong anodic polarization. Identification of the key factor that controls the rate of OE delamination is of paramount importance to achieve long-term stable operation of SOECs. Here we show from thousands of hours of testing data that the exchange current density (i°) of OE can be leveraged as a predictor for the rate of delamination. To obtain i°, we apply DC-biased electrochemical impedance spectroscopy on a three-electrode symmetrical cell to measure polarization resistance (R p ) of OE as a function of current density (i) and time (t). The collected R p -i-t raw data are then converted to overpotential ( η )-i-t, from which i° is extracted from the “low-field” approximation. An analytical relationship between i° and time-to-delamination (TTD) is further established from the established i°-i-t relationship. We show that under a constant anodic polarization current density i, the greater the ratio i/i°, the faster the delamination. Therefore, we conclude that the exchange current density of an OE, i°, can be used to predict the rate of OE degradation in solid oxide-ion electrolyzers.

25 ENERGY STORAGE↗
DOE OSTI2020

Cutting the Deployment Costs of Physics-Based MPC in Buildings by Simulation-Based Imitation Learning

It has been shown that model predictive control (MPC) is a promising solution for energy-efficient building operations. However, the deployment of MPC in a large portion of the building stock has not been possible partially because of high installation costs. Every building is unique and requires a tailored MPC solution. The best performing solutions are often based on physics-based modeling, which is, however, computationally expensive and requires dedicated software. A promising direction that tackles this problem is to train a neural network-based optimal control policy to imitate the behavior of physics-based MPC from the simulation data generated offline. The neural networks give control actions that closely approximate those produced by physics-based MPC, but with a fraction of the computational and memory requirements and without the need for licensed software. The main advantage of the proposed approach stems from simple evaluation at execution time, leading to low computational footprints and easy deployment on embedded HW platforms. In the case study, we present the energy savings potential of physics-based MPC applied to an office building in Belgium. We demonstrate how neural network approximators can be used to cut the implementation and maintenance costs of MPC deployment without compromising performance. We also critically assess the presented approach by pointing out the remaining challenges and open research questions.

Drgona, Jan↗
DOE OSTI2021

U.S. Industrial and Commercial Motor System Market Assessment Report Volume 1: Characteristics of the Installed Base

Motor systems are an integral part of our industrial and commercial facilities. They provide the motive force behind the fans, pumps, compressors, chillers, and conveyors in these facilities. Given their centrality to any facility’s operations, they are a critical energy end use to understand, particularly when developing technologies and policies to meet sustainability goals, improve productivity, and enhance resilience. In the late 1990s, the U.S. Department of Energy (DOE) conducted two seminal studies to better understand the installed stock and energy savings opportunities of industrial and commercial motor systems. In the industrial sector, The United States Industrial Electric Motor Systems Market Opportunities Assessment used primary data collected through onsite assessments and led to a greater understanding of the installed base of motor systems, their characteristics, and the opportunities for energy savings (U.S. Department of Energy, 2002). Notable findings included the following: Industrial motor systems consumed 679 billion kilowatt-hours (kWh) in 1994, representing 23 percent of U.S. electricity consumption. Cost-effective energy efficiency measures could result in 62-104 billion kWh of energy savings annually. Sixty-two percent of the energy savings potential was from fan, pump, and compressor end use equipment. Nearly half of motor system electricity consumption was attributable to approximately 3,500 facilities, or 1.5 percent of U.S. manufacturing facilities. Opportunities for Energy Savings in the Residential and Commercial Sectors with High Efficiency Electric Motors provided an evaluation of the installed stock of motor driven equipment in U.S. commercial and residential buildings and opportunities for utilization of high efficiency motors and variable speed technologies (Arthur D. Little, 1999). Notable findings included the following: Commercial motor systems consumed 343 billion kWh in 1995, with refrigeration and space conditioning constituting 93 percent of the total. Cost-effective energy efficiency measures could result in 51 billion kWh of energy savings annually. Due in no small part to these seminal studies, motor system technologies and usage characteristics have changed drastically since the late 1990s. Greater awareness of cost-effective strategies for reducing motor system electricity consumption have been developed and deployed. This includes several software tools, literature, and utility and government programs promoting energy efficiency improvements in motor driven systems. Additionally, several rounds of energy efficiency standards have been enacted, resulting in improved installed motor efficiency. The cost of variable speed drives has dropped substantially, and combined with utility rebate programs, has led to their greater adoption. Further, since these results were published, the U.S. manufacturing sector has undergone a massive transformation. Due to global competition, some sectors have relocated operations overseas. Others have brought operations onshore to avail low cost and abundant natural gas. Additionally, automation and robotics have pervaded the entire sector. Consequently, these two reports likely do not represent the current state of motor driven systems in U.S. industrial and commercial facilities. As cited in recent studies, the lack of current information on motor system electricity consumption and use characteristics limits the ability to conduct analysis on energy savings potential, develop technologies to address energy and productivity gaps, and develop programs to promote energy efficiency practices and technologies for motor systems (International Energy Agency, 2007; UNIDO, 2010; McKane and Hasanbeigi, 2011; Waide and Brunner, 2011). Specifically, the lack of information affects a range of stakeholders: Governments must rely on outdated information when setting research agendas, developing policies, and designing energy efficiency programs and offerings. Utilities and energy efficiency programs cannot identify the current market needs or potential impact when designing rebate and energy efficiency programs. Electric grid planners cannot identify motor system usage characteristics when developing plans to support the resilience of the electric grid. Manufacturers of motors, motor driven equipment, and drives are hampered when developing technologies to meet the needs of their market. Motor system end users are limited in their ability to identify energy saving opportunities within their own facilities because they do not have reliable benchmark information. In response to the lack of current information and analysis on industrial and commercial motor systems, the DOE initiated an update to these two studies. Launched in 2016 and led by Lawrence Berkeley National Laboratory (LBNL), the Motor System Market Assessment (MSMA) provides an updated, more comprehensive assessment of the installed stock of motor systems in both the industrial and commercial sectors, a review of the supply chains supporting motor and drives in the U.S., and the performance improvement opportunity available from using best available technologies and maintenance and operation practices. The outcomes of the MSMA are documented in three U.S. Industrial and Commercial Motor System Market Assessment reports, with this report being the first listed: 1. Volume 1: Characteristics of the Installed Base (this report) documents the findings on the installed base of motor systems in the U.S. industrial and commercial sectors. Quantification of energy savings potential is not documented in this report but in Volume 3. 2. Volume 2: Motors and Drives Supply Chain Review reviews the state of supply chains for motors and drives installed in U.S. industrial and commercial facilities, focusing on advanced motor and drive technologies and their constituent materials. 3. Volume 3: Energy Savings Opportunity analyzes the energy performance improvement opportunity for the installed base of U.S. industrial and commercial motor systems. This report has been prepared as a reference for motor system stakeholders. It provides factual information as could be best determined by the assessment results and avoids speculating on any findings.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Fuel Assembly and Irradiation Parametric Study for Extended-Enrichment and High-Burnup Light-Water Reactor Spent Nuclear Fuel in Dry Storage Casks and Transportation Packages

There is an increased interest in operating commercial light-water reactors (LWRs) in the United States with improved economics that would result from longer fuel cycle lengths, fewer and shorter refueling outages, and fewer fuel assemblies requiring storage at the back end of the fuel cycle. To support this, fuel discharge burnups, as well as initial 235 U enrichments, must be higher than those used in current commercial LWRs. The typical upper limit considered for assembly average burnup in this report is 75 gigawatt-days (GWd) per metric ton of uranium (MTU), as opposed to the current typical upper bound of approximately 62 GWd/MTU. The upper limit considered for initial 235 U enrichment is 8 weight percent (8 wt %), as opposed to the current regulatory limit of 5 wt %. The enrichment range from 5 to 8 wt % is referred to in this report as extended enrichment . To investigate the effect of high burnup and extended enrichment conditions on dose rates and burnup credit for dry storage casks and transportation packages, a fuel assembly and irradiation parametric study was performed. The conclusions from this study will assist U.S. Nuclear Regulatory Commission staff in reviewing applications for dry storage casks and transportation packages that contain high-burnup and extended enrichment fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling and Analysis of a Thermophotovoltaic Integrated Self-Powered Furnace

This work investigates the energy efficiency and carbon reduction potential of self-powered residential building heating equipment using a thermodynamic modeling approach. An integrated thermophotovoltaic power module and residential scale furnace system (40,000 Btu/h) were modeled and studied in detail to assess the influence of different design configurations on primary energy efficiency. Operational characteristics such as total power generation, electrical efficiency, and heat recovery were examined in a self-powered system configuration. A sensitivity analysis was conducted to determine the influence of the electric grid’s carbon dioxide footprint (carbon intensity) and the cost of electricity on the environmental, as well as the economic, benefit associated with the self-powered configuration. Compared with a traditional furnace powered by an electric grid at a carbon intensity of 0.5 kg CO 2eq /kWh EL , the self-powered furnace was shown to decrease the annual carbon dioxide emissions by approximately 550 kg (~75% reduction), while also saving more than USD 200 in utility expenses, annually. Additionally, the carbon emission reduction potential of blending different concentrations of hydrogen in natural gas fuel was also studied.

14 SOLAR ENERGY↗
DOE OSTI2020

Water and Ammonia Desorption from 5A Zeolite

Demonstrate adsorption and desorption characteristics of water and ammonia on 5A zeolite during regeneration conditions using a Thermogravimetric analyzer (TGA) on samples previously evaluated by residual gas analyzer (RGA)/mass balance. Molecular sieve (zeolite) adsorbents have the following characteristics: Crystalline alumino-silicates containing pores or 'cages' that adsorb water or other molecules. Polar compound are adsorbed with high loading even at very low concentrations of the fluid. A 5A molecular sieve was chosen for its ability to adsorb both water and ammonia (NH{sub 3}) from a carrier gas stream. Preliminary RGA/mass balance testing indicated: 5A zeolite captured both NH{sub 3} and water as anticipated based on molecular size; Reloading showed limited success after zeolite heated beyond manufacturer recommended operating conditions: Steam dealumination theorized cause; Tests were reprioritized to focus on NH{sub 3} removal and end of life potential holdup of NH{sub 3} and water on 5A zeolite material. TGA measures a sample's weight and volatile composition as it is heated or cooled in a furnace. TGA heats sample to release volatile materials, (change in mass is measured). MS detects/measures low levels of impurities. Perform TGA desorption of zeolite samples taken at the following conditions tested with RGA (RGA results not shown here): Initial Bakeout conditions: hydrogen carrier gas, slowly raise temp to approximately 300 deg. C until desired moisture content <<100 ppm is achieved. Water loading (re-loading): argon carrier gas through bubbler until breach is observed. Ammonia loading (re-loading): hydrogen carrier gas through system. Regenerations/End of life conditions: hydrogen carrier gas, raise temp to approximately 300 deg. C, then slowly to 450 deg. C (end of life). After end of life testing less than 4% of the total mass was residual ammonia or water at bakeout conditions. Peak water and ammonia removal occurred around 275 deg. C and 150 deg. C respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Multiphysics Time-Integration for Turbulent Combustion at the Exascale

Turbulent reacting flow systems are often modeled with coupled time-dependent partial differential equations (PDEs). Solving such equations can easily tax the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This generic operator-splitting strategy leads to a set of ordinary differential equations (ODEs) that need to be solved as part of an "outer-loop" time-stepping approach. In many combustion applications, the ODEs to be solved can be very stiff, exhibiting timescales that span many orders of magnitude. The SUNDIALS library provides a plethora of robust time integration algorithms for solving these ODEs on exascale-capable computing hardware, yet for many complex applications (such multicomponent fuels or emissions predictions), the chemical models remain too complex to solve using reasonable resources. The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of the simulations. In this talk, I will discuss the use of the SUDIALS library of ODE solvers together with automatic code generation tools to solve complex turbulent reacting flow problems using QSSA models.

chemistry↗
DOE OSTI2024

Helicity evolution at small x : quark to gluon and gluon to quark transition operators

We include the quark to gluon and gluon to quark shock-wave transition operators into the small Bjorken-x evolution equations for helicity in the flavor-singlet channel derived earlier in [1,2,3]. While such transitions do not affect the large-N c version of the evolution equations for helicity, the large-N c & N f equations are affected. (N c and N f are the numbers of quark colors and flavors, respectively.) We derive the corresponding corrected large-N c & N f equations for the polarized dipole amplitudes contributing to the flavor-singlet quark and gluon helicity distributions in the double-logarithmic approximation (DLA), resumming powers of α s ln 2 (1/x) with αs the strong coupling constant. We solve these equations iteratively and extract the polarized splitting functions up to four loops. We show that our splitting functions agree with the fixed-order perturbative calculations up to and including the existing three-loops results [4,5,6,7]. Similar to the large-Nc helicity evolution in the shock-wave approach [8], our large-N c & N f small-x splitting functions agree with those obtained in the infrared evolution equations framework from [9, 10] up to three loops, but appear to slightly disagree at four loops.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Identifying microstructures susceptible to pulverization in commercially irradiated high burnup UO2 under LOCA conditions

High burnup fuel fragmentation (HBFF) has been a concern in the nuclear industry for many years. When UO2 reaches and exceeds pellet average burnups around 55 GWd/tU, microstructural changes in the fuel pellet begin to occur that render the pellet susceptible to fine fragmentation under loss-of-coolant accident (LOCA) conditions. During a LOCA event, the cladding balloons and bursts, potentially releasing part of these fine fuel fragments into the reactor pressure vessel. This process is known as fuel fragmentation, relocation, and dispersal (FFRD). Presently, the US nuclear industry is developing a safety basis for FFRD with the goal of increasing burnups and pressurized water reactor cycle lengths. Increased cycle lengths would push the burnup limits of the fuel rods past the known threshold for HBFF susceptibility. In this work, the microstructural impact on HBFF was investigated by comparing the microstructures of as-irradiated rod segments to their post-LOCA tested counterparts. This investigation found that varied operational histories result in different microstructural evolutions across the pellet radius, which impacts the fragmentation behavior and radial location of the fragmentation. It is noted that regions with a high bubble density and a high density of grain boundaries fragment under LOCA-relevant conditions. In addition to the fragmentation that has been previously reported at the periphery of the fuel, fragmentation was also noted in the dark zone toward the fuel center. Analysis of the power histories of the fuel samples suggests that the dark zone forms in the central regions of the pellet between temperatures of approximately 740 and 960 °C.

McKinney, Casey [ORNL] (ORCID:0000000335383614)↗