IntroductionThe concept of using continuous energy Monte Carlo simulations for critical systems as a precise tool for burnup calculations is well known in the fi eld of research and development of fuel cycle. The beginning-of-step approximation of neutron fl ux over the time-step is often applied by many burnup codes for its simplicity. This approach assumes constant fl ux/power profi le during entire step, which can cause problems such as spatial instability of reaction rates in the simulated system [1]. Spatial oscillations of fl ux and xenon concentration appear for relatively short steps, thus affecting the core's equilibrium. Regarding fuel cycle analysis of prismatic HTGR, the spatial fl ux distribution of the core varies strongly in the vicinity of compensation rods, which are being slowly withdrawn in order to compensate the reactivity loss. In order to handle this problem properly, a better fl ux normalization procedure -so-called bridge scheme -was developed and applied in the study of PuMA project [2]. An optimal model of time-step is necessary to account for the non-constant system's behavior as well as to provide stability of burnup. Abstract. In this paper, we compare the methodology of different time-step models in the context of Monte Carlo burnup calculations for nuclear reactors. We discuss the differences between staircase step model, slope model, bridge scheme and stochastic implicit Euler method proposed in literature. We focus on the spatial stability of depletion procedure and put additional emphasis on the problem of normalization of neutron source strength. Considered methodology has been implemented in our continuous energy Monte Carlo burnup code (MCB5). The burnup simulations have been performed using the simplifi ed high temperature gas-cooled reactor (HTGR) system with and without modeling of control rod withdrawal. Useful conclusions have been formulated on the basis of results.