Influence of gas injection on phase inversion in an oil–water flow through a vertical tube

An experimental study has been made of the influence of gas injection on the phase inversion between oil and water flowing through a vertical tube. Particular attention was paid to the influence on the critical concentration of oil and water where phase inversion occurs and on the pressure drop increase over the tube during phase inversion. By using different types of gas injectors also the influence of the bubble size of the injected gas on the phase inversion was studied. It was found that gas injection does not significantly change the critical concentration, but the influence on the pressure drop is considerable. For mixture velocities larger than 1 m/s, the pressure drop over the tube increases with decreasing bubble size and at inversion can become even larger than the pressure drop during the flow of oil and water without gas injection.     

Air–water flow in a vertical pipe: experimental study of air bubbles in the vicinity of the wall

This study deals with the influence of bubbles on a vertical air–water pipe flow, for gas-lift applications. The effect of changing the bubble size is of particular interest as it has been shown to affect the pressure drop over the pipe. Local measurements on the bubbles characteristics in the wall region were performed, using standard techniques, such as high-speed video recording and optical fibre probe, and more specific techniques, such as two-phase hot film anemometry for the wall shear stress and conductivity measurement for the thickness of the liquid film at the wall. The injection of macroscopic air bubbles in a pipe flow was shown to increase the wall shear stress. Bubbles travelling close to the wall create a periodic perturbation. The injection of small bubbles amplifies this effect, because they tend to move in the wall region; hence, more bubbles are travelling close to the wall. A simple analysis based on a two-fluid set of equations emphasised the importance of the local gas fraction fluctuations on the wall shear stress.     

Effects of the gas phase molecular weight and bubble size on effervescent atomization

The effervescent atomization from an industrial Coker feed nozzle is compared for two different gas densities (air and mixed gas of 81.4 vol.% helium/18.6 vol.% nitrogen) at equivalent operating temperatures. The application is to observe the similarity of lab tests using air at 20 °C to the industrial process using steam at 300-400 °C. The effects of operating conditions, such as gas to liquid mass ratio, mixing pressure and void fraction on the flow regime, bubble size, and droplet size distribution were also examined in this study. The experiments were performed using mixtures of water with air or mixed gas, which resulted in gas to liquid mass ratios ranging from 1% to 4%.Stroboscopic back scattered imagery (SBSI) indicates that the average bubble size inside the nozzle conduit is similar when air and water are used as the process fluids, when compared to the case when mixed gas and water are used as the process fluids. Under similar conditions, the Phase Doppler Particle Anemometer (PDPA) data indicate that the droplet size in the spray is similar when using either mixed gas or air as the atomization gas.Experimental results obtained by high-speed video shadowgraphy (HSVS) indicate that the flow pattern inside the nozzle feeding conduit was slug flow with a tendency to attain annular flow with increased air to liquid mass ratios. Thus, from the experimental results it is evident that the smaller molecular weight of the mixed gas versus air (8.4 versus 29) does not significantly reduce the bubble (<±10% difference) and droplet size (<±1.5% difference), indicating a weak dependence of the gas phase density on two-phase atomization. This confirms that laboratory experiments on effervescent nozzles using air have reliable similarity to systems that use high temperature steam for the gas phase.     

Measurement of multiple gas-bubble velocities in gas–liquid flows using hot-film anemometry

A dual-probe hot-film anemometry technique has been developed to measure multiple gas-bubble velocities corresponding to different gas-bubble size groups in air–water flows. A data reduction scheme using wavelet analysis combined with a phase-detection technique is used to discriminate the hot-film anemometer output signals into signals corresponding to different bubble size groups. The phase and bubble size discrimination is based on the magnitude and time derivative of the signal, and the streamwise length of the gas bubbles. A cross-correlation between the discriminated signals from the two probes yields an average time difference of arrival of the gas bubbles at the two sensor locations. The velocities are estimated from the distance between the sensors and the time difference of arrival. The mean bubble size is estimated from the chord length distribution. Measurements performed in vertical-up air–water slug flow show the technique to be a viable method for obtaining bubble velocity and size information. The velocity measurements from the hot-film anemometry are corroborated using a high-speed quantitative flow visualization system. Received: 22 December 1999/Accepted: 8 May 2001     

