Co-current flow effects on a rising Taylor bubble

Chemical Engineering Research and Design, 2005

A n experimental study is presented concerning the transition in the velocity of individual Taylor bubbles in vertical co-current liquid flow. Velocities of individual Taylor bubbles rising in co-current liquids (kinematic viscosities from 10 26 to 5.7 Â 10 26 m 2 s 21 ) in acrylic columns of 22 mm, 32 mm and 52 mm internal diameter were measured for a wide range of Reynolds number of the flowing liquid using two nonintrusive experimental techniques. The measuring section was located at 6.0 m from the gas injection. The operating conditions used correspond to inertial controlled regime. The data showed an unexpected feature of the bubble motion: the velocity coefficient C changes even when the flow regime in the liquid ahead the bubble is still laminar, i.e., the transition in the bubble velocity starts at liquid Reynolds numbers much lower than 2100. Additional experiments, employing PIV measurements, showed a developed laminar liquid flow ahead the bubble nose. Based on a dimensional analysis, the most important dimensionless numbers for the phenomena were identified and, after processing all data, an empirical correlation was established to predict the velocity coefficient C for a large range of operation conditions. This information is very important for vertical two-phase slug flow modelling.

This study is an experimental investigation of the interactions between the bubbles, the coherent motion and the viscous drag in a Taylor Couette flow, for the outer cylinder at rest. The cylinder radius ratio η is 0.9. Bubbles are injected through a needle at the bottom of the apparatus inside the gap. Different bubble sizes are investigated (ratio between the bubble size and the gap width 0.05 and 0.12) for very small void fraction (α≤0.012). Different flow regimes are studied corresponding to Reynolds number Re based on the gap width and the velocity of the inner cylinder ranging from 400 to 20000. For these Re values, Taylor vortices are persistent leading to an axial periodicity of the flow. PIV measurements of the liquid flow features, bubble tracking in a meridian plane and viscous torque of the inner cylinder measurements are performed. This study provides a first evidence of the link between the bubble localisation, the Taylor vortices and viscous torque modifications. Bubbles are attracted towards the inner cylinder, due to the rotation of the cylinder. For small buoyancy effect, bubbles are trapped and induce a decrease in the outflow intensity, thus leading to an increase of the viscous torque. When buoyancy induced bubble motion, by comparison to the coherent motion of the liquid is increased, a decrease in the viscous torque is suspected.

Experiments in Fluids, 2001

An experimental study on the interaction between Taylor bubbles rising through a co-current¯owing liquid in a vertical tube with 32 mm of internal diameter is reported. The¯ow pattern in the bubble's wake was turbulent and the¯ow regime in the liquid slug was either turbulent or laminar. When the¯ow regime in the liquid slug is turbulent (i) the minimum distance between bubbles above which there is no interaction is 5D-6D; (ii) the bubble's rising velocity is in excellent agreement with the Nicklin relation; (iii) the experimental values of the bubble length compare well with theoretical predictions (Barnea 1990); (iv) the distance between consecutive bubbles varied from 13D to 16D and is insensitive to the liquid Reynolds number. When the¯ow regime in the liquid slug is laminar (i) the wake length is about 5D-6D; (ii) the minimum distance between bubbles above which there is no interaction is higher than 25D; (iii) the bubble's rising velocity is signi®cantly smaller than theoretical predictions. These results were explained in the light of the ®ndings of Pinto et al. (1998) on coalescence of two Taylor bubbles rising through a co-current liquid.

Physics of Fluids, 2013

International Journal of Multiphase Flow, 2007

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.

Physics of Fluids, 2015

An experimental investigation of the interactions between bubbles, coherent motion, and viscous drag in a Taylor-Couette flow with the outer cylinder at rest is presented. The cylinder radii ratio η is 0.91. Bubbles are injected inside the gap through a needle at the bottom of the apparatus. Different bubbles sizes are investigated (ratio between the bubble diameter and the gap width ranges from 0.05 to 0.125) for very small void fraction (α ≤ 0.23%). Different flow regimes are studied corresponding to Reynolds number Re based on the gap width and velocity of the inner cylinder, ranging from 6 × 10 2 to 2 × 10 4 . Regarding these Re values, Taylor vortices are persistent leading to an axial periodicity of the flow. A detailed characterization of the vortices is performed for the single-phase flow. The experiment also develops bubbles tracking in a meridian plane and viscous torque of the inner cylinder measurements. The findings of this study show evidence of the link between bubbles localisation, Taylor vortices, and viscous torque modifications. We also highlight two regimes of viscous torque modification and various types of bubbles arrangements, depending on their size and on the Reynolds number. Bubbles can have a sliding and wavering motion near the inner cylinder and be either captured by the Taylor vortices or by the outflow areas near the inner cylinder. For small buoyancy effect, bubbles are trapped, leading to an increase of the viscous torque. When buoyancy induced bubbles motion is increased by comparison to the coherent motion of the liquid, a decrease in the viscous torque is rather observed. The type of bubble arrangement is parameterized by the two dimensionless parameters C and H introduced by Climent et al.

