Hydrodynamic stability of a bubble column with a bottom-mounted point air source

Bubble column is widely used in both industrial and environmental applications. In this study, we examine the flow dynamics and stability of a bubble column driven by a point air source centrally mounted at the bottom using Phase Doppler anemometry (PDA). The model cylindrical bubble column had an inner diameter of 152 mm and was filled with the liquid to about 1 m height, above the point air source, which was made of a 30-mm diameter perforated air stone. The bubble diameters were within the range of 400–1300 μm. A customized setup was developed for accurate PDA measurements of the two phases, and detailed turbulent characteristics of the liquid phase velocity, bubble diameter, bubble velocity and the slip velocity were collected throughout the column. The comprehensiveness of the data set enabled a close examination of the hydrodynamic stability inside the column. Measurements were taken at three different air rates, namely 0.13, 0.25 and 0.38 L/min (corresponding to average gas volume fractions of 0.0065, 0.0138 and 0.0197, respectively). The results illustrated a large-scale coherent liquid circulation pattern inside the column. The circulation pattern in the upper column was relatively steady, while the pattern in the lower column was strongly unsteady with the probability density functions (pdf) for both the liquid and bubble velocities showing distinct twin peaks. An analysis based on the determination of the bubble drag forces and transversal lift forces is performed by decomposing the twin-peaked pdfs into two separated Gaussian distributions, one for the upward flow due to the bubble rises and the other for the downward flow due to circulation. Through the decomposition, a stability criterion can then be established by choosing the local bubble size as the representative length scale for the turbulent eddies inside the column. The analysis with the criterion illustrates why a steady circulation pattern was achieved in the upper column, and at the same time shows that the instability at the bottom column was induced by the low frequency meandering of the bubble swarm.     

CFD simulation and experimental measurement of gas holdup and liquid interstitial velocity in internal loop airlift reactor

This paper documents experiments and CFD simulations of the hydrodynamics of our two-phase (water, air) laboratory internal loop airlift reactor (40 l). The experiments and simulations were aimed at obtaining global flow characteristics (gas holdup and liquid interstitial velocity in the riser and in the downcomer) in our particular airlift configurations. The experiments and simulations were done for three different riser tubes with variable length and diameter. Gas (air) superficial velocities in riser were in range from 1 to 7.5 cm/s. Up to three circulation regimes were experimentally observed (no bubbles in downcomer, bubbles in downcomer but not circulating, and finally the circulating regime). The primary goal was to test our CFD simulation setup using only standard closures for interphase forces and turbulence, and assuming constant bubble size is able to capture global characteristics of the flow for our experimental airlift configurations for the three circulation regimes, and if the simulation setup could be later used for obtaining the global characteristic for modified geometries of our original airlift design or for different fluids. The CFD simulations were done in commercial code Fluent 6.3 using algebraic slip mixture multiphase model. The secondary goal was to test the sensitivity of the simulation results to different closures for the drag coefficient and the resulting bubble slip velocity and also for the turbulence. In addition to the simulations done in Fluent, simulation results using different code (CFX 12.1) and different model (full Euler–Euler) are also presented in this paper. The experimental measurements of liquid interstitial velocity in the riser and in the downcomer were done by evaluating the response to the injection of a sulphuric acid solution measured with pH probes. The gas holdup in the riser and downcomer was measured with the U-tube manometer. The results showed that the simulation setup works quite well when there are no bubbles present in the downcomer, and that the sensitivity to the drag closure is rather low in this case. The agreement was getting worse with the increase of gas holdup in the downcomer. The use of different multiphase model in the different code (CFX) gave almost the same results as the Fluent simulations.     

Use of models for lift, wall and turbulent dispersion forces acting on bubbles for poly-disperse flows

