Development of co-current air–water flow in a vertical pipe

Measurements of the cross-sectional distribution of the gas fraction and bubble size distributions were conducted in a vertical pipe with an inner diameter of 51.2 mm and a length of about 3 m for air/water bubbly and slug flow regimes. The use of a wire-mesh sensor obtained a high resolution of the gas fraction data in space as well as in time. From this data, time averaged values for the two-dimensional gas fraction profiles were decomposed into a large number of bubble size classes. This allowed the extraction of the radial gas fraction profiles for a given range of bubble sizes as well as data for local bubble size distributions. The structure of the flow can be characterized by such data. The measurements were performed for up to 10 different inlet lengths and for about 100 combinations of gas and liquid volume flow rates. The data is very useful for the development and validation of meso-scale models to account for the forces acting on a bubble in a shear liquid flow and models for bubble coalescence and break-up. Such models are necessary for the validation of CFD codes for the simulation of bubbly flows.     

Local gas and liquid phase velocity measurement in a miniature stirred vessel using PIV combined with a new image processing algorithm

A method which combines standard two-dimensional particle image velocimetry (PIV) with a new image processing algorithm has been developed to measure the average local gas bubble velocities, as well as the local velocities of the liquid phase, within small stirred vessel reactors. The technique was applied to measurements in a gas–liquid high throughput experimentation (HTE) vessel of 45 mm diameter, but it is equally suited to measurements in larger scale reactors. For the measurement of liquid velocities, 3 μm latex seeding particles were used. For gas velocity measurements, a separate experiment was conducted which involved doping the liquid phase with fluorescent Rhodamine dye to allow the gas–liquid interfaces to be identified. The analysis of raw PIV images enabled the detection of bubbles within the laser plane, their differentiation from obscuring bubbles in front of the laser plane, and their use in lieu of tracer particles for gas velocity analysis using cross-correlation methods. The accuracy of the technique was verified by measuring the velocity of a bubble rising in a vertical glass column. The new method enabled detailed velocity fields of both phases to be obtained in an air–water system. The overall flow patterns obtained showed a good qualitative agreement with previous work in large scale vessels. The downward liquid velocities above the impeller were greatly reduced by the addition of the gas, and significant differences between the flow patterns of the two-phases were observed.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

Pressure-oscillation defoaming for viscoelastic fluid

The enhancement of bubble rising velocity was experimentally investigated by mechanically applying an oscillating pressure to a single small air bubble (e.g., 1 mm³) in a viscoelastic fluid. For shear-thinning fluids, the cyclic change in bubble diameter induced by the oscillating pressure generates a continuous strong local flow near the bubble surface. Consequently, the apparent liquid viscosity is reduced and the bubble rising velocity increases by 400 times or more compared to the case without oscillating pressure. However, for a Newtonian fluid, almost no effect was observed with oscillating pressure. Time-series data of the longitudinal and horizontal bubble diameters were obtained experimentally using a stroboscope and a video system, and these data were used to estimate the local shear rate and the local shear viscosity. The increase in bubble rising velocity estimated from the shear viscosity behavior agreed well with the experimental data. Additionally, a periodic change in the bubble shape from a sphere at the maximum bubble size to a cusped shape at the minimum bubble size was observed under strong oscillatory pressure.     

Bubble dynamics and interactions with a pair of micro pillars in tandem

This study investigates flow patterns and bubble dynamics of two-phase flow around two 100 μm diameter circular pillars in tandem, which were entrenched inside a horizontal micro channel. Bubble velocity, trajectory, size, and void fraction were measured using a high speed camera and analyzed using a particle tracking velocimetry method. A range of gas and liquid superficial velocities were tested, resulting in different bubbly flow patterns, which were consistent with previous studies. These flow patterns were altered as they interacted with the pillars. Depending on the relative transverse location of bubbles to the pillars, and through bubble–bubble interaction, the flow sometimes returned to its original state. It was also determined that the pillars altered both the bubble trajectory and void fraction, especially in the pillars region.     

On the liquid turbulence energy spectra in two-phase bubbly flow in a large diameter vertical pipe

