Improved Void Fraction Correlation for Two‐Phase Flow in Large‐Diameter Annuli

A flow pattern‐independent void fraction correlation for gas‐liquid two‐phase flow in vertical large‐diameter annuli is established. Two equations are proposed for the parameters of a drift flux model‐based correlation: the distribution parameter and the drift flux velocity. These equations are expressed as a function of two‐phase flow variables including void fraction, fluid properties, pipe geometry, and phase flow rates. Experiments were performed to study the void fraction of vertical air‐water two‐phase flow in large‐diameter annuli. The obtained experimental data along with the literature data of Caetano are used to verify the performance of the proposed void fraction correlation. The accuracy of this correlation is compared with nineteen frequently used correlations in literature. The proposed correlation was found to predict the void fraction consistently with a better accuracy.     

Phase Distribution Measurements during Transient Bubbly Two‐Phase Flow in a Vertical Pipe

An experimental procedure for investigating transient bubble flow for an adiabatic air/water system with vertical upward flow in a pipe is presented. The results of the measured local transient two‐phase flow parameters are shown along a pipe length of approx. four meters. From the measured radial phase distributions under steady and under transient conditions one can draw conclusions about the interfacial forces. Here, the effects indicate the action of forces such as a transverse lift force and a time dependent force like the virtual mass force during the transient. For modelling the transverse lift force which seems to play a dominant role for that flow regime the formulation of Zun was chosen and it was implemented into the commercial CFD‐Code Fluent Release 4.4.4 via user‐defined subroutines. Finally, results from the simulation of the steady states of start and end conditions of an experimental measured transient are shown.     

Void Fraction Measurement for Two‐Phase Flow Using Electrical Resistance Tomography

Measurement of void fraction of two‐phase flows remains a challenging area. In this paper the application of an electrical resistance tomography (ERT) system for this purpose has been studied. A new approach through the direct use of the voltage data measured by the ERT system is presented. The measured voltage data are first compressed through a feature extraction, and a polynomial regression procedure is followed to obtain the relationship between the void fraction and the feature extracted. Both simulation and experiment are carried out to verify the approach. The methodology of the new approach, simulation and experimental results are presented in the paper.     

Modeling of Two‐Phase Flashing Flow of Multicomponent Mixtures in Large Diameter Pipes

A two‐phase flashing flow model is developed to predict the distributions of pressure, temperature, velocity and evaporation rate in a transfer line, which is a typical example of a two‐phase flow pipe in the petrochemical industry. The model is proposed based on the pressure drop model and the multi‐stage flash model. The results indicate that pressure drop, temperature drop, and change of evaporation rate mainly occur in the transition section and the junction site of the transfer line. The predictions of the model have been tested with reliable field data and the good agreement obtained may lead to a better understanding of the two‐phase flashing flow phenomenon, as well as demonstrating the feasibility of applying the model into the design and optimization of pipelines.     

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.     

Large‐Eddy Simulation of Oil‐Water Annular Flow in Eccentric Vertical Pipes

One of the difficulties related to oil exploration is transporting heavy oil since its high viscosity causes high‐pressure drop and energy consumption. In order to save energy, the core annular flow (CAF) can be applied where a two‐phase annular flow occurs, with peripheral water flowing offering a reduction in energy expenditure. The multiphase flow was studied experimentally in a simple purpose‐built unit. To theoretically handle the CAF, computational fluid dynamics simulations were done with the commercial package Ansys Fluent. The flow was considered turbulent, isothermal, incompressible, and 3D, and both stationary and transient cases were evaluated. The volume‐of‐fluid model was adopted for the multiphase system, and water/oil interface and turbulence phenomena were well predicted.     

垂直管内气液两相弹状流含气率研究

In order to investigate the influence of the entrance effect on the spatial distribution of phases, the experiments on gas-liquid two-phase slug flow in a vertical pipe of 0.03m ID were carried out by using optical probes and an EKTAPRO 1000 high speed motion analyzer. It demonstrates that the radial profile of slug flow void fraction is parabolic. Influenced by the falling liquid film, the radial profile curve of liquid slug void fraction in the wake region is also parabolic. Since fully turbulent velocity distribution is built up in the developed region,the void fraction profile in this region is the saddle type. At given superficial liquid velocity, the liquid slug void fraction increases with gas velocity. The radial profiles of liquid slug void fraction at different axial locations are all saddle curves, but void fraction is obviously high around the centerline in the entrance region. The nearer the measuring station is from the entrance, the farther the peak location is away from the wall.     

