Hydrodynamic Modelling of Internal Loop Airlift Reactor Applying Drift‐Flux Model in Bubbly Flow Regime

A new model for the liquid circulation rates in airlift reactor (ALR) is presented. The model is based on the energy balance for the flow loop (riser, turn riser‐downcomer, downcomer, and turn downcomer‐riser) coupled with a drift flux theory of two‐phase flow gas‐liquid system, considering a bubbly flow regime. The predicted values of the liquid circulation rates by the developed model are compared with experimental results performed in a 22 dm³ internal loop airlift reactor and with the results obtained in the literatures. The proposed model predicted the experimental results very well. Slip velocity relationship based on the drift flux model was proposed; including the gas holdup, bubble size and the liquid physical properties. The predicted slip velocity was similar to that obtained from the literature. The study revealed that appropriate arrangements of internal bioreactor parts can positively influence the liquid circulation velocity at the same energy consumption. The proposed models are useful in the design; scale up and characterization of the internal loop airlift reactors, and provides a direct method of predicting hydrodynamic behaviour in gas‐liquid airlift reactors.     

3‐D Simulations of an Internal Airlift Loop Reactor using a Steady Two‐Fluid Model

Three‐dimensional (3‐D) simulations of an internal airlift loop reactor in a cylindrical reference frame are presented, which are based on a two‐fluid model with a revised k‐? turbulence model for two‐phase bubbly flow. A steady state formulation is used with the purpose of time saving for cases with superficial gas velocity values as high as 0.12 m/s. Special 3‐D treatment of the boundary conditions at the axis is undertaken to allow asymmetric gas‐liquid flow. The simulation results are compared to the experimental data on average gas holdup, average liquid velocity in the riser and the downcomer, and good agreement is observed. The turbulent dispersion in the present two‐fluid model has a strong effect on the gas holdup distribution and wall‐peaking behavior is predicted. The CFD code developed has the potential to be applied as a tool for scaling up loop reactors.     

CFD simulation and experimental measurement of gas holdup and liquid interstitial velocity in internal loop airlift reactor

This paper documents experiments and CFD simulations of the hydrodynamics of our two-phase (water, air) laboratory internal loop airlift reactor (40 l). The experiments and simulations were aimed at obtaining global flow characteristics (gas holdup and liquid interstitial velocity in the riser and in the downcomer) in our particular airlift configurations. The experiments and simulations were done for three different riser tubes with variable length and diameter. Gas (air) superficial velocities in riser were in range from 1 to 7.5 cm/s. Up to three circulation regimes were experimentally observed (no bubbles in downcomer, bubbles in downcomer but not circulating, and finally the circulating regime). The primary goal was to test our CFD simulation setup using only standard closures for interphase forces and turbulence, and assuming constant bubble size is able to capture global characteristics of the flow for our experimental airlift configurations for the three circulation regimes, and if the simulation setup could be later used for obtaining the global characteristic for modified geometries of our original airlift design or for different fluids. The CFD simulations were done in commercial code Fluent 6.3 using algebraic slip mixture multiphase model. The secondary goal was to test the sensitivity of the simulation results to different closures for the drag coefficient and the resulting bubble slip velocity and also for the turbulence. In addition to the simulations done in Fluent, simulation results using different code (CFX 12.1) and different model (full Euler–Euler) are also presented in this paper. The experimental measurements of liquid interstitial velocity in the riser and in the downcomer were done by evaluating the response to the injection of a sulphuric acid solution measured with pH probes. The gas holdup in the riser and downcomer was measured with the U-tube manometer. The results showed that the simulation setup works quite well when there are no bubbles present in the downcomer, and that the sensitivity to the drag closure is rather low in this case. The agreement was getting worse with the increase of gas holdup in the downcomer. The use of different multiphase model in the different code (CFX) gave almost the same results as the Fluent simulations.     

Experimental Study and CFD Simulation of Hydrodynamic Behaviours in an External Loop Airlift Slurry Reactor

The local hydrodynamic behaviours in an external loop airlift slurry reactor, including the gas holdup, bubble rise velocity, bubble size, were measured with a fibre optic probe. The liquid circulation velocity was measured with an ultrasound Doppler velocimetry. Two‐dimensional simulations were carried out in the framework of Two‐Fluid formulation coupled with a k‐? turbulence model. The lateral forces and interphase turbulence were taken into account and good agreement between the experimental and simulation results was obtained. The simulations show that the lateral forces and interphase turbulence have noticeable influence and should be included in the CFD model.     

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     

