The characteristics of gas/liquid flow in large risers at high pressures

Although most of the work reported on two-phase flows are limited to small pipe diameters, two-phase flow in large risers are increasingly being encountered in the petroleum and nuclear industries. In the present work, a wire mesh sensor was employed to obtain void fraction and bubble size distribution data and visualizations of steam/water flow in a large vertical pipe (194 mm in diameter) at 46 bar. For comparison purposes, measurements were made at similar phase velocities and physical properties to a dataset for nitrogen/naphtha flow in a similar-sized riser. There exist significant differences between both sets of data. Churn-turbulent flow is observed in the present work instead of slug flow, and this differs from the intermittent and semi-annular flow patterns reported for nitrogen/naphtha data. The mean void fraction of the nitrogen/naphtha data is higher than that of the present steam/water data due to the differences in purity in the liquid phases. Furthermore, core peak distributions are observed for the present work in contrast to the flatter profiles deduced for the nitrogen/naphtha using a power law relationship.     

Flow regime transition criteria for two-phase flow at reduced gravity conditions

Flow regime transition criteria are of practical importance for two-phase flow analyses at reduced gravity conditions. Here, flow regime transition criteria which take the frictional pressure loss effect into account were studied in detail. Criteria at reduced gravity conditions were developed by extending an existing model from normal gravity to reduced gravity conditions. A comparison of the newly developed flow regime transition criteria model with various experimental datasets taken at microgravity conditions showed satisfactory agreement. Sample computations of the model were performed at various gravity conditions, such as 0.196, 1.62, 3.71 and 9.81 m/s² corresponding to micro-gravity and lunar, Martian and Earth surface gravity, respectively. It was found that the effect of gravity on bubbly–slug and slug–annular (churn) transitions in a two-phase flow system was more pronounced at low liquid flow conditions, whereas the gravity effect could be ignored at high mixture volumetric flux conditions. While for the annular flow transitions due to flow reversal and onset of droplet entrainment, higher superficial gas velocity was obtained at higher gravity level.     

An investigation of pressure transients in pipelines with two-phase bubbly flow

This paper presents a two-dimensional model for the analysis of the pressure transient of a two-phase homogeneous bubbly mixture flowing in a pipeline and the numerical integration using the centre implicit method (CIM). Experiments were conducted to confirm the proposed sonic speed equation of an air–water mixture for an air concentration of less than 1%. The 2D CIM model is compared with the method of characteristics (MoC) for a two-phase bubbly flow in a pipeline. The comparisons show that the proposed 2D CIM model generally gives good agreement with the method of characteristics.     

Fluctuation characteristics of two-phase flow splitting at a vertical impacting T-junction

In order to investigate the fluctuation characteristics of two-phase flow splitting at a T-junction, particular attention was paid on Churn flow which had the strongest fluctuation comparing with bubble flow and annular flow. The main tube of the T-junction was vertical and the two branches were horizontal. All three pipes connecting to the junction were of 15 mm inner diameter. A statistical analysis based on Root Mean Square (RMS) was applied to temporal differential pressure signals and gas flow rate signals. The Power Spectral Density (PSD) was also employed to reveal their peculiar features in frequency domain as well. The effects of the extraction flow ratio and the gas and liquid superficial velocity upstream on fluctuation characteristics of gas-liquid two-phase flow splitting at the T-junction were investigated in detail. It is found that there is a wide fluctuation in both differential pressure and gas flow rate downstream at every extraction ratio (W₃/W₁) and the fluctuation intensity increases as W₃/W₁ increasing. It is also made clear that increasing either water superficial velocity or gas superficial velocity in inlet causes fluctuation to become more intensive.     

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.     

Phase split of nitrogen/non-Newtonian fluid two-phase flow at a micro-T-junction

An experimental investigation has been undertaken to understand the phase split of nitrogen gas/non-Newtonian liquid two-phase flow passing through a 0.5 mm T-junction that oriented horizontally. Four different liquids, including water and aqueous solutions of carboxymethyl cellulose (CMC) with different mass concentrations of 0.1, 0.2 and 0.3 wt%, were employed. Rheology experiments showed that different from water, CMC solutions in this study are pseudoplastic non-Newtonian fluid whose viscosity decreases with increasing the shear rate. The inlet flow patterns were observed to be slug flow, slug–annular flow and annular flow. The fraction of liquid taken off at the side arm for nitrogen gas/non-Newtonian liquid systems is found to be higher than that for nitrogen gas/Newtonian liquid systems in all inlet flow patterns. In addition, with increasing the pseudoplasticity of the liquid phase, the side arm liquid taken off increases, but the increasing degree varies with each flow pattern. For annular flow, the increasing degree is much greater than those for slug and slug–annular flows.     

Characteristics of the critical heat flux for downward flow in a vertical tube at low flow rate and low pressure conditions

