A Jacobian‐free Newton–Krylov method for thermalhydraulics simulations

Summary The current paper is focused on investigating a Jacobian-free Newton–Krylov (JFNK) method to obtain a fully implicit solution for two-phase flows. In the JFNK formulation, the Jacobian matrix is not directly evaluated, potentially leading to major computational savings compared with a simple Newton's solver. The objectives of the present paper are as follows: (i) application of the JFNK method to two-fluid models; (ii) investigation of the advantages and disadvantages of the fully implicit JFNK method compared with commonly used explicit formulations and implicit Newton–Krylov calculations using the determination of the Jacobian matrix; and (iii) comparison of the numerical predictions with those obtained by the Canadian Algorithm for Thermaulhydraulics Network Analysis 4. Two well-known benchmarks are considered, the water faucet and the oscillating manometer. An isentropic two-fluid model is selected. Time discretization is performed using a backward Euler scheme. A Crank–Nicolson scheme is also implemented to check the effect of temporal discretization on the predictions. Advection Upstream Splitting Method+ is applied to the convective fluxes. The source terms are discretized using a central differencing scheme. One explicit and two implicit formulations, one with Newton's solver with the Jacobian matrix and one with JFNK, are implemented. A detailed grid and model parameter sensitivity analysis is performed. For both cases, the JFNK predictions are in good agreement with the analytical solutions and explicit profiles. Further, stable results can be achieved using high CFL numbers up to 200 with a suitable choice of JFNK parameters. The computational time is significantly reduced by JFNK compared with the calculations requiring the determination of the Jacobian matrix. Copyright © 2015 John Wiley & Sons, Ltd.