Effects of Orifice Orientation and Gas-Liquid Flow Pattern on Initial Bubble Size

In many gasliquid processes, the initial bubble size is determined by a series of operation parameters along with the sparger design and gasliquid flow pattern. Bubble formation models for variant gasliquid flow pat terns have been developed based on force balance. The effects of the orientation of gasliquid flow, gas velocity, liquid velocity and orifice diameter on the initial bubble size have been clarified. In ambient airwater system, thesultable gasllquid flow pattern is important to obtain smaller bubbles under the low velocity liquid crossflow con ditions with stainless steel spargers. Among the four types of gasliquid flow patterns discussed, the horizontal orifice in a vertically upward liquid flow produces the smallest initial bubbles. However the orientation effects of gas and liquid flow are found tobe insgnifican whenliq.uid velocity is.higher than. 3.2 m;sa or theorifice diameter is small enough.     

Shape oscillations and path transition of bubbles rising in a model bubble column

We investigate experimentally the occurrence of shape oscillations accompanied by path transition of periodically produced air bubbles rising in water. Within the period of bubble formation, the induced velocity is measured to examine bubble-liquid and bubble-bubble interactions. The flow is produced in a small-scale bubble column with square-shaped cross section. A capillary aerator produces bubbles of size 3.4 mm at a frequency of 5 Hz. Measuring techniques employed are high-speed imaging to capture bubble shape oscillations and path geometry, and laser-Doppler anemometry (LDA) to measure the velocity in the liquid near the rising bubbles. The experimentally obtained bubble shape data are expanded in Legendre polynomials. The results show the occurrence of oscillations by the periodicity of the expansion coefficients in space. Significant shape oscillations accompanied by path transition are observed as the second-mode oscillation frequency converges to the frequency of the initial shape oscillations. The mean velocity field in the water obtained by LDA agrees well with potential theory. An analysis of the decay of the induced flow shows that there is no interaction between the flow fields of two succeeding 3.4 mm bubbles in the rectilinear path when the bubble production frequency is lower than 7.4 Hz.     

PIV study of bubble rising behavior

Bubble behaviors are studied in a rectangular bubble column using particle image velocimetry (PIV). Bubbles rise in a chain through a stagnant liquid. As liquid viscosity reduces, bubble rising trajectory changes from one-dimensional to three-dimensional. This transition is due to different bubble wake structures. The bubble shapes also show dissimilar characteristics in liquid of different viscosities. The instantaneous liquid flow fields measured by PIV show the diversities as bubbles rise in different paths. Based on the experimental data obtained in this study, a correlation is proposed to calculate the terminal velocity of bubbles in a chain.     

Flow in the nose region and annular film around a Taylor bubble rising through vertical columns of stagnant and flowing Newtonian liquids

The flow in the nose region and in the annular film around individual Taylor bubbles rising through stagnant and co-current vertical columns of liquid were studied, employing particle image velocimetry (PIV) and pulsed shadowgraphy techniques (PST) at the same time. The combined techniques enabled simultaneous determination of the bubble shape and the velocity profiles in the liquid film. Experiments were performed with water and aqueous glycerol solutions in a wide range of viscosities , in an acrylic column of 32 mm ID.Values for the distance ahead of the nose in which the flow is disturbed by the presence of the bubble are presented for the conditions studied. The bubble shapes in the nose region are compared with Dumitrescu's shape for potential flow. The velocity profiles show that after the nose region the liquid begins to accelerate downwards, and at a certain distance from the bubble nose the velocity profile and the liquid film thickness stabilise. The liquid film acquires characteristics of a free-falling film. Values of the developing length and film thickness are reported for the experimental conditions studied. Average velocity profiles in the fully developed film are also presented. A critical Reynolds number of around 80 (based on the mean absolute velocity in the liquid film and on the film thickness) is reported for the transition from laminar to turbulent regime. Shear stress profiles (in the fully developed film) are also provided.The data reported are relevant for the validation of numerical codes in slug flow.     

Characteristics of flow fields in the gas–liquid mini-bubble columns with particle image velocimetry measurements

