Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Heat transfer measurements and correlations for air–water flow of different flow patterns in a horizontal pipe

Heat transfer coefficients were measured and new correlations were developed for two-phase, two-component (air and water) heat transfer in a horizontal pipe for different flow patterns. Flow patterns were observed in a transparent circular pipe using an air–water mixture. Visual identification of the flow patterns was supplemented with photographic data, and the results were plotted on the flow regime map proposed by Taitel and Dukler and agreed quite well with each other. A two-phase heat transfer experimental setup was built for this study and a total of 150 two-phase heat transfer data with different flow patterns were obtained under a uniform wall heat flux boundary condition. For these data, the superficial Reynolds number ranged from 640 to 35,500 for the liquid and from 540 to 21,200 for the gas. Our previously developed robust two-phase heat transfer correlation for a vertical pipe with modified constants predicted the horizontal pipe air–water heat transfer experimental data with very good accuracy. Overall the proposed correlations predicted the data with a mean deviation of 1.0% and an rms deviation of 12%.     

On the phase inversion process in an oil–water pipe flow

An experimental study of the phase inversion process in an oil–water flow through a pipe was carried out. Special attention was paid to the critical concentration of the dispersed phase at which phase inversion occurs and to the change in morphological structures during inversion. To that purpose two different types of experiments were performed: (1) continuous experiments during which the dispersed phase fraction was gradually increased and (2) direct experiments whereby oil and water were injected simultaneously into the pipe at a certain concentration. During the experiments detailed pictures were taken of the phase inversion process and simultaneously the electrical conductivity of the mixture was measured to determine which liquid formed the continuous phase and which the dispersed phase. Also the pressure gradient over several parts of the pipe was measured.     

Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m³) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe.     

Heat transfer to air–water annular flow in a horizontal pipe

Two-phase air–water flow and heat transfer in a 25 mm internal diameter horizontal pipe were investigated experimentally. The water superficial velocity varied from 24.2 m/s to 41.5 m/s and the air superficial velocity varied from 0.02 m/s to 0.09 m/s. The aim of the study was to determine the heat transfer coefficient and its connection to flow pattern and liquid film thickness. The flow patterns were visualized using a high speed video camera, and the film thickness was measured by the conductive tomography technique. The heat transfer coefficient was calculated from the temperature measurements using the infrared thermography method. It was found that the heat transfer coefficient at the bottom of the pipe is up to three times higher than that at the top, and becomes more uniform around the pipe for higher air flow-rates. Correlations on local and average Nusselt number were obtained and compared to results reported in the literature. The behavior of local heat transfer coefficient was analyzed and the role of film thickness and flow pattern was clarified.     

Effect of secondary flow on droplet distribution and deposition in horizontal annular pipe flow

In horizontal annular dispersed pipe flow the liquid film at the bottom is thicker and rougher than at the top of the pipe. A turbulent pipe flow experiencing a variation of roughness along the pipe wall will show a secondary flow. Such secondary flow, consisting of two counter-rotating cells in the cross-section of the tube, can change the distribution of the droplets inside the pipe and their deposition at the wall. Here, we compare the behaviour of the droplets (dispersed phase) with and without secondary flow, using large-eddy simulations. It is shown that the presence of secondary flow increases the droplet concentration in the core of the pipe and the droplet deposition-rate at the top of the pipe.     

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines. Friction factors play an important role in the accurate calculation of pressure drop. Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall, gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions. In this paper by modification of a friction model available in the literature, an improved semiempirical model is proposed. The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved. Pressure gradient and velocity profiles are validated against experimental data. Using the improved model, the pressure gradient deviation from experiments diminishes by about 3%; the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.     

On the liquid turbulence energy spectra in two-phase bubbly flow in a large diameter vertical pipe

The liquid turbulence structure of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. A dual optical probe was used to measure the bubble characteristics, while the liquid turbulence was measured using hot-film anemometry. Experiments were performed at two liquid superficial velocities of 0.2 and 0.68 m/s for gas superficial velocities in the range of 0–0.18 m/s, corresponding to an area averaged void fraction up to 13.6%. In general, there is an increase in the liquid turbulence energy when the bubbles are introduced into the liquid flow. The increase in the energy mainly occurs over a range of length scales that are on the order of the bubble diameter. A suppression of the turbulence was observed close to the wall at very low void fraction flows. Initially, the suppression occurs in the low wave number range and then extends to higher wave numbers as the suppression is increased.     

Measurements of cross-sectional instantaneous phase distribution in gas–liquid pipe flow

