Sedimentation of Two Side-by-Side Heavy Particles of Different Density in a Shear-Thinning Fluid with Viscoelastic Properties

Particle sedimentation has widely existed in nature and engineering fields, and most carrier fluids are non-Newtonian. Recently, the manipulation of a settling particle in liquid has been a topic of high interest to those involved in engineered processes such as composite materials, pharmaceutical manufacture, chemistry and the petroleum industry. Compared with Newtonian fluid, the viscosity of non-Newtonian fluid is closely related to the shear rate, leading to a single settling particle having different dynamic behaviors. In this article, the trajectories and velocities of two side-by-side particles of different densities (heavy and light) settling in a shear-thinning fluid with viscoelastic property were studied, as well as that for the corresponding single settling particle. Regardless of the difference in the particle density, the results show the two-way coupling interaction between the two side-by-side settling particles. As opposed to a single settling particle, the wake of the heavier particle can clearly attract or rebound the light particle due to the shear-thinning or viscoelastic property of the fluid. Regarding the trajectories of the light particle, three basic path types were found: (i) the light particle is first attracted and then repelled by the wake of the heavy one; (ii) the light particle approaches and then largely traces within the path of the heavy one in the limited field of view; (iii) the light particle is first slightly shifted away from its original position and then returns to this initial position. In addition to this, due to the existence of a corridor of reduced viscosity and negative wake generated by the viscoelastic property, the settling velocity of a light particle can exceed the terminal velocity of a single particle of the same density. On the other hand, the sedimentation of the light particle can induce the distinguishable transverse migration of the heavy one.

A velocity corrected unresolved CFD-DEM coupled method to reproduce wake effects at moderate Reynolds number

Purpose The purpose of this paper is to develop a corrected unresolved CFD-DEM method that can reproduce the wake effects in modeling particulate flows at moderate Reynolds number. Design/methodology/approach First, the velocity field in the wake behind a settling particle is numerically investigated by a resolved method, in which the finite volume method (FVM) is applied to model the fluid flow, discrete element method (DEM) is applied to simulate the motion of particles and immersed boundary method (IBM) is used to tackle fluid solid interaction. Second, an analytical scaling law is given, which can effectively describe the velocity field in the wake behind the settling particle at low and middle Reynolds numbers. Third, this analytical expression is incorporated into unresolved modeling to correct the relative velocity between the particle and its surrounding fluid and enable the influence of the wake of the particle on its neighboring particles. Findings Two numerical examples, the sedimentation of dual particles, a list of particles and even more particles are provided to show the effectiveness of the presented velocity corrected unresolved method (VCUM). It is found that, in both examples simulated with VCUM, the relative positions of the particles changed, and drafting & kissing phenomenon and particle clustering phenomenon were clearly observed. Practical implications The developed VCUM can be highly beneficial for modeling industrial particulate flows with DKT and particle clustering phenomena. Originality/value VCUM innovatively incorporates the wake effects into unresolved CFD-DEM method. It improves the computational accuracy of conventional unresolved methods with comparable results from resolved modeling, while the computational cost is greatly reduced.

Inertial drag on a sphere settling in a stratified fluid

We compute the drag force on a sphere settling slowly in a quiescent, linearly stratified fluid. Stratification can significantly enhance the drag experienced by the settling particle. The magnitude of this effect depends on whether fluid-density transport around the settling particle is due to diffusion, to advection by the disturbance flow caused by the particle or due to both. It therefore matters how efficiently the fluid disturbance is convected away from the particle by fluid-inertial terms. When these terms dominate, the Oseen drag force must be recovered. We compute by perturbation theory how the Oseen drag is modified by diffusion and stratification. Our results are in good agreement with recent direct numerical simulation studies of the problem.

Experimental Investigation of Flow Field Past a Spherical Particle Settling in Viscoelastic Fluids Using Particle Image Velocimetry Technique

Experimental investigation of flow field past a spherical particle settling in viscoelastic fluids using particle image shadowgraphy techniques studies have shown that the settling velocity of particles in viscoelastic fluids decreased significantly with the increasing elasticity of the fluids. However, our understanding of how and why the change in fluid elasticity influences the particle settling velocity are not yet fully developed. An experimental study, therefore, has been conducted to understand the reasons behind why the settling velocity of the particles decrease with the increasing fluid elasticity. The main objectives were: (i) to investigate the fluid flow field behind the settling particle by using particle image velocity (PIV) technique; (ii) to understand the changes caused by the elasticity of the fluid on the flow field past the settling particle; (iii) more specifically, to determine how the fluid velocity profile and the resultant drag forces acting on the settling particle change with the increasing fluid elasticity. Two different viscoelastic fluids were formulated by mixing 3 grades of HPAM polymer (MWs: 500,000; 8,000,000; 20,000,000; concentrations: 0.09% and 0.1%wt). The fluids were designed to have almost identical shear viscosity but significantly different elastic properties. The shear viscosity and elasticity of the fluids were determined by performing shear viscosity and frequency sweep oscillatory measurements, respectively. The settling velocities of the spherical particles in viscoelastic polymer fluids were measured by using particle image shadowgraph technique. The fluid flow field behind the settling particle was determined by using the PIV technique. Results of the PIV measurements demonstrated that negative wakes were present in viscoelastic fluids. The stagnation point (i.e. the point where the velocity becomes zero and above that the fluid starts moving in the direction opposite to the particle movement) was closer to the particle settling in the higher elasticity fluid than that in the lower elasticity fluid. The velocity of the fluid in the recirculation region was higher for the flow of the fluid with higher elasticity. The presence of negative wakes having fast moving fluid in the reverse direction near the settling particle possibly creates an additional drag force (acting on the particle in the direction opposite the particle movement), which would eventually slow down the settling particle. Knowledge of the settling behavior of particles is indispensable to design and optimize numerous industrial operations such as cuttings transport in oil and gas well drilling and proppant transport in hydraulic fracturing. In this study, by conducting experiments under controlled conditions, we were able to show how the change in fluid elasticity influenced the particle settling velocity. The results from this fundamental study can be used for development of optimum drilling and fracturing fluid formulations for effective transport of cuttings and proppants.

A numerical study of the dynamics of a particle settling at moderate Reynolds numbers in a linearly stratified fluid

In this paper, the transient settling dynamics of a spherical particle sedimenting in a linearly stratified fluid is investigated by performing fully resolved direct numerical simulations. The settling behaviour is quantified for different values of Reynolds, Froude and Prandtl numbers. It is demonstrated that the transient settling dynamics is correlated to the induced Lagrangian drift of flow around the settling particle. A simplified model is provided to predict the maximum velocity of the settling particle in linearly stratified fluids. The peak velocity can be followed by the oscillation of the settling velocity and the particle can even reverse its direction of motion before reaching to its neutrally buoyant level. The frequency of oscillation of settling velocity scales with the Brunt–Väisälä frequency and the motion of the particle can lead to the formation of secondary and tertiary vortices following the primary vortex.