Motion of a magnetic flow follower in two‐phase flow — application to the study of airlift reactor hydrodynamics

A low‐cost and simple magnetic particle tracer method was adapted to characterize the hydrodynamic behavior of an internal‐ and an external‐loop airlift reactor (ALR). The residence time distribution of three magnetic particles differing in diameter (5.5, 11.0 and 21.2 mm) and with a density very close to that of water was measured in individual reactor sections. The measured data were analyzed and used to determine the velocity of the liquid phase. Validation of the experimental results for liquid velocity was done by means of the data obtained by an independent reference method. Furthermore, analysis of the differences found in the settling velocity of the particle in single‐liquid and gas‐liquid phases was carried out, using a simplified 3D momentum transfer model. The model considering particle‐bubble interaction forces resulting from changes in the liquid velocity field due to bubble motion was able to predict satisfactorily the increase in the particle settling velocity in the homogeneous bubbly regime. The effective drag coefficient in two‐phase flow was found to be directly dependent on particle Reynolds number to the power of ? 2 but independent of gas flow‐rate for all particle diameters studied. Based on the experimental and theoretical investigations, the valid exact formulation of the effective buoyancy force necessary for the calculation of the correct particle settling velocity in two‐phase flow was done. In addition, recommendations concerning the use of flow‐following particles in internal‐loop ALRs for liquid velocity measurements are presented. Copyright © 2006 Society of Chemical Industry     

Gas holdup,bubble size distribution,and mass transfer in an airlift reactor with ceramic membrane and perforated plate distributor

The airlift reactor is one of the most commonly used gas–liquid two-phase reactors in chemical and biological processes. The objective of this study is to generate different-sized bubbles in an internal loop airlift reactor and characterize the behaviours of the bubbly flows. The bubble size, gas holdup, liquid circulation velocity, and the volumetric mass transfer coefficient of gas–liquid two-phase co-current flow in an internal loop airlift reactor equipped with a ceramic membrane module (CMM) and a perforated-plate distributor (PPD) are measured. Experimental results show that CMM can generate small bubbles with Sauter mean diameter d₃₂ less than 2.5 mm. As the liquid inlet velocity increases, the bubble size decreases and the gas holdup increases. In contrast, PPD can generate large bubbles with 4 mm < d₃₂ < 10 mm. The bubble size and liquid circulation velocity increase as the superficial gas velocity increases. Multiscale bubbles with 0.5 mm < d₃₂ < 10 mm can be generated by the CMM and PPD together. The volumetric mass transfer coefficient k_La of the multiscale bubbles is 0.033–0.062 s⁻¹, while that of small bubbles is 0.011–0.057 s⁻¹. Under the same flow rate of oxygen, the k_La of the multiscale bubbles increases by up to 160% in comparison to that of the small bubbles. Finally, empirical correlations for k_La are obtained.     

Effect of Surface Tension on Hydrodynamics and Mass Transfer Coefficient in Airlift Reactors

The hydrodynamics and volumetric mass transfer coefficient play important roles in the design and scale-up of airlift reactors. The effect of surface tension on hydrodynamics and volumetric mass transfer coefficient in internal loop airlift reactors was investigated. With reduction of the surface tension of the fluid, the hydrodynamic parameters raised, namely, gas phase holdup, flow regime transition point, and interfacial area, whereas the bubble diameter as well as the liquid velocity decreased and the volumetric mass transfer coefficient increased. Empirical correlations are proposed for gas-phase holdup and volumetric mass transfer coefficient in terms of dimensionless numbers and can be applied in the design of airlift reactors.     

Numerical simulation of the dynamic flow behavior in a bubble column: A study of closures for turbulence and interface forces

