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.     

Void fraction and pressure drop during external upward two-phase crossflow in tube bundles – part I: Experimental investigation

The present paper is the part I of a broad study concerning void fraction and pressure drop for air-water upward external flow across tube bundles. Experimental results were obtained for liquid and gas superficial velocities ranging from 0.02 to 1.50 m/s and 0.20 to 10.00 m/s, respectively. Void fraction measurements were performed for bubbly flow using a capacitive probe. The test section consisted of a triangular tube bundle counting with 19 mm OD tube and transverse pitch of 24 mm. Initially, the paper describes the test facility and the data regression and experimental procedures. Then, the pressure drop and void fraction measurements are validated based on tests for single-phase flow and quiescent liquid conditions, respectively. Finally, the experimental data are presented and analyzed. In the second part of this study (Part II), a literature review on predictive methods for void fraction and pressure drop is presented. Additionally, these methods are compared with the database presented in Part I and new predictive methods for void fraction and frictional pressure drop are proposed.     

Mass flow rate measurements in gas–liquid flows by means of a venturi or orifice plate coupled to a void fraction sensor

Two-phase flow measurements were carried out using a resistive void fraction meter coupled to a venturi or orifice plate. The measurement system used to estimate the liquid and gas mass flow rates was evaluated using an air–water experimental facility. Experiments included upward vertical and horizontal flow, annular, bubbly, churn and slug patterns, void fraction ranging from 2% to 85%, water flow rate up to 4000 kg/h, air flow rate up to 50 kg/h, and quality up to almost 10%. The fractional root mean square (RMS) deviation of the two-phase mass flow rate in upward vertical flow through a venturi plate is 6.8% using the correlation of Chisholm (D. Chisholm, Pressure gradients during the flow of incompressible two-phase mixtures through pipes, venturis and orifice plates, British Chemical Engineering 12 (9) (1967) 454–457). For the orifice plate, the RMS deviation of the vertical flow is 5.5% using the correlation of Zhang et al. (H.J. Zhang, W.T. Yue, Z.Y. Huang, Investigation of oil–air two-phase mass flow rate measurement using venturi and void fraction sensor, Journal of Zhejiang University Science 6A (6) (2005) 601–606). The results show that the flow direction has no significant influence on the meters in relation to the pressure drop in the experimental operation range. Quality and slip ratio analyses were also performed. The results show a mean slip ratio lower than 1.1, when bubbly and slug flow patterns are encountered for mean void fractions lower than 70%.     

One-group interfacial area transport equation and its sink and source terms in narrow rectangular channel

The characteristics of two-phase flow in a narrow rectangular channel are expected to be different from those in other channel geometries, because of the significant restriction of the bubble shape which, consequently, may affect the heat removal by boiling under various operating conditions. The objective of this study is to develop an interfacial area transport equation with the sink and source terms being properly modeled for the gas–liquid two-phase flow in a narrow rectangular channel. By taking into account the crushed characteristics of the bubbles a new one-group interfacial area transport equation was derived for the two-phase flow in a narrow rectangular channel. The random collisions between bubbles and the impacts of turbulent eddies with bubbles were modeled for the bubble coalescence and breakup respectively in the two-phase flow in a narrow rectangular channel. The newly-developed one-group interfacial area transport equation with the derived sink and source terms was evaluated by using the area-averaged flow parameters of vertical upwardly-moving adiabatic air–water two-phase flows measured in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm. The flow conditions of the data set covered spherical bubbly, crushed pancake bubbly, crushed cap-bubbly and crushed slug flow regimes and their superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. Good agreement with the average relative deviation of 9.98% was obtained between the predicted and measured interfacial area concentrations in this study.     

Effect of bubble size on void fraction fluctuations in dispersed bubble flows

It is known that bubble size affects seriously the average void fraction in bubbly flows where buoyant velocities vary considerably with bubble size. On the contrary, there is no systematic literature report about bubble size effects on the intensity and frequency of void fraction fluctuations around the average void fraction. This work aims to provide such information. An electrical impedance technique is employed along with non-intrusive ring electrodes to register void fraction fluctuations down to 10⁻⁵. Bubble size fluctuations are estimated from high resolution optical images. Experiments are conducted in co-current upward dispersed bubble flow inside a 21 mm tube with average bubble size between ∼50 and ∼700 μm. Water and blood simulant are used as test liquids with velocity from ∼3 to ∼30 cm s⁻¹. The above resemble conditions of Decompression Sickness (DCS) in the bloodstream of human vena cava. It is found that the intensity and frequency of void fraction fluctuations vary appreciably with bubble size at constant gas and liquid flow rates. Moreover, these variations are not random but scale with bubble size. As a first step to quantify this effect, an empirical expression is derived that relates average bubble size to the ratio standard deviation/average value of void fraction.     

