Two‐Phase Bubbly Flow Structure in Large‐Diameter Vertical Pipes

The two‐phase flow structure of an air‐water, bubbly, upward flow in a 20 cm diameter pipe is presented with particular emphasis on the local interfacial area concentration. The radial distribution of void fraction, bubble velocity, bubble size, bubble frequency, and interfacial area concentration were measured using a local dual‐optical probe. The experimental results showed that the saddle‐type distribution of void fraction and interfacial area concentration, which are common for bubbly flow in small diameter pipes, only appeared in the present experiments under conditions of very low area‐averaged void fraction (<?> < 0.04). The values for the interfacial area concentration were higher in large diameter pipes when compared with data obtained under the same flow conditions in small pipes. The area‐averaged void fraction data were correlated using the drift‐flux model.     

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.     

Experimental Investigation on the Holdup Distribution of Oil‐Water Two‐Phase Flow in Horizontal Parallel Tubes

An experimental investigation was conducted to study the holdup distribution of oil and water two‐phase flow in two parallel tubes with unequal tube diameter. Tests were performed using white oil (of viscosity 52 mPa s and density 860 kg/m³) and tap water as liquid phases at room temperature and atmospheric outlet pressure. Measurements were taken of water flow rates from 0.5 to 12.5 m³/h and input oil volume fractions from 3 to 94 %. Results showed that there were different flow pattern maps between the run and bypass tubes when oil‐water two‐phase flow is found in the parallel tubes. At low input fluid flow rates, a large deviation could be found on the average oil holdup between the bypass and the run tubes. However, with increased input oil fraction at constant water flow rate, the holdup at the bypass tube became close to that at the run tube. Furthermore, experimental data showed that there was no significant variation in flow pattern and holdup between the run and main tubes. In order to calculate the holdup in the form of segregated flow, the drift flux model has been used here.     

Influence of Gas Injection on the In‐Situ Oil Fraction of an Oil‐Water Flow in Horizontal Pipes

In this work, an experimental study was made on gas injection into an oil‐water flow in horizontal pipes with two unequal pipe diameters. Special attention was given to the influence of gas injection on the average in‐situ oil fraction. Measurements were made for input water flow rates of 1.25–5 m³/h, input oil flow rates of 0–8 m³/h and input gas flow rates of 0–9 m³/h. It was found that gas injection has a considerable influence on the in‐situ oil fraction. In general, a small increase in the rate of air injection leads to greatly decreasing in‐situ oil fractions. The in‐situ oil fraction with gas injection decreases to a greater extent than that without gas injection, at the same input liquid flow rates. At a given input water flow rate, the value of the in‐situ oil fraction in the pipe with the larger diameter is higher than that in the pipe with the smaller diameter. Furthermore, the drift flux models were extended to predict the average in‐situ fractions of the oil phase in the intermittent three‐phase flow regimes. A good agreement is obtained between theory and data, especially for the in‐situ oil fraction range of 0.2–1.0.     

Characteristics of air–water two‐phase flow in a 3‐mm smooth tube

Two‐phase flow pattern and friction characteristics for an air–water system in a 3.17 mm smooth tube are reported in this study. The range of mass flux is between 50 and 700 kg/m²s. The experimental data show that the two‐phase friction multipliers are strongly related to the flow pattern. For a stratified‐wavy flow pattern, a mass‐flux dependence of the two‐phase multipliers is seen. For a non‐stratified flow pattern, the two‐phase frictional multipliers are comparatively independent of mass flux. Correlations of the frictional multipliers are developed for stratified and non‐stratified flow. To use the appropriate correlation in different regime, a simple criterion is proposed.     

Equilibrium and non‐equilibrium gas–liquid two phase flow in long and short pipelines following a rupture

