Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance simulation”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluating the Impact of Tritium Permeation Membrane Performance and Direct Internal Recycling on Fusion Fuel Cycle Efficiency Using TMAP8

An efficient fuel cycle is vital to sustainable and cost-effective energy generation in fusion systems. Since tritium is not widely available, fusion systems must breed their own tritium for sustainable fusion deuterium-tritium reactions. An inefficient fuel cycle increases the tritium inventory needed for operations, which increases costs, constraints on tritium management systems, and safety concerns. A fuel cycle model is a powerful tool for understanding tritium inventories and flow rates across all systems in the fuel cycle. By simplifying the technical details into time-dependent tritium flow rates and inventories, the model can simulate the entire fuel cycle with high computational efficiency, even for technologies that are still under development. It can therefore quantify the impact of new tritium management technologies on fuel cycle efficiency. To evaluate the impact of key components on reducing tritium inventory, we are using and expanding an existing fuel cycle models based on latest advancements in fuel cycle research. The new model integrates Tritium Permeation Membrane (TPM) and Direct Internal Recycling (DIR) to enhance tritium transport from blanket breeders and plasma exhaust. These fuel cycle models are implemented in TMAP8 (Tritium Migration Analysis Program, version 8), a MOOSE-based open-source application designed to provide cutting-edge capabilities for tritium transport and fuel cycle modeling. The study aims to demonstrate the extensibility of existing fuel cycle modeling capability in TMAP8 and to offer a proof-of-principle design for future fusion plant systems. The presentation will cover the performance of fuel cycle modeling capabilities available in TMAP8, highlight advancements in fuel cycle research, and present a sensitivity analysis of these models. The results underline potential approaches and technology solutions to lower tritium inventory requirements, highlighting their role in shaping the future of fusion energy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗
DOE OSTI2022

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

BISON validation of FeCrAl cladding mechanical failure during simulated reactivity-initiated accident conditions*

Here, a reactivity-initiated accident (RIA) is a postulated design basis accident in light water reactors (LWRs) in which a rapid reactivity insertion induces a fission rate increase and a fuel pellet temperature rise. During RIA, the fuel pellet thermally expands and may cause pellet-cladding mechanical interaction (PCMI). Separate effects PCMI tests were performed on C26M FeCrAl cladding tube samples, introducing biaxial stress via a well understood modified burst test (MBT) system. A high-speed camera in the MBT system captured the projections, covering a 360° view of the cladding deformation and enabling a digital image correlation (DIC) method to quantify the surface strains with high fidelity. Representative hot zero-power RIA mechanical loading conditions were applied to the sample, and the test duration ranged from 20 to 500 ms at an average temperature of 573 K. BISON finite element–based fuel performance code modeling was performed against the high-fidelity DIC data produced from the MBTs. Validation calculations were conducted with 2D models and systematically compared with test data of the burst time, burst pressure, burst hoop strains, and hoop strain rates. Based on the behaviors from the separate effects test, BISON calculations satisfactorily predicted the cladding deformation behaviors. Sensitivity analysis was conducted to identify highly influential mechanical properties responsible for the cladding failure behavior during the MBT experiments. The results highlight the significance of the cladding's mechanical strength in governing cladding strain, followed by the significance of the mechanical interactions between the pellet and the cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗
DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Improving stability of ALD ZrN thin film coatings over U-Mo dispersion fuel

