Hydrodynamic characterization of a needle sparger rectangular bubble column: Homogeneous flow, static bubble plume and oscillating bubble plume

In this work a detailed experimental hydrodynamic characterization of a needle sparger rectangular bubble column has been performed. The liquid velocity profiles and bubble plume oscillation frequency have been measured by means of laser Doppler anemometry (LDA), and the bubble velocity map by particle image velocimetry (PIV). In this way, the influence of the superficial gas velocity, liquid height and aeration pattern on the column flow structure was analysed. A highly uniform upward flow structure with down flow near the walls was obtained by means of a full-length aeration pattern. This flow structure was preserved even for high gas fractions values. The partial-length aeration patterns with the aerated zone (defined as the aerated width divided by the column width) larger than 0.7 provide a bubble plume and two pure liquid vortical structures in the column bottom, although they are static in nature. With aerated zones lower than 0.6, an oscillating bubble plume is obtained. A non-dimensional analysis of bubble plume oscillation frequency shows a dependence of bubble plume behaviour with the aerated zone. In this way, two different types of bubble plume oscillations, namely confined bubble plume oscillation and free bubble plume oscillation, are introduced and analysed.     

Study of bubble induced flow structure using PIV

An experimental investigation of the flow structure induced by a chain of gas bubbles was carried out in a rectangular bubble column using particle image velocimetry (PIV). It is observed that the bubble rising trajectory changes from one dimension to three dimension as liquid viscosity reduces. The variation of bubble rising trajectory associates with the alternation of bubble motions—with or without oscillatory and rotational motion depending the bubble rising trajectory is 3-D or 1-D. The different behaviors of gas bubbles introduce various instantaneous and averaged liquid flow structures. In general, complex fluid velocity fields present in liquid system of low viscosity where free vortex, cross flow, and irregular circular flow can be observed. The liquid pseudo-turbulence measured in terms of turbulence intensity and Reynolds stress is more intense in liquid of low viscosity. The turbulence is also enhanced by the frequency of bubble formation.     

LDA measurements of liquid velocities in a refractive index matched packed bubble column

LDA has been used to measure liquid velocities in a small-scale bubble column, internal diameter of 50 mm, packed with glass Raschig rings, 10 and 15 mm. A mixture of benzyl-alcohol and ethyl alcohol was index matched against the packing material. A method to separate the signals from liquid and bubbles was developed. It was found that the axial time-averaged liquid velocity was lower than that obtained in empty bubble columns, and that both the time-averaged liquid velocity and the RMS value of the liquid increased for the larger packing size.     

Mass transfer in bubble column for industrial conditions—effects of organic medium, gas and liquid flow rates and column design

Most of available gas-liquid mass transfer data in bubble column have been obtained in aqueous media and in liquid batch conditions, contrary to industrial chemical reactor conditions. This work provides new data more relevant for industrial conditions, including comparison of water and organic media, effects of large liquid and gas velocities, perforated plates and sparger hole diameter.The usual dynamic O₂ methods for mass transfer investigation were not convenient in this work (cyclohexane, liquid circulation). Steady-state mass transfer of CO₂ in an absorption-desorption loop has been quantified by IR spectrometry. Using a simple RTD characterization, mass transfer efficiency and k_La have been calculated in a wide range of experimental conditions.Due to large column height and gas velocity, mass transfer efficiency is high, ranging between 40% and 90%. k_La values stand between 0.015 and and depend mainly on superficial gas velocity. No significant effects of column design and media have been shown. At last, using both global and local hydrodynamics data, mass transfer connection with hydrodynamics has been investigated through k_La/ε_G and k_La/a.     

Investigation of the unsteady two-phase flow with small bubbles in a model bubble column using phase-Doppler anemometry

