Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced Reactor Analysis”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermal Integration of Advanced Nuclear Reactors with a Reference Refinery, Methanol Synthesis, and a Wood Pulp Plant (Rev.1)

The present report is intended to provide process flow diagrams (PFDs) and energy and mass balance data sheets for a U.S. industrial sector subset with which nuclear heat and power could be integrated—a subset that includes the oil refining, methanol and pulp and paper industries. Coupling options for integrating nuclear energy into these industries are quantitatively outlined for reference systems, and future work will extend this analysis in greater detail. Opportunities for integrating small modular nuclear reactors (SMNRs) were investigated for each of the industrial process configurations. Aspen HYSYS and Cycle-Tempo models for a high-temperature gas-cooled reactor were developed to evaluate the proposed integration. This introductory evaluation provides a general description and assessment of the operating principles, reactor coolant core outlet temperature, and reactor size to be integrated with industry. The industrial processes of oil refining and the production of methanol, pulp and paper were simulated by using Aspen HYSYS, Aspen Plus, and the PRELIM (Petroleum Refinery Life Cycle Inventory Model) tool to develop process details. Cycle-Tempo models then extend the process modeling results to obtain net energy demands (e.g., heat, steam, and electricity) when accounting for process steam and waste heat recovery. This information is intended to foster the analysis of integrating an SMNR to decarbonize industrial facilities. The SMNR would provide reliable, competitive, and sustainable clean energy while reducing carbon emissions and other environmental impacts, such as water withdrawals, consumption, and contamination. The refining industry, exhibited in Figure ES1, is a leading consumer of fossil -fuel-based heat, power, and hydrogen in the U.S. industrial sector, generating over 164 million metric tons (MMT) of CO 2 emissions in 2023. The overall mass and energy pertaining to a generalized complex refinery in the United States is reflected in Figure ES1, along with energy metrics regarding integration with a nuclear power plant (NPP). Data sheets were developed to indicate the energy requirements for the overall refinery and each refinery process. The data sheet for the overall refinery is shown in Table ES2.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Tracking the phase changes in micelle-based NiGa nanocatalysts for methanol synthesis under activation and working conditions

The hydrogenation of CO 2 into high energy density fuels such as methanol, where the required H 2 is obtained from renewable sources, is of utmost importance for a sustainable society. In recent years, NiGa alloys have attracted attention as promising catalyst material systems for the hydrogenation of CO 2 into methanol at ambient pressures. They thus represent an energy-saving alternative to the Cu-based catalysts employed in today's catalytic industry that require high pressures for the CO 2 hydrogenation. However, the underlying reaction mechanisms for the NiGa system are still under debate. One of the challenges here is to unravel the evolution and coexistence of the different species in the heterogeneous NiGa catalyst system under activation and reaction conditions. To shed light on their evolution under activation in H 2 and their catalytic roles under CO 2 hydrogenation working conditions on well-defined Ni 3 Ga 1 and Ni 5 Ga 3 nanoparticle (NP) catalysts, we employed a multi-probe approach in this study. It included advanced machine learning-based analysis of operando X-ray absorption spectroscopy data combined with operando powder X-ray diffraction and near ambient pressure X-ray photoelectron spectroscopy measurements, as well as reactivity studies using bed-packed mass flow reactors. In addition, we employed atomic force microscopy and scanning transmission electron microscopy for structural characterization. Under H 2 activation at 1 bar total pressure, we concluded the formation of metallic Ni, starting for Ni 3 Ga 1 at 300 °C, and for Ni 5 Ga 3 at 400 °C. At higher temperatures, the formation of NiGa alloys follows. The α'-Ni 3 Ga 1 alloy phase is predominantly formed for the Ni 3 Ga 1 NPs, while the coexistence of α'-Ni 3 Ga 1 , δ-Ni 5 Ga 3 and Ga 2 O 3 phases is observed for the Ni 5 Ga 3 NPs after the H 2 activation. The formation of the Ga 2 O 3 phase also results in the presence of excess metallic Ni. Under CO 2 hydrogenation reaction conditions, Ga partially oxidizes again to form a Ga 2 O 3 -rich particle shell for both NP compositions, yet, to a larger extent for the Ni 3 Ga 1 NPs, which, in turn, feature a higher amount of excess Ni. We reveal that metallic Ni is responsible for the high selectivity of the Ni 3 Ga 1 NPs towards the production of methane in our catalytic tests. In contrast, the Ni 5 Ga 3 NPs display a strong selectivity toward methanol production (>92%), more than one order of magnitude higher than that for the Ni 3 Ga 1 NPs, which we ascribe to the presence of the δ-Ni 5 Ga 3 phase.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Siting Analyses for Elementl Power (Final CRADA Report)

