Kinetic theory of colloidal suspensions: morphology, rheology, and migration

Smoluchowski kinetic equation governing the time evolution of the pair correlation function of rigid sphericalparticles suspended in a Newtonian fluid is extended to include particle migration. The extended kinetic equation takes into account three types of forces acting on the suspended particles: a direct force generated by an interparticle potential, hydrodynamic force mediated by the host fluid, and the Faxén-type forces bringing about the across-the-streamline particle migration. For suspensions subjected to externally imposed flows, the kinetic equation is solved numerically by the proper generalized decomposition method. The imposed flow investigated inthe numerical illustrations is the Poiseuille flow. Numerical solutions provide the morphology (the pair correlation function), the rheology (the stress tensor), and the particle migration.     

Wall effects for motion of spheres in power-law fluids

The steady motion of spheres representing particles inside tubes filled with different fluids has been investigated using both a finite-element and a finite-volume method. The rheology of the fluids has been modelled by the power-law able to describe the shear-thinning (pseudoplastic) behaviour of a series of polymer solutions. New results have been obtained for a series of tube/sphere diameter ratios in order to investigate the wall effects on the drag exerted by the fluid on the sphere. The results agree well with previous simulations for an unbounded medium (infinite diameter ratio). Experimental investigations have also been carried out and simulated, and the results compare favourably with the experiments. The present simulations revealed the convergence of the drag coefficient to a constant value independent of tube-to-sphere diameter ratio when the power-law index approaches zero.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Slow Motion of an Assemblage of Porous Spherical Shells Relative to a Fluid

The body-force-driven motion of a homogeneous distribution of spherically symmetric porous shells in an incompressible Newtonian fluid and the fluid flow through a bed of these shell particles are investigated analytically. The effect of the hydrodynamic interaction among the porous shell particles is taken into account by employing a cell-model representation. In the limit of small Reynolds number, the Stokes and Brinkman equations are solved for the flow field around a single particle in a unit cell, and the drag force acting on the particle by the fluid is obtained in closed forms. For a suspension of porous spherical shells, the mobility of the particles decreases or the hydrodynamic interaction among the particles increases monotonically with a decrease in the permeability of the porous shells. The effect of particle interactions on the creeping motion of porous spherical shells relative to a fluid can be quite significant in some situations. In the limiting cases, the analytical solution describing the drag force or mobility for a suspension of porous spherical shells reduces to those for suspensions of impermeable solid spheres and of porous spheres. The particle-interaction behavior for a suspension of porous spherical shells with a relatively low permeability may be approximated by that of permeable spheres when the porous shells are sufficiently thick.     

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

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

Pipe flow with solid particles fixed in space

A group of solid particles were hung by slender rods in a pipe to make a model of two-phase flow of coarse particles. Pressure gradient and velocities were measured for different types of the models. The drag on the particles (spheres) were obtained from measurements of pressure gradient with some assumptions. The results are summarized as follows. (1) Mean velocities of fluid are lower in the central part of the pipe than in the circumferential part. Turbulence is remarkably increased by particles. The spectrum distribution of turbulent velocity becomes flatter. These results are similar to the gas-solid flow of coarse particles in a vertical pipe. (2) At a large Reynolds number, the drag coefficient per one sphere in the group is larger than that of a single isolated sphere in a uniform flow. When the spheres are arranged along the same line in the longitudinal direction, the drag coefficient becomes smaller as the longitudinal distance between the spheres is shortened.     

Freely rising light solid spheres

This paper examines the behavior of spheres rising freely in a Newtonian fluid when the ratio between the density of the spheres and that of the surrounding fluid is about 0.02. High-speed imaging is used to reconstruct three-dimensional trajectories of the rising spheres. From the analysis of the trajectories the magnitudes of the drag and lift forces exerted by the surrounding fluid are deduced. It is argued that the two main contributions to the drag force are (i) a viscous drag that may be estimated from the standard drag curve by evaluating the Reynolds number using the actual value of the velocity, and (ii) an inertial drag that arises essentially by the same mechanisms that cause the lift-induced drag familiar from wing theory. Estimates of both contributions, the latter using visualizations of the wakes of the spheres, give a favorable agreement with the measured drag forces. These findings are closely related to recent numerical results of in the literature on the forces experienced by oblate ellipsoidal bubbles rising in quiescent water.     

垂直湍流液-固流中大颗粒的相对速度

通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式．用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数．实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。     

Development and validation of a reduced order history force model

The history force model accounts for temporal development in fluid gradients in the viscous region surrounding a particle in point particle methods. The calculation of the history force typically requires storing and using relative velocity information during the life time of the particle. For a large number of particles integrated over large times, history force calculation can become prohibitively expensive. The current work presents a new modeling approach to calculate the history force in which a decay function is applied to a stored cumulative value of the history force. The proposed formulation is equivalent to applying the same function obtained from a constant acceleration assumption to a running average of the acceleration within the memory time of the particle. The new force model is validated with experimental measurements of settling spheres at Reynolds numbers ranging from around one to a few hundreds and at density ratios from 1.2 to about 9.32. More validation work was carried-out with experimental measurements of oscillating spheres at different frequencies and amplitudes, as well as bouncing spheres at different Reynolds numbers and density ratios. The model shows very good agreement with the experiments of settling spheres and reasonable/good agreement with oscillating and bouncing sphere experiments. The proposed model significantly reduces the computational resources required to calculate the history force especially when large number of particles need to be integrated over long times.     

A sharp interface Cartesian grid method for viscous simulation of shocked particle-laden flows

