MCS/TH1D analysis of VERA whole-core multi-cycle depletion problems

Abstract

This paper presents the verification and validation elements of the UNIST in-house Monte Carlo code, MCS, for the multi-cycle and multi-physics analyses of high-fidelity, large-scale commercial pressurized water reactors (PWRs). Analysis on the neutronic performance with thermal/hydraulic (T/H) feedback is the key to detecting the complex behavior of an operating nuclear power reactor. The MCS solutions with T/H feedback, TH1D, of the consortium for advanced simulation of light water reactors (CASL) virtual environment for reactor applications (VERA) core physics benchmark progression problems 6 and 7 showed excellent agreement in eigenvalues, temperature and power profiles with the MC21/COBRA-IE, MC21/CTF and VERA-CS solutions for the single assembly and whole core of Watts Bar Nuclear 1 (WBN1) Cycle 1 under beginning-of-cycle and hot-full-power condition. Furthermore, the core depletion analysis is one of the most compelling advances for reactor analysis. Therefore, this work is significantly focused on the nuclide depletion simulation coupled with TH1D of the first two cycles of WBN1 to address the VERA core physics benchmark problems 9 and 10. MCS is one of the few Monte Carlo codes that have the capability of depletion calculation for both WBN1 Cycles 1 and 2. The accuracy of MCS simulation of WBN1 Cycle 1 is within 40 ppm and 30 ppm in critical boron concentration (CBC) for all burnup points, compared with the measured data and VERA-CS solutions, respectively. To demonstrate the multi-cycle refueling capability of MCS, the WBN1 Cycle 2 is simulated and compared with the solutions of VERA-CS only, because of the lack of available measured data. MCS shows excellent agreement compared with VERA-CS within 30 ppm in CBC, and the average bias for the entire Cycle 2 is approximately 20 ppm. These results provide confidence in MCS’s capability in high-fidelity, multi-cycle calculations of the practical PWR core.