Bubble size distribution in bubble columns

The bubble size distributions are measured for the air-water system as a function of air velocity at room temperature in two bubble columns. High speed cinephotography and fiber optic probe techniques are used to measure the bubble size. Our limited measurements suggest that the bubble size may be independent of gas velocity in the range 3.6 to 9.2 cm/s and may be dependent on column diameter with smaller bubbles for narrower columns. The bubble size appears to be smaller at the column wall than at distances away from the wall.     

Population balance modeling of bubble size distributions in a direct-contact evaporator using a sparger model

The study of bubble size distributions in direct-contact evaporators was addressed both theoretically and experimentally. Recently developed models for calculating bubble coalescence and breakage frequencies in isothermal bubble columns were adapted to the population balance equation using the bubble mass as the internal coordinate which was discretized using an expansion of the number density function by impulse functions. A sparger model was developed based on experimental data for a non-coalescing system and using bubble formation models for isothermal and non-isothermal conditions. Bubble size distributions in a direct-contact evaporator operating in the quasi-steady-state regime for four different gas superficial velocities, including the homogeneous and heterogeneous regimes, together with the sparger model, were used for estimating the three empirical parameters from the population balance model, which were observed to be functions of the gas superficial velocity. In all cases considered, the population balance model fitted the experimental data rather well and the regressed parameters exhibit the physically expected behavior with changes in the gas superficial velocity.     

Effects of sparger geometry on the mechanism of flow pattern transition in a bubble column

Electrical resistance tomography (ERT) is used to measure void fraction wave characteristics and to identify flow pattern in a bubble column reactor (0.24 m diameter, 2.75 m height). The effects of column pressure and superficial gas velocities for different sparger geometry and for different flow pattern have been investigated. The ERT sensor can distinguish the void fraction disturbances in different flow regimes with a good clarity. The holdup derived from ERT is in good agreement with the hold-up values measured by pressure transmitters. Different flow regimes have been identified based on void fraction properties and wall pressure fluctuations. The spectral analysis of ERT measurements yields quantitative information, such as a characteristic time and a characteristic frequency of void fraction waves, which are closely related to flow structure in the prevailing regime. The experimental observations are compared with the literature.     

Modeling of the evolution with length of bubble size distributions in bubble columns

Many of the existing methods, for the determination of the specific interfacial area in bubble columns, consider the column in a dynamic equilibrium between bubble coalescence and breaking-up. The aim of this work is to study if this consideration can be considered true for low superficial gas velocities. Two existing models have been chosen from literature in order to predict the break-up [Wang, T., Wang, J., Jin, Y., 2003. A novel theoretical breakup kernel function for bubbles/droplets in a turbulent flow. Chemical Engineering Science 58, 4629-4637] and the coalescence [Lehr, F., Millies, M., Mewes, D., 2002. Bubble size distributions and flow fields in bubble columns. A.I.Ch.E. Journal 48, 2426] rates. In order to confirm the validity of the models, predictions were compared with experimental results obtained by image analysis. Several simulations were performed for different superficial gas velocities and initial bubble size distributions. The column length needed to reach dynamic equilibrium was calculated for each simulation. The results show that the necessary length to reach the dynamic equilibrium does not depend on the shape of the initial distribution, but essentially on its Sauter mean diameter. The necessary length to reach the dynamic equilibrium is very important for low superficial gas velocities. The assumption that the entire column is in dynamic equilibrium is in general not valid. Therefore, the initial Sauter mean diameter and the total column length are important parameters for the determination of the specific interfacial area.     

Bubble size modeling approach for the simulation of bubble columns

The constant bubble size modeling approach (CBSM) and variable bubble size modeling approach (VBSM) are frequently employed in Eulerian–Eulerian simulation of bubble columns. However, the accuracy of CBSM is limited while the computational efficiency of VBSM needs to be improved. This work aims to develop method for bubble size modeling which has high computational efficiency and accuracy in the simulation of bubble columns. The distribution of bubble sizes is represented by a series of discrete points, and the percentage of bubbles with various sizes at gas inlet is determined by the results of computational fluid dynamics (CFD)–population balance model (PBM) simulations, whereas the influence of bubble coalescence and breakup is neglected. The simulated results of a 0.15 m diameter bubble column suggest that the developed method has high computational speed and can achieve similar accuracy as CFD–PBM modeling. Furthermore, the convergence issues caused by solving population balance equations are addressed.     

Spargers for controlled bubble size by means of the multiple slot disperser (MSD)

