Abstract
We propose a new solution method to the two-phase two-fluid six-equation model, which involves six conservation equations. A Roe-type numerical flux is formulated based on a very structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the eigenvalue and eigenvector can be obtained analytically. An explicit Roe-type numerical solver is constructed based on the analytical eigenvalue and eigenvector. A critical feature of the method is that the formulation of our solver does not depend on the form of the equation of state. The proposed method is applicable to realistic two-phase problems. We applied the numerical solver to the BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark. Considering simplified physical models, the solutions are in very good agreement with solutions from either existing codes and experiment data. The numerical solver using analytical eigenvalue and eigenvector is shown to be stable and robust.
| Original language | English (US) |
|---|---|
| State | Published - Jan 1 2017 |
| Event | 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics, NURETH 2017 - Xi'an, Shaanxi, China Duration: Sep 3 2017 → Sep 8 2017 |
Other
| Other | 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics, NURETH 2017 |
|---|---|
| Country | China |
| City | Xi'an, Shaanxi |
| Period | 9/3/17 → 9/8/17 |
Fingerprint
Keywords
- BFBT
- Reactor safety
- Riemann solver
- Two-phase flow
ASJC Scopus subject areas
- Nuclear Energy and Engineering
- Instrumentation
Cite this
On the approximate riemann solver for the two-phase two-fluid six equation model and application to real system. / Hu, Guojun; Kozlowski, Tomasz.
2017. Paper presented at 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics, NURETH 2017, Xi'an, Shaanxi, China.Research output: Contribution to conference › Paper
}
TY - CONF
T1 - On the approximate riemann solver for the two-phase two-fluid six equation model and application to real system
AU - Hu, Guojun
AU - Kozlowski, Tomasz
PY - 2017/1/1
Y1 - 2017/1/1
N2 - We propose a new solution method to the two-phase two-fluid six-equation model, which involves six conservation equations. A Roe-type numerical flux is formulated based on a very structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the eigenvalue and eigenvector can be obtained analytically. An explicit Roe-type numerical solver is constructed based on the analytical eigenvalue and eigenvector. A critical feature of the method is that the formulation of our solver does not depend on the form of the equation of state. The proposed method is applicable to realistic two-phase problems. We applied the numerical solver to the BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark. Considering simplified physical models, the solutions are in very good agreement with solutions from either existing codes and experiment data. The numerical solver using analytical eigenvalue and eigenvector is shown to be stable and robust.
AB - We propose a new solution method to the two-phase two-fluid six-equation model, which involves six conservation equations. A Roe-type numerical flux is formulated based on a very structured Jacobian matrix. We formulate the Jacobian matrix with arbitrary equation of state and simplify the Jacobian matrix with the help of a few auxiliary variables, e.g. isentropic speed of sound. Because the Jacobian matrix is very structured, the eigenvalue and eigenvector can be obtained analytically. An explicit Roe-type numerical solver is constructed based on the analytical eigenvalue and eigenvector. A critical feature of the method is that the formulation of our solver does not depend on the form of the equation of state. The proposed method is applicable to realistic two-phase problems. We applied the numerical solver to the BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark. Considering simplified physical models, the solutions are in very good agreement with solutions from either existing codes and experiment data. The numerical solver using analytical eigenvalue and eigenvector is shown to be stable and robust.
KW - BFBT
KW - Reactor safety
KW - Riemann solver
KW - Two-phase flow
UR - http://www.scopus.com/inward/record.url?scp=85052380856&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85052380856&partnerID=8YFLogxK
M3 - Paper
AN - SCOPUS:85052380856
ER -