A smoothed particle hydrodynamics-based fluid model with a spatially dependent viscosity: application to flow of a suspension with a non-Newtonian fluid matrix

A smoothed particle hydrodynamics approach is utilized to model a non-Newtonian fluid with a spatially varying viscosity. In the limit of constant viscosity, this approach recovers an earlier model for Newtonian fluids of Español and Revenga (Phys Rev E 67:026705, 2003). Results are compared with numerical solutions of the general Navier–Strokes equation using the “regularized” Bingham model of Papanastasiou (J Rheol 31:385–404, 1987) that has a shear-rate-dependent viscosity. As an application of this model, the effect of having a non-Newtonian fluid matrix, with a shear-rate-dependent viscosity in a moderately dense suspension, is examined. Simulation results are then compared with experiments on mono-size silica spheres in a shear-thinning fluid and for sand in a calcium carbonate paste. Excellent agreement is found between simulation and experiment. These results indicate that measurements of the shear viscosity of simple shear-rate-dependent non-Newtonian fluids may be used in simulation to predict the viscosity of concentrated suspensions having the same matrix fluid.     

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.     

The slow motion of two particles symmetrically placed about the axis of a circular cylinder in a direction perpendicular to their line of centers

Numerical values are provided which enable the trictional force and torque on the reference sphere to be computed for the particular case when two spheres move in a direction perpendicular to their line of centers symmetrically placed about the axis of a circular cylinder. Results for this motion are also expressed in terms of the ratio of frictional forces experienced by 1) a body of arbitrary shape in a bounded fluid with another particle and 2) the body now moving alone with the same speed and orientation in the same but unbounded fluid. The computation furnishes the interaction and wall corrections correct to the first order in the ratios of characteristic particle dimension to characteristic distance of the particle from another object. The theoretical results are compared with experimental data and found to be in excellent agreement.     

The flow of closely fitting particles in tapered tubes

The flow of rigid spheres, truncated cones and elastic incompressible spheres in tapered tubes is investigated assuming that the Reynolds equation is valid in the fluid and the linear theory of elasticity is applicable in the solid. It is shown that leading terms in the asymptotic expansion of pressure drop in terms of minimum fluid film thickness for neutrally buoyant rigid spheres and truncated cones are of higher order of magnitude compared to the corresponding terms for the flow of these particles in circular cylindrical tubes. The effect of taper angle on pressure drop is reduced in the case of soft elastic particles because of particle deformations and significant velocities at the particle surface.     

An elastohydrodynamic theory for the rheology of concentrated suspensions of deformable particles

The equations which govern thin films of a Newtonian liquid confined between deformable solid surfaces are applied to the regions of near contact in a concentrated suspension of deformable particles.For the case of slightly deformable elastic particles, one obtains the socalled “elastohydrodynamic” equations of lubrication theory.The appropriate asymptotic solution of these equations yields estimates for the viscosity, of a form proposed earlier by Frankel and Acrivos [1] for rigid particles, as well as a relaxation time for a suspension of near spheres. The present method, which goes beyond the dissipation calculation of Frankel and Acrivos to a derivation of the full stress tensor, yields the same form of dependence of viscosity on particle concentration. However, there is an as yet unexplained difference between the methods in the value of a numerical coefficient determined by the assumed packing of the spheres.While further work is needed on the kinetic theory for fluid suspensions, the methods employed here for the derivation of the stress tensor should have direct utility for certain solid dispersions, where it is possible to specify a priori the particle-packing in the system.     

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.     

Dynamic simulation of non-spherical particulate suspensions

