Flooding in vertical two-phase flow

The phenomenon of flooding limits the stability of a liquid film falling along the walls of a tube in which a gas is flowing upward. This paper describes a series of experiments to observe and photograph the flooding process. Flooding in a vertical perspex tube was observed by viewing axially along the tube from the top, and flooding on the walls of a stainless steel rod supported centrally within a perspex tube was observed by viewing normally through the tube walls. The resulting observations show that flooding occurs as a result of the formation and motion of a large disturbance wave on the surface of the liquid film. A simple theoretical treatment of flooding is presented, based on the observations of the flooding disturbance. The resulting equations enable the calculation of the gas velocity and the liquid wave height at flooding, and the predictions are in good agreement with experimental observations reported here and elsewhere.     

Behaviour of liquid films and flooding in counter-current two-phase flow. Part 2. Flow in annuli and rod bundles

Experiments are described on the gas velocity at the onset of flooding and the maximum height of the wavy liquid film flowing downwards on a rod surface. On the basis of a simple analysis for a large amplitude wave on the liquid film, a flooding condition relating the maximum wave height to the gas velocity at the onset of flooding is derived. The values predicted by this condition show a good agreement with the measured results.An equivalent diameter of the channel is defined for the flooding velocity. Applying this diameter, the present data for annuli and rod bundles are well correlated by the same empirical equation as that for flow in circular tubes presented previously.     

On the motion of long bubbles in vertical tubes

Potential flow theory has been applied to study the shape and speed of an infinitely long bubble rising through flowing liquid in a vertical tube. In particular, the combined effects of surface tension and externally forced liquid motion are examined. An analytical formula for the bubble rise velocity in stagnant liquid is proposed, and shown to be in good agreement with experimental data for all values of surface tension. Numerical solutions for the bubble velocity in upward flowing liquid are obtained for laminar and turbulent velocity profiles. Approximate expressions for the bubble velocity, where the effects of liquid motion and surface tension are incorporated through the Reynolds and inverse Eötwos, are proposed and compared with experimental data. The predicted changes in bubble shape have, to a large extent, been confirmed through comparisons with photographic evidence for a wide range of parameters.     

Effect of diameter on two-phase pressure drop in narrow tubes

The effect of tube diameter on two-phase frictional pressure drop was investigated in circular tubes with inner diameters of 0.6, 1.2, 1.7, 2.6 and 3.4 mm using air and water. The gas and liquid superficial velocity ranges were 0.01-50 m/s and 0.01-3 m/s, respectively. The gas and liquid flow rates were measured and the two-phase flow pattern images were recorded using high-speed CMOS camera. Unique flow patterns were observed for smaller tube diameters. Pressure drop was measured and compared with various existing models such as homogeneous model and Lockhart-Martinelli model. It appears that the dominant effect of surface tension shrinking the flow stratification in the annular regime is important. It was found that existing models are inadequate in predicting the pressure drop for all the flow regimes visualized. Based on the analysis of present experimental frictional pressure drop data a correlation is proposed for predicting Chisholm parameter “C” in slug annular flow pattern. For all other flow regimes Chisholm’s original correlation appears to be adequate except the bubbly flow regime where homogeneous model works well. The modification results in overall mean deviation of pressure drop within 25% for all tube diameters considered. This approach of flow regime based modification of liquid gas interaction parameter appears to be the key to pressure drop prediction in narrow tubes.     

Experimental study on hydrodynamic characteristics of upward gas–liquid slug flow

To clarify the impacts of the hydrodynamic boundary layer and the diffusion boundary layer in the near wall zone on gas–liquid two-phase flow induced corrosion in pipelines, the hydrodynamic characteristics of fully developed gas–liquid slug flow in an upward tube are investigated with limiting diffusion current probes, conductivity probes and digital high-speed video system. The Taylor bubble and the falling liquid film characteristics are studied, the effects of various factors are examined, and the experimental results are compared with the data and models available in literature. The length of Taylor bubble, the local void fraction of the slug unit and the liquid slug, the shear stress and mass transfer coefficient in the near wall zone, are all increased with the increase of superficial gas velocity and decreased with the increase of superficial liquid velocity, whereas the length of liquid slug and the liquid slug frequency are changed contrarily. The alternate wall shear stress due to upward gas–liquid slug flow is considered to be one of the major causes for the corrosion production film fatigue cracking. A normalized formula for mass transfer coefficient is obtained based on the experimental data.     

