INVISCID LIQUID CIRCULATION IN BUBBLE COLUMNS

For circulation in axi-symmetric (cylindrical) bubble columns, the recently developed mathematical model^25,26 has been used along with the criterion of minimum circulation strength to determine the height of each circulation cell in a tall column. This is then used to derive a theoretical expression, first of its kind, for gas hold-up inside a bubble column. The predictions of this equation as well as the equation derived here for axial liquid velocity at column axis have been compared with available data and the comparison is found to be excellent for both the variables. An explicit relation is derived for the average liquid circulation velocity. The model is also used to derive an expression for liquid axial dispersion coefficient which compares almost exactly with Deckwer et al.'s⁴ correlation.

For circulation in two-dimensional bubble columns a new mathematical model is developed. The predictions of bubble envelope shape and bubble envelope area compare well with published data. The predictions of number of circulation cells in the horizontal direction also compare well with published data.     

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.     

Liquid phase mixing in bubble columns with Newtonian and non-Newtonian fluids

A hydrodynamic model for the liquid phase in bubble columns is developed. The proposed model for fully developed turbulent Newtonian and non-Newtonian fluids is based on an energy balance and the mixing length theory. The predictions of the model are in reasonably good agreement with data on liquid velocity at the column axis and the axial dispersion coefficient. The liquid velocity data in an inverted conical bottom gas—liquid column contactor have also been measured. They are correlated by the proposed model reasonably well.     

Mass transfer in bubble columns packed with motionless mixers

The cocurrent upward mode was employed to absorb pure oxygen into water in bubble columns packed with Koch (Sulzer) motionless mixers. The liquid-side volumetric mass transfer coefficient, K_La, in the packed bubble column was found to be always larger than that in the unpacked bubble column. In the range of liquid velocities from 6.7 cm/sec to 39.9 cm/sec, the value of K_La in the packed bubble column increased with the increasing liquid velocity while that in the unpacked bubble column was almost independent of the liquid velocity. The equation of the formK_La= mν_l^β? was successfully adopted to correlate the K_La data.     

CFD simulation of bubble columns incorporating population balance modeling

A computational fluid dynamics (CFD)-code has been developed using finite volume method in Eulerian framework for the simulation of axisymmetric steady state flows in bubble columns. The population balance equation for bubble number density has been included in the CFD code. The fixed pivot method of Kumar and Ramkrishna [1996. On the solution of population balance equations by discretization—I. A fixed pivot technique. Chemical Engineering Science 51, 1311-1332] has been used to discretize the population balance equation. The turbulence in the liquid phase has been modeled by a k-ε model. The novel feature of the framework is that it includes the size-specific bubble velocities obtained by assuming mechanical equilibrium for each bubble and hence it is a generalized multi-fluid model. With appropriate closures for the drag and lift forces, it allows for different velocities for bubbles of different sizes and hence the proper spatial distributions of bubbles are predicted. Accordingly the proper distributions of gas hold-up, liquid circulation velocities and turbulence intensities in the column are predicted. A survey of the literature shows that the algebraic manipulations of either bubble coalescence or break-up rate were mainly guided by the need to obtain the equilibrium bubble size distributions in the column. The model of Prince and Blanch [1990. Bubble coalescence and break-up in air-sparged bubble columns. A.I.Ch.E. Journal 36, 1485-1499] is known to overpredict the bubble collision frequencies in bubble columns. It has been modified to incorporate the effect of gas phase dispersion number. The predictions of the model are in good agreement with the experimental data of Bhole et al. [2006. Laser Doppler anemometer measurements in bubble column: effect of sparger. Industrial & Engineering Chemistry Research 45, 9201-9207] obtained using Laser Doppler anemometry. Comparison of simulation results with the experimental measurements of Sanyal et al. [1999. Numerical simulation of gas-liquid dynamics in cylindrical bubble column reactors. Chemical Engineering Science 54, 5071-5083] and Olmos et al. [2001. Numerical simulation of multiphase flow in bubble column reactors: influence of bubble coalescence and breakup. Chemical Engineering Science 56, 6359-6365] also show a good agreement for liquid velocity and gas hold-up profiles.     

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.     

湍动浆态床流体力学研究(Ⅱ)轴向浆料速度的径向分布

对高气速、大塔径条件下的湍动浆态床轴向速度进行了实验测定,提出了简化的一维流体力学模型。实验结果表明：固含率对浆料速度的径向分布影响不明显;采用中心速度和塔径作为参考量,则不同条件下的量纲1流速分布相似;高气速下中心速度随塔径的变化可近似用Nottenkamper关联式描述。提出了一维流体力学模型预测浆态床液速分布,模型计算结果与实验数据符合较好,模型较好地反映了浆态床反应器内液速分布随表观气速、固含率及塔径变化的规律。     

Mathematical modeling of gas-liquid mass transfer rate in bubble columns operated in the heterogeneous regime

