Abstract A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

Molten salt reactors (MSRs) are a class of next-generation nuclear reactors that have received recent industrial and research interest. A generalized species transport solver was implemented in the Virtual Environment for Reactor Applications (VERA) computing suite to extend this tool to analyze liquid fueled MSRs. This core simulator has been extended to model the transport of fission product gases into a collection of circulating gas bubbles with the purpose of removing the gases. This paper presents the governing species transport equation, along with various nuclear source terms. Development of the source term for phase migration is discussed, along with a simplified interfacial area tracking method. Finally, a case study on a simplified MSR loop is presented in which modeling parameters were varied to assess their impact on gas removal. The steady state results show that parameters such as bubble diameter, gas injection rate and mass transfer coefficient have a low to moderate effect on the fraction of xenon in the core region. Removal efficiency has the greatest effect on the fraction in the core region. After the pump bowl, bubble diameter has a minor effect on the fraction of xenon in the gas void. These results point out that increasing parameters such as mass transfer coefficient, gas injection rate, and removal efficiency drives the xenon into the circulating gas void, while decreasing bubble diameter also drives xenon into the gas void by increasing interfacial area. (c) 2021 Elsevier Ltd. All rights reserved.

approach. Furthermore, an ionic size dependent activity coefficient model is presented. This book also contains research into quaternary phase diagrams for LiCl-KCl-CsCl-RECl3 molten salt systems. This book should prove an excellent introduction to anyone pursuing a career in experimental high temperature electrochemistry and molten salt chemistry.

Matrix exponential methods have long been utilized for isotopic depletion in nuclear fuel calculations. In this paper we discuss the development of such methods in addition to species transport for liquid fueled molten salt reactors (MSRs). Conventional nuclear reactors work with fixed fuel assemblies in which fission products and fissile material do not transport throughout the core. Liquid fueled molten salt reactors work in a much different way, allowing for material to transport throughout the primary reactor loop. Because of this, fission product transport must be taken into account. The set of partial differential equations that apply are discretized into systems of first order ordinary differential equations (ODEs). The exact solution to the set of ODEs is herein being estimated using the matrix exponential method known as the Chebychev Rational Approximation Method (CRAM).

The objective of the integration roadmap is to describe the necessary physics required and a plan for modeling and simulation approach for predicting mass accountancy in molten salt reactors on an engineering scale. The thermophysical properties and the underlying thermodynamics are fundamental inputs. Therefore, the modeling will span length scales from first principles calculations to the engineering scale. The intention is to predict where material accumulates in a reactor core and loop and to understand perturbations on the systems level, for example the downstream effects from a turbine failure.

New multiphysics coupling capabilities for molten salt reactor (MSR) analysis have been developed in the Virtual Environment for Reactor Applications (VERA). This development consisted of two main efforts. First, a generic species transport module was added in the CTF code, which is the thermal-hydraulics (TH) code for VERA. This module uses the velocity fields for which CTF solves during the TH calculation to transport species through the core and around the primary loop. Additionally, a gas sparging model has been added to CTF to model the movement of certain species, namely, fission products such as xenon gas, to transport between the molten salt and gas bubbles present in the salt. The second effort in this development was coupling this capability to VERA's neutron transport code MPACT. This effort focused on coupling the detailed TH transport models in CTF to MPACT to account for feedback effects in the neutron transport calculations. Finally, the thermochemistry code Thermochimica has also been coupled to VERA. Thermochimica performs pointwise calculations for chemical potential and Gibbs free energy and determines what phases are produced by the temperature, pressure, and elemental concentrations at different locations in the primary loop.These capabilities are demonstrated using a model of the Molten Salt Reactor Experiment (MSRE). This reactor operated at Oak Ridge National Laboratory in the 1960s, providing sources of experimental data that were used to develop the model. Various combinations of species were modeled using VERA's new multiphysics coupling capabilities. Species distributions and reactivity effects behaved as expected for the MSRE model, demonstrating that the coupling is behaving correctly and causing appropriate feedback. The results of these calculations show the potential for VERA to be used for a wide variety of MSR analyses.