水气两相流中稀疏气泡流动速度场的数字图像测量初探

由空压机提供的气体通过—排微小直径的喷嘴进入静止水体，形成水气两相流流场。在单相PIV和PTV技术的基础上，研究稀疏气液两相流情况下气泡的速度场分布。PIV算法采用快速傅立叶互相关分析法，而PTV算法需要获得每幅图像中每个气泡的形心，根据连续图像中的粒子对，计算速度。用PIV和PTV两种算法处理求出气泡的速度并对两种方法进行比较，其最终研究成果可应用于流体及多相流的流量测技术，提高我们进行低密度气液两相流相关研究的测量水平。同时为水气两相流的数值分析和理论研究提供流场测试的数据。     

Development of co-current air–water flow in a vertical pipe

Measurements of the cross-sectional distribution of the gas fraction and bubble size distributions were conducted in a vertical pipe with an inner diameter of 51.2 mm and a length of about 3 m for air/water bubbly and slug flow regimes. The use of a wire-mesh sensor obtained a high resolution of the gas fraction data in space as well as in time. From this data, time averaged values for the two-dimensional gas fraction profiles were decomposed into a large number of bubble size classes. This allowed the extraction of the radial gas fraction profiles for a given range of bubble sizes as well as data for local bubble size distributions. The structure of the flow can be characterized by such data. The measurements were performed for up to 10 different inlet lengths and for about 100 combinations of gas and liquid volume flow rates. The data is very useful for the development and validation of meso-scale models to account for the forces acting on a bubble in a shear liquid flow and models for bubble coalescence and break-up. Such models are necessary for the validation of CFD codes for the simulation of bubbly flows.     

Local gas and liquid phase velocity measurement in a miniature stirred vessel using PIV combined with a new image processing algorithm

A method which combines standard two-dimensional particle image velocimetry (PIV) with a new image processing algorithm has been developed to measure the average local gas bubble velocities, as well as the local velocities of the liquid phase, within small stirred vessel reactors. The technique was applied to measurements in a gas–liquid high throughput experimentation (HTE) vessel of 45 mm diameter, but it is equally suited to measurements in larger scale reactors. For the measurement of liquid velocities, 3 μm latex seeding particles were used. For gas velocity measurements, a separate experiment was conducted which involved doping the liquid phase with fluorescent Rhodamine dye to allow the gas–liquid interfaces to be identified. The analysis of raw PIV images enabled the detection of bubbles within the laser plane, their differentiation from obscuring bubbles in front of the laser plane, and their use in lieu of tracer particles for gas velocity analysis using cross-correlation methods. The accuracy of the technique was verified by measuring the velocity of a bubble rising in a vertical glass column. The new method enabled detailed velocity fields of both phases to be obtained in an air–water system. The overall flow patterns obtained showed a good qualitative agreement with previous work in large scale vessels. The downward liquid velocities above the impeller were greatly reduced by the addition of the gas, and significant differences between the flow patterns of the two-phases were observed.     

Modeling of bubble behaviors and size distribution in a slab continuous casting mold

Population balance equations combined with Eulerian–Eulerian two-phase model are employed to predict the polydispersed bubbly flow inside the slab continuous-casting mold. The class method, realized by the MUltiple-SIze- Group (MUSIG) model, alongside with suitable bubble breakage and coalescence kernels is adopted. A two-way momentum transfer mechanism model combines the bubble induced turbulence model and various interfacial forces including drag, lift, virtual mass, wall lubrication, and turbulent dispersion are incorporated in the model. A 1/4th scaled water model of the slab continuous-casting mold was built to measure and investigate the bubble behavior and size distribution. A high speed video system was used to visualize the bubble behavior, and a digital image processing technique was used to measure the mean bubble diameter along the width of the mold. Predictions by previous mono-size model and MUSIG model are compared and validated against experimental data obtained from the water model. Effects of the water flow rate and gas flow rate on the mean bubble size were also investigated. Close agreements by MUSIG model were achieved for the gas volume fraction, liquid flow pattern, bubble breakage and coalescence, and local bubble Sauter mean diameter against observations and measurements of water model experiments.     

