Some peculiarities of two-phase flow in nozzles

Calculations are conducted for unidimensional two-phase flow in nozzles for a wide range of particle concentrations and dimensions. It is established that there exists a maximum in loss of specific momentum due to a lag in particle velocity and temperature relative to the gas. The results obtained are compared with calculations using linearized theory as well as with experimental data. The agreement between calculation and experiment is noted. Equilibrium flow of a two-phase mixture with solidification of liquid particles is considered. The presence of an anomalous flow region is established, where in the model of an ideal unidimensional equilibrium flow in a nozzle with discharge into a vacuum the presence of two successively located minimum sections is necessary.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–57, May–June, 1973.     

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.     

Dynamics of Brownian particles in a turbulent channel flow

Turbulent channel flows with suspended particles are investigated by means of numerical simulations. The fluid velocity is computed by large eddy simulation. Motion of small graphite particles with diameter of 0.01–10 m, corresponding to the Schmidt number, Sc, of 2.87 × 10²–6.22 × 10⁶ and the particle relaxation time in wall unit, _p⁺, of 9.79 × 10^–5–4.51, is computed by Lagrangian particle tracking. Relation between the particle relaxation time and the computed deposition velocity is found to be in good agreement with an empirical relation. The statistics of the particle motion in the vicinity of the wall are studied. Clear differences are found in dynamical behavior of particles with different sizes. Medium size particles show a strong dependence on the structure of the fluid flow, while small and large particles are considerably less sensitive.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

Drag models for simulating gas–solid flow in the turbulent fluidization of FCC particles

This paper examines the suitability of various drag models for predicting the hydrodynamics of the turbulent fluidization of FCC particles on the Fluent V6.2 platform. The drag models included those of Syamlal–O’Brien, Gidaspow, modified Syamlal–O’Brien, and McKeen. Comparison between experimental data and simulated results showed that the Syamlal–O’Brien, Gidaspow, and modified Syamlal–O’Brien drag models highly overestimated gas–solid momentum exchange and could not predict the formation of dense phase in the fluidized bed, while the McKeen drag model could not capture the dilute characteristics due to underestimation of drag force. The standard Gidaspow drag model was then modified by adopting the effective particle cluster diameter to account for particle clusters, which was, however, proved inapplicable for FCC particle turbulent fluidization. A four-zone drag model (dense phase, sub-dense phase, sub-dilute phase and dilute phase) was finally proposed to calculate the gas–solid exchange coefficient in the turbulent fluidization of FCC particles, and was validated by satisfactory agreement between prediction and experiment.     

Experimental and computational studies on flow behavior of gas–solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally. Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695-0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

Particle deposition on a channel wall in a turbulent disperse flow with gradient

Kinetic ideas about the motion of a set of particles (droplets) in a turbulent gas flow with gradient are used to derive a Fokker-Planck equation for the case of sufficiently large particles (more than few microns). This equation describes the process in which they are deposited on the wall of a channel. Satisfactory agreement has been obtained between the numerical solution to this equation for the deposition rate and the experimental data published in the literature. Under the assumption that the parameters of the carrier gaseous flow vary fairly slowly, a generalized equation is derived for particle diffusion in turbulent flow. This takes into account the intensity gradient of transverse pulsations in the velocity of the carrier gaseous flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 57–63, July–August, 1985.     

A Two-Dimensional Network Model to Simulate Permeability Decrease Under Hydrodynamic Effect of Particle Release and Capture

A model was developed to simulate permeability decrease induced by hydrodynamic effects when injecting a fluid in a reservoir with respect to particle release and capture mechanisms and the parameters of the fluid–rock system. The kinetics of particle release and capture were integrated after computing the initial permeability of the porous medium with a square lattice of a two–dimensional network model. The rate of particle release is related to the difference between a microscopic velocity of the fluid and a critical velocity. The permeability decrease shows a direct link to the reduction of pore throat radii by three mechanisms of particle capture: straining and particle accumulation through direct interception or diffusion. Comparison between the simulations and the experimental results shows that the model reproduces the physics of the permeability decrease phenomenon, although the values are overestimated. The difference between the two sets of results can be explained by the fact that the simulations are realized at constant pressure whereas the experiments are realized at constant flow rate, and that re–entrainment of the trapped particles was not taken into account in the model.     

Steady-state particle distributions in shear flows of colloids and finely dispersed suspensions

On the basis of an analysis of the pseudoturbulent motion of both the suspended particles and the carrier fluid, the normal stress components in the dispersed phase are obtained for the problem of inclined confined flows of finely dispersed suspensions and colloids. These hydrodynamic pulsations are due to the shear and the work done by the average relative flow of the fluid phase on random concentration fluctuations of the disperse system because of the substantial slip of the phases of the suspension under gravity. The momentum conservation equations for the particles are obtained with allowance for the angle of inclination of the flow to the vertical and on the basis of these equations the suspension capacity of the flow as a function of the angle of inclination, particle size, Galileo number and other parameters is illustrated.Ekaterinburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 78–84, January–February, 1996.     

