Weakly nonlinear stability analysis of a falling film with countercurrent gas flow

Weakly nonlinear stability analysis of a falling film with countercurrent gas–liquid flow has been investigated. A normal mode approach and the method of multiple scales are employed to carry out the linear and nonlinear stability solutions for the film flow system. The results show that both supercritical stability and subcritical instability are possible for a film flow system when the gas flows in the countercurrent direction. The stability characteristics of the film flow system are strongly influenced by the effects of interfacial shear stress when the gas flows in the countercurrent direction. The effect of countercurrent gas flow in a falling film is to stabilize the film flow system.     

Macro-to-microchannel transition in two-phase flow: Part 1 - Two-phase flow patterns and film thickness measurements

The classification of macroscale, mesoscale and microscale channels with respect to two-phase processes is still an open question. The main objective of this study focuses on investigating the macro-to-microscale transition during flow boiling in small scale channels of three different sizes with three different refrigerants over a range of saturation conditions to investigate the effects of channel confinement on two-phase flow patterns and liquid film stratification in a single circular horizontal channel (Part 2 covers the flow boiling heat transfer and critical heat flux). This paper presents the experimental two-phase flow pattern transition data together with a top/bottom liquid film thickness comparison for refrigerants R134a, R236fa and R245fa during flow boiling in small channels of 1.03, 2.20 and 3.04 mm diameter. Based on this work, an improved flow pattern map has been proposed by determining the flow patterns transitions existing under different conditions including the transition to macroscale slug/plug flow at a confinement number of Co ≈ 0.3-0.4. From the top/bottom liquid film thickness comparison results, it was observed that the gravity forces are fully suppressed and overcome by the surface tension and shear forces when the confinement number approaches 1, Co ≈ 1. Thus, as a new approximate rule, the lower threshold of macroscale flow is Co = 0.3-0.4 while the upper threshold of symmetric microscale flow is Co ≈ 1 with a transition (or mesoscale) region in-between.     

Measurement of the liquid film thickness in micro tube slug flow

Slug flow is one of the representative flow regimes of two-phase flow in micro tubes. It is well known that the thin liquid film formed between the tube wall and the vapor bubble plays an important role in micro tube heat transfer. In the present study, experiments are carried out to clarify the effects of parameters that affect the formation of the thin liquid film in micro tube two-phase flow. Laser focus displacement meter is used to measure the thickness of the thin liquid film. Air, ethanol, water and FC-40 are used as working fluids. Circular tubes with five different diameters, D = 0.3, 0.5, 0.7, 1.0 and 1.3 mm, are used. It is confirmed that the liquid film thickness is determined only by capillary number and the effect of inertia force is negligible at small capillary numbers. However, the effect of inertia force cannot be neglected as capillary number increases. At relatively high capillary numbers, liquid film thickness takes a minimum value against Reynolds number. The effects of bubble length, liquid slug length and gravity on the liquid film thickness are also investigated. Experimental correlation for the initial liquid film thickness based on capillary number, Reynolds number and Weber number is proposed.     

Algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow, focusing in particular on momentum and heat transfer (so-called ‘convective boiling’) through the annular liquid film. In contrast with single-phase wall-bounded flow theory, shear-driven annular liquid films are assumed here to behave as fluid-bounded flows, mostly interacting with the shearing gas-entrained droplets core flow. Besides providing velocity and temperature profiles through the liquid film, the turbulence model proposed here predicts key parameters such as the average liquid film thickness, the void fraction and the convective boiling heat transfer coefficient with accuracies comparable or better than those of leading design correlations. This turbulence model is part of a unified annular flow modeling suite that includes methods to predict the entrained liquid fraction and the axial frictional pressure gradient. The underlying heat transfer database covers nine fluids (water, two hydrocarbons and six refrigerants) for vertical and horizontal tubes of 1.03-14.4 mm diameter and pressures of 0.1-7.2 MPa. Importantly, this study shows that there appears to be no macro-to-microscale transition when it comes to annular flow. Simply better physical modeling is required to span this range.     

Behaviour of liquid films and flooding in counter-current two-phase flow—Part 1. Flow in circular tubes

Experimental results are presented on the flooding gas velocity in tubes over a wide range of parameters—tube diameter, tube length, liquid flow rate, liquid viscosity and surface tension. The flooding phenomenon is caused by interaction between the waves on the liquid film and the upward gas stream. By measuring variation of the maximum height of the wavy liquid films with an increase of the gas flow rate, the complicated effects of tube length and surface tension on flooding are revealed. The data of the flooding velocity are empirically correlated in termes of nondimensional groups for each tube length.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

Liquid film circumferential asymmetry prediction in horizontal annular two-phase flow

