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

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

[1]  C. Kelley.  Solving Nonlinear Equations with Newton's Method , 1987 .

[2]  Vikrant Aute,et al.  A Review on Direct Two-Phase, Phase Change Flow Simulation Methods and their Applications , 2012 .

[3]  H. Deconinck,et al.  Design principles for bounded higher-order convection schemes - a unified approach , 2007, J. Comput. Phys..

[4]  Desmond J. Higham,et al.  Mean-Square and Asymptotic Stability of the Stochastic Theta Method , 2000, SIAM J. Numer. Anal..

[5]  Xiaoying Zhang,et al.  Verification of a Higher-Order Finite Difference Scheme for the One-Dimensional Two-Fluid Model , 2013 .

[6]  J. Castro,et al.  A finite-volume/level-set method for simulating two-phase flows on unstructured grids , 2014 .

[7]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[8]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[9]  Tomasz Kozlowski,et al.  Development of High-Resolution Total Variation Diminishing Scheme for Linear Hyperbolic Problems , 2015, J. Comput. Eng..

[10]  Vincent A. Mousseau,et al.  Implicitly balanced solution of the two-phase flow equations coupled to nonlinear heat conduction , 2004 .

[11]  M. Ishii,et al.  Interfacial drag force in one-dimensional two-fluid model , 2012 .

[12]  Jacques Magnaudet,et al.  An interface-capturing method for incompressible two-phase flows. Validation and application to bubble dynamics , 2007 .

[13]  H. Qin.  COMPARISONS OF PREDICTIONS WITH NUMERICAL BENCHMARK TEST NO. 2.1: FAUCET FLOW , 1992 .

[14]  Ritesh Kumar,et al.  A high resolution total variation diminishing scheme for hyperbolic conservation law and related problems , 2006, Appl. Math. Comput..

[15]  Zbigniew Koszela,et al.  RELAP5 Capability to Predict Pressure Wave Propagation Phenomena in Single- and Two-Phase Flow Conditions , 2012 .

[16]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .