基于GPU加速的三维堆芯物理程序STORK的开发与验证  

Development and Validation of 3D Core Physics Code STORK Based on GPU Acceleration

作　　者：俞陆林 杨高升 陈国华 卑华 蒋校丰 高明敏 王涛 YU Lulin;YANG Gaosheng;CHEN Guohua;BEI Hua;JIANG Xiaofeng;GAO Mingmin;WANG Tao(Shanghai NuStar Nuclear Power Technology Co.,Ltd.,Shanghai 201103,China;Nuclear Power Operation Research Institute,Shanghai 200126,China)

机构地区：[1]上海核星核电科技有限公司,上海201103 [2]核电运行研究(上海)有限公司,上海200126

出　　处：《原子能科学技术》2024年第3期662-671,共10页Atomic Energy Science and Technology

摘　　要：基于小型多GPU计算平台,采用二维全堆逐层特征线方法(MOC)和三维逐棒(pin-by-pin)三阶简化球谐函数方法(SP3方法)相耦合的方式开发了堆芯三维输运中子学计算程序STORK。在方法论方面,首先通过对堆芯各轴向层的二维MOC输运计算在线产生栅元均匀化截面以及超级均匀化修正因子(SPH因子),然后采用SP3方法进行pin-by-pin三维堆芯计算。在程序开发方面,采用了CUDA、C++和Python的混合编程,且所有计算模块都基于CUDA/C++开发,并进行了大量的性能优化。通过对C5G7三维插棒基准题和VERA基准题的验证表明,与国际上同类中子学计算软件相比,基于CPU/GPU异构系统开发的STORK程序在计算效率和计算成本方面都具有明显优势。A 3D neutron transport computational code,STORK,has been developed based on a small-scale multi-GPU computing platform,utilizing the coupled approach of the two-dimensional full-core layer-by-layer transport calculation by the method of characteristics(MOC)and the 3D pin-by-pin simplified P3(SP3)calculation.In this code,firstly,the core was layered according to the axial characteristics and the two-dimensional multi-group(69-group)transport equation was solved by MOC method(with fully reflective boundary conditions in the axial direction)for each axial layer.Secondly,utilizing the results from 2D MOC calculations,based on the equivalent homogenization theory and the super-homogenization(SPH)technology,the heterogenous cells were homogenized,which produced the few-group homogenous cross sections as well as SPH factors.Finally,the 3D whole-core pin-by-pin SP3 calculation was carried out to obtain cell flux and power distribution.Moreover,the constructive solid geometry(CSG)was applied to enhance the complex geometric modeling capability in STORK.A combination of the enhanced neutron flow method and the equivalence theory was used to perform resonance calculations and a pre-produced table of resonance interference factors was adopted to handle the resonance interference effects.During 2D transport calculation,a two-level unstructured coarse mesh finite difference method was applied to accelerate the convergence of the MOC calculation.In the 3D pin-by-pin calculation,the 3D SP3 equations were solved by the transverse integration technique and the nodal expansion method with group transverse-integrated neutron fluxes approximated by the parabola expansion in the radial direction and by semi-analytical expansion in the axial direction.In terms of code development,a hybrid programming of CUDA,C++and Python was adopted,and all the computational modules were developed based on CUDA/C++with a large number of performance optimizations,so that 2D MOC calculations at each layer of the core could be carried out on multiple GPU

关 键 词：中子输运 图形处理器 特征线方法 在线均匀化 逐棒 SP3 超级均匀化方法 

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