Particle-level simulation has been employed to investigate rheology and microstructure of non-spherical particulate suspensions in a simple shear flow. Non-spherical particles in Newtonian fluids are modeled as three-dimensional clusters of neutrally buoyant, non-Brownian spheres linked together by Hookean-type constraint force. Rotne–Prager correction to velocity disturbance has been employed to account for far-field hydrodynamic interactions. An isolated rod-like particle in simple shear flow exhibits a periodic orientation distribution, commonly referred to as Jeffery orbit. Lubrication-like repulsive potential between clusters have been included in simulation of rod-like suspensions at various aspect ratios over dilute to semi-dilute volume fractions. Shear viscosity evaluated by orientation distribution qualitatively agrees with one obtained by direct computation of shear stress.     

Linear stability of a two-dimensional uniform fluidized bed

Basic equations in a two-dimensional fluidized bed are constructed for the particle and the fluid phases, and linear stability to two-dimensional disturbances of the volume fractions and the velocities of both phases is analyzed. The diffusion of particles and an effective viscosity in the particle phase are considered. It was found that the inertia term due to the average fluid velocity is responsible for the instability, while the particle diffusion and the effective particle viscosity suppress the growth of disturbances. It was also found that the most unstable state has a vertical wavenumber vector.     

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.     

Maximising the electrorheological effect for bidisperse nanofluids from the electrostatic force between two particles

A multipole re-expansion solution for two nonidentical dielectric spheres in a parallel electric field is used to determine the critical ratio of particle radii which leads to the strongest force of attraction between the spheres at various interstices and under varying dielectric properties. These critical ratios provide genuine optimal dimensions, in the sense that the force of attraction decreases for both increasing and decreasing ratios. Numerical results are compared with experimental results from the literature and discussed from the perspective of the impact on the design of electrorheological nanofluids.     

Finite element analysis of the axially symmetric motion of an incompressible viscous fluid in a spherical annulus

This paper presents results obtained by employing a modified Galerkin finite element method to analyse the steady state flow of a fluid contained between two concentric, rotating spheres. The spheres are assumed to be rigid and the cavity region between the spheres is filled with an incompressible, viscous, Newtonian fluid. The inner sphere is constrained to rotate about a vertical axis with a prescribed angular velocity, while the outer sphere is fixed. Results for the circumferential function Ω, streamfunction ψ, vorticity function ζ and inner boundary torque T₁ are presented for Reynolds numbers Re ? 2000 and radius ratios 0.1 ? α ? 0.9. The method proved effective for obtaining results for a wide range of radius ratios (0.1 ? α ? 0.9) and Reynolds numbers (0 ? Re ? 2000). Previous investigators who employed the finite difference method experienced difficulties in obtaining results for cases with radius ratios α ? 0.2, except for small Reynolds numbers (Re ? 100). Results for Ω, Ψ, ζ and T₁ obtained in this study for radius ratios 0.8 ≤ α ≤ 0.9 verified the development of Taylor vortices reported by other investigators. The research indicates that the method may be useful for analysing other non-linear fluid flow problems.     

Motion of a spherical particle in a viscous nonisothermal fluid

The Stokes and Hadamard-Riabouchinsky formulas are generalized to the case of steady motion of a solid spherical particle or drop in an incompressible fluid whose viscosity depends exponentially on the temperature. It is shown that for finite temperature differences between the surface of the particle and the region far from it the drag is determined by an effective viscosity with value close to the geometric mean of the viscosity on the surface of the particle and far from it.     

Numerical investigation of dynamical behavior of tethered rigid spheres in supersonic flow

The dynamical behavior of two tethered rigid spheres in a supersonic flow is numerically investigated. The tethered lengths and radius ratios of the two spheres are different. The two spheres, which are centroid axially aligned initially, are held stationary first, then released, and subsequently let fly freely in a supersonic flow. The mean qualities of the system and the qualities of the bigger sphere are considered and compared with the situations without the tether. In the separation process, six types of motion caused by the spheres, tether, and fluid interaction are found. The results show that the mean x-velocity of the system changes in a different manner for different radius ratios, and the x-velocity of the bigger sphere is uniformly reduced but through different mechanisms.     

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

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulating surface tension with smoothed particle hydrodynamics