Drop size measurements for annular two-phase flow in a 20 mm diameter vertical tube

As part of a study on the effect of tube diameter on the mean drop size and liquid film flow rate in annular two-phase flow, data was obtained for the vertical upflow of an air-water system in a 20 mm internal diameter tube, held at a pressure of 1.5 bar and ambient temperature. This complements data taken in earlier experiments on 10 and 32 mm tubes. Increases in the superficial gas velocity caused reductions in the mean drop size whilst increasing the liquid mass flux in all but the lowest gas velocity case, caused the drop size to rise. Comparisons were made between the current drop size data and that from a 10 mm and 32 mm internal diameter tube, for similar conditions of temperature and pressure. The current drop size measurements, which fall between those from earlier work, confirm the dependence of drop size on tube diameter. The performance of several drop size correlations have been tested. Because the correlations do not account for the influence of tube diameter, they fail to predict the drop size data accurately. The influence of gas and liquid flow rate on the measured film flow rate show trends similar to those seen in data from the 10 mm and 32 mm diameter tubes. Models, to calculate the entrained liquid mass flux were tested; good predictions were given.     

Nonlinearwaves in a liquid film-gas flow system

The instability and regular nonlinear waves in the film of a heavy viscous liquid flowing along the wall of a round tube and interacting with a gas flow are investigated. The solutions for the wave film flows are numerically obtained in the regimes from free flow-down in a counter-current gas stream to cocurrent upward flow of the film and the gas at fairly large gas velocities. Continuous transition from the counter-current to the cocurrent flow via the state with a maximum amplitude of nonlinear waves and zero values of the liquid flow rate and the phase velocity is investigated. The Kapitsa-Shkadov method is used to reduce a boundary value problem to a system of evolutionary equations for the local values of the layer thickness and the liquid flow rate.     

A film model for the prediction of flooding and flow reversal for gas-liquid flow in vertical tubes

A model is presented which demonstrates that the process of flooding and flow reversal can be explained on the basis of a film mechanism. The model predicts well the gas flow rate at which flooding and flow reversal begins and ends for a given liquid flow rate and the presence of a hysteresis loop between flooding and flow reversal. The predictions of the theory are in satisfactory agreement with experimental flooding data.     

Interpretation of atomization rates of the liquid film in gas-liquid annular flow

For vertical gas-liquid annular flow the fraction of the liquid in the gas is controlled by the rate of atomization of the liquid film flowing along the wall and the rate of deposition of droplets entrained in the gas. Measurements of the rate of atomization are interpreted by a Kelvin-Helmholtz mechanism. Small wavelets on the liquid film are visualized to be entrained when wave-induced variations in the gas pressure cannot be counterbalanced by surface tension effects.     

Numerical study of bubble rise and interaction in a viscous liquid

In this article, the flow instabilities during the rise of a single bubble in a narrow vertical tube are studied using a transient two-dimensional/axisymmetric model. To predict the shape of the bubble deformation, the Navier-Stokes equations in addition to an advection equation for liquid volume fraction are solved. A modified volume-of-fluid technique based on Youngs' algorithm is used to track the bubble deformation. To validate the model, the results of simulations for terminal rise velocity and bubble shape are compared with those of the experiments. The effect of different parameters such as initial bubble radius, channel height, liquid viscosity and surface tension on the shape and rise velocity of the bubble is investigated.     

Flow pattern and flow characteristics for counter-current two-phase flow in a vertical round tube with wire-coil inserts

Flow pattern, void fraction and slug rise velocity on counter-current two-phase flow in a vertical round tube with wire-coil inserts are experimentally studied. Flow pattern and slug rise velocity are measured visually with a video camera. The void fraction is measured by the quick-closing valve method. Four kinds of coils with different coil pitches and coil diameters are used as inserts. The presence of wire-coil inserts induces disturbance into gas and liquid flows so that the shape and motion of gas slug or bubbles in a wire-coil inserted tube are quite different from those observed in a smooth tube without insert. The bubbly flow occurs in the low gas superficial velocity region in the wire-coil inserted tube, while the slug or churn/annular flow only appears in the smooth tube without insert over the all test range. The measured slug rise velocity in the wire-coil inserted tube is higher than that in the smooth tube. With modified mean flow velocity calculated with core area, the slug rise velocity in wire-coil tube inserted is in good agreement with Nicklin's correlation. The void fraction in a wire-coil inserted tube is lower than that in a smooth tube in the range of high gas superficial velocities. By introducing a simple assumption on considering the effective flowing area, the measured void fractions in a wire-coil inserted tube are in relatively good agreement with the predicted result based on the drift flux model proposed by others with the correlation for slug rise velocity given by others when the coil pitch is dense.     

