EPJ Web Conf.
Volume 247, 2021PHYSOR2020 – International Conference on Physics of Reactors: Transition to a Scalable Nuclear Future
|Number of page(s)||8|
|Section||Verification & Validation|
|Published online||22 February 2021|
CODE VERIFICATION OF MPACT USING GANAPOL CRITICAL ROD BENCHMARK
1 Department of Nuclear Engineering and Radiological Sciences, University of Michigan 2355 Bonisteel Blvd, Ann Arbor, MI 48109, USA
2 Department of Nuclear Engineering and Radiological Sciences, University of Michigan 2355 Bonisteel Blvd, Ann Arbor, MI 48109, USA
3 Oak Ridge National Laboratory 1 Bethel Valley Road, Oak Ridge, TN 37830, USA
Published online: 22 February 2021
This work is dedicated to the code verification of MPACT, which is developed under the Consortium for Advanced Simulation of Light Water Reactors by the University of Michigan and Oak Ridge National Laboratory, where the numerical solution is compared to the reference solution of a benchmark problem with a known analytical solution. In this work, Benchmark Problem 3.4 in Barry Ganapol’s benchmark book was chosen as an MOC code verification test problem. Problem 3.4 is a bare cylinder of infinite height, which is an excellent benchmark problem for 2D MOC. To ensure that this benchmark problem exercised the same code as typically used by MPACT, the bare rod configuration was surrounded by a bounding box filled with a non-scattering material. To avoid implementing a critical rod search in the MPACT code, the MPACT analysis was performed using cross sections that yielded the given c, the average number of secondary neutrons per collision, and a rod radius that was the corresponding critical rod radius. MPACT agreed with all cases to within a few pcm. The convergence behavior was studied. The results show a 2nd order radial convergence, consistent with flat-source approximation. The convergence curves with respect to ray spacing and polar angle quadrature set order were also obtained. The other quantity of interest tabulated for Problem 3.4 was the radial distribution of the scalar flux. Two configurations were analyzed, and the resultant radial flux profiles agreed very well with the tabulated results. The verification of the production neutronics code MPACT has been augmented by the addition of the analytical solutions for an infinite cylinder from the Ganapol benchmark book. These test cases can be included in the regression suite for MPACT.
Key words: Code verification / MPACT / Ganapol benchmark / critical rods / convergence
© The Authors, published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Initial download of the metrics may take a while.