In this paper, a multi-scale approach is followed to study gas-liquid mass transfer in bubble columns. First, a single bubble of equivalent diameter d is considered. Its morphology and its gas to liquid relative velocity are related to the bubble diameter through the use of known correlations. Then, the gas-liquid mass transfer between the bubble and the surrounding liquid is studied theoretically. An equation describing the transport of the transferred species in the viscous boundary layer around the bubble is solved. In a second step, a bubble column of 6-10 m height is studied experimentally. The gas phase in the column is characterized experimentally by means of a gammametric technique. Finally, the two studies are linked, yielding a 1D mathematical model able to predict the gas-liquid mass transfer rate in a bubble column operated in the heterogeneous regime.     

Sectionalised bubble column: Gas hold-up and wall side solid—Liquid mass transfer coefficient

The effect of gas sparger type, clear liquid height, liquid viscosity and addition of electrolyte on fractional gas hold-up in a 0.38 m i.d. sectionalised bubble column (SBC) was studied for superficial gas velocities ranging from 0.03 to 0.15 m s⁻¹. A study of the wall side solid–liquid mass transfer coefficient, k_SL, has been made by considering the dissolution of copper in acidic dichromate solutions. The variation of k_SL with axial distance in a section for different gas velocities has been explained. A comparison between the performance of sectionalised bubble columns and conventional bubble columns has been presented. In some cases, the SBC may prove to be superior to bubble columns in practice.     

CFD simulations of bubble column reactors: 1D, 2D and 3D approach

CFD simulations have been carried out for the predictions of flow pattern in bubble column reactors using 1D, 2D and 3D k-ε models. An attempt has been made to develop a complete correspondence between the operation of a real column and the simulation. Attention has been focused on the cylindrical bubble columns because of their widespread applications in the industry. All the models showed good agreement with the experimental data for axial liquid velocity and the fractional gas hold-up profiles. However, as regards to eddy diffusivity, only the 3D model predictions agree closely with the experimental data.The CFD model has been extended for the estimation of an axial dispersion coefficient (D_L) using 1D, 2D and 3D models. Excellent agreement was found only between the experimental values and the 3D predictions. The 1D and 2D simulations, however, yielded D_L values, which were lower by 25-50%. For this, a mechanistic explanation has been provided.     

The effect of gas‐liquid counter‐current operation on gas hold‐up in bubble columns using electrical resistance tomography

BACKGROUND: In order to improve the performance of a counter‐current bubble column, radial variations of the gas hold‐ups and mean hold‐ups were investigated in a 0.160 m i.d. bubble column using electrical resistance tomography with two axial locations (Plane 1 and Plane 2). In all experiments the liquid phase was tap water and the gas phase air. The superficial gas velocity was varied from 0.02 to 0.25 m s^?1, and the liquid velocity varied from 0 to 0.01 m s^?1. The effect of liquid velocity on the distribution of mean hold‐ups and radial gas hold‐ups is discussed. RESULTS: The gas hold‐up profile in a gas–liquid counter‐current bubble column was determined by electrical resistance tomography. The liquid velocity slightly influences the mean hold‐up and radial hold‐up distribution under the selected operating conditions and the liquid flow improves the transition gas velocity from a homogeneous regime to a heterogeneous regime. Meanwhile, the radial gas hold‐up profiles are steeper at the central region of the column with increasing gas velocity. Moreover, the gas hold‐up in the centre of the column becomes steeper with increasing liquid velocity. CONCLUSIONS: The value of mean gas hold‐ups slightly increases with increasing downward liquid velocity, and more than mean gas hold‐ups in batch and co‐current operation. According to the experimental results, an empirical correlation for the centreline gas hold‐up is obtained based on the effects of gas velocity, liquid velocity, and ratio of axial height to column diameter. The values calculated in this way are in close agreement with experimental data, and compare with literature data on gas hold‐ups at the centre of the column. Copyright © 2010 Society of Chemical Industry     

Hydrodynamics of a Tapered Bubble Column

The hydrodynamic behavior of a single‐stage tapered bubble column using an air‐water two‐phase system has been studied. The experimental results indicate that the holdup increases with increasing superficial gas velocity and bubble slip velocity, while it remains constant with increasing superficial liquid velocity. The gas flow rate has a subtle effect on pressure drop owing to the dynamic pressure recovery stemming from the increase in flow area in the axial direction. The results further suggest that the tapered bubble column shows higher holdup with lower energy dissipation than the conventional bubble columns under similar hydrodynamic conditions. The experimental values of the holdup are in excellent agreement with the well‐known Akita and Yoshida correlation available in the existing literature. Also, the performance of the tapered system has been shown to be much better than that of conventional columns under similar conditions in water/alkaline scrubbing of fly ash and SO₂ either alone or in combination.     

A structured catalytic bubble column reactor: hydrodynamics and mixing studies