This report summarizes siting evaluation assistance provided to Elementl Power under CRADA/NFE-23-09638 for suitability of advanced nuclear technologies to meet siting criteria from the Nuclear Regulatory Commission (NRC) and associated guidance documents including the Electric Power Research Institute (EPRI) siting guide and other proprietary datasets. Elementl Power provided sites for reactor siting evaluations using the Oak Ridge – Siting Analysis for power Generation Expansion (OR-SAGE) tool. Individual data packages for each site are provided separately.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Fabrication of yttrium hydride for high-temperature moderator application

We report given the superior thermal stability and highly attainable hydrogen density, yttrium hydride is an excellent high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. Large-scale, crack-free, bulk yttrium hydride is in high demand; however, fabrication of yttrium hydride is challenging and has not been demonstrated for nearly half century. The associated challenges are discussed herein. In response to these challenges, a hydriding system was designed and constructed at Oak Ridge National Laboratory and was used to successfully fabricate crack-free yttrium hydride in complex geometries at large scales. This was accomplished by precisely controlling the hydrogen’s partial pressure and the retort temperature, which was informed by the well-established thermodynamic properties of the binary H–Y system. Hydrogen content in as-fabricated hydride was determined by weight change, an approach which was considered reliable based on the use of ultra–high-purity yttrium, the absence of oxide phases up to levels detectable using x-ray diffraction (XRD), and the significant weight gain. Hydrogen distribution along one yttrium hydride rod was evaluated with XRD analysis on materials extracted from different locations on the rod. The results indicated a relatively homogeneous hydrogen distribution along the hydride rod, with <3% uncertainty in the fraction of the δ-phase hydride. In addition, significant efforts are being dedicated to establish a complete database summarizing the thermomechanical and physical properties of as-fabricated yttrium hydride and the irradiation response to facilitate its deployment as a high-temperature moderator in advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

AMTS update for STTR

Advanced nuclear reactor systems currently in development include designs that utilize higher energy neutron spectrums, higher temperatures, and corrosive environments that are quite different from the current fleet of nuclear power plants in use today. As a result, new structural materials must be developed and qualified for use to support these advanced systems. For example, research is currently being conducted on a new class of materials known as multi-principle-element-alloys. These alloys show some promise to be able to meet the structural needs of these challenging reactor systems. However, due to the vast available compositional space, and the limited amount of research conducted thus far, Idaho National Laboratory (INL)sought a way to advanced mechanical testing capability for highly irradiated structural materials to increase both throughput and higher fidelity data analysis through incorporation of digital image correlation. The Automated Mechanical Testing System accomplishes these tasks by automating the testing and incorporating digital image correlation to provide higher fidelity data. INL has now received this system and is currently performing functional testing and training. INL is also performing some modifications to the system to allow for installation into a remote shielded enclosure where testing of irradiated structural materials can be performed. After testing, training, and modifications of the system are complete, INL will install the system into a shielded enclosure and begin using the system to test irradiated structural materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A New Approach to Predict Hydrogeological Parameters Using Shear Waves from the Multichannel Analysis of Surface Waves Method