Sparging technology is crucial in the dispersion of gases in liquids. In this work, it was demonstrated that an effective, controllable sparger can be made by assembling an array of flat parallel slot-nozzles which may offer new options for sparger design and operation. To illustrate the technique, a compact, ‘multiple slot disperser’ (MSD) having a slot width of and a total slot length of 1.26 m was assembled from a series of 5-mm-thick graphite plates. In water containing a low concentration of frother, a dense three-dimensional bubble plume was produced. The MSD generated consistently narrow bubble size distributions with well-defined median bubble diameters in the range 2.6-3.1 mm, equivalent to a range of gas flow rate from 11 to 26 std l/min. The bubble sizes were readily predictable from established single slot bubble size correlations. The method of construction also allows for simple maintenance, repair and replacement of individual components as needed.     

Design and selection of sparger for bubble column reactor. Part II: Optimum sparger type and design

Bubble columns are widely used for conducting gas–liquid and gas–liquid–solid mass transfer/chemical reactions. Sparger is the most important accessory because it decides the bubble size/rise velocity distribution. These, in turn, govern the radial and axial hold-up profiles, the liquid phase flow pattern and hence the performance of bubble columns. In particular, the sparger design is critical if the aspect ratio is low and the sparger design dominates the performance of the bubble column. However, systematic procedure for the selection of sparger design and type are not available in the published literature. This is the specific objective of the present work. In Part I, the performance of different spargers, including the newly developed wheel type of sparger is discussed. Thus the important considerations required for the sparger design are highlighted. The bubble column used in the manufacture of hydrogen peroxide has been considered as a case for illustration.     

The influence of electrolytes on gas hold-up and regime transition in bubble columns

Experiments were conducted in a 0.12-m-in-diameter bubble column to investigate the effect of electrolytes on gas hold-up (ε) and on the regime transition point in bubble columns. Air was used as the dispersed phase and aqueous solutions of three different salts (NaCl, Na₂SO₄ and NaI), as well as double-distilled water, were utilised as the continuous phase, varying the gas superficial velocity (u_G) in the range 0-0.26 m/s. The ε×u_G curves were a function of both the chemical nature and the concentration of the electrolytes. However, similar ε×u_G profiles were obtained regardless of the electrolyte for a given ratio between the concentration in the solution and the critical concentration of the electrolyte for bubble coalescence. This ratio therefore presents itself as a promising modelling parameter to account for the chemical nature of electrolytes. The gas hold-up data were employed to compute the regime transition point according to two different methods, evidencing its non-linear dependence on the concentration of electrolytes in the liquid.     

Bubble columns with fine pore sparger operating in the pseudo-homogeneous regime: Gas hold up prediction and a criterion for the transition to the heterogeneous regime

New experimental data concerning the gas holdup in bubble columns equipped with porous sparger were acquired. The effect of liquid properties and sparger characteristic (i.e., pore size, dimensions) on gas holdup at the pseudo-homogeneous regime has been studied and a correlation regarding the prediction of the transition point from the pseudo-homogeneous to the heterogeneous regime has been proposed and found to be in good agreement with available data. Moreover, a previously proposed correlation [Mouza, A.A., Dalakoglou, G.K., Paras, S.V., 2005. Effect of liquid properties on the performance of bubble column reactors with fine pore spargers. Chemical Engineering Science 60(5), 1465-1475], for the prediction of gas holdup at the homogeneous regime for this type of equipment, has been modified to take into account the effect of the mean pore diameter and it is also found to be in good agreement with published data.     

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.     

Correction of the penetration theory based on mass-transfer data from bubble columns operated in the homogeneous regime under high pressure

A new correction term was developed which allows the classical penetration theory to be applied successfully to k_La data obtained from oblate ellipsoidal bubbles formed in bubble columns operated in the homogeneous regime at various pressures . The correction factor is a function of both the Eötvös number Eo and dimensionless gas density ratio. The new correlation was compared with literature k_La data in 18 pure organic liquids, 14 adjusted liquid mixtures and tap water. In some of the liquids (tetralin, xylene and ethanol) not only air but also other gases (nitrogen, helium and hydrogen) were used. In total, 263 experimental k_La points are fitted with an average relative error of 10.4%.In the theoretical approach for the k_La prediction, the gas-liquid contact time (used in the penetration theory) is defined as the ratio of bubble surface to the rate of surface formation. All further calculations are based on the geometrical characteristics (bubble length and height) of an oblate ellipsoidal bubble. It was found that the new correction factor f_c gradually reduces with the increase of both superficial gas velocity u_G and gas density ρ_G (operating pressure P).     

Bubble Size Distribution in Oil‐Based Bubble Columns

