Abstract:The traditional fixed-point iteration method is typically used for neutronics/thermal-hydralics coupling problems in most nuclear safety analysis codes. But the fixed-point iteration method has a tendency to fail to be used in computing the coupling problems due to slow convergence rate in some cases and even no convergence, and thus resulting in a limited efficiency, especially for the tight-coupling and fast-transient problems. In addition, for the reactor thermal-hydraulic calculation, the traditional finite difference or volume method (FDM or FVM) is used. However, both FDM and FVM require fine mesh size to achieve the desired precision and thus also result in a limited efficiency for the large scale problems. In this paper, to ensure the accuracy, efficiency and convergence for large-scale and complicated coupling problems, the new methods-NEM_JFNK are successfully developed to simultaneously solve the neutronics/thermal-hydralics coupling problems by combining the advantage of the efficient coarse nodal expansion method (NEM) and Jacobian-Free Newton-Krylov method (JFNK). The NEM has been widely used in the reactor physics analysis due to its high efficiency and accuracy in the reactor physics analysis, and it has proved to be superior to FDM and FVM. To improve the efficiency and accuracy for the large scale problems, the NEM is first extended to thermal hydraulic problems from the reactor physics calculation. Then all the governing equations of the neutronics/thermal-hydralics coupling problems can be discretized by the NEM and all the variables can be solved on the coarse meshes so that the size of coupling problems is greatly reduced. To ensure the high accuracy for the coupling problems on the coarse meshes, the high-order coefficients in NEM are successfully transferred between the coupling terms by our research. After that, to ensure the convergence of complicated coupling problems, JFNK based on the NEMs needs to be developed. However, the researches of JFNK based on the NEM in reactor analysis are less and the existing JFNK methods are mostly based on FVM or FDM or the finite element method. In this paper, the NEM discretization equations are successfully integrated into the framework of JFNK through the special treatment and the NEM_JFNK with linear-based preconditioner named LP_NEM_JFNK is also successfully developed. In addition, to take advantage of the existing code and avoid the construction of residual formulations, the non-residual construction NEM_JFNK named NRC_NEM_JFNK is presented and the black-box coupling method is achieved by NRC_NEM_JFNK so that the existing codes only need the simple modification to achieve the combination of the NEM and JFNK. Numerical results show LP_NEM_JFNK and NRC_NEM_JFNK outperform traditional fixed-point iteration method in the sense of convergence rate and efficiency. Further studies are needed to extend the NEM_JFNK method to the multi-dimensional neutronic/thermal hydraulic coupling problems in the high temperature gas-cooled reactor.

A parallel in energy Krylov solution algorithm for the variational nodal method

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

Eigenvalue Solvers for Modeling Nuclear Reactors on Leadership Class Machines

An Efficient Parallel Krylov-Schur Method for Eigen-Analysis of Large-Scale Power Systems

Application of Finite Difference Jacobian Based Newton-Krylov Method for Coupled Neutronics/Thermal Conduction Problems of Nuclear Reactor

Simultaneous Approach for Simulation of A High-Temperature Gas-Cooled Reactor

A Parallel Function Evaluation Approach for Solution to Large-Scale Equation-Oriented Models.

Application of Nonlinear Krylov Solvers for Conjugate Heat Transfer Simulations of Electrical Battery Packs

Parallel Energy-Stable Solver for a Coupled Allen--Cahn and Cahn--Hilliard System

On the Use of the Jacobian-Free Newton Krylov Method to Generate One-Group Discontinuity and Super Homogenization Factors for Full-Core Neutron Diffusion Simulations

Research On Nodal Expansion Method For Transient Convection Diffusion Equation

Jacobian-Free Newton-Krylov Nodal Expansion Methods with Physics-Based Preconditioner and Local Elimination for Three-Dimensional and Multigroup K-Eigenvalue Problems

A Fully Implicit Parallel Solver for MHD Instabilities in a Tokamak

Towards Robust Solvers for Nuclear Fusion Simulations Using JOREK: A Numerical Analysis Perspective

A Parallel Coloring Newton-Krylov Method for Multiphysics Coupling System in Nuclear Reactors

Efficient Simultaneous Solution of Multi-Physics Multi-Scale Nonlinear Coupled System in HTR Reactor Based on Nonlinear Elimination Method

Jacobian-Free Newton-Krylov Based on Nodal Expansion Method for Neutronic-Thermal Hydraulic Coupling Problem

Jacobian-free Newton Krylov Nodal Expansion Method in Three-Dimensional Cylindrical Coordinates

A Second Moment Method for k-Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

Computing Eigenvalue Sensitivity Coefficients to Nuclear Data by Adjoint Superhistory Method and Adjoint Wielandt Method Implemented in RMC Code

A Modified JFNK with Line Search Method for Solving K-Eigenvalue Neutronics Problems with Thermal-Hydraulics Feedback