流水线并行JFNK方法及在中子k本征值问题中的应用  

Pipelined Parallel JFNK Method and its Application in Neutron k Eigenvalue Problem

作　　者：刘礼勋 张汉 邬颖杰 郭炯[1] 李富[1] Liu Lixun;Zhang Han;Wu Yingjie;Guo Jiong;Li Fu(Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education,Institute of Nuclear and New Energy Technology,Tsinghua University,Beijing,100084,China)

机构地区：[1]清华大学核能与新能源技术研究院,先进反应堆工程与安全教育部重点实验室,北京100084

出　　处：《核动力工程》2023年第5期15-22,共8页Nuclear Power Engineering

基　　金：北京市自然科学基金资助项目(1212012)。

摘　　要：JFNK(Jacobian-free Newton-Krylov)方法是求解中子k本征值和反应堆多物理场耦合等非线性问题的高效加速方法,其中的Krylov迭代常用广义极小残差法(简称GMRES)。并行JFNK方法是实现更大规模问题求解的必要手段,其核心是解决GMRES中Gram-Schmidt (简称GS)正交化过程集合通信多、并行效率低的问题。本文以三维中子k本征值问题为研究对象,开发了基于消息传递接口并行编程模型和空间区域分解技术的并行JFNK方法。针对GS正交化过程并行可扩展性差的问题,分析讨论了流水线方法,以提高并行JFNK的并行效率,并对比了采用经典GS正交化的并行JFNK、采用修正GS正交化的并行JFNK和采用流水线方法的并行JFNK的计算时间和并行效率。最后选用三维扩散基准题IAEA-3D进行了数值测试,测试结果表明采用流水线方法的并行JFNK并行效率显著高于使用经典或修正GS正交化的并行JFNK,且收敛性未受影响。JFNK(Jacobian-free Newton-Krylov) is an efficient acceleration method for solving nonlinear problems such as neutron k eigenvalue and coupling of multiple physical fields in the reactor,and the generalized minimum residual(GMRES) algorithm is commonly used in Krylov iteration.The parallel JFNK method is a necessary means for solving large-scale problems,and its main problem lies in the low parallel efficiency of Gram-Schmidt(GS) orthogonalization procedure in GMRES,which causes massive collective communications.In this paper,the parallel JFNK method based on the parallel programming model of message passing interface and spatial domain decomposition technology is developed for the three-dimensional neutron k eigenvalue problem.Aiming at the poor parallel scalability of GS orthogonalization procedure,the pipelined method is studied to improve the parallel efficiency of parallel JFNK.Then,the computation time and parallel efficiency of parallel JFNK using classical GS orthogonalization,modified GS orthogonalization and pipelined method are compared.The IAEA-3D three-dimensional diffusion benchmark problem was used for numerical test.The results show that the parallel efficiency of pipelined parallel JFNK is significantly superior to that using classical or modified GS orthogonalization,and the convergence of pipelined parallel JFNK is not affected.

关 键 词：流水线方法 并行JFNK 消息传递接口(MPI) 集合通信 非阻塞通信 并行效率 

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