Hydrodynamic Characterization and Mass Transfer Analysis of an In-Line Multi-Jets Ozone Contactor

This research study investigates mixing and ozone mass transfer characteristics of a pilot-scale in-line multi-jets ozone contacting system. The hydrodynamic characteristics of the contactor were studied using a two-dimensional laser flow map particle image velocimetry coupled with planar laser induced fluorescence (PIV/PLIF). The PIV/PLIF system provided a combination of simultaneous whole-field velocity and concentration data in two-phase flows for different operating conditions. All measurements were conducted under a total liquid flow rate of about 10 L/s with gas flow rate ranging from 0.05 to 0.4 L/s. The gas was introduced to the system through a series of side stream injectors. The side injectors were tested under opposing and alternating modes. A mass transfer study was also conducted to estimate the overall mass transfer coefficient under the same operational conditions used for the hydrodynamics study. It was found that for the same number of jets (i.e., same gas flow rate) the liquid dispersion (D_L) was higher when alternating jets were used. Higher ozone mass transfer rates were observed when using opposing jet compared to the same number of alternating jets.     

Oxygen mass transfer coefficient in bubble column slurry reactor with ultrafine suspended particles and neural network prediction

The gas–liquid volumetric mass transfer coefficient was determined by the dynamic oxygen absorption technique using a polarographic dissolved oxygen probe and the gas–liquid interfacial area was measured using dual‐tip conductivity probes in a bubble column slurry reactor at ambient temperature and normal pressure. The solid particles used were ultrafine hollow glass microspheres with a mean diameter of 8.624 µm. The effects of various axial locations (height–diameter ratio = 1–12), superficial gas velocity (u_G = 0.011–0.085 m/s) and solid concentration (ε_S = 0–30 wt.%) on the gas–liquid volumetric mass transfer coefficient k_La_L and liquid‐side mass transfer coefficient k_L were discussed in detail in the range of operating variables investigated. Empirical correlations by dimensional analysis were obtained and feed‐forward back propagation neural network models were employed to predict the gas–liquid volumetric mass transfer coefficient and liquid‐side mass transfer coefficient for an air–water–hollow glass microspheres system in a commercial‐scale bubble column slurry reactor. © 2012 Canadian Society for Chemical Engineering     

Heat transfer in a continuous bubble column

The effective thermal conductivity k_e was measured for a continuous bubble column operating within the bubbly flow regime. k_e was found to be independent of liquid flow rate but strongly dependent on gas flow rate and physical properties of the liquid phase over the ranges 0.1>u_G>4 cm/s, 0>u_L>3.4 cm/s, 0.00096> μ >0.0028 kg/m μ s and 0.047 >σ>0.072 N/m.     

Revisiting Basics of kLa Dependency on Aeration in Bubble Columns: a Is Surprisingly Stable

A comprehensive experimental characterization of a small-scale bubble column bioreactor (60 mL) is presented. Bubble size distribution (BSD), gas holdup, and k_La were determined for different types of liquids, relevant fermentation conditions and superficial gas velocities u_G. The specific interfacial area a and liquid mass transfer coefficient k_L have been identified independent of each other to unravel their individual impact on k_La. Results show that increasing u_G leads to larger bubbles and higher gas holdup. As both parameters influence a in opposite ways, no increase of a with u_G is found. Furthermore, k_L increases with increasing bubble size outlining that improved oxygen transfer is not the result of higher a but of risen k_L instead. The results build the foundation for further simulative investigations.     

Prediction of The Kolmogorov Entropy Derived from CT Data in a Bubble Column Operated under The Transition Regime and Ambient Pressure

The Kolmogorov entropy (KE) algorithm was successfully applied to single source γ‐ray Computed Tomography (CT) data measured by three scintillation detectors in a 0.162 m‐ID bubble column equipped with a perforated plate distributor (163 holes × ?? 1.32 · 10^–3 m). The aerated liquid height was set at 1.8 m. Dried air was used as a gas phase, while Therminol LT (ρ_L = 886 kg m^–3, μ_L = 0.88 · 10^–3 Pa s, σ = 17 · 10^–3 N m^–1) was used as a liquid phase. At ambient pressure, the superficial gas velocity, u_G, was increased stepwise with an increment of 0.01 m s^–1 up to 0.2 m s^–1. Based on the sudden changes in the KE values, the boundaries of the following five regimes were successfully identified: dispersed bubble regime (u_G < 0.02 m s^–1), first transition regime (0.02 ≤ u_G < 0.08 m s^–1), second transition regime (0.08 ≤ u_G < 0.1 m s^–1), coalesced bubble regime consisting of four regions (called 4‐region flow; 0.1 ≤ u_G < 0.12 m s^–1), and coalesced bubble regime consisting of three regions (called 3‐region flow; u_G > 0.12 m s^–1). The KE values derived from three scintillation detectors in the first transition regime were successfully correlated to both bubble frequency and bubble impact. The latter was found to be inversely proportional to the bubble Froude number. The KE model implies that the bubble size in this particular flow regime is a weak function of the orifice Reynolds number (d_b = 7.1 · 10^–3Re₀^–0.05).     