Void fraction and pressure waves in a transient horizontal slug flow

The void fraction and the pressure waves in an air–water mixture flowing in the slug regime are experimentally investigated in a horizontal line. The test section is made of a transparent Plexiglas pipe with 26 mm ID and 26.24 m long, operating at ambient temperature and pressure. The flow induced transients are made by quickly changing the air or the water inlet velocity. The test grid has four operational points. This choice allows one to create expansion and compression waves due to the changes to the gas or to the liquid. Each experimental run is repeated 100 times to extract an ensemble average capable of filtering out the intrinsic flow intermittence and disclosing the void fraction and pressure waves’ features. The slug flow properties such as the bubble nose translational velocity, the lengths of liquid film underneath the bubble and the liquid slug are also measured. The objective of the work is two-fold: access the main characteristics of the void fraction and pressure waves and disclose the mechanics of the transient slug flow as described through the changes of the slug flow properties.     

Numerical research on micro diesel spray characteristics under ultra-high injection pressure by Large Eddy Simulation (LES)

The Large Eddy Simulation model was introduced to study the micro spray characteristics under ultra-high injection pressure (>220 MPa). EFS8400 spray test platform was set up to verify the accuracy of the numerical model. The mechanisms of micro spray characteristics were studied intensively under different injection pressures (180 MPa, 240 MPa) and nozzle diameters (0.1 mm, 0.16 mm). The results indicated that the micro turbulence vortex structures can be captured, especially in the liquid spray core area. Large Eddy Simulation model combined with the small grid size of 0.25 mm show a huge advantage in studying the micro spray characteristics under ultra-high injection pressure; The turbulence vorticity and spray velocity for injection pressure of 240 MPa are more intensive than that of 180 MPa, and also the ultra-high injection pressure can contribute to strong turbulence disturbance between spray and surrounding air, which is helpful to improve the quality of spray; The spray velocity field extended wider for the diameter of 0.16 mm, and also the values of velocity in the spray center is higher than that of the diameter of 0.1 mm; The entrainment vortex appeared at the edge of the large velocity gradient between spray and surrounding air, and the higher velocity gradient for ultra-high injection pressure (240 MPa) between the spray and air is easier to increase the generation of entrainment vortex in the downstream of the spray, which can significantly increase the quality of spray and atomization.     

Development of a robust image processing technique for bubbly flow measurement in a narrow rectangular channel

This paper presents a robust image processing technique for bubbly flow measurement over a wide range of void fractions. The proposed algorithm combines geometrical, optical and topological information recorded with high speed cameras to separate and reconstruct the overlapping bubbles. The common difficulties such as overlapping, irregular bubble shape, surface deformation and large clustering in digital image processing are solved by combining different information based on a preset decision table and flow chart. Test with synthetic bubble images is performed to evaluate the reliability of the algorithm and quantify the uncertainty of the data. The result shows that the proposed algorithm can accurately measure bubbly flows with void fraction up to 18% for large bubbles. Four runs of bubbly flow images in a 30 mm  ×  10 mm rectangular channel are then recorded by three high speed cameras. The area-averaged void fraction of these test runs range from 2.4% to 9.1%. The axial and lateral distributions of bubble number density are obtained by the present algorithm for studying the characteristics of these flows.     

Influence of sudden contractions on in situ volume fractions for oil–water flows in horizontal pipes

Oil–water two-phase flow experiments were conducted in horizontal ducts made of Plexiglas® to determine the in situ oil fraction (holdup) by means of the closing valves technique, using mineral oil (viscosity: 0.838 Pa s at 20 °C; density: 890 kg m⁻³) and tap water. The ducts present sudden contractions from 50 mm to 40 mm i.d. and from 50 mm to 30 mm i.d., with contraction ratios of 0.64 and 0.36, respectively. About 200–320 tests were performed by varying the flow rates of the phases. Flow patterns were investigated for both the up- and downstream pipe in order to assess whether relevant variations of the flow patterns across the sudden contraction take place. Data were then compared with predictions of a specific correlation for oil–water flow and some correlations for gas–water flow. A drift-flux model was also applied to determine the distribution parameter.     

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.     

