ضياء عبد الحسين جاسم الغالبي

We study the sedimentation of finite-size particles in quiescent wall-bounded Newtonian and shear-thinning fluids. The problem is studied numerically by means of direct numerical simulations with the presence of the particles accounted for with an immersed boundary method. The suspended phase consists of Non-Brownian rigid spherical particles with particle to fluid density ratio $\rho_p/\rho_f=1.5$; three different solid volume fractions $\Phi= 1\%$, $5\%$ and $20\% $ are considered. The Archimedes number is kept constant to $Ar=36$ for all shear-thinning fluids, while it is changed to $Ar=97$ for the Newtonian fluid to reproduce the same terminal velocity of a single particle sedimenting in the shear-thinning fluid. We show that the mean settling velocities decrease with the particle concentration as a consequence of the hindering effect and that the mean settling speed is always larger in the shear thinning fluid than in the Newtonian one. This is due to the decrease of the mean viscosity of the fluid which leads to a lower drag force acting on the particles. We show that particles tend to form aggregates in the middle of the channel in a shear-thinning fluid, preferentially positioning in the wake of neighboring particles or aside them, resulting in lower levels of fluid velocity fluctuation in the gravity direction than in a Newtonian fluid.

We perform fully Eulerian numerical simulations of an initially spherical hyperelastic particle suspended in a Newtonian pressure-driven flow in a cylindrical straight pipe. We study the full particle migration and deformation for different Reynolds numbers and for various levels of particle elasticity, to disentangle the interplay of inertia and elasticity on the particle focusing. We observe that the particle deforms and undergoes a lateral displacement while traveling downstream through the pipe, finally focusing at the pipe centerline. We note that the migration dynamics and the final equilibrium position are almost independent of the Reynolds number, while they strongly depend on the particle elasticity; in particular, the migration is faster as the elasticity increases (i.e., the particle is more deformable), with the particle reaching the final equilibrium position at the centerline in shorter times. Our simulations show that the results are not affected by the particle initial conditions, position, and velocity. Finally, we explain the particle migration by computing the total force acting on the particle and its different components, viscous and elastic.

Solid or deformable particles suspended in a viscous fluid are of scientific and technological interest in a broad range of applications. Pyroclastic flows from volcanoes, sedimentation flows in river bed, food industries, oil-well drilling, as well as blood flow in the human body and the motion of suspended micro-organisms in water (like plankton) are among the possible examples. Often, in these particulate flows, the carrier fluid might exhibit an inelastic or a visco-elastic non-Newtonian behavior. Understanding the behavior of these suspensions is a very difficult task. Indeed, the complexities and challenges of multiphase flows are mainly due to the large number of governing parameters such as the physical properties of the particles (e.g., shape, size, stiffness, density difference with suspended fluid, solid volume fraction), the large set of interactions among particles and the properties of the carrier fluid (Newtonian or non-Newtonian); variations of each of these parameters may provide substantial quantitative and qualitative changes in the behavior of the suspension and affect the overall dynamics in several and sometimes surprising ways. The aim of this work is therefore to provide a deeper understanding of the behavior of particle suspensions in laminar Newtonian and non-Newtonian (inelastic and/or visco-elastic) fluid flows for a wide range of different parameters. To this purpose, particle-resolved direct numerical simulations of spherical particles are performed, using an efficient and accurate numerical tool. The code is based on the Immersed Boundary Method (IBM) for the fluid-solid interactions with lubrication, friction and collision models for the close range particle-particle (particle-wall) interactions. Both inelastic (Carreau and power-law), and visco-elastic models (Oldroyd-B and Giesekus) are employed to investigate separately the shear-thinning, shear-thickening, viscoelastic and combined shear-thinning visco-elastic features of the most commonly encountered non-Newtonian fluids. Moreover, a fully Eulerian numerical algorithm based on the one-continuum formulation is used to examine the case of an hyper-elastic neo-Hookean deformable particle suspended in a Newtonian flows. Firstly, we have investigated suspensions of solid spheres in Newtonian, shear thinning and shear thickening fluids in the simple shear flow created by two walls moving in opposite directions, considering various solid volume fractions and particle Reynolds numbers, thus including inertial effects. The results show that that the non-dimensional relative viscosity of of the suspension and the mean value of the local shear-rate can be well predicted by homogenization theory, more accurately for lower particle concentrations. Moreover, we show that in the presence of inertia, the effective viscosity of these suspensions deviates from that of Stokesian suspensions. We also examine the role of fluid elasticity, shear-thinning and combined shear-thinning visco-elastic effects on the simple linear Couette shear flow of neutrally-buoyant rigid spherical particles. It is found that the effective viscosity grows monotonically with the solid volume fraction and that all the Non-Newtonian cases exhibit a lower effective viscosity than the Newtonian ones; in addition, we show that elastic effects dominate at low elasticity whereas shear thinning is predominant at high applied shear rates. These variations in the effective viscosity are mainly due to changes in the particle-induced shear stress component. We then study the settling of spherical particles in quiescent wall-bounded Newtonian and shear-thinning fluids at three different solid volume fractions. We find that the mean settling velocities decrease with the particle concentration as a consequence of the hindering effect and thatthe mean settling speed is always larger in the shear thinning fluid than in the Newtonian one, due to the reduction of the local fluid viscosity around the particles which leads to a lower drag force acting on the particles. Finally, the inertial migration of hyper-elastic deformable particle in laminar pipe flows is also investigated. We consider different flow rates and various levels of particle elasticity. We observe that the particle deforms and experiences a lateral movement while traveling downstream through the pipe, always finding a stable position at the pipe centerline.

We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ( 0.1⩽Φ⩽0.4 ) and particle Reynolds numbers ( 0.1⩽Rep⩽10 ). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal Φ . This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. (J. Rheol., vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the Rep . The wide spectrum of the local shear rate and its dependency on Φ and Rep point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, O(Rep)∼10 .