Notes on the implementation of a fully-implicit numerical scheme for a two-phase three-field flow model

Frepoli, Cesare; Mahaffy, J.H.; Ohkawa, Katsuhiro

doi:10.1016/s0029-5493(03)00159-6

Cited by 49 publications

(19 citation statements)

References 25 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The JFNK method is robust and efficient even without implementing the analytical Jacobian matrix, which is typically required for traditional Newton's methods (Frepoli et al, 2003). The Jacobian-free feature yields significant savings in developing a fully implicit system analysis code.…”

Section: Discussionmentioning

confidence: 99%

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

Zou

Zhao

Zhang

2016

Annals of Nuclear Energy

View full text Add to dashboard Cite

a b s t r a c tThis paper 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.

show abstract

Section: Discussionmentioning

confidence: 99%

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

Zou

Zhao

Zhang

2016

Annals of Nuclear Energy

View full text Add to dashboard Cite

show abstract

“…A transient two-fluid, three-field model [16,17] is adopted in the CUPID code to model the thermal-hydraulics of PWRs, where the two fluids are vapor and liquid, and the three fields are vapor, continuous liquid and droplet fields. A major application of the CUPID code is the analysis of transient two-phase flow with a resolution ranging from the CFD scale to the component scale.…”

Section: Mathematical Models For Two-phase Flowsmentioning

confidence: 99%

Multi-Scale Thermal-Hydraulic Analysis of PWRS Using the Cupid Code

Yoon¹,

Cho²,

Lee³

et al. 2012

Nuclear Engineering and Technology

View full text Add to dashboard Cite

“…The three fields include a continuous liquid, an entrained liquid, and a vapor field [15,17]. The three-field formulation of a two-phase flow is a straightforward extension of the two-fluid single pressure model [13], in which the local instant formulation of a two-phase flow is time-averaged.…”

Section: Basic Governing Equationsmentioning

confidence: 99%