Numerical Stability of Existing Monte Carlo Burnup Codes in Cycle Calculations of Critical Reactors
2009 (English)In: Nuclear science and engineering, ISSN 0029-5639, Vol. 162, no 3, 307-311 p.Article in journal (Refereed) Published
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.
Place, publisher, year, edition, pages
2009. Vol. 162, no 3, 307-311 p.
Burn up; Burnup calculation; Critical reactors; Midpoint method; Modified Euler method; MONTE CARLO; Numerical stabilities; Spatial oscillations; Step size; Time step
IdentifiersURN: urn:nbn:se:kth:diva-14050ISI: 000267621400008ScopusID: 2-s2.0-68649086454OAI: oai:DiVA.org:kth-14050DiVA: diva2:329376
QC 201007092010-07-092010-07-092010-07-16Bibliographically approved