Research on stability of moving jet containing dense suspended solid particles

Ｎｏｍｅｎｃｌａｔｕｒｅαｊ:ｐｈａｓｅｊｖｏｌｕｍｅｆｒａｃｔｉｏｎ;ｕｊ:ｐｈａｓｅｊｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎａｘｉａｌｖｅｌｏｃｉｔｙ;Ｕｏ:ｍｅａｎｘ＿ａｘｉａｌｖｅｌｏｃｉｔｙｏｆｔｈｅｅｘｉｔ;Ｕ∞ :ｖｅｌｏｃｉｔｙｏｆｊｅｔｔｉｎｇｄｅｖｉｃｅ;ΔＵ :ｊｕｍｐｉｎｇｖｅｌｏｃｉｔｙｏｆｍｏｖｉｎｇｊｅｔ,ｃｈａｒａｃｔｅｒｉｓｔｉｃｖｅｌｏｃｉｔｙｏｆｆｌｏｗ＿ｆｉｅｌｄ ,ΔＵ =Ｕｏ-Ｕ∞ ;ｐ:ｐｒｅｓｓｕｒｅ;Ｐ :ｅｉｇｅｎｐｒｅｓ…     

Quantification of dispersed phase concentration using light sheet imaging methods

With the prevalence of particle image velocimetry (PIV) as a quantitative tool for fluid mechanics diagnostics, its application for analyzing complicated multiphase flows has been steadily increasing over the last several decades. While the primary issue in using PIV for multiphase flows is in separating the information of the phases for independent analysis with a minimum of spurious “cross-talk,” an equally crucial but often overlooked point is in the accurate quantitative measurement of the dispersed phase concentration. Accurate concentration measurement is important due to the fact that the dispersed phase is often heterogeneously distributed in both space and time, either due to a non-uniformity of the source of particulates (such as a spray nozzle or sediment boundary) or due to inertial migration of the particles even from originally homogeneous spatial distributions. In the current work, we examine the effects of light sheet profile distortion and attenuation by tracer seeding particles, as well as reflected light from local wall boundaries on the effective light sheet thickness. The effective thickness is critical for concentration measurements, as it dictates the dispersed phase detection volume. A direct calibration method is demonstrated to measure the effective light sheet thickness in a water/glass bead system, which shows that systematic bias errors on the order of 30% can result if the reflective bed condition is not accounted for, and the errors can be as high as 50% or more if a single-point measure of the sheet width is used.     

Motion of a cloud of particles under gravity and its sedimentation on a flat horizontal surface

On the basis of the nonstationary two-dimensional equations of the mechanics of heterogeneous media a numerical investigation is made of the motion of a cloud of particles under the influence of gravity in an unbounded medium and the interaction of the descending cloud of particles with a flat horizontal surface. Depending on the degree of hydrodynamic interaction between the particles, there are different regimes of motion of the cloud during sedimentation; the change in its spatial configuration determined by the large-scale vortex motion of the carrier medium is determined. The concentration distributions of the particles on the sedimentation surface are obtained. A coefficient of dispersal is introduced for the cloud of particles on the flat horizontal surface, and its dependence on the concentration of the particles, their diameters, and other parameters is investigated.     

On hindered settling of particles of different sizes

The flow of a non-dilute fluid suspension is considered in which the dispersed phase consists of particles or droplets of different sizes. A phenomenological two-phase flow theory is formulated for both continuous and discrete distributions of particle sizes and illustrated by considering the batch settling of such a mixture. The volume fractions and particle distribution functions are determined, as well as the composition of the sedimentary layer.     

Effect of Interphase Heat Transfer on the Stability of a Turbulent Multicomponent Flow in a Vortex

The parameters of an axisymmetric turbulent two-phase swirling flow of a viscous heat-conducting gas containing a liquid dispersed phase in the presence of water vapor condensation on the particles are calculated. For the dispersed phase, a model taking into account the variation of the vapor concentration and the particle size due to condensation or evaporation is proposed. The distributions of the parameters of the basic unperturbed flow obtained numerically are used in the numerical solution of the linear problem of hydrodynamic stability within the time-dependent formulation. The parameters of small-amplitude harmonic perturbations propagating along the vortex axis are investigated in the linear formulation. A significant effect of heat release in the gas due to water vapor condensation on the parameters of the neutral perturbations and the neutral-stability curves is detected.     

Rotation of dispersed particles in a vortex field