This study considers the prediction of the degree of asymmetry in the circumferential distribution of the liquid film in the tube cross section of horizontal annular gas–liquid two-phase flow, endemic of the lower region of this flow regime near the stratified-wavy flow transition boundary. Focusing on disturbance waves as the predominant mechanism for transporting the liquid in the annular film from the bottom to the top of the tube to counterbalance the draining effect of gravity, a new prediction method for the degree of asymmetry in the annular liquid film is proposed that outperforms existing correlations. Flow pattern maps for horizontal gas–liquid two-phase flow of frequent use in the design of evaporators and condensers can thus be explicitly updated to account for both symmetric and asymmetric annular flows. The underlying experimental database contains 184 measured liquid film circumferential profiles, corresponding to 1276 local liquid film thickness measurements collected from 15 different literature studies for tube diameters from 8.15 mm to 95.3 mm.     

Critical friction factor modeling of horizontal annular base film thickness

Measurements of liquid base film thickness distribution have been obtained for 206 horizontal annular two-phase (air–water) flow conditions in 8.8 mm, 15.1 mm, and 26.3 mm ID tubes. It is found that the trends in base film thickness measurement do not match trends in the literature for average film thickness, which considers waves and base film together. An iterative critical friction factor model is used to model circumferentially-averaged base film thickness; an explicit, empirical correlation is also provided. Asymmetry is well-correlated by a modified Froude number based on the correlated base film thickness and the gas mass flux. The iterative model is also extended to estimate the critical film flow rate.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

Effect of density ratio on critical heat flux in closed end vertical tubes

An experimental study was conducted to measure the critical axial heat flux in countercurrent two phase flow of liquid and its vapour in a closed-end vertical tube. The experimental results for four different fluids; carbon-tetrachloride, normal hexane, ethyl alcohol and water, were reduced to give a correlation for evaluation of the flooding critical heat flux. Results of other investigators on vertical heated tubes and vertical thermosyphons were also reduced and compared with the present experimental results. The effect of the density ratio of liquid to its vapour on the critical heat flux was shown. The length to diameter ratio of the test section was shown to have an influence on the flooding critical heat flux and was included in the correlation obtained.     

Primary instabilities of liquid film flow sheared by turbulent gas stream

Evolution of excited waves on a viscous liquid film has been investigated experimentally for the annular gas–liquid flow in a vertical tube. For the first time the dispersion relations are obtained experimentally for linear waves on liquid film surface in the presence of turbulent gas flow. Both cocurrent and countercurrent flow regimes are investigated. As an example of comparison with theory, the experimental data are compared to the results of calculations based on the Benjamin quasi-laminar model for turbulent gas flow. The calculation results are found to be in good agreement with experiments for moderate values of film Reynolds number.     

Wave structure in the radial film flow with a circular hydraulic jump

 A circular hydraulic jump is commonly seen when a circular liquid jet impinges on a horizontal plate. Measurements of the film thickness, jump radius and the wave structure for various jet Reynolds numbers are reported. Film thickness measurements are made using an electrical contact method for regions both upstream and downstream of the jump over circular plates without a barrier at the edge. The jump radius and the separation bubble length are measured for various flow rates, plate edge conditions, and radii. Flow visualization using high-speed photography is used to study wave structure and transition. Waves on the jet amplify in the film region upstream of the jump. At high flow rates, the waves amplify enough to cause three-dimensional breakdown and what seems like transition to turbulence. This surface wave induced transition is different from the traditional route and can be exploited to enhance heat and mass transfer rates. Received: 25 April 2000/Accepted: 1 June 2001     

Two-phase flow patterns in short horizontal rectangular microchannels

The two-phase flow in a short horizontal channel of a rectangular cross-section with the height of 100–500 µm and width of 9–40 mm was studied experimentally. The use of the Schliren and fluorescent methods made it possible to reveal the flow of liquid in the channel and to determine its characteristics quantitatively. The features of the churn, jet and drop flow patterns were studied in details. Two particular regimes that can be distinguished represent formation of immobile drops on the channel walls because of the liquid film or liquid bridges breakage and appearance of mobile drops due to the two-phase flow instabilities. It is found out that formation of various two-phase flow patterns and transitions between them are determined by instabilities of the liquid–gas flow in the side parts of a channel. Frontal instability has been observed during the liquid–gas interaction in the region of liquid output from the nozzle. It is shown that a change in the height and width of the horizontal channels has a substantial effect on the boundaries between the flow regimes. One of the results is that the region of the churn regime increases significantly with decreasing thickness of the channel.     

Wave frequency of falling liquid films and the effect on reflux condensation in vertical tubes

Visual investigations of effective wave frequency, structure and formation of isothermal and fully developed falling liquid films inside vertical tubes are presented without and with countercurrent flow of vapour. These are supplemented by novel transmissivity measurements leading to considerable improvement of the established wave shape and liquid film structure classifications. In addition, reflux condensation heat transfer data evaluated for the limiting case of zero shear stress are represented in terms of the wave factor. These are correlated with Reynolds and Kapitza numbers and they are interrelated with the observed wave characteristics.     