Closure laws are needed for the qualification of CFD codes for two-phase flows. In case of bubbly and slug flow, forces acting on the bubbles usually model the momentum transfer between the phases. Several models for such forces can be found in Literature. They show, that these forces depend on the liquid flow field as well as on the size and the shape of the bubbles. A validation of consistent sets of bubble force models for poly-disperse flows is given, basing on a detailed experimental database for vertical pipe flows, which contains data on the radial distribution of bubbles of different size as well as local bubble size distributions. A one-dimensional (1D) solver provides velocity profiles and bubble distributions in radial direction. It considers a large number of bubble size classes and is used for the comparison with the experiments. The simplified model was checked against the results of full 3D simulations done by the commercial code CFX-5.7 for simplified monodisperse cases. The effects of the number of bubbles classes as well as the effect of the lateral extension of the bubbles were analyzed. For the validation of bubble force models measured bubble size distributions were taken as an input for the calculation. On basis of the assumption of an equilibrium of the lateral bubble forces, radial volume fraction profiles were calculated separately for each bubble class. In the result of the validation of different models for the bubble forces, a set of Tomiyama lift and wall force, deformation force and Favre averaged turbulent dispersion force was found to provide the best agreement with the experimental data. Some discrepancies remain at high liquid superficial velocities.     

鼓泡塔内气液两相湍流实验研究

介绍了研究鼓泡塔气液两相流的实验装置、实验方法。液相用激光多普勒测速技术(LDV)测量,气相用粒子示踪测速技术(PIV)测量。实验表明,轴向液相速度的径向分布呈塔中心峰值、壁面附近倒流形式,且与气相表观速度大小有关,当液相表观速度一定时,随气相表观速度增大而愈加陡峭,返混也剧烈。当表观液速与表观气速之比小于19．6时,返混区总是存在,且返混区大小与高度有关：当表观液遣与表观气速之比大于19．6时,返混消失,含气率分布由塔中心峰值转向壁面峰值。径向液相速度既与气相表现速度有关又与位置高度有关,在塔底部呈现负值,这意味着向塔轴心方向流动。随着塔高增加。流动方向逐渐转变为向塔壁方向,且又有明显的峰值。     

Single-phase flow model development for macroscopic liquid flow evaluation in gas–liquid reactors, by computational fluid dynamics

This paper proposes a single-phase flow model to simulate the flow induced in a liquid by the injection of gas dispersed in the form of a bubble curtain. It aims at predicting macroscopic liquid flow and mixing time. This single-phase flow model is developed as an alternative to two-phase flow models. The model is based on the assumption that the liquid flow is induced by a density imbalance between the bulk zone and the bubble curtain zone. The density in the bulk is set to the water density while the density in the bubble curtain corresponds to the air–water mixture density and is assessed by numerical simulations, thanks to an iterative procedure. Only the knowledge of the injected air flow rate and the bubble liquid relative velocity is required. The single-phase flow model is applied to assess the liquid flow and the mixing in open quarries having a complex geometry. The liquid velocities and the flow structure in the open quarries simulated with the single-phase flow model are in good agreement with those predicted by numerical simulations based on a two-phase flow model.     

Time-dependent gas holdup variation in an air-water bubble column

Time-dependent gas holdup variation in a two-phase bubble column is reported with air and tap water as the working fluids. The results indicate that time-dependent gas holdup is closely related to the water, whose quality is unsteady and changes, not only during the two-phase flow, but also during idle periods. The significance and characteristics of the time-dependent gas holdup variation are influenced by the bubble column operation mode (cocurrent or semi-batch), the sparger orientation, the superficial gas velocity, and the superficial liquid velocity. It is proposed that a volatile substance (VS), which exists in the water in very small concentrations and inhibits bubble coalescence, evaporates during column operation and results in a time-dependent gas holdup. The influence of bubble column operation mode, sparger orientation, superficial gas velocity, and superficial liquid velocity on the time-dependent gas holdup variation are explained based on their effects on bubble size, bubble contacting frequency and mixing intensity. This work reveals that regular tap water may cause significant reproducibility problems in experimental studies of air-water two-phase flows.     

