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.     

Bubble column research in Japan

Many experimental studies on the bubble column have been reported by Japanese researchers since around 1960. They include studies of bubble behaviour, bubble size distribution, transition from the homogeneous bubbly flow regime to the heterogeneous liquid circulation regime, liquid velocity distribution, longitudinal liquid mixing, hydrodynamic modelling, the gas holdup, and the volumetric coefficient of gas-liquid mass tranfer k_La. Studies covered various modified bubble columns, such as the airlift reactor with an external or internal loop, the packed bubble column, and others. Performance of three-phase bubble columns, which deal with suspensions or emulsions, and their use as bioreactors or chemical reactors were also studied.     

Bubble shape,gas flow and gas–liquid mass transfer in pulp fibre suspensions

Gas–liquid mass transfer in pulp fibre suspensions in a batch‐operated bubble column is explained by observations of bubble size and shape made in a 2D column. Two pulp fibre suspensions (hardwood and softwood kraft) were studied over a range of suspension mass concentrations and gas flow rates. For a given gas flow rate, bubble size was found to increase as suspension concentration increased, moving from smaller spherical/elliptical bubbles to larger spherical‐capped/dimpled‐elliptical bubbles. At relatively low mass concentrations (C_m = 2–3% for the softwood and C_m ? 7% for the hardwood pulp) distinct bubbles were no longer observed in the suspension. Instead, a network of channels formed through which gas flowed. In the bubble column, the volumetric gas–liquid mass transfer rate, k_La, decreased with increasing suspension concentration. From the 2D studies, this occurred as bubble size and rise velocity increased, which would decrease overall bubble surface area and gas holdup in the column. A minimum in k_La occurred between C_m = 2% and 4% which depended on pulp type and was reached near the mass concentration where the flow channels first formed.     

Mass Transfer Characteristics of Syngas Components in Slurry System at Industrial Conditions

The gas‐liquid mass transfer behavior of syngas components, H₂ and CO, has been studied in a three‐phase bubble column reactor at industrial conditions. The influences of the main operating conditions, such as temperature, pressure, superficial gas velocity and solid concentration, have been studied systematically. The volumetric liquid‐side mass transfer coefficient k_La is obtained by measuring the dissolution rate of H₂ and CO. The gas holdup and the bubble size distribution in the reactor are measured by an optical fiber technique, the specific gas‐liquid interfacial area aand the liquid‐side mass transfer coefficient k_L are calculated based on the experimental measurements. Empirical correlations are proposed to predict k_L and a values for H₂ and CO in liquid paraffin/solid particles slurry bubble column reactors.     

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.     

Gas Holdup and Volumetric Mass Transfer Coefficient in a Slurry Bubble Column

The gas holdup, ?, and volumetric mass transfer coefficient, k_La, were measured in a 0.051 m diameter glass column with ethanol as the liquid phase and cobalt catalyst as the solid phase in concentrations of 1.0 and 3.8 vol.‐%. The superficial gas velocity U was varied in the range from 0 to 0.11 m/s, spanning both the homogeneous and heterogeneous flow regimes. Experimental results show that increasing catalyst concentration decreases the gas holdup to a significant extent. The volumetric mass transfer coefficient, k_La, closely follows the trend in gas holdup. Above a superficial gas velocity of 0.04 m/s the value of k_La/? was found to be practically independent of slurry concentration and the gas velocity U; the value of this parameter is found to be about 0.45 s^–1. Our studies provide a simple method for the estimation of k_La in industrial‐size bubble column slurry reactors.     

Study of Hydrodynamics,Mixing and Gas‐Liquid Mass Transfer in a Split‐Rectangular Airlift Reactor