Drift Flux Modelling for a Two‐Phase System in a Flotation Column

Using different two‐phase systems, a laboratory flotation column was operated over a wide range of bubble size from 250 to 1100 µm. The drift flux analysis for the tests was used to obtain a new correlation to relate the characteristic exponent m in the drift flux equation j_gf = U_tα_g(1 — α_g)^m, to bubble Reynolds number, where j_gf is the drift flux, U_t is the terminal velocity of an individual bubble, α_g is the gas volume fraction, and m is an exponent that depends on flow conditions. m = 20.26 + 1.89 Re_b/4.38 + Re_b.     

竖直大圆管内界面面积浓度分布特性

采用光纤探针测量方法对垂直上升管中空气-水两相流动的界面面积浓度（IAC）分布规律进行了研究。实验选用圆管直径为100mm和50mm(分别属于大管径和过渡管径范围)。气相、液相表观速度的范围分别为 0.01～1.0 m?s-1和 0～1.0 m?s-1。首先通过高速摄影获取的影像数据对光纤探针法的测量精度进行评价,并在此基础上对光纤测量方法进行了标定,然后利用实验获得了IAC、索特平均直径（气泡直径）径向分布情况。通过对比分析两个圆管管内IAC、气泡直径分布特点,得到了管径的尺度效应对IAC及气泡行为的影响规律。利用获得的IAC数据,对几类IAC计算模型进行了评价,并基于Hibiki-Ishii（2002）模型,通过实验数据拟合,得到了更高计算精度的IAC计算关联式。     

Flow pattern transition in gas‐liquid downflow through narrow vertical tubes

The present report studies on the flow pattern transitions during vertical air water downflow through millichannels (0.83 ≤ Eötvös no. ≤ 20.63). Four basic flow patterns namely falling film flow, slug flow, bubbly flow, and annular flow are observed in the range of experimental conditions studied and their range of existence has been noted to vary with tube diameter and phase velocities. Based on experimental observations, phenomenological models are proposed to predict the transition boundaries between adjacent patterns. These have been validated with experimental flow pattern maps from the present experiments. Thus the study formalizes procedure for developing a generalized flow pattern map for gas‐liquid downflow in narrow tubes. © 2016 American Institute of Chemical Engineers AIChE J, 63: 792–800, 2017     

Hydrodynamic Modelling of Internal Loop Airlift Reactor Applying Drift‐Flux Model in Bubbly Flow Regime

A new model for the liquid circulation rates in airlift reactor (ALR) is presented. The model is based on the energy balance for the flow loop (riser, turn riser‐downcomer, downcomer, and turn downcomer‐riser) coupled with a drift flux theory of two‐phase flow gas‐liquid system, considering a bubbly flow regime. The predicted values of the liquid circulation rates by the developed model are compared with experimental results performed in a 22 dm³ internal loop airlift reactor and with the results obtained in the literatures. The proposed model predicted the experimental results very well. Slip velocity relationship based on the drift flux model was proposed; including the gas holdup, bubble size and the liquid physical properties. The predicted slip velocity was similar to that obtained from the literature. The study revealed that appropriate arrangements of internal bioreactor parts can positively influence the liquid circulation velocity at the same energy consumption. The proposed models are useful in the design; scale up and characterization of the internal loop airlift reactors, and provides a direct method of predicting hydrodynamic behaviour in gas‐liquid airlift reactors.     

Explicit Critical Pressure Ratio for Choked Two‐Phase Homogeneous Nozzle Flow

The throat‐to‐stagnation critical pressure ratio for a frictionless and adiabatic nozzle flow of a homogeneous, nonflashing two‐phase mixture can only be expressed as the numerical solution of a transcendental equation. A simple, physically plausible approximation is herein proposed, which fits well over the whole range of mass flow qualities.     

Generation of Two‐Phase Air‐Water Flow with Fine Microbubbles

Aiming at the development of a low‐cost technology for multipurpose water and surface treatment in the chemical industry and beyond, using microbubbles, a novel scheme of liquid‐gas interaction within a specially designed bubble generator was tested. Its efficiency for the production of microbubbles with a size distribution in the micron range is confirmed. The basic element of the device is a vortex chamber with water supply through tangential ducts, while the gas (air) is introduced in a highly turbulent swirling flow of water in radial direction through the orifice in the gas supply duct, located on the chamber axis. Bubble diameters, bubble velocities in the pipe flow and effect of the output pressure on the bubble size distribution were studied.     

