One-dimensional interfacial area transport of vertical upward bubbly flow in narrow rectangular channel

The design and safety analysis for miniature heat exchangers, the cooling system of high performance microelectronics, research nuclear reactors, fusion reactors and the cooling system of the spallation neutron source targets requires the knowledge of the gas–liquid two-phase flow in a narrow rectangular channel. In this study, flow measurements of vertical upward air–water flows in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm were performed at seven axial locations by using the imaging processing technique. The local frictional pressure loss gradients were also measured by a differential pressure cell. In the experiment, the 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. The developing two-phase flow was characterized by the significant axial changes of the local flow parameters due to the bubble coalescence and breakup in the tested flow conditions. The existing two-phase frictional multiplier correlations such as Chisholm, 1967, Mishima et al., 1993 and Lee and Lee (2001) were verified to give a good prediction for the measured two-phase frictional multiplier. The predictions of the drift-flux model with the rectangular channel distribution parameter correlation of Ishii (1977) and several existing drift velocity correlations of Ishii, 1977, Hibiki and Ishii, 2003 and Jones and Zuber (1979) agreed well with the measured void fractions and gas velocities. The interfacial area concentration (IAC) model of Hibiki and Ishii (2002) was modified by taking the channel width as the system length scale and the modified IAC model could predict the IAC and Sauter mean diameter acceptably.     

Benchmarking of the one-dimensional one-group interfacial area transport equation for reduced-gravity bubbly flows

In order to properly design and safely operate two-phase flow systems, especially those deployed on future space missions, it is necessary to have accurate predictive capabilities. The application of a novel predictive method, the interfacial area transport equation (IATE), to dynamically predict the change of interfacial area concentration for reduced-gravity two-phase flows is described in this paper. Fluid particle interaction mechanisms such as coalescence and breakup that are present in reduced-gravity two-phase flows have been studied experimentally as reported in a previous paper by the current authors [Vasavada et al., 2007]. These mechanisms represent the source and sink terms in the IATE and their mechanistic models are benchmarked using experimental data obtained in a 25 mm inner diameter ground-based test section wherein reduced-gravity conditions were simulated. The comparison of the predictions from the model against experimental data shows good agreement. It has been found that, in contrast to the hypothesis extended in the literature, the wake entrainment based coalescence mechanism is present in reduced-gravity two-phase flows and in some cases is more important than coalescence due to random collision. Physics based arguments are extended to support this conclusion.     

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.     

Experimental study on near wall transport characteristics of slug flow in a vertical pipe

In this work, the wall shear stress and the mass transfer coefficient of the gas–liquid two-phase upward slug flow in a vertical pipe are investigated experimentally, using limiting diffusion current probes and digital high-speed video system. In experiments, the instantaneous and averaged characteristics of wall shear stress and mass transfer coefficient are concerned. The experimental results are compared with the numerical results in previous paper of the authors. Both experiment and numerical simulation show that the superficial gas and liquid velocities have an obvious influence on the instantaneous characteristics of the two profiles. The mass transfer coefficient has characteristics similar to the wall shear stress. The instantaneous wall shear stress and mass transfer coefficient profiles have the periodicity of slug flow. The averaged wall shear stress and mass transfer coefficient increase with increased superficial gas velocity. However, there is inconsistency in the variation trends of the averaged wall shear stress and mass transfer coefficient with superficial liquid velocity between experimental result and numerical simulation result, which can be attributed to the difference in flow condition. Moreover, the Taylor bubble length is also another impacting factor. The experimental and numerical results all shows that the product scale can not be damaged directly by the flow movement of slug flow. In fact, the alternative forces and fluctuations with high frequency acting on the pipe wall due to slug flow is the main cause for the slug flow enhanced CO₂ corrosion process.     

On the numerical study of bubbly flow created by ventilated cavity in vertical pipe