The liquid turbulence structure of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. A dual optical probe was used to measure the bubble characteristics, while the liquid turbulence was measured using hot-film anemometry. Experiments were performed at two liquid superficial velocities of 0.2 and 0.68 m/s for gas superficial velocities in the range of 0–0.18 m/s, corresponding to an area averaged void fraction up to 13.6%. In general, there is an increase in the liquid turbulence energy when the bubbles are introduced into the liquid flow. The increase in the energy mainly occurs over a range of length scales that are on the order of the bubble diameter. A suppression of the turbulence was observed close to the wall at very low void fraction flows. Initially, the suppression occurs in the low wave number range and then extends to higher wave numbers as the suppression is increased.     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Investigation of bubble breakup and coalescence in a packed-bed reactor – Part 1: A comparative study of bubble breakup and coalescence models

In a packed-bed reactor a comparative study of bubble breakup and coalescence models has been investigated to study bubble size distributions as a function of the axial location. The bubble size distributions are obtained by solving population balance equations that describe gas–liquid interactions. Each combination of bubble breakup and coalescence models is examined under two inlet flow conditions: (1) predominant bubble breakup flow and (2) predominant bubble coalescence flow. The resulting bubble size distributions, breakup and coalescence rates estimated by individual models, are qualitatively compared to each other. The change of bubble size distributions along the axial direction is also described with medians. The medians resulting from CFD analyses are compared against the experimental data. Since the predictions estimated by CFD analyses with the existing bubble breakup and coalescence models do not agree with the experimental data, a new bubble breakup and coalescence model that takes account of the geometry effects is required to describe gas–liquid interactions in a packed-bed reactor.     

Void fraction,bubble velocity and bubble size in two-phase flow

Experiments were performed in atmosphereic vertical air-water flows, for void fractions between 0.25 and 0.75 (cross-sectional averages) and superficial liquid velocities of 1.3, 1.7 and 2.1 m/s. Local values of void fraction and bubble velocity as well as the bubble diameter were measured by means of a resistivity probe technique. Reliable values were obtained for the local void fraction over the entire range $\text{0 ≤ α ≤ 1}$. The void fraction profiles appeared to have a local maximum at the pipe center, local maxima close to the wall were obviously absent. The resistivity probes are shown to measure the velocity of the interface between the conducting and nonconducting phases, which equals the gas velocity only for low void fractions. The measured data for void fraction and bubble velocity were correlated by means of power law distribution functions, with exponents given by a function of the cross-sectionally averaged void fraction. The Sauter mean diameters for the bubble size spectra found, agree reasonably well with diameters predicted by a theoretical model based on the energy dissipation in the flow.     

Modeling of bubble behaviors and size distribution in a slab continuous casting mold

Population balance equations combined with Eulerian–Eulerian two-phase model are employed to predict the polydispersed bubbly flow inside the slab continuous-casting mold. The class method, realized by the MUltiple-SIze- Group (MUSIG) model, alongside with suitable bubble breakage and coalescence kernels is adopted. A two-way momentum transfer mechanism model combines the bubble induced turbulence model and various interfacial forces including drag, lift, virtual mass, wall lubrication, and turbulent dispersion are incorporated in the model. A 1/4th scaled water model of the slab continuous-casting mold was built to measure and investigate the bubble behavior and size distribution. A high speed video system was used to visualize the bubble behavior, and a digital image processing technique was used to measure the mean bubble diameter along the width of the mold. Predictions by previous mono-size model and MUSIG model are compared and validated against experimental data obtained from the water model. Effects of the water flow rate and gas flow rate on the mean bubble size were also investigated. Close agreements by MUSIG model were achieved for the gas volume fraction, liquid flow pattern, bubble breakage and coalescence, and local bubble Sauter mean diameter against observations and measurements of water model experiments.     

Bubble formation in co-fed gas–liquid flows in a rotor-stator spinning disc reactor

The gas–liquid flow in a rotor-stator spinning disc reactor, with co-feeding of gas and liquid, is studied for high gas volumetric throughflow rates and high gas/liquid volumetric flow ratios. High speed imaging and spectral analysis of pressure drop signals are employed to analyse the flow. Two mechanisms of bubble formation are observed, one due to gas overpressure leading to large irregular bubbles, and one due to liquid turbulent vortices leading to small, well-defined bubbles. The two mechanisms lead to three distinct gas dispersion regimes, distinguished by their characteristic oscillations in pressure drop. At low rotational Reynolds numbers (Re_ω < 0.4 · 10⁶), in the gas spillover regime, the gas is dispersed as large bubbles only. Above this critical Re_ω, small bubbles are sheared off as well, thus forming a heterogeneous dispersion. At sufficiently high Re_ω, depending on the gas flow rate, the gas is homogeneously dispersed as small bubbles. The maximum gas flow that can be dispersed as small bubbles is linearly proportional to the local energy dissipation rate. The understanding of the bubble formation mechanisms and pressure signature allows prediction and detection of the prevailing hydrodynamic regime in scaled up spinning disc reactors and for different reaction fluids.     

