Structures in gas–liquid churn flow in a large diameter vertical pipe

Gas–Liquid two phase co-current flow in a vertical riser with an internal diameter of 127 mm was investigated in the churn flow pattern. This paper presents detailed experimental data obtained using a Wire Mesh Sensor. It shows that the most obvious features of the flow are huge waves travelling on the liquid film. Wisps, large tendrils of liquid and the product of incomplete atomisation, which had previously detected in smaller diameter pipes, have also been found in the larger diameter pipe employed here. The output of the Wire Mesh Sensor has been used to determine the overall void fraction. When examined within a drift flux framework, it shows a distribution coefficient of ∼1, in contrast to data for lower gas flow rates. Film thickness time series extracted from the Wire Mesh Sensor output have been examined and the trends of mean film thickness, that of the base film and the wave peaks are presented and discussed. The occurrence of wisps and their frequencies have been quantified.     

Prediction of the transition from stratified to slug flow or roll-waves in gas–liquid horizontal pipes

In stratified gas–liquid horizontal pipe flow, growing long wavelength waves may reach the top of the pipe and form a slug flow, or evolve into roll-waves. At certain flow conditions, slugs may grow to become extremely long, e.g. 500 pipe diameter. The existence of long slugs may cause operational upsets and a reduction in the flow efficiency. Therefore, predicting the flow conditions at which the long slugs appear contributes to a better design and management of the flow to maximize the flow efficiency.     

Interfacial area transport across a 90° vertical-upward elbow in air–water bubbly two-phase flow

This study develops a one-group interfacial area transport equation (IATE) for vertical-upward-to-horizontal air–water bubbly two-phase flows through a 90° elbow with a non-dimensional centerline radius of curvature of three. In order to develop the model, an extensive database is established by acquiring local two-phase flow parameters using a four-sensor conductivity probe upstream and downstream of the elbow. The data show there exist three characteristic regions in void distribution, including a bimodal-to-bimodal region, a bimodal-to-single-peaked region, and a developed horizontal flow region with void accumulated at the top of the pipe cross-section. Using the database, the preliminary dissipation length model developed by Yadav et al. (2014b) is improved by including the transition region near the exit of the elbow in addition to the dissipation region. To close the IATE model, the bubble velocity advection term and bubble interaction terms in the IATE are correlated with the parameter characterizing the “elbow-strength”. The two-phase pressure drop across the elbow is modeled using the modified Lockhart–Martinelli correlation which takes into account the minor loss effect. The closed IATE model is implemented to predict interfacial area transport in vertical-upward-to-horizontal two-phase flow. It is found that the developed model is capable of predicting interfacial area concentration with an average percent difference of less than ±6%.     

Friction factor correlations for gas–liquid/liquid–liquid flows with drag-reducing polymers in horizontal pipes

In the present paper, two correlations have been developed to predict the effect of DRP on friction factor of two-phase flow for any pipe diameter. The correlations have been verified using published experimental results of DRP added to air–liquid annular flow and for oil–water flows with any flow pattern at the asymptotic state (maximum drag reduction). Such correlations are not available in literature and considered to be very useful in predicting the drag reduction using DRP and in understanding the most significant parameters that affecting the mechanism of drag reduction by polymers in multiphase flow.     

A growth model for dynamic slugs in gas–liquid horizontal pipes

Long liquid slugs, with sizes reaching 500 pipe diameters or more, may form in gas–liquid horizontal pipe flow at intermediate liquid loadings. Such slugs cause serious operational upsets due to the strong fluctuations in flow supply and pressure. Therefore, predicting the transition from short (hydrodynamic) to long slug flow regimes may play a significant role in preventing or reducing the negative effects caused by the long slugs.     

A compressible two-layer model for transient gas–liquid flows in pipes

