Selecting Burnup Algorithms in OpenMC Using the Calculated Benchmark of LEU Assembly and MOX Fuel

OpenMC is a state-of-the-art Monte Carlo neutron transport simulation code that uses the Python programming language as an API. OpenMC supports eight burnout simulation algorithms. This study presents the results of choosing an integration method for modeling the burnup of fuel assemblies with burnable poisons for WWER-1000 reactors. Burnout simulation results from OpenMC were compared with those reported in the OECD benchmark. 8 different numerical integrators can be used to model burnout in OpenMC code: PI, CE/CM, LE/QI, CE/LI, CF4, EPC-RK4, SI-CE/LI, SI-LE/QI. The test results showed that the SI-CE/LI, SI-LE/QI integrators require significantly more time to calculate one burnup step than the others with the same accuracy, so they were excluded from further consideration. The PI integrator showed low integration accuracy at the same burnup steps with other integrators. However, PI has a high performance compared to other integrators, and as the integration step decreases, it converges to one solution, which can be chosen as a reference for assessing the quality of other integrators. Based on the results obtained using the fine step PI integrator, it was decided to use the CE/LI integrator for further work. The results obtained with CE/LI were compared with those obtained with the VVER-1000 LEU and MOX benchmark for codes: MCU, TVS-M, WIMS8A, HELIOS, MULTICELL and showed good agreement. Thus, we can conclude the applicability of the CE/LI integrator as part of OpenMC for modeling the burnup of fuel assemblies containing burnable poisons. During the work, the resources of the high-performance computer center of the National Research Nuclear University MEPhI were used.

1. Romano P.K., Horelik N.E., Herman B.R., Nelson A.G., Forget B., Smith, K. Openmc: A state-of-the-art monte carlo code for research and development. Annals of nuclear Energy, 2015, V. 82, рр. 90-97.

2. Dufek J., Kotlyar D., Shwageraus E. The stochastic implicit euler method–a stable coupling scheme for monte carlo burnup calculations. Annals of nuclear Energy, 2013, V. 60, рр. 295-300.

3. Iserles A., Munthe-Kaas H.Z., Nørsett S.P., Zanna A. Lie-group methods. Acta numerica. Acta Numer, 2000, V. 9, рр. 215-365.

4. Isotalo A.E., Aarnio P. Higher order methods for burnup calculations with bateman solutions. Annals of nuclear Energy, 2011, V. 38, рр. 1987-1995.

5. Josey. Development and analysis of high order neutron transport-depletion coupling algorithms, Ph.D. thesis, (Massachusetts Institute of Technology, 2017).

6. Богданович, Р.Б. Полная энергия деления в зависимости от глубины выгорания топлива для ВВЭР-1000 / Р.Б. Богданович Р.Б., А.С. Герасимов, Г.В. Тихомиров // Журнал физики: Серия конференции. – 2020. – № 1439.

7. Thilagam L., Sunil Sunny C., Jagannathan V., A VVER-1000 LEU and MOX assembly computational benchmark analysis using the lattice burnup code EXCEL.Annals of nuclear Energy, 2009, V. 36, рр. 505-519.

8. NEA/NSC/DOC. A VVER-1000 LEU and MOX Assembly Computational Benchmark. Nuclear Energy Agency, 2002.

9. Park, J., Khassenov, A., Kim, W., Choi, S., Lee, D. Comparative analysis of vera depletion benchmark through consistent code-to-code comparison.Annals of nuclear Energy, 2019, V. 124, рр. 385-398.

10. Muir, R. E. M. and D. W. The NJOY Nuclear Data Processing System, 1994.

11. Leppänen, J. Serpent–a continuous-energy monte carlo reactor physics burnup calculation code. VTT Technical Research Centre of Finland. Vol. 4.

12. Jiankai Yu a, B. F. Verification of depletion capability of OpenMC using VERA depletion benchmark. Annals of nuclear Energy, 2022, V. 170.

13. Девятко, Ю.Н. Моделирование распределения температуры и выгорания в уран-гадолиниевом твэле ВВЭР / Ю.Н. Девятко, В.В. Новиков, О.В. Хомяков // Физика атомного ядра. – 2018. – № 81 – С. 1257-1275.