Global hydrodynamic characteristics, liquid mixing and gas‐liquid mass transfer for a 63 L split‐rectangular airlift reactor were studied. Correlations for gas holdup and overall liquid circulation velocity were derived for the air‐water system as a function of the specific power input; these were compared to data and correlations for reactor volumes between 4.7 L and 4600 L. A partial recirculation of small bubbles in the riser was observed when U_gr > 0.03 m/s, which was attributed to the use of a single‐orifice nozzle as the gas phase distributor. The dimensionless mixing time and the overall axial dispersion coefficient were nearly constant for the range of gas flow rates studied. However, values of K_L/d_B were greater than those reported in previous studies and this is caused by the partial recirculation of the gas phase in the riser. While scale effects remain slight, the use of a gas distributor favouring this partial recirculation seems adequate for mass transfer in split‐rectangular airlift reactors.     

Hydrodynamics and mass transfer behavior in multiple draft tube airlift contactors

The draft tube configuration significantly affected the performance of an airlift contactor. The multiple draft tube configuration was demonstrated to give a better gas-liquid mass transfer when compared with a conventional one-draft-tube system. The airlift with a larger number of draft tubes allowed a higher level of bubble entrainment, which rendered a high downcomer gas holdup. This resulted in a higher overall gas holdup in the contactor. Liquid velocity was also enhanced by increasing the number of draft tubes. The ratio between downcomer and riser cross sectional areas, A _d /A _r , had a great effect on the system performance, where a larger A _d /A _r led to a lower downcomer liquid velocity and smaller quantity of gas bubbles being dragged into the downcomer. This resulted in low gas holdup, and consequently, low gas-liquid interfacial mass transfer area, which led to a reduction in the overall volumetric mass transfer coefficient. The presence of salinity in the system drastically reduced the bubble size and subsequently led to an enhancement of gas entrainment within the system. As a result, higher gas holdups and gas-liquid interfacial area were observed, and hence, a higher rate of gas-liquid mass transfer was obtained.     

Experimental Study on Bubble Rising and Descending Velocity Distribution in a Slurry Bubble Column Reactor

To determine bubble rising and descending velocity simultaneously, a BVW‐2 four‐channel conductivity probe bubble parameters apparatus and its analysis are used in gas‐liquid and gas‐liquid‐solid bubble columns. The column is 100 mm in internal diameter and 1500 mm in height. The solid particles used are glass beads with an average diameter of 17.82 μm, representing typical particle size for catalytic slurry reactors. The effects of superficial gas velocity (1.0 cm/s ≤ U_g ≤ 6.4 cm/s), solid holdup (0 % ≤ ?_s ≤ 30 %), and radial location (r/R = 0, 0.4, and 0.7) on bubble velocity distributions are determined. It is found that increasing U_g can increase the velocity of bubbles but do not exert much influence on bubble velocity distribution. Solid holdup mainly affects the distribution of bubble velocity while the radial direction affects bubble velocity distribution only slightly. The ratio of descending bubbles to rising bubbles increases from the bubble column center to the wall. It can be proved experimentally that large bubbles do not always rise faster than small bubbles at higher U_g (for example 6.4 cm/s).     

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.     

CFD analysis of two‐phase turbulent flow in internal airlift reactors

The hydrodynamic performance of three internal airlift reactor configurations was studied by the Eulerian–Eulerian k–ε model for a two‐phase turbulent flow. Comparative evaluation of different drag and lift force coefficient models in terms of liquid velocity in the riser and downcomer and gas holdup in the riser was highlighted. Drag correlations as a function of Eötvös number performed better results in comparison to the drag expressions related to Reynolds number. However, the drag correlation as a function of both Reynolds and Eötvös numbers fitted well with experimental results for the riser gas holdup and downcomer liquid velocity in configurations I and II. Positive lift coefficients increase the liquid velocity and decrease the riser gas holdup, while opposite results were obtained for negative values. By studying the effects of bubble size and their shape, the smaller bubbles provide a lower liquid velocity and a gas holdup. The effects of bubble‐induced turbulence and other non‐drag closure models such as turbulent dispersion and added mass forces were analysed. The gas velocity and gas holdup distributions, liquid velocity in the riser and downcomer, vectors of velocity magnitude and streamlines for liquid phase, the dynamics of gas holdup distribution and turbulent viscosity at different superficial gas velocities for different reactor configurations were computed. The effects of various geometrical parameters such as the draft tube clearance and the ratio of the riser to the downcomer cross‐sectional area on liquid velocities in the riser and the downcomer, the gas velocity and the gas holdup were explored. © 2011 Canadian Society for Chemical Engineering     