Experimental investigation of the influence of column scale,gas density and liquid properties on gas holdup in bubble columns

Measurements of gas holdups in bubble columns of 0.16, 0.30 and 0.33 m diameter were carried out. These columns were operated in co-current flow of gas and liquid phases and in semibatch mode. The column of 0.33 m diameter was operated at elevated pressures of up to 3.6 MPa. Nitrogen was employed as the gas phase and deionized water, aqueous solutions of ethanol and acetone and pure acetone and cumene as the liquid phase. The effects of differing liquid properties, gas density (due to elevated pressure), temperature, column diameter and superficial liquid velocity on gas holdup were studied. The gas holdup measurements were utilized by differential pressure measurements at different positions along the height of the bubble columns which allowed for the identification of axial gas holdup profiles. A decrease of gas holdup with increasing column diameter and an increase of gas holdup with increasing pressure was observed. The effect of a slightly decreasing gas holdup with increasing liquid velocity was found to exist at smaller column diameters. The use of organic solvents as the liquid phase resulted in a significant increase in gas holdup compared to deionized water. It is found that published gas holdup models are mostly unable to predict the results obtained in this study.     

Computational study of aeration for wastewater treatment via ventilated pump-turbine

Large eddy simulations were performed on a modular pump-turbine to study oxygen dissolution inside the draft tube. Air injection was applied over the runner cone surface during turbine operation. Data regarding bubble size, void fraction and interfacial area concentration were presented to understand their influence on oxygen dissolution. Transient single phase and multiphase flow simulations were carried out to investigate the influence of air injection and dissolution within the flow field and turbine performance. Multiphase simulations were conducted by using the mixture multiphase model. The mathematical modeling of oxygen dissolution employed was validated by comparing predicted oxygen dissolution against experimental measurements performed by Zhou et al. (2013). The averaged dissolved oxygen concentration in the range of 1.2–1.4 mg/l was obtained; which is sufficient for an active aerobic microorganism activity for wastewater treatment processes. Dissolution efficiency and the amount of averaged dissolved oxygen inside the draft tube were sensitive to the inlet bubble size. The efficiency of the dissolution increases strongly as the inlet bubble size was reduced. The obtained results revealed that vortex suppression was achieved through air admission within multiphase flow simulation. Moreover, the power generation of the turbine was hardly influenced by the aeration through the runner cone.     

An advanced LIF-PLV system for analysing the hydrodynamics in a laboratory bubble column at higher void fractions

Bubble columns are widely used in the chemical industry and biotechnology. Flow and turbulence in such an apparatus are induced by the bubble rise, and the bubble behaviour is strongly affected by swarm effects (i.e. the interaction between bubbles). For analysing the bubble swarm behaviour and simultaneously evaluating the flow structure and bubble-induced turbulence, a bubble column of 140 mm diameter and a height of 650 mm or 1,400 mm (initial water level) were considered. The bubble column was aerated with relatively fine bubbles having a mean size between about 0.5 and 4.0 mm. The gas hold-up was varied in the range between 0.5 and 19%. A two-phase pulsed-light velocimetry (PLV) system was developed to evaluate instantaneous flow fields of both rising bubbles and the continuous phase. The measurement of the liquid velocities in the bubble swarm was achieved by adding fluorescing seed particles. Images of bubbles and fluorescing tracer particles were acquired by two CCD cameras. Hence, the images from tracers and bubbles were easily separated by optical interference filters with a bandwidth corresponding to the emitting wavelength of the fluorescing tracer particles and the wavelength of the applied Nd-YAG pulsed laser, respectively. To improve the phase separation of the system, the CCD cameras were additionally placed in a non-perpendicular arrangement with respect to the light sheet. The acquired images were evaluated with the minimum-quadratic-difference algorithm. The potential of this technique for the analysis of bubbly flows with higher void fraction was explored. In order to obtain averaged velocity maps of bubble and fluid within the entire column, about 1,000 image pairs were recorded and evaluated for each phase. In addition, turbulence intensities of the fluid were deduced from the measurements. The turbulence properties were used to characterise bubble-induced turbulence for various bubble mean diameters and gas hold-ups. Moreover, the determination of the average bubble slip velocity within the bubble swarm was possible.     

