Compressibility induced bubble size variation in bubble column reactors: Simulations by the CFD–PBE

Bubble column reactors can be simulated by the two fluid model (TFM) coupled with the population balance equation (PBE). For the large industrial bubble columns, the compressibility due to the pressure difference may introduce notable bubble size variation. In order to address the compressibility effect, the PBE should be reformulated and coupled with the compressible TFM. In this work, the PBE with a compressibility term was formulated from single bubble dynamics, the mean Sauter diameters predicted by the compressible TFM coupled with the PBE were compared with the analytical solutions obtained by the ideal gas law. It was proven that the mesoscale formulations presented in this work were physically consistent with the macroscale modeling. It can be used to simulate large industrial plants when the compressibility induced bubble size variation is important.     

CFD simulation of gas–liquid flow in a high-pressure bubble column with a modified population balance model

In this study, based on the Luo bubble coalescence model, a model correction factor C_e for pressures according to the literature experimental results was introduced in the bubble coalescence efficiency term. Then, a coupled modified population balance model (PBM) with computational fluid dynamics (CFD) was used to simulate a high-pressure bubble column. The simulation results with and without C_e were compared with the experimental data. The modified CFD-PBM coupled model was used to investigate its applicability to broader experimental conditions. These results showed that the modified CFD-PBM coupled model can predict the hydrodynamic behaviors under various operating conditions.     

气液逆流鼓泡塔中的气含率与液速分布和数值模拟

分别采用线－线电导探针和背靠背式皮托和测定了气液逆流鼓泡塔中局部气含率和液体轴向速度的径向和轴向分布,讨论了操作条件对气含率与液体轴向速度分布的影响,利用双流体模型对塔内两相湍流流动进行了数值模拟。     

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.     

Numerical and technological properties of bubble column bioreactors for aerobic processes

This paper describes the method of modelling and numerical simulation of bubble column bioreactors in which the process of aerobic biodegradation of carbonaceous substrate occurs. Such bioreactors belong to systems with distributed state variables. Determining steady states of such objects results in solving non-linear boundary-value problems.The bioreactors with axial dispersion and piston flow with biomass recirculation were analysed. The effectiveness of selected algorithms used to determine the steady states of such bioreactors were compared; the numerical properties of mathematical models of analysed bioreactors were specified, the branches of steady states in bioreactors with axial dispersion and piston flow were determined and the application areas of such bioreactors were defined.     

Time-space-property least squares spectral method for population balance problems

Population Balance Equations (PBE) describe complex processes where the accurate prediction of the dispersed phase plays a major role for the overall behavior of the system. However, the solution of this type of equation is computationally intensive. In this paper, a least squares spectral method for solving the population balance in d=1,2,3 physical domains, p=1,… property domain and time is discussed. In the presented approach, the solution is found as the function that minimizes the residual function in the least squares sense. The capability of the method is shown by using some model problems related to the advective transport of the density function which describes the dispersed phase.     

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.     

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.     

A numerical method for solving the transient multidimensional population balance equation using an Euler-Lagrange formulation

A numerical method was developed to solve the population balance equation for transient multidimensional problems including particle-particle interactions. The population balance equation was written in a mixed Euler-Lagrange formulation which was solved using the discretization method that represents the number density function by impulse functions, an operator splitting method and a remeshing procedure for the internal variable that conserves the mass and the number of particles.This method was successfully tested against analytical and semi-analytical solutions for pure breakage, pure coalescence, breakage and coalescence, pure advection, advection with absorption, advection with binary uniform breakage and with constant or linear absorption. The method was also applied to a free-boundary transient one-dimensional gas-phase model in a bubble column reactor with simplified hydrodynamics. Accurate solutions were obtained for several simulation conditions for the bubble column, including gas absorption, bubble breakage, bubble coalescence and variable gas density effects. The results showed that the numerical method is adequate and robust for solving transient population balance problems with spatial dependence and particle-particle interactions.     

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.     

Modelling local bubble size distributions in agitated vessels

