Democratic People’s Republic of Korea (DPRK) has produced weapon-grade plutonium in a graphite-moderated experimental reactor at the Yongbyon nuclear facilities. The amount of plutonium produced can be estimated using the Graphite Isotope Ratio Method (GIRM), even without considering specific operational histories. However, the result depends to some degree on the operational cycle length. Moreover, an optimal cycle length can maximize the number of nuclear weapons made from the plutonium produced. For conservatism, it should be assumed that the target reactor was operated with an optimal cycle length. This study investigated the optimal cycle length using which the Calder Hall MAGNOX reactor can achieve the maximum annual production of nuclear weapons. The results show that lower enrichment fuel produced a greater number of critical plutonium spheres with a shorter optimal cycle length. Specifically, depleted uranium (0.69wt%) produced 5.561 critical plutonium spheres annually with optimal cycle lengths of 251 effective full power days. This research is crucial for understanding DPRK’s potential for nuclear weapon production and highlights the importance of reactor operational strategy in maximizing the production of weapons-grade plutonium in MAGNOX reactors.

Graphite Isotope Ratio Method (GIRM) can be used to estimate plutonium production in a graphite-moderated reactor. This study presents verification results for the GIRM combined with a 3-D polynomial regression function to estimate cumulative plutonium production in a graphite-moderated reactor. Using the 3-D Monte-Carlo method, verification was done by comparing the cumulative plutonium production with the GIRM. The GIRM can estimate plutonium production for specific sampling points using a function that is based on an isotope ratio of impurity elements. In this study, the 10B/11B isotope ratio was chosen and calculated for sampling points. Then, 3-D polynomial regression was used to derive a function that represents a whole core cumulative plutonium production map. To verify the accuracy of the GIRM with polynomial regression, the reference value of plutonium production was calculated using a Monte-Carlo code, MCS, up to 4250 days of depletion. Moreover, the amount of plutonium produced in certain axial layers and fuel pins at 1250, 2250, and 3250 days of depletion was obtained and used for additional verification. As a result, the difference in the total cumulative plutonium production based on the MCS and GIRM results was found below 3.1% with regard to the root mean square (RMS) error.