Numerical simulation of multiphase flow in bubble column reactors. Influence of bubble coalescence and break-up

Population balance equations have been combined to a classical hydrodynamic Euler/Euler simulation to investigate the operation of a cylindrical bubble column. The MUSIG (mutiple-size-group) model implemented in the CFX 4.3 commercial software has been used. Hydrodynamic experimental variables, i.e. local axial liquid velocity and local gas hold-up, have been compared to the corresponding calculated values, showing a quite good agreement, except for the gas hold-up when the column is no more operating in the homogeneous regime. Bubble sizes have been investigated, showing that two domains of superficial gas velocities can be distinguished. In the first domain, coalescence occurs predominantly, Sauter diameter increases with the superficial gas velocity, bubble size distribution is narrow and Sauter diameter is continuously evolving along the column axis. In the second domain, break-up becomes more intensive and compensates coalescence, bubble size distribution becomes wider, since more small bubbles are formed, an equilibrium Sauter diameter appears when the superficial gas velocity increases. Furthermore an equilibrium Sauter diameter appears along the column axis, and it can be noticed that this phenomenon appears lower in the column when the gas flow rate is increased. In these two domains the characteristics of the bubbles are typical of those of the homogeneous and transition regimes.     

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.     

Theoretical prediction of flow regime transition in bubble columns by the population balance model

Theoretical prediction of flow regime transition in bubble columns was studied based on the bubble size distribution by the population balance model (PBM). Models for bubble coalescence and breakup due to different mechanisms, including coalescence due to turbulent eddies, coalescence due to different bubble rise velocities, coalescence due to bubble wake entrainment, breakup due to eddy collision and breakup due to large bubble instability, were proposed. Simulation results showed that at relatively low superficial gas velocities, bubble coalescence and breakup were relatively weak and the bubble size was small and had a narrow distribution; with an increase in the superficial gas velocity, large bubbles began to form due to bubble coalescence, resulting in a much wider bubble size distribution. The regime transition was predicted to occur when the volume fraction of small bubbles sharply decreased. The predicted transition superficial gas velocity was about 4 cm/s for the air-water system, in accordance with the values obtained from experimental approaches.     

CFD simulation of gas-liquid contacting in tubular reactors

Gas-liquid contacting in tubular reactors was simulated using an Eulerian-Eulerian CFD approach in which accurate interphase momentum closure relations are incorporated, bubble-induced turbulence is accounted for, and population balance equations are used to describe bubble breakage and coalescence. The ability of two breakup kernels (Luo, H., Svendsen, H.F., 1996. Theoretical model for drop and bubble breakup in turbulent dispersions. A.I.Ch.E. Journal 42, 1225-1233; Lehr, F., Millies, M., Mewes, D., 2002. Bubble size distributions and flow fields in bubble columns. A.I.Ch.E. Journal 48, 2426-2443) and three coalescence kernels (Prince, M.J., Blanch, H.W., 1990. Bubble coalescence and breakup in air sparged bubble columns. A.I.Ch.E. Journal 36, 1485-1499; Luo, H., 1993. Coalescence, breakup and liquid recirculation in bubble column reactors. Ph.D. Thesis, Norwegian University of Science and Technology, Trondheim; Lehr, F., Millies, M., Mewes, D., 2002. Bubble size distributions and flow fields in bubble columns. A.I.Ch.E. Journal 48, 2426-2443) to accurately predict several flow parameters in pipe flow was tested.Good agreement between simulation and experimental results (radial profiles of gas holdup, turbulence intensity, and local Sauter bubble diameter) was achieved without the use of empirically derived relationships (such as Drift flux) by adjusting a single parameter which accounts for the deviation in the coalescence behaviour of tap water from that of pure water. The approach adopted in this investigation may thus be applicable to more complex hydrodynamic situations such as those encountered in mechanically agitated tanks and the need for extensive experimental testing may be replaced by single measurement of the effect interfacial properties have on coalescence rates.     