The two‐phase flow following the blowdown of pipeline carrying flashing liquid is numerically investigated by using thermodynamic equilibrium and non‐equilibrium models. Model equations are solved numerically by the finite volume method. The values of fluxes at cell boundaries are obtained by AUSM+‐up. To obtain proper values for the coefficients of dissipation, both single phase liquid and two phase shock tube problems are investigated. The transient release from the pressurized pipeline is studied for two cases of long and short pipes. Comparison of the predictions against experimental data reveals non‐equilibrium model performs a little better than equilibrium model in the prediction of temporal variations of pressure and void fraction of the long pipe. However, equilibrium model totally overestimates pressure and void fraction of the short pipe. The relative error of equilibrium model in the prediction of pressure variation with time exceeds 50% and it is 20% for non‐equilibrium model. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3214–3223, 2017     

Hydrodynamics of a Multi‐Stage Internal Loop Airlift Reactor

The local hydrodynamic properties in a multi‐stage internal loop airlift reactor were investigated in this study. The gas‐liquid two‐phase flow hydrodynamic properties, including gas holdup, bubble velocity, bubble diameter, and liquid circulation velocity at various stages were measured by dual electrical resistivity probes and conductivity cells. Detailed studies on the gas holdup, bubble velocity, bubble diameter, and liquid circulation velocity were conducted with respect to various values of superficial gas. The Zuber and Findlay drift flux model was used to represent the variation of slip velocity with total gas‐liquid velocity at various stages and the model fits the data well.     

Prediction of Layer‐Inversion Behavior of Down‐Flow Binary Solid‐Liquid Fluidized Beds

The layer‐inversion behavior of down‐flow binary solid‐liquid fluidized beds is predicted using the property‐averaging approach. The binary pair in this case consists of a larger solid species which is also heavier than its smaller counterpart, while both are lighter than the fluidizing medium. The model is based on using the generalized Richardson‐Zaki correlation for evaluation of the bed void fraction wherein mean values of particle properties are used. However, unlike the maximum bulk density condition for the conventional up‐flow binary solid fluidized bed, the model is based on a minimum bulk density condition for occurrence of layer inversion. This is due to the fact that the volume contraction phenomenon associated with the mixing of unequal solid species leads to a decrease in bulk density of the bed. Model predictions are also compared using the limited data available in the literature. Predictions are consistent with the observed mixing behavior.     

Influence of void fraction calculation on fidelity of CFD‐DEM simulation of gas‐solid bubbling fluidized beds

The correct calculation of cell void fraction is pivotal in accurate simulation of two‐phase flows using a computational fluid dynamics‐discrete element method (CFD‐DEM) approach. Two classical approaches for void fraction calculations (i.e., particle centroid method or PCM and analytical approach) were examined, and the accuracy of these methodologies in predicting the particle‐fluid flow characteristics of bubbling fluidized beds was investigated. It was found that there is a critical cell size (3.82 particle diameters) beyond which the PCM can achieve the same numerical stability and prediction accuracy as those of the analytical approach. There is also a critical cell size (1/19.3 domain size) below which meso‐scale flow structures are resolved. Moreover, a lower limit of cell size (1.63 particle diameters) was identified to satisfy the assumptions of CFD‐DEM governing equations. A reference map for selecting the ideal computational cell size and the suitable approach for void fraction calculation was subsequently developed. © 2014 American Institute of Chemical Engineers AIChE J, 60: 2000–2018, 2014     

Void fraction in vertical gas-liquid slug flow: Influence of liquid slug content

In this study we develop a model for computing the mean void fraction and the liquid slug void fraction in vertical upward gas-liquid intermittent flow. A new model for the rate of gas entrained from the Taylor bubble to the liquid slug is formulated. It uses the work done by the pressure force at the rear of the Taylor bubble. Then an iterative approach is employed for equating the gas entrainment flux and the gas flux obtained via conservation equations. Model predictions are compared with experimental data. The developed iterative method is found to provide reasonable quantitative predictions of the entrainment flux and of the void fraction at low and moderate liquid slug void fraction conditions. However, with an increased liquid slug void fraction experimental data indicate that the flow in the liquid slug transits to churn-heterogeneous bubbly flow thus gas entrainment flux tends to zero. Considering this effect in the iterative model significantly improved the predictions for large liquid slug void fraction conditions.     