The unsteady two-phase flow of water laden with small air bubbles in a model bubble column is investigated experimentally. Phase-Doppler anemometry (PDA) is used for measuring the velocities of water and bubbles. The measured sizes of reflecting tracers in the water and of the air bubbles are used to discriminate between water and bubble data. The investigations are focussed on the unsteady behaviour of the flow and on the interaction between the two phases. The measurement of relative (slip) velocities between bubbles and water reveals information about the dynamic behaviour of the two-phase system under the action of buoyancy on the disperse phase. The evaluation of time series of bubble velocities yields insight into typical frequencies at which the flow fluctuates. It is shown that, at all locations in the flow field, the velocity probability density functions of bubbles and liquid can be described by two superimposed Gaussian functions. The bubbles belonging to the two Gaussians exhibit different slip velocities. The probability for the occurrence of bubble collisions is quantified on the basis of the PDA data.     

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.     

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.     

A gas holdup model for cocurrent air-water-fiber bubble columns

A gas holdup model is developed for cocurrent air-water-fiber bubble column flows using the drift-flux model. The model coefficients are estimated using a nonlinear least square method and systematically acquired experimental data. The model correlates gas holdup with superficial gas and liquid velocity, and fiber type and mass fraction. The model reproduces most experimental data within ±10% error and all but 3 of the 3839 experimental data points within ±15% error. It also accurately predicts air-water bubble column gas holdup data; these data were not used in estimating the model coefficients. The physical implications of the model coefficients are also discussed.     

Effect of fiber length distribution on gas holdup in a cocurrent air-water-fiber bubble column

An experimental investigation is reported on the effect of fiber length distribution on gas holdup in a cocurrent air-water-fiber bubble column. Different combinations of 1 and 3 mm Rayon fibers are used to simulate different fiber length distributions. At a constant total fiber mass fraction, gas holdup generally decreases with increasing mass fraction of the 3 mm Rayon fiber while other conditions remain constant. Crowding factors estimated using four different methods (N_c=N_c,A, , N_c,L, and N_c,M) and the parameters and are tested on their performance to quantify the overall effects of fiber mass fraction and fiber length and its distribution on gas holdup. and provide the best characterization of the fiber effects on gas holdup in the cocurrent air-water-fiber bubble column. The crowding factor estimated using the model-based average fiber length (N_c,M) also provides a good characterization and is better than the other crowding factor definitions.     

Model validation for low and high superficial gas velocity bubble column flows

In the present work, gas-liquid flow dynamics in a bubble column are simulated with CFDLib using an Eulerian-Eulerian ensemble-averaging method in a two-dimensional Cartesian system. The two-phase flow simulations are compared to experimental measurements of a rectangular bubble column performed by Mudde et al. [1997. Role of coherent structures on Reynolds stresses in a 2-D bubble column. A.I.Ch.E. Journal 43, 913-926] and a cylindrical bubble column performed by Rampure et al. [2003. Modeling of gas-liquid/gas-liquid-solid flows in bubble columns: experiments and CFD simulations. The Canadian Journal of Chemical Engineering 81, 692-706] for low and high superficial gas velocities, respectively. The objectives are to obtain grid-independent numerical solutions using CFDLib to reconcile unphysical results observed using FLUENT with increasing grid resolutions [Law, D., Battaglia, F., Heindel, T.J., 2006. Numerical simulations of gas-liquid flow dynamics in bubble columns. In: Proceedings of the ASME Fluids Engineering Division, IMECE2006-13544, Chicago, IL], and to validate computational fluid dynamics (CFD) simulations with experimental data to demonstrate the use of numerical simulations as a viable design tool for gas-liquid bubble column flows. Numerical predictions are presented for the local time-averaged liquid velocity and gas fraction at various axial heights as a function of horizontal or radial position. The effects of grid resolution, bubble pressure (BP) model, and drag coefficient models on the numerical predictions are examined. The BP model is hypothesized to account for bubble stability, thus providing physical solutions.     

Experimental study of mass transfer in a dense bubble swarm

We consider the liquid-side mass transfer coefficient k_L in a dense bubble swarm for a wide range of gas volume fraction (0.45%≤α_G≤16.5%). The study is performed for an air–water system in a square column. Bubble size, shape and velocity have been measured for different gas flow rates by means of a high speed camera. Gas volume fraction and bubble velocity have also been measured by a dual-tip optical probe. Both of these measurements show that the bubble vertical velocity decreases when increasing α_G in agreement with previous investigations. The mass transfer is measured from the time evolution of the dissolved oxygen concentration, which is obtained by the gassing-out method. The mass transfer coefficient is found to be very close to that of a single bubble provided the bubble Reynolds number is based on the average equivalent diameter 〈d_eq〉 and the vertical slip velocity 〈V_z〉.     

