Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

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.     

Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column

A new approach for simulating the formation of a froth layer in a slurry bubble column is proposed. Froth is considered a separate phase, comprised of a mixture of gas, liquid, and solid. The simulation was carried out using commercial flow simulation software (FIRE v2014) for particle sizes of 60–150 μm at solid concentrations of 0–40 vol%, and superficial gas velocities of 0.02–0.034 m/s in a slurry bubble column with a hydraulic diameter of 0.2 m and height of 1.2 m. Modelling calculations were conducted using a Eulerian–Eulerian multiphase approach with k–ε turbulence. The population balance equations for bubble breakup, bubble coalescence rate, and the interfacial exchange of mass and momentum were included in the computational fluid dynamics code by writing subroutines in Fortran to track the number density of different bubble sizes. Flow structure, radial gas holdup, and Sauter mean bubble diameter distributions at different column heights were predicted in the pulp zone, while froth volume fraction and density were predicted in the froth zone. The model was validated using available experimental data, and the predicted and experimental results showed reasonable agreement. To demonstrate the effect of increasing solid concentration on the coalescence rate, a solid-effect multiplier in the coalescence efficiency equation was used. The solid-effect multiplier decreased with increasing slurry concentration, causing an increase in bubble coalescence efficiency. A slight decrease in the coalescence efficiency was also observed owing to increasing particle size, which led to a decrease in Sauter mean bubble diameter. The froth volume fraction increased with solid concentration. These results provide an improved understanding of the dynamics of slurry bubble reactors in the presence of hydrophilic particles.     

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.     

Gravity-driven azimuthal flow of a layer of thixotropic fluid on the inner surface of a horizontal tube

The gravity-driven azimuthal flow of a layer of thixotropic paint on the inner surface of a horizontal tube is studied, considering surface tension effects. Using the lubrication theory, it was shown that a non-linear, fourth-order partial differential equation governs the time evolution of the paint layer thickness distribution along the azimuthal coordinate. Three parameters arise in the analysis, namely, the Bond number and two rheology-related parameters. The governing equation is integrated via a second-order accurate finite-difference scheme. The results showed that, in situations where the capillary force dominates (Bo < 1), displacement of the paint after application is very slow. For situations where the gravity force dominates (Bo > 1), an undulation on the interface arises near the tube bottom at sufficiently large times.     

Motion of singles bubbles moving under a slightly inclined surface through stationary liquids

The influence of the liquid properties on the dynamical bubble shape and on the bubble motion has been investigated for bubbles moving under a downward facing inclined surface. The Morton number Mo varied from 2.59 × 10⁻¹¹ to 2.52 × 10⁺⁰¹. The Bond number Bo covered the range from 10 to 150 and the surface inclination angle θ was varied from 2° to 6°. To cover the wide range of Mo, several liquids such as glycerine, propanediol, water and isopropanol were used. The results have shown that the relation Fr = Fr(Bo, Mo, θ) is not adequate to describe the bubble motion, where Fr is the terminal Froude number. The choice of the terminal Reynolds number Re as the dependent parameter, allowed the clarification of the role of the Morton number on the bubble motion. At a given Bond number, the bubble Reynolds number decreases monotonously with the Morton number. Furthermore, an empirical correlation Re = Re(Bo, Mo, θ) is given that can be readily used in the mathematical modelling of bubble laden flows under solids.     

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.     

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.     

Experimental investigation on flow asymmetry in solid entrance region of a square circulating fluidized bed

To study the influence of back feeding particles on gas-solid flow in the riser, this paper investigated the flow asymmetry in the solid entrance region of a fluidized bed by particle concentration/velocity measurements in a cold square circulating fluidized beds （CFB）. The pressure drop distribution along the riser and the saturation carrying capacity of gas for Geldart-B type particles were first analyzed. Under the condition of u0 = 4 m/s and Gs = 21 kg/（m^2 s）, the back feeding particles were found to penetrate the lean gas-solid flow near the entrance （rear） wall before reaching the opposite （front） wall, thus leading to a relatively denser region near the front wall in the bottom bed. Higher solid circulation rate （u0 =4 m/s, Gs = 33 kg/（m^2 s）） resulted in a higher particle concentration in the riser. However the back feeding particles with higher momentum increased the asymmetry of the particle concentration/velocity profile in the solid entrance region. Lower air velocity （u0 =3.2 m/s） and Gs =21 kg/（m2 s）, beyond the saturation carrying capacity of gas, induced an S-shaped axial solid distribution with a denser bottom zone. This limited the penetration of the back feeding particles and forced the flnidizing air to flow in the central region, thus leading to a higher solid holdup near the rear wall. Under the conditions of uo = 4 m/s and Gs = 21 kg/（m^2 s）, addition of coarse particles （dp= 1145 μm） into the bed made the radial distribution of solids more symmetrical.     

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.     

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.     

Oil mist transport process in a long pipeline on turbulent flow transition region

Internal gas velocity fluctuations and their effects on the mist diffusion process were examined in a long horizontal pipe to understand oil mist transportation, particularly in the laminar-to-turbulent flow transition region. Three hot-wire anemometers and aerosol concentration monitors were used to deduce these effects as the two-phase mist flow gradually developed in the stream-wise direction. We found significant axial mist diffusion at Reynolds numbers (Re) < 1000 because of passive scalar transport by Poiseuille flow. However, this diffusion was restricted by the non-zero inertia of the mist at a Stokes number, O(10⁻⁵), relying on the Brownian motion of the mist. At Re > 2400, a sharp mist waveform was maintained by a turbulent flow with active radial mixing. New data were obtained within the range of 1000 < Re < 2400, which cannot be explained by interpolation between the above-mentioned two states. The mist concentration displays multiple temporal peaks at Re < 2000 owing to perturbations of localized turbulence as well as radial anisotropy as being conveyed more than 2000-diameters in distance. This behavior is caused by intermittent disturbances induced by the pipe wall roughness, which sharply distorts the wall-aligned laminar mist layer left by parabolic axial stretching of local laminar flow.     