Liquid holdup in concentric annuli during cocurrent gas‐liquid upflow

In the recent paper, an in‐depth investigation of liquid holdup during air‐water upflow through concentric annuli has been reported. The liquid holdup has been determined experimentally for the bubbly, slug and churn flow regimes. The drift flux model has been adopted for the theoretical estimation of holdup in the bubbly, dispersed bubbly and slug flow regimes. The pronounced effect of flow regime on this parameter as observed from experiments has been incorporated in the model by adopting different values of U₀, n and C₀. The asymmetry of the Taylor bubbles has been incorporated in the slug flow regime. The theoretical predictions exhibit a good agreement with the experimental data of the present work and that available in literature (Caetano et al., 1989b). The Hughmark's correlation is observed to correlate the churn flow data of the present work reasonably well.     

Characterization of Air‐Water Two‐Phase Vertical Flow by Using Electrical Resistance Imaging

The characterization of air‐water two‐phase vertical flow in a 12 m flow loop with 1.5 m of vertical section is studied by using electrical resistance tomography (ERT). By applying a fast data collection to a dual‐plane ERT sensor and an iterative image reconstruction algorithm, relevant information is gathered for implementation of flow characteristics, particularly for flow regime recognition. A cross‐correlation method is also used to interpret the velocity distribution of the gas phase on the cross section. The paper demonstrates that ERT can now be deployed routinely for velocity measurements and this capability will increase as faster measurement systems evolve.     

Druckabfall bei der Strömung von Gas/Dampf-Flüssigkeits-Gemischen in Rohren

Pressure drop during gas/vapour-liquid flow in pipes . Pressure drop in two phase flow is considered to be made up of terms for geodetical elevation or depression, acceleration, and friction. The geodetical and momentum pressure drop are discussed and reasonable correlations are presented, along with their limitations and range of application. Various relationships are available for calculating the technically important friction component. However, such predictions are still relatively inexact. If purely empirical correlations are neglected, all the remaining predictive schemes can be related to three basic physical models. These models as well as the commonly accepted and most reliable relationships are discussed and classified according to type. Furthermore, their scope is defined and the accuracy of prediction systematically compared with the aid of a newly set-up data bank. The extensive literature data consulted refer to single component two phase flow of water and various refrigerants and several two-component systems in horizontal and vertical unheated straight pipes under industrially relevant flow conditions. Finally, the accuracy of the prediction of some generally accepted void correlations is dealt with; here too, numerous published void fraction data have been gathered and checked.     

Frictional Pressure Drop of Air Non‐Newtonian Liquid Flow through Helical Coils in Horizontal Orientation

Investigations have been carried out to evaluate the two‐phase frictional pressure drop for air non‐Newtonian liquid flow through helical coils in horizontal orientation. The experiments performed using 36 different helical coils and 4 different concentrations of sodium salt of carboxymethyl—cellulose (SCMC) as non‐Newtonian liquids. The effects of air and liquid flow rate, coil diameter, helix angle and liquid properties‐ on two‐phase frictional pressure drop have been discussed. An attempt has been made to fit the experimental two‐phase frictional pressure drop data by the Lockhart and Martinelli, Chem. Eng. Prog. 45 , 39–48 (1949) correlation and the modified Lockhart‐Martinelli correlation as presented by different authors. In another approach, friction factor method was adopted to correlate the experimental data by dimensional analysis. The correlation developed predicts the two‐phase frictional pressure drop with acceptable statistical accuracy.     

Predicting the Viscosity of Non‐Newtonian Water‐in‐Oil Emulsions Based on the Lederer Model

The accurate prediction of the viscosity of emulsions is highly important for oil well exploitation. Commonly used models for predicting the viscosity of water‐in‐oil (W/O) emulsions composed by two or three factors cannot always fit well the viscosity of W/O emulsions, especially in the case of non‐Newtonian W/O emulsions. An innovative and comprehensive method for predicting the viscosity of such emulsions was developed based on the Lederer, Arrhenius, and Einstein models, using experimental data. Compared with the commonly applied W/O emulsion viscosity models, the proposed method considers more factors, including temperature, volume fraction of water, shear rate, and viscosity of the continuous (oil) and dispersed phase (water). Numerous published data points were collected from the literature to verify the accuracy and reliability of the method. The calculation results prove the high accuracy of the model.     

