Numerical Simulations of Subcritical Reactor Kinetics in Thermal Hydraulic Transient Phases

A subcritical reactor driven by a linear proton accelerator has been considered as a nuclear waste incinerator at Korea Atomic Energy Research Institute(KAERI). Since the multiplication factor of a subcritical reactor is less than unity, to compensate exponentially decreasing fission neutrons from spallation reactions are essentially required for operating the reactor in its steady state. furthermore, the profile of accelerator beam currents is very important in controlling a subcritical reactor, because the reactor power varies in accordance of the profile of external neutrons. We have developed a code system to find numerical solutions of reactor kinetics equations, which are the simplest dynamic model for controlling reactors. In a due course of our previous numerical study of point kinetics equations for critical reactors, however, we learned that the same code system can be used in studying dynamic behavior of the subcritical reactor. Our major motivation of this paper is to investigate responses of subcritical reactors for small changes in thermal hydraulic parameters. Building a thermal hydraulic model for the subcritical reactor dynamics, we performed numerical simulations for dynamic responses of the reactor based on point kinetics equations with a source term. Linearizing a set of coupled differential equations for reactor responses, we focus our research interest on dynamic responses of the reactor to variations of the thermal hydraulic parameters in transient phases.