Efficient calculation of the history force at finite Reynolds numbers

Well-acknowledged problems associated with modeling the history force in large, many-particle simulations are related to the need to store and integrate over the entire lifetime of the particle. To address this concern, a computationally efficient method for calculating the history force (the “window model”) was developed based on the assumption of weak changes in acceleration in the recent relevant history of the particle. This assumption leads to the design of a model with a truncated integration interval which requires storage of and integration over a much shorter period of the particle’s history compared to other history force models. The truncation of the integration window can yield more than an order of magnitude savings in CPU time. In a related study, the two empirical coefficients of the Mei & Adrian history force kernel have been optimized (based on comparison with experimental data for falling particles) to give improved predictions of the data. Both the new history force kernel and the window model have been investigated for a large range of experimental data yielding, to the authors’ knowledge, the most extensive comparison yet conducted. For falling particles, the new history force kernel shows good predictions for particle Reynolds numbers ranging from 9 to 853 and density ratios from 1.17 to 9.32. Good predictions were also obtained using the window model when changes in particle relative acceleration over the window period were modest. For particles under forced oscillating in a quiescent fluid, the history kernel was generally reasonable but did not predict the peak forces well in all cases. This may be explained by noting that the assumption of a t⁻² long-time dependence for the finite Reynolds number history force kernel may become invalid during rapid deceleration and wake ingestion (which can lead to exponential or t⁻¹ behavior). However, the finite Reynolds number kernel gives better predictions in all cases than those made using the Basset history force. The window model was only reasonable for the oscillating particle cases when the changes in the relative particle acceleration over the integration window were small.     

Dynamics of dual-particles settling under gravity

As a first step towards understanding particle–particle interaction in fluid flows, the motion of two spherical particles settling in close proximity under gravity in Newtonian fluids was investigated experimentally for particle Reynolds numbers ranging from 0.01 to 2000. It was observed that particles repel each other for Re>0.1 and that the separation distance of settling particles is Reynolds number dependent. At lower Reynolds numbers, i.e. for Re<0.1, particles settling under gravity do not separate.The orientation preference of two spherical particles was found to be Reynolds number dependent. At higher Reynolds numbers, the line connecting the centres of the two particles is always horizontal, regardless of the way the two particles are launched. At lower Reynolds numbers, however, the particle centreline tends to tilt to an arbitrary angle, even of the two particles are launched in the horizontal plane. Because of the tilt, a side migration of the two particles was found to exist. A linear theory was developed to estimate the side migration velocity. It was found that the maximum side migration velocity is approximately 6% of the vertical settling velocity, in good agreement with the experimental results.Counter-rotating spinning of the two particles was observed and measured in the range of Re=0–10. Using the linear model, it is possible to estimate the influence of the tilt angle on the rate of rotation at low Reynolds numbers. Dual particles settle faster than a single particle at small Reynolds numbers but not at higher Reynolds numbers, because of particle separation. The variation of particle settling velocity with Reynolds number is presented. An equation which can be used to estimate the influence of tilt angle on particle settling velocity at low Reynolds number is also derived.     

Dynamic simulation of sedimentation of solid particles in an Oldroyd-B fluid

In this paper we present a two-dimensional numerical study of the viscoelastic effects on the sedimentation of particles in the presence of solid walls or another particle. The Navier-Stokes equations coupled with an Oldroyd-B model are solved using a finite-element method with the EVSS formalism, and the particles are moved according to their equations of motion. In a vertical channel filled with a viscoelastic fluid, a particle settling very close to one side wall experiences a repulsion from the wall; a particle farther away from the wall is attracted toward it. Thus a settling particle will approach an eccentric equilibrium position, which depends on the Reynolds and Deborah numbers. Two particles settling one on top of the other attract and form a doublet if their initial separation is not too large. Two particles settling side by side approach each other and the doublet also rotates till the line of centers is aligned with the direction of sedimentation. The particle-particle interactions are in qualitative agreement with experimental observations, while the wall repulsion has not been documented in experiments. The driving force for lateral migrations is shown to correlate with the pressure distribution on the particle's surface. As a rule, viscoelasticity affects the motion of particles by modifying the pressure distribution on their surface. The direct contribution of viscoelastic normal stresses to the force and torque is not important.     

Size segregation and particle velocity fluctuations in settling concentrated suspensions

