The determination of the bulk stress in a suspension of spherical particles to order c 2
TL;DR: In this article, an exact formula 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.8.
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 . 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.
Citations
More filters
TL;DR: In this paper, a simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure.
Abstract: Recent advances in theory and experimen- tation motivate a thorough reassessment of the physics of debris flows. Analyses of flows of dry, granular solids and solid-fluid mixtures provide a foundation for a com- prehensive debris flow theory, and experiments provide data that reveal the strengths and limitations of theoret- ical models. Both debris flow materials and dry granular materials can sustain shear stresses while remaining stat- ic; both can deform in a slow, tranquil mode character- ized by enduring, frictional grain contacts; and both can flow in a more rapid, agitated mode characterized by brief, inelastic grain collisions. In debris flows, however, pore fluid that is highly viscous and nearly incompress- ible, composed of water with suspended silt and clay, can strongly mediate intergranular friction and collisions. Grain friction, grain collisions, and viscous fluid flow may transfer significant momentum simultaneously. Both the vibrational kinetic energy of solid grains (mea- sured by a quantity termed the granular temperature) and the pressure of the intervening pore fluid facilitate motion of grains past one another, thereby enhancing debris flow mobility. Granular temperature arises from conversion of flow translational energy to grain vibra- tional energy, a process that depends on shear rates, grain properties, boundary conditions, and the ambient fluid viscosity and pressure. Pore fluid pressures that exceed static equilibrium pressures result from local or global debris contraction. Like larger, natural debris flows, experimental debris flows of ;10 m 3 of poorly sorted, water-saturated sediment invariably move as an unsteady surge or series of surges. Measurements at the base of experimental flows show that coarse-grained surge fronts have little or no pore fluid pressure. In contrast, finer-grained, thoroughly saturated debris be- hind surge fronts is nearly liquefied by high pore pres- sure, which persists owing to the great compressibility and moderate permeability of the debris. Realistic mod- els of debris flows therefore require equations that sim- ulate inertial motion of surges in which high-resistance fronts dominated by solid forces impede the motion of low-resistance tails more strongly influenced by fluid forces. Furthermore, because debris flows characteristi- cally originate as nearly rigid sediment masses, trans- form at least partly to liquefied flows, and then trans- form again to nearly rigid deposits, acceptable models must simulate an evolution of material behavior without invoking preternatural changes in material properties. A simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure. These equations can describe a spectrum of debris flow behav- iors intermediate between those of wet rock avalanches and sediment-laden water floods. With appropriate pore pressure distributions the equations yield numerical so- lutions that successfully predict unsteady, nonuniform motion of experimental debris flows.
2,426 citations
Cites background from "The determination of the bulk stres..."
...…1967] have deduced equations to predict the effective viscosity of concentrated suspensions of fine spheres, other investigations have explained why no such equation can be expected to work well for the full range of yfines and all conceivable flow fields [Batchelor and Green, 1972; Acrivos, 1993]....
[...]
TL;DR: In this paper, the authors review the analysis of composite materials from the applied mechanics and engineering science point of view, including elasticity, thermal expansion, moisture swelling, viscoelasticity, conductivity, static strength, and fatigue failure.
Abstract: The purpose of the present survey is to review the analysis of composite materials from the applied mechanics and engineering science point of view. The subjects under consideration will be analysis of the following properties of various kinds of composite materials: elasticity, thermal expansion, moisture swelling, viscoelasticity, conductivity (which includes, by mathematical analogy, dielectrics, magnetics, and diffusion) static strength, and fatigue failure.
2,216 citations
TL;DR: In this article, the effect of Brownian motion on the probability density of the separation vector of rigid spherical particles in a dilute suspension is investigated and an explicit expression for this leading approximation is constructed in terms of hydrodynamic interactions between pairs of particles.
Abstract: The effect of Brownian motion of particles in a statistically homogeneous suspension is to tend to make uniform the joint probability density functions for the relative positions of particles, in opposition to the tendency of a deforming motion of the suspension to make some particle configurations more common. This smoothing process of Brownian motion can be represented by the action of coupled or interactive steady ‘thermodynamic’ forces on the particles, which have two effects relevant to the bulk stress in the suspension. Firstly, the system of thermodynamic forces on particles makes a direct contribution to the bulk stress; and, secondly, thermodynamic forces change the statistical properties of the relative positions of particles and so affect the bulk stress indirectly. These two effects are analysed for a suspension of rigid spherical particles. In the case of a dilute suspension both the direct and indirect contributions to the bulk stress due to Brownian motion are of order o2, where o([Lt ] 1) is the volume fraction of the particles, and an explicit expression for this leading approximation is constructed in terms of hydrodynamic interactions between pairs of particles. The differential equation representing the effects of the bulk deforming motion and the Brownian motion on the probability density of the separation vector of particle pairs in a dilute suspension is also investigated, and is solved numerically for the case of relatively strong Brownian motion. The suspension has approximately isotropic structure in this case, regardless of the nature of the bulk flow, and the effective viscosity representing the stress system to order ϕ2 is found to be
\[
\mu^{*} = \mu(1+2.5\phi + 6.2\phi^2).
\]
The value of the coefficient of o2 for steady pure straining motion in the case of weak Brownian motion is known to be 7[sdot ]6, which indicates a small degree of ‘strain thickening’ in the o2-term.
1,956 citations
TL;DR: This work reviews many significant developments over the past decade of the lattice-Boltzmann method and discusses higherorder boundary conditions and the simulation of microchannel flow with finite Knudsen number.
Abstract: With its roots in kinetic theory and the cellular automaton concept, the lattice-Boltzmann (LB) equation can be used to obtain continuum flow quantities from simple and local update rules based on particle interactions. The simplicity of formulation and its versatility explain the rapid expansion of the LB method to applications in complex and multiscale flows. We review many significant developments over the past decade with specific examples. Some of the most active developments include the entropic LB method and the application of the LB method to turbulent flow, multiphase flow, and deformable particle and fiber suspensions. Hybrid methods based on the combination of the Eulerian lattice with a Lagrangian grid system for the simulation of moving deformable boundaries show promise for more efficient applications to a broader class of problems. We also discuss higherorder boundary conditions and the simulation of microchannel flow with finite Knudsen number. Additionally, the remarkable scalability of the LB method for parallel processing is shown with examples. Teraflop simulations with the LB method are routine, and there is no doubt that this method will be one of the first candidates for petaflop computational fluid dynamics in the near future.
1,585 citations
01 Jan 1974
1,090 citations
References
More filters
TL;DR: In this article, it was shown that if the radius of the suspended drops or the velocity of distortion of the fluid are small, surface tension may be expected to keep them nearly spherical, and in that case Einstein's analysis may be extended so as to include the case of liquid drops.
Abstract: The viscosity of a fluid in which small solid spheres are suspended has been studied by Einstein as a problem in theoretical hydrodynamics. Einstein’s paper gave rise to many experimental researches on the viscosity of fluids containing solid particles, and it soon became clear that though complete agreement with the theory might be expected when the particles are true sphered, some modification is necessary when the particles are flattened or elongated. The theory of such systems was developed by G. B. Jeffery, who calculated the motion of ellipsoidal particles in a viscous fluid and their effect on the mean viscosity. Some of his conclusions have been verified by observation. So far no one seems to have extended Einstein’s work to liquids containing small drops of another liquid in suspension. The difficulties in the way of a complete theory when solid particles are replaced by fluid drops are almost insuperable, partly because the correct boundary conditions are not known, and partly because a fluid drop would deform under the combined action of viscous forces and surface tension. Even if the boundary conditions were known to be those commonly used in hydrodynamical theory, the calculation of the shape of the deformed drop would be exceedingly difficult. When the radius of the suspended drops or the velocity of distortion of the fluid are small, surface tension may be expected to keep them nearly spherical, and in that case Einstein’s analysis may be extended so as to include the case of liquid drops.
1,708 citations
TL;DR: In this paper, the authors consider the properties of the bulk stress in a suspension of non-spherical particles, on which a couple (but no force) may be imposed by external means, immersed in a Newtonian fluid.
Abstract: The purpose of the paper is to consider in general terms the properties of the bulk stress in a suspension of non-spherical particles, on which a couple (but no force) may be imposed by external means, immersed in a Newtonian fluid. The stress is sought in terms of the instantaneous particle orientations, and the problem of determining these orientations from the history of the motion is not considered. The bulk stress and bulk velocity gradient in the suspension are defined as averages over an ensemble of realizations, these averages being equal to integrals over a suitably chosen volume of ambient fluid and particles together when the suspension is statistically homogeneous. Without restriction on the type of particle or the concentration or the Reynolds number of the motion, the contribution to the bulk stress due to the presence of the particles is expressed in terms of integrals involving the stress and velocity over the surfaces of particles together with volume integrals not involving the stress. The antisymmetric part of this bulk stress is equal to half the total couple imposed on the particles per unit volume of the suspension. When the Reynolds number of the relative motion near one particle is small, a suspension of couple-free particles of constant shape is quasi-Newtonian; i.e. the dependence of the bulk stress on bulk velocity gradient is linear. Two significant features of a suspension of non-spherical particles are (1) that this linear relation is not of the Newtonian form and (2) that the effect of exerting a couple on the particles is not confined to the generation of an antisymmetrical part of the bulk stress tensor. The role of surface tension at the particle boundaries is described.In the case of a dilute suspension the contributions to the bulk stress from the various particles are independent, and the contributions arising from the bulk rate of strain and from the imposed couple are independent for each particle. Each particle acts effectively as a force doublet (i.e. equal and opposite adjoining ‘Stokeslets’) whose tensor strength determines the disturbance flow far from the particle and whose symmetrical and antisymmetrical parts are designated as a stresslet and a couplet. The couplet strength is determined wholly by the externally imposed couple on the particle; but the stresslet strength depends both on the bulk rate of strain and, for a non-spherical particle, on the rate of rotation of the particle relative to the fluid resulting from the imposed couple. The general properties of the stress system in a dilute suspension are illustrated by the specific and complete results which may be obtained for rigid ellipsoidal particles by use of the work by Jeffery (1922).
1,428 citations
TL;DR: In this article, the authors considered a large number of identical small rigid spheres with random positions which are falling through Newtonian fluid under gravity and determined the mean value of the velocity of a sphere (U).
Abstract: The dispersion considered consists of a large number of identical small rigid spheres with random positions which are falling through Newtonian fluid under gravity. The volume fraction of the spheres (c) is small compared with unity. The dispersion is statistically homogeneous, and the axes of reference are chosen so that the mean volume flux across any stationary surface is zero. The problem is to determine the mean value of the velocity of a sphere (U). In §3 there is described a systematic and rigorous procedure which overcomes the familiar difficulty presented by the occurrence of divergent integrals, essentially by the choice of a quantity V whose mean value can be found exactly and which has the same long-range dependence on the position of a second sphere as U so that the mean of U – V can be expressed in terms of an absolutely convergent integral. The result is that, correct to order c, the mean value of U is U0(1 – 6.55 c), where U0, is the velocity of a single sphere in unbounded fluid. The only assumption made in the calculation is that the centres of spheres in the dispersion take with equal probability all positions such that no two spheres overlap; arguments are given in support of this assumption, which is expected to be valid only when the spheres are identical. Calculations which assume a simple regular arrangement of the spheres or which adopt a cell model of the hydrodynamic interactions give the quite different result that the change in the mean speed of fall is proportional to , for reasons which are made clear.The general procedure described here is expected to be applicable to other problems concerned with the effect of particle interactions on the average properties of dispersions with small volume fraction of the particles.
1,158 citations
TL;DR: In this article, the authors provide a systematic and explicit description of the interaction between two rigid spheres that are relevant in a calculation of the mean stress in a suspension of spherical particles subjected to bulk deformation.
Abstract: Two rigid spheres of radii a and b are immersed in infinite fluid whose velocity at infinity is a linear function of position. No external force or couple acts on the spheres, and the effect of inertia forces on the motion of the fluid and the spheres is neglected. The purpose of the paper is to provide a systematic and explicit description of those aspects of the interaction between the two spheres that are relevant in a calculation of the mean stress in a suspension of spherical particles subjected to bulk deformation. The most relevant aspects are the relative velocity of the two sphere centres (V) and the force dipole strengths of the two spheres (S′ij, S″ij), as functions of the vector r separating the two centres.It is shown that V, S′ij and S″ij depend linearly on the rate of strain at infinity and can be represented in terms of several scalar parameters which are functions of r/a and b/a alone. These scalar functions provide a framework for the expression of the many results previously obtained for particular linear ambient flows or for particular values of r/a or of b/a. Some new results are established for the asymptotic forms of the functions both for r/(a + b) [Gt ] 1 and for values of r − (a + b) small compared with a and b. A reasonably complete numerical description of the interaction of two rigid spheres of equal size is assembled, the main deficiency being accurate values of the scalar functions describing the force dipole strength of a sphere in the intermediate range of sphere separations.In the case of steady simple shearing motion at infinity, some of the trajectories of one sphere centre relative to another are closed, a fact which has consequences for the rheological problem. These closed forms are described analytically, and also numerically in the case b/a = 1.
697 citations