Effect of Delayed Energy Release on Power Normalization in Reactor Depletion Calculations
Reactor depletion calculations typically assume that the thermal power of the reactor is directly proportional to the instantaneous fission rate within the core. However, approximately 7% of the sensible energy released per fission comes from the radioactive decay of fission products, which occurs alter the fission event. Thus, the assumed direct relationship between reactor power and fission rate implies that all fission products have reached secular equilibrium within the core. In situations where fission product concentrations have not reached equilibrium (e.g., shortly after reactor startup or a change in reactor power during operation), the delayed energy release due to fission product decay creates a nonlinear relationship between reactor power and fission rate. This non-linear relationship, in turn, affects the neutron flux levels and fuel consumption rate required to maintain the desired thermal power output for the reactor. For example, during initial startup of a commercial light water reactor the apparent energy release per fission is ∼7{sup %} lower than the typical value of ∼200 MeV/fission due to the potential energy stored in radioactive fission products. As a consequence, the flux level (and, therefore, uranium consumption rate) in the reactor must be correspondingly higher to preserve the desired reactor power. As the reactor operates and fission product inventories increase, the decay heat produced by the reactor will rise, producing a corresponding increase in the apparent energy release per fission. The characteristic time-dependent behavior of energy release due to the decay of fission products (i.e., decay heat) following reactor startup, shutdown, or changes in reactor power is well-known and widely-used in reactor safety analyses. However, to our knowledge, the effects of time-dependent decay heating rates on reactor depletion calculations have not been thoroughly investigated. In this paper, we present a simple methodology for explicitly accounting for delayed energy release due to the decay of fission products in reactor depletion calculations. In the proposed method, the balance equations for fission product decay energy are incorporated into the depletion coefficient matrix, using a strategy similar to that described by Isotalo in 2016. The modified system of differential equations can be solved via a matrix exponential approximation (or other suitable numerical method) to yield solutions for both nuclide densities and decay energy emission densities as a function of time. Numerical results from the method are presented for simplified (one-group point (0-D) and 1-D) representations of a commercial pressurized water reactor (PWR), in order to estimate the effect of delayed energy release on fuel consumption and flux level immediately following reactor startup.