American Nuclear Society

Home / Publications / Journals / Nuclear Science and Engineering / Volume 162 / Number 3

Numerical Stability of Existing Monte Carlo Burnup Codes in Cycle Calculations of Critical Reactors

Jan Dufek, J. Eduard Hoogenboom

Nuclear Science and Engineering / Volume 162 / Number 3 / July 2009 / Pages 307-311

Technical Note /

We show that major existing Monte Carlo burnup codes are numerically unstable in cycle calculations of critical reactors; spatial oscillations of the neutron flux can be observed even when relatively small time steps are used. This is caused by using the explicit Euler or midpoint method that appear to be numerically unstable with the step sizes common in cycle calculations. More stable methods that are common in deterministic burnup calculations, like the modified Euler method, can easily be introduced into the Monte Carlo burnup codes.

Questions or comments about the site? Contact the ANS Webmaster.