Numerical simulations of the bubbly flow in two square cross-sectioned bubble columns were conducted with the commercial CFD package CFX-4.4. The effect of the model constant used in the sub-grid scale (SGS) model, C_S, as well as the interfacial closures for the drag, lift and virtual mass forces were investigated. Furthermore, the performance of three models [Pfleger, D., Becker, S., 2001. Modeling and simulation of the dynamic flow behavior in a bubble column. Chemical Engineering Science, 56, 1737-1747; Sato, Y., Sekoguchi, K.,1975. Liquid velocity distribution in two-phase bubble flow. International Journal of Multiphase Flow 2, 79-95; Troshko, A.A., Hassan, Y.A., 2001. A two-equation turbulence model of turbulent bubbly flows. International Journal of Multiphase Flow 27, 1965-2000] to account for the bubble-induced turbulence in the k-ε model was assessed. All simulation results were compared with experimental data for the mean and fluctuating liquid and gas velocities. It is shown that the simulation results with C_S=0.08 and 0.10 agree well with the measurements. When C_S is increased, the effective viscosity increases and subsequently the bubble plume becomes less dynamic. All three bubble-induced turbulence models could produce good solutions for the time-averaged velocity. The models of Troshko and Hassan and Pfleger and Becker reproduce the dynamics of the bubbly flow in a more accurate way than the model of Sato and Sekoguchi. Based on the comparison of the results obtained for two columns with different aspect ratio (H/D=3 and H/D=6), it was found that the model of Pfleger and Becker performs better than the model of Troshko and Hassan, while the model of Sato and Sekoguchi performs the worst. It was observed that the interfacial closure model proposed by Tomiyama [2004. Drag, lift and virtual mass forces acting on a single bubble. Third International Symposium on Two-Phase Flow Modeling and Experimentation, Pisa, Italy, 22-24 September] performs better for the taller column. With the drag coefficient proposed by Tomiyama, the predicted slip velocity agrees well with the experimental data in both columns. The virtual mass force has a small influence on the investigated bubbly flow characteristics. However, the lift force strongly influences the bubble plume dynamics and consequently determines the shape of the vertical velocity profile. In a taller column, the lift coefficient following from the model of Tomiyama produces the best results.     

气升式环流反应器的理论研究进展

概述了气升式环流反应器的流体力学模型,对工程实用性较强的气液两流体模型和气泡导致的湍动模型进行了分类,并对常见相间作用力模型和多相流求解技术进行了总结和评价. 指出将Favre平均的两流体模型和采用Sato模型考虑气泡导致湍动的k-e模型相结合是现阶段气升式环流反应器设计和放大的有力工具. 在分析了气液两相流理论现状的基础上,指出了其存在的问题和进一步的研究方向.     

Hydrodynamics of a Multi‐Stage Internal Loop Airlift Reactor

The local hydrodynamic properties in a multi‐stage internal loop airlift reactor were investigated in this study. The gas‐liquid two‐phase flow hydrodynamic properties, including gas holdup, bubble velocity, bubble diameter, and liquid circulation velocity at various stages were measured by dual electrical resistivity probes and conductivity cells. Detailed studies on the gas holdup, bubble velocity, bubble diameter, and liquid circulation velocity were conducted with respect to various values of superficial gas. The Zuber and Findlay drift flux model was used to represent the variation of slip velocity with total gas‐liquid velocity at various stages and the model fits the data well.     

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.     

Solid effects on hydrodynamics and heat transfer in an external loop airlift reactor

The effect of a solid presence on global hydrodynamic parameters and heat transfer in an external loop airlift reactor has been experimentally investigated. Results obtained in both two- and three-phase flow are presented in this study. Two different external loop airlift reactor sizes have been used and local hydrodynamic characteristics including local gas hold-up and bubble velocity have been obtained in two-phase flow. Optical and ultrasound probes have been used to obtain this information, respectively. It was found that an increase of solid hold-up leads to a decrease of liquid velocity and heat transfer coefficient. Measured in a two- and three-phase reactor using a horizontal-heating probe, a correlation of the average gas hold-up and heat transfer coefficient is proposed. Correlation parameters are identified in homogeneous and heterogeneous flow regimes, which have been derived from the gas slip velocity concept. The experimental liquid velocity and gas hold-up in the riser have been represented in a satisfactory way by a hydrodynamic model, either in the absence or in the presence of solid particles.     

A modified pseudopotential for a lattice Boltzmann simulation of bubbly flow

The pseudopotential in Shan and Chen-type multiphase models was investigated and modified based on a virial equation of state with newly proposed parameters. This modified pseudopotential was used in a lattice Boltzmann model and shown to be suitable for simulating sufficiently large gas–liquid density ratios with good numerical stability and only small spurious velocities. The spurious velocity was reduced by reducing the pseudo-sound speed by the use of suitable parameters. The multicomponent multiphase model based on this modified pseudopotential can be used in bubbly flow simulations. Bubble rise behavior was simulated using a 3D multicomponent and multiphase model with a high density ratio. The predicted terminal velocity and drag coefficient of a single bubble agreed well with those calculated from empirical correlations. The drag coefficient of bubbles in the homogenous regime decreased with increased gas holdup. A new relationship between the bubble drag coefficient and gas holdup in the homogenous regime was proposed.     