An expression is obtained for the angular velocity of a spherical dispersed particle in a viscous fluid in an external vortex field with an harmonic time dependence. This expression is then used for investigating a system of two rotating dispersed particles whose rotation is the result of the interaction of the particles in the field of an incident sound wave. It is found that such a system possesses a rather interesting nontrivial property: under certain conditions it has a resonant frequency at which the rotation of the particles relative to the fluid is most intense.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 186–188, July–August, 1992.     

Interaction forces of particles dispersed in an electrolyte

Expressions are derived for the forces acting in a disperse medium in the presence of interaction of the double layers surrounding particles or drops of the dispersed phase when the potential of the dispersed particles is small. It is found that the force produced by the presence of double layers is proportional to the concentration gradient of the dispersed particles. It is shown that this force is comparable with the force produced by Brownian motion of the particles and may even exceed it. The equations of motion for the dispersed phase are derived with allowance for the convective terms, the pressure gradient, and the forces caused by Brownian motion and the presence of the double layers. A generalized Fick's law is obtained with effective diffusion coefficient. The equilibrium distribution of the particle concentration in a uniformly rotating cylinder is found with allowance for the interaction of the double layers.Translated from Izvestiya Akademii Nauk SSSR, Meklianika Zhidkosti i Gaza, No. 5, pp. 98–102, September–October, 1984.     

Thermocapillary convection in a liquid layer with a quasi-point heat source in the substrate

By the method of tracer particles the velocity field of thermocapillary convection in a thin layer of silicone oil, excited by a quasi-point heat source in the rigid substrate, is investigated as a function of the layer thickness, the temperature of the heater, and the liquid viscosity. The vertical velocity distributions are plotted in several cross-sections at different distances from the vortex axis. A novel method of measuring the profile of the thermocapillary depression, based on mirror reflection from the free liquid surface of radiation scattered by a tracer particle, is proposed. The central segment of the profile of the thermocapillary depression is obtained for different values of the layer thickness, the liquid viscosity and the heater temperature.     

Representation of classifier distributions in terms of hypergeometric functions

This paper derives alternative analytical expressions for classifier product distributions in terms of Gauss hypergeometric function, 2F1, by considering feed distribution defined in terms of Gates-Gaudin-Schumann function and efficiency curve defined in terms of a logistic function. It is shown that classifier distributions under dispersed conditions of classification pivot at a common size and the distributions are difference similar. The paper also addresses an inverse problem of classifier distributions wherein the feed distribution and efficiency curve are identified from the measured product distributions without needing to know the solid flow split of particles to any of the product streams.     

Direct numerical simulation of modulation of isotropic turbulence by poly‐dispersed particles

A direct numerical simulation technique based on two‐way coupling is presented to study a particle‐laden, decaying isotropic turbulent flow. Physical characteristics of turbulence modulation because of the mono‐dispersed (i.e., particles with single Stokes number) and poly‐dispersed particles (i.e., particles with more than one Stokes number) were investigated. A scale dependent effective viscosity that summarizes the aspects of the interaction between the velocity field and particles is defined in the study. Particles of Stokes number (St) 3.2,6.4 and 12.8 were used in performing the simulations. Poly‐dispersed particles were acquired by mixing particles of two different Stokes numbers at a time. As a whole, decay of turbulence because of the poly‐dispersed particles is observed to be larger than that of the decay of turbulence because of the mono‐dispersed particles. Simulations of poly‐dispersed particle indicate nonlinear characteristics in the modification of the temporal evolution of turbulence energy and dissipation. The scale dependent effective viscosity, which correlates with the energy spectrum plot, indicates that the decay of turbulence is mostly observed at the intermediate scales of turbulence. The effective viscosity for the simulations of the poly‐dispersed particles was calculated to be higher than that of the simulations of the mono‐dispersed particles. Copyright © 2011 John Wiley & Sons, Ltd.     

水击驻波场中分散相颗粒的运动分析

根据水击在管段内形成的驻波场现象,分析了流体内分散相颗粒受到的驻波作用力;运用李雅普诺夫稳定判据研究了颗粒积聚与分离的机理;考虑到颗粒运动方程的严重刚性而很难进行数值求解,采用相空间和非对称分析方法获得了分散相颗粒的运动轨迹近似解,并进行了实验验证。结果表明:水击驻波场中分散相颗粒的受力方程中惯性项对颗粒初始运动速率的影响不可忽略;在水击驻波波节的±λ/4范围内,分散相颗粒经过一定的时间会发生积聚,其运动速度呈对称分布,最大速度出现在3λ/8位置处;随着分散相颗粒粒径和密度等物性参数以及水击驻波的频率和连续相初始速度的增大,颗粒达到平衡位置的时间呈减小趋势,且连续相的初始速度对颗粒到达波节时间的影响显著。     

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.     