The characteristics of the critical heat flux (CHF) for downward flow were studied experimentally with an Inconel 600 circular tube test section in a water test loop at low-flow rate (0 200 kg/m²s) and low-pressure (0.1 0.7 MPa) conditions. The attention was given to the effects of upstream conditions—upper plenum and inlet throttling. Two totally different kinds of CHF behaviors were observed. It seems appropriate to interpret them as flooding-type CHF and dryout in annular flow. The CHF in downward flow may vary from extremely unstable flow CHF as low as near the flooding CHF value to stable flow CHF as high as that of upflow, depending on the upstream conditions of the test section. The CHF correlation by Mishima and that by Weber were proposed for the presentation of the lower and upper limits of the CHF for downward flow in a vertical tube at low-flow rate and low-pressure conditions.     

Evaluation analysis of prediction methods for two-phase flow pressure drop in mini-channels

Two-thousand and ninety-two data of two-phase flow pressure drop were collected from 18 published papers of which the working fluids include R123, R134a, R22, R236ea, R245fa, R404a, R407C, R410a, R507, CO₂, water and air. The hydraulic diameter ranges from 0.506 to 12 mm; Re_l from 10 to 37,000, and Re_g from 3 to 4 × 10⁵. Eleven correlations and models for calculating the two-phase frictional pressure drop were evaluated based upon these data. The results show that the accuracy of the Lockhart–Martinelli method, Mishima and Hibiki correlation, Zhang and Mishima correlation and Lee and Mudawar correlation in the laminar region is very close to each other, while the Muller-Steinhagen and Heck correlation is the best among the evaluated correlations in the turbulent region. A modified Chisholm correlation was proposed, which is better than all of the evaluated correlations in the turbulent region and its mean relative error is about 29%. For refrigerants only, the new correlation and Muller-Steinhagen and Heck correlation are very close to each other and give better agreement than the other evaluated correlations.     

Flow-induced vibration in two-phase flow with wire coil inserts

Information of flow-induced vibration (FIV) in two-phase flow with wire coil inserts at atmospheric pressure, is presented in this study. FIV was measured in an upward vertical tube for four different wire coil inserts using an air–water mixture as process fluid. Vibration increased along with mass flux and quality. The narrower wire coils produced more vibration. The FIV prediction correlation for two-phase flow with wire coil inserts was experimentally developed with coefficient correlation value of 0.956.     

Surface wettability effect on flow pattern and pressure drop in adiabatic two-phase flows in rectangular microchannels with T-junction mixer

Wettability is an important parameter in micro-scale flow patterns. Previous research has usually been conducted in conventional microtubes due to limitations of visualizing flow patterns and fabricating microchannels. However, most microchannels in practical applications have rectangular shape. Furthermore, pressure drop is closely related with flow pattern. Hence, we studied water liquid and nitrogen gas flows in rectangular microchannels with different wettabilities. The rectangular glass microchannels were fabricated from photosensitive glass, whose surface is hydrophilic. The surface of one was silanized using octadecyl-trichloro-silane (OTS) to prepare a hydrophobic microchannel. The two-phase flow pattern was visualized with a high-speed camera and a long distance microscope. The frictional pressure drop in the microchannel was measured directly with embedded pressure ports. The flow pattern and pressure drop in the hydrophobic microchannel were totally different from those in the hydrophilic microchannel. Finally, the two-phase frictional pressure drop was analyzed based on the flow patterns of different wettabilities.     

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle flows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained. Received 14 June 1996 / Accepted 19 October 1996     

Recognition of gas–liquid two-phase flow patterns based on improved local binary pattern operator

A new method to pattern recognition of gas–liquid two-phase flow regimes based on improved local binary pattern (LBP) operator is proposed in this paper. Five statistic features are computed using the texture pattern matrix obtained from the improved LBP. The support vector machine and back-propagation neural network are trained to flow pattern recognition of five typical gas–liquid flow regimes. Experimental results demonstrate that the proposed method has achieved better recognition accuracy rates than others. It can provide reliable reference for other indirect measurement used to analyze flow patterns by its physical objectivity.     

Similarity in stationary motions of gas and two-phase medium with incompressible component

In this paper we suggest the transformation between the equations for a perfect gas and the equations describing in one-velocity approach the two-phase medium with any volume occupied by the incompressible phase. It is proved that the motion of a two-phase medium in the transformed coordinate system is similar with certain accuracy to that of a perfect gas. It means that the solutions obtained for perfect gas can be used to solve wave problems for media with incompressible component. There is no necessity directly to solve the problem for medium with incompressible component, and it is only sufficient to transform the known solution of the similar problem for a homogeneous medium. Thus, the solutions of many hydrodynamic problems for multi-component media with incompressible phase can be obtained without solving the original set of equations. The scope for the suggested transformation is demonstrated by reference to the strong explosion in a two-phase medium.     

Numerical calculations for two-phase jet flow of a mixture of pulverized-coal and gas

