Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

Accounting for finite-size effects in simulations of disperse particle-laden flows

A numerical formulation for Eulerian–Lagrangian simulations of particle-laden flows in complex geometries is developed. The formulation accounts for the finite-size of the dispersed phase. Similar to the commonly used point-particle formulation, the dispersed particles are treated as point-sources, and the forces acting on the particles are modeled through drag and lift correlations. In addition to the inter-phase momentum exchange, the presence of particles affects the fluid phase continuity and momentum equations through the displaced fluid volume. Three flow configurations are considered in order to study the effect of finite particle size on the overall flowfield: (a) gravitational settling, (b) fluidization by a gaseous jet, and (c) fluidization by lift in a channel. The finite-size formulation is compared to point-particle representations, which do not account for the effect of finite-size. It is shown that the fluid displaced by the particles plays an important role in predicting the correct behavior of particle motion. The results suggest that the standard point-particle approach should be modified to account for finite particle size, in simulations of particle-laden flows.     

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.     

Mean and fluctuating components of drag and lift forces on an isolated finite-sized particle in turbulence

We perform fully resolved direct numerical simulations of an isolated particle subjected to free-stream turbulence in order to investigate the effect of turbulence on the drag and lift forces at the level of a single particle, following Bagchi and Balachandar’s work (Bagchi and Balachandar in Phys Fluids 15:3496–3513, 2003). The particle Reynolds numbers based on the mean relative particle velocity and the particle diameter are Re?=?100, 250 and 350, which covers three different regimes of wake evolution in a uniform flow: steady axisymmetric wake, steady planar symmetric wake, and unsteady planar symmetric vortex shedding. At each particle Reynolds number, the turbulent intensity is 5–10% of the mean relative particle velocity, and the corresponding diameter of the particle is comparable to or larger than the Kolmogorov scale. The simulation results show that standard drag values determined from uniform flow simulations can accurately predict the drag force if the turbulence intensity is sufficiently weak (5% or less compared to the mean relative velocity). However, it is shown that for finite-sized particles, flow non-uniformity, which is usually neglected in the case of the small particles, can play an important role in determining the forces as the relative turbulence intensity becomes large. The influence of flow non-uniformity on drag force could be qualitatively similar to the Faxen correction. In addition, finite-sized particles at sufficient Reynolds number are inherently subjected to stochastic forces arising from their self-induced vortex shedding in addition to lift force arising from the local ambient flow properties (vorticity and strain rate). The effect of rotational and strain rate of the ambient turbulence seen by the particle on the lift force is explored based on the conditional averaging using the generalized representation of the quasi-steady force proposed by Bagchi and Balachandar (J Fluid Mech 481:105–148, 2003). From the present study, it is shown that at Re?=?100, the lift force is mainly influenced by the surrounding turbulence, but at Re = 250 and 350, the lift force is affected by the wake structure as well as the surrounding turbulence. Thus, for a finite-sized particle of sufficient Reynolds number supporting self-induced vortex shedding, the lift force will not be completely correlated with the ambient flow. Therefore, it appears that in order to reliably predict the motion of a finite-sized particle in turbulence, it is important to incorporate both a deterministic component and a stochastic component in the force model. The best deterministic contribution is given by the conditional average. The influence of ambient turbulence at the scale of the particle, which are not accounted for in the deterministic contribution, can be considered in stochastic manner. In the modeling of lift force, additional stochastic contribution arising from self-induced vortex shedding must also be included.     

关于颗粒悬浮机理和悬浮动的讨论

从分析气体分子的悬浮和静水中Ｂｒｏｗｎ微粒的悬浮之机理出发，论述了重力场中粒子（分子、微粒等）的悬浮不一定需要其它外力，粒子本身的任何形式的无规则运动，达到一定强度后都能使粒子弥散悬浮．河流中的泥沙颗粒和气（水）力输送管道中的颗粒的悬浮也主要靠颗粒物的无规则运动．作用于颗粒的升力和其它力可改变颗粒悬浮沿高度的分布，但仅用这些力（若无任何无规则运动）无法解释颗粒的弥散悬浮状态．讨论了颗粒对流动阻力的双重作用：支持颗粒悬浮的湍流脉动因引入颗粒而削弱，这是颗粒的减阻作用；颗粒增阻的一个主要机制是，流体给予颗粒的水平动量在颗粒一壁面碰撞中不断地损失．用悬浮动概念解释颗粒引起的增阻是不正确的．     

