Bubble mechanisms and characteristics at pore scale in a packed-bed reactor

An image processing technique was used to study dominant bubble mechanisms in a two-dimensional packed-bed at pore level under the bubbly flow regime. Bubble breakup and coalescence were identified as dominant mechanisms using a large number of image samples. Two types of coalescence mechanisms were identified that occur due to compression and deceleration associated with the bubbles and three breakup mechanisms were identified that are result of liquid shear force, bubble acceleration, and bubble impact. Data on various two-phase parameters, such as local void fraction, bubble velocity, size, number, and shape were obtained from the images. Results indicated that when a flow regime changed from bubbly to either trickling or pulsing flow, the number of average sized bubbles significantly decreased and the shape of the majority of the bubbles was no longer spherical. Although a mean bubble velocity of all sized bubbles was uniform for given gas and liquid superficial velocities, individual bubble velocities were quite different depending on the bubble location in the pore. The present bubble size distributions were compared with previous studies and the results on bubble size are in general agreement.     

Pressure drop in monolith reactors

A theoretical model for the computation of pressure drop in bubble-train flow inside capillaries of square cross-section was developed. The model is based on three contributions: hydrostatics, viscous pressure drop, and capillary pressure drop. Capillary pressure drop is related to the shape of the fronts and ends of the bubbles. The model does not include entrance or exit effects, has no adjustable parameters, and agrees very well with available experimental data.

For a given set of flow parameters, bubble velocity and liquid slug average velocity are computed as a function of gas and liquid superficial velocities. The length of the unit cell determines the number of bubbles inside the capillary for a given flow situation. The model requires experimental information of average bubble lengths to compute the length of a unit cell consisting of a bubble and a liquid slug.

The three pressure contributions for a unit capillary length are linear functions of the number of bubbles inside the capillary. The length of the bubbles in bubble-train flows is a critical parameter in the computation of pressure drop.     

Viscous co‐current downward Taylor flow in a square mini‐channel

This article presents a computational study of the co‐current downward Taylor flow of gas bubbles in a viscous liquid within a square channel of 1 mm hydraulic diameter. The three‐dimensional numerical simulations are performed with an in‐house computer code, which is based on the volume‐of‐fluid method with interface reconstruction. The computed (always axi‐symmetric) bubble shapes are validated by experimental flow visualizations for varying capillary number. The evaluation of the numerical results for a series of simulations reveals the dependence of the bubble diameter and the interfacial area per unit volume on the capillary number. Correlations between bubble velocity and total superficial velocity are also provided. The present results are useful to estimate the values of the bubble diameter, the liquid film thickness and the interfacial area per unit volume from given values of the gas and liquid superficial velocities. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

矩形截面螺旋通道内气液两相流局部含气率分布实验研究

应用电导探针测量技术,对矩形截面螺旋通道内气液两相流局部含气率进行实验研究。在不同的气相折算速度下,应用电导探针测量了弹状流弹单元的长度,并与可视化方法进行对比,验证了电导探针的可靠性,并为信号处理选择合适的阈值。分别在泡状流、弹状流及环状流三种流型的条件下,分析了气相与液相折算速度对局部含气率分布的影响。实验结果发现,螺旋通道气液两相局部含气率呈非对称的抛物线形分布,这种非对称性受流型和液相折算速度的影响。     

Local hydrodynamics modeling of a gas–liquid–solid three-phase bubble column

A three-dimensional (3D) transient model was developed to simulate the local hydrodynamics of a gas–liquid–solid three-phase bubble column using the computational fluid dynamic method, where the multiple size group model was adopted to determine the size distribution of the gas bubbles. Model simulation results, such as the local time-averaged gas holdups and axial liquid velocities, were validated by experimental measurements under varied operating conditions, e.g., superficial gas velocities and initial solid loadings at different locations in the three-phase bubble column. Furthermore, the local transient hydrodynamic characteristics, such as gas holdups, liquid velocities, and solid holdups, as well as gas bubble size distribution were predicted reasonably by the developed model for the dynamic behaviors of the three-phase bubble column. © 2007 American Institute of Chemical Engineers AIChE J, 2007     

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.     

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.     

鼓泡塔内气液两相湍流实验研究

