scholarly journals Effect of confinement in wall-bounded non-colloidal suspensions

2016 ◽  
Vol 799 ◽  
pp. 100-127 ◽  
Author(s):  
Stany Gallier ◽  
Elisabeth Lemaire ◽  
Laurent Lobry ◽  
Francois Peters

This paper presents three-dimensional numerical simulations of non-colloidal dense suspensions in a wall-bounded shear flow at zero Reynolds number. Simulations rely on a fictitious domain method with a detailed modelling of particle–particle and wall–particle lubrication forces, as well as contact forces including particle roughness and friction. This study emphasizes the effect of walls on the structure, velocity and rheology of a moderately confined suspension (channel gap to particle radius ratio of 20) for a volume fraction range $0.1\leqslant {\it\phi}\leqslant 0.5$. The wall region shows particle layers with a hexagonal structure. The size of this layered zone depends on volume fraction and is only weakly affected by friction. This structure implies a wall slip which is in good accordance with empirical models. Simulations show that this wall slip can be mitigated by reducing particle roughness. For ${\it\phi}\lessapprox 0.4$, wall-induced layering has a moderate impact on the viscosity and second normal stress difference $N_{2}$. Conversely, it significantly alters the first normal stress difference $N_{1}$ and can result in positive $N_{1}$, in better agreement with some experiments. Friction enhances this effect, which is shown to be due to a substantial decrease in the contact normal stress $|{\it\Sigma}_{xx}^{c}|$ (where $x$ is the velocity direction) because of particle layering in the wall region.

2011 ◽  
Vol 686 ◽  
pp. 26-39 ◽  
Author(s):  
Étienne Couturier ◽  
François Boyer ◽  
Olivier Pouliquen ◽  
Élisabeth Guazzelli

AbstractWe measure the second normal-stress difference in suspensions of non-Brownian neutrally buoyant rigid spheres dispersed in a Newtonian fluid. We use a method inspired by Wineman & Pipkin (Acta Mechanica, vol. 2, 1966, pp. 104–115) and Tanner (Trans. Soc. Rheol., vol. 14, 1970, pp. 483–507), which relies on the examination of the shape of the suspension free surface in a tilted trough flow. The second normal-stress difference is found to be negative and linear in shear stress. The ratio of the second normal-stress difference to shear stress increases with increasing volume fraction. A clear behavioural change exhibiting a strong (approximately linear) growth in the magnitude of this ratio with volume fraction is seen above a volume fraction of 0.22. By comparing our results with previous data obtained for the same batch of spheres by Boyer, Pouliquen & Guazzeli (J. Fluid Mech., 2011, doi:10.1017/jfm.2011.272), the ratio of the first normal-stress difference to the shear stress is estimated and its magnitude is found to be very small.


2014 ◽  
Vol 758 ◽  
pp. 486-507 ◽  
Author(s):  
Braden Snook ◽  
Levi M. Davidson ◽  
Jason E. Butler ◽  
Olivier Pouliquen ◽  
Élisabeth Guazzelli

AbstractMeasurements of normal stress differences are reported for suspensions of rigid, non-Brownian fibres for concentrations of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}nL^2d=1.5\text {--}3$ and aspect ratios of $L/d=11\text {--}32$, where $n$ is the number of fibres per unit volume, $L$ is the fibre length and $d$ is the diameter. The first and second normal stress differences are determined experimentally from measuring the deformation in the free surface in a tilted trough and in a Weissenberg rheometer. Simulations are performed as well, and the hydrodynamic and contact contributions to the normal stresses are calculated. The experiments and simulations indicate that the second normal stress difference is negative and that its magnitude increases as the concentration is raised and the aspect ratio is lowered. The first normal stress difference is positive and its magnitude is approximately twice that of the second normal stress difference. Simulation results indicate that, for the concentrations and aspect ratios studied, contact forces between fibres form the dominant contribution to the normal stress differences.


2015 ◽  
Vol 781 ◽  
pp. 506-549 ◽  
Author(s):  
Mohsen Daghooghi ◽  
Iman Borazjani

We investigate the rheological properties of a suspension of neutrally buoyant rigid ellipsoids by fluid–structure interaction simulations of a particle in a periodic domain under simple shear using the curvilinear immersed-boundary (CURVIB) method along with a quaternion–angular velocity technique to calculate the dynamics of the particle’s motion. We calculate all the different terms of particle stress for the first time for non-spherical particles, i.e. in addition to the stresslet, we calculate the acceleration and Reynolds stress, which are typically ignored in previous similar works. Furthermore, we derive analytical expressions for all these terms to verify the numerical results and deduce the effect of inertia by comparing our numerical results with the analytical solution. The effect of particle Reynolds number ($\mathit{Re}$), volume fraction (${\it\phi}$), and the shape of particles has been studied on all mechanisms of stress generation, the intrinsic viscosity, and normal stress differences of the suspension for the range$0.008\leqslant {\it\phi}\leqslant 0.112$and$0.01\leqslant \mathit{Re}\leqslant 10.0$. We found that inertia increases the shear and the second normal difference of the stresslet (dominant term of the particle stress), and decreases the first normal difference that is generated due to the strain field. The contribution of acceleration stress to the total stress is found to be important in the second normal stress difference, with a cycle-average comparable to the stresslet component. We also discovered that the contribution of Reynolds stress in the first normal stress difference becomes important even when inertia is as low as$\mathit{Re}\sim O(0.1)$, and its value can be even greater than the stresslet when inertia increases, i.e. Reynolds stresses cannot be ignored for non-spherical particles. For concentrations in the range from dilute to semi-dilute, the effect of inertia on the intrinsic viscosity of a suspension is found to be comparable to the volume fraction. Furthermore, our calculations show that for a dilute concentration and the low-inertia regime ($\mathit{Re}<1.0$), the intrinsic viscosity of a suspension consisting of ellipsoids with an aspect ratio of five can be 20 % higher than its Stokesian analytical value.


