Rheology of concentrated hard-sphere suspensions Complex behaviour in a simple system Wilson Poon School of Physics & Astronomy, The University of Edinburgh Rut Besseling Lucio Isa ( ETH, Zürich) Khoa Pham ( CSIRO) Pierre Ballesta ( Crete) Jon Singer (??) Peter Pusey Dimitris Vlassopoullos (Crete) George Petekedis (Crete) Eric Weeks (Atlanta)
What am I doing here?
New tool: single-particle, real time imaging of concentrated suspensions under flow used to study a zoo of complex behaviour, much of which is yet to be theoretically explained.
Experimental hard sphere system PMMA 1µm 10nm Stabilised by chemically-grafted polymer brush Rheology textbook: Shear-thinning fluid Yield-stress fluid fluid ~ 12% polydispersity or more f + c crystal 0.494 0.545 ~ 0.58 Glass φ
Popular constitutive equations for yield-stress suspensions Bingham, E. C. Fluidity and Plasticity. New York: McGraw-Hill, 1922. Herschel, W.H.; Bulkley, R. (1926), Konsistenzmessungen von Gummi- Benzollösungen, Kolloid Zeitschrift 39: 291 300 Casson N. (1959), A flow equation for pigment-oil suspensions of the printing ink type. In: Mill CC, editor. Rheology of disperse systems. London, U.K.: Pergamon.
Chocolate + molten cocoa butter + shear ( conching ) = 50% suspension of sugar in oil
PMMA suspensions Fits H-B Slope = 1 MCT can predict such behaviour ab initio Mike Cates talk
Measuring the constitutive relation is not totally straightforward Cone: rough Plate: Filled = rough Open = smooth No difference when φ < 0.58
Bulk rheology: φ = 0.64 is a H-B yield-stress fluid HB fluid in pipe: plug flow; plug diminishes with flow rate Huilgol & You, JNFM (2005)
Od Flow rate Plug size (fraction of pipe size)
Fluorescent, indexmatched PMMA φ ~ 0.64, 2R ~ 2.6 µm Particle tracking under shear: Decompose image into strips Shift to maximise correlation between successive strips Average v of strips Subtract average motion Track (Crocker-Grier etc.) Constant P 20 d 170 frames/s Add back average motion Fast confocal
1. No dependence of shear zone size on flow speed 2. In any case Od ~ 0.02 Plug size should be 0 Plug size (fraction of pipe size) Od
1. Pouliquen: granular flow [PRE 1996] Yielding due to stress fluctuations ( σ) spreading from walls Coulomb friction dominates Flow speed stress friction yield stress constant shear zone (width controlled by strength of σ) 2. Bocquet [Nature 2008, PRL 2009] ( Dhont): non-local fluidity (STZ) Colin s talk
Confinement-induced flow instabilities ~20D ~40D ~30D
Certainly doesn t tell the whole story for the pipe flow of concentrated hard sphere suspensions how about simpler, rheometric, geometries?
3D confocal imaging in rheometric geometries ~ 5 slices/particle 0.1 s / particle Vertical slice ACIS (2009) Home-built parallel-plate shear cell or AR2000 rheometer cone-plate
Genuine slip next to smooth plate, not shear localization r = 2.5 mm (in 20mm-radius plate)
But shear localization does occur in rough/rough cone/plate 1000 100 σ m (Pa) 10 1 0.1. ϕ 63% 61% 59% 52% H-B bulk rheology?! 1E-3 0.1 10. γ a (s -1 ) 1.0 30 0.5 0.1 0.05 0.01 exponential fit v / v cone 1.0 0.5 v / v cone 0.5 0.0 0.0 0.2 0.4 0.6 0.8 1.0 z / z gap 0.0 0.0 0.2 0.4 0.6 0.8 1.0 z / z gap
Bocquet [Nature 2008, PRL 2009] ( Dhont): non-local fluidity True bulk flow curve Finite size effects?? σ m (Pa) 1000 100 10. 1 0.1 1E-3 0.1 10. γ a (s -1 ) ϕ 63% 61% 59% 52% v / v cone 1.0 0.5 30 0.5 0.1 0.05 0.01 exponential fit 0.0 0.0 0.2 0.4 0.6 0.8 1.0 z / z gap
Parallel-plate shear cell (rough walls) Analyse local regions with linear velocity gradient Measured shear rate > average imposed shear rate
Particle dynamics as a function of local shear rate x, velocity ~ x, z plane z, gradient 15 µm 1 µm 15 µm y,z plane Substract affine component of shear-induced motion ~ x,y plane t = 800 s ~ x, y plane Significant heterogeneity in shear-induced motions Shear-induced cage breakout long-range diffusion in all directions
Structural relaxation Cf. Time-density superposition of quiescent HS colloids at φ < φ g. (van Megen) Time-shear superposition Cf. Theory (Fuchs & Cates; Berthier) F s ( Q m, t ) 1.0. γ = 0. 0.8. 0.6 10 3 γ ( s -1 ) 0.31 0.4 0.45 0.93 0.2 1.48 τ α 2.82 0.0 1 10 100 1000 t (s) 1.0 0.8 f s 0.6 Exponentiality cf. Varnik L-J simulations 0.4 0.2 0.0 0.1 1 10 t / τ α
Scaling with shear rate
0.01., D = 0.9 D 0 (τ B γ ) 0.8 D y (µm 2 /s) 1E-3 1E-4 a=0.85 µm (ϕ~0.62) a=0.43 µm (ϕ~0.6) 1E-4 1E-3 0.01 0.1. γ (s -1 ) 0.1 D y (µm 2 /s) 0.01 1E-3 1E-4 no shear a=0.85 µm φ = 0.62, n=0.8 φ = 0.55, n=0.77 1E-5 1E-5 1E-4 1E-3 0.01 0.1 1. γ (s -1 )
Increasing stress Increasing φ Kobelev PRE 2005 Schweizer 2007
?
Extreme shear thickening Norman Wagner s talk
Conclusions Concentrated suspensions of hard spheres show complex flow behaviour Real-time single-particle resolution imaging helps to reveal (and perhaps elucidate!) this zoo of complexity Method can be extended to other systems That s all, folks!