高速斜入水和水平入水气炮水箱实验系统

为开展模型高速斜入水和水中高速航行的水流场实验研究,研制了立式和卧式气炮与水箱组合的实验系统。通过快速阀和活塞阀控制气炮激发和驱动状态,一级气炮采用高压空气直接驱动弹托和模型,二级气炮采用高压空气驱动重活塞压缩使集气腔中产生高压气体,再驱动弹托和模型达到预定速度。通过调节水箱和发射管角度,使高速模型斜入水或水平入水。其中,立式可变发射角二级气炮可发射质量1～1000 g的模型至2500 m/s最大速度,卧式一级气炮可发射质量1～100 kg的模型至300 m/s最大速度。和小气室、高燃气压力火药驱动方式相比,新型气炮采用大体积、中低驱动压力气室,高压气体更接近等熵膨胀做功,调节高压气体压力,能较好地满足模型质量和速度的宽范围要求。结合光反射通断法测速、高速摄影和阴影流场显示等测量技术,得到立式气炮压缩管重活塞运动速度、压缩管末端压力时间曲线和模型倾斜与水平入水的流场阴影图像。结果表明:重活塞速度在膜片破裂前和理论计算值符合较好,但破膜后差异较大。立式气炮流场阴影图像反映了模型斜入水产生的空中和水中激波以及在气水界面的反射激波、空泡形成和侧向气水界面的破碎与飞溅等现象。从卧式气炮的模型水平入水阴影图像提取气泡轮廓,清楚地看出尾部气泡气水界面的波动和失稳。和商业计算软件Fluent计算结果相比,空泡上游区域基本重合,但尾流区域强湍流导致两者存在明显差异。和水洞实验相比,气炮水箱实验系统近真实地再现高速入水过程伴随的冲击和动态空化等物理现象和模型尺度效应。     

CFD modelling of most probable bubble nucleation rate from binary mixture with estimation of components’ mole fraction in critical cluster

The employment of different mathematical models to address specifically for the bubble nucleation rates of water vapour and dissolved air molecules is essential as the physics for them to form bubble nuclei is different. The available methods to calculate bubble nucleation rate in binary mixture such as density functional theory are complicated to be coupled along with computational fluid dynamics (CFD) approach. In addition, effect of dissolved gas concentration was neglected in most study for the prediction of bubble nucleation rates. The most probable bubble nucleation rate for the water vapour and dissolved air mixture in a 2D quasi-stable flow across a cavitating nozzle in current work was estimated via the statistical mean of all possible bubble nucleation rates of the mixture (different mole fractions of water vapour and dissolved air) and the corresponding number of molecules in critical cluster. Theoretically, the bubble nucleation rate is greatly dependent on components’ mole fraction in a critical cluster. Hence, the dissolved gas concentration effect was included in current work. Besides, the possible bubble nucleation rates were predicted based on the calculated number of molecules required to form a critical cluster. The estimation of components’ mole fraction in critical cluster for water vapour and dissolved air mixture was obtained by coupling the enhanced classical nucleation theory and CFD approach. In addition, the distribution of bubble nuclei of water vapour and dissolved air mixture could be predicted via the utilisation of population balance model.     

