Application of PIV to velocity measurements in a liquid film flowing down an inclined cylinder

PIV technique is applied for measurements of instant velocity distributions in a liquid film flowing down an inclined tube in the form of a wavy rivulet. An application of special optical calibration is applied to correct distortion effects caused by the curvature of the interface. A vortex flow of liquid is observed inside a wave hump in the reference system moving with wave phase velocity. Conditionally averaged profiles of longitudinal and transverse components of liquid velocity are obtained for different cross-sections of developed non-linear waves. It is shown that the increase in wave amplitude slightly changes the location of the vortex center. The analysis of modification of vortex motion character due to wavy flow conditions, such as tube inclination angle, film Reynolds number, wave excitation frequency, is fulfilled.     

Numerical analysis on flow dynamics and heat transfer of falling liquid films with interfacial waves

Flow dynamics and heat transfer of falling liquid films with interfacial waves flowing on a vertical plate have been studied with originally proposed numerical simulation method. To discretize basic equations a staggered grid fixed on a physical space is employed. A small amplitude disturbance generated at inflow boundary develops to a solitary wave which consists of a large amplitude roll wave and small amplitude capillary waves. Instantaneous streamwise velocity profiles at the wave crest and trough are very different from a laminar flow. A circulation flow occurs in the roll wave and it affects temperature distributions, especially the strong effect is observed for high Prandtl number liquids. The interfacial wave enhances the heat transfer by two kinds of effects which are a film thinning effect and a convection effect. The dominating effect depends on the Prandtl number. Received on 23 December 1998     

Interfacial instabilities in a stratified flow of two superposed fluids

Here we shall present a linear stability analysis of a laminar, stratified flow of two superposed fluids which are a clear liquid and a suspension of solid particles. The investigation is based upon the assumption that the concentration remains constant within the suspension layer. Even for moderate flow-rates the base-state results for a shear induced resuspension flow justify the latter assumption. The numerical solutions display the existence of two different branches that contribute to convective instability: long and short waves which coexist in a certain range of parameters. Also, a range exists where the flow is absolutely unstable. That means a convectively unstable resuspension flow can be only observed for Reynolds numbers larger than a lower, critical Reynolds number but still smaller than a second critical Reynolds number. For flow rates which give rise to a Reynolds number larger than the second critical Reynolds number, the flow is absolutely unstable. In some cases, however, there exists a third bound beyond that the flow is convectively unstable again. Experiments show the same phenomena: for small flow-rates short waves were usually observed but occasionally also the coexistence of short and long waves. These findings are qualitatively in good agreement with the linear stability analysis. Larger flow-rates in the range of the second critical Reynolds number yield strong interfacial waves with wave breaking and detached particles. In this range, the measured flow-parameters, like the resuspension height and the pressure drop are far beyond the theoretical results. Evidently, a further increase of the Reynolds number indicates the transition to a less wavy interface. Finally, the linear stability analysis also predicts interfacial waves in the case of relatively small suspension heights. These results are in accordance with measurements for ripple-type instabilities as they occur under laminar and viscous conditions for a mono-layer of particles.     

Experimental study on axial development of liquid film in vertical upward annular two-phase flow

Accurate measurements of the interfacial wave structure of upward annular two-phase flow in a vertical pipe were performed using a laser focus displacement meter (LFD). The purpose of this study was to clarify the effectiveness of the LFD for obtaining detailed information on the interfacial displacement of a liquid film in annular two-phase flow and to investigate the effect of axial distance from the air–water inlet on the phenomena. Adiabatic upward annular air–water flow experiments were conducted using a 3 m long, 11 mm ID pipe. Measurements of interfacial waves were conducted at 21 axial locations, spaced 110 mm apart in the pipe. The axial distances from the inlet (z) normalized by the pipe diameter (D) varied over z/D = 50–250. Data were collected for predetermined gas and liquid flow conditions and for Reynolds numbers ranging from Re_G = 31,800 to 98,300 for the gas phase and Re_L = 1050 to 9430 for the liquid phase. Using the LFD, we obtained such local properties as the minimum thickness, maximum thickness, and passing frequency of the waves. The maximum film thickness and passing frequency of disturbance waves decreased gradually, with some oscillations, as flow developed. The flow development, i.e., decreasing film thickness and passing frequency, persisted until the end of the pipe, which means that the flow might never reach the fully developed state. The minimum film thickness decreased with flow development and with increasing gas flow rate. These results are discussed, taking into account the buffer layer calculated from Karman’s three-layer model. A correlation is proposed between the minimum film thickness obtained in relation to the interfacial shear stress and the Reynolds number of the liquid.     