The dispersion of bubbles into a down-liquid flow in a vertical pipe is investigated. At low flow rates, the intended design of a swarm of discrete bubbles is achieved. At high flow rates, a ventilated cavity is nonetheless formed, which is attached close to the gas sparger. Behind this ventilated cavity, three different flow regimes characterize the complex bubbly flow field downstream of the down-liquid flow: vortex region with high void fraction, transitional region and pipe flow region. In this study, a numerical model that solved the entire development of the gas–liquid flow including the extended single-phase liquid region upstream to the wall-jet and recirculating-vortex zones in order to allow a more realistic determination of the boundary conditions of the down-liquid flow was adopted. Coupling with the Eulerian–Eulerian two-fluid model to solve the respective gas and liquid phases, a population balance model was also applied to predict the bubble size distribution in the wake right below the cavity base as well as further downstream in the transitional and fully-developed pipe flow regions. The numerical model was evaluated by comparing the numerical results against the data derived from theoretical, numerical and experimental approaches. Prediction of the Sauter mean bubble diameter distributions by the population balance approach at different axial locations confirmed the dominance of breakage due to the high turbulent intensity below the ventilated cavity which led to the generation of small gas bubbles at high void fraction. Further downstream, the coalescence effect dominated leading to merging of the small bubbles to form bigger bubbles.     

A DNS study of laminar bubbly flows in a vertical channel

Direct numerical simulations are used to examine laminar bubbly flows in vertical channels. For equal size nearly spherical bubbles the results show that at steady state the number density of bubbles in the center of the channel is always such that the fluid mixture there is in hydrostatic equilibrium. For upflow, excess bubbles are pushed to the walls, forming a bubble rich wall-layer, one bubble diameter thick. For downflow, bubbles are drawn into the channel center, leading to a wall-layer devoid of bubbles, of a thickness determined by how much the void fraction in the center of the channel must be increased to reach hydrostatic equilibrium. The void fraction profile can be predicted analytically using a very simple model and the model also gives the velocity profile for the downflow case. For the upflow, however, the velocity increase across the wall-layer must be obtained from the simulations. The slip velocity of the bubbles in the channel core and the velocity fluctuations are predicted reasonably well by results for homogeneous flows.     

Void fraction and pressure fluctuations of bubbly flow in a vertical annular channel

In this investigation some hydrodynamic characteristics of two phase, two component, air water bubbly flow in a vertical annulus were studied. In particular, the void fraction profiles, and the pressure fluctuations were measured by the electrical resistivity probe and a capacitive type differential transducer respectively. These measurements were assessed under various system parameters, viz the air and water flux, the perforation ratio (Area of holes/channel cross sectional area) and the dimensionless axial distance. In addition, the pressure drop calculated from the void fraction measurements was in very good agreement with the corresponding one measured by the pressure transducers.List of symbols D _eq equivalent diameter of the annular channel (m) - j flux (discharge/channel cross sectional area) (m/s) - m mass flow rate (kg/s) - P pressure (Pa) - AP static pressure difference along the test section (Pa) - P pressure fluctuations (Pa) - P ^* dimensionless pressure (P _m/P _S.P. ) - P dimensionless pressure fluctuations (P _max /P _T.P. ) - r radius (m) - z axial distance (m) Greek symbols void fraction - dimensionless axial distance (Z/Dimeq) - perforation ratio (area of holes/channel cross sectional area) - density (kg/m3) - time (s) - dimensionless radial distance (r–r _i )/(r _o-r _i ) Suffix g gas - i inner - L liquid - m mean - Max Maximum - O outer - S.P. single-phase - T.P. two-phase     

On thermodynamic closures for two-phase flow with interfacial area concentration transport equation

The present work is a part of a modelling of forest fires fighting by aerial means. In this paper, we study different kind of closures for modelling two-phase flows with an almost “infinite range” of scales. Since theories like homogenization are not, in this case, relevant for obtaining the equivalent medium equations, the averaging method has been preferred. The variables are averaged by convolution with a smooth kernel with compact support, as the equations are non-linear, new quantities are defined in order to obtain the equations satisfied by averaged quantities; the entropy production is determined and closures or phenomenological equations are obtained using the second principle of thermodynamics. Main features of this work are, firstly a derivation in this framework of a balance equation for the interfacial area concentration and secondly, since this introduces a new unclosed variable: the mean velocity of interfaces, extended irreversible thermodynamics is used to obtain the general form of the appropriate closures equations.     

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.     

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.     

Experimental study of gas-liquid slug flow in a small-diameter vertical pipe