For near-surface contaminant characterization, the accurate prediction of hydrogeological parameters in anisotropic and heterogeneous environments has been a challenge since the last decades. However, recent advances in near-surface geophysics have facilitated the use of geophysical data for hydrogeological characterization in the last few years. A pseudo 3-D high resolution P-wave shallow seismic reflection survey was performed at the P Reactor Area, Savannah River Site, South Carolina in order to delineate and predict migration pathways of a large contaminant plume including trichloroethylene. This contaminant plume originates from the northwest section of the reactor facility that is located within the Upper Atlantic Coastal Plain. The data were collected with 40 Hz geophones, an accelerated weight-drop as seismic source and 1 m receiver spacing with near- and far-offsets of 0.5 and 119.5 m, respectively. In such areas with near-surface contaminants, a detailed subsurface characterization of the vadose zone hydraulic parameters is very important. Indeed, an inexpensive method of deriving such parameters by the use of seismic reflection surveys is beneficial, and our approach uses the relationship between seismic velocity and hydrogeological parameters together with empirical observations relating porosity to permeability and hydraulic conductivity. Shear wave velocity ( V s ) profiles were estimated from surface wave dispersion analysis of the seismic reflection data and were subsequently used to derive hydraulic parameters such as porosity, permeability, and hydraulic conductivity. Additional geophysical data including core samples, vertical seismic profiling, surface electrical resistivity tomography, natural gamma and electrical resistivity logs allowed for a robust assessment of the validity and geological significance of the estimated V s and hydrogeological models. The results demonstrate the usefulness of this approach for the upper 15 m of shallow unconsolidated sediments even though the survey design parameters were not optimal for surface wave analysis due to the higher than desired frequency geophones.

Engineering↗
DOE OSTI2022

Post-irradiation characterization of a high burnup mixed oxide fuel rod with minor actinides

The Advanced Fuels Campaign performed a series of irradiation tests of minor actinide-bearing mixed oxide fuel (MA-MOX), the so-called AFC-2C&D experiments, to investigate the transmutation of long-lived transuranic actinide isotopes contained in spent nuclear fuel via fast reactor technology at burnups exceeding 10 % fission of initial metallic atoms. This manuscript reports the test results derived from one of the five MA-MOX rodlets taken to higher burnup in the AFC-2D irradiation. This includes both non-destructive investigations, such as gamma and neutron spectrometry, and destructive investigations, such as fission gas release, ceramography, and chemical burnup analysis. In addition, the microstructure of the fuel was investigated using advanced electron microscopy techniques including electron backscatter diffraction (EBSD) and transmission electron microscopy (TEM). Here, it was observed with EBSD that the pellet had subdivision of the grains and the TEM observed migration of cladding material into the 5 metal precipitates in the fuel which could have been from the higher than desired oxygen/metal ratio. The TEM also showed an enrichment of Cr in fuel clad chemical interaction (FCCI) layer.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Microstructure of irradiated AGR TRISO particle buffer layers as measured by X-ray computed tomography

Shrinkage of the initially low-density buffer layer in tristructural isotropic (TRISO) coated fuel particles during irradiation is a well-known phenomenon with potential implications for fission product and actinide transport, as well as potential pyrocarbon fracture that in some cases has been observed to impact particle performance. During post-irradiation examination, the buffer layer's structure is commonly determined using 2D microscopy of a particle cross section or a series of particle cross sections at staggered depths. Although these methods provide a general idea of the irradiated buffer microstructure, they do not provide a full picture of the TRISO coating microstructure. By contrast, x-ray computed tomography (XCT) provides full 3D imaging of the TRISO particle. Particles from various compacts irradiated for the Advanced Gas Reactor Fuel Development and Qualification Program were imaged using XCT. Some of these particles were selected because of abnormal fission product inventories (e.g., low 137 Cs), whereas others were randomly selected from the center of the fission product inventory distribution. Qualitative and quantitative analysis techniques were applied to the randomly selected particles as representatives of typical TRISO behavior to study the post-irradiation structure of the buffer layer. These results showed that, while separation of the buffer and inner pyrolytic carbon layers was a common behavior during the radiation-induced shrinkage of the buffer, a portion of the original buffer/inner pyrolytic carbon interface remained intact throughout irradiation in nearly all cases. Furthermore, these results also indicated clear trends in the degree of buffer densification with irradiation temperature and fluence.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling the structure of molten salts: a review of high-energy x-ray diffraction (HEXRD) studies