环形喷管喷口气泡演化的实验研究

水下气泡的发展演化及气泡动力学行为是气液两相动力学的基础理论与水下射流应用的重要基础. 环形喷管/喷口形成的气泡及气体射流具有其不同于圆孔实心射流的特殊表现与规律机制,随着同心筒破水发射等特殊应用的出现,环形喷口气体射流/泡流的基础现象观测和机制分析成为迫切的需求. 基于环形喷管的设计和水下射流条件的分析,设计建立了一套环形喷管水箱实验系统,对水下环形喷管喷口气泡发展演化过程进行了初步的实验研究. 为观测研究气体通过环形喷管气泡生长发展过程,在较低压力、较低流速下,采用高速摄影仪记录气泡生长及发展演化过程. 结合对气泡发展演化过程的图像处理与分析,研究分析了环形喷口气泡形成区制、气泡生长过程形态发展特点、以及气泡形成时间及气泡体积变化特点. 研究表明:在本实验气体流量范围内(50.8~237.3 dm3/min),环形喷口气泡发展演化过程呈现较为明显的三周期区制,前泡尾流影响是环形气泡呈三周期区制的主要原因;不同周期内的气泡形成时间具有较稳定规律,并受到流量影响;气泡生长过程中有较为明显的下沉、回升特征;气泡表面张力、液体惯性与流动的共同作用,造成了典型的气泡顶部坍塌现象.      

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work, co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 mm in diameter. Each tube includes two 10 m long pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly flow, stratified flow, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high-speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood–Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/non-Newtonian liquid.     

Experimental investigation of three-phase flow in a vertical pipe: Local characteristics of the gas phase for gas-lift conditions

Laboratory experiments have been performed on the flow of oil, water and air through a vertical pipe in order to study the gas-lift technique for oil–water flows. Special attention was paid to the phase inversion phenomenon, by which the continuous phase switches to the dispersed phase and vice versa. By using different types of gas injectors the influence of the bubble size of the injected air on the efficiency of the gas-lift technique (in particular at the point of phase inversion) was studied. Also the gas and liquid mixture velocities were varied. The air bubbles were detected by means of optical fibre probes. Local measurements of the time-averaged gas volume fraction, bubble size and bubble velocity were carried out, as well as pressure measurements.     

On the numerical study of bubbly flow created by ventilated cavity in vertical pipe

The dispersion of bubbles into a down-liquid flow in a vertical pipe is investigated. At low flow rates, the intended design of a swarm of discrete bubbles is achieved. At high flow rates, a ventilated cavity is nonetheless formed, which is attached close to the gas sparger. Behind this ventilated cavity, three different flow regimes characterize the complex bubbly flow field downstream of the down-liquid flow: vortex region with high void fraction, transitional region and pipe flow region. In this study, a numerical model that solved the entire development of the gas–liquid flow including the extended single-phase liquid region upstream to the wall-jet and recirculating-vortex zones in order to allow a more realistic determination of the boundary conditions of the down-liquid flow was adopted. Coupling with the Eulerian–Eulerian two-fluid model to solve the respective gas and liquid phases, a population balance model was also applied to predict the bubble size distribution in the wake right below the cavity base as well as further downstream in the transitional and fully-developed pipe flow regions. The numerical model was evaluated by comparing the numerical results against the data derived from theoretical, numerical and experimental approaches. Prediction of the Sauter mean bubble diameter distributions by the population balance approach at different axial locations confirmed the dominance of breakage due to the high turbulent intensity below the ventilated cavity which led to the generation of small gas bubbles at high void fraction. Further downstream, the coalescence effect dominated leading to merging of the small bubbles to form bigger bubbles.     

Experimental investigation of single helium bubbles rising in FLiNaK molten salt