介绍了研究鼓泡塔气液两相流的实验装置、实验方法。液相用激光多普勒测速技术(LDV)测量,气相用粒子示踪测速技术(PIV)测量。实验表明,轴向液相速度的径向分布呈塔中心峰值、壁面附近倒流形式,且与气相表观速度大小有关,当液相表观速度一定时,随气相表观速度增大而愈加陡峭,返混也剧烈。当表观液速与表观气速之比小于19．6时,返混区总是存在,且返混区大小与高度有关：当表观液遣与表观气速之比大于19．6时,返混消失,含气率分布由塔中心峰值转向壁面峰值。径向液相速度既与气相表现速度有关又与位置高度有关,在塔底部呈现负值,这意味着向塔轴心方向流动。随着塔高增加。流动方向逐渐转变为向塔壁方向,且又有明显的峰值。     

加温加压下CFD-PBM耦合模型空气-水两相流数值模拟研究

计算流体力学与群体平衡模型（CFD-PBM）结合可有效地模拟鼓泡塔内流体行为,较准确地预测流场、相含率以及局部气泡尺寸分布。以直径100 mm、高1.3 m的加温加压鼓泡塔为模拟对象,在系统压力为1 MPa、表观气速为0.08~0.24 m/s、温度为30~160℃条件下系统地考察了空气-水体系的表观气速、温度以及固含率对平均气含率、大小气泡气含率、气泡直径和气泡尺寸分布等参数的影响。结果表明,平均气含率的模拟结果和实验值在10%的误差范围内吻合较好;温度的变化主要影响了塔内气泡的聚并和破碎,并用聚并破碎的机理解释了温度对其流体行为的影响。     

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.     

Generation of uniform small bubbles and hydrodynamic characterization of a bubble column with high pressure jet sparger

Systems generating uniform small bubbles are used in many mineral processing and chemical operations. We investigated the generation of smaller bubbles by using a two fluid jet system. Gas holdup results are reported in terms of the effect of superficial gas and liquid velocities in relation to the pressure in a bubble column with a water jet sparger. Experiments were conducted with hydrostatic head of 80 cm, 100 cm, and 120 cm in the bubble column. The gas velocity varied from 0.122 to 1.22 cm/s, and water flow rate from 33.3 to 333 cm³/s. Experiments were conducted at pressures of 2 atms., 3 atms. 4 atms. and 5 atms., and bubble sizes were measured by a digital camera (bubble compared to a reference wire inside the bubble column). Results show that the gas holdup increases with the pressure and superficial gas velocities; and at pressures of 2, 3, 4 and 5 atms., the gas holdup increases by 8.75%, 9.166%, 10% and 10%, respectively. The maximum gas holdup of 16.4% was observed at a liquid level of 80 cm and pressure of 4 atms. Optimum conditions for generating smaller bubbles with larger gas holdup are increased liquid flow rate, low liquid level, and high gas pressure. Experimental results also indicate that the column operates in both the homogeneous and heterogeneous regimes of gas-liquid flow.     

Bubble formation in flowing liquid

Initial bubbles in flowing liquid from a nozzle were observed from two mutually perpendicular directions. Two nozzles of 0.086 cm and 0.305 cm in diameter were used. The gas flow rate and the superficial liquid velocity ranged from 0.33 cm³/s to 36.2 cm³/s and from 0 cm/s to 154.9 cm/s, respectively. The bubble size formed in flowing liquid decreased with decreasing gas flow rate and with increasing superficial liquid velocity. Three types of bubble formation, i.e. single bubbles, coalescent bubbles and gas jets, were observed depending upon the gas rate and the liquid velocity. Two empirical equations of the bubble sizes are given.     

EFFECT OF STATIC LIQUID HEIGHT ON GAS-LIQUID MASS TRANSFER IN A DRAFT-TUBE BUBBLE COLUMN AND THREE-PHASE FLUIDIZED BED

Experiments are performed under batch-liquid operating conditions to investigate the effect of static liquid height on the gas-liquid mass transfer coefficient (K_La) in a draft-tube bubble column (DTBC) and a draft-tube three-phase fluidized bed (DTFB). In addition, the effects of column diameter, gas-distributor, and draft-tube diameter are studied. The results indicate that for a given system with a porous plate gas-distributor at low superficial gas velocities (<70 m/hr), increasing static liquid height decreases K_La. At high gas velocities, K_La is independent of the static liquid height. For systems with a perforated gas-distributor, there is no effect of static liquid height on K_La. The formation of small dispersed bubbles at low gas velocities in the porous plate distributor system accounts for the considerably high K_La values and the observed effect of liquid height. On the other hand, the formation of large spherical-cap bubbles and the bubble coalescence at high gas velocities reduce the performance of the porous plate distributor system to that of the perforated one.     

