Effects of hydrophilic particles on bubbly flow in slurry bubble column

To investigate the effects of hydrophilic particles on slurry bubble flows in a bubble column, distributions of the local gas holdup and the bubble frequency are measured using an electric conductivity probe. Particles are made of silica and their diameter is 100 μm. The particle volumetric concentration C_S is varied from 0 to 0.40. The measured data imply that the presence of particles promotes bubble coalescence. The film drainage time for two coalescing bubbles in a quasi two-dimensional bubble flow in a small vessel is also measured to quantitatively evaluate the particle effect on coalescence. A particle-effect multiplier is introduced into a coalescence efficiency model by taking into account the data of film drainage time and is implemented into a multi-fluid model. The main conclusions obtained are as follows: (1) the local gas holdup and bubble frequency in slurry bubble flows decrease with increasing the particle concentration, (2) the hydrophilic particles enhance bubble coalescence and the enhancement saturates at C_S ≃ 0.45, (3) the particle effect on coalescence is well accounted for by introducing the particle-effect multiplier to the film drainage time, and (4) the multi-fluid model can give good predictions for the distribution of the local gas holdup in the slurry bubble column.     

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.     

Origin of regime transition to turbulent flow in bubble column: Orifice- and column-induced transitions

The possible events during bubble formation on an orifice were investigated using a rectangular bubble column (30 cm × 30 cm × 100 cm). The gas flow rate through a single orifice was adjusted from 0.1 dm³/min to 5.0 dm³/min covering a high flow rate regime. At the high gas flow rate, the bubble formation process was complicated by diverse events, such as wake effect, channeling, and orifice-induced turbulent flow. The detachment period could be used to discern the bubble formation steps because it was strongly affected by the above events. The bubble size distribution around the orifice was also analyzed to gain a clearer understanding of the bubble formation process. Above the rate of 3.0 dm³/min through a single orifice, the detachment period converged to a value of 25 ms irrespective of the orifice diameter. The bubble size distribution also showed little difference in this range of gas flow rate. This could be explained by the development of turbulent flow around the orifice. A 0.15 m in-diameter bubble column was tested to investigate the effect of orifice-induced turbulent flow on the regime transition in which the homogeneous flow regime is converted into the heterogeneous flow regime in the column. Obvious distinction between the orifice- and column-induced transitions was observed.     

The effect of aspect ratio in counter-current gas-liquid bubble columns: Experimental results and gas holdup correlations

It is generally admitted that the gas holdup is independent of the column dimensions and gas sparger design if three criteria are satisfied: the diameter of the bubble column is larger than 0.15 m, gas sparger openings are larger than 1–2 mm and the aspect ratio is larger than 5. This paper contributes to the existing discussion; in particular, the effect of the aspect ratio (within the range 1–15) in a counter-current gas-liquid bubble column has been experimentally studied and a new gas holdup correlation to estimate the influence of aspect ratio, operation mode and working fluid on the gas holdup has been proposed. The bubble column, equipped with a spider gas sparger, is 5.3 m in height, has an inner diameter of 0.24 m; gas superficial velocities in the range of 0.004–0.23 m/s have been considered, and, for the runs with water moving counter-currently to the gas phase, the liquid has been recirculated at a superficial velocity of −0.0846 m/s. Filtered air has been used as the gaseous phase in all the experiments, while the liquid phase has included tap water and different aqueous solutions of sodium chloride as electrolyte. Gas holdup measurements have been used to investigate the flow regime transitions and the global bubble column hydrodynamics. The counter-current mode has turned out to increase the gas holdup and destabilize the homogeneous flow regime; the presence of electrolytes has resulted in increasing the gas holdup and stabilizing the homogeneous flow regime; the aspect ratio, up to a critical value, has turned out to decrease the gas holdup and destabilize the homogeneous flow regime. The critical value of the aspect ratio ranged between 5 and 10, depending on the bubble column operation (i.e., batch or counter-current modes) and liquid phase properties. Since no correlation has been found in the literature that can correctly predict the gas holdup under the investigated conditions, a new scheme of gas holdup correlation has been proposed. Starting from considerations concerning the flow regime transition, corrective parameters are included in the gas holdup correlation to account for the effect of the changes introduced by the aspect ratio, operation mode and working fluid. The proposed correlation has been found to predict fairly well the present experimental data as well as previously published gas holdup data.     

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.     

Axisymmetric simulation of the interaction of a rising bubble with a rigid surface in viscous flow