Bubble–wall interaction and two-phase flow parameters on a full-scale boat boundary layer

Full scale bubbly flow experiments were performed on a 6 m flat bottom survey boat, measuring the void fraction, bubble velocity and size distributions as the bubbles naturally entrained at the bow of the boat interact with the boat’s boundary layer. Double-tip sapphire optical probes capable of measuring bubbles down to 50 μm in diameter were specifically designed and built for this experiment. The probes were positioned under the hull at the bow near the bubble entrainment region and at the stern at the exit of the bottom flat plate. Motorized positioners were used to vary the probe distance to the wall from 0 to 50 mm. The experiments were performed in fresh water (Coralville Lake, IA) and salt water (Panama City Beach, FL), at varying velocities with most data analysis performed at 10, 14 and 18 knots. The results indicate that the bubbles interact significantly with the boundary layer. At low velocity in fresh water, bubble accumulation under the hull and coalescence are evident by the presence of large bubbles at the stern. At high speeds bubble breakup dominates and very small bubbles are produced near the wall. It is also observed that salt water inhibits coalescence, even at low boat speeds. The void fraction increases with speed beyond 10 knots and peaks near the wall. Bubble velocities show slip with the wall at all speeds and exhibit large RMS fluctuations, increasing near the wall.     

Experimental study on local characteristics of oil–water dispersed flow in a vertical pipe

The local flow characteristics of oil–water dispersed flow in a vertical upward pipe were studied experimentally. The inner diameter and length of the test section are 40 mm and 3800 mm, respectively. A double-sensor conductivity probe was used to measure the local interfacial parameters, including interfacial area concentration, oil phase fraction, interfacial velocity, and oil drops Sauter mean diameter. The water flow rates varied from 0.12 m/s to 0.89 m/s, while the oil flow rates ranged from 0.024 m/s to 0.198 m/s. Typical radial profiles of interfacial area concentration, oil phase fraction, interfacial velocity, and oil drops Sauter mean diameter are presented. An interesting phenomenon is that the local and cross-section-averaged interfacial area concentrations display concave change with water flow rate under constant oil flow rate. The physical mechanism of such a variation is discussed in details.     

Experimental study of lean flammability limits of methane/hydrogen/air mixtures in tubes of different diameters

Lean limit flames in methane/hydrogen/air mixtures propagating in tubes of internal diameters (ID) of 6.0, 8.9, 12.3, 18.4, 25.2, 35.0, and 50.2 mm have been experimentally studied. The flames propagated upward from the open bottom end of the tube to the closed upper end. The content of hydrogen in the fuel gas has been varied in the range 0–40 mol%. Lean flammability limits have been determined; flame shapes recorded and the visible speed of flame propagation measured. Most of the observed limit flames in tubes with diameters in the range of 8.9–18.4 mm had enclosed shape, and could be characterized as distorted or spherical flame balls. The tendency was observed for mixtures with higher hydrogen content to form smaller size, more uniform flame balls in a wider range of tube diameters. At hydrogen content of 20% or more in the fuel gas, limit flames in largest diameters (35.0 mm and 50.2 mm ID) tubes had small, compared to the tube diameter, size and were “lens”-shaped. “Regular” open-front lean limit flames were observed only for the smallest diameters (6.0 mm and 8.9 mm) and largest diameters (35.0 and 50.2 mm ID), and only for methane/air and (90% CH₄ + 10% H₂)/air mixtures, except for 6 mm ID tube in which all limit flames had open front. In all experiments, except for the lean limit flames in methane/air and (90% CH₄ + 10% H₂)/air mixtures in the 8.9 mm ID tube, and all limit flames in 6.0 mm ID tube, visible flame speeds very weakly depended on the hydrogen content in the fuel gas and were close to- or below the theoretical estimate of the speed of a rising hot bubble. This observation suggests that the buoyancy is the major factor which determines the visible flame speed for studied limit flames, except that last mentioned. A decrease of the lean flammability limit value with decreasing the tube diameter was observed for methane/air and (90% CH₄ + 10% H₂)/air mixtures for tubes having internal diameters in the range of 18.4–50.2 mm. This effect has been attributed to the stronger combined effect of the preferential diffusion and flame stretch in narrower tubes for flames which resemble rising bubble.     

