共查询到20条相似文献,搜索用时 31 毫秒
1.
《International Journal of Multiphase Flow》2005,31(7):824-842
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 (de ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (de ⩾ 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. 相似文献
2.
《International Journal of Multiphase Flow》1999,24(8):1383-1406
An experimental investigation on flow around an oscillating bubble and solid ellipsoid with a flat bottom was conducted. A single air bubble (equivalent diameter De=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 (De∼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. 相似文献
3.
《Particuology》2018
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. 相似文献
4.
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 CT = 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 JL = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond JG = 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. 相似文献
5.
Paulo R. de Souza Mendes Flávio H. Marchesini Priscilla R. Varges 《ournal of non Newtonian Fluid Mechanics》2011,166(17-18):1004-1011
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. 相似文献
6.
《International Journal of Multiphase Flow》2005,31(10-11):1063-1096
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. 相似文献
7.
Zhengyang Wang Shaozeng Sun Hao Chen Qigang Deng Guangbo Zhao Shaohua Wu 《Particuology》2009,7(6):483-490
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. 相似文献
8.
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 CS 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 CS ≃ 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. 相似文献
9.
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−11 to 2.52 × 10+01. 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. 相似文献
10.
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 dm3/min to 5.0 dm3/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 dm3/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. 相似文献
11.
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/m3 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. 相似文献
12.
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. 相似文献
13.
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−5), 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. 相似文献
14.
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. 相似文献
15.
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/m2s, heat fluxes from 2.3 kW/m2 to 250 kW/m2 at Tsat = 31 °C and ΔTsub 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. 相似文献
16.
《Experimental Thermal and Fluid Science》2009,33(1):132-139
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 × 104, 2 × 104 and 3 × 104) 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. 相似文献
17.
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. 相似文献
18.
Numerical investigation is made for three-dimensional fluid flow and convective heat transfer from an array of solid and perforated fins that are mounted on a flat plate. Incompressible air as working fluid is modeled using Navier–Stokes equations and RNG based k ? ? turbulent model is used to predict turbulent flow parameters. Temperature field inside the fins is obtained by solving Fourier’s conduction equation. The conjugate differential equations for both solid and gas phase are solved simultaneously by finite volume procedure using SIMPLE algorithm. Perforations such as small channels with square cross section are arranged streamwise along the fin’s length and their numbers varied from 1 to 3. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 104 to 4 × 104 based on the fin length and Prandtl number is taken Pr = 0.71. Numerical computations are validated with experimental studies of the previous investigators and good agreements were observed. Results show that fins with longitudinal pores, have remarkable heat transfer enhancement in addition to the considerable reduction in weight by comparison with solid fins. 相似文献
19.
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 (mf) of CMC in the solution is varied in the range 0.1% ⩽ mf ⩽ 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 mf 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 mf, 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 mf. 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. 相似文献
20.
A circular water jet (Re = 1.6 × 105; We = 8.8 × 103) plunging at shallow angles (θ ≈ 12.5°) into a quiescent pool is investigated computationally and experimentally. A surprising finding from the work is that cavities, of the order of jet diameter, are formed periodically in the impact location, even though the impinging flow is smooth and completely devoid of such a periodicity. Computational prediction of these frequencies was compared with experimental findings, yielding excellent agreement. The region in the vicinity of the impact is characterized by strong churning due to splashing and formation of air cavities. Measured velocity profiles indicate a concentration of momentum beneath the free surface slightly beyond the impact location (X/Dj ≈ 14), with a subsequent shift towards the free surface further downstream of this point (X/Dj ≈ 30). This shift is due primarily to the action of buoyancy on the cavity/bubble population. Comparisons of the mean velocity profile between simulations and experiments are performed, yielding good agreement, with the exception of the relatively small churning flow region. Further downstream (X/Dj ≳ 40), the flow develops mostly due to diffusion and the location of peak velocity coincides with the free surface. In this region, the free surface acts as an adiabatic boundary and restricts momentum diffusion, causing the peak velocity to occur at the free surface. 相似文献