Simulation and Experimental Validation of Two‐Phase Flow in an Aerosol Counter‐Flow Reactor using Computational Fluid Dynamics

A simulation of the hydrodynamic behavior of an aerosol‐counter flow reactor was conducted using an Euler‐Lagrange method. The simulation results were then verified with experiments. The process simulated was a separation process required during the production of biodiesel (fatty acid methyl ester). In this process, the liquid ester/glycerol phases are continuously injected through a hollow cone nozzle with an overpressure of 10⁶ Pa into the reactor, operated at 15000 Pa. The liquid is atomized because of the pressure drop and a liquid particle spray is generated with an inlet velocity of 44.72 m/s. Water vapor of temperature 333 K is injected tangentially through two side, gas inlets with an inlet velocity of 1.2 m/s. Excess methanol is subjected to a mass transfer from the liquid phase into the gas phase, which is withdrawn through the head of the reactor and condensed in an external condenser unit. The stripping of the methanol off the liquid leads to a sharp interface between the glycerol and the ester phase, which can then be easily separated by gravity or pumping. The gas velocity field, pressure field and the liquid particle trajectories were calculated successfully. Simulated dwell time distribution curves were derived and analyzed with the open‐open vessel dispersion model. Experimental dwell time distribution curves were measured, analyzed with the open‐open vessel dispersion model, and compared with the simulated curves. A good consistency between simulated and measured Bodenstein numbers was achieved, but 25 % of the simulated particles exited at the reactor's head, contrary to experimental observations. The difference between simulated and measured dwell times was within one order of magnitude.     

Accuracy of Safety Valve Two‐Phase Mass Flow Capacity Sizing

The method specified in ISO/CD 4126–10 exhibits the highest accuracy compared to other models used in industry and academia. At the same time it allows for an oversizing of the necessary relief area under all test conditions. The average or statistical reproductive accuracy is characterized by an unacceptable logarithmic scatter of about 80 %.     

Phase Inversion and Associated Phenomena in Oil‐Water Vertical Pipeline Flow

Phase inversion and its associated phenomena are experimentally investigated in co‐current upward and downward oil‐water flow in a vertical stainless steel test section (38 mm I.D.). Oil (ρ_o=828 kg/m³, µ_o=5.5 mPa s) and tap water are used as test fluids. Two inversion routes (w/o to o/w and o/w to w/o) are followed in experiments where either the mixture velocity is kept constant and the dispersed phase fraction is increased (type I experiments), or the continuous phase flow rate is kept constant and that of the dispersed phase is increased (type II experiments). By monitoring phase continuity at the pipe centre and at the wall it was found that phase inversion does not happen simultaneously at all locations in the pipe cross‐section. In type I experiments, the velocity ratios (U_o/U_w) where complete inversion appeared acquired the same constant value in both flow directions, although the phase inversion points, based on input phase fractions, were different. In contrast to previous results in horizontal flows, frictional pressure gradient was found to be minimum at the phase inversion point. The interfacial energies of the two dispersions before and after phase inversion, calculated from the measured drop sizes, were found to be different in contrast to the previously suggested criterion of equal energies for the appearance of the phenomenon. In type II experiments the phase inversion point was found to depend on mixture velocity for low and medium velocities but not for high ones. In all cases studied an ambivalent region, commonly reported for inversion in stirred vessels, was not observed.     

Distribution of Two‐Phase Flow across 2.7‐mm Vertically Split U‐Type Junctions

A visual observation of the two‐phase flow across vertically split U‐type junctions and its flow redistribution inside two 2.7‐mm diameter smooth tubes with curvature ratios (2R/D) of 3 and 7, respectively, are reported. The range of mass flux is between 100 and 700 kg/m²s and quality (x) ranges from 0.001 to 0.5. The ratio of liquid distribution between the upper and lower outlet legs is related to the inlet flow pattern, but its influence is reduced at higher mass flux. The difference in liquid flow rates in the lower and upper legs is significantly affected by gravity at a small inlet mass flux, but this difference becomes less profound when the inlet mass flux is increased. The difference between the liquid flux in the upper and lower leg is reduced for the smaller curvature radius due to the reduced effect of gravity between the upper and lower legs. However, there is no consistent trend of gas flow distribution across the U‐type junction as compared to liquid flow distribution. The air mass flux in the upper and lower legs always increase with an increase in both gas quality and the total mass flux.     

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.