Effect of gravity on axial development of bubbly flow at low liquid Reynolds number

In view of the great importance of two geometrical parameters such as void fraction and interfacial area concentration to the accurate two-phase flow analysis at microgravity conditions, axial developments of flow parameters such as void fraction, interfacial area concentration, bubble Sauter mean diameter, and bubble number density were measured in bubbly flow at microgravity and low liquid Reynolds number conditions where the gravity effect on the flow parameters were pronounced. A total of seven data sets were acquired in the flow range of the void fraction from 1.01% to 3.36% and the liquid Reynolds number from 1,400 to 4,750. The measurements were also performed in the similar flow range at normal gravity conditions. The transport mechanisms of the flow parameters are discussed in detail based on the data measured at normal and microgravity conditions, and the drift-flux model developed at microgravity conditions are compared with the measured data.An erratum to this article can be found at     

An experimental study of two-phase flow in simulated reduced-gravity condition: dispersed droplet to slug flow transition and slug flow

The results from an experimental study of reduced-gravity two-phase flows are reported in this paper. The experiments were conducted in simulated reduced-gravity conditions in a ground-based test facility with a circular test section of 25 mm inner diameter. The flow conditions for which data were acquired lie in the dispersed droplet to slug flow transition and slug flow regime. Local data were acquired for 17 different flow conditions at three axial locations. The acquired data complement and extend those discussed in an earlier paper by the authors (Vasavada et al. in, Exp Fluids 43: 53–75, 2007). The radial profiles and axial changes in the local data are analyzed and discussed in this paper. The area-averaged data, in conjunction with the local data, are discussed to highlight important interaction mechanisms occurring between fluid particles, i.e., drops. The data clearly show the effect of progressive coalescence leading to formation of slug drops. Furthermore, the shape of slug drops in reduced-gravity conditions was observed to be different from that in normal-gravity case. The analyses presented here show the presence of drop coalescence mechanisms that lead to the formation of slug drops and transition from dispersed droplet flow to the slug flow regime. The most likely causes of the coalescence mechanism are random collision of drops driven by turbulence eddies in the continuous phase and wake entrainment of smaller drops that follow preceding larger drops in the wake region. Data from flow conditions in which the breakup mechanism due to impact of turbulent eddies on drops illustrate the disintegration mechanism.     

Bubbly-to-cap bubbly flow transition in a long-26 m vertical large diameter pipe at low liquid flow rate

The concurrent upward two-phase flow of air and water in a long vertical large diameter pipe with an inner diameter (D) of 200 mm and a height (z) of 26 m (z/D = 130) was investigated experimentally at low superficial liquid velocities from 0.05009 to 0.3121 m/s and the superficial gas velocities from 0.01779 to 0.5069 m/s. The resultant void fractions range from 0.03579 to 0.4059. According to the observations using a high speed video camera, the flow regimes of bubbly, developing cap bubbly and fully-developed cap bubbly flows prevailed in the flows. The developing cap bubbly flow appeared as a flow regime transition from bubbly to fully-developed cap bubble flow in the vertical large diameter pipe. The developing cap bubbly flow changes gradually and lasts for a long time period and a wide axial region in the flow direction, in contrast to a sudden transition from bubbly to slug flows in a small diameter pipe. The analysis in this study showed that the flow regime transition depends not only on the void fraction but also on the axial distance in the flow and the pipe diameter. The axial flow development brings about the transition to happen in a lower void fraction flow and the increase of pipe diameter causes the transition to happen in a higher void fraction flow. The measured void fraction showed an N-shaped axial changing manner that the void fraction increases monotonously with axial position in the bubbly flow, decreases non-monotonously with axial position in the developing cap bubbly flow, and increases monotonously again with axial position in the fully-developed cap bubbly flow. The temporary void fraction decrease phenomenon in the transition region from bubbly to cap bubbly flow can be attributed to the formation of medium to large cap bubbles and their gradual growth into the maximum size of cap bubble and/or cluster of large cap bubbles in the developing cap bubbly flow. In order to predict the N-shaped axial void fraction changing behaviors in the flow regime transition from bubbly to cap bubbly flow, the existing 12 drift flux correlation sets for large diameter pipes are reviewed and their predictabilities are studied against the present experimental data. Although some drift flux correlation sets, such as those of Clark and Flemmer (1986) and Hibiki and Ishii (2003), can predict the present experimental data with reasonable average relative deviations, no drift flux correlation set for distribution parameter and drift velocity can give a reliable prediction for the observed N-shaped axial void fraction changing behaviors in the region from bubbly to cap bubbly flow in a vertical large diameter pipe.     