Atomic layer deposition (ALD) of ZrN is a candidate technology for coating U-Mo dispersion fuel as a diffusion barrier coating. During the early development stages of the coating the ZrN deposition showed to be mechanically unstable and ultimately resulted in spalling. Based on experimental outcomes it was found that instabillity can be eliminated through the introduction of a thin amorphous (a) Al 2 O 3 interlayer coating deposited in between the ~1-μm thick ZrN and the U-Mo. substrate. To elucidate these findings simulations were performed with density functional theory (DFT) to measure work of adhesions at different interfaces while finite element modeling (FEM) was performed to measure the residual stress distribution. DFT indicated that the ZrN coating deposited over (a)-Al 2 O 3 is ~2.5 times stronger when compared to direct deposition over U or UO 2 substrate. While calculations from FEM recognized; (1) large stress concentrations can originate from the irregularly distributed native surface oxides (UO 2 ) over U-Mo and (2) stress concentrations can be reduced if those surface UO 2 can be modified to a uniform layer (achieved after application of a 8 nm (a)-Al 2 O 3 interlayer). These results provided explanations and confirmed the role played by interlayer in eliminating the original mechanical instability of the ZrN.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Evaluation of Sb-Nd and Te-Nd phases within the U-Zr fuel matrix and their interactions with HT9 alloy

Antimony (Sb) and tellurium (Te) were investigated as potential additives for U-10Zr (wt.%) metallic fuel to limit the fuel-cladding chemical interaction (FCCI) with HT-9 alloy. Neodymium (Nd) was utilized to simulate the formation of lanthanide-based solid fission products which are known to play a detrimental role in FCCI. Fuel alloys of U-Zr-Sb-Nd and U-Zr-Te-Nd were evaluated in their annealed condition and compared against their as-cast conditions. Isothermal diffusion couple experiments were performed between U-Zr-Nd, U-Zr-Sb-Nd, and U-Zr-Te-Nd against the cladding alloy HT9 to evaluate the effectiveness of the additives to stabilize Nd within the fuel alloys, as well as investigate the interaction regions that form between the different fuel alloys and HT9. Further, SbNd and Sb 3 Nd 4 , and TeNd are found to be the primary neodymium-based phases formed in the U-Zr-Sb-Nd and U-Zr-Te-Nd alloys, respectively. The zirconium-based phase, Zr 2 Sb, is also found to form within the former alloy. All phases were found to remain stable through the diffusion experiments and exhibited no interaction with HT9 constituent elements. Preferential interaction between Nd with additivities Te and Sb compared to constituting elements in HT9 was further verified based on density functional theory (DFT) calculated enthalpy of mixing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multi-Physics Simulations of Molten Chloride Fast Reactor using Nek5000 and PROTEUS-NODAL

This report documents the FY20 work to demonstrate a specific multi-physics capability within NEAMS tools to model molten chloride fast reactors (MCFR). Specifically, the CFD code Nek5000 and the neutronics code PROTEUS-NODAL were coupled off-line to calculate the delayed neutron precursor distribution within the flowing fuel salts in an MCFR core. An MCFR benchmark problem and a 3D test problem based on a conceptual MSR design were successfully solved through the coupled simulation of Nek5000 and PROTEUS-NODAL. In the coupled simulations for these test problems, Nek5000 provided high-quality flow field information, and PROTEUS-NODAL performed the neutron diffusion calculation along with the precursor drift model. The impact of precursor drift on the core status was also investigated. Additionally, the precursor drift model implemented in PROTEUS-NODAL was verified by making use of the reference precursor distribution generated from Nek5000. Results of these test problems suggest that the multi-physics capability of Nek5000 and PROTEUS-NODAL was properly developed for tracking delayed neutron precursor drift.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Front-End Engineering Design (FEED) Study for a Carbon Capture Plant Retrofit to a Natural Gas-Fired Gas Turbine Combined Cycle Power Plant

A comprehensive front-end engineering design (FEED) study has been undertaken for a post-combustion capture (PCC) unit located at Panda’s Sherman natural gas–combined cycle (NGCC) power plant in Sherman, Texas. This is described in a full and unredacted FEED study report with all supporting documents, numbering over 150, also publicly available.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Front-End Engineering Design (FEED) Study for a Carbon Capture Plant Retrofit to a Natural Gas-Fired Gas Turbine Combined Cycle Power Plant (2x2x1 Duct-Fired 758-MWe Facility with F Class Turbines)