A Cartesian grid-based sharp interface method is presented for viscous simulations of shocked particle-laden flows. The moving solid–fluid interfaces are represented using level sets. A moving least-squares reconstruction is developed to apply the no-slip boundary condition at solid–fluid interfaces and to supply viscous stresses to the fluid. The algorithms developed in this paper are benchmarked against similarity solutions for the boundary layer over a fixed flat plate and against numerical solutions for moving interface problems such as shock-induced lift-off of a cylinder in a channel. The framework is extended to 3D and applied to calculate low Reynolds number steady supersonic flow over a sphere. Viscous simulation of the interaction of a particle cloud with an incident planar shock is demonstrated; the average drag on the particles and the vorticity field in the cloud are compared to the inviscid case to elucidate the effects of viscosity on momentum transfer between the particle and fluid phases. The methods developed will be useful for obtaining accurate momentum and heat transfer closure models for macro-scale shocked particulate flow applications such as blast waves and dust explosions.     

非球形颗粒两相流的数值模拟研究进展

非球形颗粒两相流是多相流的重要研究方向之一, 常见于自然界及工业生产过程中. 不同于球形颗粒, 由于非球形颗粒形状的各向异性, 除了颗粒平动行为, 还需要考虑颗粒的转动与取向行为, 颗粒的取向与转动行为会影响颗粒所受的力和力矩. 为了准确模拟非球形颗粒的运动行为, 目前非球形颗粒两相流的数值模拟研究主要基于欧拉?拉格朗日的求解框架展开, 常见的非球形颗粒两相流数值模拟方法主要包括点颗粒法与全分辨颗粒法. 本文将对这两类方法进行介绍, 同时会全面介绍非球形颗粒两相流研究的基础理论模型, 并系统总结非球形颗粒在简单基本流和复杂湍流中的研究进展, 包括对于非球形颗粒在湍流中的取向与转动行为机理, 以及颗粒对湍流减阻调制作用的研究. 最后, 本文提出了非球形颗粒两相流研究存在的问题及未来研究方向.      

细长椭球体在水中自由下落的运动特性实验研究

固体颗粒在液体中的运动现象在日常生活和工业应用领域广泛存在, 其中因蕴含着丰富的流体力学现象而受到学者们的广泛关注. 本文通过实验研究了细长椭球体在水中受浮力影响的下落特性. 实验中采用带有两台相互垂直的高速摄像机和光源组成的运动跟踪平台并结合荧光染色技术对细长椭球体下落过程中的运动轨迹和尾涡结构进行研究. 文中选用的细长椭球体与环境流体的密度比为1.2, 其长短轴比范围为2$\sim $10, 相应的阿基米德数范围为400$\sim$1400, 对应实现的终态雷诺数范围为120$\sim$1350. 实验过程中我们观察到细长椭球体在水中下落过程中产生的5种典型路径, 分别为: 小振幅不规则运动、小振幅高频振荡运动、大振幅低频振荡运动、高度非线性运动以及直线运动, 并得到了对应的速度振荡以及倾斜角的演化规律. 进一步地, 分析了细长椭球体运动过程中受到的阻力系数与雷诺数之间的关系. 随后采用荧光可视化技术清晰获得了颗粒下落过程中的尾涡结构特性, 并结合颗粒的运动状态详细分析了涡脱落过程对颗粒运动状态转捩的影响. 最后, 通过对比前人关于圆柱体下落的运动特性的相关结果, 获得了细长椭球体和细长圆柱体运动特性之间的异同点以及其潜在的物理机理.      

Linear stability of particle laden flows: the influence of added mass,fluid acceleration and Basset history force

Both modal and non-modal linear stability analysis of a channel flow laden with particles is presented. The particles are assumed spherical and solid and their presence modelled using two-way coupling, with Stokes drag, added mass and fluid acceleration as coupling terms. When the particles considered have a density ratio of order one, all three terms become important. To account for the volume and mass of the particles, a modified Reynolds number is defined. Particles lighter than the fluid decrease the critical Reynolds number for modal stability, whereas heavier particles may increase the critical Reynolds number. Most effect is found when the Stokes number defined with the instability time scale is of order one. Non-modal analysis shows that the generation of streamwise streaks is the most dominant disturbance-growth mechanism also in flows laden with particles: the transient growth of the total system is enhanced proportionally to the particle mass fraction, as observed previously in flows laden with heavy particles. When studying the fluid disturbance energy alone, the optimal growth hardly changes. We also show that the Basset history force has a negligible effect on stability. The inclusion of the extra interaction terms does not show any large modifications of the subcritical instabilities in wall-bounded shear flows.     

CFD-DEM simulation of the hole cleaning process in a deviated well drilling: The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid–particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle−particle, particle−wall, and particle−drill pipe are taken into account with the Hertz–Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the multi-sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid–particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid–solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.     

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.      

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.     

Slow motion of a slip spherical particle perpendicular to two plane walls

The problem of the quasisteady motion of a spherical fluid or solid particle with a slip-flow surface in a viscous fluid perpendicular to two parallel plane walls at an arbitrary position between them is investigated theoretically in the limit of small Reynolds number. To solve the axisymmetric Stokes equation for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both circular cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the plane walls by the Hankel transform and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle and of the relative separation distances between the particle and the confining walls. For the motions of a spherical particle normal to a single plane wall and of a no-slip sphere perpendicular to two plane walls, our drag results are in good agreement with the available solutions in the literature for all relative particle-to-wall spacings. The boundary-corrected drag force acting on the particle in general increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given geometry, but there are exceptions. For a specified wall-to-wall spacing, the drag force is minimal when the particle is situated midway between the two plane walls and increases monotonically when it approaches either of the walls. The boundary effect on the particle motion normal to two plane walls is found to be significant and much stronger than that parallel to them.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.