Atomization of Liquids by Two‐phase Gas‐liquid Flow through a Plain‐orifice Nozzle: Flow Regimes inside the Nozzle

Liquids or suspensions are divided into sprays of small droplets by atomization of two‐phase gas‐liquid mixtures. In this way either an equal distribution of the droplets or the generation of large surface areas of the liquid phase are accomplished, leading to increased heat‐ and mass‐transfer. The spatial and time dependency of the mean droplet diameter is a function of the total pressure upstream of the nozzle, the volumetric flow rate of the liquid and the gas, as well as on the flow regime in the nozzle. Thus the radial and axial profile of the void fraction inside the nozzle are measured with an electrical measurement technique. In addition, the flow in the nozzle is imaged by a high‐speed camera. Three flow regimes are identified. These are bubbly flow, plug flow and annular flow. A continuous flow of the emitting spray is observed for bubbly flow and annular flow only. The distribution of the dispersed bubble phase is given by ratio of the isothermic compression energy needed to pressurize the gas mass flow rate from atmospheric pressure up to the total pressure in front of the nozzle, and the potential energy of the supplied liquid mass flow rate.     

垂直圆管内液氮流动沸腾的理论模型及数值模拟

分析了液氮流动沸腾过程中气液两相间动量、能量以及质量的传输规律,建立了相应的理论模型,新模型重点修正了界面面积浓度和气泡挣脱直径的计算式;采用新建立的理论模型作为封闭方程对CFX-4.3中内建的双流体模型进行了修正,并采用修正后的双流体模型模拟了液氮在垂直圆管内的流动沸腾过程.数值模拟的结果与文献中的实验数据吻合较好,证明了本文所建模型的合理性.通过数值模拟发现,两相流参数分布的不均匀性对液氮流动沸腾过程中的热质传输特性有重要影响.     

Drag reduction phenomenon in viscous oil‐water dispersed pipe flow: Experimental investigation and phenomenological modeling

An experimental study on drag‐reduction phenomenon in dispersed oil‐water flow has been performed in a 26‐mm‐i.d. Twelve meter long horizontal glass pipe. The flow was characterized using a novel wire‐mesh sensor based on capacitance measurements and high‐speed video recording. New two‐phase pressure gradient, volume fraction, and phase distribution data have been used in the analysis. Drag reduction and slip ratio were detected at oil volume fractions between 10 and 45% and high mixture Reynolds numbers, and with water as the dominant phase. Phase‐fraction distribution diagrams and cross‐sectional imaging of the flow suggested the presence of a higher amount of water near to the pipe wall. Based on that, a phenomenology for explaining drag reduction in dispersed flow in a flow situation where slip ratio is significant is proposed. A simple phenomenological model is developed and the agreement between model predictions and data, including data from the literature, is encouraging. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Voidage characteristics and prediction of bed expansion in liquid-solid inverse fluidized bed

Studies on voidage fluctuations, axial voidage profile and bed expansion are carried out by measuring the local void fraction using particles of wide ranging characteristics in liquid-solid inverse fluidized bed. The quality of fluidization is elucidated by the local voidage fluctuations. The RMS voidage fluctuation depicts a maximum with respect to average bed void fraction and increases with increase in Archimedes number. The fluidization quality has been quantified using average normalized RMS voidage fluctuation in terms of Transition number. The axial void fraction is almost uniform throughout the bed except for particles with size distribution. All the literature and present experimental data on bed expansion are unified in terms of Richardson and Zaki equation using experimental terminal velocities. A new correlation is proposed for predicting the wall effect corrected experimental terminal velocities, as a substitute for standard drag equation. The bed expansion data are also predicted using the drift flux model.