基于GPU并行的二维时空中子动力学MOC程序开发及验证  

Development and Verification of 2D Kinetic MOC Code Based on GPU

作　　者：邹航 梁亮 张乾 宋佩涛 赵强 ZOU Hang;LIANG Liang;ZHANG Qian;SONG Peitao;ZHAO Qiang(Fundamental Science on Nuclear Safety and Simulation Technology Laboratory,Harbin Engineering University,Harbin 150001,China;Nuclear Energy Creative Power Co.,Ltd.,Xi’an 710077,China;China Institute for Radiation Protection,Taiyuan 030006,China)

机构地区：[1]哈尔滨工程大学核安全与仿真技术国防重点学科实验室,黑龙江哈尔滨150001 [2]西安核创能源科技有限公司,陕西西安710077 [3]中国辐射防护研究院,山西太原030006

出　　处：《原子能科学技术》2022年第1期87-95,共9页Atomic Energy Science and Technology

基　　金：国家自然科学基金(12105063);核反应堆系统设计技术重点实验室运行基金(HT-KFKT-02-2019004);核数据重点实验室稳定支持基金(JCKY2021201C154);中核集团“青年英才”项目;黑龙江省自然科学基金(LH2020A001)。

摘　　要：目前特征线方法(MOC)被广泛应用于反应堆精细中子输运计算。为提高基于MOC方法的时空中子动力学输运计算效率,本文开发了ALPHA程序的动力学计算模块,实现了基于GPU并行的二维精细动力学输运计算。同时,实现了基于GPU并行的CMFD加速计算,并对TWIGL基准题和MINI-CORE基准题进行验证。数值结果显示,基于GPU并行的中子动力学计算方法能保证良好的计算精度,且具有明显的加速效果。Time-dependent,whole-core calculation with high-fidelity pin-resolved details serves an important role in multi-physics reactor applications.However,the limitation of the current CPU has become a sore point.In order to maintain the high accuracy and reduce the computational burden,it is necessary to develop the advanced technology to speed-up the transient calculation.Recently,heterogeneous computing has been increasingly widespread used with the high-performance computing(HPC)systems,which is widely equipped in HPC clusters.Compared to CPU,current GPU has much higher rate of FLOPs/s(floating-point operations per second),larger memory bandwidth,and lower energy consumption per FLOP.These features leverage the parallel compute engine in GPUs to solve many complex computational problems in a more efficient way.Method of Characteristics(MOC)has been widely used in GPU-based whole-core transport calculation for its feature of natural massive parallelization.Based on the above points,a GPU-based 2D MOC transient fixed source problem(TFSP)solver is implemented on lattice physics code ALPHA(advanced lattice physics code based on heterogeneous architecture).In this paper,the full implicit method(FIM)were adopted to solve the TFSP.The Jacobi transport sweep algorithm was introduced in MOC TFSP solver.The transport sweep over energy groups and polar angles were performed for each ray segment during the ray tracing.To accelerate the convergence of the TFSP,the pin-based coarse mesh finite difference(CMFD)was implemented in the TFSP solver.In order to parallelize the CMFD TFSP solver on GPU,the red-black ordering algorithm was introduced in the CMFD TFSP solver.The results of two problem are presented including the TWIGL benchmark problem and the 2D MINI-CORE benchmark problem to verify the capability of the GPU-based TFSP solver.Numerical results demonstrate that the TFSP solver built in ALPHA has the desired accuracy.The percent difference of core power history between ALPHA and the reference are less than 0.5%.Compared wi

关 键 词：MOC 时空中子动力学 CMFD 输运计算 GPU加速 

分 类 号：TL329[核科学技术—核技术及应用]