Pressure drop prediction in annular two-phase flow in macroscale tubes and channels

A new prediction method for the frictional pressure drop in annular two-phase flow is presented. This new prediction method focuses on the aerodynamic interaction between the liquid film and the gas core in annular flows, and explicitly takes into account the asymmetric liquid film distribution in the tube cross section induced by the action of gravity in horizontal tubes operated at low mass fluxes. The underlying experimental database contains 6291 data points from the literature with 13 fluid combinations (both single-component saturated fluids such as water, carbon dioxide and refrigerants R12, R22, R134a, R245fa, R410a, R1234ze, and two-component fluids such as water-argon, water-nitrogen, alcohol-argon, water plus alcohol-argon and water-air), vertical and horizontal tubes and annuli with diameters from 3 mm to 25 mm, and both adiabatic and evaporating flow conditions. The new prediction method is very simple to implement and use, is physically based and outperforms existing pressure drop correlations (mean absolute error of 12.9%, and 7 points out of 10 captured to within ±15%).     

Effects of fluid properties and pipe diameter on two-phase flow patterns in horizontal lines

Extensive new data have been obtained on the transitions between two-phase flow patterns during co-current gas liquid flow in horizontal lines. Fluid properties were varied in a systematic manner to determine the effects of liquid viscosity, liquid density, interfacial tension and gas density. Line sizes varied from 1.2 to 5 cm for most of the tests. Visual observations were supplemented by an analysis of pressure drop fluctuations and hence the present data are believed to be less subjective than most past observations.

The transition data from the present tests, as well as available literature data, were compared to the most frequently used transition line correlations. In almost all cases serious deficiencies were observed. Revised dimensionless correlations which fit present data, and those previously available, are presented.     

Wave behavior in horizontal annular air–water flow

In this study, non-intrusive pressure drop, liquid base film thickness distribution, and wave behavior measurements have been obtained for 206 horizontal annular two-phase (air–water) flow conditions in 8.8, 15.1, and 26.3 mm ID tubes. Reliable wave velocity measurements are available for 185 of these flow conditions, while 131 flow conditions allow for reliable wave frequency measurements. The wave velocity is found to be predicted to within 9% by gas friction velocity and 6% by an optimized correlation of similar structure. Wave frequency can also be predicted with a simple correlation to within 5% for the 8.8 and 15.1 mm tubes, but a separate relation is required to explain 26.3 mm frequency data. The differences in wave behavior between the annular and wavy-annular/wavy regimes are also discussed.     

Evaporation of single component viscous liquids in a metal falling film evaporator

Heat transfer for single component falling film evaporation has been investigated in a stainless steel single tube falling film evaporator. The tube had a heated length of 2,500 mm. Propylene glycol and cyclohexanol have been used as evaporating media. Liquid film running down the tube, is formed on the inner side of the tube. For the distribution of liquid two different devices were examined. Process equipment was operated in pump-around mode with the distillate being condensed and recycled. Results show that none of the available correlations for heat transfer in falling film evaporation is able to describe the results qualitatively as well as quantitatively. Using different film distribution devices, a significant influence of the Reynolds number for the transition from laminar to turbulent flow is seen. However, differences between experimental results and correlations in literature cannot be explained only by usage of different film distributions, in particular when the correlation is based on measurements with a different tube length. A model approach is presented for cyclohexanol as evaporating medium with a flat weir as film distributor.     

Forced convective flow and heat transfer of upward cocurrent air–water slug flow in vertical plain and swirl tubes

This experimental study comparatively examined the two-phase flow structures, pressured drops and heat transfer performances for the cocurrent air–water slug flows in the vertical tubes with and without the spiky twisted tape insert. The two-phase flow structures in the plain and swirl tubes were imaged using the computerized high frame-rate videography with the Taylor bubble velocity measured. Superficial liquid Reynolds number (Re_L) and air-to-water mass flow ratio (AW), which were respectively in the ranges of 4000–10000 and 0.003–0.02 were selected as the controlling parameters to specify the flow condition and derive the heat transfer correlations. Tube-wise averaged void fraction and Taylor bubble velocity were well correlated by the modified drift flux models for both plain and swirl tubes at the slug flow condition. A set of selected data obtained from the plain and swirl tubes was comparatively examined to highlight the impacts of the spiky twisted tape on the air–water interfacial structure and the pressure drop and heat transfer performances. Empirical heat transfer correlations that permitted the evaluation of individual and interdependent Re_L and AW impacts on heat transfer in the developed flow regions of the plain and swirl tubes at the slug flow condition were derived.     

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work, co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 mm in diameter. Each tube includes two 10 m long pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly flow, stratified flow, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high-speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood–Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/non-Newtonian liquid.