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.     

Heat transfer metrology issues in two‐phase bubble column reactors

The metrology and the impact of various parameters and operating conditions on the bulk‐to‐tube heat transfer coefficients in two‐phase bubble columns are investigated on a small‐scale mock‐up. It is shown that (1) quasi‐adiabatic conditions can be reached in the column; (2) the bulk‐to‐tube heat transfer coefficients for each U‐tube downward and upward sections may or may not differ significantly, depending on the way uncertainty of the measurements is estimated; (3) using the different measurements and uncertainty estimates for given conditions, a mean heat transfer coefficient over all tubes is estimated within ±5%. The consequences for bulk‐to‐tube heat transfer coefficient prediction in a larger column are discussed.     

Hydrokinetic optimization of commercial scale slurry bubble column reactor

An in‐depth numerical study has been carried out to investigate a high‐pressure commercial scale (2–8 m diameter, 30–40 m in height) slurry bubble column reactor. Typical superficial gas velocities are in the range of 0.5–3 m/s, and overall vapor holdups are in the range of 0.45–0.85. The study revealed that steady compartmental reaction models do not match plant data when reaction time constants are fast. Also, off‐the‐shelf commercial computational fluid dynamics codes do not produce useful information about a reactive column of this scale without first validating the model using data “anchors” from full‐scale operational columns. Important measures include both transient and time‐averaged profiles, integrals, and extrema of vapor holdup and reactants. Reactor designs based on this study show both improved productivity and product quality, allowing record production from existing plants along with substantial capital scope reduction for new plants. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Liquid velocity in a high‐solids‐loading three‐phase external‐loop airlift reactor

BACKGROUND: Airlift reactors are of interest for many different processes, especially for three‐phase systems. In this study the behavior of a high‐loading three‐phase external‐loop airlift reactor was examined. In particular, the effect of parameters such as airflow rate (riser superficial gas velocities between 0.003 and 0.017 m s^?1), solids loading (up to 50%, v/v) on liquid circulation velocity in the air‐water‐alginate beads system as a crucial hydrodynamic parameter was studied. RESULTS: It was observed that increase of the airflow rate resulted in increase of the liquid velocity in the system. The same result but less pronounced was observed by introducing small amounts of solid particles up to 7.5% v/v. However, further introduction of solids caused decrease of the liquid velocity. Laminar regime for the liquid circulation was observed for low gas velocities. Minimum gas velocities for recirculation initiation in the reactor were determined for all solid loadings and linear dependence on the solid content was found. Gas holdups for the three‐phase system were larger than for the two‐phase system in all experiments. A simple model for predicting the liquid circulation velocity in the three‐phase system with high solid loading of low‐density particles was developed. This model is based on the viscosity of integrated medium (solid + liquid) which is a new aspect to analyze this phenomenon. CONCLUSIONS: The developed model shows very good agreement with the experimental results for all solid loadings. It also includes the influence of reactor geometry on the liquid circulation velocity thus enabling optimization. Copyright © 2012 Society of Chemical Industry     

CFD simulations of hydrodynamic/thermal coupling phenomena in a bubble column with internals

CFD simulations have been carried out in a full three‐dimensional, unsteady, Eulerian framework to simulate hydrodynamic/thermal coupling in a bubble column with internals. A first part of the study, dedicated to the hydrodynamic/thermal coupling in liquid single‐phase flows, showed that assuming constant wall temperature on the internals constitutes a reasonable approximation in lieu of comprehensive simulations encompassing shell flow and coolant flow together. A second part dealing with the hydrodynamics of gas–liquid flows in a bubble column with internals showed that a RNG k–ε turbulence model formulation accounting for gas‐induced turbulence was a relevant choice. The last part used these conclusions to build a hydrodynamic/thermal coupling model of a gas–liquid flow in a bubble column with internals. With a per‐phase RNG k–ε turbulence model and assuming constant wall temperature, it was possible to simulate heat transfer phenomena consistent with experimentally measured heat transfer coefficients. © 2010 American Institute of Chemical Engineers AIChE J, 2010     

Bubble splitting under gas–liquid–liquid three‐phase flow in a double T‐junction microchannel

