Jacobian-free Newton Krylov coarse mesh finite difference algorithm based on high-order nodal expansion method for three-dimensional nuclear reactor pin-by-pin multiphysics coupled models

Abstract

The nuclear reactor core is a large-scale, complicated and tight-coupled nonlinear multiphysics coupled system including neutron transport, heat conduction of nuclear fuel pins, coolant flow, heat transfer and so on. To meet the demands of the advanced nuclear reactor analysis and design, there is an urgent need to reduce the computational costs and improve the convergence rate for the three-dimensional (3D) nuclear reactor core multiphysics coupled models, especially for a high-fidelity pin-by-pin (at the level of fuel pin by fuel pin) full core simulation. In this work, we develop CNJ, an efficient Jacobian free Newton Krylov (JFNK) coarse mesh finite difference (CMFD) algorithm based on the nodal expansion method (NEM) to simultaneously solve the complicated pin-by-pin neutronics-thermal hydraulic (N-TH) coupled models. The high-order NEM can improve the numerical accuracy of the CMFD method on a coarse mesh and the JFNK method can ensure the rapid convergence and high efficiency on large-scale and nonlinear coupled discrete systems. By combining the CMFD, NEM and JFNK methods and making full use of their respective advantages, the CNJ algorithm can obtain high accuracy and efficiency even on a coarse mesh to greatly reduce the number of solution variables and the computational cost of the 3D coupled models. Finally, two representative and complicated 3D pin-by-pin N-TH coupled models are analyzed to evaluate the numerical accuracy and efficiency advantage of CNJ in detailed.