The interaction between a rising deformable gas bubble and a solid wall in viscous liquids is investigated by direct numerical simulation via an arbitrary-Lagrangian–Eulerian (ALE) approach. The flow field is assumed to be axisymmetric. The bubble is driven by gravity only and the motion of the gas inside the bubble is neglected. Deformation of the bubble is tracked by a moving triangular mesh and the liquid motion is obtained by solving the Navier–Stokes equations in a finite element framework. To understand the mechanisms of bubble deformation as it interacts with the wall, the interaction process is studied as a function of two dimensionless parameters, namely, the Morton number (Mo) and Bond number (Bo). We study the range of Bo and Mo from (2, 6.5 × 10⁻⁶) to (16, 0.1). The film drainage process is also considered in this study. It is shown that the deformation of a bubble interacting with a solid wall can be classified into three modes depending on the values of Mo and Bo.     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

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.     

Turbulent flow past a bubble and an ellipsoid using shadow-image and PIV techniques

An experimental investigation on flow around an oscillating bubble and solid ellipsoid with a flat bottom was conducted. A single air bubble (equivalent diameter D_e=9.12 mm) was attached to a small disk (∼1 mm) at the end of a needle and suspended across a vertical square channel (100 mm) by wire wherein water flowed downward at a constant flowrate. The solid ellipsoid (D_e∼9.1 mm) was suspended across the square channel in the same manner. The equivalent diameter-based Reynolds and Eotvos number range, 1950<Re<2250 and 11<Eo<11.5, placed the bubble in the ‘wobbly’ regime while the flow in its wake was turbulent. A constant flowrate and one bubble size was used such that flow in the wake was turbulent. Velocity measurements of the flow field around the bubble or solid were made using a one CCD camera Digital Particle Image Velocimetry (DPIV) system enhanced by Laser Induced Fluorescence (LIF). The shape of the bubble or solid was simultaneously recorded along with the velocity using a second CCD camera and an Infrared Shadow Technique (IST). In this way both the flow-field and the boundary of the bubble (solid) were measured. The velocity vector plots of flow around and in the wake of a bubble/solid, supplemented by profiles and contours of the average and root-mean-square velocities, vorticity, Reynolds stress and turbulent kinetic energy, revealed differences in the wake flow structure behind a bubble and solid. One of the significant differences was in the inherent, oscillatory motion of the bubble which not only produced vorticity in the near-wake, but as a result of apparent vorticity stretching distributed the turbulent kinetic energy associated with this flow more uniformly on its wake, in contrast to the solid.     

CFD simulation of a gas–solid fluidized bed with two vertical jets

A computational fluid dynamics (CFD) model is used to investigate the hydrodynamics of a gas–solid fluidized bed with two vertical jets. Sand particles with a density of 2660 kg/m³ and a diameter of 5.0 × 10^?4 m are employed as the solid phase. Numerical computation is carried out in a 0.57 m × 1.00 m two-dimensional bed using a commercial CFD code, CFX 4.4, together with user-defined Fortran subroutines. The applicability of the CFD model is validated by predicting the bed pressure drop in a bubbling fluidized bed, and the jet detachment time and equivalent bubble diameter in a fluidized bed with a single jet. Subsequently, the model is used to explore the hydrodynamics of two vertical jets in a fluidized bed. The computational results reveal three flow patterns, isolated, merged and transitional jets, depending on the nozzle separation distance and jet gas velocity and influencing significantly the solid circulation pattern. The jet penetration depth is found to increase with increasing jet gas velocity, and can be predicted reasonably well by the correlations of Hong et al. (2003) for isolated jets and of Yang and Keairns (1979) for interacting jets.     

Experimental investigation of the influence of column scale,gas density and liquid properties on gas holdup in bubble columns

Measurements of gas holdups in bubble columns of 0.16, 0.30 and 0.33 m diameter were carried out. These columns were operated in co-current flow of gas and liquid phases and in semibatch mode. The column of 0.33 m diameter was operated at elevated pressures of up to 3.6 MPa. Nitrogen was employed as the gas phase and deionized water, aqueous solutions of ethanol and acetone and pure acetone and cumene as the liquid phase. The effects of differing liquid properties, gas density (due to elevated pressure), temperature, column diameter and superficial liquid velocity on gas holdup were studied. The gas holdup measurements were utilized by differential pressure measurements at different positions along the height of the bubble columns which allowed for the identification of axial gas holdup profiles. A decrease of gas holdup with increasing column diameter and an increase of gas holdup with increasing pressure was observed. The effect of a slightly decreasing gas holdup with increasing liquid velocity was found to exist at smaller column diameters. The use of organic solvents as the liquid phase resulted in a significant increase in gas holdup compared to deionized water. It is found that published gas holdup models are mostly unable to predict the results obtained in this study.     

Hot-film anemometry in bubbly flow I: bubble–probe interaction

