On the Simulation of Multiphase Flow Containing Hydrate Particles

Author(s):  
Jian-Kui Zhao ◽  
Jing Gong

The article presents a mathematical model to simulate steady multiphase flow containing hydrate particles. The approach is based on the one-dimensional two-fluid model in which the two phases consist of the gas and the mixture of the liquid and hydrate particles. That is to say the hydrate particle and liquid phase are treated as pseudo-fluid. When a compositional model is used to simulate two-phase flow, it is required to couple mass, momentum, energy equations and equation of state, then the pressure and temperature and mass transfer between gas and liquid phase along the pipeline can be simulated. In the steady condition, it’s assumed that the composition of fluid is unchangeable along the pipeline and the flow can be described by pressure and temperature. In this paper, calculation method of gas-liquid two-phase flow which required a double iterative procedure of pressure and temperature respectively was improved. Liquid holdup and pressure drop were calculated by momentum equation. Enthalpy balance equation was substituted by explicit formulation of temperature calculation which meant that the loop of temperature was not required. So the calculation speed was enhanced. With strict flash calculation the following items were determined: the amount of hydrate; phase number; the location that hydrate appeared; flow rate and molar component of gas phase and liquid phase. Then thermodynamic quantities were calculated with proper relational expressions.

Author(s):  
Shouxu Qiao ◽  
Wenyi Zhong ◽  
Sijia Hao ◽  
Peiyao Qi ◽  
Sichao Tan

Abstract The present study investigates the air-water two-phase flow across a 90-degree vertical-upward elbow with the computational fluid dynamics (CFD) simulation. The Eulerian-Eulerian two-fluid model and the Multi Size Group (MUSIG) model are used to predict the development of the detailed interfacial structures between the two phases. The axial development of the void fraction and the interfacial area concentration are investigated and benchmarked with the experimental data measured using the four-sensor conductivity probe. It is concluded that CFD simulation can predict the characteristics distributions of void fraction and interfacial area concentration and their development downstream of the elbow. The double-peaked void fraction distribution is found to be caused by the secondary flow induced by the elbow. The liquid phase on the outer curvature moves to the inner curvature and forms a double counter rotating vortex, entraining the bubbles to form a double-peaked distribution. The elbow effects become dissipated between 33 and 63 hydraulic diameters. The simulation results of liquid-phase and gas-phase parameters can be used to develop the theoretical two-phase flow models for the elbow region.


2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano

2019 ◽  
Vol 53 (5) ◽  
pp. 1763-1795 ◽  
Author(s):  
Khaled Saleh

This article is the first of two in which we develop a relaxation finite volume scheme for the convective part of the multiphase flow models introduced in the series of papers (Hérard, C.R. Math. 354 (2016) 954–959; Hérard, Math. Comput. Modell. 45 (2007) 732–755; Boukili and Hérard, ESAIM: M2AN 53 (2019) 1031–1059). In the present article we focus on barotropic flows where in each phase the pressure is a given function of the density. The case of general equations of state will be the purpose of the second article. We show how it is possible to extend the relaxation scheme designed in Coquel et al. (ESAIM: M2AN 48 (2013) 165–206) for the barotropic Baer–Nunziato two phase flow model to the multiphase flow model with N – where N is arbitrarily large – phases. The obtained scheme inherits the main properties of the relaxation scheme designed for the Baer–Nunziato two phase flow model. It applies to general barotropic equations of state. It is able to cope with arbitrarily small values of the statistical phase fractions. The approximated phase fractions and phase densities are proven to remain positive and a fully discrete energy inequality is also proven under a classical CFL condition. For N = 3, the relaxation scheme is compared with Rusanov’s scheme, which is the only numerical scheme presently available for the three phase flow model (see Boukili and Hérard, ESAIM: M2AN 53 (2019) 1031–1059). For the same level of refinement, the relaxation scheme is shown to be much more accurate than Rusanov’s scheme, and for a given level of approximation error, the relaxation scheme is shown to perform much better in terms of computational cost than Rusanov’s scheme. Moreover, contrary to Rusanov’s scheme which develops strong oscillations when approximating vanishing phase solutions, the numerical results show that the relaxation scheme remains stable in such regimes.