A practical population balance model was used to evaluate the bubble size distribution in a bubble column. In addition, the bubble size distribution in the bubble column was measured at different gas velocities by photography and analysis of the pictures. Four types of liquid, i.e., water and three petroleum‐based liquids, were used in the experiments. The gas phase was air. It was found that the existing models in the literature are not able to satisfactorily predict the experimentally measured bubble size distribution. The model can be corrected by applying a correction factor to the energy dissipation rate. The corrected model fits the experimental bubble size distribution considerably better than the existing models. The variation of this correction factor is reported for different systems at different gas velocities.     

Liquid phase axial mixing in bubble columns operated at high pressures

Liquid phase axial mixing was measured in a 100 mm i.d. bubble column operated in the pressure range of 0.1-0.5 MPa. Water, ethanol and 1-butanol were used as the liquid phase and nitrogen as the gas phase. The temperature and superficial gas velocity were varied in the range of 298-323 K and 0.01-0.21 m/s, respectively. The axial dispersion coefficient increased with an increase in the gas density due to pressure. The temperature had surprisingly a small effect. A CFD model was developed for the prediction of flow pattern in terms of mean velocity and eddy diffusivity profiles. The model was further extended for the prediction of residence time distribution and hence the axial dispersion coefficient (D_L). The predictions of axial dispersion coefficient agree favorably with all the experimental data collected in this work as well as published in the literature. The model was extended for different gas-liquid systems. The predicted values of axial dispersion coefficient were found to agree very well with all the experimental data.     

Liquid-phase backmixing in bubble columns, structured by introduction of partition plates

Installation of perforated sieve plates into a bubble column has the effect of introducing structure into an otherwise chaotic hydrodynamic behaviour. In this study, we focus on the reduction of backmixing of the liquid phase in compartmentalised bubble columns. Liquid-phase residence time distribution (RTD) measurements were carried out in bubble columns with diameters D_T=0.10, 0.15 and 0.38 m with air–water system operating at superficial gas velocities of U_G=0.05–0.4 m/s. Partition sieve plates with open areas of 18.6 and 30.7% were used in the studies. The measured data on RTD were interpreted in terms of an axial dispersion model extended to allow for liquid interchange between compartments. The interchange velocity was found to be strongly dependent on the free area of the plates but practically independent of the column diameter.     

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.     

Bubble characteristics in shallow bubble column reactors

The bubble characteristics have been investigated in an air–water bubble column with shallow bed heights. The effect of bed height, location and the presence of solids on the bubble size, bubble rise velocity and overall and sectional gas holdup are studied over a range of superficial gas velocities. Optimal shallow bed operation relies on the combined entrance and exit effects at the distributor and the liquid bed surface. The gas holdup is found to decrease with an increase in H/D ratio but the effect is diminishing at high H/D ratios. A H/D ratio of 2–4 is found to be suitable for shallow bed operation. The presence of solids causes the formation of larger bubbles at the distributor and the effect is diminishing as the gas velocity is increased.     

The important role of local dispersed phase hold-ups for the calculation of three-phase bubble columns

A computational fluid dynamics model is used to calculate a three-phase (air-water-solid particles) flow in a bubble column. The calculation of multi-phase flows is significantly influenced by the formulation of the inter-phase drag and the modelling of the turbulence. Both are influenced by the dispersed phases. The k-ε turbulence model extended with terms accounting for the bubble-induced turbulence is used to calculate the eddy viscosity of the liquid phase. Bubble-bubble and particle-particle interactions are considered as well as a direct momentum transfer between the two dispersed phases bubbles and solid particles. The local volume fractions of the dispersed phases are considered for the calculation of the drag coefficients between the dispersed phases and the continuous phase. Measured local gas and solid hold-ups as well as measured liquid velocities are compared with the corresponding calculated results. The measured and the calculated results show good agreement.     

Analysis of bubble behaviors in bubble columns using electrical resistance tomography

The distribution of gas holdup, the rise velocity of gas bubble swarm and the Sauter mean bubble size are estimated with a small diameter laboratory scale bubble column using electrical resistance tomography (ERT). The theory of gas disengagement based on ERT methods has been developed for estimations of bubble size and bubble rise velocity. The gas holdups of large bubble swarm and small bubble swarm, the distribution of both bubble size are derived through the analysis of gas disengagement based on the differences of the rise velocity of bubble swarm at the cross-section imaged by electrical resistance tomography. Experimental results are in very good agreement with correlations and conventional estimation obtained using pressure transmitter methods. The proposed methodology can be also used as an analysis tool for quantifying and optimizing the performance of other types of complex reaction systems.