热毛细对流速度场测试研究

为了探索微重力状态下表面张力驱动流速度场的测试技术，本文在地面进行了模拟实验，对位于热壁下流体中气泡周围的表面张力驱动流进行了研究，并用ＰＩＶ技术测量了流场的速度分布．     

A novel numerical scheme for the investigation of surface tension effects on growth and collapse stages of cavitation bubbles

In the present study the effects of surface tension on the growth and collapse stages of cavitation bubbles are studied individually for both spherical and nonspherical bubbles. The Gilmore equation is used to simulate the spherical bubble dynamics by considering mass diffusion and heat transfer. For the collapse stage near a rigid boundary, the Navier–Stokes and energy equations are used to simulate the flow domain, and the VOF method is adopted to track the interface between the gas and the liquid phases. Simulations are divided into two cases. In the first case, the collapse stage alone is considered in both spherical and nonspherical situations with different conditions of bubble radius and surface tension. According to the results, surface tension has no significant effects on the flow pattern and collapse rate. In the second case, both the growth and collapse stages of bubbles with different initial radii and surface tensions are considered. In this case surface tension affects the growth stage considerably and, as a result, the jet velocity and collapse time decrease with increasing surface tension coefficient. This effect is more significant for bubbles with smaller radii.     

Numerical study on the effect of initial flow velocity on liquid film thickness of accelerated slug flow in a micro tube

Numerical simulation of air–water slug flows accelerated from steady states with different initial velocities in a micro tube is conducted. It is shown that the liquid film formed between the gas bubble and the wall in an accelerated flow is significantly thinner than that in a steady flow at the same instantaneous capillary number. Specifically, the liquid film thickness is kept almost unchanged just after the onset of acceleration, and then gradually increases and eventually converges to that of an accelerated flow from zero initial velocity. Due to the flow acceleration, the Stokes layer is generated from the wall, and the instant velocity profile can be given by superposition of the Stokes layer and the initial parabolic velocity profile of a steady flow. It is found that the velocity profile inside a liquid slug away from the bubble can be well predicted by the analytical solution of a single-phase flow with acceleration. The change of the velocity profile in an accelerated flow changes the balance between the inertia, surface tension and viscous forces around the meniscus region, and thus the resultant liquid film thickness. By introducing the displacement thickness, the existing correlation for liquid film thickness in a steady flow (Han and Shikazono, 2009) is extended so that it can be applied to a flow with acceleration from an arbitrary initial velocity. It is demonstrated that the proposed correlation can predict liquid film thickness at Re < 4600 within the range of ±10% accuracy.     

Direct numerical simulation of a liquid sheet in a compressible gas stream in axisymmetric and planar configurations

A thin liquid sheet present in the shear layer of a compressible gas jet is investigated using an Eulerian approach with mixed-fluid treatment for the governing equations describing the gas–liquid two-phase flow system, where the gas is treated as fully compressible and the liquid as incompressible. The effects of different topological configurations, surface tension, gas pressure and liquid sheet thickness on the flow development of the gas–liquid two-phase flow system have been examined by direct solution of the compressible Navier–Stokes equations using highly accurate numerical schemes. The interface dynamics are captured using volume of fluid and continuum surface force models. The simulations show that the dispersion of the liquid sheet is dominated by vortical structures formed at the jet shear layer due to the Kelvin–Helmholtz instability. The axisymmetric case is less vortical than its planar counterpart that exhibits formation of larger vortical structures and larger liquid dispersion. It has been identified that the vorticity development and the liquid dispersion in a planar configuration are increased at the absence of surface tension, which when present, tends to oppose the development of the Kelvin–Helmholtz instability. An opposite trend was observed for an axisymmetric configuration where surface tension tends to promote the development of vorticity. An increase in vorticity development and liquid dispersion was observed for increased liquid sheet thickness, while a decreasing trend was observed for higher gas pressure. Therefore surface tension, liquid sheet thickness and gas pressure factors all affect the flow vorticity which consequently affects the dispersion of the liquid.      

Effects of a countercurrent gas flow on falling-film mode transitions between horizontal tubes