Characteristics of developing vertical bubbly flow under normal and microgravity conditions

The axial development of the void fraction profile, interfacial area concentration and Sauter mean bubble diameter of adiabatic nitrogen-water bubbly flows in a 9 mm-diameter pipe were measured using stereo image processing under normal and microgravity conditions. The flow measurements were performed at four axial locations (axial distance from the inlet, z normalized by the pipe diameter, D, z/D = 5, 20, 40 and 60) and with various flows: superficial gas velocity of 0.00840-0.0298 m/s, and superficial liquid velocity of 0.138-0.914 m/s. The effect of gravity on radial distribution of bubbles and the axial development of two-phase flow parameters is discussed in detail based on the obtained database and visual observation. Following Serizawa-Kataoka’s phase distribution pattern criteria under normal gravity conditions, the phase distribution pattern map was developed. Similar to normal gravity two-phase flows, wall, core and intermediate void peak patterns are observed under microgravity conditions but a transition void distribution pattern is not observed in the current experimental conditions. The data obtained in the current experiment are expected to contribute to the benchmarking of CFD simulation of phase distribution pattern and interfacial area concentration in forced convective pipe flow under microgravity conditions.     

Interfacial structure of air-water two-phase flow in a relatively large pipe

Miniaturized four-sensor conductivity probes are used to study flow structure development in air-water bubbly flow, cap-bubbly flow, and transition to slug flow. The measurements are performed at three different elevations in a vertical round pipe with an inner diameter of 101.6 mm. The time-averaged local void fraction, interfacial velocity, and bubble number frequency are measured by the conductivity probes. Also, the interfacial area concentration and averaged bubble Sauter mean diameter are obtained. A detailed representation of the flow structure is revealed by investigating the acquired data. Furthermore, comparisons of the data at three elevations demonstrate the development of the interfacial structure along the flow direction due to bubble interactions and hydrodynamic effects. This may provide the community with a better knowledge about two-phase flow in a relatively large pipe. In addition, these data can also serve as an experimental database for investigation of the interfacial area transport in large-pipe two-phase flow. Published online: 19 November 2002 This work was performed under the auspices of the U.S. Nuclear Regulatory Commission through the Institute of Thermal-hydraulics.     

Damping of tubes due to internal two-phase flow

Two-phase internal flow is present in many piping system components. Although two-phase damping is known to be a significant constituent of the total damping, the energy dissipation mechanisms that govern two-phase damping are not well understood. In this paper, damping of three different clamped–clamped tubes subjected to two-phase air–water internal axial flow is investigated. Experimental data are reported, showing a strong dependence of two-phase damping on void fraction, flow velocity and flow regime. Data-points plotted on two-phase flow pattern maps indicate that damping is greater in a bubbly flow regime. The two-phase damping ratio reaches a maximum value at the highest void fraction before the transition to a churn flow regime. An analytical model that relates the two-phase damping ratio to the interface surface area is proposed. The model is based on rigid spherical bubbles in cubic elementary flow volumes. The analytical results are well correlated with the experiments.     

Visualization of flow boiling of liquid nitrogen in a vertical mini-tube

The flow boiling patterns of liquid nitrogen in a vertical mini-tube with an inner diameter of 1.931 mm are visualized with a high-speed digital camera. The superficial gas and liquid velocities are in the ranges of 0.01–26.5 m/s and 0.01–1.2 m/s, respectively. Four typical flow patterns, namely, bubbly, slug, churn and annular flow are observed. Some interesting scenes about the entrainment and liquid droplet deposition in the churn and annular flow, and the flow reversal with the indication of negative pressure drop, are also presented. Based on the visualization, the two-phase flow regime maps are obtained. Compared with the flow regime maps for gas–water flow in tubes with similar hydraulic diameters, the region of slug flow in the present study reduces significantly. Correspondingly, the transition boundary from the bubbly flow to slug flow shifts to higher superficial gas velocity, and that of churn to annular flow moves to lower superficial gas velocity. Moreover, time-averaged void fraction is calculated by quantitative image-digitizing technique and compared with various prediction models. Finally, three kinds of oscillations with long-period and large-amplitude are found, possible explanation for the oscillations is given by comparing the instantaneous flow images with the data of pressure, mass flux and temperature recorded synchronously.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

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.     

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.     

