Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code-to-code benchmark”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Assessment of the CTF subchannel code for modeling a large-break loss-of-coolant accident reflood transient

With increased industry interest in extending reactor operating cycles, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has been investigating the behavior of high-burnup fuel during design basis accidents such as the large-break loss-of-coolant accident (LBLOCA) with consideration for risk of fuel fragmentation, relocation, and dispersal (FFRD). As part of that activity, the NEAMS subchannel thermal/ hydraulics (T/H) code, CTF, is being used for modeling of LBLOCA and to determine the impact of subchannel resolution on results. Although CTF includes a wide range of models for LBLOCA conditions, the code has not been used for this application while maintained at Oak Ridge National Laboratory (ORNL) until now. Therefore, here, in this work, a preliminary assessment of several of these models was performed using openly available reflood experimental data from the Flooding Experiments in Blocked Arrays (FEBA) tests. One coarse mesh and one fine mesh model were set up in CTF for high and low flooding rate tests performed in the unblocked FEBA facility. A coarse TRACE model was set up to be as consistent as possible with the coarse CTF model to allow for code-to-code benchmarking. The assessment shows a tendency of the codes to over-predict peak cladding temperature (PCT) near the top of the bundle and to quench early. Advanced spacer grid models were shown to improve upper bundle predictions in CTF. The resolved CTF model over-predicted PCT by a larger degree in the center channels in the low-flooding rate test, and it is believed that the radiative heat transfer model, which was not used in this study, may be needed to correct this over-prediction. Finally, this work demonstrates the importance of the droplet model in determining quench time and vapor temperature and PCT prediction, which necessitates a more in-depth validation of these models in the future.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

SMR safety through HTTF modeling and benchmark efforts for code validation for gas-cooled reactor applications

Accurate modeling and simulation tools for thermal-hydraulics calculations are a key element needed to design and license new advanced reactors including Small Modular Reactors (SMR) and Microreactors. Uncertainties in modeling and simulation can have significant safety and economic implications. The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled prismatic-block nuclear reactors. High-quality measurement data is available from the HTTF that is suitable for a thermal-hydraulics code validation benchmark for gas-cooled reactor simulations. Here, this paper summarizes individual HTTF modeling efforts to date for tool validation at Idaho National Laboratory (INL), Argonne National Laboratory (ANL), Oregon State University (OSU) and Canadian Nuclear Laboratories (CNL) using system thermal-hydraulics codes, Computational Fluid Dynamics (CFD) codes and system-CFD code couplings. Also, the paper introduces the ongoing OECD Nuclear Energy Agency (NEA) High Temperature Gas Reactor Thermal-Hydraulics (HTGR T/H) benchmark that allows for better comparisons of results between different international modeling teams. The benchmark provides well defined computational problems that include code-to-code comparisons and comparisons to measured data. These problems provide an avenue for quantifying accuracy and identifying sources of uncertainty in thermal-hydraulics calculations, including in measured thermophysical properties, as part of validation for gas-cooled reactor simulation tools.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Second-order wave excitation forces in WEC-Sim/MOST: Implementation, experimental validation, and code-to-code comparison

Accurate prediction of second-order hydrodynamic loads is essential for floating bodies, including floating offshore wind turbines, wave energy converters, and hybrid wind–wave platforms. These nonlinear effects, arising from both sum- and difference-frequency forcing, are critical for capturing key response characteristics but remain challenging to model efficiently. In this work, we extend the open-source Wave Energy Converter Simulator / MATLAB for Offshore Simulation Tool by implementing second-order wave excitation forces, supporting both the full Quadratic Transfer Function formulation and the Newman approximation. The full Quadratic Transfer Function method is used for all code-to-code comparisons and experimental validation, while the Newman approximation is provided as a computationally lighter alternative. To benchmark the new capability, we perform a code-to-code comparison with OpenFAST and OrcaFlex. We then validate the enhanced model using wave-tank measurements of a 1:96 scale DeepCwind semi-submersible, showing that second-order effects are required to reproduce platform motions. The implementation employs a computationally efficient pre-computation strategy for second-order wave excitation forces, reducing simulation cost while maintaining engineering accuracy. Overall, this work advances the tool as an open-source and versatile tool for modelling floating offshore renewable-energy systems requiring second-order hydrodynamic fidelity.

17 WIND ENERGY↗
DOE OSTI2024

Feasibility and performance of the staged Z-pinch: A one-dimensional study with FLASH and MACH2

Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (>300) but achieve fuel ion temperatures of many kilo-electron volts. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Code Benchmark of the HTTF Pressurized Conduction Cooldown Test Using SAM