2000 ◽  
Vol 407 ◽  
pp. 167-200 ◽  
Author(s):  
DAVID R. FOSS ◽  
JOHN F. BRADY

The non-equilibrium behaviour of concentrated colloidal dispersions is studied using Stokesian Dynamics, a molecular-dynamics-like simulation technique for analysing suspensions of particles immersed in a Newtonian fluid. The simulations are of a monodisperse suspension of Brownian hard spheres in simple shear flow as a function of the Péclet number, Pe, which measures the relative importance of hydrodynamic and Brownian forces, over a range of volume fraction 0.316 [les ] ϕ [les ] 0.49. For Pe < 10, Brownian motion dominates the behaviour, the suspension remains well-dispersed, and the viscosity shear thins. The first normal stress difference is positive and the second negative. At higher Pe, hydrodynamics dominate resulting in an increase in the long-time self-diffusivity and the viscosity. The first normal stress difference changes sign when hydrodynamics dominate. Simulation results are shown to agree well with both theory and experiment.


2015 ◽  
Vol 786 ◽  
pp. 110-127 ◽  
Author(s):  
D. Matsunaga ◽  
Y. Imai ◽  
T. Yamaguchi ◽  
T. Ishikawa

We present a numerical analysis of the rheology of a dense suspension of spherical capsules in simple shear flow in the Stokes flow regime. The behaviour of neo-Hookean capsules is simulated for a volume fraction up to${\it\phi}=0.4$by graphics processing unit computing based on the boundary element method with a multipole expansion. To describe the specific viscosity using a polynomial equation of the volume fraction, the coefficients of the equation are calculated by least-squares fitting. The results suggest that the effect of higher-order terms is much smaller for capsule suspensions than rigid sphere suspensions; for example,$O({\it\phi}^{3})$terms account for only 8 % of the specific viscosity even at${\it\phi}=0.4$for capillary numbers$Ca\geqslant 0.1$. We also investigate the relationship between the deformation and orientation of the capsules and the suspension rheology. When the volume fraction increases, the deformation of the capsules increases while the orientation angle of the capsules with respect to the flow direction decreases. Therefore, both the specific viscosity and the normal stress difference increase with volume fraction due to the increased deformation, whereas the decreased orientation angle suppresses the specific viscosity, but amplifies the normal stress difference.


This paper describes a series of experiments in which the three material functions of steady viscometric flows were measured for a given polyisobutene solution. A number of instruments and measuring techniques were used in order to check the experimental method. The shear stress was determined from the torque transmitted by the fluid in a cone-and-plate apparatus and in Couette flow between concentric cylinders. The results obtained from these measurements were in good agreement with each other. The primary normal-stress difference was determined from the normal force acting on the plate of a cone-and-plate apparatus, and from stress-optical measurements on Couette flow between concentric cylinders. These results are in good agreement with each other. Detailed measurements of the distribution f Permanent address: Fluid Mechanics Research Institute, University of Essex, Colchester, Essex. of the normal stress acting on the plate of the cone-and-plate apparatus were made for three cone angles and for two boundary configurations at the rim of the apparatus: from these results a combination of the primary and the secondary normal-stress differences was deduced, thereby making possible the computation of the secondary normal-stress difference. When the normal stress acting on a rigid surface is measured by means of a hole leading to a pressure transducer the results are in error by an amount roughly proportional to the primary normal-stress difference of the fluid (cf. Kaye, Lodge & Vale 1968). In the present experiments this error was determined from measurements of the distribution of the normal stress acting on the plates of a plate-and-plate apparatus, together with the assumption that the error is a function only of the shear rate at the position o the hole in the undisturbed viscometric flow. The values of the measuring error thus obtained are in goo agreement with measurements made in Gouette flow between concentric cylinders. The secondary normal-stress difference, P2, was measured in a number of different ways. From the results it is suggested that the methods of Jackson & Kaye and of Marsh & Pearson may be imprecise and, in particular, may yield incorrect values for P2- A new, direct, method of estimating P2, suggested by Higashitani & Pritchard (1971) and outlined in appendix A, may provide a more convenient means of determining P2.


Sign in / Sign up

Export Citation Format

Share Document