Two-phase damping and interface surface area in tubes with vertical internal flow

Two-phase flow is common in the nuclear industry. It is a potential source of vibration in piping systems. In this paper, two-phase damping in the bubbly flow regime is related to the interface surface area and, therefore, to flow configuration. Experiments were performed with a vertical tube clamped at both ends. First, gas bubbles of controlled geometry were simulated with glass spheres let to settle in stagnant water. Second, air was injected in stagnant alcohol to generate a uniform and measurable bubble flow. In both cases, the two-phase damping ratio is correlated to the number of bubbles (or spheres). Two-phase damping is directly related to the interface surface area, based on a spherical bubble model. Further experiments were carried out on tubes with internal two-phase air–water flows. A strong dependence of two-phase damping on flow parameters in the bubbly flow regime is observed. A series of photographs attests to the fact that two-phase damping in bubbly flow increases for a larger number of bubbles, and for smaller bubbles. It is highest immediately prior to the transition from bubbly flow to slug or churn flow regimes. Beyond the transition, damping decreases. It is also shown that two-phase damping increases with the tube diameter.     

On the development of an objective flow regime indicator

A high intensity dual beam X-ray system was designed and constructed to make chordal-average void fraction measurements. This X-ray system employed a DC excited tube filament, full wave rectification and high voltage filtering to produce a stable photon source. The large photon flux produced by the X-ray system allowed analog linearization of the signal.A series of chordal-average void fraction measurements were made and used to generate probability density functions (PDF) and power spectral density (PSD) functions. The first four moments associated with these distributions were studied as possible flow regime indicators.The moments of the PDF indicated the various flow regime transitions. The moments of the PSD also show some flow regime transition information, but were sensitive to liquid phase velocity. The PDF variance, or second moment about the mean, was found to be the best indicator of flow regime. A variance of 0.04 appear to adequately discriminate between the bubbly, slug and annular flow regimes for low pressure air/water flow in a 2.54 cm I.D. vertical tube.     

Liquid-phase axial mixing in two-phase horizontal pipe flow

The liquid-phase axial dispersion coefficient and volume-averaged fractional phase hold-ups have been measured in two-phase horizontal pipe flow. Radioactive ^99mTc—technetium-99 metastable—(as an aqueous solution of sodium pertechnate) was used as a tracer. The pulse technique with two-point measurement was employed. Superficial gas (air) and liquid (water) velocities were varied in the range 20–2300 and 30–800mm/s, respectively. The flow regimes covered were bubbly, elongated bubbly, stratified, wavy and slug. Experiments were also performed using single-phase pipe flow. The liquid-phase dispersion coefficient has been shown to depend upon the flow regime and the superficial gas and liquid velocities.     

Spatial distribution of void fraction within a liquid slug and some other related slug parameters

The structure of vertical upward slug flow in a pipe is studied. The distribution of the phases in the Taylor bubble zone and the liquid slug zone is investigated by simultaneous measurements with two optical fiber probes. In the Taylor bubble zone the shape of the Taylor bubble and the distribution of the bubble length is reported. In the liquid slug region, the distribution of the void fraction is obtained over a dense grid in both the axial and radial directions. These experimental results shed some light on the hydrodynamics of the two-phase slug flow, in particular regarding the production of the dispersed bubbles and their distribution along the liquid slug.     

Microgravity flow regime data and analysis

To utilize the advantageous properties of two-phase flow in microgravity applications, the knowledge base of two-phase flow phenomena must be extended to include the effects of gravity. In the experiment described, data regarding the behavior of two-phase flow in a conduit under microgravity conditions (essentially zero gravity) are explored. Of particular interest, knowledge of the void fraction of the gas and liquid in a conduit is necessary to develop models for heat and mass transfer, pressure drop, and wall shear. An experiment was conducted under reduced gravity conditions to collect data by means of a capacitance void fraction sensor and high speed visual imagery. Independent parameters were varied to map the flow regime regions. These independent parameters include gas and liquid volumetric flow rates and saturation pressures. Void fraction measurements were taken at a rate of 100 Hz with six sensors at two locations along the conduit. Further, statistical parameters were developed from the void fraction measurements. Statistical parameters such as variance, signal-to-noise ratio, half height value, and linear area difference were calculated and found to have characteristics allowing flow regime identification.     