BUBBLE MEASUREMENTS IN A GAS-LIQUID JET

Measurements have been carried out in the developing and fully developed regions of a free, axisymmetric, isothermal, air-water, bubbly jet. Three experiments have been conducted at a fixed jet-exit Reynolds number and gas superficial velocity using three different bubble injector assemblies producing bubbles of moderately different average sizes and size distributions. The volume fraction of the bubbly jet flow examined in this study is low and the resulting dispersed flow is dilute. A one-component Phase-Doppler Velocimetry system has been employed to measure bubble size and velocity non-intrusively. Visual data collected simultaneously with the light-scattering measurements were analyzed with the aid of image processing and used to verify the trends portrayed by the light-scattering measurements and to determine average bubble size. Our measurements show that, even in the dilute flow examined here, differences in initial bubble-size and size distribution can influence the RMS velocity fluctuations of the bubbles, particularly in the jet development region. The average bubble velocities are less sensitive. Evidence that the development pattern of the jet near the exit is affected by the presence of the bubbles is also presented. Near the exit of the jet, bubbles are shown to be ejected laterally outside the jet due to the significant lift force caused by the high velocity gradient in the axisymmetric shear layer. The observed sensitivity of the bubble flow to size-related parameters and initial conditions in this dilute case, indicates that discrepancies in previous measurements of dispersed, bubbly flows could be attributed to different size characteristics and/or initial conditions.     

Destabilisation of homogeneous bubbly flow in an annular gap bubble column

Experimental results are presented to show that there are very significant differences in the mean gas void fractions measured in an open tube and a annular gap bubble column, when operated at the same gas superficial velocity, using a porous sparger. Measurements were carried out in a vertical 0.102 m internal diameter column, with a range of concentric inner tubes to form an annular gap, giving diameter ratios from 0.25 to 0.69; gas superficial velocities in the range 0.014–0.200 m/s were investigated. The mean gas void fraction decreases with increasing ratio of the inner to outer diameter of the annular gap column and the transition to heterogeneous flow occurs at lower gas superficial velocities and lower void fractions. Two reasons are proposed and validated by experimental investigations: (1) the presence of the inner tube causes large bubbles to form near the sparger, which destabilise the homogeneous bubbly flow and reduce the mean void fraction; this was confirmed by deliberately injecting large bubbles into a homogeneous dispersion of smaller bubbles, and (2) the shape of the void fraction profiles changes with gap geometry and this affects the distribution parameter in the drift‐flux model. Both of these effects serve to reduce the mean gas void fraction in an annular gap bubble column compared to an open tube at the same gas superficial velocity.     

Heat transfer and bubble dynamics in a three-phase inverse fluidized bed

In this study, surface-to-bed heat transfer experiments were performed to gain insight on heat transfer and hydrodynamics in a three-phase inverse fluidized bed. Air, tap water or 0.5 wt.% aqueous ethanol, and polypropylene were, respectively, the gas, liquid and solid phases. The solid loading was varied from 0 to 30 vol.%, and the gas and liquid superficial velocities from 2 to 50 mm/s and 0 to 21 mm/s, respectively. Visual observations were associated with measured phase holdups and instantaneous heat transfer coefficients. Larger gas velocities lead to an increase in bubble size due to the transition to the coalesced bubble flow regime. The greater turbulence induced by the larger bubbles increases the average heat transfer coefficient. On the other hand, adding ethanol reduces the heat transfer coefficient. Solid concentrations up to ∼13 vol.% increase the average heat transfer coefficient whereas higher solid concentrations tend to lower it. The distribution of instantaneous heat transfer coefficient peak height is wider at higher gas and liquid velocities while the addition of a surfactant narrows it. Gas holdups and average heat transfer coefficients are both compared with existing correlations, which are then adjusted for a better fit.     

Local gas–liquid slip velocity distribution in bubble columns and its relationship with heat transfer

