Strong and weak blow-up of the viscous dissipation rates for concentrated suspensions

We study the overall dissipation rate of highly concentrated non-colloidal suspensions of rigid neutrally buoyant particles in a Newtonian fluid. This suspension is confined to a finite size container, subject to shear or extensional boundary conditions at the walls of the container. The corresponding dissipation rates determine the effective shear viscosity μ* and the extensional effective viscosity λ*. We use recently developed discrete network approximation techniques to obtain discrete forms for the overall dissipation rates, and analyse their asymptotics in the limit when the characteristic interparticle distance goes to zero. The focus is on the finite size and particle wall effects in spatially disordered arrays. Use of the network approximation allows us to study the dependence of μ* and λ* on variable distances between neighbouring particles in such arrays. Our analysis, carried out for a two-dimensional model, can be characterized as global because it goes beyond the local analysis of flow in a single gap between two particles and takes into account hydrodynamic interactions in the entire particle array. The principal conclusion in the paper is that, in general, asymptotic formulae for μ* and λ* obtained by global analysis are different from the formulae obtained from local analysis. In particular, we show that the leading term in the asymptotics of μ* is of lower order than suggested by the local analysis (weak blow-up), while the order of the leading term in the asymptotics of λ* depends on the geometry of the particle array (either weak or strong blow-up). We obtain geometric conditions on a random particle array under which the asymptotic order of λ* coincides with the order of the local dissipation in a gap between two neighbouring particles, and show that these conditions are generic. We also provide an example of a uniformly closely packed particle array for which the leading term in the asymptotics of λ* degenerates (weak blow-up).

[1]  de Kruif CG,et al.  Linear viscoelastic behavior of dense hard-sphere dispersions. , 1989, Physical review. A, General physics.

[2]  John F. Brady,et al.  Rheology and microstructure in concentrated noncolloidal suspensions , 2002 .

[3]  Joseph B. Keller,et al.  Effective viscosity of a periodic suspension , 1984, Journal of Fluid Mechanics.

[4]  Mihail C. Roco,et al.  Slurry flow : principles and practice , 1991 .

[5]  N. Bakhvalov,et al.  Homogenisation: Averaging Processes in Periodic Media , 1989 .

[6]  G. Batchelor,et al.  The determination of the bulk stress in a suspension of spherical particles to order c2 , 1972, Journal of Fluid Mechanics.

[7]  E. S. Palencia Non-Homogeneous Media and Vibration Theory , 1980 .

[8]  J. Stickel,et al.  FLUID MECHANICS AND RHEOLOGY OF DENSE SUSPENSIONS , 2001 .

[9]  John F. Brady,et al.  Microstructure of strongly sheared suspensions and its impact on rheology and diffusion , 1997, Journal of Fluid Mechanics.

[10]  Sangtae Kim,et al.  Microhydrodynamics: Principles and Selected Applications , 1991 .

[11]  John F. Brady,et al.  Accelerated Stokesian Dynamics simulations , 2001, Journal of Fluid Mechanics.

[12]  M. Denn,et al.  On the viscosity of a concentrated suspension of solid spheres , 1985 .

[13]  D. S. Pearson,et al.  Viscoelastic behavior of concentrated spherical suspensions , 1994 .

[14]  A. Einstein Eine neue Bestimmung der Moleküldimensionen , 1905 .

[15]  J. Lumley,et al.  Mechanics of non-Newtonian fluids , 1978 .

[16]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .

[17]  L. Berlyand,et al.  Fictitious Fluid Approach and Anomalous Blow-up of the Dissipation Rate in a 2D Model of Concentrated Suspensions , 2006, math/0608671.

[18]  Leonid Berlyand,et al.  Network Approximation in the Limit of¶Small Interparticle Distance of the¶Effective Properties of a High-Contrast¶Random Dispersed Composite , 2001 .

[19]  R. Chhabra,et al.  Transport Processes In Bubbles, Drops And Particles , 1991 .

[20]  Alexander Panchenko,et al.  Network Approximation for Effective Viscosity of Concentrated Suspensions with Complex Geometry , 2005, SIAM J. Math. Anal..

[21]  Herbert Edelsbrunner,et al.  Triangulations and meshes in computational geometry , 2000, Acta Numerica.

[22]  Robin Ball,et al.  A simulation technique for many spheres in quasi-static motion under frame-invariant pair drag and Brownian forces , 1997 .

[23]  L. Silbert,et al.  Colloidal microdynamics : pair-drag simulations of model-concentrated aggregated systems , 1997 .

[24]  V. Zhikov,et al.  Homogenization of Differential Operators and Integral Functionals , 1994 .