Experimental and numerical investigations of scale-up effects on the hydrodynamics of slurry bubble columns

Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid concentrations(0–30 vol%). Radial profiles of time-averaged gas holdup, axial liquid velocity, and turbulent kinetic energy were measured by using in-house developed conductivity probes and Pavlov tubes. Effects of column diameter, superficial gas velocity, and solid concentration were investigated in a wide range of operating conditions. Experimental results indicated that the average gas holdup remarkably increases with superficial gas velocity, and the radial profiles of investigated flow properties become steeper at high superficial gas velocities. The axial liquid velocities significantly increase with the growth of the column size, whereas the gas holdup was slightly affected. The presence of solid in bubble columns would inhibit the breakage of bubbles, which results in an increase in bubble rise velocity and a decrease in gas holdup, but time-averaged axial liquid velocities remain almost the same as that of the hollow column. Furthermore, a 2-D axisymmetric k–ε model was used to simulate heterogeneous bubbly flow using commercial code FLUENT 6.2. The lateral lift force and the turbulent diffusion force were introduced for the determination of gas holdup profiles and the effects of solid concentration were considered as the variation of average bubble diameter in the model. Results predicted by the CFD simulation showed good agreement with experimental data.     

Multicompartment hydrodynamic model for slurry bubble columns

A core-annulus multicompartment two-dimensional two-bubble class model accounting for slurry recirculation and coupled with catalyst transport was developed as a part and parcel of the analysis of the behavior of slurry bubble column reactors at high gas throughputs corresponding to the churn turbulent flow regime. The model analyzed the contributions of bubble-induced turbulence closures, bubble coalescence and breakup phenomena, and catalyst axial distribution as the resultant of sedimentation, advection via liquid-solid slip, per-compartment axial dispersion and core-annulus lateral exchange of catalyst by bubble-induced turbulence. The model was also used to analyze the effects of catalyst loading, gas density and superficial velocity, and column diameter and vessel aspect ratio on the hydrodynamics of slurry bubble column reactors, namely, the per-compartment phase holdups and interstitial velocities, pressure gradient, bubble coalescence and break-up rates, and loci of velocity inversion for the gas and slurry profiles.     

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.     

Interaction between partitioning porous plate and rising bubbles in a trayed bubble column

In a trayed bubble column, the structure of the partitioning plate plays an important role on the bubble behavior. This study examined the effect of the opening ratio and pore size of the plate on the bubble break-up frequency and bubble size distribution. The sieve tray was used as the partitioning plate. The opening ratio was closely related to gas cap development. The stagnation of bubble flow and a gas cap were observed with an opening ratio less than 48.5%. The gas cap increased with decreasing opening ratio and increasing superficial gas velocity. The main effect of the sieve tray could be categorized into the additional drag force and bubble break-up depending on the sieve pore size. When the sieve pore size was smaller than the Sauter diameter of the bubble swarm, the movement of rising bubbles was interrupted by the drag force applied by the surrounding mesh lines. On the other hand, when the sieve pore size was larger than the Sauter diameter, the bubbles were affected by the additional bubble break-up. After the bubbles penetrated the sieve tray, the bubble size distribution shifted to a smaller one and the Sauter diameter decreased.     

Effect of antifoam agents on bubble characteristics in bubble columns based on acoustic sound measurements

In this paper, the effect of antifoam agents on bubble characteristics in bubble columns is studied. Specifically, the bubble characteristics of air in tap water are compared to those of air in 5% and 10% antifoam solutions. Bubble characteristics such as gas holdup, bubble diameter, bubble-size distribution, and damping ratio were investigated at various superficial gas velocities. These properties were deduced from the acoustic sound measurement. The study revealed that the addition of antifoam chemicals reduces the overall gas holdup and increases the average bubble diameter. The bubble-size distribution in tap water is found to be homogeneous while in antifoam solutions to be heterogeneous. It is also found that at low gas velocities the damping ratio for antifoam solutions is higher than that for tap water, while at high gas velocities the damping ratio is not affected. The results affirm that acoustic probes are excellent measuring tools over classical tools at moderate gas velocities.     

