DOI QR코드

DOI QR Code

Fission source convergence diagnosis in Monte Carlo eigenvalue calculations by skewness and kurtosis estimation methods

  • Received : 2024.01.27
  • Accepted : 2024.04.18
  • Published : 2024.09.25

Abstract

In this study, a skewness estimation method (SEM) and kurtosis estimation method (KEM) are introduced to determine the number of inactive cycles in Monte Carlo eigenvalue calculations. The SEM and KEM can determine the number of inactive cycles on the basis that fully converged fission source distributions may follow normal distributions without asymmetry or outliers. Two convergence criteria values and a minimum cycle length for the SEM and KEM were determined from skewness and kurtosis analyses of the AGN-201K benchmark and 1D slab problems. The SEM and KEM were then applied to two OECD/NEA slow convergence benchmark problems to evaluate the performance and reliability of the developed methods. Results confirmed that the SEM and KEM provide appropriate and effective convergence cycles when compared to other methods and fission source density fraction trends. Also, the determined criterion value of 0.5 for both ε1 and ε2 was concluded to be reasonable. The SEM and KEM can be utilized as a new approach for determining the number of inactive cycles and judging whether Monte Carlo tally values are fully converged. In the near future, the methods will be applied to various practical problems to further examine their performance and reliability, and optimization will be performed for the convergence criteria and other parameters as well as for improvement of the methodology for practical usage.

Keywords

Acknowledgement

This work was supported by the National Research Foundation of Korea (NRF) grant funded by theKorean government (MSIT) (NRF2022M2C7A1A01053531).

References

  1. L.L. Carter, N.J. McCormick, Source convergence in Monte Carlo calculations, Nucl. Sci. Eng. 36 (1969) 438.
  2. T. Kitada, T. Takeda, Effective convergence of fission source distribution in Monte Carlo simulation, J. Nucl. Sci. Technol. 38 (2001) 324-329.
  3. T. Ueki, F.B. Brown, Autocorrelation and dominance ratio in Monte Carlo criticality calculations, Nucl. Sci. Eng. 145 (2003) 279-290.
  4. Y. Naito, J. Yang, The sandwich method for determining source convergence in Monte Carlo calculation, J. Nucl. Sci. Technol. (Tokyo, Jpn.) 41 (5) (2004) 559-568.
  5. T. Ueki, F.B. Brown, Stationarity modelling and informatics-based diagnostics in Monte Carlo criticality calculation, Nucl. Sci. Eng. 149 (2005) 38-50.
  6. H.J. Shim, C.H. Kim, Stopping criteria of inactive cycle Monte Carlo calculations, Nucl. Sci. Eng. 157 (2007) 132-141.
  7. S.H. Kim, et al., A new approach to check and diagnose the fission source convergence in Monte Carlo criticality calculations, Nucl. Sci. Eng. 178 (2014) 29-41.
  8. T.A. Jones, Skewness and Kurtosis as criteria of normality in observed frequency distributions, J. Sediment. Res. 39 (1969) 1622-1623.
  9. D.N. Joanes, C.A. Gill, Comparing measures of sample skewness and kurtosis, The Statistician 47 (1) (1998) 183-189.
  10. H.J. Park, et al., Preliminary benchmark evaluation of AGN-201K educational and research reactor, in: Proceedings of the Reactor Physics Asia 2023 (RPHA2023) Conference, 2023. Gyeongju, Korea, October 24-26.
  11. T. Endo, A. Yamamoto, K. Sakata, Theoretical prediction on underestimation of statistical uncertainty for fission rate tally in Monte Carlo calculation, in: PHYSOR2014 - the Role of Reactor Physics toward a Sustainable Future, the Westin Miyako, 2014. Kyoyo, Japan, September 28 - October 3.
  12. R.N. Blomquist, M. Armishaw, et al., Source convergence in criticality safety analyses phase I: results for four test problems, Nuclear Science ISBN 92-64-02304-6, OECD 2006 NEA No. (2006) 5431.
  13. H.J. Shim, et al., McCARD: montel Carlo code for advanced reactor design and analysis, Nucl. Eng. Technol. 44 (2012) 161-176.
  14. R.N. Blomquist, et al., OECD/NEA source convergence benchmark program: overview and summary of results, in: ICNC 2003 - 7th International Conference on Nuclear Criticality Safety, 2003. Tokai-mura, Japan, October 20-24.