A liquid film falling between horizontal tubes is known to take the form of droplets, jets or sheets, depending on the liquid flow rate; the form of the flow is the so-called “falling-film mode”. Although previously neglected in studies of mode transition, a countercurrent gas flow often exists in falling-film heat exchangers, and its effect on the liquid flow might be important: it could impact the flow regime, lead to local “dryout,” and decrease the heat transfer rate. Experiments are conducted to explore the effects of a countercurrent gas flow and liquid feeding length on falling-film mode transitions for a liquid flowing over horizontal tubes. The effects on mode transition are shown to depend on fluid properties and are explained in terms of unsteadiness and film thickness. In general, transition hysteresis is reduced with an increasing gas velocity. A correlation is developed to predict the countercurrent gas flow effects on falling-film mode transitions. The liquid feeding length can affect mode transitions in quiescent surroundings and when a countercurrent gas flow imposed.     

Bubble generation and transport in a microfluidic device with high aspect ratio

Bubble generation and transport in a micro-device composed of a micro-T-junction and a following serpentine micro-channel was experimentally investigated. It has a rectangular cross-sectional with an aspect ratio of 7.425. Air and water were used as gas and liquid, respectively. Mixtures of water–glycerol and water–Tween-20 were used to study the effects of liquid viscosity and surface tension. Compared with previous T-junction bubble generation, the liquid and gas inlets orientation was switched in this work. The continuous flow was driven from the perpendicular channel and the dispersion flow was from the main channel. It shows that the break-up process has three periodic steps under certain operating conditions. The dimensionless bubble length L/w in the micro-channel with high aspect ratio is much larger than that in square microchannels. A correlation is proposed to correlate L/w with liquid flow rate J_L, gas flow rate J_G, and liquid viscosity μ_L. Surface tension σ can change the bubble shape but almost does not affect the bubble length in this fast break-up process. Additionally, a long bubble may be broken up at the corners at the same time because the locations of gas and liquid are exchanged relative to the concave and convex portions of an elbow after a turn which may result in the change of fluid velocities and gas–liquid pressure drop.     

内分液结构调控竖直管间歇流型建模研究

内分液流型调控管依靠微尺度网孔阻气通液的毛细力学特性,调控气液两相间歇流型以实现传热强化.基于Lockhart-Martinelli 分相模型以及Zuber-Findlay 漂移流动模型,建立描述内分液竖直管内流体动力特性的一维数学模型. 采用模型求解实验工况,计算结果与实验结果误差均在20% 以内. 计算发现,液速对流动现象起决定作用,而气速影响通过丝网的渗透程度. 在定性分析基础上,采用三角立方插值与最小二乘B 样条拟合获得了流动特性与气速、液速的定量函数关系. 据此得出结论,当Re_l < 693 7 时,一定出现第1 类工况;当Re_l > 693 7,且Re_g < 67 时,可能会出现第2 类工况,此时较低的气速会促进第2 类工况的出现. 根据建立的模型与拟合关系式可实现内分液调控管的优化设计.      

Propagation of pressure waves in two-phase flow

We deal with a pressure wave of finite amplitude propagating in a gas and liquid medium or in the fluid in an elastic tube. We study the effects of pipe elasticity on the propagation velocity of the pressure wave. Pressure waves of finite amplitude progressing in the two-phase flow are treated considering the void fraction change due to pressure rise. The propagation velocity of the two-phase shock wave is also investigated, and the behavior of the reflection of the pressure wave at the rigid wall is analyzed and compared to that in a pure gas or liquid. The results are compared to experimental data of a pressure wave propagating in the two-phase flow in a vertical shock tube.     

Flooding in two-phase counter-current flows: Numerical investigation of the gas–liquid wavy interface using the Navier–Stokes equations

This paper is devoted to a theoretical analysis of counter-current gas–liquid wavy film flow between vertical plates. We consider two-dimensional nonlinear waves on the interface over a wide variation of parameters. We use the Navier–Stokes equations in their full statement to describe the liquid phase hydrodynamics. For the gas phase equations, we use the Benjamin-Miles approach where the liquid phase is a small disturbance for the turbulent gas flow. We find a region of the superficial velocity where we have two solutions at one set of the problem parameters and where the flooding takes place. We calculate the flooding dependences on the gas/liquid physical properties, on the liquid Reynolds number and on the distance between the plates. These computations allow us to present the correlation for the onset of flooding that based on the fundamental equations and principles.