As a new type of gas–liquid microreactors, the gas–liquid mini-bubble column has potential applications. However, few studies on the flow fields in the mini-bubble column can be found at present. In this work, particle image velocimetry (PIV) was used to visually study the velocity fields, vorticity fields and bubble dynamics in the gas–liquid mini-bubble columns with column inner diameters of 1–3 mm and mini-bubble diameters ranged from 0.7 to 1.3 mm. It is found that with the increase of superficial liquid velocity, bubbles rose from almost straight line to Z-shaped or S-shaped trajectory, and the bubble trajectory changed from one-dimension to three-dimension; when the bubble velocity changed, the bubble size and gas holdup decreased; bubble terminal velocity was controlled by bubble buoyancy and flow resistance, and increased slightly with bubble coalescence. These findings may provide basic reference for the design and scale-up of such a mini-bubble column reactor.     

Experimental Study on Bubble Rising and Descending Velocity Distribution in a Slurry Bubble Column Reactor

To determine bubble rising and descending velocity simultaneously, a BVW‐2 four‐channel conductivity probe bubble parameters apparatus and its analysis are used in gas‐liquid and gas‐liquid‐solid bubble columns. The column is 100 mm in internal diameter and 1500 mm in height. The solid particles used are glass beads with an average diameter of 17.82 μm, representing typical particle size for catalytic slurry reactors. The effects of superficial gas velocity (1.0 cm/s ≤ U_g ≤ 6.4 cm/s), solid holdup (0 % ≤ ?_s ≤ 30 %), and radial location (r/R = 0, 0.4, and 0.7) on bubble velocity distributions are determined. It is found that increasing U_g can increase the velocity of bubbles but do not exert much influence on bubble velocity distribution. Solid holdup mainly affects the distribution of bubble velocity while the radial direction affects bubble velocity distribution only slightly. The ratio of descending bubbles to rising bubbles increases from the bubble column center to the wall. It can be proved experimentally that large bubbles do not always rise faster than small bubbles at higher U_g (for example 6.4 cm/s).     

固含率对沸腾床反应器气泡行为的影响

采用动态气体逸出法,在高7.0 m、直径0.3 m的有机玻璃塔中研究了固含率对沸腾床反应器内气泡行为特性的影响。在表观气速2.16~21.62 cm/s和固含率9.8%~39.0%（体积分数）范围内测定了反应器内的总气含率、大小气泡含率、大小气泡上升速度及其尺寸等参数。结果表明：总气含率随着表观气速的增大而增大,随着固含率的增大而减小。随着表观气速的增大,大气泡含率、大气泡直径及其上升速度均呈增大趋势;小气泡含率明显增大,但小气泡上升速度和直径趋于减小。随着固含率的增大,大气泡含率略有降低,但大气泡直径及其上升速度都明显增大;当固含率超过19.5%（体积分数）后,小气泡上升速度几乎下降为0;当固含率达到29.3%（体积分数）时,小气泡基本消失。     

Flow patterns in the wake of a Taylor bubble rising through vertical columns of stagnant and flowing Newtonian liquids: An experimental study

The flow in the wake and near-wake regions of individual Taylor bubbles rising through stagnant and co-current vertical columns of Newtonian liquids was studied, employing simultaneously particle image velocimetry (PIV) and pulsed shadowgraphy techniques (PST). Experiments were made with water and aqueous glycerol solutions covering a wide range of viscosities , in an acrylic column of 32 mm ID.Different wake structures (laminar, transitional and turbulent) are identified, in both stagnant and co-current flow conditions. In stagnant liquids, the wake flow pattern is only dependent on the dimensionless group N_f. The different types of wakes obtained are in accordance with the critical N_f numbers proposed in previous works. For co-current flow conditions, the flow patterns in the wake depend on the Reynolds number based on the relative (to the bubble) average velocity of the upward liquid flow, the laminar-transitional and transitional-turbulent limits being for the first time experimentally determined.The wake flow patterns are quantified by means of instantaneous and average flow fields. Values for the wake length and wake volume are also presented and compare well with correlations found in literature. Study of the flow in the near-wake zone enabled determination of the distance needed to recover the undisturbed liquid velocity profile.The detailed study of the flow in the wake and near-wake regions is an important contribution to better understanding the interaction and coalescence mechanisms between Taylor bubbles.The data reported are relevant to the validation of numerical simulation codes in the vertical slug flow regime.     

黏性流体中单气泡的运动特性