A comprehensive front-end engineering design (FEED) study has been undertaken by Bechtel National Inc. (Bechtel) for locating a post-combustion capture and compression (PCC) unit at Panda’s Sherman natural gas–combined cycle (NGCC) power plant in Sherman, Texas. This is described in the unredacted FEED Study report (Attachment 1) with all supporting documents, numbering over 150. The Study Report is publicly available. Sizing of the PCC plant is based on treating an amount of flue gas equivalent to that produced when generating 420 MW, which is approximately 68% of the total flue gas emitted by the NGCC power plant operating at guarantee condition with duct burners off. A reduced power plant capacity factor was used for sizing the PCC plant because the gas turbines at the site often operate at reduced load due to the high penetration of renewable power in the ERCOT region. The cost of carbon capture is primarily driven by capital cost (and therefore is highly sensitive to capacity factor). Sizing the capture unit so that when used it is nearly always operating at full capacity is critical to the economic viability of the proposed investment.

03 NATURAL GAS↗
DOE OSTI2023

Transient Efficiency, Flexibility, and Reliability Optimization of Coal-Fired Power Plants - Final Program Review

This program developed an advanced model-based monitoring and model-predictive control algorithms for a coal fired power plant (CFPP), and deployed these algorithms in a real-time platform to demonstrate performance benefits for transient flexibility and plant operation efficiency. More specifically, the objectives were successfully achieved through a combination of (i) developing a high-fidelity transient plant model in Apros, which was used as a high-fidelity plant simulation between $100-50\% TMCR$ where TMCR denotes the turbine maximum continuous rating, i.e., baseload, (ii) developing a very fast physics-based reduced-order model (ROM) of the plant, which ran more than $100\times$ faster than real-time, enabling its use as real-time embedded model for model-based estimation (MBE) and model predictive control (MPC) (iii) implementing a real-time MBE based on ROM using a robust unscented Kalman filter (UKF) to continuously tune the ROM to match the measurements from high-fidelity Apros plant model despite significant plant-model mismatch, and thus, obtain a Digital Twin of the plant (iv) designing and implementing a real-time MPC with dual objectives of transient plant load tracking with high ramp rates and minimizing coal consumption, i.e., improving plant efficiency in the baseload-partload operation range of $100-50\% TMCR$. Each key element above was developed and tested individually, and has been reported in corresponding Topical Reports. Finally, all the individual elements were integrated in an overall closed-loop system, that was successfully tested in desktop Simulink test harness simulations with ROM or high-fidelity model as the plant. Thereafter, the Simulink implementation was used to auto-generate C-code and deploy as real-time Docker microservice containers in Linux, and validate that they can run in real-time in the hardware-in-the loop (HIL) setup and produce the same results as in Simulink. The results of the integrated simulation tests in Simulink as well as the real-time HIL deployment are documented in this final report, showing good load tracking for load ramps at $3-4\%/min$ ramp rates, and achieving up to $5.5\%$ reduction in coal relative to baseline operation at $50\% TMCR$ load. The desktop and HIL simulations show successful performance of the overall model based estimation and control solution and achieve the key objectives of the program for flexible, efficient and reliable operation of sub-critical coal fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

System code evaluation of near-term accident tolerant claddings during pressurized water reactor station blackout accidents

