Hostname: page-component-7c8c6479df-94d59 Total loading time: 0 Render date: 2024-03-26T19:50:19.122Z Has data issue: false hasContentIssue false

The determination of the bulk stress in a suspension of spherical particles to order c2

Published online by Cambridge University Press:  29 March 2006

G. K. Batchelor
Affiliation:
Department of Applied Mathematics and Theoretical Physics, University of Cambridge
J. T. Green
Affiliation:
Department of Applied Mathematics and Theoretical Physics, University of Cambridge

Abstract

An exact formula is obtained for the term of order c2 in the expression for the bulk stress in a suspension of force-free spherical particles in Newtonian ambient fluid, where c is the volume fraction of the spheres and c [Lt ] 1. The particles may be of different sizes, and composed of either solid or fluid of arbitrary viscosity. The method of derivation circumvents the familiar obstacle, of non-absolutely convergent integrals representing the effect of all pair interactions in which one specified particle takes part, by the judicious use of a certain quantity which is affected by the presence of distant particles in a similar way and whose mean value is known exactly. The bulk stress is in general of non-Newtonian form and depends on the statistical properties of the suspension which in turn are dependent on the type of bulk flow.

The formula contains two functions which are parameters of the flow field due to two spherical particles immersed in fluid in which the velocity gradient is uniform at infinity. One of them, p(r, t), represents the probability density for the vector r separating the centres of the two particles. The variation of p(r, t) for a moving material point in r-space due to hydrodynamic action is found in terms of a function q(r), and this gives p(r, t) explicitly over the whole of the region of r-space occupied by trajectories of one particle centre relative to another which come from infinity. In a region of closed trajectories, steady-state hydrodynamic action alone does not determine the relation between the values of p (r, t) for different material points. The function q(r) is singular when the spheres touch, and the contribution of nearly-touching spheres to the bulk stress is evidently important. Approximate numerical values of all the relevant functions are presented for the case of rigid spherical particles of uniform size.

In the case of steady pure straining motion of the suspension, all trajectories in r-space come from infinity, the suspension has isotropic structure and the stress behaviour can be represented (to order c2) in terms of an effective viscosity ${\mathop\mu\limits^{*}}$. It is estimated from the available numerical data that for a suspension of identical rigid spherical particles \[ {\mathop\mu\limits^{*}}/\mu = 1 + 2.5c + 7.6c^2, \] the error bounds on the coefficient of c2 being about ∓ 0.8. In the important case of steady simple shearing motion, there is a region of closed trajectories of one sphere centre relative to another, of infinite volume. The stress system is here not of Newtonian form, and numerical results are not obtainable until the probability, density p(r, t) can be made determinate in the region of closed trajectories by the introduction of some additional physical process, such as three-sphere encounters or Brownian motion, or by the assumption of some particular initial state.

In the analogous problem for an incompressible solid suspension it may be appropriate to assume that for many methods of manufacture p(r, t) is uniform over the accessible part of r-space, in which event the solid suspension has ‘Newtonian’ elastic behaviour and the ratio of the effective shear modulus to that of the matrix is estimated to be 1 + 2·5c + 5·2c2 for a suspension of identical rigid spheres.

Type
Research Article
Copyright
© 1972 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Batchelor, G. K. 1970 The stress system in a suspension of force-free particles. J. Fluid Mech. 41, 545.Google Scholar
Batchelor, G. K 1972 Sedimentation in a dilute suspension of spheres J. Fluid Mech. 52, 245.Google Scholar
Batchelor, G. K. & Green, J. T. 1972 The hydrodynamic interaction of two small freely-moving spheres in a linear flow field. J. Fluid Mech. 56, 375.Google Scholar
Cox, R. G. & Brenner, H. 1971 The rheology of a suspension of particles in a Newtonian fluid. Chem. Engng Sci. 26, 65.Google Scholar
Darabaner, C. L. & Mason, S. G. 1967 Particle motions in sheared suspensions. Xxii: Interactions of rigid spheres. (Experimental.) Rheol. Acta 6, 273.CrossRefGoogle Scholar
Einstein, A. 1906 Eine neue Bestimmung der Moleküldimensionen. Annln. Phys. 19, 289 (and 34, 591).Google Scholar
Hashin, Z. 1964a Bounds for viscosity coefficients of fluid mixtures by variational methods. Proceedings of I.U.T.A.M. Symposium on ‘Second-order effects in Elasticity, Plasticity and Fluid Dynamics’, Haija, 1962, edited by M. Reiner and D. Abir. Pergamon Press.
Hashin, Z. 1964b Theory of mechanical behaviour of heterogeneous media. Appl. Mech. Rev. 17, 1.Google Scholar
Keller, J. B., Rubenfeld, L. A. & Molyneux, J. E. 1967 Extremum principles for slow viscous flows with applications to suspensions. J. Fluid Mech. 30, 97.Google Scholar
LIN, C. J., LEE, K. J. & Satecer, N. F. 1970 Slow motion of two spheres in a shear field. J. Fluid Mech. 43, 35.Google Scholar
Peterson, J. M. & Fixman, M. 1963 Viscosity of polymer solutions. J. Chern. Phys. 39, 2516.Google Scholar
Rutgers, R. 1962a Relative viscosity of suspensions of rigid spheres in Newtonian liquids. Rheol. Acta 2, 202.Google Scholar
Rutuers, R. 1962b Relative viscosity and concentration. Rheol. Acta. 2, 305.Google Scholar
Taylor, G. I. 1932 The viscosity of a fluid containing small drops of another fluid. Proc. Roy. Soc. A 138, 41.Google Scholar
Walpole, L. J. 1972 The elastic behaviour of a suspension of spherical particles. Quart. J. Mech. Appl. Math. 25, 153.Google Scholar