Research on fully-implicit numerical algorithm of two-fluid three-field two-phase flow model based on JFNK method

FAN Jie,WU Hexin,GOU Junli,SHAN Jianqiang
DOI: https://doi.org/10.11889/j.0253-3219.2022.hjs.45.040605
2022-01-01
Nuclear Techniques
Abstract:BackgroundMost nuclear reactor safety analysis codes using the two-fluid three-field two-phase flow model are based on semi-implicit numerical algorithm. The stability of the numerical scheme is limited to the Courant value. Few codes based on two-fluid three-field model are solved by classic Newton iteration method, which needs a Jacobian Matrix. The writing of the Jacobian matrix is difficult.PurposeIn order to improve the stability of the numerical scheme and avoid writing the Jacobian matrix, JFNK (Jacobian-free Newton-Krylov) method is adopted to fully-implicitly solve the two-fluid three-field two-phase flow model.MethodsThe conservation equations of the two-fluid three-field two-phase flow model are discretized based on the staggered grid and finite volume difference. The equation set is solved by JFNK method. When one phase or two phases of the three-field model aren't existed, the fraction of the non-existent phase is assumed as a tiny value, which can avoid the Jacobian matrix and precondition matrix being singular.ResultsBy simulating the Ransom water faucet test, the accuracy of the numerical algorithm is verified. The numerical treatments for the non-existent phase are feasible. For the simulations of Barotolimei sub-cooled boiling experiments and dryout experiments, the performances of the developed numerical algorithm for the single-phase flow, two-phase flow and the transition from the single-phase flow to the two-phase flow are validated.ConclusionsBy simulating the water faucet test and pipe heat transfer experiments, the fully-implicit numerical algorithm based on JFNK method for two-fluid three-field two-phase flow model is successfully achieved. It lays a foundation for the advanced nuclear reactor safety analysis code based on the two-fluid three-field two-phase flow model and fully-implicit numerical algorithm.
What problem does this paper attempt to address?