In this paper, Euler-Lagrange type equations are used to describe the jet flow of a mixture of pulverized-coal and gas, which is an unsteady axisymmetric two-phase flow. By means of the finite-difference method, the coal particle's distribution, velocity and trajectory in the flow field are obtained. The coal particles are represented by a finite number of computational particles. Each particle's diameter is randomly assigned according to a given distribution. The states of the computational particles are different from each other. Turbulence is accounted for in a stochastic model. Explicit time-splitting scheme is used to calculate the strongly coupling interphase term. The numerical results are reasonable. The comparison between the numerical results and the experiment data for the case of the oil droplet injection shows good agreement. This numerical technique can be extended to the calculation of other two-phase flows of dilute particles or a droplet system. Mr. Mei Renwei also participated in the work of this paper.     

Optimal systems of Lie subalgebras for a two-phase mass flow

We apply the Lie symmetry method to a two-phase mass flow model (Pudasaini, 2012 [18]) and construct one-, two- and three-dimensional optimal systems of Lie subalgebras corresponding to the non-linear PDEs. As an optimal system contains structurally important information about different types of invariant solutions, it provides precise insights into all possible invariant solutions emerging from infinitesimal symmetries. We use the optimal system of one-dimensional Lie subalgebras to reduce the two-phase mass flow model to other systems of PDEs. Using the fact that the Lie bracket contains information about further reduction, we further reduce to systems of ODEs and PDEs. We solve a system numerically and present results for different physical and Lie parameters. Simulations reveal fluid and solid dynamics are distinctly sensitive to different Lie parameters, whereas both phases are influenced by the solid and the fluid pressure parameters. Higher pressure gradients result in higher flow velocities and lower flow heights. Fluid velocities dominate solid velocities, but the solid heights are higher than the fluid heights. Results provide an overall picture of the physical process, and the coupled dynamics of the solid and fluid phase velocities and the flow heights. These are physically meaningful results in sheared inclined channel flow of coupled two-phase mixture. This confirms the consistency of the obtained similarity solutions and potential applicability of the models and the constructed optimal systems.     

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.     

An experimental investigation of roll waves in high pressure two-phase inclined pipe flows

This paper examines the effects of small upward inclinations on the formation of roll waves and the properties of fully developed roll waves at high pressure conditions. A total of 984 experiments were conducted at six positive pipe inclinations θ = 0.00°, 0.10°, 0.25°, 1.00°, 2.50° and 5.00° using a 25 m long 10 cm i.d. pipe. Sulfur hexafluoride (SF₆) was used at 8 bara giving a gas density of 50 kg/m³. Two independent mechanisms for the formation of roll waves were identified; (1) interaction between 2D shallow water waves and (2) a visible long wavelength instability of the stratified layer. Viscous long wavelength linear stability analysis predicted the critical liquid flow rate and liquid height for the initiation of roll waves when roll waves were formed due to the second mechanism. A simple equation from shallow water wave theory agreed with measurements for critical liquid flow rate when roll waves were formed due to the first mechanism. Shallow water wave speed agreed with critical wave speeds at transition and nonlinear wave speeds for fully developed roll waves in certain cases. The increase in interfacial friction due to the presence of large waves was compared with models from the literature.     

Rheological behaviour of two-phase flow of solid particles in a gas

The present investigation was concerned with the rheological behaviour of dilute suspensions of solid particles in a gas in a vertical cocurrent flow moving upwards. Starting from the experimentally determined dependence of the pressure drop on the concentration of solid particles and the Reynolds number of the carrier medium in the steady flow region, the rheological parameters were estimated using pseudo-shear diagrams. Air was the carrier medium and the dispersed phase was one of six fractions of polypropylene powder and five fractions of glass ballotini. The results show that the investigated two-phase systems have pseudoplastic character which becomes more pronounced with increases in concentration, equivalent diameter and density of solid particles in the flowing suspension. C _d coefficient of particle resistance - d _e equivalent diameter of particles - D column diameter - Fr Froude number - g gravitational acceleration - K rheological parameter - L length - n rheological parameter - p _t pressure drop due to friction - p _m total pressure drop - p _ag pressure drop due to acceleration of the gas phase - p _as pressure drop due to acceleration of the solid phase - p _g hydrostatic pressure of the gas phase - p _s specific effective weight of the dispersed phase - r radius - Re Reynolds number - Re _p Reynolds number of a particle - Re _G generalized Reynolds number - Re _G1 generalized Reynolds number relating to the end of the laminar flow region - Re _G2 generalized Reynolds number relating to the beginning of the turbulent flow region - w _z axial component of velocity - u _t steady free-fall velocity of a single particle - w average velocity - w _g average velocity of the gas phase - w _s average velocity of the dispersed phase of solid particles - relative mass fraction of solid particles - x _s volume fraction of solid particles - _g coefficient of pressure drop due to friction - µ dynamic viscosity - _g density of the gas phase - _m density of the suspension - _s density of solid particles - _ds density of the dispersed phase - _w shear stress at the wall