We investigate the sedimentation of concentrated suspensions at low Reynolds numbers to study collective particle effects on local particle velocity fluctuations and size segregation effects. Experiments are carried out with polymethylmetacrylate (PMMA) spheres of two different mean diameters (190 and 25 μm) suspended in a hydrophobic index-matched fluid. Spatial repartitions of both small and large spheres and velocity fluctuations of particles are measured using fluorescently labelled PMMA spheres and a particle image velocimetry method. We also report measurements of the interstitial fluid pressure during settling. Experiments show that size segregation effects can occur during the sedimentation of concentrated suspensions of either quasi-monodisperse or bidisperse spheres. Size segregation is correlated to the organisation of the sedimentation velocity field into vortex-like structures of finite size. A loss of size segregation together with a significant decrease of the fluid pressure gradient in the bulk suspension is observed when the size of vortex-like structures gets on the order of the container size. However, the emergence of channels through the settling zone prevents a complete loss of size segregation in very concentrated suspensions.     

颗粒在黏弹性流体扩张和收缩流中沉降的格子Boltzmann模拟分析

对于Oldroyd-B型黏弹性流体,本文应用格子Boltzmann方法(LBM),实现了流体在二维1:3扩张流道及3:1收缩流道中流动的数值模拟,获得了黏弹性流体在扩张和收缩流道中的流场分布．结合颗粒的受力和运动规则,基于点源颗粒模型,数值分析了颗粒在扩张流和收缩流中的沉降过程和特征,讨论了颗粒相对质量和起始位置以及雷诺数Re和威森伯格数Wi对颗粒沉降特征的影响．结果表明,颗粒相对质量和起始位置以及Re对颗粒沉降轨迹和落点影响较大,而Wi的影响则较小．     

DRAG FORCE IN DENSE GAS—PARTICLE TWO—PHASE FLOW

Numerical simulations of flow over a stationary particle in a dense gas-particle two-phase flow have been carried out for small Reynolds numbers (less than 100). In order to study the influence of the particles interaction on the drag force, three particle arrangements have been tested: a single particle, two particles placed in the flow direction and many particles located regularly in the flow field. The Navier-Stokes equations are discretized in the three-dimensional space using finite volume method. For the first and second cases, the numerical results agree reasonably well with the data in literature. For the third case, i.e., the multiparticle case, the influence of the particle volume fraction and Reynolds numbers on the drag force has been investigated. The results show that the computational values of the drag ratio agree approximately with the published results at higher Reynolds numbers (from 34.2 to 68.4), but there is a large difference between them at small Reynolds numbers. The project supported by the Special Funds for Major State Basis Research Projects in China (G19990222).     

Acceleration of heavy particles in isotropic turbulence

The paper concerns the effect of particle inertia on acceleration statistics. A simple analytical model for predicting the acceleration of heavy particles suspended in an isotropic homogeneous turbulent flow field is developed. This model is capable of describing the influence of both Stokes and Reynolds numbers on the particle acceleration variance. Comparisons of model predictions with numerical simulations are presented.     

Second-order finite difference solutions for the flow between rotating concentric spheres

This paper describes a second-order method to calculate approximate solutions to flow of viscous incompressible fluid between rotating concentric spheres. The governing partial differential equations are presented in the stream–vorticity formulation and are written as a series of second-order equations. The technique employed makes use of second-order approximations for all terms in the governing equations and is dependent upon the direction of flow at a given point. This upwind technique has allowed us to generate approximate solutions with larger Reynolds numbers than has generally been possible for second and higher-order techniques. Solutions have been obtained with Reynolds numbers as large as 3000 and with grids as fine as a 40 × 40 mesh. Results are displayed in the form of level curves for both the stream and vorticity functions. A dimensionless quantity related to the torque acting on both spheres has been calculated from the approximate solution and compared with other results. Results with smaller Reynolds numbers such as 100 and 1000 are in excellent agreement with other published results.     

The sedimentation of polydisperse suspensions in vessels having inclined walls

A theory is presented for describing the sedimentation of polydisperse suspensions in two-dimensional channels having walls that are inclined to the vertical. The theory assumes that the flow is laminar and that the suspension consists of spherical beads having small particle Reynolds numbers. The suspension may consist of either $\text{N}$ distinct species of particles or of a continuum of particle sizes and densities. For the sake of simplicity, the analysis is mostly confined to the case in which the hindered settling velocity of each particle is given by its Stokes settling velocity multiplied by a function of the total local solids concentration. Under these conditions, results are developed that are useful for the design of either batch or continuous settling devices. Experimental observations were found to be in good agreement with the predictions of the present theory.     

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%.     

Oldroyd-B黏弹性液滴弹跳行为的改进SPH模拟