The interaction between a bubble, which is rising in a descending water flow, and a hot-film anemometer was experimentally investigated using stereoscopic high-speed imaging. The mean downward water velocity varied from 0 up to 0.15 m/s, i.e., relatively low, allowing for an extended bubble–probe interaction. Moreover, the direction of the water causes the wake of the probe to play a role before the bubble touches the probe. The equivalent bubble radii were 0.4–2.8 mm and the bubble velocities relative to the probe ranged from 0.04 to 0.38 m/s. Image processing techniques were applied to reconstruct the bubbles’ path, shape, and orientation during the interaction process. As a result, three types of interactions were found, namely penetrating, bouncing, and splitting interactions. The image sequences were compared with the corresponding time series of the hot-film anemometer. From the time series the type of interaction cannot be deduced, at least not for the analyzed flow situation. Furthermore, we demonstrated that the residence time estimate from the hot-film data is systematically biased in our type of experiments. Finally, it was found that the velocity of a bubble may be altered considerably due to the interaction.     

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 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.     

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.     

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.     

UDV measurements of single bubble rising in a liquid metal Galinstan with a transverse magnetic field

Liquid metal is an important type of energy transport carrier in nuclear reactors, such as in accelerator-driven sub-critical systems, fusion reactors and spallation neutron source devices. It is necessary to conduct research for bubbles rising in a liquid metal under different magnetic field intensities. The Perspex container is positioned concentrically inside a transverse magnetic field, which provides a homogeneous DC longitudinal magnetic field that passes through the fluid district. The coils are supplied with maximum field strength of 1.97 T. The equivalent diameter of the bubble is 3.1–5.6 mm. The Ultrasonic Doppler Velocimetry (UDV) method is used to evaluate the internal flow velocity of opaque liquid metals. Research shows that the influence of the Lorenz force on the bubble ascension velocity is not simply positive or negative. The magnetic field inhibits the ascension velocity of small bubbles with diameters of 3.1 mm and 3.4 mm. The terminal velocity for large bubbles with diameters of 4.57 mm, 5.15 mm and 5.6 mm is higher under a weak magnetic field than without a magnetic field. The positive effect happens under strong magnetic intensity. The target is to obtain the hydro-dynamical relationships between the terminal velocity, drag coefficient, the Eötvös number, Reynolds number, and Stuart number in a strong magnetic field using a multiple regression method to reveal that the mechanism of the induced current's restraining influence determines the ascension velocity of the bubble in viscous electric liquids with a strong magnetic field.     

Vertical stability of bubble chain: Multiscale approach

Linear stability is investigated of a uniform chain of equal spherical gas bubbles rising vertically in unbounded stagnant liquid at Reynolds number Re = 50–200 and bubble spacing s > 2.6 bubble radii. The equilibrium bubble positions are questioned for their stability with respect to small displacements in the vertical direction, parallel to the chain motion. The transverse displacements are not considered, and the chain is assumed to be laterally stable. The bubbles are subjected to three kinds of forces: buoyant, viscous, inviscid. The viscous and inviscid forces have both pairwise (local) and distant (nonlocal) components. The pairwise forces are expressed by the leading-order formulas known from the literature. The distant forces are expressed as a linear superposition of the pairwise forces taken over several farther neighbours. The stability problem is addressed on three different length scales corresponding to: discrete chain (microscale), continuous chain (mesoscale), bubbly chain flow (macroscale). The relevant governing equations are derived for each scale. The microscale equations are a set of ODE’s, the Newton force laws for the individual discrete bubbles. The mesoscale equation is a PDE for bubbles continuously distributed along a line, obtained by taking the continuum limit of the microscale equations. The macroscale equations are two PDEs, the mass and momentum conservation equations, for an ensemble of noninteracting mesoscale chains rising in parallel. This transparent two-step process (micro → meso → macro) is an alternative to the usual one-step averaging, in obtaining the macroscale equations from microscale information. Here, the scale-up methodology is demonstrated for 1D motion of bubbles, but it can be used for behaviour of 2D and 3D lattices of bubbles, drops, and solids.It is found that the uniform equilibrium spacing results from a balance between the attractive and repulsive forces. On all three length scales, the equilibrium is stabilized by the viscous drag force, and destabilized by the viscous shielding force (shielding instability). The inviscid forces are stability neutral and generate conservative oscillations and concentration waves. The stability region in the parameter plane s − Re is determined for each length scale. The stable region is relatively small on the microscale, larger on the mesoscale, and shrinks to zero on the macroscale where the bubbly chain flow is inherently unstable.The shielding instability is expected to occur typically in intermediate Re flows where the vertical bubble interactions dominate over the horizontal interactions. This new kind of instability is studied here in a great detail, likely for the first time. Its relation to the elasticity properties of bubbly suspension on different length scales is discussed too. The shielding force takes the form of a negative bulk modulus of elasticity of the bubbly mixture.