装配水平筛板的气升式反应器中气液传质特性的数值模拟

利用Turbulent–Lehr组合模型对装配水平筛板的气升式反应器进行了计算流体力学(CFD)模拟,研究水平筛板对气含率、气泡直径、体积传质系数(kLa)和气液流速的影响。结果表明,筛板对气相的囤积作用和对液相的阻碍作用增加了反应器的整体气含率。筛板对气相的二次均布作用减弱了筛板和液面之间区域的气泡聚并过程,筛板筛孔对气泡的破碎作用产生了大量小于初始直径的气泡,增加了气泡比表面积(a);筛板对液相的阻碍作用提高了筛板附近的气–液相流动速度差,从而提高了该区域的液膜传质系数(kL),强化了反应器内的气液传质效果。     

Studies on the hydrodynamic behavior and mass transfer in a down‐flow jet loop reactor with a coaxial draft tube

The effects of certain pertinent parameters such as gas and liquid flow rates and nozzle position on the behavior of a down‐flow jet loop reactor (DJR) have been studied. The mean residence times of gas and liquid phases and the gas holdup within the reactor have been measured. In addition, the overall volumetric mass transfer coefficient, and the influence of the gas flow rate and the position of the nozzle inside the draft tube on the latter has been determined. Correlations have been presented for the gas holdup and k_La which take into account the length of the draft tube and the nozzle immersion height. The k_La values obtained at different power per unit volume (P/V) values in the DJR used in the present study compare favorably with data presented for stirred tanks and bubble columns in the literature. The liquid residence time distribution (RTD) within the reactor has been studied by tracer analysis for various operating conditions and nozzle immersion height and the results are indicative of the high mixing intensities that can be obtained in such reactions. © 2001 Society of Chemical Industry     

Gas–liquid mass transfer in a new three-phase stirred airlift reactor

The gas–liquid mass transfer performances of a novel three-phase reactor involving both airlift and mechanical stirring have been tested using aqueous solutions of glucose. Stirring in addition to classical airlift leads to an importance increase of k_La. The absolute increase depends mainly on the stirrer speed are not on the gas velocity. A slight effect of the solid loading with a maximum at about 2% (w/v) was observed. Two correlations that show the influence of physical parameters are proposed for both water and glucose solutions. © 1998 SCI     

Volumetric mass transfer coefficients in an internal loop airlift reactor with low‐density particles

The influence of organic additives (propanol, benzoic acid, isoamyl alcohol and carboxymethylcellulose) on the volumetric mass transfer coefficient, k_La, in an internal loop airlift reactor with low‐density particles (nylon‐6 and polystyrene) was investigated. The k_La values increased with increase in superficial gas velocity, U_sg, and decreased with increase in solid loading. A draft tube to reactor diameter ratio, D_R/D, of 0.4 gave maximum k_La values. The addition of benzoic acid and propanol increased the k_La values owing to their coalescence inhibiting characteristics. The addition of isoamyl alcohol decreased k_La, owing to the formation of rigid bubbles and recirculation of small bubbles having a low oxygen content. The k_La values decreased with increase in the concentration of the non‐Newtonian fluid carboxymethylcellulose (CMC). The proposed correlations predicted the experimental data well. Copyright © 2006 Society of Chemical Industry     

COMPREHENSIVE STUDY OF THE GAS HOLD UP AND BUBBLE SIZE DISTRIBUTION IN HIGHLY VISCOUS LIQUIDS II. CMC SOLUTIONS

Following properties of short bubble columns employing CMC solutions (1.0, 1.4 and 2.0%) and perforated plates (0.5, 1.0 and 3.0 mm hole diameters) were determined: relative mean gas hold up, E_G , bubble size distribution, Sauter bubble diameter, d_s , and the specific geometrical bubble surface areas, a′ and a″ due to the “intermediate to large bubbles.” The a″ values were compared with the corresponding volumetric mass transfer coefficients, k_La's and the mass transfer coefficients, k_L , were estimated. The properties of these systems were investigated as function of the superficial gas velocity, W_SG , CMC concentration and aerator type.     

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     