Numerical modelling of the entrainment of particles in inviscid supersonic flow

The interaction between particles situated in close proximity and moving at supersonic speeds is investigated computationally. The simplest case of the motion of a single particle travelling behind a lead particle is used to elucidate the role of aerodynamic forces in the motion of a group of particles. The effect of the following parameters on the drag and lift forces acting on each of two particles of equal diameter in proximity is investigated: the free-stream Mach number, and the axial and lateral displacements of the trailing particle. The two-dimensional flow field is numerically simulated using an unsteady Euler CFD code to find the steady-state drag and lift coefficients for both particles. Three static zones of aerodynamic influence in the wake of the lead particle are identified, which are denoted as the entrainment, lateral attraction, and ejection zones. A non-dimensional representation of the zones of influence is given. It is shown that the dynamic entrainment of particles can occur even when the path of the trailing particle originates outside the entrainment and lateral attraction zones.     

On the interaction between two fixed spherical particles

The variation of the drag (C_D) and lift coefficients (C_L) of two fixed solid spherical particles placed at different positions relative each other is studied. Simulations are carried out for particle Reynolds numbers of 50, 100 and 200 and the particle position is defined by the angle between the line connecting the centers of the particles and the free-stream direction (α) and the separation distance (d₀) between the particles. The flow around the particles is simulated using two different methods; the Lattice Boltzmann Method (LBM), using two different computational codes, and a conventional finite difference approach, where the Volume of Solid Method (VOS) is used to represent the particles. Comparisons with available numerical and experimental data show that both methods can be used to accurately resolve the flow field around particles and calculate the forces the particles are subjected to. Independent of the Reynolds number, the largest change in drag, as compared to the single particle case, occurs for particles placed in tandem formation. Compared to a single particle, the drag reduction for the secondary particle in tandem arrangement is as high as 60%, 70% and 80% for Re = 50, 100 and 200, respectively. The development of the recirculation zone is found to have a significant influence on the drag force. Depending on the flow situation in-between the particles for various particle arrangements, attraction and repulsion forces are detected due to low and high pressure regions, respectively. The results show that the inter-particle forces are not negligible even under very dilute conditions.     

Comparison of the RANS and PDF methods for air-particle flows

Two simulation methods, namely Reynolds-Averaged Navier–Stokes (RANS) equations, and Probability Distribution Function (PDF) are currently widely used for the modeling of multiphase flows. These two approaches are supplemented with appropriate closure equations that take into account all the pertinent forces and interaction effects on the solid particles, such as: particle–turbulence interactions; turbulence modulation; particle–particle interactions; particle–wall interactions; gravitation, drag and lift forces. The two methods have been used in order to simulate the turbulent particulate flow in upward pipes. The flow domain in all cases was a cylindrical pipe and the computations were carried for upward pipe flow. Monodisperse as well as polydisperse mixtures of particles have been considered. In general, the average velocity results obtained from the two methods are in close agreement, because the methods predict well the average velocity distribution of the carrier fluid as well as the solids. Thus, the differences in the average axial velocities predicted by the methods are not substantial. Differences in the turbulence intensity are more significant. A comparison of the numerical results obtained shows the relative importance of retaining the diffusion terms in both the axial and radial directions in the RANS method. Also the comparisons of the results show the relative effect of the lift forces in the distribution of solid particles.     

明渠层流近床面颗粒拖曳力研究

泥沙颗粒受到的拖曳力是泥沙运动的主要驱动力,而当前应用于计算流体力学-离散颗粒法(CFD-DPM)耦合模型进行水沙运动模拟的泥沙颗粒拖曳力公式均没有考虑明渠流底床边壁作用的影响。求解不可压缩Navier-Stokes方程,对明渠层流不同雷诺数条件下床面附近不同高度处颗粒所受拖曳力进行了模拟,根据模拟结果变化规律,提出了综合考虑床面和水流惯性对标准拖曳力影响的修正拖曳力计算公式。与常用的单颗粒标准拖曳力公式和考虑遮蔽效应的多颗粒拖曳力公式相比,采用本文修正公式得到的水沙作用力更接近高精度数值解,应用于CFD-DPM输沙模拟获得的输沙结果与输沙率公式结果一致,应用分析表明输沙模拟应当采用粗糙底床边界。     

