稳态及瞬态中子输运模型的JFNK求解方法研究

Jacobian-Free Newton Krylov Methods for Steady and Transient Neutron Transport Models

  • 摘要: 三维全堆芯pin-by-pin中子输运模型的高效加速方法是核反应堆高精度计算的重点和难点。本文有效融合课题组开发的并行多维离散纵坐标(SN)中子输运程序comeSn和Jacobian-Free Newton Krylov(JFNK)通用求解框架comeJFNK的高效并行特性、鲁棒性和强收敛性,开发了一套三维稳态及瞬态中子输运模型的JFNK并行求解程序comeSn_JFNK。为了提高计算效率,选择中子标通量密度(而非中子角通量密度)作为JFNK全局求解变量,并利用基于空间区域并行的KBA输运扫描方法和物理预处理方法,分别构建了稳态及瞬态模型的JFNK统一残差计算模型。计算结果表明,comeSn_JFNK相比于comeSn,计算效率具有显著优势,对于三维pin-by-pin稳态KAIST-3A算例,加速比为10倍以上;对于栅元均匀化的二维七群瞬态C5G7-TD2系列基准算例,加速比约为30倍。

     

    Abstract: There is an urgent need to reduce the computational costs and improve the convergence rate for the three-dimensional (3D) high-fidelity pin-by-pin full core simulation. Therefore, the efficient and robust acceleration method of 3D large-scale pin-by-pin neutron transport models is a primary objective of high-resolution reactor calculations. In this paper, comeSn_JFNK, an efficient unified parallel solver was developed for 3D steady-state and transient pin-by-pin neutron transport models. The comeSn_JFNK solver integrated the parallel discrete ordinate (SN) neutron transport code comeSn into the parallel computational framework comeJFNK of Jacobian-Free Newton Krylov (JFNK). The comeSn code and comeJFNK framework were developed by the Virtual Reactor Coupling Analysis Team (VRCAT) at Huazhong University of Science and Technology (HUST). Therefore, comeSn_JFNK took advantages of the fast and robust convergence of the parallel JFNK framework and the high accuracy and efficiency of the SN method based on the KBA algorithm. To further improve computational efficiency of comeSn_JFNK, in the JFNK solution, neutron scalar fluxes instead of neutron angular fluxes were chosen as the global solution variables. This can reduce the number of variables and minimize the computing scale. To rapidly construct the unified JFNK residuals, the parallel KBA transport sweep methods and physics-based preconditioning techniques were utilized to improve the computational efficiency. To rapidly construct the unified JFNK residuals, the parallel KBA transport sweep methods and physics-based preconditioning techniques were utilized. This way of constructing residuals, for both steady-state and transient models, can further improve the computational efficiency. Finally, the detailed analysis of the computational accuracy and acceleration characteristics of comeSn_JFNK was presented by solving the pin-by-pin steady-state KAIST-3A and homogenized-pin transient C5G7-TD2 benchmark cases. The relative errors of radial average power density, effective multiplication factor and relative powers as a function of time were shown in this paper. There are almost the same numerical solutions between the parallel comeSn code and the parallel solver called comSn_JFNK. Numerical results also show that the parallel solver called comeSn_JFNK can achieve an acceleration of over 10 times for KAIST-3A benchmark problems and approximately 30 times for C5G7-TD2 cases compared to the original parallel comeSn code using the traditional source iteration/power iteration methods. These solutions indicate that the JFNK methods offer significant acceleration for solving both steady-state and transient SN neutron transport models compared to traditional source iteration/power iteration methods. In summary, the JFNK method provides reliable computational accuracy and high computational efficiency, which demonstrates the potential and advantages to accelerate neutron transport solutions. It also establishes a foundation for efficient simultaneous solution of more complicated pin-by-pin neutron transport and thermal-hydraulic coupling problems in the nuclear reactor cores using the unified JFNK framework in the coupling multiphysics environment (COME) developed by VRCAT at HUST.

     

/

返回文章
返回