skip to main content

DOE PAGESDOE PAGES

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

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 ofmore » 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).« less
Authors:
ORCiD logo [1] ;  [1] ;  [1]
  1. Idaho National Lab. (INL), Idaho Falls, ID (United States)
Publication Date:
Report Number(s):
INL/JOU-17-40865
Journal ID: ISSN 0149-1970; PII: S0149197017301750
Grant/Contract Number:
AC07-05ID14517
Type:
Accepted Manuscript
Journal Name:
Progress in Nuclear Energy
Additional Journal Information:
Journal Volume: 100; Journal ID: ISSN 0149-1970
Publisher:
Elsevier
Research Org:
Idaho National Lab. (INL), Idaho Falls, ID (United States)
Sponsoring Org:
USDOE Office of Nuclear Energy (NE)
Country of Publication:
United States
Language:
English
Subject:
22 GENERAL STUDIES OF NUCLEAR REACTORS; Newton-Krylov method; drift-flux model; high-order numerical methods
OSTI Identifier:
1375246

Zou, Ling, Zhao, Haihua, and Zhang, Hongbin. Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model. United States: N. p., Web. doi:10.1016/j.pnucene.2017.07.008.
Zou, Ling, Zhao, Haihua, & Zhang, Hongbin. Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model. United States. doi:10.1016/j.pnucene.2017.07.008.
Zou, Ling, Zhao, Haihua, and Zhang, Hongbin. 2017. "Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model". United States. doi:10.1016/j.pnucene.2017.07.008. https://www.osti.gov/servlets/purl/1375246.
@article{osti_1375246,
title = {Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model},
author = {Zou, Ling and Zhao, Haihua and Zhang, Hongbin},
abstractNote = {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).},
doi = {10.1016/j.pnucene.2017.07.008},
journal = {Progress in Nuclear Energy},
number = ,
volume = 100,
place = {United States},
year = {2017},
month = {8}
}