High energy x-ray diffraction (HEXRD) has emerged as a pivotal technique for examining the structural properties of molten salts, which play vital roles in various high-temperature applications, such as energy storage, nuclear reactors, and metallurgy. This review provides a comprehensive analysis of recent progress in HEXRD studies of various molten salt systems, including molten halide salts, their mixtures, and salts with molecular anions, focusing on the insights into their atomic-scale arrangements. Additionally, we discuss methodological advancements that have improved the resolution and accuracy of HEXRD measurements, as well as recent developments in computational simulations, such as machine learning interatomic potentials, for interpreting results. This review aims to serve as a valuable resource for researchers in this field, offering a detailed overview of the progress of HEXRD studies of molten salts and suggesting pathways for further exploration.

Blockmon, Avery L. [Argonne National Laboratory (A↗
DOE OSTI2025

Investigation of radiated-power for low aspect ratio fusion plasmas

Investigating radiated power in fusion plasmas is of utmost importance to understand the effect of undesired impurities, such as metals, in present devices, or also desired impurities, such as noble gasses, to purposefully radiate a large fraction of power in future devices. These studies are especially important for high Z impurities which will play a crucial role in future generation fusion pilot plants. In this work, we have developed a power radiation analysis module, which is used to investigate 2D distributions of impurity densities and radiated power asymmetries caused due to both plasma rotation and the cooling rate dependence on temperature profile for the cases of experimental NSTX plasmas and designed scenarios for the spherical tokamak advanced reactor (STAR), both being low aspect ratio tokamaks. Two different atomic databases have been tested during this work to study their impact on the radiated power distribution, especially for high Z impurities. Also in PRAM, self-consistent calculations of two-dimensional electron, main ion, and impurity ion densities are carried out using one-dimensional input density, temperature, and rotation profiles. In the case of NSTX, discharges with high rotation of ∼ 170 km s, measured with charge exchange recombination spectroscopy, have been investigated. Rotation-induced charge separation, leading to an electrostatic potential, is calculated iteratively to a self-consistent solution while testing high Z impurities to observe any 2D asymmetry in the core radiated power due to centrifugal forces. The STAR design, being much larger (R = 4 m), is projected to have a much lower rotation, and is shown to have low rotation-induced asymmetries, on the order of ten percent or less, between the low field and high field sides. However, another effect not due to rotation but to the dependence of impurity cooling rates on temperature can lead to radiation peaking off-axis, near the plasma edge. This effect is noticeable for argon in NSTX, for example, but can also be enhanced for certain impurities at much higher temperatures projected for STAR (T e0 ~ 32 keV), for example for undesired tungsten or possibly desired xenon.

NSTX↗
DOE OSTI2020

Development of a Fast Spectrum Test Loop at the Advanced Test Reactor

The United States currently has no infrastructure to provide a fast neutron flux for research purposes. With the increased interest in Sodium Fast Reactor (SFR) technology, this type of infrastructure is necessary for further development of key components such as Versatile Test Reactor. The purpose of this research is to characterize the thermal hydraulics of a test loop capable of providing SFR like operating conditions. This includes showcasing feasibility in mechanical design, reaching a high fast to thermal neutron flux ratio ideally near 15, an inlet temperature of 600 K, a temperature increase of 50 K from inlet to outlet, and test chamber radial temperature profiles similar to SFR fuel rods. Existing infrastructure must be used to provide the neutron flux for the test loop, and it was decided that the Advanced Test Reactor (ATR) was the optimal choice due to its availability and previous work conducted similar to this project. To achieve the goals of this project, several mechanical designs were put forward with the most feasible design being chosen. The neutronics and thermal hydraulics analysis were conducted in tandem, with an iterative process used to maximize efficiency of both systems. Analysis of the neutronics system was done using Attila and MCNP while the thermal hydraulic analysis was done using RELAP5. This document outlines the thermal hydraulic analysis, though it includes the results from the mechanical design and neutronics analysis. The loop design includes a primary system with molten sodium, a set of “booster fuel plates” to increase the fast flux within the test section, a water system to cool the booster fuel plates, a helium system to separate the sodium and water systems, and a thermal neutron shield to decrease the thermal flux. Each of these systems is modeled in the neutronics and thermal hydraulic analyses. The RELAP5 model showcased the feasibility of the design proposed by this report. For all designs tested, the desired 50 K temperature increase from inlet to outlet was reached. The fuel rod radial temperature profile was lower than was expected, likely due to a lower rod thermal power than expected. A parametric study was also conducted to characterize the system with three fuel pins. This study showed the effect of thermal power, fluid pressure, fluid velocity, and inlet temperature on the system. ©Copyright by Anton Higgins December 11, 2019 All Rights Reserved

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Mid-Year Progress Update on ORNL Support for Developing the Guidance for Microreactor Manufacturing Licenses

Factory fueling and assembly, multi-site operation, and the associated transportation of advanced nuclear systems present both new opportunities and challenges for microreactor deployment. These areas have not been demonstrated under Nuclear Regulatory Commission rules and regulations. A goal of the Systems Integration and Analysis technical area under the US Department of Energy Microreactor Program for the current fiscal year is to articulate these challenges and provide industry recommendations. This mid-year progress report describes the progress and some of the contributions to this effort by the Advanced Reactor Engineering and Development section at Oak Ridge National Laboratory. The content provided here is preliminary and may change before being incorporated into a later report, which is expected to be completed in September 2022.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Updating and Refining BISON TRISO Assessment

Several advanced reactor concepts consider tristructural isotropic (TRISO) fuel in their designs. To support the advanced reactor industry, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program within the U.S. Department of Energy (DOE) is developing modeling and simulation capabilities within the BISON fuel performance code. Assessment of these capabilities against available experi- mental data, such as the Advanced Gas Reactor (AGR) program, is essential in instilling confidence in the code predictions. Associated uncertainties in model inadequacy, experimental measurements, and manufacturing are also of importance. This report discusses the updating and standardization of the existing assessment suite in BISON for TRISO fuel performance analysis. The cases considered the AGR-1, AGR-2, and AGR-3/4 experimental programs. The report also includes a preliminary investigation into the sources of model inadequacy associated with silver (Ag) release fractions using the AGR-2 Compact 2-1-2 as a case study. The results indicate that uncertainties in the operating temperatures supplied to the model and the form of those temperature conditions (spatially averaged versus spatially resolved) both affect Ag release predictions. This suggests that previously observed model inadequacy can be attributed to both physical models and the operating conditions supplied to them. These findings will be leveraged to guide development of assessment cases for additional experiments and refinement of physical models for behaviors that potentially influence Ag release, such as palladium penetration.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A sensitivity analysis to predict the neutronics behavior of samples irradiated in the VTR rabbit system

We report a low-order neutronics model is developed to carry out hundreds of simulations efficiently and investigate the neutronics behavior of samples being irradiated in a test reactor setting under different geometrical constraints. The low-order model allowed for simulations that yield the expected neutronics behavior of any irradiated sample in any environment and allows for the calculation of highly accurate spatially averaged statistics and idealized spatial distributions in the neutron flux. Several benchmarks are performed to evaluate the performance and limitations of the low-order model revealing many important findings. The low-order model predicted the LHGR in the EBR-II driver fuel to within 2.34% by only simulating the fuel rod by itself, which served as a validation for the model. Sensitivity studies investigated 3% enriched UO 2 and U-10Zr being irradiated in the Versatile Test Reactor rabbit system. The analyses investigated a range of combinations of 15 radii and 5 heights for each sample in the rabbit system. Similar data sets are also provided for irradiations in the Advanced Test Reactor’s B-10 irradiation position, which is a thermal neutron spectrum environment. Generalized fits and fit coefficients are obtained for sample heating, reaction rate densities, and local multiplication rate characteristics, allowing the predictions of the neutronics behavior of the samples based on their geometrical constraints. The analyses and fits laid the groundwork for developing a user-end Multiphysics analysis framework to assist and accelerate irradiation experiment design and optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