Gas–aqueous liquid–oil three‐phase flow was generated in a microchannel with a double T‐junction. Under the squeezing of the dispersed aqueous phase at the second T‐junction (T2), the splitting of bubbles generated from the first T‐junction (T1) was investigated. During the bubble splitting process, the upstream gas–oil two‐phase flow and the aqueous phase flow at T2 fluctuate in opposite phases, resulting in either independent or synchronous relationship between the instantaneous downstream and upstream bubble velocities depending on the operating conditions. Compared with two‐phase flow, the modified capillary number and the ratio of the upstream velocity to the aqueous phase velocity were introduced to predict the bubble breakup time. The critical bubble breakup length and size laws of daughter bubbles/slugs were thereby proposed. These results provide an important guideline for designing microchannel structures for a precise manipulation of gas–liquid–liquid three‐phase flow which finds potential applications among others in chemical synthesis. © 2017 American Institute of Chemical Engineers AIChE J, 63: 376–388, 2018     

CFD simulations to portray the bubble distribution and the hydrodynamics in an annulus sparged air‐lift bioreactor

Air‐lift bioreactor (ABR) are applied in the pharmaceutical industry for those processes of low‐oxygen demand. The characteristics of the air bubbles in an ABR are important since they influence the overall performance of the bioreactor. In this article, two‐dimensional CFD simulations are carried out to portray the flow characteristics in an annulus‐sparged ABR with a low‐aeration rate, which is a scaled geometry of a 200 m³ annulus‐sparged ABR for the L ‐ascorbic acid manufacturing. VOF model is applied firstly to capture the flow behaviour of the bubbles in the bioreactor, and the distributions of the bubbles and local mean bubble sizes are approximately obtained. With the local mean bubble size distribution, Eulerian model is adopted to investigate the performance of the bioreactor. The liquid velocity and the gas holdup are further discussed. The results show the existence of the local back mixing as well as the mal‐distributions of the velocity and the gas holdup along the radial direction in the annulus‐sparged ABR.     

Studies on Liquid—Gas and Three‐Phase Fluidized Beds with Pulsating Air Flows

The effects of air‐flow pulsation and water and air flowrates on the hydrodynamics of liquid—gas and three‐phase fluidized beds containing 3‐mm glass beads have been studied in a 90‐mm i.d. column. Under steady‐flow conditions, both types of bed contained a relatively large number of small bubbles. With a pulsing air flow, however, a smaller number of much larger bubbles or slugs were formed. This was attributed to different mechanisms of bubble formation at the distributor. Variations in phase holdup were explained in terms of the effects of the operating parameters on the bubble characteristics.     

Gas phase flow in bubble columns: A convective phenomenon

The axial dispersion model has been commonly used to describe gas phase flow in bubble columns. Scatter in dispersion coefficients reported to date may be a result of the misuse of the axial dispersion model when a convective model would be more appropriate. Using simple tests with radioactive tracer response curve moments, convective and dispersive behaviours are differentiated. A convective model is presented. The model fits both tracer response curves and mean gas velocity well in both the bubbly and churn turbulent flow regimes, and may be used as a technique to calculate bubble rise velocity distributions.     

Hydrodynamic feedback on bubble breakup at a T‐junction within an asymmetric loop

Bubble breakup at a microfluidic T‐junction by taking into consideration the hydrodynamic feedback at the downstream channels is presented. Experiments are conducted in square microchannels with 400 μm in width. The splitting ratio of the bubble size in the bifurcations varies nonmonotonically with the flow rate ratio of gas/liquid phases, and it is also affected by the liquid viscosity. A critical size of the mother bubble determines the variation trend of the splitting ratio of bubble size with flow rates of both phases and the liquid viscosity, which is related to the different breakup mechanisms for long and short bubbles at the junction and the different additional resistances induced by long and short bubbles in downstream channels. A theoretical model is proposed to predict the tailoring size of bubbles at the T‐junction by taking into account of the additional resistance in the presence of bubbles in downstream channels. © 2014 American Institute of Chemical Engineers AIChE J, 60: 1920–1929, 2014     

Characteristics of a countercurrent reciprocating plate bubble column. I. Holdup,pressure drop and bubble diameter

The hydrodynamics of countercurrent air/water flow in a 5 cm diameter reciprocating plate bubble column have been studied; the plates contained 14 mm diameter perforations and had a fractional open area of 0.57. The ranges of superficial velocities of air and water were respectively 0-0.99 cm/s and 0-3.95 cm/s. The stroke was in most cases 4.5 cm and the reciprocation frequency was in the range 0–5 Hz. The pressure drops were measured in the absence of reciprocation for single phase and two phase flow conditions. Pressure fluctuations and time-averaged pressure drops were measured with plate reciprocation under single and two-phase conditions. The results were described in terms of the simple quasi-steady state model; the effective orifice coefficients of the perforations were within the range 0.4 to 0.97 depending on the reciprocation conditions. The Sauter mean diameters of the bubbles decreased with agitation; they were about twice the values predicted from an earlier correlation developed for liquid-liquid systems. The gas holdups were also substantially greater than predicted from correlations based on liquid-liquid systems. Both these effects were explained as due to the tendency for bubbles to cluster in the plate region.     

