skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model

Journal Article · · Progress in Nuclear Energy

This study concerns the application and solver robustness of the Newton-Krylov method in solving two-phase flow drift-flux model problems using high-order numerical schemes. In our previous studies, the Newton-Krylov method has been proven as a promising solver for two-phase flow drift-flux model problems. However, these studies were limited to use first-order numerical schemes only. Moreover, the previous approach to treating the drift-flux closure correlations was later revealed to cause deteriorated solver convergence performance, when the mesh was highly refined, and also when higher-order numerical schemes were employed. In this study, a second-order spatial discretization scheme that has been tested with two-fluid two-phase flow model was extended to solve drift-flux model problems. In order to improve solver robustness, and therefore efficiency, a new approach was proposed to treating the mean drift velocity of the gas phase as a primary nonlinear variable to the equation system. With this new approach, significant improvement in solver robustness was achieved. With highly refined mesh, the proposed treatment along with the Newton-Krylov solver were extensively tested with two-phase flow problems that cover a wide range of thermal-hydraulics conditions. Satisfactory convergence performances were observed for all test cases. Numerical verification was then performed in the form of mesh convergence studies, from which expected orders of accuracy were obtained for both the first-order and the second-order spatial discretization schemes. Finally, the drift-flux model, along with numerical methods presented, were validated with three sets of flow boiling experiments that cover different flow channel geometries (round tube, rectangular tube, and rod bundle), and a wide range of test conditions (pressure, mass flux, wall heat flux, inlet subcooling and outlet void fraction).

Research Organization:
Idaho National Lab. (INL), Idaho Falls, ID (United States)
Sponsoring Organization:
USDOE Office of Nuclear Energy (NE)
Grant/Contract Number:
AC07-05ID14517
OSTI ID:
1375246
Alternate ID(s):
OSTI ID: 1549808
Report Number(s):
INL/JOU-17-40865; PII: S0149197017301750
Journal Information:
Progress in Nuclear Energy, Vol. 100; ISSN 0149-1970
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

References (20)

Development and assessment of system analysis code, TASS/SMR for integral reactor, SMART journal March 2012
Applicability of TASS/SMR using drift flux model for SMART LOCA analysis journal September 2013
A new numerical method for solution of boiling flow using combination of SIMPLE and Jacobian-free Newton-Krylov algorithms journal March 2017
An advanced new fully implicit numerical method for two-phase flow subchannel analysis based on the Drift Flux Model journal October 2017
Benchmarking a sub-channel program based on a drift-flux model with 8×8 NUPEC BWR rod bundle journal August 2013
Sub-channel analysis of 8×8 and 9×9 BWR fuel assemblies with different two-phase flow models journal December 2013
Subchannel analysis in BWR fuel bundles journal January 1985
Jacobian-free Newton–Krylov methods: a survey of approaches and applications journal January 2004
Implicitly balanced solution of the two-phase flow equations coupled to nonlinear heat conduction journal October 2004
A Fully Implicit Hybrid Solution Method for a Two-Phase Thermal-Hydraulic Model journal May 2005
Study on the stability behaviour of two-phase natural circulation systems using a four-equation drift flux model journal February 2007
A numerical technique for analysis of transient two-phase flow in a vertical tube using the Drift Flux Model journal January 2012
Numerical study of nuclear coupled two-phase flow instability in natural circulation system under low pressure and low quality journal December 2011
Applications of high-resolution spatial discretization scheme and Jacobian-free Newton–Krylov method in two-phase flow problems journal September 2015
Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method journal January 2016
Application of Jacobian-free Newton–Krylov method in implicitly solving two-fluid six-equation two-phase flow problems: Implementation, validation and benchmark journal April 2016
Implicitly solving phase appearance and disappearance problems using two-fluid six-equation model journal April 2016
Solving phase appearance/disappearance two-phase flow problems with high resolution staggered grid and fully implicit schemes by the Jacobian-free Newton–Krylov Method journal April 2016
Numerical study on the Welander oscillatory natural circulation problem using high-order numerical methods journal January 2017
Average Volumetric Concentration in Two-Phase Flow Systems journal November 1965

Cited By (1)

A Generalized Reduced-Order Dynamic Model for Two-Phase Flow in Pipes journal June 2019