Impinging-Jet Ozone Bubble Column Modeling: Hydrodynamics,Gas Hold-up,Bubble Characteristics,and Ozone Mass Transfer

A transient back flow cell model was used to model the hydrodynamic behaviour of an impinging-jet ozone bubble column. A steady-state back flow cell model was developed to analyze the dissolved ozone concentration profiles measured in the bubble column. The column-average overall mass transfer coefficient, k_La (s^?1), was found to be dependent on the superficial gas and liquid velocities, u_G (m.s^?1) and u_L (m.s^?1), respectively, as follows: k_La?=?55.58 · u_G ^1.26· u_L ^0.08 . The specific interfacial area, a (m^?1), was determined as a = 3.61 × 10³ · u_G ^0.902 · u_L ^?0.038 by measuring the gas hold-up (ε _G?=?4.67 · u_G ^1.11 · u_L ^?0.05 ) and Sauter mean diameter, d_S (mm), of the bubbles (d_S?=?7.78 · u_G ^0.207 · u_L ^{? 0.008} ). The local mass transfer coefficient, k_L (m.s^?1), was then determined to be: k_L?=?15.40 · u_G ^0.354 · u_L ^0.118 .     

Effect of geometrical parameters of draft tubes and clear liquid height on gas holdup in a bubble column for gas dispersion into tubes

The effects of the geometrical parameters of draft tubes and the clear liquid height on the average gas holdup E_G in a 0.16 m I.D. bubble column for gas dispersion into the tubes were experimentally studied in an airtap water system. The gas holdup depended on the superficial gas velocity U(ING), the kinds of gas spargers, the diameter and length of the draft tubes, clearance C_b between the lower end of the draft tube and the bottom of the bubble column, and the clear liquid height H_L. E_G increased with decreasing hole diameter of the gas sparger at a small gas velocity U_G, but did not depend on the kinds of gas spargers at a large U_G. E_G decreased with increasing clear liquid height H_L. The effect of H_L on E_G was well expressed by the modified three-region model. The experimental data of E_G were correlated.     

Mass transfer properties of bubble columns suspending immobilized glucose oxidase gel beads for gluconic acid production

The volumetric gas-liquid oxygen transfer coefficient, k_L a, and the liquid–solid coefficient, k_S, were measured in a 6.7 L external loop airlift bubble column (ELBC), a 2.5 L internal loop airlift (ILBC) and a 2.5 L normal bubble column (NBC) by the steady state method proposed previously using the oxidation of glucose with air catalyzed by glucose oxidase, GO. For an improved and simultaneous determination of k_L a and k_S, GO was entrapped in calcium alginate gel beads together with fine palladium particles instead of catalase to decompose the hydrogen peroxide produced. The gas holdup, ?_G, in each type of bubble column and the liquid circulation velocity, u_L, governing ?_G in the ELBC were also measured to correlate the data on k_La according to the previous correlations proposed for a larger scale of the ELBC, ILBC and NBC. The data on k_L a, k_S, ?_G and u_L (only for the ELBC) in the reaction system were compared to each other for the three types of bubble columns. The results are well predicted by the previous correlations.     

Multiple effects of operating variables on heat transfer in three-phase slurry bubble columns

Characteristics of heat transfer were investigated in pressurized slurry bubble column reactors whose diameter was either 0.051, 0.076, 0.102 or 0.152 m (ID) and 1.5 m in height, respectively. Effects of gas velocity (U _G ), solid contents (S _C ), pressure (P), liquid viscosity (μ _L ) and column diameter (D) on the heat transfer coefficient (h) between the immersed vertical heater and the column were determined. Multiple effects such as UG and D, P and D, μ _L and D, and S _C and D on the value of heat transfer coefficient were discussed. Temperature fluctuations were also measured and analyzed by adapting chaos theory, which was used to explain the effects of operating variables on the heat transfer in the column. The heat transfer coefficient increased with increasing gas velocity, pressure or solid content in the slurry phase, but decreased with increasing liquid viscosity or column diameter. The decrease trend of h with increasing column diameter was somewhat sensitive when the gas velocity was relatively high (U _G ⩾12 cm/s). The effects of column diameter on the h value became almost linear when the operating pressure (P=4−10 kg _f /cm²), liquid viscosity (μ _L =20−38 mPa·s) or solid content in the slurry phase (S _C =10−20 wt%) was relatively high and gas velocity was relatively low, within these experimental conditions. The heat transfer coefficient was well correlated in terms of dimensionless groups as well as operating variables.     

