skip to main content


Title: Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method

This study presents a numerical investigation on using the Jacobian-free Newton–Krylov (JFNK) method to solve the two-phase flow four-equation drift flux model with realistic constitutive correlations (‘closure models’). The drift flux model is based on Isshi and his collaborators’ work. Additional constitutive correlations for vertical channel flow, such as two-phase flow pressure drop, flow regime map, wall boiling and interfacial heat transfer models, were taken from the RELAP5-3D Code Manual and included to complete the model. The staggered grid finite volume method and fully implicit backward Euler method was used for the spatial discretization and time integration schemes, respectively. The Jacobian-free Newton–Krylov method shows no difficulty in solving the two-phase flow drift flux model with a discrete flow regime map. In addition to the Jacobian-free approach, the preconditioning matrix is obtained by using the default finite differencing method provided in the PETSc package, and consequently the labor-intensive implementation of complex analytical Jacobian matrix is avoided. Extensive and successful numerical verification and validation have been performed to prove the correct implementation of the models and methods. Code-to-code comparison with RELAP5-3D has further demonstrated the successful implementation of the drift flux model.
 [1] ;  [1] ;  [1]
  1. Idaho National Lab. (INL), Idaho Falls, ID (United States)
Publication Date:
Report Number(s):
Journal ID: ISSN 0306-4549; PII: S0306454915003977; TRN: US1700088
Grant/Contract Number:
Accepted Manuscript
Journal Name:
Annals of Nuclear Energy (Oxford)
Additional Journal Information:
Journal Name: Annals of Nuclear Energy (Oxford); Journal Volume: 87; Journal Issue: P2; Journal ID: ISSN 0306-4549
Research Org:
Idaho National Lab. (INL), Idaho Falls, ID (United States)
Sponsoring Org:
Country of Publication:
United States
42 ENGINEERING; 97 MATHEMATICS AND COMPUTING; drift flux model; Jacobian-free Newton–Krylov method; verification; validation
OSTI Identifier:
Alternate Identifier(s):
OSTI ID: 1396942