有限元SN中子输运模拟的区域分解并行  

Parallelization of Finite Element SN Neutron Transport Simulation Based on Domain Decomposition

作　　者：郭海兵[1] 黄洪文[1] 马纪敏[1] 丁文杰[1] GUO Haibing;HUANG Hongwen;MA Jimin;DING Wenjie(Institute of Nuclear Physics and Chemistry,China Academy of Engineering Ph ysics,Mianyang 621900,China)

机构地区：[1]中国工程物理研究院核物理与化学研究所,四川绵阳621900

出　　处：《原子能科学技术》2020年第6期1074-1084,共11页Atomic Energy Science and Technology

基　　金：中国工程物理研究院战略科技指南项目资助(TCGH0213);国家自然科学基金青年科学基金项目资助(11405155)。

摘　　要：确定论中子输运方法具有计算速度快、可获取物理量的精细场分布、可高效多物理耦合等优点,随着有限元方法在中子输运模拟中的应用,复杂几何结构、大尺度下的屏蔽问题和临界问题都能得到高保真建模和分析.离散纵标(SN)法是求解中子输运方程的有效数值方法,基于OpenMP并行机制对各独立离散方向进行并行求解,可提高SN输运模拟的计算速度,但并行规模较有限.对几何空间进行区域分解并采用MPI并行机制,可实现大规模并行扩展,进而实现对大型问题的高精度快速求解.本文采用并行自适应非结构网格应用框架JAUMIN进行区域分解和进程间通信,通过并行SN扫描实现了自主有限元输运程序ENTER的高效并行,完成正确性检验后在天河Ⅱ号超级计算机上使用1440个CPU核完成了1.43×10^7网格单元、2.81×10^9自由度规模问题的测试,计算时间约7.4h.表明该程序具备了有效模拟大型复杂结构中子输运问题的能力,具有一定工程应用价值.The deterministic neutron transport method has the advantage of fast calcula-tion speed,and is capable of giving fine field distribution of physical parameters and simulating multi-physics coupling efficiently.With the application of finite element method in neutron transport calculation,high fidelity modeling and simulation can be performed for complex geometric structures and large scale radiation shielding and criti-cality calculation.The discrete ordinate(SN)algorithm is an efficient method in solving.the Boltzmann neutron transport equation,and the parallelization of independent ordi-nate directions can be easily realized by the OpenMP scheme,but the speedup effect is limited by the number of discrete ordinate directions and the computation resources on single computer node or workstation.In order to improve the calculation speed further,the geometry domain decomposition method and the MPI parallel computing scheme can be applied to extend the parallel computing capability.The domain decomposition and communication between processes were accomplished by utilizing the parallel adaptive unstructured meshes applications infrastructure JAUMIN.And the parallelization of the formerly developed finite element SN neutron transport code ENTER was realized by implementing the parallel SN sweeping algorithm.After verifying this parallel code on a workstation,several large scale problems on Tianhe-Ⅱ supercomputer with 1440 CPUs were tested at most.The largest number of elements and degree of freedom were approximately 1.43× 10^7 and 2.81× 10^9,respectively,which cost simulation time of 7.4h.The newly developed parallel neutron transport code ENTER has the capability of efficiently simulating large scale and complex geometry problems.

关 键 词：中子输运 有限元方法 区域分解 MPI并行计算 JAUMIN 

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