OpenMC – это современный код моделирования процессов переноса нейтронов методом Монте-Карло, который использует в качестве программного интерфейса (API) язык программирования Python. OpenMC поддерживает восемь алгоритмов моделирования выгорания. В этом исследовании представлены результаты выбора метода интегрирования при моделировании выгорания ТВС с выгорающими поглотителями для реакторов ВВЭР-1000. Результаты моделирования выгорания, полученные по OpenMC, сравнивались результатами, представленными в бенчмарк OECD. Для моделирования выгорания в коде OpenMC можно использовать 8 различных численных интеграторов: PI, CE/CM, LE/QI, CE/LI, CF4, EPC-RK4, SI-CE/LI, SI-LE/QI. Результаты тестов показали, что интеграторам SI-CE/LI, SI-LE/QI требуется значительно больше времени для расчета одного шага по выгоранию, чем остальным при одинаковой точности, поэтому они были исключены из дальнейшего рассмотрения. Интегратор PI показал низкую точность интегрирования при одинаковых шагах по выгоранию с другими интеграторами. Однако PI обладает высоким быстродействием, в сравнении с другими интеграторами, и по мере уменьшения шага интегрирования обладает сходимостью к одному решению, которое может быть выбрано как реперное для оценки качества других интеграторов. На основе результатов, полученных с использованием интегратора PI с мелким шагом, было принято решение использовать для дальнейшей работы интегратор CE/LI. Результаты, полученные с помощью CE/LI, были сравнены с результатами, полученными в бенчмарке VVER-1000 LEU and MOX по кодам: MCU, TVS-M, WIMS8A, HELIOS, MULTICELL, и показали хорошее совпадение. Таким образом можно сделать вывод применимости интегратора CE/LI в составе OpenMC для моделирования выгорания ТВС, содержащих выгорающие поглотители. При проведении работ были использованы ресурсы высокопроизводительного вычислительного центра НИЯУ МИФИ
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.
The authors declare that there are no conflicts of interest present.