Comparison of Hydrodynamics and Mass Transfer in Airlift and Bubble Column Reactors Using CFD

Computational Fluid Dynamics (CFD) is used to compare the hydrodynamics and mass transfer of an internal airlift reactor with that of a bubble column reactor, operating with an air/water system in the homogeneous bubble flow regime. The liquid circulation velocities are significantly higher in the airlift configuration than in bubble columns, leading to significantly lower gas holdups. Within the riser of the airlift, the gas and liquid phases are virtually in plug flow, whereas in bubble columns the gas and liquid phases follow parabolic velocity distributions. When compared at the same superficial gas velocity, the volumetric mass transfer coefficient, k_La, for an airlift is significantly lower than that for a bubble column. However, when the results are compared at the same values of gas holdup, the values of k_La are practically identical.     

Mass transfer characteristics of low H/D bubble columns

The mass transfer characteristics of 0.2, 0.6 and 1.0 m diameter bubble columns having a low height to diameter ratio (0.6 < H/D < 4) and operated at low superficial gas velocities (0.01 < V_G < 0.08 m/s) were investigated. Different types of spargers were used to study their effect on the column performance. The values of effective interfacial area, a , and volumetric mass transfer coefficient, k_L a , were measured by using chemical methods. The values of a and k_L a were found to vary from 40 to 420 m²/m³ of clear liquid volume and from 0.01 to 0.16 s^?1, respectively, in the range of V_G, and V_L covered in this investigation. The value of the liquid-side mass transfer coefficient, k_L, was found to vary from 3 × 10^?4 to 7 × 10⁴ m/s. The effect of the physical properties of the system on the values of a was also investigated. The height to diameter ratio and the column diameter did not have significant effect on the values of gas holdup, a and k_L a . It was found that the sparger design is not of critical importance, provided multipoint/multiorifice gas spargers are used. The comparative performance of bubble columns having low H/D with horizontal sparged contactors and tall bubble columns has been considered.     

Hydrodynamics and mass transfer in an upward venturi/bubble column combination

The hydrodynamics and mass transfer characteristics of a venturi/bubble column combination were studied at high liquid superficial velocities of up to 0.35 m/s. The gas hold-up was increased by 50% to 150% and the overall volumetric mass transfer coefficient was tripled when the venturi was used as “gas distributor” instead of a porous distributor. A correlation of the overall volumetric mass transfer coefficient (K_La) with the gas hold-up, valid for gas hold-ups as high as 0.3, was proposed for the cylindrical bubble column section. The energy consumption per mole of oxygen transferred was lower than with most distributors and the oxygen transfer rate per unit of reactor volume was higher than in a bubble column with a porous distributor. The venturi/bubble column combination is a compact and efficient system which does not have the operating problems of systems which require internals.     

EFFECT OF STATIC LIQUID HEIGHT ON GAS-LIQUID MASS TRANSFER IN A DRAFT-TUBE BUBBLE COLUMN AND THREE-PHASE FLUIDIZED BED

Experiments are performed under batch-liquid operating conditions to investigate the effect of static liquid height on the gas-liquid mass transfer coefficient (K_La) in a draft-tube bubble column (DTBC) and a draft-tube three-phase fluidized bed (DTFB). In addition, the effects of column diameter, gas-distributor, and draft-tube diameter are studied. The results indicate that for a given system with a porous plate gas-distributor at low superficial gas velocities (<70 m/hr), increasing static liquid height decreases K_La. At high gas velocities, K_La is independent of the static liquid height. For systems with a perforated gas-distributor, there is no effect of static liquid height on K_La. The formation of small dispersed bubbles at low gas velocities in the porous plate distributor system accounts for the considerably high K_La values and the observed effect of liquid height. On the other hand, the formation of large spherical-cap bubbles and the bubble coalescence at high gas velocities reduce the performance of the porous plate distributor system to that of the perforated one.     

Local bubble size distribution,gas–liquid interfacial areas and gas holdups in an up-flow ejector

Particle image velocimetry (PIV) was used to measure local bubble size distributions (BSD), gas–liquid interfacial areas and gas holdups in an up-flow ejector, based on the water–air system with different liquid and gas flow rates under the presence/absence of the swirl body. The results show that the bubble flow patterns are different whether to add the swirl body into the nozzle, especially at low gas flow rate because the bubbles formed “bubble chain” in the ejector with swirl. The mean bubble sizes D₃₂ of the two are both related to the pressure drop between import and export, gas ratio and liquid flow. The interfacial area and D₃₂ are both mainly dependent on the local gas holdups. The mean bubble sizes in the absence of swirl body are smaller than that in the presence of swirl under different operating conditions. The gas holdups and interfacial area are larger with swirl than those without swirl. With the increase of the gas fraction, the differences of D₃₂, a_t and e_G become smaller.     