Photography and capillary suction probe were used to measure local bubble size distributions (BSDs) from Rushton turbine agitated (14/200 L) air-tap water and CO₂-n-butanol dispersions. A multiblock stirred tank model with population balances (PBs) for bubbles was created to describe local BSDs in agitated vessels. Unknown parameters in breakage and coalescence models were adjusted by comparing the predicted and measured local BSDs. The BSDs from both investigated systems and varying vessel-operating conditions were included simultaneously to the fitting. The adjusted models were incorporated to MUSIG PB model in CFX-5.7 and tested for the laboratory stirred tanks. The multiblock model showed to be an optimal trade-off between the accuracy and CPU time for the investigation of gas-liquid hydrodynamics and validation of closure models. As a result of fitting, the adjusted model seems to describe local BSDs more accurately in agitated vessels than the model of Lehr et al. [2002. Bubble-size distributions and flow fields in bubble columns. A.I.Ch.E. Journal 48, 2426-2443], which has been successful in bubble column studies. This shows that phenomenological breakage and coalescence closures need experimental validation for various flow environments.     

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.     

Hydrodynamics of air–kerosene bubble column under elevated pressure in homogeneous flow regime

A multiphase computational fluid dynamics(CFD) model coupled with the population balance equation(PBE) was developed in a homogeneous air–kerosene bubble column under elevated pressure(P). The specific pressure drop(DP/L), gas holdup(a_G), and Sauter mean diameter(d_(32)) were experimentally measured in the bubble column with 1.8 m height and 0.1 m inner diameter, which was operated at a superficial gas velocity of 12.3 mm·s~(-1), and P = 1–35 bar(1 bar = 10~5 Pa). A modified drag coefficient model was proposed to consider the effect of bubble swarm and pressure on hydrodynamics of the bubble column.The Luo breakage model was modified to account for liquid density, viscosity, surface tension and gas density. The DP/L, a_G, and d_(32) obtained from the CFD model were compared with experimental data,and the gas density-dependent parameters of the CFD model were identified. With increasing P from 1 to 35 bar, the aGvaried from 5.4% to 7.2% and the d_(32) decreased from 2.3 to 1.5 mm. The CFD-PBE model is applicable to predict hydrodynamics of pressurized bubble columns for gas–organic liquid in the homogeneous regime.     

Three-dimensional mathematical modeling of dispersed two-phase flow using class method of population balance in bubble columns

Computational fluid dynamics (CFD) simulations of bubble columns have received recently much attention and several multiphase models have been developed, tested, and validated through comparison with experimental data. In this work, we propose a model for two-phase flows at high phase fractions. The inter-phase forces (drag, lift and virtual mass) with different closure terms are used and coupled with a classes method (CM) for population balance. This in order to predict bubble’s size distribution in the column which results of break-up and coalescence of bubbles. Since these mechanisms result greatly of turbulence, a dispersed k– turbulent model is used.The results are compared to experimental data available from the literature using a mean bubble diameter approach and CM approach and the appropriate formulations for inter-phase forces in order to predict the flow are highlighted.The above models are implemented using the open source package OpenFoam.     

Lift force on bubbles in a bubble column reactor: Experimental analysis

The lift force acting on bubbles in a swarm has been estimated by analyzing the instantaneous velocity-time data obtained using LDA in a cylindrical bubble column. Phase distinction was achieved through the multiresolution analysis of the velocity-time data. Several important issues related to the transverse motion of bubbles subjected to a shear field have been discussed quantitatively. The actually measured bubble sizes, the respective slip velocity values in transverse and axial directions and the local shear rates (γ) enabled the verification of known formulations for the lift coefficient (C_L) for bubbles. At many locations in the column the radial flux of the gas phase by turbulent dispersion and the radial slip were estimated. The radially inward movement of bubbles from low to high axial velocity (from column wall to center, i.e., C_L<0) was observed at most of the measurement locations. The local lift coefficient was estimated using the transverse drag force and the values support the results from the material balance approach. The estimated C_L values showed a wide variation over the column cross-section.