Measurement of the aerodynamic forces on a small particle attached to a wall

Using operating principles similar to that applied in atomic force microscopes, we have developed a novel measuring method to study the aerodynamic forces, in particular the lift and drag force, acting on a small particle attached to a wall and immersed in a linear shear flow. Results thus far have shown that the system is capable of measuring both the minute aerodynamic lift and drag forces that a particle experiences as a result of the flow.C. Muthanna has also published under the name C. M. Kolera

---

History force on a sphere in a weak linear shear flow

A numerical study of history forces acting on a spherical particle in a linear shear flow, over a range of finite Re, is presented. In each of the cases considered, the particle undergoes rapid acceleration from Re₁ to Re₂ over a short-time period. After acceleration, the particle is maintained at Re₂ in order to allow for clean extraction of drag and lift kernels. Good agreement is observed between current drag kernel results and previous investigations. Furthermore, ambient shear is found to have little influence on the drag kernel. The lift kernel is observed to be oscillatory, which translates to a non-monotonic change in lift force to the final steady state. In addition, strong dependence on the start and end conditions of acceleration is observed. Unlike drag, the lift history kernel scales linearly with Reynolds number and shear rate. This behavior is consistent with a short-time inviscid evolution. A simple expression for the lift history kernel is presented.     

Effect of Shape on Inertial Particle Dynamics in a Channel Flow

Particle dynamics in a channel flow are investigated using large eddy simulation and a Lagrangian particle tracking technique. Following validation of single-phase flow predictions against DNS results, fluid velocities are subsequently used to study the behaviour of particles of differing shape assuming one-way coupling between the fluid and the particles. The influence of shape- and orientation-dependent drag and lift forces on both the translational and rotational motion of the particles is accounted for to ensure accurate representation of the flow dynamics of non-spherical particles. The size of the particles studied was obtained based on an equivalent-volume sphere, and differing shapes were modelled using super-quadratic ellipsoid forms by varying their aspect ratio, with their orientation predicted using the incidence angle between the particle relative velocity and the particle principal axis. Results are presented for spherical, needle- and platelet-like particles at a number of different boundary layer locations along the wall-normal direction within the channel. The time evolution and probability density function of selected particle translational and rotational properties show a clear distinction between the behaviour of the various particles types, and indicate the significance of particle shape when modelling many practically relevant flows.     

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.     

Drag of a flat wedge in a two-phase flow

The drag of a flat wedge in a subsonic two-phase flow is investigated. In contrast to earlier work of Balanin and Zlobin [1] particular attention is devoted to the influence of the particle size. Detailed investigations are made of the dependences of the forces and aerodynamic coefficients on the concentration of the solid phase, the opening angle of the wedge, and the particle size. It is established that the drag coefficients depend on the particle size only for particles with diameters less than 30 um.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 177–180, March–April, 1982.     

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.      

Drag forces of interacting spheres in power-law fluids

Drag forces of interacting particles suspended in power-law fluid flows were investigated in this study. The drag forces of interacting spheres were directly measured by using a micro-force measuring system. The tested particles include a pair of interacting spheres in tandem and individual spheres in a cubic matrix of multi-sphere in flows with the particle Reynolds number from 0.7 to 23. Aqueous carboxymethycellulose (CMC) solutions and glycerin solutions were used as the fluid media in which the interacting spheres were suspended. The range of power-law index varied from 0.6 to 1.0. In conjunction to the drag force measurements, the flow patterns and velocity fields of power-law flows over a pair of interacting spheres were also obtained from the laser assisted flow visualization and numerical simulation.

Both experimental and computational results suggest that, while the drag force of an isolated sphere depends on the power-index, the drag coefficient ratio of an interacting sphere is independent from the power-law index but strongly depends on the separation distance and the particle Reynolds number. Our study also shows that the drag force of a particle in an assemblage is strongly positions dependent, with a maximum difference up to 38%.     

Hydrodynamic and thermal fields analysis in gas-solid two-phase flow