This work is dedicated to the modeling of gas–liquid flows in pipes. As a first step, a new two-layer model is proposed to deal with the stratified regime. The starting point is the isentropic Euler set of equations for each phase where the classical hydrostatic assumption is made for the liquid. The main difference with the models issued from the classical literature is that the liquid as well as the gas is assumed compressible. In that framework, an averaging process results in a five-equation system where the hydrostatic constraint has been used to define the interfacial pressure. Closure laws for the interfacial velocity and source terms such as mass and momentum transfer are provided following an entropy inequality. The resulting model is hyperbolic with non-conservative terms. Therefore, regarding the homogeneous part of the system, the definition and uniqueness of jump conditions is studied carefully and acquired. The nature of characteristic fields and the corresponding Riemann invariants are also detailed. Thus, one may build analytical solutions for the Riemann problem. In addition, positivity is obtained for heights and densities. The overall derivation deals with gas–liquid flows through rectangular channels, circular pipes with variable cross section and includes vapor–liquid flows.     

Review on vertical gas–liquid slug flow

Vertical slug flow is characterized by the rise of long bullet-shaped gas bubbles with a diameter almost matching that of the tube - Taylor bubbles. Liquid slugs separate consecutive Taylor bubbles, which may interact and coalesce if the distance between them is small. Slug flow has numerous industrial applications, being also observed on physiological and geological systems. In spite of the contribution of the development of non-intrusive experimental techniques to a deeper understanding of slug flow features, the complexity of this flow pattern requires the combined use of numerical approaches to overcome some of the optical problems reported in experimental methods, and other limitations related to the flow aperiodic behavior.The need to systematize the large amount of data published on the subject and to understand the limitations of the techniques employed constitutes the motivation for this review. In the present work, literature on vertical gas–liquid slug flow, with Newtonian fluids, from 1943 to 2015, covering theoretical, experimental and numerical approaches, is reviewed. Focus is given to single and trains of Taylor bubbles rising through stagnant and co-current liquids.It should be emphasized, however, that further research still needs to be conducted in some particular areas, namely the hydrodynamics of the liquid film surrounding the Taylor bubbles, the interaction between consecutive bubbles, and a more detailed approach to the flow of Taylor bubbles through co-current liquids.     

Algebraic turbulence modeling in adiabatic gas–liquid annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic gas–liquid annular two-phase flow. After reviewing the existing literature, two new algebraic turbulence models are proposed for both the liquid film and the droplet laden gas core of annular two-phase flow. Both turbulence models are calibrated with experimental data taken from the open literature and their performance critically assessed. Although the proposed turbulence models reproduce the key parameters of annular flow well (average liquid film thickness and pressure gradient) and the predicted velocity profiles for the core flow compare favorably with available core flow velocity measurements, a more accurate experimental database is required to further improve the models accuracy and range of applicability.     

Prediction of the entrained liquid fraction in vertical annular gas–liquid two-phase flow

The study considers the prediction of the entrained liquid fraction in adiabatic gas–liquid annular two-phase flow in vertical pipes. Nine empirical correlations have been tested against an experimental data bank drawn together in this study containing 1504 points for 8 different gas–liquid combinations and 19 different tube diameters from 5.00 mm to 57.1 mm. The correlation of Sawant, Ishii and Mishima and the one of Oliemans, Pots and Trompé were found to best reproduce the available data. A new correlating approach, derived from both physical intuition and dimensional analysis and capable of providing further physical insight into the liquid film atomization process, was proposed and worked better than any of the existing methods. This new correlation is based on the core flow Weber number that is also a controlling dimensionless group in determining the wall shear stress and associated frictional pressure gradient of annular flows.     

Residual mass and flow regimes for the immiscible liquid–liquid displacement in a plane channel

The motion of two immiscible liquids in a plane channel is analyzed for the case in which the flow conditions and the interactions between the liquids and the solid surface maintain the displaced fluid attached to the wall. The Galerkin Finite Element Method is used to compute the velocity field and the configuration of the interface between the two fluids. We compare the residual mass fraction left on the wall with its two counterparts in capillary tubes, namely residual mass fraction and dimensionless layer thickness of the displaced fluid. The main result of this comparison was that although there is a qualitative similarity concerning the layer thickness between the two cases, the residual fraction of mass presented an important difference, showing that when the aspect ratio of the capillary passage is large there is an increase in the displacement efficiency. The thickness of the displaced liquid film attached to the channel walls is a function of the capillary number (Ca) and the viscosity ratio (N_μ). A map of streamlines in the Cartesian space (Ca, N_μ) with the different flow regimes of the problem is presented. We also showed that we can adapt the available analytical results obtained for gas-displacement in capillary tubes to the plane channel case, for low values of Ca.     