The High Temperature Test Facility (HTTF) at Oregon State University is an integral system test facility to simulate postulated reactor transients of prismatic high-temperature gas-cooled reactors(HTGRs). A series of test campaigns was launched, providing abundant test data that could be used to benchmark reactor system analysis codes like the System Analysis Module (SAM). In this study, a SAM model of the facility is developed based on the two-dimensional (2D) ring model approach. All components including the ceramic matrix, graphite heaters, helium coolant channels, core barrel, upcomer, pressure vessel, and reactor cavity cooling system are modeled as concentric cylindrical rings. The model is used to simulate one of the benchmark problems-Pressurized Conduction Cooldown (PCC)-within the scope of the Organisation for Economic Co-operation and Development Nuclear Energy Agency International HTTF Benchmark. The simulations consist of two parts. In the first part, operating and boundary conditions as well as thermophysical properties of materials are specified for the benchmark problem. In this work, results from the first part will be used in code-to-code comparison. In the second part, the SAM model is used to simulate Test PG-27, which is the first PCC test carried out in the HTTF, with only two of the ten heater banks activated. The results in the second part are used for code-to-data comparison. Because the helium coolant flow rate is not measured in this facility, it is estimated using the input power and inlet/outlet coolant temperatures. Additionally, radial heat flow in the ceramic blocks is complicated by hundreds of cylindrical coolant channels and heater rods embedded in them. As such, it is necessary to deduce an effective thermal conductivity for the ceramic to analyze the core thermal behavior. SAM predictions of the helium coolant and ceramic temperatures are compared with test data measured in three equivalent sectors. Overall, the SAM results agree reasonably well with test data within the variation of data among the three sectors, which demonstrates SAM's capability in capturing transient effects in HTGR using the simplified 2D ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Development and validation of fully open-source R2S shutdown dose rate capabilities in OpenMC *

We present the first fully open-source capabilities for shutdown dose rate (SDR) calculations of fusion energy facilities based on the Rigorous 2-Step (R2S) methodology. These capabilities have been implemented in the OpenMC Monte Carlo particle transport code, building on its existing capabilities while also leveraging new features that have been added to the code to support SDR calculations, such as decay photon source generation. Each of the individual physics components in the R2S workflow—neutron transport, activation, decay photon source generation, and photon transport—have been verified through code-to-code comparisons with MCNP6.2 and FISPACT-II 4.0. These comparisons generally demonstrate excellent agreement between codes for each of the physics components. The full cell-based R2S workflow was validated by performing a simulation of the first experimental campaign from the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For short cooling times, the dose calculated by OpenMC agrees with the experimental measurements within the stated experimental uncertainties. For longer cooling times, an overprediction of the shutdown dose was observed relative to experiment, which is consistent with previous studies in the literature. Altogether, these features constitute a combination of capabilities in a single, open-source codebase to provide the fusion community with a readily-accessible option for SDR calculations and a platform for rapidly analyzing the performance of fusion technology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A code-to-code benchmark for magneto-convection in a horizontal duct

Liquid metals and magnetic fields are used in many technical applications such as metallurgy, crystal growth and nuclear fusion reactors. When an electrically conducting fluid moves in a magnetic environment, electric currents and electromagnetic forces are generated that affect velocity and pressure losses in the flow. These magnetohydrodynamic (MHD) interactions have to be investigated to optimize the engineering processes. The characteristics of MHD flows depend on the geometrical configuration, the strength of the applied magnetic field, the electrical properties of fluid and structural materials and the thermal conditions. In the so-called blankets for fusion reactors, where liquid metals are used to breed the plasma fuel component tritium and to extract the generated heat, magneto-convective flows play a crucial role in determining heat and mass transfer. Therefore, the availability of numerical codes to simulate this type of flow is mandatory and their validation is a necessary step to guarantee the reliability of the results. For that reason, a benchmark problem has been defined to simulate liquid metal flows in a horizontal rectangular duct heated from below and exposed to a non-uniform magnetic field. Results obtained by five research groups using different codes are compared.

benchmark↗
DOE OSTI2024

Wind Turbine Aeroelastic Stability in OpenFAST

Wind turbines are growing in size and increasingly suffer from aeroelastic instabilities. Unfortunately, numerical models often show inconsistent results during verification studies. We address this gap by first introducing novel linearization capabilities within the open-source aero-hydro-servo-elastic framework OpenFAST. Next, a code-to-code benchmark study is presented that compares modal parameters between OpenFAST and HAWCStab2 for a land-based version of the International Energy Agency 15-MW reference wind turbine modeled with quasi-steady aerodynamics. The two solvers are in strong agreement except for discrepancies in the second rotor flapwise modes. The differences are attributed to the torsional flexibility of the tower, which is assumed torsionally stiff in the OpenFAST model. Work is ongoing to close this modeling gap. The aeroelastic stability of a low-specific-power land-based wind turbine is also investigated. The impact of design choices is discussed, high-lighting how narrow the margins are between a stable design and an unstable design.

