Benjamin R. Betzler

High-Fidelity Modeling and Simulation for a High Flux Isotope Reactor Low-Enriched Uranium Core Design

Abstract A high-fidelity model of the High Flux Isotope Reactor (HFIR) with a low-enriched uranium (LEU) fuel design and a representative experiment loading has been developed to serve as a new reference model for LEU conversion studies. With the exception of the fuel elements, this HFIR LEU model is completely consistent with the current highly enriched uranium HFIR model. Results obtained with the new LEU model provide a baseline for analysis of alternate LEU fuel designs and further optimization studies. The newly developed HFIR LEU model has an explicit representation of the HFIR-specific involute fuel plate geometry, including the within-plate fuel meat contouring, and a detailed geometry model of the fuel element side plates. Such high-fidelity models are necessary to accurately account for the self-shielding from 238U and the depletion of absorber materials present in the side plates. In addition, a method was developed to account for fuel swelling in the high-density LEU fuel plates during the depletion simulation. Calculated time-dependent metrics for the HFIR LEU model include fission rate and cumulative fission density distributions, flux and reaction rates for relevant experiment locations, point kinetics data, and reactivity coefficients.

Two-Dimensional Neutronic and Fuel Cycle Analysis of the Transatomic Power Molten Salt Reactor

This status report presents the results from the first phase of the collaboration between Transatomic Power Corporation (TAP) and Oak Ridge National Laboratory (ORNL) to provide neutronic and fuel cycle analysis of the TAP core design through the Department of Energy Gateway for Accelerated Innovation in Nuclear, Nuclear Energy Voucher program. The TAP design is a molten salt reactor using movable moderator rods to shift the neutron spectrum in the core from mostly epithermal at beginning of life to thermal at end of life. Additional developments in the ChemTriton modeling and simulation tool provide the critical moderator-to-fuel ratio searches and time-dependent parameters necessary to simulate the continuously changing physics in this complex system. Results from simulations with these tools show agreement with TAP-calculated performance metrics for core lifetime, discharge burnup, and salt volume fraction, verifying the viability of reducing actinide waste production with this design. Additional analyses of time step sizes, mass feed rates and enrichments, and isotopic removals provide additional information to make informed design decisions. This work further demonstrates capabilities of ORNL modeling and simulation tools for analysis of molten salt reactor designs and strongly positions this effort for the upcoming three-dimensional core analysis.

Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method

Abstract For a nuclear system in which the entire -eigenvalue spectrum is known, eigenfunction expansion yields the time-dependent flux response to any arbitrary source. Applications in which this response is of interest include pulsed-neutron experiments, accelerator-driven subcritical systems, and fast burst reactors, where a steady-state assumption used in neutron transport is invalid for characterizing the time-dependent flux. To obtain the -eigenvalue spectrum, the transition rate matrix method (TRMM) tallies transition rates describing neutron behavior in a discretized position-direction-energy phase space using Monte Carlo. Interpretation of the resulting Markov process transition rate matrix as the operator in the adjoint -eigenvalue problem provides an avenue for determining a large finite set of eigenvalues and eigenfunctions of a nuclear system. Results from the TRMM are verified using analytic solutions, time-dependent Monte Carlo simulations, and modal expansion from diffusion theory. For simplified infinite-medium and one-dimensional geometries, the TRMM accurately calculates eigenvalues, eigenfunctions, and eigenfunction expansion solutions. Applications and comparisons to measurements are made for the small fast burst reactor CALIBAN and the Fort St. Vrain high-temperature gas-cooled reactor. For large three-dimensional geometries, discretization of the large position-energy-direction phase space limits the accuracy of eigenfunction expansion solutions using the TRMM, but it can still generate a fair estimate of the fundamental eigenvalue and eigenfunction. These results show that the TRMM generates an accurate estimate of a large number of eigenvalues. This is not possible with existing Monte Carlo–based methods.

Calculating alpha Eigenvalues in a Continuous-Energy Infinite Medium with Monte Carlo

The {alpha} eigenvalue has implications for time-dependent problems where the system is sub- or supercritical. We present methods and results from calculating the {alpha}-eigenvalue spectrum for a continuous-energy infinite medium with a simplified Monte Carlo transport code. We formulate the {alpha}-eigenvalue problem, detail the Monte Carlo code physics, and provide verification and results. We have a method for calculating the {alpha}-eigenvalue spectrum in a continuous-energy infinite-medium. The continuous-time Markov process described by the transition rate matrix provides a way of obtaining the {alpha}-eigenvalue spectrum and kinetic modes. These are useful for the approximation of the time dependence of the system.