Buoyancy-driven motion of bubbles in square channels

The motion of air bubbles in square capillaries moving under the influence of gravity is studied over a range of Reynolds numbers. The steady shapes and velocities of the bubbles as a function of the bubble size are determined experimentally at moderate Bond and capillary numbers. Bubbles are nearly spherical at lower bubble volumes and become prolate losing their fore and aft symmetry at larger bubble volumes. At higher Weber numbers, a reentrant cavity develops at the rear of bubble. The critical Weber number at which this shape transition occurs lies between 0.89 and 1.38. At small Weber numbers, the terminal velocity of bubbles increases monotonically with bubble volume and eventually reaches a plateau value, which is independent of the bubble size. At higher Weber numbers, a maxima develops in the velocity-volume curve at moderate bubble sizes which grows in magnitude as the Weber number increases. Even at small bubble volumes with nearly spherical shape, the terminal velocity of the bubbles is less than the Hadamard-Rybczynski velocity due to the wall drag. The speed and the maximum bubble width for air bubbles rising in a square channel is higher than that of an air bubble rising in a circular channel with the same hydraulic diameter. The experimental data compares well with predicted trends in the viscous and inertial limits for long bubbles.     

Shape oscillations and path transition of bubbles rising in a model bubble column

We investigate experimentally the occurrence of shape oscillations accompanied by path transition of periodically produced air bubbles rising in water. Within the period of bubble formation, the induced velocity is measured to examine bubble-liquid and bubble-bubble interactions. The flow is produced in a small-scale bubble column with square-shaped cross section. A capillary aerator produces bubbles of size 3.4 mm at a frequency of 5 Hz. Measuring techniques employed are high-speed imaging to capture bubble shape oscillations and path geometry, and laser-Doppler anemometry (LDA) to measure the velocity in the liquid near the rising bubbles. The experimentally obtained bubble shape data are expanded in Legendre polynomials. The results show the occurrence of oscillations by the periodicity of the expansion coefficients in space. Significant shape oscillations accompanied by path transition are observed as the second-mode oscillation frequency converges to the frequency of the initial shape oscillations. The mean velocity field in the water obtained by LDA agrees well with potential theory. An analysis of the decay of the induced flow shows that there is no interaction between the flow fields of two succeeding 3.4 mm bubbles in the rectilinear path when the bubble production frequency is lower than 7.4 Hz.     

Configurations and dynamics of single air/alcohol gas-liquid compound drops in vegetable oil

Gas-liquid compound drops can be used as a low energy liquid-in-liquid disperser, mixer and transporter in place of the energy intensive mechanical agitator. This paper presents an experimental investigation of single air/methanol and air/ethanol gas-liquid compound drops (alcohol-floated bubbles) in sunflower oil at room temperature. Dispersed methanol (ethanol) droplets rise and dissolve slowly in vegetable oil. The air bubbles can float such droplets and increase their rise velocities. The shape of alcohol droplet on the bubble surface was changed to be an alcohol film covering around the bubble and was continuously thrown into the sunflower oil in the form of a toroidal-shaped film, which increases a huge contact area between the alcohol and the oil. A double tube nozzle was used to inject the air bubble (internal tube) and the alcohols (external tube) into the oil. The images showing the configurations of the compound drops during their rise are shown and discussed. The experimental results of the rise velocity and the drag coefficient of the compound drops are presented and compared to those of the single air bubbles and the single alcohol droplets. Rise velocities and drag coefficients predicted by theory of immobile and mobile surface air bubble in sunflower oil are also included for comparison.     

Experimental and numerical investigation of liquid circulation induced by a bubble plume in a baffled tank