The knowledge of the local gas–liquid slip velocity distribution can offer a better understanding for the complex transport phenomena in bubble columns. In this work, CFD–PBM simulations are carried out to investigate the effect of superficial gas velocities, axial positions, and scale of bubble columns on the time-averaged radial profiles of gas–liquid slip velocities. Furthermore, the relationship between local slip velocities and local heat transfer coefficients in pilot-scale bubble columns at superficial gas velocities of 0.05 m/s, 0.20 m/s, and 0.35 m/s is studied. The results indicate that the slip velocities decrease with the increase of r/R (r-radial position, R-column radius), while increase with increasing superficial gas velocities in general. In the fully developed region, the axial positions have small impact on the local slip velocities. A strong linear relation between heat transfer coefficients and slip velocities in the fully flow developed region is observed.     

Use of models for lift, wall and turbulent dispersion forces acting on bubbles for poly-disperse flows

Closure laws are needed for the qualification of CFD codes for two-phase flows. In case of bubbly and slug flow, forces acting on the bubbles usually model the momentum transfer between the phases. Several models for such forces can be found in Literature. They show, that these forces depend on the liquid flow field as well as on the size and the shape of the bubbles. A validation of consistent sets of bubble force models for poly-disperse flows is given, basing on a detailed experimental database for vertical pipe flows, which contains data on the radial distribution of bubbles of different size as well as local bubble size distributions. A one-dimensional (1D) solver provides velocity profiles and bubble distributions in radial direction. It considers a large number of bubble size classes and is used for the comparison with the experiments. The simplified model was checked against the results of full 3D simulations done by the commercial code CFX-5.7 for simplified monodisperse cases. The effects of the number of bubbles classes as well as the effect of the lateral extension of the bubbles were analyzed. For the validation of bubble force models measured bubble size distributions were taken as an input for the calculation. On basis of the assumption of an equilibrium of the lateral bubble forces, radial volume fraction profiles were calculated separately for each bubble class. In the result of the validation of different models for the bubble forces, a set of Tomiyama lift and wall force, deformation force and Favre averaged turbulent dispersion force was found to provide the best agreement with the experimental data. Some discrepancies remain at high liquid superficial velocities.     

CFD Simulations of a Bubble Column Operating in the Homogeneous and Heterogeneous Flow Regimes

Bubble columns are operated either in the homogeneous or heterogeneous flow regime. In the homogeneous flow regime, the bubbles are nearly uniform in size and shape. In the heterogeneous flow regime, a distribution of bubble sizes exists. In this paper, a CFD model is developed to describe the hydrodynamics of bubble columns operating in either of the two flow regimes. The heterogeneous flow regime is assumed to consist of two bubble classes: “small” and “large” bubbles. For the air‐water system, appropriate drag relations are suggested for these two bubble classes. Interactions between both bubble populations and the liquid are taken into account in terms of momentum exchange, or drag‐, coefficients, which differ for the “small” and “large” bubbles. Direct interactions between the large and small bubble phases are ignored. The turbulence in the liquid phase is described using the k‐ϵ model. For a 0.1 m diameter column operating with the air‐water system, CFD simulations have been carried out for superficial gas velocities, U, in the range 0.006–0.08 m/s, spanning both regimes. These simulations reveal some of the characteristic features of homogeneous and heterogeneous flow regimes, and of regime transition.     

Flow behaviors of multiscale bubbles with porous ceramic membrane distributors: Visualization and modeling

The distribution of gas-liquid two-phase flow is one of the significant effects on heterogeneous catalytic reactions. Ceramic membrane gas distributors (CMGD) were applied in improving gas-liquid distribution, and flow behavior of gas as dispersed phase in liquid phase was visualized via high-speed photograph. The average diameters of multiscale bubbles were measured and modeled ranging from 10⁻⁵ to 10⁻² m. The coalescence and trajectory of bubbles during rising process were observed, and two typical trajectories straight and spiral types were tracked. To inhibit coalescence of bubbles during rising process, internals manufactured by 3D printing were installed in the channel of ceramic membrane. The average bubble size of CMGD decreases 12% from 392 to 345 μm compared with that of the original CMGD. The CMGD with internals enhances the heterogeneous catalytic reaction performance via providing large quantity of stabile multiscale bubbles which could match the porous structure of catalyst.