A quasi-static Monte Carlo algorithm for the simulation of sub-prompt critical transients

Abstract

A pseudo-time-dependent Monte Carlo algorithm suitable for simulating thermal reactor transients with timescales on the order of seconds is presented. An initial experimental implementation has been written into the Monte Carlo Application ToolKit (MCATK) developed at Los Alamos National Laboratory. The algorithm is essentially a combination of the fully time-dependent Monte Carlo algorithm, including the explicit simulation of delayed neutron precursors, and a quasi-static algorithm that does not use the point-kinetics equations. The fully time-dependent calculation is performed for some number of inactive and active micro time steps. The dynamic reactor time constant and delayed neutron precursor production rates are calculated directly from the active cycles, rather than from point-kinetics parameters. These quantities are then used to advance the neutron population and delayed neutron precursors forward in time across a macro time step. The method is applied to the C5G7-TD benchmark operating in steady state and in transient mode.