Bubble induced liquid circulation is important in applications such as bubble columns and air-lift reactors. In this work, we describe an experimental and numerical investigation of liquid circulation induced by a bubble plume in a tank partitioned by a baffle. The baffle divides the tank into two compartments. Liquid can flow from one compartment to the other through openings at the top and the bottom of the baffle. Gas (air) was injected in the riser section in the form of bubbles at one corner of the tank. The temporal and spatial variation of velocity field in the liquid as a function of the gas flow rate was measured using particle image velocimetry (PIV). At a constant gas flow rate, the liquid flow field is unsteady due to the interaction with the bubbles. The time scales associated with the velocity-time series and the bubble plume thickness variation were calculated. The time averaged-velocity field was used to quantify the variation of the liquid circulation rate with gas flow rate. The turbulence in the liquid was measured in terms of turbulent intensities. These were calculated from the experimental data and were observed to be less than 3 cm/s. A 2-d Euler-Euler two-fluid model with buoyancy and drag as the interaction terms was used to simulate the flow. The parameters chosen for the simulations were selected from literature. It is shown that inclusion of turbulence model such as k-ε is necessary to capture the overall flow behavior. Good agreement was observed between experimentally obtained velocity profiles and the recirculation rates with the simulation results.     

Inviscid bubble formation on porous plates and sieve plates

We experimentally and theoretically investigate the formation of bubbles on wetted and non-wetted sieves and porous plates, which are submerged in an inviscid liquid. It is assumed that the process of bubble formation occurs as if each bubble were alone on the plate. The interaction between orifices is modeled by estimating the gas flow field caused by a growing bubble. Depending on the material below the plate, a growing bubble influences the formation of bubbles at other orifices either in its vicinity or over the whole plate, and bubbling sites are regularly or irregularly distributed over the plate. Values for the volume fluxes at which the bubbling behavior changes are given, e.g., the minimum gas volume flux at which a disperser becomes effective or the volume flux at which growing bubbles densely populate the plates.For an example system of air in water we discuss the four cases of wetted sieve, non-wetted sieve, wetted porous plate and non-wetted porous plate. We present data for a non-wetted porous plate and find fair agreement between theoretical predictions and experiment.     

Application of the double optic probe technique to distorted tumbling bubbles in aqueous or organic liquid

The optic probe technique is widely used to investigate bubble reactors. To derive values of bubble local velocities and bubble local sizes, a specific signal treatment is usually applied under severe assumptions for bubble path and shape. However, in most industrial reactors, bubble motion is chaotic and no common shape can be assumed.In this work, the reliability of the signal treatment associated with the optic probe technique is examined for distorted and tumbling bubbles. A double-tip optic probe is settled in a glass tank and the rise of bubbles is filmed simultaneously. Several trains of bubbles are studied, interactions between bubbles being gradually increased.Referring to image analysis, several ways to derive mean bubble velocities from optic probe data have been compared. Crenels from front tip and rear tip raw signals are associated and individual bubble velocities are derived. Nevertheless, complete velocity distributions are difficult to obtain, as they depend on the choice of the time within which the bubble is searched on the second tip. Using a simpler approach it is shown that the most probable velocity, calculated through the raw signals inter-correlation, is a correct estimation of the average bubble velocity.Concerning bubble size, bubble chord distributions show too high values due to bubble distortion and deviation. A simplified estimation of bubble mean Sauter diameter, using the most reliable measurements only (i.e., local gas hold-up, local mean bubbling frequency, and most probable bubble velocity), was tested for highly distorted bubbles; this method was validated both in water and cyclohexane.     

Double-sensor conductivity probe applied in a flotation system