Calculation of the momentum and heat transfer in turbulent gas-particle jet flows

The equations for the second moments of the dispersed-phase velocity and temperature fluctuations are used for calculating gas-suspension jet flows within the framework of the Euler approach. The advantages of introducing the equations for the second moments of the particle velocity fluctuations has previously been quite convincingly demonstrated with reference to the calculation of two-phase channel boundary flows [9–11]. The flows considered below have a low solid particle volume concentration, so that interparticle collisions can be neglected and, consequently, the stochastic motion of the particles is determined exclusively by their involvement in the fluctuating motion of the carrier flow. In addition to the equations for the turbulent energy of the gas and its dissipation, the calculation scheme includes the equations for the turbulent energy and turbulent heat transfer of the solid phase; however, the model constructed does not contain additional empirical constants associated with the presence of the particles in the flow.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.3, pp. 69–80, May–June, 1992.     

Effect of particle size on nonlinear wave propagation in a fluidized bed

The effect of particle size (Archimedes number) on the propagation of a kinematic particle concentration wave in a fluidized bed is investigated. The dependence of the characteristic wave velocity on the porosity of the bed (particle concentration) and the Archimedes number (or the Reynolds number for flow past individual particles of the dispersed phase) is determined. The evolution of a nonlinear perturbation of the bed porosity is investigated and the formation of discontinuities in the concentration of the dispersed phase is studied in relation to the particle size (Archimedes number). It is shown, in particular, that, as distinct from a bed of small particles, in a bed of large particles with quadratic interphase interaction only compression discontinuities can be formed. The results obtained can be used to analyze the formation of inhomogeneities (slugs and bubbles) in a fluidized bed in relation to the particle size.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 96–100, May–June, 1987.     

A particle suspension model: an unstructured finite-volume implementation

In this paper, a modified particle temperature model for concentrated suspensions is proposed, which allows for the shear-induced migration of particles. The migration is modelled by a convection–diffusion equation, derived from the particle mass and momentum conservation. The model is implemented in an unstructured finite volume method and is utilized to investigate the shear-induced particle migration in channel flow. The profiles and the evolution of the velocity, concentration and particle temperature along the channel are presented. The entrance lengths needed to reach a fully developed profile of the corresponding field variables are also checked against different averaged concentrations and different relative particle radii. Comparison with available experimental data is made whenever possible.     

Vibration induced flow in hoppers: DEM 2D polygon model

A two-dimensional discrete element model （DEM） simulation of cohesive polygonal particles has been developed to assess the benefit of point source vibration to induce flow in wedge-shaped hoppers. The particle-particle interaction model used is based on a multi-contact principle. The first part of the study investigated particle discharge under gravity without vibration to determine the critical orifice size （Bc） to just sustain flow as a function of particle shape. It is shown that polygonal-shaped particles need a larger orifice than circular particles. It is also shown that Bc decreases as the number of particle vertices increases. Addition of circular particles promotes flow of polygons in a linear manner. The second part of the study showed that vibration could enhance flow, effectively reducing Bc. The model demonstrated the importance of vibrator location （height）, consistent with previous continuum model results, and vibration amplitude in enhancing flow.     

Motion of inhomogeneities of a developed fluidized bed at small reynolds numbers

The instability of a fluidized system in which the particles are uniformly distributed in space [1–3] leads to the development of local inhomogeneities in the internal structure, these taking the form of more or less stable formations of packets of particles [4]. In accordance with the existing ideas based on experimental data [5–8, 13], the particle concentration within a packet may vary in a wide range from very small values (10^–2–10^–3 [8]) for bubbles to the concentration of the unfluidized bed for bunches of particles in a nearly closely packed state. The paper considers the steady disturbed motion of the fluid and solid phases near an ascending or descending packet of particles in a developed fluidized bed. It is assumed that the motion of the solid phase corresponds to a creeping flow of viscous fluid, and the viscosity of the fluidizing agent is taken into account only in the terms that describe the interphase interaction. The velocity fields and pressure distributions of the phases inside and outside a packet are determined. If the particle concentration within a packet tends to zero, the solution describes the slow motion of a bubble in a fluidized bed. The results of the paper are compared with results obtained earlier for the model of ideal fluids [9] and Batchelor's model [10], in which the fluidized bed is treated in a simplified form as a viscous quasihomogeneous continuum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 57–65, July–August, 1984.     

Charging of disperse particles in two-phase media with unipolar charge