Experimental investigation of upward gas-liquid slug flow in a vertical pipe in 15 mm ID has been carried out. The electrochemical method which permits the determination of the value and direction of instantaneous wall shear stress as well as the mean and fluctuating components of the liquid velocity is used for measurements. It is shown that the change of the sign of the velocity near the wall usually occurs at the moment of slug passage; the time-averaged wall shear stress at low liquid velocities is significantly lower than the value obtained by means of common prediction methods. The results of measuring of the local void fraction. liquid velocity and components of liquid velocity fluctuations are presented. The time-dependent behavior of the instantaneous hydrodynamic characteristics is described.     

Bubble and liquid turbulence characteristics of bubbly flow in a large diameter vertical pipe

The bubble and liquid turbulence characteristics of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. The bubble characteristics were measured using a dual optical probe, while the liquid-phase turbulence was measured using hot-film anemometry. Measurements were performed at six liquid superficial velocities in the range of 0.2–0.68 m/s and gas superficial velocity from 0.005 to 0.18 m/s, corresponding to an area average void fraction from 1.2% to 15.4%. At low void fraction flow, the radial void fraction distribution showed a wall peak which changed to a core peak profile as the void fraction was increased. The liquid average velocity and the turbulence intensities were less uniform in the core region of the pipe as the void fraction profile changed from a wall to a core peak. In general, there is an increase in the turbulence intensities when the bubbles are introduced into the flow. However, a turbulence suppression was observed close to the wall at high liquid superficial velocities for low void fractions up to about 1.6%. The net radial interfacial force on the bubbles was estimated from the momentum equations using the measured profiles. The radial migration of the bubbles in the core region of the pipe, which determines the shape of the void profile, was related to the balance between the turbulent dispersion and the lift forces. The ratio between these forces was characterized by a dimensionless group that includes the area averaged Eötvös number, slip ratio, and the ratio between the apparent added kinetic energy to the actual kinetic energy of the liquid. A non-dimensional map based on this dimensionless group and the force ratio is proposed to distinguish the conditions under which a wall or core peak void profile occurs in bubbly flows.     

Turbulence structure of air-water bubbly flow—III. transport properties

An experimental study of heat and bubble transport in turbulent air-water bubbly flow was carried out by means of tracer techniques. Helium tracer gas concentration data and temperature distributions were used to extract bubble and heat diffusivity information. The results indicated that the turbulent velocity components of the liquid phase play a predominant role in the turbulent transport process. A systematic increase of diffusivity of heat, ?_H, with quality and water velocity was observed. An empirical correlation for the diffusivity ratio ?_H,TP/?_H,SP is presented. The Péclet number, u′_covbar|dφ, for bubble dispersion can be approximated by 2.0, independent of the flow variables. The bubble-to-heat diffusivity ratio, φ/?_H, approaches unity with increasing quality and water velocity. Momentum transport is also discussed, based on a mixing length theory.     

Separated flow modelling and interfacial transport phenomena

A physical explanation of the different types of interfacial waves that appear in stratified and annular gas-liquid flows is presented. The role of waves in affecting process performance is discussed. Particular attention is paid to interfacial drag, gas absorption, the initiation of slug flow and atomization.This paper was presented at the Shell Conference on Computational Fluid Dynamics in Apeldoorn, December, 1989.     

Experimental investigation of heat transfer in vertical upward and downward supercritical CO2 flow in a circular tube

An experimental investigation of turbulent heat transfer in vertical upward and downward supercritical CO₂ flow was conducted in a circular tube with an inner diameter of 4.5 mm. The experiments were performed for bulk fluid temperatures from 29 to 115 °C, pressures from 74.6 to 102.6 bar, local wall heat fluxes from 38 to 234 kW/m², and mass fluxes from 208 to 874 kg/m² s. At a moderate wall heat flux and low mass flux, the wall temperature had a noticeable peak value for vertical upward flow, but increased monotonically along the flow direction without a peak value for downward flow. The ratios of the experimental Nusselt number to the value obtained from a reference correlation were compared with Bo^* and q⁺ distributions to observe the buoyancy and flow-acceleration effects on heat transfer. In the experimental range of this study, the flow acceleration predominantly affected the heat-transfer phenomena. Based on an analysis of the shear-stress distribution in the turbulent boundary layer and the significant variation of the specific heat across the turbulent boundary layer, a new heat-transfer correlation for vertical upward and downward flow of supercritical pressurized fluid was developed; this correlation agreed with various experimental datasets within ±30%.     

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.     

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.