17 WIND ENERGY↗
DOE OSTI2022

Cell Dancoff-based embedded self-shielding capability for doubly heterogeneous particulate fuels in SCALE/Polaris

Polaris is a two-dimensional (2D) transport lattice physics capability in the SCALE code system for light-water reactor (LWR) analysis. SCALE/Polaris is being extended to treat doubly heterogeneous particulate fuels in LWRs with accident-tolerant fuels and prismatic high-temperature gas-cooled reactor (HTGR) fuel. Recently, a Dancoff-based Wigner-Seitz approximation capability based on the embedded self-shielding method (DWA-ESSM) was implemented into Polaris for efficient computation. A new double heterogeneity treatment capability based on DWA-ESSM in conjunction with the Hebert's collision probability method for double heterogeneity was developed and implemented for Polaris. The new capability was verified by performing benchmark calculations for the prismatic HTGR fuel compact problems with various design parameters through a code-to-code comparison between the Polaris and continuous-energy Monte Carlo results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 1: Pin and lattice

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. This paper focuses on analysis of the predictions for pin and lattice designs representative of HB/HE core reloads from the Westinghouse advanced lattice code PARAGON2, the Monte-Carlo code SERPENT2 and MPACT, the deterministic neutronic component of VERA. A companion paper presented at this conference analyzes results obtained with the Westinghouse ANC9 core simulator, using lattice data provided by PARAGON2 with the NEXUS cross-section methodology, and MPACT for 2D and 3D core configurations incorporating HB/HE fuel reloads to achieve cycle length extension to 24-month in high-power density PWRs. The results of the pin and lattice benchmark presented in this paper show that the neutronic predictions from PARAGON2 for HB/HE fueled PWRs are in remarkable agreement with predictions from both SERPENT2 and MPACT including for depletion to high fuel burnup values, and employing a set of challenging lattices which feature combined use of burnable absorbers to obtain the required reactivity hold-down to support PWR transition to HB/HE fuel and 24-month operational cycles. As corroborated by this analysis and the analysis contained in the companion paper, the Westinghouse NEXUS-based PWR core analysis package can be confidently applied by industry for design analysis of PWR core reloads adopting HB/HE fuel and properly support the industry commercialization plans for in-reactor deployment. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 2: 2D and 3D core simulations

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. The benchmark consisted of a progression in complexity from modeling single 2D pin cells, through 2D lattices and finally to whole core 2D and then 3D scenarios. The accompanying paper (Part 1) presented the results for a variety of pin cell and lattices representative of high energy cores and established the ability of the MPACT code to accurately match the predictions of the Monte-Carlo code SERPENT2 such that the former may be used as a surrogate reference code in spatially complex whole core modeling applications where using Monte-Carlo methods is somewhat impractical. This paper focuses on comparisons at the whole core level using a representative high power density 3-loop Westinghouse NSSS core operating for 24 month cycles employing high enriched fuel with use of combined burnable absorbers to hold down excess reactivity. The evolution of the core reactivity and power distribution characteristics throughout the cycle are calculated using the Westinghouse NEXUS-based code suite and compared to the results using MPACT. The results of this whole core benchmark confirm the conclusions that the Westinghouse code suite accurately models HB/HE fuel designs, with very good prediction of core reactivity and power distribution throughout the cycle. The industry can be confident that the NEXUS code system can be deployed to support PWR reload safety and operational analyses for cores featuring HB/HE fuel and 24-month cycle operation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Benchmark for Fuel Shuffling and Depletion for Pebble-Bed Reactors