Ocean Engineering, 2018

Slug flow is one of the main flow regimes encountered in multiphase flow systems especially in oil and gas production systems. In the present study, the rise of single Taylor bubble through vertical stagnant Newtonian liquid is investigated by performing complete dimensionless treatment followed by an order of magnitude analysis of the terms of equations of motion. Based on this analysis, it is concluded that Froude, Eötvös and Reynolds numbers are the sole physical parameters influencing the dimensionless slug flow equations. Using the guidelines of the order of magnitude analysis, computational fluid dynamics simulation is carried out to investigate the dynamics of Taylor bubbles in vertical pipe using the volume-of-fluid (VOF) method. Good agreement with previous experimental data and models available in the literature is established confirming that the density ratio, viscosity ratio and the initial ratio of bubble size to pipe diameter (/) have minimal effect on the main hydrodynamic features of slug flow. Based on the developed results, correlations for the terminal velocity of the Taylor bubble and the dimensionless wall shear stress are proposed showing the significance of these main dimensionless parameters and support other important theoretical and experimental work available in the literature.

Chemical Engineering Science, 2006

The flow in the wake and near-wake regions of individual Taylor bubbles rising through stagnant and co-current vertical columns of Newtonian liquids was studied, employing simultaneously particle image velocimetry (PIV) and pulsed shadowgraphy techniques (PST). Experiments were made with water and aqueous glycerol solutions covering a wide range of viscosities (1 × 10 −3 Pa s < < 1.5 Pa s), in an acrylic column of 32 mm ID.

International Journal of Multiphase Flow, 2021

When a long gas bubble travels in a horizontal liquid-filled channel of circular crosssection, a liquid film is formed between the bubble and the channel wall. At low Reynolds and Bond numbers, inertial and buoyancy effects are negligible, and the liquid film thickness is a function of the capillary number only. However, as the tube diameter is increased to the millimetre scale, both buoyancy and inertial forces may become significant. We present the results of a systematic analysis of the bubble shape, inclination, and liquid film thickness for a wide range of capillary, Bond, and Reynolds numbers, namely 0.024 ≤ Ca l ≤ 0.051, 0.11 ≤ Bo ≤ 3.5, and 1 ≤ Re l ≤ 750. Three-dimensional numerical simulations of the flow are performed by employing the Volume-Of-Fluid method implemented in OpenFOAM. In agreement with previous studies, we observe that buoyancy lifts the bubble above the channel axis, making the top liquid film thinner, and thickening the bottom film. As the Bond number approaches unity, the cross-sectional shape of the bubble deviates significantly from a circular shape, due to flattening of the bottom meniscus. The simulations demonstrate the existence of a cross-stream film flow that drains liquid out of the top film and drives it towards the bottom film region. This drainage flow causes inclination of the bubble, with a larger inclination angle along the bottom plane of the bubble than the top. As buoyancy becomes even more significant, draining flows become less effective and the bubble inclination reduces. A theoretical model for the liquid film thickness and bubble speed is proposed embedding dependencies on both capillary and Bond numbers, which shows good agreement with the reported numerical results. Inertial forces tend to shrink the bubble cross-section and further lift the bubble above the channel centreline, so that the bottom film thickness increases significantly with the Reynolds number, whereas the top film thickness is less sensitive to it.

Asia-Pacific Journal of Chemical Engineering, 2017

The coalescence and rising behavior of co-axial and lateral bubbles in viscous fluid is reported for nonambient conditions in line with actual operating conditions of industrial bubble column. The volume of fluid with the continuum surface force (VOF-CSF) model embedded in FLUENT CFD package is used. Three main case studies, namely under ambient and aqueous condition (μ* = σ* = 1), under reduced viscosity (μ* = 0.1) and surface tension (σ* = 0.1) are investigated and compared using two co-axial bubble and three lateral bubble rise configuration models. The latter two cases represent operating conditions at high temperature or pressure. Details of co-axial bubble (or a pair of bubbles rising in a vertical line) and three lateral bubble rise and coalescence characteristics are presented.