### 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*. Paper presented at 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics, NURETH 2017, Xi'an, Shaanxi, China.

**On the approximate riemann solver for the two-phase two-fluid six equation model and application to real system.** / Hu, Guojun; Kozlowski, Tomasz.

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 -