This paper reports the results of a comprehensive experimental study of the hydrodynamics and mixing in two bubble column reactors of 0.1 and 0.24 m in diameter with KATAPAK-S^® as packing material. Total gas hold up and axial dispersion coefficients were measured in the structured bubble columns and the values were compared with experimental results obtained in the same work with empty bubble columns. The results reveal that the gas hold up in structured bubble columns is practically the same as in empty bubble columns when compared at the same superficial gas velocity based on open area available for gas–liquid dispersion. The presence of the structured elements in the bubble column reactor reduces the liquid phase backmixing by one order of magnitude.     

Turbulent circulation in bubble columns from eddy viscosity distributions of single-phase pipe flow

The exact prediction of the flow structure in bubble columns is important for their design and scale-up. This paper proposes a theoretical model for the liquid circulation on the basis of Reynolds equation of motion and the eddy viscosity distribution of single phase pipe flow, and presents a derivation of an analytical equation for the axial liquid velocity profiles that is fast and easy to use. The model shows a strong analogy for the eddy viscosity between multiphase and single-phase systems, and how both eddy and molecular viscosities affect the flow. The velocity profiles calculated from the model are shown to agree well with reported experimental data for both low and high viscosity fluids.     

Axial dispersion in a rectangular bubble column

This paper presents the results of an experimental study on the gas holdup and the liquid phase axial dispersion coefficient in a narrow packed and unpacked rectangular bubble column. In both cases the gas and liquid flow rates were varied and the data were obtained by employing standard tracer technique. The gas holdup and the axial dispersion coefficient for both the packed and unpacked columns were found to be dependent on the gas and liquid flow rates. For given gas and liquid velocities and a given packing size in the case of the packed column, the rectangular column gave significantly higher dispersion coefficients than a cylindrical column of the equivalent cross sectional area. This result agrees very well with the one predicted by the velocity distribution model. The correlations for the Peclet number, the axial dispersion coefficient, and the fluid holdup for both the unpacked and packed bubble columns are presented.     

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.     

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.     

基于气泡群相间作用力模型的加压鼓泡塔流体力学模拟

目前,多数文献报道了冷态加压湍动鼓泡塔内流动特征,并且通过实验数据回归相关经验关联式。然而,此类关联式适用范围有限,难以直接外推到工业鼓泡塔反应器条件。因此,在FLUENT平台上建立了基于气泡群相间作用力的、动态二维加压鼓泡塔计算流体力学模型。通过数值模拟考察了操作压力为0.5~2.0 MPa,表观气速为0.20~0.31 m·s^-1,内径0.3 m鼓泡塔内流场特性参数分布,并且与冷态实验数据进行比较。结果表明,采用修正后的气泡群曳力模型、径向力平衡模型以及壁面润滑力模型描述气泡群相间作用力,能够较为准确地反映平均气含率和气含率径向分布随操作压力和表观气速变化的规律。     

Hydrodynamic and mass transfer characteristics of bubble and packed bubble columns with downcomer

The hydrodynamic and mass transfer characteristics of bubble and packed bubble columns with downcomer were investigated. The contactor consisted of two concentric columns of 0.11 and 0.2 m i.d., with the annulus acting as the downcomer. The packing used in this investigation was standard 16 mm stainless steel Pall rings. The superficial gas and liquid velocities, V_G and V_L, were varied from 0.01 to 0.09 and 1 × 10^?3 to 8.8 × 10^?3 m s^?1 respectively. Two flow patterns, namely the bubble and pulse flows were observed in the packed bubble column with downcomer, as shown by a flow map. The liquid circulation velocity in both the contactors was observed to be constant throughout the ranges of V_G and V_L covered in this work. The effect of liquid viscosity (0.8 to 9.5 mPa ? s) and surface tension (45 to 72 mN m^?1) on the flow pattern, liquid circulation, gas hold-up and pressure drop was investigated. The pressure drop characteristics across the two contactors have been compared with those across a bubble column. Values of the effective interfacial area, a, and the volumetric mass transfer coefficient, k_L a, were measured by using chemical methods. Values of a as high as 180 and 700 m^?1 and k_L a as high as 0.075 and 0.22 s^?1, in the bubble and packed bubble columns with downcomer, respectively, were obtained. The values of true liquid-side mass transfer coefficient, k_L, were found to be independent of V_G and were of the order of 5.5 × 10^?4 and 3.5 × 10^?4 m s^?1, respectively, in the two contactors.     

多级鼓泡塔内液体速度分布的实验研究

液体循环流动是多级鼓泡塔重要流体力学特征之一,文中在内径为282 mm,高2000 mm的鼓泡塔内,采用不同类型的筛板将普通鼓泡塔分割成双级气液鼓泡塔.采用Pavlov管测液速的方法考察了不同筛板、不同表观气速下该鼓泡塔中上下二侧的液体速度分布.根据实验结果得出了液体速度在塔中心处最大,且与表观气速有关,随着表观气速的...