Practical methods for GPU-based whole-core Monte Carlo depletion calculation
0
조회
0
다운로드
0
인용
0
참고문헌
키워드
초록
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
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
참고문헌 (27)
- Tara M. Pandya, "Implementation, capabilities, and benchmarking of Shift, a massively parallel Monte Carlo radiation transport code", Journal of Computational Physics, 308, pp.239~272 (2016), Elsevier BV[DOI]
- Paul K. Romano, "The OpenMC Monte Carlo particle transport code", Annals of Nuclear Energy, 51, pp.274~281 (2013), Elsevier BV[DOI]
- Steven P. Hamilton, "Multigroup Monte Carlo on GPUs: Comparison of history- and event-based algorithms", Annals of Nuclear Energy, 113, pp.506~518 (2018), Elsevier BV[DOI]
- Steven P. Hamilton, "Continuous-energy Monte Carlo neutron transport on GPUs in the Shift code", Annals of Nuclear Energy, 128, pp.236~247 (2019), Elsevier BV[DOI]
- J.R. Tramm, "Toward portable GPU acceleration of the OpenMC Monte Carlo particle transport code", International Conference on Physics of Reactors, PHYSOR 2022, pp.2734~2743 (2022)
- M.S. Mckinley, "Status of LLNL Monte Carlo transport codes on sierra GPUs", International Conference on Mathematics Computational Methods Applied to Nuclear Science and Engineering, M&C, 2019, pp.2160~2165 (2019)
- N. Choi, "Initial Development of PRAGMA-A GPU-Based Continuous Energy Monte Carlo Code for Practical Applications", Transactions of the Korean Nuclear Society Autumn Meeting (2019)
- N. Choi, "Analytic Treatment of Intra-Fuel-Rod Temperature Distributions in the GPU-Based Continuous Energy Monte Carlo Code PRAGMA", Transactions of the American Nuclear Society, 122, pp.721~724 (2020), AMNS[DOI]
- 최남재, "Domain decomposition for GPU-Based continuous energy Monte Carlo power reactor calculation", Nuclear Engineering and Technology, 52(11), pp.2667~2677 (2020), 한국원자력학회[DOI]
- N. Choi, "Algorithms and performance of GPU-based continuous energy Monte Carlo code PRAGMA", Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo, SNA+MC 2020, pp.373~380 (2020)
- Gregory G. Davidson, "Nuclide depletion capabilities in the Shift Monte Carlo code", Annals of Nuclear Energy, 114, pp.259~276 (2018), Elsevier BV[DOI]
- W. Haeck, "An Optimum Approach to Monte Carlo Burnup", Nuclear Science and Engineering, 156(2), pp.180~196 (2007), Informa UK Limited[DOI]
- E. Brun, "Systematic uncertainty due to statistics in Monte Carlo burnup codes : application to a simple benchmark with TRIPOLI-4-D", Progress in NUCLEAR SCIENCE and TECHNOLOGY, 2, pp.879~885 (2011)
- J.L. Salcedo-Perez, "Hybrid tallies to improve performance in depletion Monte Carlo simulations", International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M&C, 2019, pp.927~936 (2019)
- K. M. Kim, "Initial development of depletion capability in the GPU-based Monte Carlo code PRAGMA", Transactions of the Korean Nuclear Society Autumn Meeting (2020)
- Pusa Maria, "Higher-Order Chebyshev Rational Approximation Method and Application to Burnup Equations", Nuclear Science and Engineering, 182(3), pp.297~318 (2016), Informa UK Limited[DOI]
- 심형진, "MCCARD: MONTE CARLO CODE FOR ADVANCED REACTOR DESIGN AND ANALYSIS", Nuclear Engineering and Technology, 44(2), pp.161~176 (2012), 한국원자력학회
- H.J. Shim, "Monte Carlo depletion analysis of a PWR with the MCNAP", Mathematics and Computation, Reactor Physics and Environmental Analysis in Nuclear Applications (M&C99), pp.1026~1035 (1999)
- User's Manual for the ORIGEN2 Computer Code, "A.G. Croff", ~ (1980)
- HELIOS Methods, "R. J. Stamm’ler", Studsvik Scandpower, ~ (1998)
- J.R. Askew, "General description of the lattice code WIMS", J. Br. Nucl. Energy Soc, 5, pp.564~585 (1966)[DOI]
- SCALE Code System 6 1 (2), "B. T. Rearden", ~ (2016)
- Y. Arafa, "Low overhead instruction latency characterization for NVIDIA GPGPUs", IEEE High Performance Extreme Computing Conference (HPEC) (2019)
- H.G. Lee, "Optimization of the GPU-based depletion solver in nTRACER", Transactions of the Korean Nuclear Society Autumn Meeting (2020)
- N. Horelik, "Benchmark for evaluation and validation of reactor simulations (BEAVRS)", International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, pp.2986~2999 (2013)
- Kyung Min Kim, "High-Performance and High-Fidelity GPU-Based Monte Carlo Solutions to the BEAVRS Benchmark", Nuclear Science and Engineering, 197(8), pp.1823~1844 (2023), Informa UK Limited[DOI]
- J. Im, "A study on the optimal use of ramp-up technique under massive particle condition with CMFD acceleration", Transactions of the Korean Nuclear Society Spring Meeting (2020)