Prediction of the interfacial shear-stress in vertical annular flow

Many improvements of the Wallis correlation for the interfacial friction in annular flow have been proposed in the literature. These improvements give in general a better fit to data, however, their physical basis is not always justified. In this work, we present a physical approach to predict the interfacial shear-stress, based on the theory on roughness in single-phase turbulent pipe flows. Using measured interfacial shear-stress data and measured data on roll waves, which provide most of the contribution to the liquid film roughness, we show that the interfacial shear-stress in vertical annular flow is in very close agreement with the theory. We show that the sand-grain roughness of the liquid film is not equal to four times the mean film thickness, as it is assumed in the Wallis correlation. Instead, the sand-grain roughness is proportional to the wave height, and the proportionality constant can be predicted accurately using the roughness density (or solidity). Furthermore, we show that our annular flow, which is in similar conditions to others in the literature, is fully rough. Hence, the bulk Reynolds number should not appear in the prediction of the interfacial friction coefficient, as is often done in the improvements of the Wallis correlation proposed in the literature.     

Planar laser-induced fluorescence (PLIF) measurements of liquid film thickness in annular flow. Part II: Analysis and comparison to models

Film thickness distributions in upward vertical air–water annular flow have been determined using planar laser-induced fluorescence (PLIF). Film thickness data are frequently used to estimate interfacial shear and pressure loss. This film roughness concept has been used in a number of models for annular flow of varying complexity. The PLIF data are presently applied to the single-zone interfacial shear correlation of Wallis; the more detailed model of Owen and Hewitt; and the two-zone (base film and waves) model of Hurlburt, Fore, and Bauer. For the present data, these models all under-predict the importance of increasing liquid flow on pressure loss and interfacial shear. Since high liquid flow rates in annular flow induce disturbance wave and entrainment activity, further modeling in these areas is advised.     

Nature of the surface heart transfer fluctuation in a hypersonic separated turbulent flow

This paper presents the results of an experimental study of the unsteady nature of a hypersonic separated turbulent flow. The nomimal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5×10⁷/m. The separated flow was generated using finite span forward facing steps. An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make multi-channel measurements of the fluctuating surface heat trtansfer within the separated flow. Conditional sampling analysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness. The compression waves converge into a single leading shock beyond the boundary layer. The shock structure is unsteady and undergoes large-scale motion in the streamwise direction. The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave. There exists a wide band of frequency of oscillations of the shock system. Most of the frequencies are in the range of 1–3 kHz. The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave. This intermittent phenomenon is considered as the consequence of the large-scale shock system oscillations. Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed. The project supported by National Natural Science Foundation of China     

Measurements of velocity,temperature and velocity fluctuation distributions in falling liquid films

The velocity, temperature and velocity fluctuation distributions within falling spindle oil films in an inclined rectangular channel were measured using hot-wire techniques and thin thermocouples. The interfacial shear was caused by cocurrent air flow.The results indicate that the liquid films are as a whole much more laminar-like than turbulent in a range of Reynolds numbers (4γ/μ) up to the experimental limit of 6000. Mixing motion occurs in the vicinity of the interface; however, the flow near the wall surface exhibits no sign of such eddy motions, as predicted by the wall law for single phase turbulent flow. Although velocity fluctuation is observed within films with interfacial shear, mean velocity profiles are approximately the same as those obtained by the laminar film prediction.     

Flow visualization of a recirculating flow by rheoscopic liquid and liquid crystal techniques