Prediction of mono-disperse gas–liquid turbulent flow in a vertical pipe

A two-fluid model in the Eulerian–Eulerian framework has been implemented for the prediction of gas volume fraction, mean phasic velocities, and the liquid phase turbulence properties for gas–liquid upward flow in a vertical pipe. The governing two-fluid transport equations are discretized using the finite volume method and a low Reynolds number $k-\varepsilon$ model is used to predict the turbulence field for the continuous liquid phase. In the present analysis, a fully developed one-dimensional flow is considered where the gas volume fraction profile is predicted using the radial force balance for the bubble phase. The current study investigates: (1) the turbulence modulation terms which represent the effect of bubbles on the liquid phase turbulence in the k–ε transport equations; (2) the role of the bubble induced turbulent viscosity compared to turbulence generated by shear; and (3) the effect of bubble size on the radial forces which results in either a center-peak or a wall-peak in the gas volume fraction profiles. The results obtained from the current simulation are generally in good agreement with the experimental data, and somewhat improved over the predictions of some previous numerical studies.     

Prediction of void fraction for gas–liquid flow in horizontal,upward and downward inclined pipes using artificial neural network

In this work, the ability of artificial neural networks (ANNs) to predict void fraction of gas–liquid two–phase flow in horizontal and inclined pipes was investigated. For this purpose, an ANN model was designed and trained using a total of 301 experimental data points reported in the literature for inclination angles between –20° and +20°. Pipe inclination angle as well as superficial Reynolds number of gas (Re_sg) and liquid (Re_sl) were chosen as input parameters of different structures of multilayer perceptron (MLP) neural networks, while the corresponding void fraction was selected as their output parameter. A hyperbolic tangent sigmoid and a linear function were employed as transfer functions of hidden and output layers, respectively, and Levenberg–Marquardt back propagation algorithm was used to train the networks. By trial–and–error method, a three–layer network with 10 neurons in the hidden layer was achieved as optimal structure of the ANN which made it possible to predict the void fraction with a high accuracy. Mean absolute percent error (MAPE) of 1.81% and coefficient of determination (R²) of 0.9976 for training data and MAPE of 1.52% and R² value of 0.9948 for testing data were obtained. Also for all data, MAPE of 1.95% and R² value of 0.9972 were calculated, and 96% data were within ±5% error band. In addition, the accuracy of the proposed ANN model was compared with the predictions from 17 void fraction correlations available in the literature for different flow patterns and horizontal and inclined flows. For all cases, the proposed ANN model gave better performance than all of the studied correlations. The results confirm the very good capability of the ANNs to predict the void fractions of gas–liquid flow in inclined pipes, regardless of flow pattern. Finally, by performing interpolation using the trained network, the void fraction values for some other conditions were predicted.     

Experimental flow pattern map,slippage and time–frequency representation of oil–water two-phase flow in horizontal small diameter pipes

We detect the flow structures of a horizontal oil–water two-phase flow in a 20 mm inner-diameter pipe using 8-channels radial mini-conductance probes. In particular, we present an experimental flow pattern map that includes 218 flow conditions and compare this map to the flow pattern transitional boundaries predicted by published models. In addition, using the Adaptive Optimal Kernel Time–Frequency Representation, we analyze the conductance fluctuating signals and characterize the flow pattern in terms of the total energy and dominant frequency. Based on the liquid holdup measurements using the quickly closing valve technology combined with three parallel-wire capacitance probes, we investigate the slip effect between the oil and water phases under various flow conditions. The results show that the flow structures greatly affect the slippage, and the slip ratio is sensitive to flow pattern variations.     

The flow pattern map of a two-phase non-Newtonian liquid–gas flow in the vertical pipe

