DOI QR코드

DOI QR Code

Practical methods for GPU-based whole-core Monte Carlo depletion calculation

  • Received : 2023.01.26
  • Accepted : 2023.04.15
  • Published : 2023.07.25

Abstract

Several practical methods for accelerating the depletion calculation in a GPU-based Monte Carlo (MC) code PRAGMA are presented including the multilevel spectral collapse method and the vectorized Chebyshev rational approximation method (CRAM). Since the generation of microscopic reaction rates for each nuclide needed for the construction of the depletion matrix of the Bateman equation requires either enormous memory access or tremendous physical memory, both of which are quite burdensome on GPUs, a new method called multilevel spectral collapse is proposed which combines two types of spectra to generate microscopic reaction rates: an ultrafine spectrum for an entire fuel pin and coarser spectra for each depletion region. Errors in reaction rates introduced by this method are mitigated by a hybrid usage of direct online reaction rate tallies for several important fissile nuclides. The linear system to appear in the solution process adopting the CRAM is solved by the Gauss-Seidel method which can be easily vectorized on GPUs. With the accelerated depletion methods, only about 10% of MC calculation time is consumed for depletion, so an accurate full core cycle depletion calculation for a commercial power reactor (BEAVRS) can be done in 16 h with 24 consumer-grade GPUs.

Keywords

Acknowledgement

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021M2D6A1048220).