A method for simulating two‐phase flows including surface tension is presented. The approach is based upon smoothed particle hydrodynamics (SPH). The fully Lagrangian nature of SPH maintains sharp fluid–fluid interfaces without employing high‐order advection schemes or explicit interface reconstruction. Several possible implementations of surface tension force are suggested and compared. The numerical stability of the method is investigated and optimal choices for numerical parameters are identified. Comparisons with a grid‐based volume of fluid method for two‐dimensional flows are excellent. The methods presented here apply to problems involving interfaces of arbitrary shape undergoing fragmentation and coalescence within a two‐phase system and readily extend to three‐dimensional problems. Boundary conditions at a solid surface, high viscosity and density ratios, and the simulation of free‐surface flows are not addressed. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamic simulation of arbitrarily shaped particles in shear flow

A new simulation framework was created for modeling the dynamics of arbitrarily shaped particles dispersed in Newtonian fluid. Theoretical complexity usually restricts suspension simulations to those for spheroids. This new simulation is loosely based on the Stokesian Dynamics method including long range hydrodynamic interaction and uses spheres as building components for greater particulates of arbitrary shape. This approach is capable of accurately reproducing the dynamics of an isolated arbitrarily shaped particle. As verification, the simulated results are compared against known results for a rod-like particle. An elongated rod-shaped structure made from linked spheres is shown to reproduce the well-known elongated ellipsoidal particle dynamics described by [Jeffrey Proc R Soc Lond A 102:161–179, 1923]. The predicted orbital period and spin rates for a fiber in shear are reproduced and compare well with theoretical prediction over a wide aspect ratio range. Predicted particle dynamics for other shaped particles are then demonstrated.     

Elongational flows of some non-colloidal suspensions

The rheology of a system must be explored not only in viscometric flows, but also in other flow classes, and so, we present some results for the axisymmetric elongational flow of non-colloidal suspensions of spheres. We compare our results with data from shear flows using the same matrices and spheres. We have experimented with non-colloidal suspensions of 40-μm diameter polystyrene spheres with volume fractions (?) varying from 0.3 to 0.5. Two matrix fluids were used—one was a near-Newtonian polydimethyl siloxane of 12 Pa-s viscosity and the other was a variant of the M1 Boger fluid sample of Sridhar which we call M1*. We did not find that the Trouton ratio for either of these fluids was 3; generally, the ratio was larger. We investigated the role of sphere roughness using spheres roughened to 5.3 % of the radius in a 50 % suspension in silicone oil and found an increase of elongational viscosity of about 65 % which is comparable with the 60 % increase in shear viscosity with roughness noted previously. For the silicone oil matrix, we found no rate effect, with very little strain-hardening. By contrast, the M1-type matrix suspensions showed strain-hardening and an increase of elongational viscosity with elongation rate.     

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.     

Studies on the axisymmetric sphere–sphere interaction problem in Newtonian and non-Newtonian fluids

In this research, experimental studies have been performed on the hydrodynamic interaction between two spheres by using particle image velocimetry and measuring the force between the spheres. To approach the system as a resistance problem, a servo-driving system was set-up by assembling a microstepping motor, a ball screw and a linear motion guide for the particle motion. Glycerin and a dilute solution of polyacrylamide in glycerin were used as Newtonian and non-Newtonian fluids, respectively. The polymer solution behaves like a Boger fluid when the concentration is 1000 ppm or less. The experimental results were compared with the asymptotic solution of Stokes equation. The result shows that fluid inertia and unsteadiness play important roles in the particle–particle interaction in the Newtonian fluid. This implies that the motion of two particles in suspension is not reversible even in the Newtonian fluid. In the non-Newtonian fluid, in addition to inertial effect, normal stress differences and viscoelasticity play important roles as expected. In dilute solutions weak shear thinning and the migration of polymer molecules in the inhomogeneous flow field also appear to affect the physics of the problem.