The present work aims to investigate numerically the flowfield and heat transfer process in gas-solid suspension in a vertical pneumatic conveying pipe. The Eulerian-Lagrangian model is used to simulate the flow of the two-phases. The gas phase is simulated based on Reynolds Average Navier-Stokes equations (RANS) with low Reynolds number k-ε model, while particle tracking procedure is used for the solid phase. An anisotropic model is used to calculate the Reynolds stresses and the turbulent Prandtl number is calculated as a function of the turbulent viscosity. The model takes into account the lift and drag forces and the effect of particle rotation as well as the particles dispersion by turbulence effect. The effects of inter-particles collisions and turbulence modulation by the solid particles, i.e. four-way coupling, are also included in the model. Comparisons between different models for turbulence modulation with experimental data are carried out to select the best model. The model is validated against published experimental data for velocities of the two phases, turbulence intensity, solids concentration, pressure drop, heat transfer rates and Nusselt number distribution. The comparisons indicate that the present model is able to predict the complex interaction between the two phases in non-isothermal gas-solid flow in the tested range. The results indicate that the particle-particle collision, turbulence dispersion and lift force play a key role in the concentration distribution. In addition, the heat transfer rate increases as the mass loading ratio increases and Nusselt number increases as the pipe diameter increases.     

Migration of settling particles in a horizontal viscous flow through a vertical slot with porous walls

Particle migration in a horizontal flow of dilute suspension through a vertical slot with porous walls is studied using the two-continua approach. The lateral migration is induced by two opposite effects: an inertial lift force due to particle settling and directed toward the slot centre-line, and a drag due to leak-off entraining particles toward the walls. An expression for the inertial lift on a settling particle in a horizontal channel flow found recently is generalized to the case of a low leak-off velocity. The evolution of an initial uniform particle concentration profile is studied within the full Lagrangian approach. Four migration regimes are found differing by the direction of particle migration and numbers of equilibrium positions. Conditions of the regime change and a critical value of dimensionless leak-off velocity for particle deposition on the walls are obtained analytically. Suspension flows with zones where the particle concentration is zero or increases infinitely, are studied numerically.     

Study of direct force measurement and characteristics on blades of Savonius rotor at static state

Experiments were conducted to measure and to clarify quantitative variations of fluid force on three kinds of setting blade configuration on a Savonius rotor at the static state in the steady uniform flow. In the experiments, by splitting both the end-disks of the blades into two half parts, the two blades were independent on each other and the one blade attached to both the half end-disks was only connected to the axis rod of the rotor, and the other blade connecting to the other half end-disks was free from the axis rod. Hence, the drag and lift forces acting on the former blade were directly measured in a short time at every phase angle and every overlap ratio of the two blades using the force balance meter. The experiment for the rotor of the diameter 160 mm and the height 160 mm was conducted in the steady wind flow of the Reynolds number, about 0.64 × 10⁵ which based on the rotor diameter and the upstream velocity. From the experiment, value of the drag coefficient of the one blade in two blades with the overlap ratio of zero was taken from about −0.3 to 1.6, and the lift coefficient of it from about −1.0 to 1.4 in a cycle. The minus value of drag only appears in the setting of one blade in two blades, however the minus value does not exist in the distribution of tangential force on every overlap. And it was shown that the force distributions on the one blade in a cycle were very different in the first and second half cycle and interfered by the other blade, and strongly dependent on the overlap ratio between the blades.     

Particle statistics in a gas–solid coaxial strongly swirling flow: A direct numerical simulation

A direct numerical simulation of a strongly coaxial swirling particle-laden flow is conducted with reference to a previous experiment. The carrier phase is simulated as a coaxial swirling flow through a short nozzle injecting into a large container. The particle phase is carried by the primary jet, and simulated in the Lagrangian approach. The drag force, slip-shear force and slip-rotation force experienced by particles are calculated. A partial validation of the results is followed. The results are analyzed in Eulerian approach focusing on the statistical behavior of particle motion. The relative importance of the drag, slip-shear and slip-rotation forces under different Stokes numbers is indicated quantitatively. The particle velocity profiles, fluctuations, Reynolds stress, and turbulence intensity are demonstrated and analyzed respectively. An important “choke” behavior for large particles within the mainstream is found and interpreted. Additionally, the patterns of particle distribution and the helical structures of particle motion under different Stokes numbers are demonstrated qualitatively and analyzed quantitatively.