Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fuel performance modeling”

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,279 records · Page 71

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

EVALUATION OF HOT CHANNEL FACTOR FOR SODIUM-COOLED FAST REACTORS WITH MULTI-PHYSICS TOOLKIT

The evaluation of hot channel factor (HCF) is of great significance to the quantification of safety margins for reactor designs. In this paper, HCFs for a sodium-cooled fast reactor (SFR) are evaluated with the Simulation-based High-efficiency Advanced Reactor Prototyping (SHARP) toolkit, which is developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign of DOE for multi-physics reactor performance and safety simulations. The high-fidelity neutronics and thermal hydraulics solvers PROTEUS and Nek5000 in the SHARP toolkit are coupled to perform the multi-physics simulations for HCF evaluation. The HCFs induced by cladding manufacturing tolerance, fissile content mal-distribution, wire orientation and uncertainties on the cladding, coolant, and fuel properties are evaluated for a reference core SFR design (AFR-100). The HCFs calculated with the SHARP toolkit are compared to legacy HCFs for similar reactor types. The comparison demonstrates the reduction or elimination of modeling uncertainties in the calculation of HCFs using high fidelity advanced modeling and simulation tools without the need of expensive experiments. Moreover, the reduction of the uncertainties on HCFs evaluation allows an increase in nominal parameters and safety margin, which in turn improves the economic competitiveness of the SFR.

Hot channel factor (HCF)↗

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Derived Dynamic Operating Reserve Requirements in High-Renewable Power Systems

Accurately forecasting wind and solar power output poses challenges for deeply decarbonized electricity systems. Grid operators must commit resources to provide reserves to ensure reliable operations in the face of forecast errors, a process which can increase fuel consumption and emissions. We apply neural network-based machine learning to expand the usefulness of median point forecast data by creating probabilistic distributions of short-term uncertainty in demand, wind, and solar forecasts that adapt to prevailing grid conditions. Machine learning derived estimates of forecast errors compare favorably to estimates based on incumbent methods. Reserves derived from machine learning are usually smaller than values derived using incumbent methods, which enables fuel savings during most hours. Machine learning reserves are generally larger than incumbent reserves during times of higher forecast error, potentially improving system reliability. Performance is tested using multi-stage production simulation modeling of the California Independent System Operator (CAISO) system. Machine learning reserves provide production cost and greenhouse gas (GHG) emission reductions of approximately 0.3% relative to historical 2019 requirements. Savings in the 2030 timeframe are highly dependent on battery storage capacity. At lower levels of battery capacity, savings of 0.4% from machine learning reserves are shown. Significant quantities of battery storage are expected to be added to meet California's resource adequacy needs and GHG reduction targets. Addition of these batteries saturate reserve needs and results in minimal within-hour balancing costs in 2030.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluation of XCT for Matrix Density Measurement of Particle Fuel Forms

Particle fuel forms generally consist of a dispersion of fuel, such as tristructural isotropic (TRISO) particles, within a refractory matrix (e.g., graphite or silicon carbide). The density of matrix materials for particle fuel forms is of interest for modeling fuel form strength and thermal properties and may be specified as a quality control parameter, depending on reactor design. Some of the uncertainty associated with traditional, manual approaches can be eliminated by performing x-ray computed tomography (XCT) on the fuel forms and applying image processing methods to generate a precise count of the number of particles. This also removes the need to include determination of particle count within each individual fuel form during fabrication. Unfortunately, reconstruction artifacts from high-Z uranium-bearing kernels prevent accurate measurement of individual particle volumes using this approach, so the use of mean particle mass and volume are still necessary for computation of average fuel form matrix density. This method of using XCT to count particles in individual fuel form for determination of average matrix density was applied to three archived compacts from the Advanced Gas Reactor Fuel Development and Qualification (AGR)-1 campaign, four archived compacts with uranium carbide/uranium oxide (UCO) TRISO from the AGR-2 campaign, and three archived UO 2 -TRISO compacts from the AGR-2 campaign. The resulting density values were compared with those previously reported, showing slight changes due to uncertainties in the previously used number of particles in each of these cylindrical, graphite matrix compacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High Temperature High Vacuum Mechanical Property Assessment of Zirconium Nuclear Fuel Cladding