On two-phase flow patterns and transition criteria in aqueous methanol and CO2 mixtures in adiabatic,rectangular microchannels

This paper presents flow map investigations of adiabatic two-phase flow in square cross-sectioned, 200 μm deep microchannels fabricated in silicon, employing laser induced fluorescence microscopy. The influence of surface tension and nozzle geometry on the flow pattern transition was investigated using two nozzle widths (orifices of 30 μm and 50 μm, respectively) and methanol–water solutions with CO₂ as the gas phase. It was found and quantified that smaller nozzle geometries and smaller liquid surface tension promote the propagation of capillary gas bubbles at lower superficial gas and liquid velocities. Within the measurement domain of superficial gas (0.01–0.625 m/s) and liquid (0.0005–0.5000 m/s) velocities, we observed dispersed bubbly, regularly ordered bubbly, wedging, slug and annular flows, thus extending the experimental knowledge base to smaller superficial liquid velocities by almost two orders of magnitude. With the help of the flow maps presented herein, we were able to characterize the observed regularly ordered bubbly flow as the transition regime between dispersed bubbly and wedging flow. The results of the present investigation are of direct relevance to the operation of small-scale direct methanol fuel cells.     

Experimental study of “laminar” bubbly flows in a vertical pipe

Measurement of bubbly two-phase flow parameters in a vertical pipe were performed. To keep the pipe Reynolds number below that for single-phase turbulent transition, a water-glycerin solution was used as the test liquid. Local void fraction and liquid velocity profiles along with the wall shear stress were measured by an electrochemical method. Experiments were made with bubbles of two different sizes. As the gas flow rate was increased, a gradual development of the liquid velocity profile from the parabolic Poiseuille flow to a flattened two-phase profile was observed. The evolution of the wall shear stress and of the velocity fluctuations were also quantified.Centre National de la Recherche Scientifique. Université Joseph Fourier, Institut National Polytechnique de Grenoble.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

Void fraction measurement of bubble and slug flow in a small channel using the multivision technique

Using the multivision technique, a new void fraction measurement method was developed for bubble and slug flow in a small channel. The multivision system was developed to obtain images of the two-phase flow in two perpendicular directions. The obtained images were processed—using image segmentation, image subtraction, Canny edge detection, binarization, and hole filling—to extract the phase boundaries and information about the bubble or slug parameters. With the extracted information, a new void fraction measurement model was developed and used to determine the void fraction of the two-phase flow. The proposed method was validated experimentally in horizontal and vertical channels with different inner diameters of 2.1, 2.9, and 4.0 mm. The proposed method of measuring the void fraction has better performance than the methods that use images acquired in only one direction, with a maximum absolute difference between the measured and reference values of less than 6%.     

Experimental investigation of horizontal forced-vibration effect on air-water two-phase flow

In order to investigate the potential seismic vibrations effect on two-phase flow in an annular channel, experimental tests with air-water two-phase flow under horizontal vibrations were carried out. A low-speed eccentric-cam vibration module capable of operating at motor speed of 45–1200 rpm (f = 0.75–20 Hz) was attached to an annular channel, which was scaled down from a prototypic BWR fuel sub-channel with inner and outer diameters of 19.1 mm and 38.1 mm, respectively. The two-phase flow was operated in the ranges of 〈j_f〉 = 0.25–1.00 m/s and 〈j_g〉 = 0.03–1.46 m/s with 27 flow conditions, and the vibration amplitudes controlled by cam eccentricity (E) were designed for the range of 0.8–22.2 mm. Ring-type impedance void meters were utilized to detect the area-averaged time-averaged void fraction under stationary and vibration conditions. A systematic experimental database was built and analyzed with effective maps in terms of flow conditions (〈j_g〉-〈j_f〉) and vibration conditions (E-f and f-a), and the potential effects were expressed by regions on the maps. In the 〈j_g〉-〈j_f〉 maps, the void fraction was found to potentially decrease under vibrations in bubbly flow regime and relatively lower liquid flow conditions, which may be explained by the increase of distribution parameter. Whereas and the void fraction may increase at the region closed to bubbly-to-slug transition boundary under vibrations, which may be explained by the changes of drift velocity due to flow regime change from bubbly to slug flows. No significant change in void fraction was found in slug flow regime under the present test conditions.