A solver for the two-phase two-fluid model based on high-resolution total variation diminishing scheme

Finite volume techniques with staggered mesh are used to develop a new numerical solver for the one-dimensional two-phase two-fluid model using a high-resolution, Total Variation Diminishing (TVD) scheme. The solver is implemented to analyze numerical benchmark problems for verification and testing its abilities to handle discontinuities and fast transients with phase change. Convergence rates are investigated by comparing numerical results to analytical solutions available in literature for the case of the faucet flow problem. The solver based on a new TVD scheme is shown to exhibit higher-order of accuracy compared to other numerical schemes. Mass errors are also examined when phase change occurs for the shock tube problem, and compared to those of the 1st-order upwind scheme implemented in the nuclear thermal-hydraulics code TRACE. The solver is shown to exhibit numerical stability when applied to problems with discontinuous solutions and results of the new solver are free of spurious oscillations.