With the resurgence of interest in molten salt reactors, there is a need for new experiments and modeling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory to be placed in the Advanced Test Reactor (ATR). One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this has on the fuel-salt properties, source term generation, and corrosion control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system, and noble metals (Mo, Tc, Te, etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in CTF (Coolant-Boiling in Rod Arrays—Two Fluids) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design. A noble metal–helium bubble mass transfer model is coupled with the thermal-hydraulic results from CTF to determine the effectiveness of this insoluble fission product (IFP) extraction method for VESIL. Amounts of IFP species extracted to the off-gas system and species distributions in VESIL after a 60-day ATR cycle are reported.

Molten salt reactors (MSRs) are a family of advanced, generation-IV, nuclear reactors that use a circulating molten alkali fluoride fuel salt as both a primary working fluid and a fuel matrix. Since the fuel salt is liquid, gases become entrained in the fuel salt flow, and these entrained gases are called circulating voids. This paper reviews some aspects of circulating voids on reactor behavior. A formal definition of the circulating void fraction is given. A time-line of Oak Ridge reports on circulating void behavior is presented. Methods to measure the circulating void fraction are described. Foam induced by gas entrained in the fuel salt is discussed. The discovery of circulating voids in the molten salt reactor experiment (MSRE) is detailed. Aspects related to bubble size and interfacial area are reviewed. The Laplace pressure is examined. The void and pressure coefficients of reactivity are examined. Some requirements for a bubble model for xenon behavior analysis are proposed. The effects of bubble surface mobility are examined. Finally, sources of gas production in the fuel salt are investigated.

Xenon behaves differently in molten salt reactors (MSRs) compared to solid fuel reactors. This behavior needs exploring due to the large reactivity effect of the 135Xe isotope, given the current interest in MSR power plant development for commercial deployment. This paper focuses on select topics in xenon transport, reviews relevant past works, and proposes specific research questions to advance the state of the art in each of the focus areas. Specifically, the paper discusses the issue of xenon solubility in MSRs, the behavior of particulates circulating in MSR fuel salt and its influence on the xenon transport, the possibility of ionization of xenon atoms which changes its effective size and thus affects its mass transport, and finally the issue of circulating void fraction and how it is measured. This work presents specific recommendations for MSR designers to research the limits of Henry's law validity, circulating particulate scrubbers, validity of mass transport coefficients in high radiation fields, and the effects of pump speed on circulating void fraction.

This article discusses the aspects of Xe-135 behavior in molten salt reactor (MSR) graphite. Models of MSR graphite are described. The related mass transfer and mass diffusion coefficients are described and means by which they can be calculated are detailed. Xenon reactivity effects are explored. A method is presented to model the internal xenon distribution within the graphite stringers.

Molten salt reactor (MSR) designs require robust multiphysics modelling capabilities to further the development of this next generation technology. Due to the unique liquid fuel design, traditional reactor physics analysis codes do not capture all of the important physics of these reactors. This work builds upon a general species transport model implemented into the multiphysics core simulator suite VERA-CS, by adding source and sink models for insoluble fission products, specifically the noble metals. Noble metal (NM) fission products are born in the fuel-salt, but do not form stable fluorides, and instead plate-out on various surfaces in the reactor fuel loop. A concentration gradient advection-diffusion based mass transfer model is implemented into the general species transport model within CTF – the subchannel thermal hydraulics code in VERA-CS – and verified using convergence tests in space and time. Two noble metal decay chains are analyzed, and their steady state bulk liquid densities and surface concentrations are calculated for a simple flow loop in CTF that is roughly representative of the Oak Ridge National Laboratory Molten Salt Reactor Experiment (MSRE) [1]. Insoluble species 99Mo decaying to 99Tc captures wall accumulation of decaying species, and insoluble 131Sb to semi-soluble 131Te to soluble 131I captures wall loss of decaying species. Results are compared with MSRE data and the possible multiphysics effects of this phenomena are identified. Future work to enhance the robustness of the physical model and couple it with neutronics to evaluate the effects of NM mass transport is discussed.

This document analyzes molten salt reactor (MSR) xenon theory and reviews the corpus of literature pertaining to it. A history of xenon behavior in MSRs is presented. Relevant literature is tabulated. Assumptions in MSR xenon theories are made explicit. The structure of graphite stringers is examined. The behavior of bubbles is investigated along with interfacial area perturbation in response to thermodynamic state perturbation. This document serves as an underpinning to further investigation into the MSR xenon theory.