Calculation of two-fluid subchannels model of pressurized water reactor: Picard Krylov method

Yuhang Zhang,Zhaofei Tian,Lei Li,Guangliang Chen,Hao Qian,Lixuan Zhang,Yuguan Jin

DOI: https://doi.org/10.1016/j.icheatmasstransfer.2024.107356

IF: 6.782

2024-03-09

International Communications in Heat and Mass Transfer

Abstract:A Picard Krylov method for solving the two-fluid subchannel analysis model has been proposed to improve the accuracy and efficiency of thermal conduction and convection coupling calculations in the pressurized water reactor core. Generally speaking, the Picard algorithm has high coupling accuracy and good robustness, but has low computational efficiency, while the Newton Krylov method has higher coupling accuracy and efficiency than the Picard iteration method, but has low ability to solve ill-conditioned matrices. In order to obtain higher convergence accuracy, robustness and calculation efficiency, the Picard Krylov method has been proposed, in which the Krylov subspace method was used to improve the local Picard iteration calculation efficiency. Considering that the pressure distribution of the pressurized water reactor is affected by both flow and heating, pressure update equations were generated from mass-energy equations processed by the Gauss elimination, and the Krylov subspace method was used to improve the calculation efficiency of pressure update equations. Finally, in order to further improve the computational efficiency of the Krylov subspace method in solving pressure update equations, a physics-based precondition matrix was proposed. Numerical experiments showed that the Picard Krylov method had good calculation accuracy and efficiency in the OECD/NRC Benchmark based on NUPEC PWR Subchannel and Bundle Test calculation.

thermodynamics,mechanics