Full-field bubbly flow velocity measurements using a multiframe particle tracking technique

The study is an examination of two-phase dispersed air bubble flow about a cylindrical conductor emitting a constant heat flux. The technique of Particle Image Velocimetry is utilized in order to obtain a full-field non-invasive measurement of the resulting bubbly flow velocity field. The employed approach utilizes a flow visualization technique in which the instantaneous velocity profile of a given flow field is determined by digitally recording particle or bubble images within the flow over multiple successive video frames and then conducting a completely computational analysis of the data. The use of particle tracking algorithms which perform a point-by-point matching of seed images from one frame to the next allows construction of particle or bubble pathlines and instantaneous velocity field. Results were initially obtained for a synthetically created flow field and a single phase liquid convective field seeded with flow-following tracer particles. The method was additionally extended to measurements within a gas/liquid system in which bubble rise velocities over a substantial two-dimensional flow area were determined in order to demonstrate the effectiveness of the developed digital data acquisition and analysis methodology.A version of this paper was presented at the 12th Symposium on Turbulence, University of Missouri-Rolla, 24–26 September, 1990     

Heat transfer during bubble formation

The flow around a gas bubble forming at a submerged orifice without vaporisation is analysed both with and without heat transfer taking place from the bubble surface. The results of experiments are given in which high-speed photography was used to study air bubbling through four different sizes of orifice into water, for the isothermal and diabatic case. A numerical solution of the dimensionless equations using Hamming's modified predictor-corrector method gave excellent agreement with experimental measurements of the bubble volume-time history for both stages of growth as well as predicting accurately the transition point between the two stages. Most of the heat transfer was found to occur almost immediately after the start of bubble growth     

Three-dimensional CFD simulation of bubble–melt two-phase flow with air injecting and melt stirring

This paper reports on progress in developing CFD simulations of gas bubble–metallic melt turbulent flows induced by a pitched-blade impeller with an inclined shaft. Foaming process of aluminum foams, in which air is injected into molten aluminum composites and the melt is mechanical stirred by the impeller, has been investigated. A two-fluid model, incorporated with the Multiple Reference Frames (MRF) method is used to predict the three-dimensional gas–liquid flow in the foaming tank, in which a stirring shaft is positioned inclined into the melt. Locally average bubble size is also predicted by additively solving a transport equation for the bubble number density function, which accounts for effects of bubble breakup and coalescence phenomena. The computed bubble sizes are compared with experimental data from our water model measurement and reasonable agreements are obtained. Further, simulated results show that the volume averaged total and local gas fractions are generally increased with rising impeller speed and gas flow rate. The local averaged bubble size increases with increasing gas flow rate and orifice diameter and decreasing liquid viscosity, and decreases also with rising rotation speed of the impeller.     

The effect of aspect ratio in counter-current gas-liquid bubble columns: Experimental results and gas holdup correlations

It is generally admitted that the gas holdup is independent of the column dimensions and gas sparger design if three criteria are satisfied: the diameter of the bubble column is larger than 0.15 m, gas sparger openings are larger than 1–2 mm and the aspect ratio is larger than 5. This paper contributes to the existing discussion; in particular, the effect of the aspect ratio (within the range 1–15) in a counter-current gas-liquid bubble column has been experimentally studied and a new gas holdup correlation to estimate the influence of aspect ratio, operation mode and working fluid on the gas holdup has been proposed. The bubble column, equipped with a spider gas sparger, is 5.3 m in height, has an inner diameter of 0.24 m; gas superficial velocities in the range of 0.004–0.23 m/s have been considered, and, for the runs with water moving counter-currently to the gas phase, the liquid has been recirculated at a superficial velocity of −0.0846 m/s. Filtered air has been used as the gaseous phase in all the experiments, while the liquid phase has included tap water and different aqueous solutions of sodium chloride as electrolyte. Gas holdup measurements have been used to investigate the flow regime transitions and the global bubble column hydrodynamics. The counter-current mode has turned out to increase the gas holdup and destabilize the homogeneous flow regime; the presence of electrolytes has resulted in increasing the gas holdup and stabilizing the homogeneous flow regime; the aspect ratio, up to a critical value, has turned out to decrease the gas holdup and destabilize the homogeneous flow regime. The critical value of the aspect ratio ranged between 5 and 10, depending on the bubble column operation (i.e., batch or counter-current modes) and liquid phase properties. Since no correlation has been found in the literature that can correctly predict the gas holdup under the investigated conditions, a new scheme of gas holdup correlation has been proposed. Starting from considerations concerning the flow regime transition, corrective parameters are included in the gas holdup correlation to account for the effect of the changes introduced by the aspect ratio, operation mode and working fluid. The proposed correlation has been found to predict fairly well the present experimental data as well as previously published gas holdup data.     

