GPU-accelerated 3D neutron diffusion code based on finite difference method
- Dept. of Engineering Physics, Tsinghua Univ. (China)
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
- Research Organization:
- American Nuclear Society, Inc., 555 N. Kensington Avenue, La Grange Park, Illinois 60526 (United States)
- OSTI ID:
- 22105600
- Resource Relation:
- Conference: PHYSOR 2012: Conference on Advances in Reactor Physics - Linking Research, Industry, and Education, Knoxville, TN (United States), 15-20 Apr 2012; Other Information: Country of input: France; 7 refs.
- Country of Publication:
- United States
- Language:
- English
Similar Records
SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation
Quantum Monte Carlo Endstation for Petascale Computing
Related Subjects
22 GENERAL STUDIES OF NUCLEAR REACTORS
ACCELERATION
ACCURACY
COMPUTER CODES
COMPUTER GRAPHICS
DIFFUSION
EXTRAPOLATION
FINITE DIFFERENCE METHOD
GHZ RANGE 01-100
IMPLEMENTATION
NEUTRON DIFFUSION EQUATION
NEUTRONS
POLYNOMIALS
PWR TYPE REACTORS
THREE-DIMENSIONAL CALCULATIONS