CFD modelling of slug flow in vertical tubes
Introduction
When gas and liquid flow in a pipe they tend to distribute themselves in a variety of configurations. These characteristic distributions of the fluid–fluid interface are called flow patterns or flow regimes. Much time and effort has been expended in determining these regimes for various pairs of fluids, channel geometries, and inclinations (Mandhane et al., 1974, Taitel and Duckler, 1976, Taitel, 1986, Barnea, 1987). For vertical co-current flow and at low gas flow rates, the flow pattern observed is bubbly. Here, the gas phase is distributed as discrete bubbles within the liquid continuum. At higher gas flow rates, some of the bubbles have nearly the same cross-sectional area as that of the channel. These bullet-shaped bubbles—sometimes referred to as ‘Taylor bubbles’ or ‘slugs’—move along and are separated by liquid plugs that may or may not contain a dispersion of smaller gas bubbles. An increase in the gas flow rate in a two-phase mixture flowing in slug flow will eventually result in a complete destruction of slug flow integrity with consequential churning or oscillatory action. At very high gas flow rates, the flow becomes annular, in which, adjacent to the channel wall. There is a liquid continuum and the core of the channel is a gas continuum (Hewitt and Hall-Taylor, 1970).
Slug flow is the most important of the two-phase flow regimes primarily because of the numerous industrial and practical applications. Some of these include buoyancy-driven fermenters, production and transportation of hydrocarbons, boiling and condensation processes in thermal power plants, and emergency cooling of nuclear reactors. Slug flow is characterised by its random intermittence and inherent unsteadiness. A fixed observer would see a quasi-periodic occurrence of long, bullet-shaped Taylor bubbles followed by liquid plugs sometimes carrying dispersed bubbles and taking on the appearance of bubbly flow in a pipe. The bubble region may take on stratified or annular configurations depending upon the tube inclination and flow conditions. In order to understand the complex features of intermittent slug flow, mainly experimental research has been conducted to study the motion of isolated Taylor bubbles in motionless and flowing liquids for various inclination angles. Here, vertical slug flow will be emphasised.
The shape and the velocity with which a single Taylor bubble ascending through a denser stagnant liquid is influenced by the forces acting on it, namely the viscous, inertial and interfacial forces. Dimensionless analysis based on Pi-theorem leads to the following dimensionless groups:where is the diameter of the tube, is Taylor bubble velocity, and are the density of the liquid and gas, respectively, and are the viscosity of the liquid and gas respectively, is the surface tension and is the length of Taylor bubble. For cylindrical bubbles, the film thickness and the bubble rise velocity are independent of the bubble length (Griffith and Wallis, 1961, Nicklin et al., 1962, Mao and Duckler, 1989, Polonsky et al., 1999). Under the assumption that inertial forces in the gas are far smaller than the inertial forces in the liquid (), can be eliminated. And if the viscosity of the gas in the bubble is neglected the following set of three dimensionless groups is sufficient to characterise the motion of a single bubble rising through a motionless liquid. These are the Eötvös number , the Morton number defined as and the Froude number defined as . The Eötvös number represents the relative significance of surface tension and buoyancy. The Morton number is sometimes referred to as the property group. The Froude number represents the ratio of inertial to gravitational forces. Other dimensionless groups could be used, e.g., Fabre and Line (1992) in their review paper on the motion of Taylor bubbles used Froude number as a unique function of Eötvös number and a dimensionless inverse viscosity number, , given by . Wallis (1969) used a set of Fr, and Archimedes number Ar (. Other and these widely used dimensionless groups can all be derived by manipulating and/or combining two or more of groups adopted in this work, e.g., and the capillary number, which is the ratio of viscous to surface tension forces. White and Beardmore (1962) described a wide spectrum of experimental results on Taylor bubbles drifting through motionless liquids in vertical tubes. The authors presented a cross plot (Fig. 6 in their paper) showing the regions in which various retarding forces may be neglected. The density of air was neglected and instead of Fr, is plotted to have a reasonable spread of data (White and Beardmore, 1962). In the region where surface tension dominates the bubble does not move at all where the hydrostatic forces are completely balanced by surface tension forces. This occurs at (Hattori, 1935, Bretherton, 1961). For inertia-controlled region when viscosity and surface tension can be neglected (, ), the bubble rise velocity is given solely in terms of Fr (Nicklin et al., 1962, White and Beardmore, 1962, Zukoski, 1966, Mao and Duckler, 1990). In the centre of the graph, the relative magnitude of all retarding forces, namely the viscous, inertial, and interfacial forces are significant. In this region the relationship between Fr, Eo, and M for vertical tubes has been presented by White and Beardmore (1962) as a graphical map which plots lines of constant M on Fr–Eo axis. Similar maps have been produced for non-vertical tubes by Wallis (1969) and Weber et al. (1986).
The description of the motion of single Taylor bubbles in flowing liquids dates back to the pioneering work of Nicklin's et al. (1962) who placed the corner stone of slug flow modelling by recognising the fact that the bubble velocity is a superimposition of two components:The second term represents the drift due to buoyancy (the bubble velocity in a stagnant liquid) and the first term refers to the transport by the mean flow, ). is a dimensionless coefficient that depends on the velocity profile ahead of the bubble, and can be seen as the ratio of the maximum to the mean velocity in the profile. Hence for turbulent flows, while for laminar pipe flow, (Nicklin et al., 1962, Collins et al., 1978, Grace and Clift, 1979, Bendiksen, 1985, Polonsky et al., 1999). The velocity field around the bubble has been measured by many researchers primarily because of its importance regarding the interaction between two Taylor bubbles. Several techniques have been adopted, e.g.: the photochromic dye activation (PDA) method (DeJesus et al., 1995, Kawaji et al., 1997; Ahmad et al., 1998), the particle image velocimetry (PIV) technique (Polonsky et al., 1999; Nogueria et al., 2000) and the laser doppler velocimetry (LDV) technique (Kvernvold et al., 1984). Measurement of wall shear stress in stagnant (Mao and Duckler, 1989) and flowing (Nakoryakov et al., 1989) water showed clearly the reversal of the flow in the liquid film. Experimental work on the propagation of Taylor bubbles in downward flow has been limited to a few studies on ascending bubbles in vertically downward flow (Martin, 1976; Polonsky et al., 1999) and one single study (to our knowledge) on inclined downward flow by Bendiksen (1984). Theoretical treatments of the problem for predicting the rise bubble velocity have been successful for the special case of vertical tubes. This is, however, restricted to the case where both viscous and surface tension effects are negligible (Dumitrescu, 1943; Davies and Taylor, 1950), or the case where viscous effects dominate (Goldsmith and Mason, 1962). Later, Bendiksen (1985) extended the theoretical approach of Dumitrescu (1943) and Davies and Taylor (1950) to account for the surface tension. Mao and Duckler, 1990, Mao and Duckler, 1991 and Clarke and Issa (1997) performed numerical simulations in order to calculate the velocity of the bubble and the velocity field in the liquid film. Propagation of Taylor bubbles in non-vertical tubes has been less well studied experimentally, and even less well theoretically, for the obvious reason of asymmetry.
Thus, studying the motion of a single Taylor bubble in stagnant and in moving liquid is essential in order to understand the intrinsically complicated nature of slug flow. The rise of a single bubble in both stagnant and flowing liquid inside vertical tubes has been extensively studied by numerous researchers, both experimentally and theoretically. All the published numerical methods to model slug flow are restricted only to vertical tubes and they assume either the shape of the bubble or a functional form for the shape (Dumitrescu, 1943, Davies and Taylor, 1950). These assumptions constrain the nature of the solution, while the approach adopted here (the volume of fluid method, VOF) lays no such a priori foundations. The solution domain in the present model not only includes the field around the bubble, as in the study of Mao and Duckler (1990), but also extends behind the bubble, allowing field information to be obtained in the wake region. The difficulty of obtaining local data is even more complicated in horizontal and inclined slug flow by the fact that the flow is asymmetric, and consequently, very few detailed data have been published in the open literature. Previous models available in literature also fail to provide a fully satisfactory mechanistic model of slug flow in horizontal configuration (Duckler and Hubbard, 1975; Taitel and Barnea, 1990). Thus, it is evident that more insight into slug flow is needed to attain a thorough understanding of the internal structure of the slug flow pattern. In this work, an attempt is made to calculate the shape and rising velocity of a single Taylor bubble in stagnant and in moving liquid in vertical tubes. The velocity field in a slug unit (Taylor bubble + liquid plug) and wall shear stress are also calculated. The intention of future work is to characterise the hydrodynamics of slug flow when the tube is tilted away from the vertical and also to investigate the effect of the leading bubble on the shape and velocity of the trailing one.
Section snippets
CFD model development
The CFD software FLUENT (Release 5.4.8, 1998) is used to simulate the motion of a single Taylor bubble rising in a motionless or flowing liquid through a vertical tube. In FLUENT, the control volume method—sometimes referred to as the finite volume method—is used to discretise the transport equations. The movement of the gas–liquid interface is tracked based on the distribution of , the volume fraction of gas in a computational cell, where in the liquid phase and in the gas phase (
Drift of Taylor bubbles in stagnant liquids
Consider a single Taylor bubble rising in stagnant liquid inside a vertical tube with velocity . The bubble may be made stationary by superimposing a downward velocity to the liquid and to the tube walls. The bubble has a round nose and fills almost the cross sectional area of the tube (Fig. 2). The liquid ahead of the bubble moves around the bubble as a thin liquid film moving downwards in the annular space between the tube wall and the bubble surface. Alongside the bubble, the liquid
Conclusions
In this paper we presented a numerical study to investigate the motion of single Taylor bubbles in vertical tubes. A complete description of the bubble propagation in both stagnant and flowing liquids was obtained. The bubble was found to have a cylindrical body with a spherical nose and a fluctuating tail. Sometimes, small bubbles were seen to be sheered off the tail due to the liquid jet coming down from the annular region that separates the bubble from the tube walls. As the bubble moves up,
Notation
Ar Archimedes number, dimensionless dimensionless coefficient, dimensionless dimensionless coefficient, dimensionless Ca capillary number, dimensionless diameter of the tube, m Eo Eötvös number, dimensionless friction factor, dimensionless liquid holdup in the liquid film, dimensionless Fr Froude number, dimensionless g acceleration due to gravity, h liquid film thickness, m length of liquid plug, m length of Taylor bubble, m m dimensionless coefficient, dimensionless M Morton number,
Acknowledgements
T. Taha is grateful to the Karim Rida Said Foundation for financial support and to Dr. David Kenning for helpful suggestions. This work is partially sponsored by EPSRC (GR/66438).
References (50)
- et al.
Falling film hydrodynamics in slug flow
Chemical Engineering Science
(1998) A unified model for predicting flow pattern transition in the whole range of pipe inclination
International Journal of Multiphase Flow
(1987)An experimental investigation of the motion of the long bubble in inclined tubes
International Journal of Multiphase Flow
(1984)On the motion of long bubbles in vertical tubes
International Journal of Multiphase Flow
(1985)- et al.
A continuum method for modeling surface tension
Journal of Computational Physics
(1992) - et al.
A numerical model of slug flow in vertical tubes
Computers and Fluids
(1997) - et al.
Dependence of slug rise velocity on tube Reynolds number in vertical gas–liquid flow
Chemical Engineering Science
(1979) - et al.
Volume of fluid (VOF) method for the dynamics of free boundaries
Journal of Computational Physics
(1981) - et al.
Investigation of flow structures in vertical slug flow
Nuclear Engineering and Design
(1997) - et al.
Velocity distribution in horizontal slug flow
International Journal of Multiphase Flow
(1984)
A flow pattern map for gas–liquid flow in horizontal pipes
International Journal of Multiphase Flow
The motion of Taylor bubbles in vertical tubes—I. A numerical simulation for the shape and the rise velocity of Taylor bubbles in stagnant and flowing liquids
Journal of Computational Physics
The motion of Taylor bubbles in vertical tubes: II. Experimental data and simulations for laminar and turbulent flow
Chemical Engineering Science
Coalescence of two gas slugs in rising in a vertical column of liquid
Chemical Engineering Science
The relation between the Taylor bubble motion and the velocity field ahead of it
International Journal of Multiphase Flow
Stability of severe slugging
International Journal of Multiphase Flow
Two-phase slug flow
Advances in Heat Transfer
On interaction between leading and trailing elongated bubbles in a vertical pipe flow
International Journal of Multiphase Flow
Experimental investigation of the velocity field induced by a Taylor bubble rising in stagnant water
International Journal of Multiphase Flow
Velocities of extended bubbles in inclined tubes
Chemical Engineering Science
The velocity of rise of single cylindrical air bubbles through liquids contained in vertical tubes
Chemical Engineering Science
The motion of long bubbles in tubes
Journal of Fluid Mechanics
An experimental study of the wake of gas slugs rising in liquids
Journal of Fluid Mechanics
The motion of a large gas bubble rising through liquid flowing in a tube
Journal of Fluid Mechanics
The mechanics of large bubble rising through extended liquids and through liquids in tubes
Proceedings of the Royal Society of London, Series A
Cited by (214)
Critical review of vertical gas-liquid slug flow: An insight to better understand flow hydrodynamics' effect on heat and mass transfer characteristics
2024, International Journal of Heat and Mass TransferAn inclusive review on structural enhancement techniques in forced condensation inside tubes
2024, International Journal of Heat and Fluid FlowFlow pattern dependent model for airlift pumps performance: Analytical simulation and experimental verification
2024, Chemical Engineering Research and DesignInvestigation of Taylor bubble dynamics in annular conduits with counter-current flow
2024, International Journal of Multiphase FlowDynamics of Taylor bubble interface in vertical turbulent counter-current flow
2023, International Journal of Multiphase Flow