Effect of glycerol addition on phase inversion in horizontal dispersed oil-water pipe flows

The effect of interfacial tension on the phase inversion process during horizontal pipe flow of an oil-aqueous solution was investigated. Interfacial tension was varied by adding small amounts of glycerol in the water phase. At these glycerol concentrations the density and viscosity of the aqueous phase changed by 1% or less. Exxsol™ D140 (5.5 mPa s, 828 kg m⁻³) was used as the oil phase. The experiments were carried out in a 38 mm ID acrylic test pipe. The phase continuity and appearance of phase inversion were investigated using conductivity (wire and ring) probes and an Electrical Resistance Tomographic (ERT) system. The ERT also provided diagrams of the phase distribution in a pipe cross section. Drop size distribution was monitored using a dual impedance probe. It was found that starting from a water continuous flow with increasing oil fraction at constant mixture velocity the mixture inverted initially in the middle of the pipe (measured at 19 mm from the top pipe wall) while a higher oil fraction was required for inversion at the top (measured at 4 mm from the top pipe wall) and finally the rest of the pipe. The addition of glycerol did not affect the phase fraction where the initial inversion occurred but caused an increase in the oil fraction needed to complete the inversion. The drop size measurements were used to explain this behaviour. Pressure drop was found to decrease with increasing oil fraction but this trend reversed when inversion spread to the pipe wall and the oil continuous phase came in contact with it.     

Analysis of the characteristic lengths in the bubble and slug flow regimes generated in a capillary T-junction

We present an analysis of the geometry of the continuous and disperse phases in the bubble and slug flow regimes in air–water mixtures generated in a capillary T-junction of 1  mm internal diameter. Bubble size dispersion is very low in the considered flow patterns. The concept of unit cell is used to identify two characteristic lengths of the two-phase flow, namely, the unit cell length and the bubble length. The relationship between these lengths and the gas and liquid superficial velocities, gas mean velocity, bubble generation frequency and volume average void fraction is analysed. We conclude that in the considered configuration the unit cell and bubble lengths can be predicted either by the ratio of the gas–liquid superficial velocities or the volume average void fraction.     

Liquid flow patterns about barbotage bubbles

This paper summarizes the results of a flow visualization study on the liquid motion around barbotage bubbles during growth and departure. Flow patterns, as well as for the first time, instantaneous velocities, are reported as a function of time and location about the bubbles. The experiments, employing the hydrogen-bubble technique and high-speed cine photography, were with: water as the liquid, air as the bubbled gas, orifice diameters of 0.116 and 0.252 cm, and different air flow rates; the two limiting cases of constant supply pressure and constant volumetric flow rate were covered. It was found that the liquid around a barbotage bubble assumes two velocity maxima, the first an outward maximum during bubble growth and the second in the opposite direction approximately at the time of bubble departure; further, liquid velocities were found to be higher close to the bubbling site. Certain differences in liquid velocities between the constant pressure and constant flow cases are explained in terms of available theoretical solutions to the bubble growth rate. Qualitative comparisons of the barbotage liquid flow patterns and those recently reported for boiling flow patterns are also presented.