Statistical analysis of the hydrodynamic pressure in the near field of compressible jets

This paper is devoted to the statistical characterization of the pressure fluctuations measured in the near field of a compressible jet at two subsonic Mach numbers, 0.6 and 0.9. The analysis is focused on the hydrodynamic pressure measured at different distances from the jet exit and analyzed at the typical frequency associated to the Kelvin–Helmholtz instability. Statistical properties are retrieved by the application of the wavelet transform to the experimental data and the computation of the wavelet scalogram around that frequency. This procedure highlights traces of events that appear intermittently in time and have variable strength. A wavelet-based event tracking procedure has been applied providing a statistical characterization of the time delay between successive events and of their energy level. On this basis, two stochastic models are proposed and validated against the experimental data in the different flow conditions     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

Effect of viscoelasticity on the rotation of a sphere in shear flow

When particles are dispersed in viscoelastic rather than Newtonian media, the hydrodynamics will be changed entailing differences in suspension rheology. The disturbance velocity profiles and stress distributions around the particle will depend on the viscoelastic material functions. Even in inertialess flows, changes in particle rotation and migration will occur. The problem of the rotation of a single spherical particle in simple shear flow in viscoelastic fluids was recently studied to understand the effects of changes in the rheological properties with both numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346] and experiments [Snijkers et al., J. Rheol. 53 (2009) 459–480]. In the simulations, different constitutive models were used to demonstrate the effects of different rheological behavior. In the experiments, fluids with different constitutive properties were chosen. In both studies a slowing down of the rotation speed of the particles was found, when compared to the Newtonian case, as elasticity increases. Surprisingly, the extent of the slowing down of the rotation rate did not depend strongly on the details of the fluid rheology, but primarily on the Weissenberg number defined as the ratio between the first normal stress difference and the shear stress.In the present work, a quantitative comparison between the experimental measurements and novel simulation results is made by considering more realistic constitutive equations as compared to the model fluids used in previous numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346]. A multimode Giesekus model with Newtonian solvent as constitutive equation is fitted to the experimentally obtained linear and nonlinear fluid properties and used to simulate the rotation of a torque-free sphere in a range of Weissenberg numbers similar to those in the experiments. A good agreement between the experimental and numerical results is obtained. The local torque and pressure distributions on the particle surface calculated by simulations are shown.     

Probabilistic approach for transport of contaminants through porous media

The channels formed between individual particles in porous media have variable dimensions and orientations. The porosity, permeability and its anisotropy exhibit random spatial distributions. The probabilistic approach can effectively describe the transport of contaminants through porous media and is analysed in this paper. Numerical results are obtained by considering (I) random dispersion coefficients without and with spatial structure, (II) random time distribution of concentration at the inlet boundary, (III) random velocity distribution in the flow field without and (IV) with variable dispersion coefficient, (V) non-linearity of the governing equation and (VI) anisotropy of the dispersion coefficient. Two methods are used for probabilistic predictions: (1) Gaussian field approach in conjunction with Monte Carlo method and (2) random walk method. The input random parameters are assumed to have normal and log-normal distributions according to available experimental data. The probability distribution functions of the contaminant concentration at different locations within the flow domain are calculated and compared with the input distributions as a function of the mean and fluctuation Peclet numbers. The one-dimensional case is analysed in detail and the illustrative numerical predictions are compared with analytical and experimental results. The extension to a two-dimensional domain is discussed in the last part of this paper.     

Numerical simulation of oscillations of a monodisperse gas-particle mixture in a nonlinear wave field

This paper presents a numerical simulation procedure for the dynamics of a monodisperse gas-particle mixture in the nonlinear wave field of an acoustic resonator using a two-temperature two-velocity model ignoring phase transitions, particle collision, and possible coagulation. It is assumed that viscosity is present only in the carrier medium described by the Navier-Stokes equations for a compressible gas. The dispersed phase is described by the equation of conservation of mass, momentum, and energy. A monotonic solution is obtained by solving the equations of motion for the carrier medium and dispersed phase in generalized moving coordinates using the explicit McCormack method with splitting in the spatial directions and a conservative correction scheme. The method can be used to study nonlinear oscillations of two-phase mixtures in the vicinity of the first three eigenfrequencies in a flat channel.