Following the Fukushima Daiichi nuclear accident in 2011, researches on Accident-tolerant fuels (ATFs) are currently of high interest in not only the nuclear industry but also governmental and international organizations. In this work, a quantitative evaluation of the performance of monolithic FeCrAl cladding and Cr-coated Zircaloy cladding has been performed for Pressurized Water Reactor (PWR) Station Blackout (SBO) accidents. A generic PWR model has been built in system thermal-hydraulics code TRACE based on the Surry Nuclear Power Station with counter-current natural circulation modelling capability for hotleg and steam generator U-tube components during the accidents. The base model results are then compared to MELCOR and RELAP simulations to verify the system component implementation in TRACE. Two PWR SBO scenarios were investigated, including: short-term SBO and long-term SBO with early reactor coolant pump (RCP) seal failure. These scenarios are defined to be very similar to the accidents studied in the State-of-the-Art Reactor Consequence Analysis (SOARCA) project. TRACE code is modified to reflect the oxidation kinetics of FeCrAl and Cr-coating. Larson-Miller creep rupture model is also implemented in TARCE using its built-in control systems to simulate the creep rupture of hotlegs. Additionally, the comparison between the TRACE models with and without the counter-current flow modeling resulted in significant difference when comparing ATF cladding to Zircaloy for short term SBO, while it marginal impacted the performance during long term SBO with RCP seal failure. For short term SBO, both ATF cladding underwent hot leg creep rupture ~20 min after Zircaloy cladding. While Zircaloy and Cr-coated cladding had generated significant amount of hydrogen gas (>10 kg) before the creep rupture event, FeCrAl cladding had only generated <0.5 kg of hydrogen gas. For long term SBO with RCP seal failure, significant hydrogen generation and fuel melting was predicted before hot leg creep rupture for the ATF cladding while providing only 10–20 min additional coping time compared to Zircaloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Source term analysis of FeCrAl accident tolerant fuel using MELCOR

It has been established that incremental improvements in beyond design basis accident performance can be achieved through accident tolerant fuel (ATF). However, they have the potential to recover margin with respect to conventional fuel and therefore enhance plant economics through uprate or cycle length increase. To realize this potential, it is necessary to quantify the reduction in source term due to use of ATF, and correspondingly how this is affected by increasing the cycle length and/or burnup. This requires development of a risk informed analysis methodology for ATF under high burnup conditions, which is being developed within LWRS. To this end, a MELCOR simulation of FeCrAl ATF in a recovered Large Break LOCA (LBLOCA) scenario has been developed, using a model based on the Zion Pressurized Water Reactor (PWR). The new user defined material capability, along with the inclusion of detailed neutronics- and depletion-derived parameters such as core power distribution, decay heat behavior, and fission product inventories, allows a more detailed simulation of FeCrAl clad material properties and behavior than has previously been possible using MELCOR. These detailed FeCrAl results were compared with a zircaloy clad model to investigate the differences between the two clad types and quantify the benefits of FeCrAl ATF with an 18-month cycle. Next, the fuel cycle with ATF was extended to 24 months, to determine whether any additional safety margin provided by FeCrAl ATF could be leveraged to implement high burnup FeCrAl-clad fuel while retaining the same operating and safety limits as current zircaloy-clad fuel. For the particular scenario analyzed, the delay to fuel failure from using ATF was of the same order as the LOCA recovery time, and hence significant in reducing fission product release, with the 24-month FeCrAl core performing better than the 18-month zircaloy core. It is noted that for other transients, the reduction in release due to using FeCrAl may be less significant. Furthermore, the material model developed here can be used in such further studies in support of determining the overall source term.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Exploring Building Retrofit Strategies Using AutoBEM Under Future Weather Scenarios

This study evaluates the long-term effectiveness of energy conservation measures (ECMs) on building energy consumption using AutoBEM, a scalable modeling framework driven by the high-resolution Model America dataset. We simulated 18,951 buildings in Flagstaff, Arizona under four climate scenarios using Future Typical Meteorological Year (fTMY) weather files for six time periods spanning from 1980 to 2099. Six ECMs were analyzed across electricity and gas usage, including HVAC fuel-switching, insulation upgrades, and infiltration control. While some measures, such as reducing space infiltration by percentage, showed minimal or even negative impact on total energy savings at the aggregate level, they proved highly effective for specific building types. Conversely, HVAC electrification offers high gas reduction but shifts demand to electricity, highlighting critical trade-offs under different climate trajectories. Building-type-specific analysis under SSP5-RCP8.5 (2080–2099) revealed significant variation in ECM performance, underscoring the need for targeted retrofit strategies. This study demonstrates the power of combining fTMY projections with large-scale simulations to inform data-driven retrofit planning.

Chowdhury, Shovan [ORNL]↗
DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