Flow pattern change in horizontal rectangular laterally ribbed ducts through alteration of the ribs thickness and pitch

Experiments were conducted on two-phase flow in laterally ribbed rectangular ducts. Air–water adiabatic flow at atmospheric pressure and room temperature was driven through a 3.6 m long rectangular ribbed test section with cross-section of 100 × 50 mm. To investigate the effect of rib thickness and pitch on flow pattern diagrams and transition boundaries, nine various rib arrangements were implemented with thicknesses of 2, 4 and 8 mm and pitches of 50, 60, and 80 mm. Unlike non-ribbed rectangular duct, lateral rib arrangement did not allow any stratified flow to occur. However wavy, plug and slug flows were parallel in both flow conditions, rib existence caused explicitly coarser pattern shapes. Increasing the rib thickness, while keeping the pitch constant, results in different flow patterns to occur as well as dramatic changes in boundaries positions and shapes. On the other hand, as pitch shifts up at a constant rib thickness, one can notice the duplication of almost identical flow patterns and their boundaries however, boundary values undergo tangible changes. Consistent attention was paid to conditions under which wavy pattern zone extends while intermittent flow zones were avoided. Studies concerned ribbed duct are of major applicable value to designing and enhancing heat transfer systems.     

Gas temperature in the trace of water droplets streamlined by hot air flow

In order to obtain the knowledge necessary for developing new effective fire extinguishing technologies, we determined experimentally the gas temperature in the trace of water droplets streamlined by hot air flow. It was important to establish how much the temperature in the droplet trace decreases and how fast it recovery to the initial temperature field after the droplet evaporation. The following parameters were varied: droplet size from 1.3 mm to 1.7 mm, velocity from 1 m/s to 5 m/s, initial airflow temperature from 473 K to 773 K, number of droplets (one or two), and the arrangement of droplets relative to the hot inflow (serial or parallel). The study proves the theoretical hypothesis about a significant influence of evaporation on the temperature in the water droplet trace. When a temperature trace of water droplets is formed, irrespective of their arrangement, the role of the evaporation process strengthens with the gas flow temperature rising. Furthermore, the study specifies typical longitudinal dimensions of the aerodynamic and temperature traces of water droplets. It has been established that when droplets are located in series and in parallel, their combined impact on the temperature and velocity of the gas flow in the medium differs rather considerably.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Tandem fluidized bed elutriator—Pneumatic classification of coal particles in a fluidized conveyer

Coal moisture control (CMC) in coking process, which reduces coal moisture before loading the coal into the coke oven, allows substantial reduction in coking energy consumption and increase in coke productivity. The technology is seeking to integrate the coal classification, thus calling it the coal classifying moisture control (CCMC), to separate the fine and coarse coal fractions in the CMC process so that the downstream coal crushing can only treat the coarse fraction. CCMC adopts a reactor that integrates a fluidized bottom section and a pneumatic conveying top section. The present work investigates the pneumatic classification behavior in a laboratory CCMC reactor with such a configuration by removing the coal fraction below a given size (e.g., 3.0 mm) from a 0 to 20.0 mm coal feed. The results show that the coal classification were dominated by the gas velocity in the top conveying section, and the required gas velocity for ensuring the maximal degree of removing a fine coal fraction could be roughly predicted by the Richardson and Zaki equation. The effect of bottom fluidization on the performance of CCMC is also examined.     

An experimental investigation of two-phase air–water flow through a horizontal circular micro-channel

This paper is a continuation of the authors’ previous work. Two-phase air–water flow experiments are performed in a horizontal circular micro-channel. The test section is made of a fused silica tube with an inner diameter of 0.15 mm and a length of 104 mm. The flow phenomena, which are liquid/unstable annular alternating flow (LUAAF), liquid/annular alternating flow (LAAF), and annular flow, are observed and recorded by a high-speed camera mounted together with a stereozoom microscope. A flow pattern map is presented in terms of the phase superficial velocities and is compared with those of other researchers obtained from different working fluids. Image analysis is performed to determine the void fraction, which increases non-linearly with increasing volumetric quality. It is revealed that the two-phase frictional multiplier data show a dependence on flow pattern rather than mass flux. Based on the present data, a new pressure drop correlation is proposed for practical applications. According to the present study, in general the data for the two-phase air–water flow characteristics are found to comply with those of working fluids other than air–water mixtures.     

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.