采用高速CCD成像技术实验测定了单个气泡在水及不同黏度甘油水溶液中的形变特性及上升终速度.结果表明,随着液相黏度的增加,小气泡均能呈良好球形,而大气泡则由类似椭球形的不规则形状逐渐向球帽形转变;气泡的周期性振荡趋于平缓,气泡上升终速度降低.提出了计算气泡上升终速度的关联方程,其预测结果与实验值吻合良好.     

Interaction of two in-line bubbles of equal size rising in viscous liquid

The interaction of bubbles is the key to understand gas–liquid bubbling flow. Two-dimensional axis-symmetry computational fluid dynamics simulations on the interactive bubbles were performed with VOF method,which was validated by experimental work. It is testified that several different bubble interactive behaviors could be acquired under different conditions. Firstly, for large bubbles(d: 4, 6, 8, 10 mm), the trailing bubble rising velocity and aspect ratio have negative correlations with liquid viscosity and surface tension. The influences of viscosity and surface tension on leading bubble are negligible. Secondly, for smaller bubbles(d: 1, 2 mm), the results are complicated. The two bubbles tend to move together due to the attractive force by the wake and the potential repulsive force. Especially for high viscous or high surface tension liquid, the bubble pairs undergo several times acceleration and deceleration. In addition, bubble deformation plays an important role during bubble interaction which cannot be neglected.     

Bubble formation in cocurrently upward flowing liquid

The effects of liquid velocity, nozzle diameter, gas chamber volume and gas flow rate on volumes, shapes and growth curves of bubbles formed at a nozzle submerged in a cocurrently upward flowing liquid in a bubble column were experimentally investigated. The bubble volume decreases with increasing liquid flow velocity. The effect of liquid flow velocity on the volume of bubble increases with an increase in the gas flow rate. To simulate bubble formation at a nozzle submerged in cocurrently upward flowing liquid, a revised non-spherical bubble formation model was proposed. Bubble volumes, bubble growth curves and shapes experimentally obtained in this study, as well as in previous experimental studies, are well predicted by the present model.     

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.     

Drag force and clustering in bubble swarms

In this article, results of detailed numerical simulations are reported meant to provide a closure relation for the drag force acting on bubbles rising in a dense swarm. The formation of clusters of bubbles in a periodic domain and the effect thereof on the rise velocity and effective drag coefficient on the bubbles are studied. Using smaller bubble sizes than presented in our earlier work, we are also able to refine our correlation for the drag coefficient acting on bubbles rising in a swarm, such that it is applicable for a large range of bubble sizes. The simulations are performed with an advanced Front‐Tracking model in which Lagrangian marker points are used to track the gas–liquid interface, while accounting for surface tension and substantial interface deformation. Simulations were performed using periodic domains to simulate rising air bubbles in water from 1.0 mm up to 6.0 mm in diameter. The effect of liquid phase viscosity was also studied to extend the range of validity of the drag correlation. For the 1.0 and 1.5 mm cases, strong horizontal clustering effects are observed. Especially, at high gas fractions, the bubbles tend to form rigid horizontal arrays, which have been shown to strongly increase the drag force acting on the bubbles in the cluster. For viscous liquids, the tendency to form horizontal clusters is lower, and even vertical clustering is observed. The bubble slip velocity was compared with the experimental results of Zenit et al., which agree very well taking into account the differences between simulations and experiments. Based on our simulations, a new drag correlation was proposed, taking into account Eötvös numbers ranging from 0.13 to 4.9, and Morton numbers in the range 3.8 ≤ ? log Mo < 6.6, and gas hold‐ups up to 40% (30% for Eo < 0.3). At lower values for ?log Mo, the Reynolds number drops to the order of unity, and the correlation overpredicts the drag coefficient, which defines the range of applicability of the currently proposed drag correlation. The correlation itself describes a linear increase of the normalized drag coefficient as a function of the gas hold‐up. The strength of linear increase is stronger at lower Eötvös numbers. © 2012 American Institute of Chemical Engineers AIChE J, 59: 1791–1800, 2013     

Study of bubble induced flow structure using PIV

