Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Unsteady flow, clusters, and bands in a model shear-thickening fluid

Shibu Saw, Matthias Grob, Annette Zippelius, and Claus Heussinger
Phys. Rev. E 101, 012602 – Published 13 January 2020
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • MODEL
  • RESULTS AND DISCUSSION
  • CONCLUSION
  • ACKNOWLEDGMENTS
    • Supplemental Material
    References

    Abstract

    We analyze the flow curves of a two-dimensional assembly of granular particles which are interacting via frictional contact forces. For packing fractions slightly below jamming, the fluid undergoes a large scale instability, implying a range of stress and strain rates where no stationary flow can exist. Whereas small systems were shown previously to exhibit hysteretic jumps between the low and high stress branches, large systems exhibit continuous shear thickening arising from averaging unsteady, spatially heterogeneous flows. The observed large scale patterns as well as their dynamics are found to depend on strain rate: At the lower end of the unstable region, force chains merge to form giant bands that span the system in the compressional direction and propagate in the dilational direction. At the upper end, we observe large scale clusters which extend along the dilational direction and propagate along the compressional direction. Both patterns, bands and clusters, come in with infinite correlation length similar to the sudden onset of system-spanning plugs in impact experiments.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    2 More
    • Received 15 May 2019
    • Revised 23 October 2019

    DOI:https://doi.org/10.1103/PhysRevE.101.012602

    ©2020 American Physical Society

    Physics Subject Headings (PhySH)

    1. Research Areas
    JammingNon-Newtonian fluidsShear flows
    1. Physical Systems
    Flow instabilityGranular flows
    1. Properties
    Shear thickening
    1. Techniques
    Molecular dynamics
    Fluid DynamicsPolymers & Soft Matter

    Authors & Affiliations

    Shibu Saw, Matthias Grob, Annette Zippelius, and Claus Heussinger

    • Institute for Theoretical Physics, Georg-August University of Göttingen, Friedrich-Hund Platz 1, 37077 Göttingen, Germany

    Article Text (Subscription Required)

    Click to Expand

    Supplemental Material (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand
    Issue

    Vol. 101, Iss. 1 — January 2020

    Reuse & Permissions
    Access Options
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      Flow curves for different system sizes N and volume fractions ϕ: ϕ=0.78,N=3600 (black), ϕ=0.7975,N=8000 (blue), and N=32000 (green). Inset: different starting configurations (ϕ=0.801,N=80000); the yield-stress branch is metastable.

      Reuse & Permissions
    • Figure 2
      Figure 2

      Left: stress σ as a function of time for N=8000, ϕ=0.7975, and γ̇=1.143×10−4; right: corresponding distribution of stress P(σ). The average over this distribution, restricted to the regions around the two peaks, gives the two data points in the blue curve of Fig. 1 that carry an error bar. The error is obtained by varying the stress value at which the distribution is split.

      Reuse & Permissions
    • Figure 3
      Figure 3

      Snapshots of local shear stress for (a) γ̇=3.297×10−5, (b) γ̇=1.143×10−4, (c) γ̇=10−3, and (d) γ̇=0.05. Particle i is colored in black (gray) whenever its local shear stress σi is above (below) average. (e) Snapshot of local shear stress for a high aspect ratio sample: Lx=10,Ly=260 (γ̇=1.732×10−4,ϕ=0.79875).

      Reuse & Permissions
    • Figure 4
      Figure 4

      Correlation function vs distance r taken along (a) the compressive and (b) dilational directions. Data are averaged over the corresponding perpendicular direction. N=80000, ϕ=0.79750.

      Reuse & Permissions
    • Figure 5
      Figure 5

      Correlation length ξ vs strain rate γ̇. Length scale extracted from the decay of the stress correlation function Cσ(r) with r taken along (a) the compressive and (b) the tensile directions.

      Reuse & Permissions
    • Figure 6
      Figure 6

      Profiles of local shear stress σ and number of contacts per particle z, taken along the (tensile) diagonal direction r. Values are averaged over the compressive direction. N=80000, ϕ=0.79875.

      Reuse & Permissions
    • Figure 7
      Figure 7

      Average contact lifetime 〈τ〉 vs strain rate γ̇. Highlighted are the strain rate γ̇l at the onset and at the upper end γ̇u of shear thickening. N=80000, ϕ=0.79750.

      Reuse & Permissions
    • Figure 8
      Figure 8

      Global pressure p (black) and global stress σ (red) as a function of strain γ. N=80000, ϕ=0.79875.

      Reuse & Permissions
    • Figure 9
      Figure 9

      (a) Time evolution of stress bands at the lower end of ST (γ̇=1.1×10−4). The x axis represents the spatial coordinate along the dilational direction, i.e., along the propagation direction of the band; values are averaged over perpendicular stripes. The y axis is time. The color code is shear stress. Four system-spanning bands move along the dilational direction but in opposite directions, collide, and move on. (b) The same event, now with particle accelerations. (c) Stress at the upper end of ST (γ̇=1.0×10−3). The x axis now represents the spatial coordinate along the compressional direction; values are averaged over perpendicular stripes.

      Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review E

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×