References

  1. T.M. Pandya, S.R. Johnson, T.M. Evans, G.G. Davidson, S.P. Hamilton, A.T. Godfrey, Implementation, capabilities, and benchmarking of Shift, a massively parallel Monte Carlo radiation transport code, J. Comput. Phys. 308 (Mar. 2016) 239-272, https://doi.org/10.1016/j.jcp.2015.12.037.
  2. P.K. Romano, B. Forget, The OpenMC Monte Carlo particle transport code, Ann. Nucl. Energy 51 (Jan. 2013) 274-281, https://doi.org/10.1016/j.anucene.2012.06.040.
  3. S.P. Hamilton, S.R. Slattery, T.M. Evans, Multigroup Monte Carlo on GPUs: comparison of history- and event-based algorithms, Ann. Nucl. Energy 113 (Mar. 2018) 506-518, https://doi.org/10.1016/j.anucene.2017.11.032.
  4. S.P. Hamilton, T.M. Evans, Continuous-energy Monte Carlo neutron transport on GPUs in the Shift code, Ann. Nucl. Energy 128 (Jun. 2019) 236-247, https://doi.org/10.1016/j.anucene.2019.01.012.
  5. J.R. Tramm, et al., Toward portable GPU acceleration of the OpenMC Monte Carlo particle transport code, in: International Conference on Physics of Reactors, PHYSOR 2022, 2022, pp. 2734-2743, https://doi.org/10.13182/PHYSOR22-37847.
  6. M.S. Mckinley, et al., Status of LLNL Monte Carlo transport codes on sierra GPUs, in: International Conference on Mathematics Computational Methods Applied to Nuclear Science and Engineering, M&C, 2019, pp. 2160-2165, 2019.
  7. N. Choi, K.M. Kim, H.G. Joo, Initial Development of PRAGMA-A GPU-Based Continuous Energy Monte Carlo Code for Practical Applications, Transactions of the Korean Nuclear Society Autumn Meeting, Oct. 2019.
  8. N. Choi, H.G. Joo, Analytic treatment of intra-fuel-rod temperature distributions in the GPU-based continuous energy Monte Carlo code PRAGMA, Trans. Am. Nucl. Soc. 122 (2020) 721-724, https://doi.org/10.13182/T122-31902.
  9. N. Choi, H.G. Joo, Domain decomposition for GPU-Based continuous energy Monte Carlo power reactor calculation, Nucl. Eng. Technol. 52 (11) (2020) 2667-2677, Nov, https://doi.org/10.1016/j.net.2020.04.024.
  10. N. Choi, K.M. Kim, J. Im, H.G. Joo, Algorithms and performance of GPU-based continuous energy Monte Carlo code PRAGMA, in: Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo, SNA+MC 2020), 2020, pp. 373-380.
  11. G.G. Davidson, et al., Nuclide depletion capabilities in the Shift Monte Carlo code, Ann. Nucl. Energy 114 (Apr. 2018) 259-276, https://doi.org/10.1016/j.anucene.2017.11.042.
  12. W. Haeck, B. Verboomen, An optimum approach to Monte Carlo burnup, Nucl. Sci. Eng. 156 (2) (2007) 180-196, https://doi.org/10.13182/NSE07-A2695.
  13. E. Brun, E. Dumonteil, F. Malvagi, Systematic uncertainty due to statistics in Monte Carlo burnup codes: application to a simple benchmark with TRIPOLI4-D, Progress in NUCLEAR SCIENCE and TECHNOLOGY 2 (2011) 879-885. https://doi.org/10.15669/pnst.2.879
  14. J.L. Salcedo-Perez, B. Forget, K. Smith, P.K. Romano, Hybrid tallies to improve performance in depletion Monte Carlo simulations, in: International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M&C, 2019, pp. 927-936, 2019.
  15. K.M. Kim, N. Choi, H.G. Lee, H.G. Joo, Initial development of depletion capability in the GPU-based Monte Carlo code PRAGMA, in: Transactions of the Korean Nuclear Society Autumn Meeting, Dec. 2020.
  16. M. Pusa, Higher-order Chebyshev rational approximation method and application to burnup equations, Nucl. Sci. Eng. 182 (3) (2016) 297-318, https://doi.org/10.13182/NSE15-26.
  17. H.J. Shim, B.S. Han, S.J. Jong, H.J. Park, C.H. Kim, McCard: Monte Carlo code for advanced reactor design and analysis, Nucl. Eng. Technol. 44 (2) (2012) 161-176, https://doi.org/10.5516/NET.01.2012.503.
  18. H.J. Shim, C.H. Kim, W.S. Park, H.K. Joo, Monte Carlo depletion analysis of a PWR with the MCNAP, in: Mathematics and Computation, Reactor Physics and Environmental Analysis in Nuclear Applications (M&C99), 1999, pp. 1026-1035.
  19. A.G. Croff, User's Manual for the ORIGEN2 Computer Code, 1980, https://doi.org/10.2172/5285077. United States.
  20. R.J. Stamm'ler, HELIOS Methods, Studsvik Scandpower, 1998.
  21. J.R. Askew, F.J. Fayers, P.B. Kemshell, General description of the lattice code WIMS, J. Br. Nucl. Energy Soc. 5 (1966) 564-585. https://www.osti.gov/biblio/4462252.
  22. B.T. Rearden, M.A. Jessee, SCALE Code System 6 1 (2) (2016), https://doi.org/10.2172/1326509. United States.
  23. Y. Arafa, A.-H. Badawy, G. Chennupati, N. Santhi, S. Eidenbenz, Low overhead instruction latency characterization for NVIDIA GPGPUs, IEEE High Performance Extreme Computing Conference (HPEC) (2019) 1-8. May 2019.
  24. H.G. Lee, H.G. Joo, Optimization of the GPU-based depletion solver in nTRACER, in: Transactions of the Korean Nuclear Society Autumn Meeting, Dec. 2020.
  25. N. Horelik, B. Herman, B. Forget, K. Smith, Benchmark for evaluation and validation of reactor simulations (BEAVRS), in: International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, 2013, pp. 2986-2999 [Online]. Available: https://www.osti.gov/biblio/22212933.
  26. K.M. Kim, J. Im, N. Choi, H.G. Lee, H.G. Joo, High-performance and high-fidelity GPU based Monte Carlo solutions to the BEAVRS benchmark, Nucl. Sci. Eng. (Jan. 2023) 1-22, https://doi.org/10.1080/00295639.2022.2148812.
  27. J. Im, N. Choi, H.G. Joo, A study on the optimal use of ramp-up technique under massive particle condition with CMFD acceleration, in: Transactions of the Korean Nuclear Society Spring Meeting, May 2020.