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 d...

Full description

Saved in:
Bibliographic Details
Main Authors: Peijun Li, Chen Hao, Ning Xu, Yanling Zhu, Yizhen Wang, Zhigang Zhang
Format: Article
Language:English
Published: Elsevier 2025-02-01
Series:Nuclear Engineering and Technology
Subjects:
Online Access:http://www.sciencedirect.com/science/article/pii/S1738573324004583
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:1738-5733