Numerical simulation of bubble columns flows: effect of different breakup and coalescence closures

Two-dimensional axisymmetric Eulerian/Eulerian simulations of two-phase (gas/liquid) transient flow were performed using a multiphase flow algorithm based on the finite-volume method. These numerical simulations cover laboratory scale bubble columns of different diameters, operated over a range of superficial gas velocities ranging from the bubbly to the churn turbulent regime. The bubble population balance equation (BPBE) is implemented in the two-fluid model that accounts for the drag force and employs the modified k-ε turbulence model in the liquid phase. Several available bubble breakup and coalescence closures are tested. Quantitative agreements between the experimental data and simulations are obtained for the time-averaged axial liquid velocity profiles, as well as for the kinetic energy profiles, only when model predicted breakup rate is increased by a factor of ten to match the coalescence rate. The calculated time-averaged gas holdup profiles deviate in shape from the measured ones and suggest that full three-dimensional simulation is needed. Implementation of BPBE leads to better agreement with data, especially in the churn-turbulent flow regime, compared to the simulation based on an estimated constant mean bubble diameter. Differences in the predicted interfacial area density, with and without BPBE implementation, are significant. The choice of bubble breakup and coalescence closure does not have a significant impact on the simulated results as long as the magnitude of breakup is increased tenfold.     

Numerical simulation of two-phase flow with bubble break-up and coalescence coupled with population balance modeling

A computational fluid dynamic (CFD) model was developed with an improved source term based on previous work by Hagesaether et al. [1] for bubble break up and bubble coalescence to carry out numerical prediction of number density of different bubble class in turbulent dispersed flow. The numerical prediction was based on two fluid models, using the Eulerian–Eulerian approach where the liquid phase was treated as a continuum and the gas phase (bubbles) was considered as a dispersed phase. Bubble–bubble interactions, such as breakage due to turbulence and coalescence due to the combined effect of turbulence and laminar shear were considered. The result shows that the radial distributions of number densities of lower bubble classes are more than its higher counterpart. The result also shows that the Sauter mean diameter increases with the increase of height up to 1 m and then become steady. Simulated results are found to be in good agreement with the experimental data.     

CHARACTERIZATION OF THE DOWNFLOW SECTION OF AN AIRLIFT COLUMN USING BUBBLE SIZE DISTRIBUTION MEASUREMENTS

Bubble size distributions in an airlift column were investigated with an emphasis on the downflow section. Measurements have been made using direct photographic techniques in conjunction with image analysis in a split cylinder airlift column. Information extracted from these measurements includes local gas hold-up, variation of Sauter mean bubble diameter with column length, and liquid circulation velocity. An air-tap water system was studied for purposes of comparison, while effects of electrolyte concentration and viscosity were studied using salt water, and two carboxymethyl cellulose (CMC) concentrations in water, respectively. The effect of energy input was studied by varying air flow rate to produce superficial velocities ranging from 2·59 cm/s to 10·36 cm/s.     

Radial Profiles of Gas Bubble Behavior in a Gas‐Liquid Bubble Column Reactor under Elevated Pressures

The effects of operating conditions on radial variation of gas holdups, bubble swarm rising velocity, and Sauter mean diameter were investigated in a bubble column reactor under elevated pressures using a conductivity probe method. Air served as gas phase and tap water as liquid phase with varying gas velocity and pressure. All three parameters increased constantly with higher superficial gas velocity. Maximum holdup, velocity, and Sauter mean diameter were found at the center of the cross section. Two different cases for Sauter mean diameter distribution were observed. The gas holdups increase continuously with higher system pressure, but decrease for bubble swarm rising velocity and Sauter mean diameter. According to experimental results, an empirical correlation of the gas holdup profiles is presented.     

Effect of liquid properties on the performance of bubble column reactors with fine pore spargers

