Improved Reaction Rate Tracking and Fission Product Yield Determinations for the Monte Carlo-Linked Depletion Capability in MCNPX

Abstract

As advanced reactor concepts challenge the accuracy of current modeling technologies, a higher-fidelity depletion calculation is necessary to model time-dependent core reactivity properly for accurate cycle length and safety margin determinations. The recent integration of CINDER90 into the MCNPX Monte Carlo radiation transport code provides a completely self-contained Monte Carlo-linked depletion capability. Two advances have been made in the latest MCNPX capability based on problems observed in prereleased versions: continuous-energy collision density tracking and adequate fission yield selection.

Prereleased versions of the MCNPX depletion code calculated the reaction rates for (n,2n), (n,3n), (n,p), and (n,α) by matching the MCNPX steady-state 63-group flux with 63-group cross sections inherent in the CINDER90 library and then collapsing to one-group collision densities for the depletion calculation. The accuracy of this procedure is therefore dictated by the adequacy of the 63-group energy structure of the cross-section set to accurately capture the spectrum of a specific model. Different types of models would therefore require different types of cross-section energy group structure. MCNPX 2.6.A eliminates this dependency by using the continuous-energy reaction rates determined during the MCNPX steady-state calculation to calculate energy-integrated collision rates to be used by CINDER90.

MCNPX 2.6.A now also determines the proper fission yield to be used by the CINDER90 code for the depletion calculation. The CINDER90 code offers a thermal, fast, and high-energy fission yield for each fissile isotope contained in the CINDER90 data file. MCNPX 2.6.A determines which fission yield to use for a specified problem by calculating the integral fission rate for the defined energy boundaries (thermal, fast, and high energy), determining which energy range contains the majority of fissions, and then selecting the appropriate fission yield for the energy range containing the majority of fissions.

The MCNPX depletion capability enables complete, relatively easy-to-use depletion calculations in a single Monte Carlo code. This study focuses on the methodology development of the two improvements described here. Further improvements are under development to enhance the usefulness of this new capability.