Gas hold-up and liquid mixing at low and intermediate gas velocities I. Air-water system

Gas hold-up and liquid phase dispersion experiments have been carried out in a 0.06 m bubble column at varying liquid and gas velocities. The results obtained show that the coefficient of liquid mixing varies with the flow regime. The isotropic turbulence theory of Baird and Rice (Chem. Eng. J., 9 (1975) 171) was used to provide dimensionally consistent correlations for the chain bubbling, bubbly and churn turbulent flow regimes. The gas hold-up was determined to increase with gas velocity in the chain bubbling and bubbly flow regimes. The results obtained from this study also show that the Froude number represents a useful criterion for mapping flow regimes in vertical bubble columns.     

用CFD研究气升式内环流生物反应器下降管中的流体力学性质

在对气升式内环流生物反应器内部流动分析基础上,全面考虑反应器下降管中气泡的并聚破碎、气液两相间相互作用和滑移等, 建立了能描述反应器下降管中复杂流动的CFD数学模型. 运用CFX-4.4对模型方程进行求解, 通过求解得到了包括气液两相速度场、局部气含率分布等详尽信息,并就液相流动速度与相应条件下的PIV测试结果进行了比较,主体流动速度的偏差在20%以下,且两者总的变化趋势一致.该模型能较好地预测反应器下降管内的复杂流场.     

Advanced PIV/LIF and shadowgraphy system to visualize flow structure in two-phase bubbly flows

Particle image velocimetry (PIV) is a promising technique to measure dispersed phase size, dispersed phase hold-up and velocity of both the phases. The current work reports measurement of the shape, size, velocity and acceleration of bubbles using shadowgraphy, and liquid velocity measurement obtained using PIV/LIF with fluorescent tracer particles. Measurements were performed in a narrow rectangular column at moderate gas hold-up (～5%) with wide variation of bubble sizes (0.1–15 mm). The liquid velocity field was subjected to 2D discrete wavelet transform (DWT) to visualize the flow structures in the bubbly flow. Further, the slip velocity of individual bubbles was obtained from the DWT filtered liquid velocity field. The results are compared with the slip velocity correlations reported in literature for single bubbles rising in quiescent water. The comparison shows the difference in slip velocity of single bubbles and bubbles rising in swarm. The scale wise decomposition obtained from DWT was also used to quantify the liquid velocity field in terms of wavenumber spectrum. The velocity and acceleration measurements are demonstrated on a single spherical cap bubble rising in quiescent water. The measurements show the potential of the 2D acceleration measurement to facilitate the estimation of unsteady drag on bubbles.     

Estimation of kLa Values in Bench‐Scale Stirred Tank Reactors with Self‐Inducing Impeller by Multiphase CFD Simulations

A multiphase computational fluid dynamics (CFD) simulation methodology is developed and proposed for the estimation of the spatial distribution of k_La values in a bench‐scale reactor equipped with a self‐inducing impeller. The importance of estimating an apparent drag coefficient, which considers the effect of turbulence on the gas bubble rising velocity, is also tackled by applying different correlations available in literature, namely, Brucato, modified Brucato, and Pinelli correlations. The spatial distribution of k_La values in the agitated vessel is found from the CFD results using Danckwert's surface renewal model. An analysis of the gas volume fraction distribution obtained from the simulations is performed in order to choose the most suitable drag model. The modified Brucato correction correlation for the drag force exhibits the best agreement with experimental data.     

二维鼓泡床内气液流动特性实验与数值模拟

采用高速摄像法测量了0.20 m×0.02 m×2.00 m拟二维床内气泡尺寸分布和流型等变化规律,结果表明,随着表观气速的增大,鼓泡床内依次呈现均匀鼓泡区、过渡区和湍动区3种形式,以气泡个数概率表示的气泡尺寸分布呈对数正态分布。以计算流体力学软件ANSYS CFX 10.0为平台,采用k-ε湍流模型和GRACE曳力模型对气液鼓泡床内流体动力学行为展开了数值模拟,其结果与实验值比较吻合。研究表明,从多相流理论出发的计算流体力学模拟方法可以用来预报鼓泡床内流型过渡等流体动力学特性。