This work is a study of the effect of liquid properties on the performance of bubble column reactors with fine pore spargers. Various liquids covering a range of surface tension and viscosity values are employed, while the gas phase is atmospheric air. A fast video technique is used for visual observations and, combined with image processing, is used for gas holdup and bubble size measurements. New data on average gas holdup values, bubble size distributions and Sauter diameters are presented and are consistent with existing physical models on coalescence/breakage. A correlation based on dimensionless groups for the prediction of gas holdup in the homogeneous regime is proposed and found to be in good agreement with available data.     

Flow regimes and mass transfer in counter-current bubble columns

Counter current bubble columns have the feature that specific gas-liquid interfacial area and gas holdup are larger than those for standard and cocurrent bubble columns. In this study, three different flow regimes, churn-turbulent flow, bubble flow and bubble down-flow, have been observed in a counter-current bubble column and correlations of gas holdup and volumetric liquid-phase mass transfer coefficient have been proposed as functions of operating variables such as the superficial velocities of gas and liquid, the gas-liquid slip velocity and the liquid properties.     

Small bubble formation via a coalescence dependent break-up mechanism

A mechanism has been elucidated for the coalescence-mediated break-up of bubbles in gas-liquid systems. Images taken from dynamic systems (a coalescence cell and laboratory-scale bubble columns) show that in some instances the coalescence of two bubbles is accompanied by the formation of a much smaller daughter bubble. Following the coalescence process an annular wave is formed due to the very rapid expansion of the hole following the instant of film rupture. As the wave moves along the length of the bubble, away from the point of rupture it causes a rippling effect which distorts the newly coalesced bubble and may result in the formation of an unstable extension. Instabilities due to the annular wave pinch off a portion of this extension, resulting in the generation of a small daughter bubble. In coalescence dominated systems the process results in the generation of significant numbers of bubbles much smaller (100- diameter) than the Sauter mean diameter (3-).     

Effect of gas distributor on hydrodynamics in shallow bubble column reactors

The effects of gas distributor on hydrodynamics in an air–water shallow bubble column reactor are investigated. Three types of distributors, namely, single nozzle, perforated plate and porous plate are being studied. The overall gas holdup, bubble size distributions and bubble rise velocity distributions are studied over a range of superficial gas velocities. The results show that single nozzle is not suitable for shallow bed operation. While perforated plate and porous plate distributors have comparable behaviour in the absence of solids, the addition of solids particles causes the two distributors to behaviour differently. The presence of solids promotes bubble coalescence for perforated plate distributor while the same inhibits bubble coalescence for porous plate distributor.     

Numerical simulations of gas-liquid mass transfer in bubble columns with a CFD-PBM coupled model

Gas-liquid mass transfer in a bubble column in both the homogeneous and heterogeneous flow regimes was studied by numerical simulations with a CFD-PBM (computation fluid dynamics-population balance model) coupled model and a gas-liquid mass transfer model. In the CFD-PBM coupled model, the gas-liquid interfacial area a is calculated from the gas holdup and bubble size distribution. In this work, multiple mechanisms for bubble coalescence, including coalescence due to turbulent eddies, different bubble rise velocities and bubble wake entrainment, and for bubble breakup due to eddy collision and instability of large bubbles were considered. Previous studies show that these considerations are crucial for proper predictions of both the homogenous and the heterogeneous flow regimes. Many parameters may affect the mass transfer coefficient, including the bubble size distribution, bubble slip velocity, turbulent energy dissipation rate and bubble coalescence and breakup. These complex factors were quantitatively counted in the CFD-PBM coupled model. For the mass transfer coefficient k_l, two typical models were compared, namely the eddy cell model in which k_l depends on the turbulent energy dissipation rate, and the slip penetration model in which k_l depends on the bubble size and bubble slip velocity. Reasonable predictions of k_la were obtained with both models in a wide range of superficial gas velocity, with only a slight modification of the model constants. The simulation results show that CFD-PBM coupled model is an efficient method for predicting the hydrodynamics, bubble size distribution, interfacial area and gas-liquid mass transfer rate in a bubble column.     

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.     

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.