Experiments were conducted in a three-dimensional lid-driven cavity flow to study the behavior of longitudinal Taylor-Görtler-like vortices. Flow visualization was accomplished by use of a rheoscopic liquid and of liquid crystals, together with laser-light and white-light sheets, respectively. Photographs of the lighted planes in the flow confirmed the existence of the vortices for a wide range of Reynolds numbers and for stable, neutrally-buoyant and buoyant global flow conditions. As usual the flow visualization revealed flow patterns not deducible by in situ measurements; the liquid crystal photographs give both flow pathlines and temperature distribution on any lighted plane.     

Co-current stratified gas–liquid downflow—Influence of the liquid flow field on interfacial structure

The transition from smooth to wavy stratified flow is studied for various pipe inclination angles with the aim to contribute to the realistic modeling and simulation of long distance two-phase pipe flow. The influence of the liquid flow field on interfacial structure is studied through local axial velocity measurements in the liquid phase in conjunction with other liquid layer characterization experiments. Observations based on fast-video recordings, are used to identify the main patterns of wave evolution. Liquid layer thickness time records are acquired using a parallel wire conductance technique from which mean layer thickness, RMS and power spectra of the fluctuations, as well as wave celerities are calculated. Laser Doppler Anemometry (LDA) is employed to investigate the flow structure in the thin liquid layer both with and without interfacial shear induced by a co-current gas flow. The results reveal the influence of the liquid flow field development on the interfacial structure. In particular, the new data of axial velocity profiles and liquid layer characteristics suggest that the onset of the interfacial waves is strongly affected by the liquid flow structure, possibly by the laminar to turbulent transition within the layer.     

Wavy flow of a liquid film in the presence of a cocurrent turbulent gas flow

Wavy downflow of viscous liquid films in the presence of a cocurrent turbulent gas flow is analyzed theoretically. The parameters of two-dimensional steady-state traveling waves are calculated for wide ranges of liquid Reynolds number and gas flow velocity. The hydrodynamic characteristics of the liquid flow are computed using the full Navier-Stokes equations. The wavy interface is regarded as a small perturbation, and the equations for the gas are linearized in the vicinity of the main turbulent flow. Various optimal film flow regimes are obtained for the calculated nonlinear waves branching from the plane-parallel flow. It is shown that for high velocities of the cocurrent gas flow, the calculated wave characteristics correspond to those of ripple waves observed in experiments.     

Liquid flow patterns about barbotage bubbles

This paper summarizes the results of a flow visualization study on the liquid motion around barbotage bubbles during growth and departure. Flow patterns, as well as for the first time, instantaneous velocities, are reported as a function of time and location about the bubbles. The experiments, employing the hydrogen-bubble technique and high-speed cine photography, were with: water as the liquid, air as the bubbled gas, orifice diameters of 0.116 and 0.252 cm, and different air flow rates; the two limiting cases of constant supply pressure and constant volumetric flow rate were covered. It was found that the liquid around a barbotage bubble assumes two velocity maxima, the first an outward maximum during bubble growth and the second in the opposite direction approximately at the time of bubble departure; further, liquid velocities were found to be higher close to the bubbling site. Certain differences in liquid velocities between the constant pressure and constant flow cases are explained in terms of available theoretical solutions to the bubble growth rate. Qualitative comparisons of the barbotage liquid flow patterns and those recently reported for boiling flow patterns are also presented.     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

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     

Local film thickness and temperature distribution measurement in wavy liquid films with a laser-induced luminescence technique

Heat transfer in falling liquid film systems is enhanced by waviness. Comprehension of the underlying kinetic phenomena requires experimental data of the temperature field with high spatiotemporal resolution. Therefore a non-invasive measuring method based on luminescence indicators is developed. It is used to determine the temperature distribution and the local film thickness simultaneously. Results are presented for the temperature distribution measurement in a laminar-wavy water film with a liquid side Reynolds number of 126 flowing down a heated plane with an inclination angle of 2° at two positions in flow direction. The measured temperature distributions are used to calculate the local heat transfer coefficient for solitary waves at two positions in flow direction.     

