Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Unsteady turbulent flow of a gas suspension in a channel under conditions of injection and forced pressure oscillations

A turbulent flow of a suspension of solid particles in a gas is considered. The suspension is located in a channel with permeable walls (the pressure at the left end face of the channel follows a sinusoidal law). The flow considered here reflects the principal features of the flow in the combustion chamber of a solid-propellant rocket motor. The unsteady flow of the gas suspension is described by using the Eulerian-Lagrangian approach. A stochastic variant of the discrete-trajectory approach is used for modeling the particle motion. The influence of the condensed phase on the turbulence characteristics and acoustic oscillations of the parameters of the working medium in the channel in the case of injection is discussed. The calculated results are compared with data obtained in a physical experiment.     

Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

Simulation on motion of particles in vortex merging process

In a two-phase flow, the vortex merging influences both the flow evolution and the particle motion. With the blobs-splitting-and-merging scheme, the vortex merging is calculated by a corrected core spreading vortex method （CCSVM）. The particle motion in the vortex merging process is calculated according to the particle kinetic model. The results indicate that the particle traces are spiral lines with the same rotation direction as the spinning vortex. The center of the particle group is in agreement with that of the merged vortex. The merging time is determined by the circulation and the initial ratio of the vortex radius and the vortex center distance. Under a certain initial condition, a stretched particle trail is generated, which is determined by the viscosity, the relative position between the particles and the vortex, and the asymmetrical circulation of the two merging vortices.     

Motion and heat and mass transfer in a disperse admixture in turbulent nonisothermal jets of a gas and a low-temperature plasma

The motion and heat and mass transfer of particles of a disperse admixture in nonisothermal jets of a gas and a low-temperature plasma are simulated with allowance for the migration mechanism of particle motion actuated by the turbophoresis force and the influence of turbulent fluctuations of the jet flow velocity on heat and mass transfer of the particle. The temperature distribution inside the particle at each time step is found by solving the equation of unsteady heat conduction. The laws of scattering of the admixture and the laws of melting and evaporation of an individual particle are studied, depending on the injection velocity and on the method of particle insertion into the jet flow. The calculated results are compared with data obtained with ignored influence of turbulent fluctuations on the motion and heat and mass transfer of the disperse phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 95–108, May–June, 2008.     

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.     

Calculation of particle concentration in the problem of aerosol aspiration into a thin-walled tube

A mathematical model of the problem of aerosol aspiration into a thin-walled tube from a moving gas is developed and particle concentration fields are calculated. In the absence of a particle effect on the gas flow, the carrier medium is calculated in the potential-flow and viscous-flow approximations for an incompressible gas, using boundary-element and finite-volume methods. For the viscous-flow model, a numerical solution is found using the FLUENT program. The particle motion equations are complemented with equations for calculating the concentration along the particle trajectories. The spatial distributions of the particle concentration near and inside the sampler are studied for different ratios of the wind and aspiration velocities and for different Stokes numbers. The effect of nonuniformity of the particle concentration distributions on the aspiration coefficient is discussed.     

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.     

Atomization of a Powder by a Swirling Flow with a Recirculation Zone

The problem of the motion of a swirling flow in an axisymmetric channel with permeable walls is investigated numerically. Various flow regimes including those with the formation of recirculation zones are obtained. The problem of atomization of a powder by a swirling flow for the purpose of obtaining a finely dispersed mixture is considered. Particle concentration distributions in the flow are calculated, the formation of characteristic deposition zones is demonstrated, and the unsteady process of particle transfer is investigated with allowance for deposition on the lateral surface of the channel.     

Modeling of Ice Accretion on the Airfoil Surface in an Air Flow Containing Ice Particles

The process of airfoil icing caused by incidence of ice crystals is considered. A physicomathematical model of motion of spheroidal crystals in the gas flow and their interaction with the body is formulated. The model takes into account the non-spherical particle orientation with respect to the velocity vector of the gas flow. It is assumed that the particle impact onto the body surface leads to partial destruction of the particle under the action of normal and tangential stresses, and some part of the particle mass remains in the vicinity of the impact point. An inviscid flow around the airfoil with a time-dependent shape is calculated by the method of approximate conformal mapping.     

Formation of particle jetting in a cylindrical shock tube

A dense, solid particle flow is numerically studied at a mesoscale level for a cylindrical shock tube problem. The shock tube consists of a central high pressure gas driver section and an annular solid powder bed with air in void regions as a driven section with its far end adjacent to ambient air. Simulations are conducted to explore the fundamental phenomena, causing clustering of particles and formation of coherent particle jet structures in such a dense solid flow. The influence of a range of parameters is investigated, including driver pressure, particle morphology, particle distribution and powder bed configuration. The results indicate that the physical mechanism responsible for this phenomenon is twofold: the driver gas jet flow induced by the shock wave as it passes through the initial gaps between the particles in the innermost layer of the powder bed, and the chaining of solid particles by inelastic collision. The particle jet forming time is determined as the time when the motion of the outermost particle layer of the powder bed is first detected. The maximum number of particle jets is bounded by the total number of particles in the innermost layer of the powder bed. The number of particle jets is mainly a function of the number of particles in the innermost layer and the mass ratio of the powder bed to the gas in the driver section, or the ratio of powder bed mass (in dimensionless form) to the pressure ratio between the driver and driven sections.     

Numerical study of the dynamics of a vapor envelope round a heated solid particle in a liquid