Prediction of dispersed phase holdup in a modified multi‐stage bubble column scrubber

This paper reports on the experimental investigation carried out to evaluate fractional dispersed phase holdup for a gas‐liquid mixture in a modified multi‐stage bubble column (with contraction and expansion disks), which has been conceived, designed and fabricated as a wet scrubber for control of air pollution; in addition it has versatile use as a gas‐liquid contactor in chemical process industries. A correlation developed for predicting fractional dispersed phase holdup has been found to be encouraging and highly significant from statistical analysis.     

Three‐component solids velocity measurements in the outlet section of a riser

Coincident (simultaneous) three‐component particle velocity measurements performed using two laser Doppler anemometry probes at the outlet section of a 9 m high cylindrical riser are for the first time presented for dilute flow conditions. Near the blinded extension of the T‐outlet a solids vortex is formed. Particle downflow along the riser wall opposite the outlet tube is observed, which is restricted to higher riser heights at higher gas flow rates. Increased velocity fluctuations are observed in the solids vortex and downflow region as well as at heights corresponding to the outlet tube. Contrary to the rest of the riser, in the downflow region time and ensemble velocity averages are not equal. Given the local bending of the streamlines, axial momentum transforms to radial and azimuthal momentum giving rise to the corresponding shear stresses. Turbulence intensity values indicate the edges of the downflow region. © 2016 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 62: 3575–3584, 2016     

Relationship between Riser and Downcomer Gas Holdups in Semi‐batch Circulating Bubble Columns

Gas holdups were measured in a circulating bubble column (CBC) using air‐water system with various additives. The liquid volume in the gas separator affected downcomer gas holdup up to volume ratio equal to 6%. Presence of surface active agents, alcohols, solid particles had little effects on the gas holdups ?_gd — a?_gr relationships. The gas holdups relationships are theoretically related by a single parameter model ?_gd = a?_gr — a + 1, and empirically by two parameter model ?_gd = a?_gr — constant. The theoretical model yielded a values ? 0.45, while the empirical model yielded a values ? 1.16. Analogy between the two models (Contreras et al., 1998) will lead to incorrect conclusions.     

Bubble characteristics of a scaled‐down three‐phase fluidized bed

Experiments were carried out in a 0.29 m diameter column to provide information on bubble characteristics relevant to an industrial three‐phase fluidized bed reactor. Using a specially built electrical conductivity probe, together with a series of computer programs to gather and filter the data, local bubble properties (rise velocities, pierced chord lengths, frequencies, and direction of travel) were measured. The data are compared to existing correlations. Since the operating conditions were within the dispersed bubble regime, uniform bubbling was anticipated. However, the actual system indicated considerable non‐uniformity, with individual bubbles, large swarms of bubbles and particles traveling at all angles. As a result, it is shown that bubble rise velocities are not necessarily a good measure of the bubble properties and must be vectorized to provide a good means of estimating gas hold‐up.     

Influence of the lift force closures on the numerical simulation of bubble plumes in a rectangular bubble column

Numerical Eulerian-Eulerian simulations of the unsteady gas-liquid flow in a centrally aerated two-dimensional bubble column were carried out in order to understand the effect of different formulations of the lift force coefficient (C_L) on the computational results. Three different values of the superficial gas velocity (U_G=2.4, 12.0 and 21.3 mm s⁻¹) that ensure the existence of different flow regimes were experimentally and computationally studied. The validation of the simulated results was based on visual observations and measurements of the global gas hold-up (ε_G) and the plume oscillation period (POP). The results presented reveal that, at U_G=12.0 and 21.3 mm s⁻¹, using C_L<0 results in under- and over-estimation of the ε_G and POP, respectively. On the other hand, taking C_L>0 does not affect the POP while it leads to increasingly higher ε_G values, which are different from those experimentally reported. At U_G=2.4 mm s⁻¹, the effect of the lift force is not so evident, although it slightly improves the prediction of experimental values. Particularly interesting is the case of C_L>0.4 at U_G=21.3 mm s⁻¹, producing a non-symmetric bubble plume oscillation. Since using Tomiyama's lift coefficient correlation does not improve the results, including the lift force into the simulation of bubble plumes is not recommended.