Numerical study on gas and liquid slugs for Taylor flow in a T-junction microchannel

The rapid development of microfabrication techniques creates new opportunities for applications of microchannel reactor technology in chemical reaction engineering. The extremely large surface-to-volume ratio and the short transport path in microchannels enhance heat and mass transfer dramatically, and hence provide many potential opportunities in chemical process development and intensification. Multiphase reactions involving gas/liquid reactants with a solid as a catalyst are ubiquitous in chemical and pharmaceutical industries. The hydrodynamics of the flow affects the reactor performance significantly; therefore it plays a prominent role in reactor design. For gas/liquid two-phase flow in a microchannel, the Taylor slug flow regime is the most commonly encountered flow pattern. The present study deals with the numerical simulation of the Taylor flow in a microchannel, particularly on gas and liquid slugs. A T-junction empty microchannel with varying cross-sectional width (0.25, 0.5, 0.75, 1, 2 and 3 mm) served as the model micro-reactor, and a finite volume based commercial computational fluid dynamics (CFD) package, FLUENT, was adopted for the numerical simulation. The gas and liquid slug lengths at various operating and fluid conditions were obtained and found to be in good agreement with the literature data. Several correlations in the T-junction microchannel were developed based on the simulation results. The slug flows for other geometries and inlet conditions were also studied.     

X-ray computed tomography in large bubble columns

X-ray computed tomography (CT) is used to explore the differences in a semi-batch bubble column operated at superficial gas velocities of U_g=3, 10, and 18 cm/s. Air-water or air-water-cellulose fiber systems comprise the multiphase flow, and the bubble column has a 32.1 cm internal diameter. A CT image of a phantom object composed of several air-filled tubes immersed in water is used to identify several characteristic features of the X-ray CT system. CT images are then compared between air-water and air-water-cellulose fiber systems. When the fiber mass fraction is 0.1%, gas holdup is slightly higher than that of the air-water system in the column center and near the column wall. In 1.0% cellulose fiber slurries, gas holdup is lower than that of air-water results at all radial positions.     

Unsteady motion of a single bubble in highly viscous liquid and empirical correlation of drag coefficient

The unsteady motion of single bubbles rising freely in a quiescent liquid with high viscosity was measured using a CCD (charge coupled device) camera. Sequences of the recorded frames were digitized and analyzed using image analysis software and the measurements of the acceleration and steady motion of bubbles were obtained. The total drag coefficient was calculated from the accelerating motion to the steady motion with the added mass force and history force included. In virtue of dimensional analysis, the total drag coefficient of single bubbles is correlated as a function of the acceleration number, Archimedes number and Reynolds number based on the equivalent bubble diameter. The proposed correlation represents very well the experimental data of the total drag force in a wide range covering both unsteady accelerating motion and steady motion. The combined added mass and history force coefficient accounting for the accelerating effect on single bubbles was evaluated and correlated.     

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.     

Identification of flow regime transitions in an annulus sparged internal loop airlift reactor based on higher order statistics and Winger trispectrum

Higher order statistics and Wigner higher order moment spectra were used to extract useful flow regime characteristics from wall pressure fluctuation signals in an annulus sparged internal loop airlift reactor. It is found that the pressure fluctuation in the airlift reactor is a typical nonlinear and non-stationary process, which exhibits different frequency characteristics depending on flow regimes. Analysis methods based on bispectrum and Wigner trispectrum are powerful tools to reveal frequency characteristics of pressure signals. To identify flow regime transitions in the reactor, two new characteristic quantities, namely average bispectrum and generalized average frequency, are defined from bispectrum and Wigner trispectrum of the pressure signal, respectively. Two flow regime transition points corresponding to three flow regimes in the reactor are successfully detected by using these two characteristic quantities.     

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.