JFNK method with a physics-based preconditioner for the fully implicit solution of one-dimensional drift-flux model in boiling two-phase flow

The Jacobian-free Newton-Krylov (JFNK) method with an efficient physics-based preconditioner is applied for the numerical solution of the one-dimensional drift-flux model with the closure constitutive equations. Additional closure correlations, including the flow pattern dependent heat transfer correlations and the flow pattern independent kinematic constitutive correlation, are used to close the governing equations of one-dimensional drift-flux model. The governing equations have been discretized using the first-order upwind method for spatial discretization and the fully implicit method for temporal discretization. An efficient physics-based preconditioner derived from the semi-implicit solution of governing equations is used in the JFNK method to improve the efficiency and numerical stability. The numerical verification and code validation have been performed for subcooled boiling two-phase flow in a vertical tube. By comparing with the other methods (JFNK method without preconditioner, the Broyden method and the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method), the preconditioning JFNK method shows the robustness and a good computational efficiency. Moreover, by comparing the numerical simulation results with the experimental results, it is found that the JFNK method with a physics-based preconditioner shows the good accuracy for the numerical simulation for one-dimensional boiling two-phase flow.