A new Monte Carlo approach for solution of the time dependent neutron transport equation based on nodal discretization to simulate the xenon oscillation with feedback

Abstract

In this paper a probabilistic methodology based on core nodalization is proposed to estimate the core power in the presence of xenon oscillation. A time-dependent Monte Carlo neutron transport code named MCSP-NOD is developed for dynamic analysis in arbitrary 3D geometries to simulate xenon oscillations as well as sub-critical condition with feedbacks. The new code is based on the approach adopted in MCNP-NOD which was previously introduced as a tool for core transient analysis using the MCNPX platform. As before, the core is divided into nodes of arbitrary dimensions, and all terms of the transport equation e.g. interaction rates, leakage ratio are estimated using the MC techniques. However, as a new option the concentration of iodine and xenon are also estimated which enables us to predict the oscillatory behavior of reactor power following poison oscillation. These quantities are then employed within the time-dependent neutron transport equation for each node independently to compute the neutron population. Simulations prove the robustness of the method.