An experimental investigation of the flow structure induced by a chain of gas bubbles was carried out in a rectangular bubble column using particle image velocimetry (PIV). It is observed that the bubble rising trajectory changes from one dimension to three dimension as liquid viscosity reduces. The variation of bubble rising trajectory associates with the alternation of bubble motions—with or without oscillatory and rotational motion depending the bubble rising trajectory is 3-D or 1-D. The different behaviors of gas bubbles introduce various instantaneous and averaged liquid flow structures. In general, complex fluid velocity fields present in liquid system of low viscosity where free vortex, cross flow, and irregular circular flow can be observed. The liquid pseudo-turbulence measured in terms of turbulence intensity and Reynolds stress is more intense in liquid of low viscosity. The turbulence is also enhanced by the frequency of bubble formation.     

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.     

非牛顿流体中的两并行运动气泡周围流场及气泡间相互作用的PIV研究

The flow fields surrounding two parallel moving bubbles rising from two identical orifices submerged in non-Newtonian fluid of carboxymethylcellulose (CMC) solution of three different mass concentration were measured experimentally by the use of particle image velocimetry (PIV). The influences of gas flowrate, solution mass concentration, orifice interval and the angle between two bubble centers line and vertical direction on the flow field surrounding bubbles were discussed respectively by analyzing the velocity vector, velocity contours as well as individual velocity components. The results show that the liquid velocity both in front of two bubbles and behind increases with gas flowrate duo to shear-thinning effect of previous bubbles, whereas decreases with the increase of CMC concentration due to the increase of drag force acting on bubbles. The effect of the orifice interval on the flow field around two moving bubbles becomes gradually obvious as the interval becomes closer. Moreover, two adjacent side-by-side bubbles repulse each other during rising, leading to the practical interval between them increased somewhat above the orifice interval. When the distance between bubbles is less than the orifice interval l0 mm, the interaction between two neighboring bubbles changed from mutual repellence to attraction with the decrease of the angle of the line of linking two bubble centers to the vertical direction.     

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.     

幂律流体中平行上升双气泡间相互作用研究

采用流体体积(VOF)法对非牛顿幂律流体中平行上升双气泡间相互作用进行了模拟研究.结果表明:剪切变稀效应及气泡周围流场结构对气泡间相互作用的贡献分别随着气泡间距增大而减小,气泡间相互作用主要受气泡间流体涡旋结构控制.当初始间距很小时,气泡平行聚并发生;随着初始间距增大,双气泡行为转变为相互影响的平行上升;当初始间距增大...     

Numerical simulation of single bubble motion in ionic liquids

In this work a numerical simulation is studied to investigate the motion of single bubble in ionic liquids using an improved volume-of-fluid (VOF) method. In the improved method, besides the gravity and surface tension, a new drag force is added to the momentum equation in order to describe the gas–liquid interaction in the ionic liquids, which possess some special properties compared with the traditional solvents. The deformation, velocity and equivalent diameter of single bubble rising in three ionic liquids, i.e., bmimBF₄, bmimPF₆ and omimBF₄, are simulated and the calculation results agree well with the experimental data. Furthermore, the detailed velocity fields and pressure fields around the bubbles are predicted with the proposed numerical simulation model. This work is important for understanding the fluid dynamic performance of bubbles in ionic liquids, and could provide a useful tool for designing a bubble column with ionic liquids as its solvents.     

Behavior of a single coherent gas bubble chain and surrounding liquid jet flow structure

The motion of a single nitrogen gas bubble chain and the structure of water flow field surrounding the chain were experimentally studied. We developed a bubble generator that can control both the bubble diameter and the generation frequency independently. Experimental conditions of bubble Reynolds number and bubble distance divided by bubble diameter were from 300 to 650 and from 6.5 to 300, respectively. We discuss the interaction effects on the motion of each bubble rising in a chain, as compared to the effects of a single rising bubble. The bubble trajectories and the surrounding water flow fields in the state of bubbles rising in a chain were investigated using a high-speed digital video camera and an analog single-lens-reflex camera. We observed two important physical phenomena. First, bubbles passed through a nearly identical path in the case of low frequency of bubble production. On the contrary, at a height of approximately 50 mm from the nozzle, the bubbles in the case of high frequency deviated and scattered from this path due to bubble-bubble interaction. Second, with higher bubble production frequency, coherent bubble chain and the characteristic structure of the surrounding water flow called “liquid jet” were observed near the nozzle. The direction of liquid jet flow differed from the bubble trajectory. We theoretically investigated the relation of coherent bubble chain and liquid jet by applying the conservation of liquid momentum.