### 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.

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 -