This report presents the mechanical characterization of a specific Zry-4 cladding batch serving as the foundation for a diverse range of fuel performance research at Oak Ridge National Laboratory (ORNL). This effort supports research needs for the U.S. Department of Energy (DOE), particularly regarding evaluating accident tolerant fuel (ATF) cladding coating concepts, expanding understanding of cladding response to loss-of-coolant accidents (LOCA) transients, refining post-critical heat flux (CHF) limits (t@T), and upcoming irradiation campaigns. The central objective was to define the baseline performance of the substrate Zircaloy-4 (Zry-4) material leveraged across ORNL Advanced Fuel Campaign (AFC) efforts through controlled high-temperature vacuum tensile testing. This work begins to address gaps in existing models where implementation based on nominal heat-treatment labels, such as stress relief annealed (SRA), often fail to capture the interplay of recovery, recrystallization, and grain growth. To quantify this, data was benchmarked against the Pacific Northwest National Laboratory (PNNL) stress strain model to determine where this material falls in comparison to assumed values for materials in the same heat treatment regime. Analysis of the tensile data revealed that this specific SRA batch exhibits a transitional microstructural state best described by an effective cold-work (CW) parameter of 0.09, diverging from the previous estimation of 0.5 for SRA materials. Additionally, comparative testing of Cr coated specimens demonstrated no distinct difference in axial strength relative to the bare substrate. This suggests that the strengthening benefits of Cr coatings observed in burst scenarios are driven by residual stress mechanisms acting solely in the hoop direction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An optimized rotation-axis model-following controller for STS orbiter attitude maneuvers

This paper describes a new controller for automatic attitude maneuvers of the Space Shuttle Orbiter. The controller incorporates a new solution algorithm for partially linearized, fuel-optimal maneuvers for RCS jet controlled spacecraft. These solutions provide model trajectories and jet commands that are easily tracked using a feed-forward, feedback controller structure. The result is a substantial performance improvement over the current Orbiter attitude maneuver controller. The paper discusses optimal maneuvers for this type of vehicle, giving detail of the solution algorithm. Closed-loop controller design and implementation is discussed. Detailed performance results are presented.

Redding, D. C.↗

BETO 2021 Peer Review: 2.1.0.100 - Biochemical Platform Analysis

The objective of this project is to perform techno-economic analysis (TEA) to guide Biochemical Platform efforts, utilizing models for purposes of setting future R&D targets and tracking performance progress against those targets. Outcomes of our work are leveraged by BETO to guide program plans, as well as by other NREL/partner projects to quantify the impact of research on key technology barriers and to prioritize future efforts. This project provides high impact and relevance through establishing “bottom-up” TEA models as a basis for understanding the technical feasibility to meet “top-down” BETO cost targets. By providing a framework to translate technical performance to cost reductions in a biorefinery, our TEA models may be leveraged to maximize the efficiency of research funding towards the most economically impactful priorities, ultimately in support of BETO’s 2030 fuel cost targets below $2.5/GGE. In order to mitigate a key risk/challenge to this project in overly-constraining our analyses to a singular technology focus or TEA metric, our approach continuously re-assesses opportunities for better optimization and alternative technology pathway options, while maintaining close interaction with other BETO analysis partners. We have made numerous recent accomplishments rooted around identifying and working with the researchers to solve key technical and TEA/LCA challenges, reflected through notable improvements in State-of-Technology (SOT) updates over prior benchmarks.

biochemical↗

Benefits of Dual Fuel Heat Pump Grid-responsive Control: A Model-based Control Optimization Approach Using Building and Equipment Co-simulation

