An innovative and efficient implementation of matrix-free Newton krylov method for neutronics/thermal-hydraulics coupling simulation

The core physical behavior of reactors is essentially the result of multi-physical fields coupling feedback. High-fidelity neutronics/thermal-hydraulics (N/TH) analysis can simulate and predict nuclear reactor core phenomena realistically, providing advanced and reliable technical means during the design and safety analysis of nuclear reactor. In this work, an efficient and robustness coupling method using power density as the coupling parameter, Matrix-Free Newton Krylov (MFNK) method, is successfully developed and innovatively implemented in HNET for high-fidelity N/TH coupling simulation. To enhance the efficiency and stability, the multi-level generalized equivalence theory-based CMFD (ML-gCMFD) iterative acceleration method and ML-gCMFD coupling acceleration method are proposed. In addition, the nonlinear preconditioning and hybrid perturbation size formula are implemented to further improve the convergence. Finally, to evaluate the numerical accuracy, convergence, efficiency and stability of MFNK method, a series of representative problems, including a three-dimensional (3D) single fuel pin problem, VERA Benchmark Problem 6, and VERA Benchmark Problem 7, are analyzed by comparing with the current N/TH coupling methods. Numerical results indicate that MFNK method can obtain strong stability, high convergence performance, and relatively high computational efficiency while ensuring high accuracy. It demonstrates that MFNK method has significant performance advantages and potential for high-fidelity N/TH coupling simulation.