Preferential concentration driven instability of sheared gas–solid suspensions

2015 ◽  
Vol 770 ◽  
pp. 85-123 ◽  
Author(s):  
M. Houssem Kasbaoui ◽  
Donald L. Koch ◽  
Ganesh Subramanian ◽  
Olivier Desjardins

We examine the linear stability of a homogeneous gas–solid suspension of small Stokes number particles, with a moderate mass loading, subject to a simple shear flow. The modulation of the gravitational force exerted on the suspension, due to preferential concentration of particles in regions of low vorticity, in response to an imposed velocity perturbation, can lead to an algebraic instability. Since the fastest growing modes have wavelengths small compared with the characteristic length scale ($U_{g}/{\it\Gamma}$) and oscillate with frequencies large compared with ${\it\Gamma}$, $U_{g}$ being the settling velocity and ${\it\Gamma}$ the shear rate, we apply the WKB method, a multiple scale technique. This analysis reveals the existence of a number density mode which travels due to the settling of the particles and a momentum mode which travels due to the cross-streamline momentum transport caused by settling. These modes are coupled at a turning point which occurs when the wavevector is nearly horizontal and the most amplified perturbations are those in which a momentum wave upstream of the turning point creates a downstream number density wave. The particle number density perturbations reach a finite, but large amplitude that persists after the wave becomes aligned with the velocity gradient. The growth of the amplitude of particle concentration and fluid velocity disturbances is characterised as a function of the wavenumber and Reynolds number ($\mathit{Re}=U_{g}^{2}/{\it\Gamma}{\it\nu}$) using both asymptotic theory and a numerical solution of the linearised equations.

Author(s):  
M. Sato ◽  
M. Tanahashi ◽  
T. Miyauchi

Direct numerical simulations of homogeneous isotropic turbulence laden with particles have been conducted to clarify the relationship between particle dispersion and coherent fine scale eddies in turbulence. Dispersion of 106 particles are analyzed for several particle Stokes numbers. The spatial distributions of particles depend on their Stokes number, and the Stokes number that causes preferential concentration of particles is closely related to the time scale of coherent fine scale eddies in turbulence. On the plane perpendicular to the rotating axes of fine scale eddies, number density of particle with particular Stokes number is low at the center of the fine scale eddy, and high in the regions with high energy dissipation rate around the eddy. The maximum number density can be observed at about 1.5 to 2.0 times the eddy radius on the major axis of the fine scale eddy.


2011 ◽  
Vol 29 (2) ◽  
pp. 201-217 ◽  
Author(s):  
Kaushik Balakrishnan ◽  
Suresh Menon

AbstractA new multiphase buoyancy-drag model is developed for the study of Rayleigh-Taylor and Richtmyer-Meshkov instabilities in dusty gases, extending on a counterpart single-phase model developed in the past by Srebro et al. (2003). This model is applied to single- and multi-mode perturbations in dusty gases and both Rayleigh-Taylor and Richtmyer-Meshkov instabilities are investigated. The amplitude for Rayleigh-Taylor growth is observed to be contained within a band, which lies within limits identified by a multiphase Atwood number that is a function of the fluid densities, particle size, and a Stokes number. The amplitude growth is subdued with (1) an increase in particle size for a fixed particle number density and with (2) an increase in particle number density for a fixed particle size. The power law index for Richtmyer-Meshkov growth under multi-mode conditions also shows dependence to the multiphase Atwood number, with the index for the bubble front linearly decreasing and then remaining constant, and increasing non-linearly for the spike front. Four new classes of problems are identified and are investigated for Rayleigh-Taylor growth under multi-mode conditions for a hybrid version of the model: (1) bubbles in a pure gas rising into a region of particles; (2) spikes in a pure gas falling into a region of particles; (3) bubbles in a region of particles rising into a pure gas; and (4) spikes in a region of particles falling into a pure gas. Whereas the bubbles accelerate for class (1) and the spikes for class (4), for classes (2) and (3), the spikes and bubbles, respectively, oscillate in a gravity wave-like phenomenon due to the buoyancy term changing sign alternatively. The spike or bubble front, as the case may be, penetrates different amounts into the dusty or pure gas for every subsequent penetration, due to drag effects. Finally, some extensions to the presently developed multiphase buoyancy-drag model are proposed for future research.


2017 ◽  
Vol 9 (1) ◽  
pp. 32-46 ◽  
Author(s):  
Javier Achury ◽  
Wolfgang Polifke

Multiple interactions may occur when a poly-disperse spray is exposed to an acoustic field. In the context of spray combustion instabilities, acoustic agglomeration, the formation of a droplet number density wave and the modulation of the droplet size distribution are interesting effects. A droplet number density wave, i.e. preferential concentration of droplets in space, may result from size-dependent, one-way momentum coupling between the acoustic field and the spray. The modulation of the droplet size distribution, which has been evidenced in the experimental work of Gurubaran and Sujith (AIAA 2008-1046), is thus a consequence of the droplet number density wave formation. In the present work, the mechanisms that produce these effects are simulated and analyzed in depth by means of computational fluid dynamics. The spray is modeled with both Lagrangian (particles mass-point approach) and Eulerian (continuous phase approach) descriptions. The particular Eulerian method used is a variant of the presumed density function method of moments, which allows to account for the effects of poly-dispersity, in particular the size-dependence of particle velocity. Both the Lagrangian and Eulerian models are validated against experimental data for spray dynamics and spray response to an acoustic field.


