Validation of 2D CFD for Two-Phase Transient Flow in a Channel and Comparison With 1D Model

Author(s):  
Stamatis Kalogerakos ◽  
Mustapha Gourma ◽  
Chris Thompson

Severe limitations of the use of three-dimensional computational fluid dynamics codes (CFD) arise when trying to simulate multiphase flow in long pipes due to time constraints. 1D codes for two-phase flow, based on two-fluid models, are fast but are known to be accurate only when the velocities are within the Kelvin-Helmholtz inviscid limit [1]. An alternative is to carry out a two-dimensional CFD simulation of a channel based on the Volume of Fluid (VOF) model. 2D CFD has a wider applicability range compared to 1D, it does not have the issue of ill-posedness and it also has better turbulence models built in. Again compared to 1D the 2D VOF model has a better interface description and wall treatment. In this paper a novel method is introduced that allows swift simulations of pipeline two-phase flow in the stratified and slug flow regime, by approximating the pipe as a channel and with a methodology that solves the problem of the interfacial velocity differences, inherent in the volume of fluid model. An initial validation using the wave growth problem has already been carried out [2]. Here a set consisting of 92 experimental cases in the slug flow regime has been simulated with 2D CFD, and the simulation results showed a good agreement with experimental results. Discussions in the paper include also the question of the range of applicability for 2D CFD, and the advantages and disadvantages compared to 3D CFD and also to 1D code based on the two-fluid model. Shear stresses are then extracted from the 2D CFD simulations and used to recalibrate the friction factors [3] used in the 1D code.

1994 ◽  
Vol 59 (12) ◽  
pp. 2595-2603
Author(s):  
Lothar Ebner ◽  
Marie Fialová

Two regions of instabilities in horizontal two-phase flow were detected. The first was found in the transition from slug to annular flow, the second between stratified and slug flow. The existence of oscillations between the slug and annular flows can explain the differences in the limitation of the slug flow in flow regime maps proposed by different authors. Coexistence of these two regimes is similar to bistable behaviour of some differential equation solutions.


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

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):  
Yoshiteru Komuro ◽  
Atsushi Kodama ◽  
Yoshiyuki Kondo ◽  
Koichi Tanimoto ◽  
Takashi Hibiki

Abstract Two-phase flows are observed in various industrial plants and piping systems. Understanding two-phase flow behaviors such as flow patterns and unsteady void fraction in horizontal and vertical pipes are crucial in improving plant safety. Notably, the flow patterns observed in a large diameter pipe (approx. 4–6 in or larger) are significantly different from those observed in a medium diameter pipe. In a vertical large diameter pipe, no slug flow is observed due to the instantaneous slug bubble breakup caused by the surface instability. Besides, in a horizontal pipe, flow regime transition from stratification of liquid and gas to slug (plug) flow that induces unsteady flow should be taken into account. From this viewpoint, it is necessary to predict the flow regime in horizontal and vertical large diameter pipes with some elbows and to evaluate the unsteady flow regime. In this study, the simulation method based on the two-fluid model is developed. The two-fluid model is considered the most accurate model because the governing equations for mass, momentum, and energy transfer are formulated for each phase. When using the two-fluid model, some constitutive equations should be given in computing the momentum transfer between gas and liquid phases. In this study, several state-of-art constitutive equations of the bubble diameter, the interfacial drag force and non-drag forces such as the lift force and the bubble-bubble collision force, are implemented in the platform of ANSYS FLUENT. The developed simulation method is validated with visualization results and force from an air-water flow at the elbow of the piping system.


2004 ◽  
Vol 45 (10) ◽  
pp. 1049-1066 ◽  
Author(s):  
Moon-Sun Chung ◽  
Seung-Kyung Pak ◽  
Keun-Shik Chang

SPE Journal ◽  
2012 ◽  
Vol 17 (03) ◽  
pp. 942-955 ◽  
Author(s):  
Mahdy Shirdel ◽  
Kamy Sepehrnoori

Summary A great deal of research has been focused on transient two-phase flow in wellbores. However, there is lack of a comprehensive two-fluid model in the literature. In this paper, we present an implementation of a pseudo-compositional, thermal, fully implicit, transient two-fluid model for two-phase flow in wellbores. In this model, we solve gas/liquid mass balance, gas/liquid momentum balance, and two-phase energy balance equations to obtain five primary variables: liquid velocity, gas velocity, pressure, holdup, and temperature. This simulator can be used as a stand-alone code or can be used in conjunction with a reservoir simulator to mimic wellbore/reservoir dynamic interactions. In our model, we consider stratified, bubbly, intermittent, and annular flow regimes using appropriate closure relations for interphase and wall-shear stress terms in the momentum equations. In our simulation, we found that the interphase and wall-shear stress terms for different flow regimes can significantly affect the model's results. In addition, the interphase momentum transfer terms mainly influence the holdup value. The outcome of this research leads to a more accurate simulation of multiphase flow in the wellbore and pipes, which can be applied to the surface facility design, well-performance optimization, and wellbore damage estimation.


1997 ◽  
Vol 63 (611) ◽  
pp. 2377-2385
Author(s):  
Kiyoshi MINEMURA ◽  
Tomomi UCHIYAMA ◽  
Katsuhiko KINOSHITA ◽  
Lin LYU ◽  
Shinji SYODA ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document