This study proposed a new approach for measuring bubble size distribution, bubble mean diameter, Sauter mean bubble diameter, and gas holdup using a double-sensor conductivity probe in an air/water two-phase system bubble column. The results for the two-phase system were compared and calibrated using analyses from bubble images taken by a digital camera from the side of the column wall. Good agreement was observed between the two techniques. The same double-sensor conductivity was used in an air/water/solids three-phase system. The conductivity probe captured the change in bubble dynamic behaviour inside the pulp phase; however, the presence of the solids made it more challenging to measure. As a result, the VisioFroth commercial package, using images taken from the top of the froth layer, could be used in conjunction with the double-sensor conductivity probe to show the dynamic evolution of mineralized bubbles from the pulp zone to the froth zone in a flotation process.     

The motion of small air bubbles in stagnant and flowing water

The motion and interaction of small air bubbles rising vertically in stagnant and flowing water has been studied. Singlet, doublet and triplet bubble configurations with bubble generation frequencies of between 0.5 and 3.0 bubbles per second were used to examine the effects of radial and axial bubble interaction. The observed bubble flow behavior and subsequent analysis has allowed an identification of both the steady state and unsteady state (accelerating) regions. The critical separation distance, within which non-steady state prevails, has been determined and the extrapolated limiting value for a Taylor bubble confirms the accepted value. Based on the experimental observations and the analysis, a prediction method to determine the details of bubble flow behavior within a limited range of the variables is proposed.     

Advanced PIV/LIF and shadowgraphy system to visualize flow structure in two-phase bubbly flows

Particle image velocimetry (PIV) is a promising technique to measure dispersed phase size, dispersed phase hold-up and velocity of both the phases. The current work reports measurement of the shape, size, velocity and acceleration of bubbles using shadowgraphy, and liquid velocity measurement obtained using PIV/LIF with fluorescent tracer particles. Measurements were performed in a narrow rectangular column at moderate gas hold-up (～5%) with wide variation of bubble sizes (0.1–15 mm). The liquid velocity field was subjected to 2D discrete wavelet transform (DWT) to visualize the flow structures in the bubbly flow. Further, the slip velocity of individual bubbles was obtained from the DWT filtered liquid velocity field. The results are compared with the slip velocity correlations reported in literature for single bubbles rising in quiescent water. The comparison shows the difference in slip velocity of single bubbles and bubbles rising in swarm. The scale wise decomposition obtained from DWT was also used to quantify the liquid velocity field in terms of wavenumber spectrum. The velocity and acceleration measurements are demonstrated on a single spherical cap bubble rising in quiescent water. The measurements show the potential of the 2D acceleration measurement to facilitate the estimation of unsteady drag on bubbles.     

Study of Air‐Liquid Flow Patterns in Hydrocyclone Enhanced by Air Bubbles

In order to improve the oil‐water separation efficiency of a hydrocyclone, a new process utilizing air bubbles has been developed to enhance separation performance. Using the two‐component phase Doppler particle analyzer (PDPA) technique, the velocities of two phases, air and liquid, and air bubble diameter were measured in a hydrocyclone. The air‐liquid mixing pump can produce 15 to 60 μm‐diameter air bubbles in water. There is an optimum air‐liquid ratio for oil‐water separation of a hydrocyclone enhanced by air bubbles. An air core occurs in the hydrocyclone when the air‐liquid ratio is more than 1 %. The velocities of air bubbles have a similar flow pattern to the water phase. The axial and tangential velocity differences of the air bubbles at different air‐liquid ratio are greater near the wall and near the core of the hydrocyclone. The measured results show that the size distribution of the air bubbles produced by the air‐liquid mixing pump is beneficial to the process where air bubbles capture oil droplets in the hydrocyclone. These studies are helpful to understand the separation mechanism of a hydrocyclone enhanced by air bubbles.     

CFD Simulations of a Bubble Column Operating in the Homogeneous and Heterogeneous Flow Regimes

