Initial Testing of the Microscopic Depletion Implementation in the MAMMOTH Reactor Physics Application
Present and new nuclear fuels that will be tested at the Transient Reactor Test (TREAT) facility will be analyzed with the MAMMOTH reactor physics application, currently under development, at Idaho National Laboratory. MAMMOTH natively couples the BISON, RELAP-7, and Rattlesnake applications within the MOOSE framework. This system allows the irradiation of fuel from beginning of life in a nuclear reactor until it is placed in TREAT for fuel testing within the same analysis mesh and, thus, retaining a very high level of resolution and fidelity. The calculation of the isotopic distribution in fuel requires the solution to the decay and transmutation equations coupled to the neutron transport equation. The Chebyshev Rational Approximation Method (CRAM) is the current state-of-the-art in the field, as was chosen to be the solver for the decay and transmutation equations. This report shows that the implementation of the CRAM solver within MAMMOTH is correct with various analytic benchmarks for decay and transmutation of nuclides. The results indicate that the solutions with CRAM order 16 achieve the level of precision of the benchmark. The CRAM solutions show little sensitivity to the time step size and consistently produce a high level of accuracy for isotopic decay for time steps of 1x10^11 years. Comparisons to DRAGON5 with 297 isotopes yield comparable results, but some differences need to be further analyzed.