2001 ◽  
Vol 123 (4) ◽  
pp. 811-818 ◽  
Author(s):  
Jun Ishimoto ◽  
Mamoru Oike ◽  
Kenjiro Kamijo

The two-dimensional characteristics of the vapor-liquid two-phase flow of liquid helium in a pipe are numerically investigated to realize the further development and high performance of new cryogenic engineering applications. First, the governing equations of the two-phase flow of liquid helium based on the unsteady thermal nonequilibrium multi-fluid model are presented and several flow characteristics are numerically calculated, taking into account the effect of superfluidity. Based on the numerical results, the two-dimensional structure of the two-phase flow of liquid helium is shown in detail, and it is also found that the phase transition of the normal fluid to the superfluid and the generation of superfluid counterflow against normal fluid flow are conspicuous in the large gas phase volume fraction region where the liquid to gas phase change actively occurs. Furthermore, it is clarified that the mechanism of the He I to He II phase transition caused by the temperature decrease is due to the deprivation of latent heat for vaporization from the liquid phase. According to these theoretical results, the fundamental characteristics of the cryogenic two-phase flow are predicted. The numerical results obtained should contribute to the realization of advanced cryogenic industrial applications.


Author(s):  
David Heinze ◽  
Thomas Schulenberg ◽  
Lars Behnke

A simulation model for the direct contact condensation of steam in subcooled water is presented that allows determination of major parameters of the process, such as the jet penetration length. Entrainment of water by the steam jet is modeled based on the Kelvin–Helmholtz and Rayleigh–Taylor instability theories. Primary atomization due to acceleration of interfacial waves and secondary atomization due to aerodynamic forces account for the initial size of entrained droplets. The resulting steam-water two-phase flow is simulated based on a one-dimensional two-fluid model. An interfacial area transport equation is used to track changes of the interfacial area density due to droplet entrainment and steam condensation. Interfacial heat and mass transfer rates during condensation are calculated using the two-resistance model. The resulting two-phase flow equations constitute a system of ordinary differential equations, which is solved by means of the explicit Runge–Kutta–Fehlberg algorithm. The simulation results are in good qualitative agreement with published experimental data over a wide range of pool temperatures and mass flow rates.


Author(s):  
F Bakhtar ◽  
H Mashmoushy ◽  
O C Jadayel

During the course of expansion of steam in turbines the fluid first supercools and then nucleates to become a two-phase mixture. The liquid phase consists of a large number of extremely small droplets which are difficult to generate except by nucleation. To reproduce turbine two-phase flow conditions requires a supply of supercooled vapour which can be achieved under blow-down conditions by the equipment employed. This paper is the third of a set describing an investigation into the performance of a cascade of rotor tip section profiles in wet steam and presents the results of the wake traverses.


Author(s):  
Marco Pellegrini ◽  
Giulia Agostinelli ◽  
Hidetoshi Okada ◽  
Masanori Naitoh

Steam condensation is characterized by a relatively large interfacial region between gas and liquid which, in computational fluid dynamic (CFD) analyses, allows the creation of a discretized domain whose average cell size is larger than the interface itself. For this reason generally one fluid model with interface tracking (e.g. volume of fluid method, VOF) is employed for its solution in CFD, since the solution of the interface requires a reasonable amount of cells, reducing the modeling efforts. However, for some particular condensation applications, requiring the computation of long transients or the steam ejected through a large number of holes, one-fluid model becomes computationally too expensive for providing engineering information, and a two-fluid model (i.e. Eulerian two-phase flow) is preferable. Eulerian two-phase flow requires the introduction of closure terms representing the interactions between the two fluids in particular, in the condensation case, drag and heat transfer. Both terms involve the description of the interaction area whose definition is different from the typical one adopted in the boiling analyses. In the present work a simple but effective formulation for the interaction area is given based on the volume fraction gradient and then applied to a validation test case of steam bubbling in various subcooling conditions. It has been shown that this method gives realistic values of bubble detachment time, bubble penetration for the cases of interest in the nuclear application and in the particular application to the Fukushima Daiichi accident.


Sign in / Sign up

Export Citation Format

Share Document