Bubble columns are operated either in the homogeneous or heterogeneous flow regime. In the homogeneous flow regime, the bubbles are nearly uniform in size and shape. In the heterogeneous flow regime, a distribution of bubble sizes exists. In this paper, a CFD model is developed to describe the hydrodynamics of bubble columns operating in either of the two flow regimes. The heterogeneous flow regime is assumed to consist of two bubble classes: “small” and “large” bubbles. For the air‐water system, appropriate drag relations are suggested for these two bubble classes. Interactions between both bubble populations and the liquid are taken into account in terms of momentum exchange, or drag‐, coefficients, which differ for the “small” and “large” bubbles. Direct interactions between the large and small bubble phases are ignored. The turbulence in the liquid phase is described using the k‐ϵ model. For a 0.1 m diameter column operating with the air‐water system, CFD simulations have been carried out for superficial gas velocities, U, in the range 0.006–0.08 m/s, spanning both regimes. These simulations reveal some of the characteristic features of homogeneous and heterogeneous flow regimes, and of regime transition.     

The influence of approach velocity on bubble coalescence

Experiments have been carried out in which a cloud of air bubbles has been prevented from rising by downflowing water in a tube. High speed photography revealed an almost complete absence of bubble coalescence. This has been attributed to the large approach velocities of bubbles in the cloud.Further experiments in which a single bubble has been allowed to coalesce with a plane air—water interface have demonstrated the effect more clearly. Two basic types of bubble coalescence have been recognised depending on the approach velocity of the bubbles. At a low approach velocity, bubble coalescence is rapid, but coalescence times are considerably increased at large approach velocities. For pure liquids, a theory is put forward which shows that at low approach velocities film rupture can occur before the approaching bubbles are brought to rest. At large approach velocities the bubbles are brought to rest before rupture occurs. In the latter case bubble bounce can occur and the total coalescence time is thereby considerably increased.Based on observed approach velocities in a stationary bubble cloud, it is suggested that large approach velocities in a bubble column may be an important factor in limiting bubble coalescence.     

Bubble swarm behaviour and gas absorption in non-newtonian fluids in sparged columns

Mean relative gas hold up, slip velocity, bubble size distribution, and volumetric mass transfer coefficient of oxygen were measured in sparged columns of highly viscous non-Newtonian fluids (CMC solutions) as a function of the gas flow rate, and CMC concentration (fluid consistency index k, and flow behaviour index n).By comparison of the measured bubble swarm velocities with those calculated by relations for single bubbles the bubble swarm behaviour was investigated. It could be shown that small bubbles in swarm have higher rising velocities than single bubbles, expecially in highly viscous media. Large single bubbles rise with high velocity due to the change of their shape caused by the swarm of the smaller bubbles. No large bubbles with spherical cap shape could be observed. The volumetric mass transfer coefficient decreases rapidly with increasing CMC-concentration.A comparison of the volumetric mass transfer coefficients with those measured in mechanically agitated vessels indicates, that the performance of sparged columns is comparable with the one of agitated vessels. Because of their lower energy requirement sparged columns are more economical than mechanically agitated vessels. It is possible to improve the performance of sparged columns by the redispersion of large bubbles in a multistage equipment.     

Bubble mechanisms and characteristics at pore scale in a packed-bed reactor

An image processing technique was used to study dominant bubble mechanisms in a two-dimensional packed-bed at pore level under the bubbly flow regime. Bubble breakup and coalescence were identified as dominant mechanisms using a large number of image samples. Two types of coalescence mechanisms were identified that occur due to compression and deceleration associated with the bubbles and three breakup mechanisms were identified that are result of liquid shear force, bubble acceleration, and bubble impact. Data on various two-phase parameters, such as local void fraction, bubble velocity, size, number, and shape were obtained from the images. Results indicated that when a flow regime changed from bubbly to either trickling or pulsing flow, the number of average sized bubbles significantly decreased and the shape of the majority of the bubbles was no longer spherical. Although a mean bubble velocity of all sized bubbles was uniform for given gas and liquid superficial velocities, individual bubble velocities were quite different depending on the bubble location in the pore. The present bubble size distributions were compared with previous studies and the results on bubble size are in general agreement.