Pebble bed reactors have specific operational characteristics when their fuel-cycle and fueling operations are considered. They are specifically distinguished by other type of nuclear reactor designs by their online fuel recycling scheme, where the fuel elements that have not yet reached discharge burnup can be reloaded and recycled continuously during normal operation. The fuel in a pebble bed reactor is not stationary and stochastically moves through the core once or several times during its lifetime, which allows them to operate without requiring a large excess reactivity hold for the burnup. However, this characteristic of pebble bed reactors introduces challenges in simulation, as each pebble can take many different trajectories through the core, its composition depends on the details of the irradiation history that is unique to its aggregated path through the core. For predicting the safety performance characteristics, such as source term, maximum fuel temperatures and fuel failure rates, etc., it is important to accurately incorporate the movement of pebbles through the core during their lifetime in a multi-physics simulation together with other phenomena. The equilibrium core analysis for pebble bed reactors are performed with multi-physics tools including fuel depletion in a multi pass reload coupled to the fuel movement. Currently, there are only a few legacy multi-physics simulation tools that can implement the pebble flow characteristics and perform equilibrium core analysis for pebble bed reactors. However, there are development efforts on-going under Department of Energy's Nuclear Energy Advanced Modelling and Simulation program and also in private industry for including these capabilities into their modelling and simulation tools. Any new development in the modelling and simulation tools needs to be validated by using tools such as experiments, analytical solutions or code-to-code benchmarks. In this work, a code-to-code benchmark for the equilibrium core analysis capability of pebble bed reactors was developed. Multiple cases were identified to capture different fuel cycle strategies that can be used in PBRs. The results of each case are presented in terms of overall equilibrium core characteristics: the discharge burnup; spatial burnup distribution; spatial isotopic distributions; axial and radial neutron flux distributions and power history of fuel elements per pass through the core for both a prototypical pebble bed High Temperature Gas-cooled Reactor and a prototypical pebble bed Fluoride-salt cooled High temperature Reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Isotopic and Fuel Lattice Parameter Trends in Extended Enrichment and Higher Burnup LWR Fuel Vol I: PWR fuel

Commercial light water reactor (LWR) operators and fuel vendors in the United States are pursuing changes to nuclear fuel that include extended enrichment (EE) and accident-tolerant fuel (ATF) designs. The term EE (8% > 235 U > 5%) is used in this report to refer to a subset of high assay low-enriched uranium (HALEU) that is considered usable in commercial US LWRs in the near term. ATF features are designed to improve fuel system performance under accident conditions. One goal of EE is to improve fuel cycle economy by enabling fuel to be depleted to higher burnup than the typical current maximum pin burnup limits (62 gigawatt-days per metric ton of uranium [GWd/MTU]). Adoption of EE, ATF, and high burnup (HBU) fuels in the US commercial fleet requires a clear understanding of the effects on core physics parameters and used fuel isotopic content, as well as confidence in the accuracy of computer code predictions over an expanded range of materials, enrichment, and burnup. A thorough understanding of the applicability and adequacy of benchmark data (e.g., criticality, decay heat, isotopic content) for computer code validation is necessary to ensure that appropriate safety margins are maintained. To prepare for and support these potential changes, the effects of EE, ATF, and HBU are being assessed for selected representative LWR fuel designs. The project is divided into phases: this report summarizes the findings of Phase 1, which focuses on the lattice physics parameter and used fuel isotopic changes for a conventional Westinghouse 17×17 pressurized water reactor (PWR) design. The primary investigation tool is the SCALE Polaris code using the SCALE 56-group Evaluated Nuclear Data File (ENDF)/B-VII.1 cross sections. The goal of the current work is to (1) identify and explain important effects of EE and HBU (reactivity, lattice physics, and isotopic effects) assuming that PWR fuel design and usage remain similar to those for current enrichment fuel, (2) provide limited code-to-code comparisons with higher order cross section libraries and/or codes, and (3) identify any apparent anomalous trends in the results for further investigation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Coupled Neutronic and Thermal Hydraulic Analysis of a Natural Circulation Based Small Modular Reactor (SMR) Using VERA-CS

As part of the work supported by a US Department of Energy (DOE) Office of Nuclear Energy Gateway for Accelerated Innovation in Nuclear FY 2020 Voucher, Holtec International subsidiary SMR LLC (hereinafter referred to as Holtec) and Oak Ridge National Laboratory entered into a cooperative research and development agreement (CRADA) to develop coupled multiphysics core models of the Holtec developed small modular reactor (SMR), SMR-160. The scope of the 1 year project was to use the DOE–developed tool set, VERA, to analyze several aspects of a representative SMR-160 core design. The goals of the project were to perform the code-to-code benchmarking and to provide Holtec with a confirmatory calculations to compare against the results with the codes used in the SMR-160 design, identify and resolve areas of discrepancy between the two models to give Holtec a basis for selecting certain tool and model options, and develop a workflow so that the VERA model could be adapted easily to any future changes in the design of the SMR-160 core. The work proceeded in four phases: 1. Developing the initial VERA model based on Holtec-provided core design information; 2. Comparing VERA results with Holtec-provided calculations for the first five cycles; 3. Developing a VERA transient model based on the prior model; and 4. Training Holtec personnel on VERA and hand-off the developed models. This report summarizes at a high level the completed tasks. Details of the tasks and results are reserved for the full report, which will not be publicly released for 5 years per the terms of the CRADA. The report is titled “Coupled Neutronic and Thermal Hydraulic Analysis of a Natural Circulation Based Small Modular Reactor (SMR) using VERA-CS” and has the ORNL technical report number ORNL/TM 2021/2303 (NFE-20-08305)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