Conventional dual fuel heat pumps lack the intelligent control mechanisms to efficiently manage the switch between heat pump and furnace, leading to sub-optimal energy usage and, in some cases, increased operating costs. To resolve this gap, this study applies optimized control on hybrid heat pumps. With a focus on equipment control strategies, we compare the performances of five spacing heating equipment, including a conventional heat pump (HP), a conventional furnace, a dual fuel heat pump (DFHP) with conventional control, a dual fuel heat pump with smart control, and a novel seamlessly fuel flexible heat pump (SFFHP). While DFHP runs on either gas or electricity at any given moment, SFFHP concurrently consumes gas and electricity by continuously optimizing the proportion of each. In this research, a co-simulation framework is developed by integrating a building envelope model with a physics-based heat pump simulation model to analyze the benefits of grid-responsive controls of DFHP and SFFHP. The model-based optimal controls adjust the operation of the heat pump and gas furnace based on utility price signals and marginal grid emission to minimize utility cost and CO 2 emissions for multiple climate zones, different utility tariffs, and marginal grid emission scenarios. Case studies in Chicago and Los Angeles demonstrate that SFFHP and DFHP, with model-based optimal control, can deliver significant reductions in peak demand, utility cost, and CO 2 emission. In Chicago, SFFHP and smart controlled DFHP yield up to 64.7% and 61.7% utility cost reduction and up to 15.7% and 8.5% CO 2 emission reduction compared to the gas furnace. In Los Angeles, SFFHP and smart controlled DFHP achieve up to 43.6% and 40.1% utility cost reduction and up to 13.8% and 14.1% CO2 emission reduction compared to conventional heat pumps. In conclusion, by leveraging the fuel flexibility nature of dual fuel heat pumps, the model-based control optimization approach makes dual fuel heat pump an attractive option for demand response programs.

Control↗

Nitric Oxide and Oxygen Air-Contamination Effects on Extinction Limits of Non-Premixed Hydrocarbon-Air Flames for a HIFiRE Scramjet

Unique nitric oxide (NO) and oxygen air-contamination effects on the extinction Flame Strength (FS) of non-premixed hydrocarbon (HC) vs. air flames are characterized for 7 gaseous HCs, using a new idealized 9.3 mm straight-tube Opposed Jet Burner (OJB) at 1 atm. FS represents a laminar strain-induced extinction limit based on cross-section-average air jet velocity, Uair, that sustains combustion of a counter jet of gaseous fuel just before extinction. Besides ethane, propane, butane, and propylene, the HCs include ethylene, methane, and a 64 mole-% ethylene / 36 % methane mixture, the writer s previously recommended gaseous surrogate fuel for HIFiRE scramjet tests. The HC vs. clean air part of the work is an extension of a May 2008 JANNAF paper that characterized surrogates for the HIFiRE project that should mimic the flameholding of reformed (thermally- or catalytically-cracked) endothermic JP-like fuels. The new FS data for 7 HCs vs. clean air are thus consolidated with the previously validated data, normalized to absolute (local) axial-input strain rates, and co-plotted on a dual kinetically dominated reactivity scale. Excellent agreement with the prior data is obtained for all 7 fuels. Detailed comparisons are also made with recently published (Univ. Va) numerical results for ethylene extinction. A 2009-revised ethylene kinetic model (Univ. Southern Cal) led to predicted limits within approx. 5 % (compared to 45 %, earlier) of this writer s 2008 (and present) ethylene FSs, and also with recent independent data (Univ. Va) obtained on a new OJB system. These +/- 5 % agreements, and a hoped-for "near-identically-performing" reduced kinetics model, would greatly enhance the capability for accurate numerical simulations of surrogate HC flameholding in scramjets. The measured air-contamination effects on normalized FS extinction limits are projected to assess ongoing Arc-Heater-induced "facility test effects" of NO production (e.g., 3 mole-%) and resultant oxygen depletion (from 21 to 19.5 %), for testing the "64/36" surrogate fuel in Langley s Arc-Heated Scramjet Test Facility for HIFiRE engine designs. The FS results show a generally small (< 4 %) "nitric oxide enhancement" effect, relative to clean air, for up to 3 % NO (freestream Mach number up to 7 in Arc Jet testing). However, a progressively large "oxygendeficiency weakening" effect develops. For 3 % NO, a net weakening of 26 % in FS is derived for the "64/36" fuel vs. air. The corresponding net weakening for pure ethylene is 20 %. A number of practical recommendations regarding facility test effects are offered.

Pellett, Gerald L.↗

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

ZrO 2 /MWCNT as a Support Platform for Acid and Metal Catalysis in Water: The Investigation of Pd-WO x -ZrO 2 /MWCNT Catalysts (Final Technical Report)

