Application of implicit Roe-type scheme and Jacobian-Free Newton-Krylov method to two-phase flow problems

Guojun Hu, Tomasz Kozlowski

Research output: Contribution to journalArticle

Abstract

A new implicit second-order accurate spatial scheme for steady-state thermal-hydraulic simulations of the two-phase two-fluid six-equation model is proposed. The new scheme is based on a Roe-type numerical flux that is formulated with the help of analytical approximate eigenvalues and eigenvectors of the two-phase system. Approximate eigenvalues and eigenvectors are obtained with a structured Jacobian matrix that is general for arbitrary Equation of State (EOS). Three issues are solved in this article: (1) A fully-implicit scheme using the backward Euler method is proposed to improve the stability and avoid the time step limit of the originally explicit scheme. The fully implicit scheme is solved with the Jacobian-Free Newton-Krylov method. (2) A second-order accurate spatial scheme is proposed by converting the first-order Roe-type numerical flux to a second-order one. The conversion is made by extending an existing procedure for single-phase flows to two-phase flows. (3) Phase appearance issue is treated with a numerical procedure, which requires little modification to the schemes. The new scheme is verified and validated with two numerical tests: the faucet flow problem and the BFBT benchmark.

Original languageEnglish (US)
Pages (from-to)180-190
Number of pages11
JournalAnnals of Nuclear Energy
Volume119
DOIs
StatePublished - Sep 2018
Externally publishedYes

Fingerprint

Newton-Raphson method
Eigenvalues and eigenfunctions
Two phase flow
Fluxes
Jacobian matrices
Equations of state
Hydraulics
Fluids
Hot Temperature

Keywords

  • Implicit scheme
  • JFNK
  • Riemann solver
  • Second-order
  • Two-phase flow

ASJC Scopus subject areas

  • Nuclear Energy and Engineering

Cite this

Application of implicit Roe-type scheme and Jacobian-Free Newton-Krylov method to two-phase flow problems. / Hu, Guojun; Kozlowski, Tomasz.

In: Annals of Nuclear Energy, Vol. 119, 09.2018, p. 180-190.

Research output: Contribution to journalArticle

@article{792cb09d867349a79460ba93d5cbdc14,
title = "Application of implicit Roe-type scheme and Jacobian-Free Newton-Krylov method to two-phase flow problems",
abstract = "A new implicit second-order accurate spatial scheme for steady-state thermal-hydraulic simulations of the two-phase two-fluid six-equation model is proposed. The new scheme is based on a Roe-type numerical flux that is formulated with the help of analytical approximate eigenvalues and eigenvectors of the two-phase system. Approximate eigenvalues and eigenvectors are obtained with a structured Jacobian matrix that is general for arbitrary Equation of State (EOS). Three issues are solved in this article: (1) A fully-implicit scheme using the backward Euler method is proposed to improve the stability and avoid the time step limit of the originally explicit scheme. The fully implicit scheme is solved with the Jacobian-Free Newton-Krylov method. (2) A second-order accurate spatial scheme is proposed by converting the first-order Roe-type numerical flux to a second-order one. The conversion is made by extending an existing procedure for single-phase flows to two-phase flows. (3) Phase appearance issue is treated with a numerical procedure, which requires little modification to the schemes. The new scheme is verified and validated with two numerical tests: the faucet flow problem and the BFBT benchmark.",
keywords = "Implicit scheme, JFNK, Riemann solver, Second-order, Two-phase flow",
author = "Guojun Hu and Tomasz Kozlowski",
year = "2018",
month = "9",
doi = "10.1016/j.anucene.2018.05.003",
language = "English (US)",
volume = "119",
pages = "180--190",
journal = "Annals of Nuclear Energy",
issn = "0306-4549",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - Application of implicit Roe-type scheme and Jacobian-Free Newton-Krylov method to two-phase flow problems

AU - Hu, Guojun

AU - Kozlowski, Tomasz

PY - 2018/9

Y1 - 2018/9

N2 - A new implicit second-order accurate spatial scheme for steady-state thermal-hydraulic simulations of the two-phase two-fluid six-equation model is proposed. The new scheme is based on a Roe-type numerical flux that is formulated with the help of analytical approximate eigenvalues and eigenvectors of the two-phase system. Approximate eigenvalues and eigenvectors are obtained with a structured Jacobian matrix that is general for arbitrary Equation of State (EOS). Three issues are solved in this article: (1) A fully-implicit scheme using the backward Euler method is proposed to improve the stability and avoid the time step limit of the originally explicit scheme. The fully implicit scheme is solved with the Jacobian-Free Newton-Krylov method. (2) A second-order accurate spatial scheme is proposed by converting the first-order Roe-type numerical flux to a second-order one. The conversion is made by extending an existing procedure for single-phase flows to two-phase flows. (3) Phase appearance issue is treated with a numerical procedure, which requires little modification to the schemes. The new scheme is verified and validated with two numerical tests: the faucet flow problem and the BFBT benchmark.

AB - A new implicit second-order accurate spatial scheme for steady-state thermal-hydraulic simulations of the two-phase two-fluid six-equation model is proposed. The new scheme is based on a Roe-type numerical flux that is formulated with the help of analytical approximate eigenvalues and eigenvectors of the two-phase system. Approximate eigenvalues and eigenvectors are obtained with a structured Jacobian matrix that is general for arbitrary Equation of State (EOS). Three issues are solved in this article: (1) A fully-implicit scheme using the backward Euler method is proposed to improve the stability and avoid the time step limit of the originally explicit scheme. The fully implicit scheme is solved with the Jacobian-Free Newton-Krylov method. (2) A second-order accurate spatial scheme is proposed by converting the first-order Roe-type numerical flux to a second-order one. The conversion is made by extending an existing procedure for single-phase flows to two-phase flows. (3) Phase appearance issue is treated with a numerical procedure, which requires little modification to the schemes. The new scheme is verified and validated with two numerical tests: the faucet flow problem and the BFBT benchmark.

KW - Implicit scheme

KW - JFNK

KW - Riemann solver

KW - Second-order

KW - Two-phase flow

UR - http://www.scopus.com/inward/record.url?scp=85046666348&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85046666348&partnerID=8YFLogxK

U2 - 10.1016/j.anucene.2018.05.003

DO - 10.1016/j.anucene.2018.05.003

M3 - Article

AN - SCOPUS:85046666348

VL - 119

SP - 180

EP - 190

JO - Annals of Nuclear Energy

JF - Annals of Nuclear Energy

SN - 0306-4549

ER -