Abstract
In the nuclear industry, a method based on a staggered grid is used in two-phase flow system codes such as RELAP, TRAC, and CATHARE. Solving the two-phase two-fluid model with this method is complicated. The objective of this article is to develop a new solver, which is mathematically consistent and algebraically simpler than existing codes. The extension of existing shock-capturing upwind schemes for single-phase flows is our way. A numerical solver with a Roe-type numerical flux is formulated based on a very well-structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix to a simple and structured form with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the characteristic polynomial of the Jacobian matrix is simple and suitable for analytical analysis. Results from the characteristic analysis of the two-phase system are consistent with well-known facts, such as the ill-posedness of the basic two-phase two-fluid model which assumes all pressure terms are equal. An explicit numerical solver, with a Roe-type numerical flux, is constructed based on the characteristic analysis. A critical feature of the method is that the formulation does not depend on the form of equation of state and the method is applicable to realistic two-phase problems. We demonstrate solver performance based on three two-phase benchmark problems: two-phase shock-tube problem, faucet flow problem, and Christensen boiling pipe problem. The solutions are in excellent agreement with analytical solutions and numerical solutions from a system code. The new solver provides essential framework for developing a more accurate and robust solver for realistic reactor safety analysis. However, improvements on the new solver is necessary for achieving a high-order accuracy and increasing the robustness.
Original language | English (US) |
---|---|
Pages (from-to) | 354-370 |
Number of pages | 17 |
Journal | Nuclear Engineering and Design |
Volume | 326 |
DOIs | |
State | Published - Jan 2018 |
Fingerprint
Keywords
- Hyperbolicity
- Realistic EOS
- Riemann solver
- Two-phase flow
ASJC Scopus subject areas
- Nuclear and High Energy Physics
- Nuclear Energy and Engineering
- Materials Science(all)
- Safety, Risk, Reliability and Quality
- Waste Management and Disposal
- Mechanical Engineering
Cite this
A roe-type numerical solver for the two-phase two-fluid six-equation model with realistic equation of state. / Hu, Guojun; Kozlowski, Tomasz.
In: Nuclear Engineering and Design, Vol. 326, 01.2018, p. 354-370.Research output: Contribution to journal › Article
}
TY - JOUR
T1 - A roe-type numerical solver for the two-phase two-fluid six-equation model with realistic equation of state
AU - Hu, Guojun
AU - Kozlowski, Tomasz
PY - 2018/1
Y1 - 2018/1
N2 - In the nuclear industry, a method based on a staggered grid is used in two-phase flow system codes such as RELAP, TRAC, and CATHARE. Solving the two-phase two-fluid model with this method is complicated. The objective of this article is to develop a new solver, which is mathematically consistent and algebraically simpler than existing codes. The extension of existing shock-capturing upwind schemes for single-phase flows is our way. A numerical solver with a Roe-type numerical flux is formulated based on a very well-structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix to a simple and structured form with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the characteristic polynomial of the Jacobian matrix is simple and suitable for analytical analysis. Results from the characteristic analysis of the two-phase system are consistent with well-known facts, such as the ill-posedness of the basic two-phase two-fluid model which assumes all pressure terms are equal. An explicit numerical solver, with a Roe-type numerical flux, is constructed based on the characteristic analysis. A critical feature of the method is that the formulation does not depend on the form of equation of state and the method is applicable to realistic two-phase problems. We demonstrate solver performance based on three two-phase benchmark problems: two-phase shock-tube problem, faucet flow problem, and Christensen boiling pipe problem. The solutions are in excellent agreement with analytical solutions and numerical solutions from a system code. The new solver provides essential framework for developing a more accurate and robust solver for realistic reactor safety analysis. However, improvements on the new solver is necessary for achieving a high-order accuracy and increasing the robustness.
AB - In the nuclear industry, a method based on a staggered grid is used in two-phase flow system codes such as RELAP, TRAC, and CATHARE. Solving the two-phase two-fluid model with this method is complicated. The objective of this article is to develop a new solver, which is mathematically consistent and algebraically simpler than existing codes. The extension of existing shock-capturing upwind schemes for single-phase flows is our way. A numerical solver with a Roe-type numerical flux is formulated based on a very well-structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix to a simple and structured form with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the characteristic polynomial of the Jacobian matrix is simple and suitable for analytical analysis. Results from the characteristic analysis of the two-phase system are consistent with well-known facts, such as the ill-posedness of the basic two-phase two-fluid model which assumes all pressure terms are equal. An explicit numerical solver, with a Roe-type numerical flux, is constructed based on the characteristic analysis. A critical feature of the method is that the formulation does not depend on the form of equation of state and the method is applicable to realistic two-phase problems. We demonstrate solver performance based on three two-phase benchmark problems: two-phase shock-tube problem, faucet flow problem, and Christensen boiling pipe problem. The solutions are in excellent agreement with analytical solutions and numerical solutions from a system code. The new solver provides essential framework for developing a more accurate and robust solver for realistic reactor safety analysis. However, improvements on the new solver is necessary for achieving a high-order accuracy and increasing the robustness.
KW - Hyperbolicity
KW - Realistic EOS
KW - Riemann solver
KW - Two-phase flow
UR - http://www.scopus.com/inward/record.url?scp=85035768168&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85035768168&partnerID=8YFLogxK
U2 - 10.1016/j.nucengdes.2017.11.023
DO - 10.1016/j.nucengdes.2017.11.023
M3 - Article
AN - SCOPUS:85035768168
VL - 326
SP - 354
EP - 370
JO - Nuclear Engineering and Design
JF - Nuclear Engineering and Design
SN - 0029-5493
ER -