A map for the determination of flow pattern for two-phase flow of gas and non-Newtonian liquid in the vertical pipe has been presented. Our own experimental data confirm applicability of such a map.     

A theoretical model for gas–liquid slug flow in down inclined ducts

Analytical solutions are obtained for flows in downwardly inclined ducts, partly filled by a liquid and containing finite amplitude moving jumps. A unified theory for both roll waves and periodic slug flows in rounded ducts of arbitrary cross-section is worked out by means of some simplifications. The article is focused on slugs: a set of equations is obtained, which predicts the transition between roll waves and slug regimes and gives access to all flow characteristics without any need of closure laws concerning either the speed of propagation or the slug length. As a result, we gain a new insight on the physical structure of slug flow. The proposed model is valid for sufficient inclination, small pressure gradient along the duct and negligible superficial tension. Owing to assumptions, only main trends and orders of magnitude observed in experiments are to be checked. In this connection the model fits most of the previously published experimental results obtained in ducts of circular cross-section: the domain of occurrence of downwardly propagating slugs is satisfactorily predicted, the limitations in drift velocity and in liquid layer thickness are demonstrated and upwardly propagating slugs are possible.     

Analysis of spatial and temporal spectra of liquid film surface in annular gas–liquid flow

Wavy structure of liquid film in annular gas–liquid flow without liquid entrainment consists of fast long-living primary waves and slow short-living secondary waves. In present paper, results of spectral analysis of this wavy structure are presented. Application of high-speed LIF technique allowed us to perform such analysis in both spatial and temporal domains. Power spectra in both domains are characterized by one-humped shape with long exponential tail. Influence of gas velocity, liquid Reynolds number, liquid viscosity and pipe diameter on frequency of the waves is investigated. When gravity effect is much lesser than the shear stress, similarity of power spectra at different gas velocities is observed. Using combination of spectral analysis and identification of characteristic lines of primary waves, frequency of generation of secondary waves by primary waves is measured.     

Characteristics of flow and liquid distribution in a gas–liquid vortex separator with multi spiral arms

Fischer–Tropsch (F–T) synthesis is an important route to achieve the clean fuel production. The performance of gas–liquid separation equipment involving in the progressive condensation and separation of light and heavy hydrocarbons in the oil-gas products has become a bottleneck restricting the smooth operation of the F–T process. In order to remove the bottleneck, a gas–liquid vortex separator with simple structure, low pressure drop and big separation capacity was designed to achieve the efficient separation between gas and droplets for a long period. The RSM (Reynolds Stress Model) and DPM (Discrete Phase Method) are employed to simulate the flow characteristics and liquid distribution in the separator. The results show that the separation efficiency is influenced by the flow field and liquid phase concentration in the annular zone. The transverse vortex at the top of spiral arm entrains the droplets with small diameter into the upper annular zone. The entrained droplets rotate upward at an angle of about 37.4°. The screw pitch between neighbor liquid threads is about 0.3 m. There is a top liquid ring in the top of annular zone, where the higher is the liquid phase concentration, the lower is the separation efficiency. It is found that by changing the operating condition and the annular zone height the vortex can be strengthened but not enlarged by the inlet velocity. The screw pitch is not affected by both inlet velocity and annular zone height. The liquid phase concentration in the top liquid ring decreases with both the increases of inlet velocity and annular zone height. The total pressure drop is almost not affected by the annular zone height but is obviously affected by the inlet velocity. When the height of annular zone is more than 940 mm, the separation efficiency is not changed. Therefore, the annular zone height of 940 mm is thought to be the most economical design.     

Investigation of two-phase flow pattern,liquid holdup and pressure drop in viscous oil–gas flow

An experimental investigation was carried out on viscous oil–gas flow characteristics in a 69 mm internal diameter pipe. Two-phase flow patterns were determined from holdup time-traces and videos of the flow field in a transparent section of the pipe, in which synthetic commercial oils (32 and 100 cP) and sulfur hexafluoride gas (SF₆) were fed at oil superficial velocities from 0.04 to 3 m/s and gas superficial velocities from 0.0075 to 3 m/s.