High-Order Numerical Methods for Compressible Two-Phase Flows

AbstractIn this article a comparison study of the numerical methods for compressible two-phase flows is presented. Although many numerical methods have been developed in recent years to deal with the jump conditions at the fluid-fluid interfaces in compressible multiphase flows, there is a lack of a detailed comparison of these methods. With this regard, the transport five equation model, the modified ghost fluid method and the cut-cell method are investigated here as the typical methods in this field. A variety of numerical experiments are conducted to examine their performance in simulating inviscid compressible two-phase flows. Numerical experiments include Richtmyer-Meshkov instability, interaction between a shock and a rectangle SF6 bubble, Rayleigh collapse of a cylindrical gas bubble in water and shock-induced bubble collapse, involving fluids with small or large density difference. Based on the numerical results, the performance of the method is assessed by the convergence order of the method with respect to interface position, mass conservation, interface resolution and computational efficiency.

Download Full-text

Numerical solution of one-dimensional transient, two-phase flows with temporal fully implicit high order schemes: Subcooled boiling in pipes

Nuclear Engineering and Design ◽

10.1016/j.nucengdes.2016.12.027 ◽

2017 ◽

Vol 313 ◽

pp. 319-329 ◽

Cited By ~ 3

Author(s):

R. López ◽

A. Lecuona ◽

J. Nogueira ◽

C. Vereda

Keyword(s):

Numerical Solution ◽

High Order ◽

Subcooled Boiling ◽

Two Phase Flows ◽

Two Phase ◽

One Dimensional ◽

High Order Schemes ◽

Fully Implicit

Download Full-text

Accurate Treatment of Material Interface Dynamics in the Calculation of One-Dimensional Two-Phase Flows by the Integral Method of Characteristics

Journal of Pressure Vessel Technology ◽

10.1115/1.3264341 ◽

1984 ◽

Vol 106 (3) ◽

pp. 261-269 ◽

Cited By ~ 1

Author(s):

Y. W. Shin ◽

A. H. Wiedermann

Keyword(s):

Numerical Methods ◽

Moving Boundary ◽

Critical Evaluation ◽

Special Technique ◽

Flow Conditions ◽

Two Phase Flows ◽

Two Phase ◽

Material Interface ◽

One Dimensional ◽

Accurate Treatment

Accurate numerical methods for treating the junction and boundary conditions needed in the transient two-phase flows of a piping network were published earlier by us; the same methods are used to formulate the treatment of the material interface as a moving boundary. The method formulated is used in a computer program to calculate sample problems designed to test the numerical methods as to their ability and the accuracy limits for calculation of the transient two-phase flows in the piping network downstream of a PWR pressurizer. Independent exact analytical solutions for the sample problems are used as the basis of a critical evaluation of the proposed numerical methods. The evaluation revealed that the proposed boundary scheme indeed generates very accurate numerical results. However, in some extreme flow conditions, numerical difficulties were experienced that eventually led to numerical instability. This paper discusses further a special technique to overcome the difficulty.

Download Full-text