1993 ◽  
Vol 7 (6) ◽  
pp. 842-851 ◽  
Author(s):  
M. Queiroz ◽  
M. P. Bonin ◽  
J. S. Shirolkar ◽  
R. W. Dawson

2016 ◽  
Vol 7 (32) ◽  
pp. 5123-5131 ◽  
Author(s):  
O. L. J. Virtanen ◽  
M. Brugnoni ◽  
M. Kather ◽  
A. Pich ◽  
W. Richtering

Many applications of poly(N-isopropylacrylamide) microgels necessitate robust control over particle size.


2017 ◽  
Vol 34 (7) ◽  
pp. 075203
Author(s):  
Rang-Yue Zhang ◽  
Yan-Hong Liu ◽  
Feng Huang ◽  
Zhao-Yang Chen ◽  
Chun-Yan Li

2017 ◽  
Vol 822 ◽  
pp. 640-663 ◽  
Author(s):  
J. L. G. Oliveira ◽  
C. W. M. van der Geld ◽  
J. G. M. Kuerten

Three-dimensional particle tracking velocimetry is applied to particle-laden turbulent pipe flows at a Reynolds number of 10 300, based on the bulk velocity and the pipe diameter, for developed fluid flow and not fully developed flow of inertial particles, which favours assessment of the radial migration of the inertial particles. Inertial particles with Stokes number ranging from 0.35 to 1.11, based on the particle relaxation time and the radial-dependent Kolmogorov time scale, and a ratio of the root-mean-square fluid velocity to the terminal velocity of order 1 have been used. Core peaking of the concentration of inertial particles in up-flow and wall peaking in down-flow have been found. The difference in mean particle and Eulerian mean liquid velocity is found to decrease to approximately zero near the wall in both flow directions. Although the carrier fluid has all of the characteristics of the corresponding turbulent single-phase flow, the Reynolds stress of the inertial particles is different near the wall in up-flow. These findings are explained from the preferential location of the inertial particles with the aid of direct numerical simulations with the point-particle approach.


2001 ◽  
Author(s):  
S. L. Chang ◽  
C. Q. Zhou ◽  
B. Golchert ◽  
M. Petrick

Abstract A typical glass furnace consists of a combustion space and a melter. The intense heat, generated from the combustion of fuel and air/oxygen in the combustion space, is transferred mainly by radiation to the melter where the melt sand and cullet (scrap glass) are melted, creating molten glass. The melter flow is a complex multi-phase flow including solid batches of sand/cullet and molten glass. Proper modeling of the flow patterns of the solid batch and liquid glass is a key to determining the glass quality and furnace efficiency. A multi-phase CFD code has been developed to simulate glass melter flow. It uses an Eulerian approach for both the solid batch and the liquid glass-melt flows. The mass, momentum, and energy conservation equations of the batch flow are used to solve for local batch particle number density, velocity, and temperature. In a similar manner, the conservation equations of mass, momentum, and energy of the glass-melt flow are used to solve for local liquid molten glass pressure, velocity, and temperature. The solid batch is melted on the top by the heat from the combustion space and from below by heat from the glass-melt flow. The heat transfer rate from the combustion space is calculated from a radiation model calculation while the heat transfer rate from the glass-melt flow to the solid batch is calculated from a model based on local particle number density and glass-melt temperature. Energy and mass are balanced between the batch and the glass-melt. Batch coverage is determined from local particle number density and velocity. A commercial-scale glass melter has been simulated at different operating/design conditions.


1991 ◽  
Vol 130 ◽  
pp. 71-74
Author(s):  
A.Z. Dolginov ◽  
N.A. Silant’ev

AbstractA new method for the calculation of kinetic coefficients is presented. This method allows us to obtain the distribution of scalar and vector fields (such as the temperature, the admixture particle number density and the magnetic field) in turbulent cosmic media with any value of S = u0т0/R0. The explicit expression for the “turbulent” diffusivity DT is obtained. In some cases DT becomes negative, implying the clustering of the admixture particles in patches (a local increase of the temperature and magnetic fields). The magnetic α-effect is considered for the case S ~ 1.


2008 ◽  
Vol 4 (1) ◽  
pp. 69-78 ◽  
Author(s):  
D. B. Kirk-Davidoff ◽  
J.-F. Lamarque

Abstract. Previous work has shown that polar stratospheric clouds (PSCs) could have acted to substantially warm high latitude regions during past warm climates such as the Eocene (55 Ma). Using a simple model of stratospheric water vapor transport and polar stratospheric cloud (PSC) formation, we investigate the dependence of PSC optical depth on tropopause temperature, cloud microphysical parameters, stratospheric overturning, and tropospheric methane. We show that PSC radiative effects can help slow removal of water from the stratosphere via self-heating. However, we also show that the ability of PSCs to have a substantial impact on climate depends strongly on the PSC particle number density and the strength of the overturning circulation. Thus even a large source of stratospheric water vapor (e.g. from methane oxidation) will not result in substantial PSC radiative effects unless PSC ice crystal number density is high compared to most current observations, and stratospheric overturning (which modulates polar stratospheric temperatures) is low. These results are supported by analysis of a series of runs of the NCAR WACCM model with methane concentrations varying up to one thousand times present levels.


Sign in / Sign up

Export Citation Format

Share Document