Gas hold‐up in bubble columns: Operation with concentrated slurries versus high viscosity liquid

The hydrodynamics of bubble columns with concentrated slurries of paraffin oil (density, ρ_L = 790 kg/m³; viscosity, μ_L = 0.0029 Pa·s; surface tension, σ = 0.028 N·m¹) containing silica particles (mean particle diameter d_p = 38 μm) has been studied in columns of three different diameters, 0.1, 0.19 and 0.38 m. With increasing particle concentration, the total gas hold‐up decreases significantly. This decrease is primarily caused by the destruction of the small bubble population. The hold‐up of large bubbles is practically independent of the slurry concentration. The measured gas hold‐up with the 36% v paraffin oil slurry shows remarkable agreement with the corresponding data obtained with Tellus oil (ρ_L = 862 kg/m³; μ_L = 0.075 Pa·s; σ = 0.028 N·m^?1) as the liquid phase. Dynamic gas disengagement experiments confirm that the gas dispersion in Tellus oil also consists predominantly of large bubbles. The large bubble hold‐up is found to decrease significantly with increasing column diameter. A model is developed for estimation of the large bubble gas hold‐up by introduction of an wake‐acceleration factor into the Davies‐Taylor‐Collins relation (Collins, 1967), describing the influence of the column diameter on the rise velocity of an isolated spherical cap bubble.     

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.     

Numerical analysis of particle mixing in a rotating fluidized bed

This paper describes the numerical analysis of particle mixing in a rotating fluidized bed (RFB). A two-dimensional discrete element method (DEM) and computational fluid dynamics (CFD) coupling model were proposed to analyze the radial particle mixing in the RFB. Spherical polyethylene particles (Geldart group B particles) were used as model particles under the assumptions that they were cohesionless and mono-disperse with their diameter of 0.5 mm.The validity of the proposed model was confirmed by the comparison between the calculated degree of particle mixing and the experimental one, which was obtained by measuring the lightness of the recorded image taken by a high-speed video camera. Effects of the operating parameters (gas velocity, centrifugal acceleration, particle bed height, and vessel radius) on the radial particle mixing rate were numerically analyzed. The radial particle mixing rate was found to be strongly affected by the bubble characteristics, especially by the bubble size. The mathematical model for the rate coefficient of particle mixing as functions of operating parameters was empirically proposed. The radial particle mixing rate in a RFB could be well correlated by the three dimensionless numbers: dimensionless acceleration (Ac), bubble Froude number (Fr_b), and dimensionless radius on the surface of particle bed (β_s).     

Treatment of gases and dust using “droplet columns”

The properties of the air-water state diagram, representing the liquid holdups according to gas velocities, in a 0.075 m diameter column are restated. After measurement of the interfacial areas and mass-transfer coefficients, the part of the diagram corresponding to high gas velocities and low liquid contents (10 < U_G < 14 m/s and 0.005 < U_L < 0.04 m/s) was chosen for the treatment of polluted gas streams. Under these conditions, it was shown that a “droplet column” is very efficient for the treatment of gases polluted by acid vapors (SO₂, HCl) and dust (iron oxide, talc, etc.). The cost of energy appeared more favorable than for classical bubble columns.     

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.     

Gas–liquid interfacial area and mass transfer coefficient in a co‐current down flow contacting column

The gas–liquid interfacial area and mass transfer coefficient for absorption of oxygen from air into water, aqueous glycerol solutions up to 1.5% (w/w) and fermentation medium containing glucose up to a 3% concentration were determined in a co‐current down flow contacting column (CDCC; 0.05 m i.d. and 0.8 m length). Experimental studies were conducted using various nozzle diameters at different gas and re‐circulation liquid rates. Specific interfacial area (a) is determined from the fractional gas hold‐up (ε_G) and the average bubble diameter (d_b). Once the interfacial area is determined, the volumetric mass transfer coefficient (k_La) is then used to evaluate the film mass transfer coefficient in the CDCC. The effects of operating conditions and liquid properties on the specific interfacial area were investigated. The values of interfacial area in air–aqueous glycerol solutions and fermentation media were found to be lower than those in the air–water system. As far as experimental conditions were concerned, the values of interfacial area obtained from this study were found to be considerably higher than those of the literature values of conventional bubble columns. The penetration theory is used to interpret the film mass transfer coefficient and results match the experimental k_L data reasonably well. Copyright © 2006 Society of Chemical Industry