Ozone Mass Transfer in the Mixing Zone of a Confined Plunging Liquid Jet Contactor

In a Confined Plunging Liquid Jet Contactor (CPLJC) a jet of liquid is introduced into an enclosed cylindrical column (downcomer) that generates fine gas bubbles that are contacted with the bulk liquid flow. The region where the liquid jet impinges the receiving liquid and expands to the wall of the downcomer is called the Mixing Zone (MZ). In the MZ, the energy of the liquid jet is dissipated by the breakup of the entrained gas into fine bubbles, and the intense recirculation of the two-phase mixture. The study presented here was undertaken to quantify the ozone-water mass transfer performance of the MZ through the determination of the volumetric mass transfer coefficient, k_La (s^?1), and to produce a model for predicting k_La based on the specific energy dissipation rate. It was found experimentally that k_La in the MZ increased with increasing superficial gas velocity. A maximum experimental k_La value of 0.84 s^?1 was achieved which compares well to other contactors used in water treatment. Such a large k_La value combined with the small volume of the reactor, favorable energy requirements and safety features of the system, suggests that the CPLJC provides an attractive alternative to conventional ozone contactors. The relatively large mass transfer rates were found to be a function of the high gas holdup and fine bubble size generated in the MZ, which results in an almost froth-like consistency. A model based on the specific energy dissipation rate of the water jet, E (kg · m^?1· s^?3), and MZ bubble size was used to predict k_La in the MZ. Using E, the number average bubble size was predicted which was then used to calculate the liquid phase mass transfer coefficient k_L. The bubble size was also used with the predicted mixing zone gas holdup to calculate the specific interfacial area, a (m^?1), which was then combined with k_L to determine a predicted value of k_La. The average deviation between experimental and predicted k_La was 6.2%.     

Gas‐Liquid Mass Transfer in Fibrous Bed Reactor with Counter‐Current Liquid Recycle

In this work, the gas‐liquid mass transfer in a lab‐scale fibrous bed reactor with liquid recycle was studied. The volumetric gas‐liquid mass transfer coefficient, k_La, is determined over a range of the superficial liquid velocity (0.0042–0.0126 m.s^–1), gas velocity (0.006–0.021 m.s^–1), surface tension (35–72 mN/m), and viscosity (1–6 mPa.s). Increasing fluid velocities and viscosity, and decreasing interfacial tension, the volumetric oxygen transfer coefficient increased. In contrast to the case of co‐current flow, the effect of gas superficial velocity was found to be more significant than the liquid superficial velocity. This behavior is explained by variation of the coalescing gas fraction and the reduction in bubble size. A correlation for k_La is proposed. The predicted values deviate within ± 15 % from the experimental values, thus, implying that the equation can be used to predict gas‐liquid mass transfer rates in fibrous bed recycle bioreactors.     

Oxygen transfer in viscous non-Newtonian liquids having yield stress in bubble columns

The oxygen transfer and hydrodynamics in viscous media having a yield stress in bubble columns operated under the slug flow regime were investigated to design an optimum bubble column fermentor for culture media having a yield stress.The gas holdup of escapable bubbles was well estimated by the equation of Nicklin et al. (Trans. Inst. Chem. Eng. 40 (1962) 61), which was modified for the viscous liquid having a yield stress by Terasaka and Tsuge (Chem. Eng. Sci. 58 (2003) 513). The volumetric oxygen transfer coefficient k_La increased with increasing superficial gas velocity and decreasing column diameter under the present conditions. To predict k_La in the non-Newtonian liquids having a yield stress under the operation in slug flow regime, the proposed correlation equation estimated relatively well the experimental k_La.To increase oxygen transfer rate, two types of novel bubble columns were compared with the standard bubble column. The partitioned bubble column presented the better performance than those of the other ones.