Charging of disperse particles with good conduction in two-phase media with unipolar charge is considered in the case when the volume concentration of the particles is low. For this, in the framework of electrohydro-dynamics [1, 2], a study is made of the charge of one perfectly conducting liquid particle in a gas (or liquid) with unipolar charge in a fairly strong electric field. The influence of the inertial and electric forces on the motion of the gas is ignored, and the velocities are found by solving the Hadamard—Rybczynski problem. We consider the axisymmetric case when the gas velocity and electric field intensity far from the particle are parallel to a straight line. The analogous problem for a solid spherical particle was solved in [3–6] (in [3], the relative motion of the gas was ignored, while in [4–6] Stokes flow around the particle was considered). The two-dimensional problem of the charge of a solid circular, perfectly conducting cylinder in an irrotational flow of gas with unipolar charge was studied in [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 108–115, November–December, 1980.We thank L. I. Sedov and V. V. Gogosov for a helpful discussion of the present work.     

Equations of electrohydrodynamics of weakly ionized aerosols with diffusion charging of dispersed phase particles

A closed model describing the motion of weakly ionized aerosols with allowance for dispersed phase particle charging processes due to ion deposition is constructed within the framework of continuum mechanics [1]. Both the general process of particle charging in a weakly ionized gas and its limiting cases, where the limiting stage of the process is the diffusion of the ions towards the particles or the reactions leading to their deposition on the particles, are investigated. Expressions are obtained for the positive and negative ion flows to a particle in a weakly ionized gas. The basic equations of electrohydro-dynamics of weakly ionized aerosols, in which the dispersed phase particle charging mechanism in question leads to the interphase transfer of elctrical charge, are formulated. Cases where the system of equations of electrohydrodynamics obtained can be simplified by replacing the differential equations of motion and charging of the dispersed phase and, moreover, the positive and negative ion balance equations by algebraic relations such as Ohm's law and Saha's equation are investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 54–60, January–February, 1986.     

Gas–solid interaction force from direct numerical simulation (DNS) of binary systems with extreme diameter ratios

Fluid–particle systems as commonly encountered in chemical, metallurgical and petroleum industries are mostly polydisperse in nature. However, the relations used to describe fluid–particle interactions are originally derived from monodisperse systems, with ad hoc modifications to account for polydispersity. In previous work it was shown that for bidisperse systems with moderate diameter ratios of 1:2 to 1:4, this approach leads to discrepancies, and a correction factor is needed. In this work we demonstrate that this correction factor also holds for more extreme diameter ratios of 1:5, 1:7 and 1:10, although the force on the large particles is slightly overestimated when using the correction factor. The main origin of the correction is that the void surrounding the large particles becomes less in case of a bidisperse mixture, as compared to a monodisperse system with the same volume fraction. We further investigated this discrepancy by calculating the volume per particle by means of Voronoi tessellation.     

The effect of two-way coupling and inter-particle collisions on turbulence modulation in a vertical channel flow

Turbulence modulation due to its interaction with dispersed solid particles in a downward fully developed channel flow was studied. The Eulerian framework was used for the gas-phase, whereas the Lagrangian approach was used for the particle-phase. The steady-state equations of conservation of mass and momentum were used for the gas-phase, and the effect of turbulence on the flow-field was included via the standard k–ε model. The particle equation of motion included the drag, the Saffman lift and the gravity forces. Turbulence dispersion effect on the particles was simulated as a continuous Gaussian random field. The effects of particles on the flow were modeled by appropriate source terms in the momentum, k and ε equations. Particle–particle collisions and particle–wall collisions were accounted for in these simulations. Gas-phase velocities and turbulence kinetic energy in the presence of 2–100% mass loadings of two particle classes (50 μm glass and 70 μm copper) were evaluated, and the results were compared with the available experimental data and earlier numerical results. The simulation results showed that when the inter-particle collisions were important and was included in the computational model, the fluid turbulence was attenuated. The level of turbulence attenuation increased with particle mass loading, particle Stokes number, and the distance from the wall. When the inter-particle collisions were negligible and/or was neglected in the model, the fluid turbulence was augmented for the range of particle sizes considered.     

Influence of the viscous boundary layer on the deposition of particles from a gas suspension flowing past a sphere

The flow of incompressible gas containing particles past bodies of simple shapes at moderate and high velocities is investigated in [1–5], in which the flow of the carrier medium is assumed to be irrotational. The estimates made in [3] for the neighborhood of the stagnation point show that it is necessary to take into account the viscous boundary layer in the case of fine particles. In the present paper, the viscous flow of a gas suspension over the front surface of a sphere at Reynolds numbers R = 10³–10⁷ is considered. It is assumed that the carrier gas is incompressible and the particle concentra ion negligibly small. The influence of the boundary layer on the particle trajectories and the deposition of the disperse phase on the surface of the sphere is investigated. It is shown that there is a wide range of flow parameters for the gas suspension in which the influence of the boundary layer is important. The limits of this range are established.Translated from Izvestiya Akademli Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–66, January–February, 1982.I thank Yu. P. Savel'ev for a helpful discussion of the work.