Oscillation and breakup of a bubble under forced vibration

Coupled shape oscillations and translational motion of an incompressible gas bubble in a vibrating liquid container is studied numerically. The bubble oscillation characteristics are mapped based on the bubble Bond number (Bo) and the ratio of the vibration amplitude of the container to the bubble diameter (A/D). At small Bo and A/D, the bubble oscillation is found to be linear with small amplitudes, and at large Bo and A/D, it is nonlinear and chaotic. This chaotic bubble oscillation is similar to those observed in two coupled nonlinear systems, here being the gas inside the bubble and its surrounding liquid. Further increases in the forcing, results in the bubble breakup due to large liquid inertia.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

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.     

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.     

Dynamics of a Taylor bubble in steady and pulsatile co-current flow of Newtonian and shear-thinning liquids in a vertical tube

A computational analysis is carried out to ascertain the effects of steady and pulsatile co-current flow, on the dynamics of an air bubble rising in a vertical tube containing water or a solution of Carboxymethylcellulose (CMC) in water. The mass fraction (m_f) of CMC in the solution is varied in the range 0.1% ⩽ m_f ⩽ 1% to accommodate zero-shear dynamic viscosities in the range 0.009–2.99 Pa-s. It was found that the transient and time-averaged velocities of Taylor bubbles are independent of the bubble size under both steady as well as pulsatile co-current flows. The lengths of the Taylor bubbles under the Newtonian conditions are found to be consistently greater than the corresponding shear-thinning non-Newtonian conditions for any given zero-shear dynamic viscosity of the liquid. In contrast to observations in stagnant liquid columns, an increase in the dynamic viscosity of the liquid (under Newtonian conditions) results in a concomitant increase in the bubble velocity, for any given co-current liquid velocity. In shear-thinning liquids, the change in the bubble velocity with an increase in m_f is found to be relatively greater at higher co-current liquid velocities. During pulsatile shear-thinning flows, distinct ripples are observed to occur on the bubble surface at higher values of m_f, the locations of which remain stationary with reference to the tube for any given pulsatile flow frequency, while the bubble propagated upwards. In such a pulsatile shear-thinning flow, a localised increase in dynamic viscosity is accompanied near each ripple, which results in a localised re-circulation region inside the bubble, unlike a single re-circulation region that occurs in Newtonian liquids, or shear-thinning liquids with low values of m_f. It is also seen that as compared to frequency, the amplitude of pulsatile flow has a greater influence on the oscillating characteristics of the rising Taylor bubble. The amplitude of oscillation in the bubble velocity increases with an increase in the CMC mass fraction, for any given value of pulsatile flow amplitude.     

Experimental study of parallel and inclined turbulent wall jets

An experimental study of the flow field in a two-dimensional wall jet has been conducted. All measurements were carried out using hot-wire anemometry. The experimental facility has a rectangular slot nozzle of high aspect ratio l/b = 100 (where l and b are the length and height slot, respectively). Mean velocities and Reynolds stresses were determined with three nozzle Reynolds numbers (Re = 1 × 10⁴, 2 × 10⁴ and 3 × 10⁴) and four different inclination angles between the wall and the flow velocity at the nozzle (β = 0°, 10°, 20° and 30°). Results indicate that all wall jets are self-preserving in the developed region. Normal to the wall two regions can be identified: one similar to a plane free jet and the other similar to a boundary layer. Downstream the interaction between these two regions creates a mixed or third region. The logarithmic region increases with the distance from the nozzle and with the Reynolds number. For the inclined wall jet, the spreading rate expressed in terms of jet half-width or maximum velocity decay with respect to the streamwise distance, asymptotes to a linear law. The streamwise locations where the jet becomes self-similar are farther from the exit than in parallel wall jet. The slope of both half-width and maximum velocity decay in the developed region are affected by both wall jet inclination angle and nozzle exit Reynolds number.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

Large-eddy simulation of thermal fatigue in a mixing tee

Temperature fluctuations occur due to thermal mixing of hot and cold streams in the T-junctions of the piping system in nuclear power plants, which may cause thermal fatigue of piping system. In this paper, three-dimensional, unsteady numerical simulations of coolant temperature fluctuations at a mixing T-junction of equal diameter pipes were performed using the large eddy simulation (LES) turbulent model. The experiments used in this paper to benchmark the simulations were performed by Hitachi Ltd. The calculated normalized mean temperatures and fluctuating temperatures are in good agreement with the measurements. The influence of the time-step ranging from 100 Hz to 1000 Hz on the numerical simulation results was explored. The simulation results indicate that all the results with different frequencies agree well with the experimental data. Finally, the attenuation of fluctuation of fluid temperature was also investigated. It is found that, drastic fluctuation occurs within the range of less than L/D = 4.0; the fluctuation of fluid temperature does not always attenuate from the pipe center to the wall due to the continuous generation of vortexes. At the top wall, the position of L/D = 1.5 has a minimum normalized mean temperature and a peak value of root-mean square temperature, whereas at the bottom wall, the position having the same characteristics is L/D = 2.0.