基于光滑粒子流体动力学SPH(Smoothed Particle Hydrodynamics)方法对Oldroyd-B黏弹性液滴撞击固壁面产生的弹跳行为进行了模拟与分析。首先,为了解决SPH模拟黏弹性自由表面流出现的张力不稳定性问题,联合粒子迁移技术提出了一种改进SPH方法。然后,对Oldroyd-B黏弹性液滴撞击固壁面产生的铺展行为进行了改进SPH模拟,与文献结果的比较验证了方法的有效性。最后,通过降低Reynolds数捕捉到了液滴的弹跳行为;并在此基础上,分析了液滴黏度比、Weissenberg数和Reynolds数对液滴弹跳行为的影响。结果表明,改进SPH方法可有效地模拟黏弹性自由表面流问题;液滴黏度比、Weissenberg数和Reynolds数对液滴最大回弹高度均有显著的影响。     

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.     

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.     

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

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

Flow induced oscillation of two rigid rectangular plates in a side-by-side configuration

The present work is an experimental study of two oscillating rigid plates placed in side-by-side configuration, hinged at their leading edges, subjected to low subsonic flow. This problem is investigated using smoke-wire flow visualization, hot-wire anemometry, and time resolved particle image velocimetry. It is found that beyond a critical Reynolds number, the plates set into oscillatory motion. This critical Reynolds number depends on the gap between the plates. It is also seen that this value of Reynolds number, at lower values of gap to thickness ratio (<7) is significantly higher than that of the single plate configuration value. The frequency and amplitude of the oscillating plates at various gaps and Reynolds numbers have been studied and compared with the characteristics of an oscillating single plate. It is also found that depending on the gap and acceleration of the free-stream, there exist two modes of oscillation - (i) in-phase and (ii) out-of-phase. For gap to thickness ratio less than 10, only in-phase oscillations take place for all values of free-stream velocity considered in the present work, whereas, when this ratio is greater than 10, the mode of oscillation depends on the initial conditions up to a certain free-stream velocity, beyond which the plates switch to in-phase mode. Smoke wire flow visualization technique along with time resolved particle image velocimetry reveal that the vorticity distributions around the plates are responsible for the initiation of the two modes of oscillations.     

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.     

The motion of rodlike particles in the pressure-driven flow between two flat plates

The motion of freely suspended rodlike particles has been observed in the pressure-driven flow between the two flat plates of a Hele Shaw flow cell at low Reynolds numbers. Data are reported for rodlike particles with aspect ratios of 12.0 suspended in a Newtonian fluid for gap thickness to particle length ratios of 3, 6, and 20; and for rodlike particles with aspect ratios between 5 and 8 in a non-Newtonian fluid (79.25 wt.% water, 20.2 wt.% glycerine, and 0.55 wt.% polyacrylamide). For the Newtonian fluid, the time-dependent orientation of the particles near and far from walls was shown to be in quantitative agreement with Jeffery's theory for ellipsoids suspended in a simple shear flow if an effective aspect ratio is calculated from the experimental period of rotation. Particles aligned with the flow direction and less than a particle half-length from a wall interacted irreversibly with the wall. For the non-Newtonian fluid, the timedependent orientation far from a wall was shown to be in qualitative agreement with Leal's theory for a second-order fluid; however, particles that were aligned with the flow direction and were near walls did not rotate.     

Überlegungen zur Scherviskosität von Suspensionen aufgrund einer Dimensionsanalyse

Dimensional analysis of the motion of solid particles suspended in a fluid phase shows that the macroscopic relative shear viscosity of suspensions generally depends not only on the volume concentration and particle shape but also on two Reynolds numbers and a dimensionless sedimentation number. These dimensionless numbers are formed using parameters characterizing the structure and motion of the suspension at the microscopic level. The analysis was based on the assumptions that the dispersed particles are rigid and sufficiently large that Brownian motion may be neglected, that the continuous fluid phase is Newtonian and that the interactions between particles and between particles and fluid phase are only hydrodynamic. The Reynolds numbers describe the influence of the inertial forces at the microscopic level, and the sedimentation number the influence of gravity. The dimensionless numbers can be neglected if their values are much smaller than one. For each of the dimensionless numbers both the shear rate and the particle size influence the shear viscosity. Thus sedimentation number is large for low shear rates, whereas the Reynolds numbers are large for high shear rates. The viscosity function for one suspension can be transformed into the viscosity function for another suspension with geometrically similar particles but of a different size. The scale-up rules are derived from the requirement that the relevant dimensionless numbers must be constant. The influence of non-hydrodynamic effects at the microscopic level on the shear viscosity can be detected by deviations from the derived scale-up rules.     

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.