A study is made of the radial motion of a vapor envelope surrounding an isolated spherical particle in an unbounded mass of liquid. It is assumed that the liquid is viscous and incompressible and that the temperature is distributed uniformly in the solid particle. A model of a calorifically perfect gas is used for the vapor phase. The same assumptions are made as in Rayleigh's formulation for the problem of the dynamics of a single bubble: that the process is spherically symmetric and that the pressure P₂ (t) in the vapor phase is homogeneous. The justification for making these assumptions in problems of the dynamics of gas, vapor, and vaporgas bubbles is discussed in [1–5]. In this paper, the collapse of the vapor layer and the boiling of the liquid on the surface of the heated particle are not considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 154–157, July–August, 1985.     

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4 :2244–2251; Int. J. Multiphase Flow 2000; 26 :1583–1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamics of particles floating on liquid flowing in a semicircular open channel

The dynamic characteristics of surface-floating particles in liquids flowing in a two-dimensional, semicircular open channel is studied experimentally. For high visibility in the experiments, relatively large particles are employed whose particle-liquid density ratio is either equal to or less than unity. Particles of different size and geometry are tested in a water-glycerin mixture. A video camera traces the pathline of each particle from which the velocity and direction of particle motion are evaluated. Liquid velocity distribution is determined by hot-film anemometry. A modified dynamics (Basset-Boussinesq-Oseen) equation is derived and numerically solved by means of a finite-difference technique to determine fluid velocity. A new dimensionless parameter is disclosed which is pertinent to both particle geometry and fluid flow conditions. It correlates particle trajectory and velocity, trajectory dispersion and fluid-particle velocity ratio.Visiting Scholar on leave from Department of Mechanical Engineering, Fukuyama University, Fukujama, Japan     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

An analysis of the chaotic motion of particles of different sizes in a gas fluidized bed

The dynamic behavior of individual particles during the mixing/segregation process of particle mixtures in a gas fluidized bed is analyzed. The analysis is based on the results generated from discrete particle simulation, with the focus on the trajectory of and forces acting on individual particles.Typical particles are selected representing three kinds of particle motion:a flotsam particle which is initially at the bottom part of the bed and finally fluidized at the top part of the bed; a jetsam particle which is initially at the top part of the bed and finally stays in the bottom de-fluidized layer of the bed; and a jetsam particle which is intermittently joining the top fluidized and bottom de-fluidized layers. The results show that the motion of a particle is chaotic at macroscopic or global scale, but can be well explained at a microscopic scale in terms of its interaction forces and contact conditions with other particles, particle-fluid interaction force, and local flow structure. They also highlight the need for establishing a suitable method to link the information generated and modeled at different time and length scales.     

Calculation of the flow of a two-phase mixture in a Laval nozzle,allowing for the turbulent diffusion of the particles

The flow of a mixture of gas and condensed particles in an axisymmetric Laval nozzle is considered. The motion of the particles is calculated in a specified field of gas flow, with due allowance for their turbulent diffusion. The results of calculations indicating the necessity of allowing for this phenomenon when considering the motion of particles toward the wall of a profiled nozzle are presented.Translated from Izvestiya Akademii Nauk SSSR, No. 2, pp. 161–165, March–April, 1973.     

Cold flow experiments in an entrained flow gasification reactor with a swirl-stabilized pulverized biofuel burner

Short particle residence time in entrained flow gasifiers demands the use of pulverized fuel particles to promote mass and heat transfer, resulting high fuel conversion rate. The pulverized biomass particles have a wide range of aspect ratios which can exhibit different dispersion behavior than that of spherical particles in hot product gas flows. This results in spatial and temporal variations in temperature distribution, the composition and the concentration of syngas and soot yield. One way to control the particle dispersion is to impart a swirling motion to the carrier gas phase. This paper investigates the dispersion behavior of biomass fuel particles in swirling flows. A two-phase particle image velocimetry technique was applied to simultaneously measure particle and gas phase velocities in turbulent isothermal flows. Post-processed PIV images showed that a poly-dispersed behavior of biomass particles with a range of particle size of 112–160 µm imposed a significant impact on the air flow pattern, causing air flow decelerated in a region of high particle concentration. Moreover, the velocity field, obtained from individually tracked biomass particles showed that the swirling motion of the carrier air flow gives arise a rapid spreading of the particles.     

Relative particle velocity in a turbulent flow

On the basis of a statistical approach using a probability density function for the coordinates of two particles in a turbulent flow, the parameters of the relative particle motion are investigated. For the functions describing particle entrainment in the turbulence, rigorous results are obtained using a 3D turbulence spectrum. A method of calculating the particle relative-velocity rate with account for particle trajectory correlation is presented. The effects of particle inertia and velocity slip on the parameters of the relative particle motion are studied. Simple approximating formulas for calculating the relative particle motion in a turbulent flow are proposed. The calculation results are compared with the data of direct numerical simulation of stochastic particle trajectories in an isotropic turbulent field.     

Migration of particles in a turbulent nonisothermic flow as a result of dependence of the viscosity of the gas on the temperature

The motion is considered of a Stokes-like spherical particle in a turbulent nonisothermic gas flow whose viscosity depends on the temperature. The field of the turbulent velocity is assumed to be homogeneous, isotropic, and steady. It is shown that if there is a mean temperature gradient in the gas, and, consequently, a heat flow due to turbulent pulsations, then there may be turbulent migration of particles in a direction collinear with the gradient of the mean temperature. The migration is due to statistical correlation of turbulent pulsations of velocity and temperature, and is not connected with the phenomenon of ordinary thermophoresis. Upon the introduction of a number of simplifying assumptions, the rate of migration is calculated in dependence on the characteristics of the particle and the flow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 53–58, November–December, 1986. The author is grateful to V. S. Galkin, V. A. Zharov, M. N. Kogan, and V. A. Sabel'nikov for discussions of the study.