Two novel complementing methods that enable experimental study of gas and liquid phases distribution in two-phase pipe flow are considered. The first measuring technique uses a wire-mesh sensor that, in addition to providing data on instantaneous phase distribution in the pipe cross-section, also allows measuring instantaneous propagation velocities of the phase interface. A novel algorithm for processing the wire-mesh sensor data is suggested to determine the instantaneous boundaries of gas–liquid interface. The second method applied here takes advantage of the existence of sharp visible boundaries between the two phases. This optical instrument is based on a borescope that is connected to a digital video camera. Laser light sheet illumination makes it possible to obtain images in the illuminated pipe cross-section only. It is demonstrated that the wire-mesh-derived results based on application of the new algorithm improve the effective spatial resolution of the instrument and are in agreement with those obtained using the borescope. Advantages and limitations of both measuring techniques for the investigations of cross-sectional instantaneous phase distribution in two-phase pipe flows are discussed.     

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.     

On the use of the stratified momentum balance for the deduction of shear stress in horizontal gas–liquid pipe flow

The use of the stratified flow momentum balance for the deduction of interfacial and liquid wall shear stresses from experimental measurements is examined. A systematic error analysis is applied to the governing equations using the principle of maximum uncertainty. A series of air–water experiments were conducted in 50 and 80 mm diameter pipes, in which gas pressure drop, liquid height and gas wall shear stress were measured. A framework for the correlation of the deduced shear stresses is proposed from the experimental measurements. The uncertainty analysis is used to show that the definition of mean liquid height does not significantly influence the overall results. The development of empirical equations based on such methods may lead to total uncertainties of up to 40%, irrespective of accuracy of the experimental data or the appropriateness of the correlating technique. Comparisons with state-of-the-art correlations for the liquid wall and interfacial friction factor data showed even larger discrepancies between measurement and prediction.     

Two-phase flow distribution of air–water annular flow in a parallel flow heat exchanger

The air and water flow distribution are experimentally studied for a round header – flat tube geometry simulating a parallel flow heat exchanger. The number of branch flat tube is 30. The effects of tube outlet direction, tube protrusion depth as well as mass flux, and quality are investigated. The flow at the header inlet is identified as annular. For the downward flow configuration, the water flow distribution is significantly affected by the tube protrusion depth. For flush-mounted configuration, most of the water flows through frontal part of the header. As the protrusion depth increases, more water is forced to the rear part of the header. The effect of mass flux or quality is qualitatively the same as that of the protrusion depth. Increase of the mass flux or quality forces the water to rear part of the header. For the upward flow configuration, however, most of the water flows through rear part of the header. The protrusion depth, mass flux, or quality does not significantly alter the flow pattern. Possible explanations are provided based on the flow visualization results. Negligible difference on the water flow distribution was observed between the parallel and the reverse flow configuration.     

Friction factor improved correlations for laminar and turbulent gas–liquid flow in horizontal pipelines

We develop improved correlations for two-phase flow friction factor that consider the effect of the relative velocity of the phases, based on a database that includes 2560 gas–liquid flow experiments in horizontal pipes. The database includes a wide range of operational conditions and fluid properties for two-phase friction factor correlations. We classify the experiments by liquid holdup ranges to obtain composite analytical expressions for two-phase friction factor vs. the Reynolds number by fitting logistic dose curves to the experimental data with. We compute the liquid holdup values used to classify the experimental data using correlations proposed previously. The Reynolds number is based on the mixture velocity and the liquid kinematic viscosity. The Fanning friction factor for gas–liquid is defined in term of the mixture velocity and density. Additionally, we sort the experimental data by flow regime and obtain the two-phase friction factor improved correlations for dispersed bubble, slug, stratified and annular flow for different holdup ranges. We report error estimates for the predicted vs. measured friction factor together with standard deviation for each correlation. The accuracy of the correlations developed in this study is compared with that of other 21 correlations and models widely available in the specialized literature. Since different authors use different definitions for friction factors and Reynolds numbers, we present comparisons of the predicted pressure drop for each and every data point in the database. In most cases our correlations predict the pressure drop with much greater accuracy than those presented by previous authors.     

Numerical study on the two-phase flow distribution in a T-junction

The main objective of this paper is to investigate the ability of a two-dimensional two-fluid computer code to predict the phase separation in a T-junction. A new semi-implicit numerical scheme is developed for solving the two-fluid model equations. Special attention is directed to the modelling of the constitutive for the interfacial friction term. Detailed distribution of void fraction, pressure and velocities are obtained for an air–water mixture in a vertical tee. Good agreement was obtained between the computer code results and the experimental data for the phase separation in the T-junction.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

PIV measurement of particle motion in spiral gas–solid two-phase flow

With a concise review on some basic and novel algorithms and methods for the techniques of particle-imaging velocimetry (PIV), the paper reports an application of the PIV techniques to the investigation of particle motion in a gas–solid two-phase spiral flow in a horizontal tube. Axial velocities of the transported particles are obtained. Some important features of particle motion governing high transportation efficiency of the spiral flow are revealed by investigating probability density distribution of particle locations in a pipe cross-section.