Abstract:The Jacobian-free Newton-Krylov (JFNK) method is a widely used and flexible numerical method for solving the neutronic/thermal-hydraulic coupling system. The main property of JFNK is that the Jacobian-vector product is evaluated approximately by finite difference, avoiding the forming and storage of Jacobian explicitly. However, the lack of an efficient preconditioner is a major bottleneck for the JFNK method, leading to poor convergence. The finite difference Jacobian-based Newton-Krylov (DJNK) method is another advanced numerical method, in which the Jacobian matrix is formed and stored explicitly. The DJNK method can provide a better preconditioner for Krylov iteration than JFNK. However, how to compute the Jacobian matrix efficiently and automatically is a key issue for the DJNK method. By fully utilizing the sparsity of the Jacobian matrix and graph coloring algorithm, the Jacobian can be computed efficiently. Unfortunately, when there are dense rows/blocks, a huge computational burden will emerge due to the lack of sparsity, resulting in the extremely poor efficiency of Jacobian computation. In this work, a Jacobian-split Newton-Krylov (JSNK) method is proposed to resolve the dense row/block problem by combining the advantages of JFNK and DJNK. The main feature of the JSNK method is to split the Jacobian matrix into sparse and dense parts. The sparse part of the Jacobian matrix is explicitly constructed using the graph coloring algorithm while for the dense part, the Jacobian-vector product is approximated by finite difference. The computational complexity of the JSNK method is analyzed and compared to the JFNK method and the DJNK method from theoretical and experimental aspects and under different meshes. A simplified two-dimensional (2-D) high-temperature gas-cooled reactor (HTR) model and a simplified 2-D pressurized water reactor model are utilized to demonstrate the superiority of the JSNK method. The numerical results show that the JSNK method successfully resolved the dense rows/blocks. More importantly, its efficiency significantly outperforms the JFNK method and the DJNK method.

Calculation of Extraction Equilibrium in Multi-Component System by Newton-Raphson Algorithm

An Undecomposed Hybrid Algorithm for Nonlinear Coupled Constitutive Relations of Rarefied Gas Dynamics

Extraction Equilibrium Algorithm for Multi-Component System with Different-Valance Ions

Unified approach for distillation column simulation and optimization based on equation-oriented algorithms

Dynamic Optimization of Multivariable Endothermic Reaction in Cascade CSTR

A Newton-Krylov Method with Dense Row Decomposition for Steady-State Neutronics/thermal-Hydraulics Problems

Models for predicting multi-component phase equilibrium ratio based on artificial neural networks

A Robust and Efficient Calculation Procedure for Determining Concentration Distribution of Multicomponent Mixtures

Simultaneous Determination of Two-Component Isotherm Parameters and Lumped Mass Transfer Coefficients in RPLC with the 0-1 Model-Inverse Method.

Accelerating Multicomponent Phase-Coexistence Calculations with Physics-informed Neural Networks

Finite Difference Jacobian Based Newton-Krylov Coupling Method for Solving Multi-Physics Nonlinear System of Nuclear Reactor

A Dynamic Forward Mixing Model for Evaluating the Mass Transfer Performances of an Extraction Column

A Dynamic Interaction Mass Transfer Model for Simulating the Mass Transfer Process in Extraction Columns

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

Calculational Methods for Determining the Distribution of Components in A Separation Cascade for Multicomponent Mixture

A projected Newton algorithm based on chemically allowed interval for chemical equilibrium computations

Splitting Approach for Solving Multi-Component Transport Models with Maxwell-Stefan-Diffusion

One-step Algorithm for Multiphase Equilibrium Calculation Using State Equations

An Innovational Jacobian-Split Newton-Krylov Method Combining the Advantages of the Jacobian-Free Newton-Krylov Method and the Finite Difference Jacobian-Based Newton-Krylov Method

Anderson acceleration. Convergence analysis and applications to equilibrium chemistry

An Algorithm for Simultaneous Chemical and Phase-Equilibrium Calculation