This final technical report covers the final renewal period, 08/01/2019 to 07/31/2020, and the no cost extension 08/01/2020 to 07/31/2021, that was required because of the pandemic interruption. The proposed research was in the area of sustainable recovery of fuels and/or chemicals from biomass; it was by design intended to address fundamental science of the catalysts and that required using a model reaction. The model reaction molecule was phenol because it is a reasonably simply molecule that can be derived from lignin, and requires both hydrogenation and dehydration catalytic reactions that can be performed in water. The needed bifunctional catalyst can be provided by a metal for hydrogenation and by an acid (preferably Brønsted acid) for the dehydration function. Therefore, phenol conversion to an alkane in liquid water at a practical temperature, e.g., 200°C, that will not require extreme pressures, captures many of the complications of biomass derived fuels so it may reveal fundamental science that is relevant to development and implementation of biomass refining.

36 MATERIALS SCIENCE↗

Model-Based Self-Tuning Multiscale Method for Combustion Control

A multi-scale representation of the combustor dynamics was used to create a self-tuning, scalable controller to suppress multiple instability modes in a liquid-fueled aero engine-derived combustor operating at engine-like conditions. Its self-tuning features designed to handle the uncertainties in the combustor dynamics and time-delays are essential for control performance and robustness. The controller was implemented to modulate a high-frequency fuel valve with feedback from dynamic pressure sensors. This scalable algorithm suppressed pressure oscillations of different instability modes by as much as 90 percent without the peak-splitting effect. The self-tuning logic guided the adjustment of controller parameters and converged quickly toward phase-lock for optimal suppression of the instabilities. The forced-response characteristics of the control model compare well with those of the test rig on both the frequency-domain and the time-domain.

Le, Dzu, K.↗

CFD-Chemical Model of One-Third Scale Demonstration of DOE Sealed Canister with ATR Fuel

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. A demonstration case for these DOE sealed canisters will be eventually performed with a one-third scale mockup that has been instrumented with thermocouples and gas concentration probes. This study seeks to model the planned canister for validation of the previously developed model, such that confidence in its long-term prediction can be increased. The deployment of the instrumented lid for online monitoring in intended for a period of 10+ years based on previous monitoring of commercial fuel storage; however, the model is still run for the previously used 50-year storage periods. This case is modeled with a G-value using a bi-linear function such that it decreases at higher dose rates to be consistent with experiments. Validation efforts of the model would likely revolve around the results for the 1st year, so the results of this timeframe are also highlighted. For dried fuel of the nominal decay heat (18W), the predicted hydrogen concentration is 0.19% after 1 year, 0.91% after 10 years, and 2.8% after 50 years, with a maximum pressure of 1.26 atm. Consistent with previous modeling, the decay heat of the fuel is the main factor that influences the results. For undried fuel in pure helium, the hydrogen concentration ranges from 0.23 to 1.2% after 1 year, 1.25-6.2% after 10 years and 8.8 to 18.84% after 50 years. For dried fuel in pure helium the hydrogen concentration ranges from 0.06% to 0.34% after 1 year, 0.35-2.11% after 10 years and 0.95 to 6.6% after 50 years. While long-term results in the presence of residual air are mostly the same, early reactions with O 2 can delay significant production of H 2 until it is consumed to form more water vapor, lowering the range to 0.14 to 1.0% after 1 year. The maximum absolute pressure that is reached across any scenario is 2.06 atm. In the event of residual air, the presence of nitric acid is possible in the range of 18-131 ppm after 1 year, 174-1180 ppm after 10 years, and 586-3500 ppm after 50 years. As long as the fuel is sufficiently dried, or of nominal decay heat, a 4% lower flammability limit of hydrogen will not be reached within a 10-year monitoring period.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

M3 milestone: Advanced contact 2021

This report describes advanced contact approaches and implementations in FY21 for use in emerging reactors with geometries that differ from LWR/EBR-II cylindrical fuel. Contact remains a challenging and crucially important component of our modeling toolkit. While significant improvements have been made through mortar and automatic differentiation, we explored advanced approaches (i.e., dual mortar) that yields better conditioned systems and will improve the robustness and performance of BISON contact problems. This effort is relevant to emerging reactor types, such as micro reactors and the corresponding complex and numerous contact surfaces. In this report, improvements of the mortar-based contact formulation are described. Benchmark problems and BISON assessment cases are included to demonstrate the improvements and advantages that are brought by the recent capability. Remaining issues are discussed and corresponding development plans are outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