Surface instabilities of thin liquid film flow on a rotating disk

 Steady flow of a liquid jet from a nozzle onto the centre of a rotating disk is studied with a streak line method to determine the superficial velocity of the spreading liquid film. Good agreement is found with an asymptotic analysis of the unperturbed flow field. Experimentally, the liquid surface is always perturbed by surface waves which appear as regular spirals, steady in the laboratory system in the low Reynolds number range. It could be shown that wave formation is very sensitive to entrance conditions. Therefore, it is assumed that wave generation is an entrance effect which acts as periodic forcing on the forming liquid film. Wave velocities outside the entrance region are measured and proved to be in good agreement with the prediction of a linear stability theory, as long as the flow rate and entrance perturbations are small. At higher flow rates or stronger disturbances, the radial development of the wave velocities takes on the characteristics predicted by nonlinear stability theories and is in qualitative agreement with experiments performed on an inclined plane. Received: 15 January 1998/Accepted: 8 June 1998     

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.     

Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air

A numerical study of the evaporation in mixed convection of a pure alcohol liquid film: ethanol and methanol was investigated. It is a turbulent liquid film falling on the internal face of a vertical tube. A laminar flow of dry air enters the vertical tube at constant temperature in the downward direction. The wall of the tube is subjected to a constant and uniform heat flux. The model solves the coupled parabolic governing equations in both phases including turbulent liquid film together with the boundary and interfacial conditions. The systems of equations obtained by using an implicit finite difference method are solved by TDMA method. A Van Driest model is adopted to simulate the turbulent liquid film flow. The influence of the inlet liquid flow, Reynolds number in the gas flow and the wall heat flux on the intensity of heat and mass transfers are examined. A comparison between the results obtained for studied alcohols and water in the same conditions is made.     

Experimental study of flow and heat transfer in rib-roughened rectangular channels

Secondary flow patterns, pressure drop and heat transfer in rib-roughened rectangular channels have been investigated. The aspect ratio of the channels is 1–8, and ribs are attached to the wide channel walls in order to set up swirling motions. The geometries tested consist of channels having cross ribs, parallel ribs, cross V-ribs, parallel V-ribs, and multiple V-ribs (Swirl Flow Tube). The flow patterns were investigated using smoke wire visualization and LDV measurements. The smoke wire experiments have been performed at Re=1100 and the LDV measurements at Re=3000 at periodic fully developed conditions. The heat transfer and pressure drop are described by j and f factors for Reynolds numbers from 500 to 15 000. The distributions of axial mean velocity and turbulent fluctuations are strongly influenced by the secondary flows. Large mean velocities and small fluctuations are found in regions where the secondary flow is directed towards a surface, while small mean velocities and large fluctuations are found in regions where the secondary flow is directed away from a surface. The Swirl Flow Tube provides a significant increase in the j factor at Reynolds numbers from 1000 to 2000, but unfortunately also an increase in the f factor. At higher Reynolds numbers, the j and f factors of the Swirl Flow Tube are of the same order of magnitude as for the other rib-roughened channels. It is found that the flow direction in a channel with parallel V-ribs has important influence on the j/f ratio. At Reynolds numbers above 4000, this channel provides the highest j/f ratio if the V-ribs are pointing upstream; while it provides the lowest j/f ratio of all rib configurations, if the V-ribs are pointing downstream.     

Characteristics of an interfacial wave in a horizontal air-water stratified flow

Interfacial wave parameters, in this case the frequency, height, velocity, and slope, were investigated experimentally in a horizontal air-water stratified flow. Experiments were conducted with a parallel wire conductance sensor and PIV visualization in a rectangular channel, of which the width and height are 40 mm and 50 mm, respectively. In the experiments, the flow condition covered the liquid Reynolds number Re_l range of 450 to 3540 and the gas Reynolds number Re_g range of 14,000 to 70,000. The results revealed that the observed wave types according to the flow conditions in the rectangular channel are similar to those in a horizontal pipe. The frequency, height, and slope of the interfacial wave show complicated tendencies according to the combination of Re_g and Re_l, which affects the coalescence and breakup of the wave. Specifically, the wave height and wave slope have opposite tendencies regarding the criterion of Re_g = 34,000. For cases in which Re_g  ≥  34,000, the interfacial drag force significantly affects the height and slope of the disturbance wave. In contrast, for Re_g < 34,000, the growth of the wave has an important effect on the wave parameters. Finally, new empirical correlations for the frequency, height, and slope of the interfacial wave were proposed for application to the development of a droplet entrainment model in a horizontal stratified flow.