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.     

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.     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

Study of two-phase bubbly flow dynamics in binary mixtures

The results of an investigation of the bubble dynamics and wave phenomena in two-component vapor-liquid mixtures are presented. These mixtures are widely used in industrial systems as heat transfer media. The effects of various additives on the wave dynamics of vapor-liquid mixtures are of particular interest. A single-velocity two-pressure model was used which takes into account both the liquid radial inertia due to medium volume changes, and the temperature distribution around the bubbles. The claim that mixture composition may have a peculiar effect on the bubble dynamics of a boiling non-ideal solution is substantiated. It is noted that the small free radial oscillation damping ratio for some binary systems lies outside the domain defined by the damping ratio of the constituents as a result of phase change diffusion effects. A criterion is proposed to identify cases of diffusion resistance responsible for the anomalous effect of component concentration on bubble behavior. The phase change delay due to diffusion results in observably higher mixture wave velocities and a smaller damping ratio than for respective single-component systems.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1992 and the AMS fonts, developed by the American Mathematical Society.     

A USM-Θ two-phase turbulence model for simulating dense gas-particle flows

A second-order moment two-phase turbulence model for simulating dense gas-particle flows (USM- model), combining the unified second-order moment two-phase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM- model is also better than the k--kp- model and the k--kp-p- model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.The project supported by the Special Funds for Major State Basic Research of China (G-1999-0222-08), the National Natural Science Foundation of China (50376004), and Ph.D. Program Foundation, Ministry of Education of China (20030007028)     

Interfacial instability in turbulent flow over a liquid film in a channel

We revisit the stability of a deformable interface that separates a fully-developed turbulent gas flow from a thin layer of laminar liquid. Although this problem has received considerable attention previously, a model that requires no fitting parameters and that uses a base-state profile that has been validated against experiments is, as yet, unavailable. Furthermore, the significance of wave-induced perturbations in turbulent stresses remains unclear. To address these outstanding issues, we investigate this problem and introduce a turbulent base-state velocity that requires specification of a flow rate or a pressure drop only; no adjustable parameters are necessary. This base state is validated extensively against available experimental data as well as the results of direct numerical simulations. In addition, the effect of perturbations in the turbulent stress distributions is investigated, and demonstrated to be small for cases wherein the liquid layer is thin. The detailed modelling of the liquid layer also elicits two unstable modes, ‘interfacial’ and ‘internal’, with the former being the more dominant of the two. We show that it is possible for interfacial roughness to reduce the growth rate of the interfacial mode in relation to that of the internal one, promoting the latter, to the status of most dangerous mode. Additionally, we introduce an approximate measure to distinguish between ‘slow’ and ‘fast’ waves, the latter being the case for ‘critical-layer’-induced instabilities; we demonstrate that for the parameter ranges studied, the large majority of the waves are ‘slow’. Finally, comparisons of our linear stability predictions are made with experimental data in terms of critical parameters for onset of wave-formation, wave speeds and wavelengths; these yield agreement within the bounds of experimental error.     

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.     

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.     

An investigation of pressure transients in pipelines with two-phase bubbly flow

This paper presents a two-dimensional model for the analysis of the pressure transient of a two-phase homogeneous bubbly mixture flowing in a pipeline and the numerical integration using the centre implicit method (CIM). Experiments were conducted to confirm the proposed sonic speed equation of an air–water mixture for an air concentration of less than 1%. The 2D CIM model is compared with the method of characteristics (MoC) for a two-phase bubbly flow in a pipeline. The comparisons show that the proposed 2D CIM model generally gives good agreement with the method of characteristics.     

On the interaction between two consecutive elongated bubbles in a vertical pipe

The motion of elongated air bubbles in a vertical pipe filled with water is studied quantitatively using video imaging of the flow and subsequent digital image processing of the recorded sequence of images. Experiments are carried out to determine the influence of the separation distance between two consecutive bubbles (liquid slug length) upon the behavior of the trailing bubble in vertical slug flow. The details of the trailing bubble acceleration and merging process are observed and the instantaneous parameters of the trailing bubble, such as its shape, velocity, acceleration, etc., are measured as a function of the separation distance. The leading bubble is found to be unaffected by the trailing elongated bubble.     

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.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

Fluctuation characteristics of two-phase flow splitting at a vertical impacting T-junction

In order to investigate the fluctuation characteristics of two-phase flow splitting at a T-junction, particular attention was paid on Churn flow which had the strongest fluctuation comparing with bubble flow and annular flow. The main tube of the T-junction was vertical and the two branches were horizontal. All three pipes connecting to the junction were of 15 mm inner diameter. A statistical analysis based on Root Mean Square (RMS) was applied to temporal differential pressure signals and gas flow rate signals. The Power Spectral Density (PSD) was also employed to reveal their peculiar features in frequency domain as well. The effects of the extraction flow ratio and the gas and liquid superficial velocity upstream on fluctuation characteristics of gas-liquid two-phase flow splitting at the T-junction were investigated in detail. It is found that there is a wide fluctuation in both differential pressure and gas flow rate downstream at every extraction ratio (W₃/W₁) and the fluctuation intensity increases as W₃/W₁ increasing. It is also made clear that increasing either water superficial velocity or gas superficial velocity in inlet causes fluctuation to become more intensive.     

Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus

In this work radial and axial flow regime development in adiabatic upward air-water two-phase flow in a vertical annulus has been investigated. Local flow regimes have been identified using conductivity probes and neural networks techniques. The inner and outer diameters of the annulus are 19.1 mm and 38.1 mm, respectively. The equivalent hydraulic diameter of the flow channel, D_H, is 19.0 mm and the total length is 4.37 m. The flow regime map includes 1080 local flow regimes identifications in 72 flow conditions within a range of 0.01 m/s < 〈j_g〉 < 30 m/s and 0.2 m/s < 〈j_f〉 < 3.5 m/s where 〈j_g〉 and 〈j_f〉 are, respectively, superficial gas and liquid velocities. The local flow regime has been classified into four categories: bubbly, cap-slug, churn-turbulent and annular flows. In order to study the radial and axial development of flow regime the measurements have been performed at five radial locations. The three axial positions correspond to z/D_H = 52, 149 and 230, where z represents the axial position. The flow regime indicator has been chosen as some statistical parameters of local bubble chord length distributions and self-organized neural networks have been used as mapping system. This information has been also used to compare the results given by the existing flow regime transition models. The local flow regime is characterized basically by the void fraction and bubble chord length. The radial development of flow regime shows partial and complete local flow regime combinations. The radial development is controlled by axial location and superficial liquid velocity. The radial flow regime transition is always initiated in the center of the flow channel and it is propagated towards the channel boundaries. The axial development of flow regime is observed in all the flow maps and it is governed by superficial liquid velocity and radial location. The prediction results of the models are compared for each flow regime transition.     

Scaling of energy spectra in weakly compressible turbulence

We find an asymptotic expression for the char-acteristic timescales of decorrelation processes in weakly compressible and isothermal turbulence. This result is used in the Eddy-Damped Quasi-Normal Markovian equation to derive the scalings of compressible energy spectra: (1) if the acoustic waves are dominant, the compressible energy spectra exhibit ?7/3 scaling; (2) if local eddy straining is dominant, the compressible energy spectra are scaled as ?3. Meanwhile, the energy spectra of incompressible components display the same scaling of ?5/3 as those in incompressible turbulence. The direct numerical simulations of weakly compressible turbulence are used to examine the scaling.     

Effect of secondary flow on droplet distribution and deposition in horizontal annular pipe flow

In horizontal annular dispersed pipe flow the liquid film at the bottom is thicker and rougher than at the top of the pipe. A turbulent pipe flow experiencing a variation of roughness along the pipe wall will show a secondary flow. Such secondary flow, consisting of two counter-rotating cells in the cross-section of the tube, can change the distribution of the droplets inside the pipe and their deposition at the wall. Here, we compare the behaviour of the droplets (dispersed phase) with and without secondary flow, using large-eddy simulations. It is shown that the presence of secondary flow increases the droplet concentration in the core of the pipe and the droplet deposition-rate at the top of the pipe.     

Particle modulations to turbulence in two-phase round jets

The particle modulations to turbulence in round jets were experimentally studied by means of two-phase velocity measurements with Phase Doppler Anemometer （PDA）. Laden with very large particles, no significant attenuations of turbulence intensities were measured in the farfields, due to small two-phase slip velocities and particle Reynolds number. The gas-phase turbulence is enhanced by particles in the near-fields, but it is significantly attenuated by the small particles in the far-fields. The smaller particles have a more profound effect on the attenuation of turbulence intensities. The enhancements or attenuations of turbulence intensities in the far-fields depends on the energy production, transport and dissipation mechanisms between the two phases, which are determined by the particle prop- erties and two-phase velocity slips. The non-dimensional parameter CTI is introduced to represent the change of turbulence intensity.     

Local flow structure beyond bubbly flow in large diameter channels

A thorough review of the available literature has revealed a significant lack of usable data regarding the transport of interfacial area in large diameter channels. This represents a concern for various industrial systems, but especially for predicting the performance of safety systems in nuclear reactor systems. In order to remedy this gap in the current experimental database a series of experiments has been performed. These experiments included the measurement of the local interfacial area concentration and other parameters using local electrical conductivity probes in pipes with diameters of 0.152 m [6 in.], 0.203 m [8 in.] and 0.304 m [12 in.]. Volumetric fluxes ranged up to 2 m/s [6.56 ft/s] for the liquid phase and 10 m/s [32.8 ft/s] for the gas phase, and two nominal pressure conditions of 180 kPa [26.1 psia] and 280 kPa [40.6 psia] were included. Gas was injected as large cap bubbles in order to provide a basis for evaluating models for cap-bubbly flow at low void fractions. Measurements were performed simultaneously at three axial locations to allow the evaluation of interfacial area transport. The resulting data provides valuable insight into the flow structure and behavior in all flow regimes other than annular flow and will serve as a valuable database for the evaluation of models for predicting the transport of interfacial area across a wide variety of flow conditions and pipe sizes.