An experimental investigation of single helium bubbles rising in a stagnant molten fluoride salt mixture was conducted in a 25.4 mm quartz tube at 600 °C and atmospheric pressure. The fluoride salt chosen for this research was the eutectic mixture of LiF-NaF-KF (46.5–11.5–42%), also known as FLiNaK. The images obtained using a high-speed camera were processed to estimate the bubble size, rising velocity, trajectory and shape. The shape of the bubble showed typical characteristics of wobbling bubble corresponding to published classifications. The trajectory of the rising bubble showed complex path oscillations and combinations of rectilinear, zig-zag and helical motion. The projected area-equivalent diameter and experimental terminal velocity varied from 5.26 to 6.23 mm and 232.35–259.57 mm/s, respectively. Several important dimensionless numbers were calculated and reported based on the experimental results, which are closely related to the phenomena involving a gas bubble rising in liquid. The measured terminal velocity was compared with the predictions of existing correlations that include the dimensionless numbers and showed reasonable agreement. The Particle Image Velocimetry (PIV) technique was applied to complement the experimental data set with high-fidelity measurements of the velocity field of the liquid molten salt surrounding the bubble. The paper presents a unique set of PIV and two-phase experimental data capturing the behavior of rising helium bubbles in molten FLiNaK and its surrounding flow field.     

Eulerian–Lagrangian 3-D simulations of unsteady two-phase gas–liquid flow in a rectangular column by considering bubble interactions

This work discusses the development of a three-dimensional Eulerian–Lagrangian CFD model for a gas–liquid flow in a rectangular column. The model resolves the time-dependent, three-dimensional motion of small gas bubbles in a liquid to simulate the dynamic characteristics of the oscillating bubble plume. Our model incorporates drag, gravity, buoyancy, lift, pressure gradient and virtual mass forces acting on a bubble rising in a liquid, and accounts for two-way momentum coupling between the phases. We use MUSIG model that provides a framework in which the population balance method together with the break up and coalescence models can be incorporated into three-dimensional CFD calculations. We use turbulent flow to describe liquid flow field. The standard κ–ε of turbulence is selected for calculating the properties of turbulent flow. The effect of aspect ratio of the column on the flow pattern, liquid velocity and gas hold-up profiles is discussed.     

Eulerian–Lagrangian simulations of unsteady gas–liquid flows in bubble columns

We studied the dynamics of gas–liquid flows in a rectangular bubble column using Eulerian–Lagrangian simulations. Three-dimensional, unsteady simulations were performed to simulate the dynamic characteristics of the oscillating bubble plume. The effect of superficial gas velocity and aerated liquid height-to-column width (H/W) ratio on the dynamic and time-averaged flow properties was studied and the simulated results were validated using wall pressure and voidage fluctuation measurements. The effect of lift force and numerical diffusion on the dynamic and time-averaged properties is discussed in detail. Further, the results obtained using the Eulerian–Lagrangian simulations were compared with the Eulerian–Eulerian simulations. The bubble scale information, which is otherwise lost in the Eulerian–Eulerian simulations, was validated using the voidage fluctuation measurements. Such experimentally validated Eulerian–Lagrangian models will be useful for the simulation of mass transfer and reactions in bubble columns.     

On the role of the lateral lift force in poly-dispersed bubbly flows

An extensive experimental database comprising air–water as well as steam-water upwards vertical pipe flows for a pressure up to 6.5 MPa was used to investigate the effect of the lateral lift force on turbulent poly-dispersed flows with medium or high gas volume fraction. It was clearly shown that the lift force plays an important role also in such flows. Several effects such as bubble coalescence and breakup as well as fast rising large bubbles which push small bubbles towards the pipe wall superpose the effect of the lift force but can be separated from this effect. The critical bubble diameter, at which the lift force changes its sign, predicted by using Tomiyama’s correlation agrees well with experimental data obtained for turbulent air–water and steam-water flows with medium and high void fraction and a broad spectrum of bubbles sizes. The values for this critical bubble diameter are confirmed by the experimental data within the frame of the uncertainty of the data. Consequences of the action of the lateral lift force on flow structures in different flow situations are discussed. From the investigations it can be concluded that the lift force including the bubble size dependent change of its sign should be considered in a proper numerical 2D or 3D-simulation on flows in which bubbles in the range of several millimeters are present.     

Elongated bubble shape in inclined air-water slug flow

The shape of elongated bubbles in upward inclined air-water slug flow was studied experimentally by quantitative measurements of the cross sectional distribution of the phases within the pipe, using a wire mesh sensor. Ensemble-averaged shapes of elongated bubbles were determined for a wide range of gas and liquid flow rates, as well as for different pipe inclination angles. The elongated bubble nose can be characterized by an annular domain where liquid is present above the gas. The effect of gas and liquid flow rates, as well as of the pipe inclination angle on the bubble shape (front and tail) is studied. A simplified theoretical model is proposed to determine the bubble front shape. The model predictions compare favorably with the experimental results.