1
|
Sun Y, Wang C, Yang J, Shi W, Pang Q, Wang Y, Li J, Hu B, Xia C. Evident structural anisotropies arising from near-zero particle asphericity in granular spherocylinder packings. Phys Rev E 2024; 110:014903. [PMID: 39161035 DOI: 10.1103/physreve.110.014903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2024] [Accepted: 07/02/2024] [Indexed: 08/21/2024]
Abstract
With magnetic resonance imaging experiments, we study packings of granular spherocylinders with merely 2% asphericity. Evident structural anisotropies across all length scales are identified. Most interestingly, the global nematic order decreases with increasing packing fraction, while the local contact anisotropy shows an opposing trend. We attribute this counterintuitive phenomenon to a competition between gravity-driven ordering aided by frictional contacts and a geometric frustration effect at the marginally jammed state. It is also surprising to notice that such slight particle asphericity can trigger non-negligible correlations between contact-level and mesoscale structures, manifested in drastically different nonaffine structural rearrangements upon compaction from that of granular spheres. These observations can help improve statistical mechanical models for the orientational order transformation of nonspherical granular particle packings, which involves complex interplays between particle shape, frictional contacts, and external force field.
Collapse
Affiliation(s)
| | | | | | | | | | - Yujie Wang
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
- Department of Physics, College of Mathematics and Physics, Chengdu University of Technology, Chengdu 610059, China
| | | | | | | |
Collapse
|
2
|
Anzivino C, Ness C, Moussa AS, Zaccone A. Shear flow of non-Brownian rod-sphere mixtures near jamming. Phys Rev E 2024; 109:L042601. [PMID: 38755845 DOI: 10.1103/physreve.109.l042601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2023] [Accepted: 04/05/2024] [Indexed: 05/18/2024]
Abstract
We use the discrete element method, taking particle contact and hydrodynamic lubrication into account, to unveil the shear rheology of suspensions of frictionless non-Brownian rods in the dense packing fraction regime. We find that, analogously to the random close packing volume fraction, the shear-driven jamming point of this system varies in a nonmonotonic fashion as a function of the rod aspect ratio. The latter strongly influences how the addition of rodlike particles affects the rheological response of a suspension of frictionless non-Brownian spheres to an external shear flow. At fixed values of the total (rods plus spheres) packing fraction, the viscosity of the suspension is reduced by the addition of "short"(≤2) rods but is instead increased by the addition of "long"(≥2) rods. A mechanistic interpretation is provided in terms of packing and excluded-volume arguments.
Collapse
Affiliation(s)
- Carmine Anzivino
- Department of Physics "A. Pontremoli", University of Milan, via Celoria 16, 20133 Milan, Italy
| | - Christopher Ness
- School of Engineering, University of Edinburgh, Edinburgh EH9 3JL, United Kingdom
| | | | - Alessio Zaccone
- Department of Physics "A. Pontremoli", University of Milan, via Celoria 16, 20133 Milan, Italy
| |
Collapse
|
3
|
Heussinger C. Start-up shear of spherocylinder packings: Effect of friction. Phys Rev E 2021; 103:052903. [PMID: 34134248 DOI: 10.1103/physreve.103.052903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Accepted: 05/03/2021] [Indexed: 11/07/2022]
Abstract
We study the response to shear deformations of packings of long spherocylindrical particles that interact via frictional forces with friction coefficient μ. The packings are produced and deformed with the help of molecular dynamics simulations combined with minimization techniques performed on a GPU. We calculate the linear shear modulus g_{∞}, which is orders of magnitude larger than the modulus g_{0} in the corresponding frictionless system. The motion of the particles responsible for these large frictional forces is governed by and increases with the length ℓ of the spherocylinders. One consequence of this motion is that the shear modulus g_{∞} approaches a finite value in the limit ℓ→∞, even though the density of the packings vanishes, ρ∝ℓ^{-2}. By way of contrast, the frictionless modulus decreases to zero, g_{0}∼ℓ^{-2}, in accordance with the behavior of density. Increasing the strain beyond a value γ_{c}∼μ, the packing strain weakens from the large frictional to the smaller frictionless modulus when contacts saturate at the Coulomb inequality and start to slide. In this regime, sliding friction contributes a "yield stress" σ_{y}=g_{∞}γ_{c} and the stress behaves as σ=σ_{y}+g_{0}γ. The interplay between static and sliding friction gives rise to hysteresis in oscillatory shear simulations.
Collapse
Affiliation(s)
- Claus Heussinger
- Institute for Theoretical Physics, Georg August University Göttingen, Friedrich Hund Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
4
|
Dong J, Trulsson M. Transition from steady shear to oscillatory shear rheology of dense suspensions. Phys Rev E 2020; 102:052605. [PMID: 33327063 DOI: 10.1103/physreve.102.052605] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Accepted: 10/20/2020] [Indexed: 11/07/2022]
Abstract
Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.
Collapse
Affiliation(s)
- Junhao Dong
- Theoretical Chemistry, Lund University, 221 00 Lund, Sweden
| | | |
Collapse
|
5
|
Keta YE, Olsson P. Translational and rotational velocities in shear-driven jamming of ellipsoidal particles. Phys Rev E 2020; 102:052905. [PMID: 33327139 DOI: 10.1103/physreve.102.052905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 11/05/2020] [Indexed: 06/12/2023]
Abstract
We study shear-driven jamming of ellipsoidal particles at zero temperature with a focus on the microscopic dynamics. We find that a change from spherical particles to ellipsoids with aspect ratio α=1.02 gives dramatic changes of the microscopic dynamics with much lower translational velocities and a new role for the rotations. Whereas the velocity difference at contacts-and thereby the dissipation-in collections of spheres is dominated by the translational velocities and reduced by the rotations, the same quantity is in collections of ellipsoids instead totally dominated by the rotational velocities. By also examining the effect of different aspect ratios we find that the examined quantities show either a peak or a change in slope at α≈1.2, which thus gives evidence for a crossover between different regions of low and high aspect ratio.
Collapse
Affiliation(s)
- Yann-Edwin Keta
- Department of Physics, Umeå University, 901 87 Umeå, Sweden
- Département de Physique, École Normale Supérieure de Lyon, 69364 Lyon Cedex 07, France
- Département de Physique, Université Claude Bernard Lyon 1, 69622 Villeurbanne Cedex, France
| | - Peter Olsson
- Department of Physics, Umeå University, 901 87 Umeå, Sweden
| |
Collapse
|
6
|
Marschall TA, Teitel S. Depletion forces in athermally sheared mixtures of frictionless disks and rods in two dimensions. Phys Rev E 2020; 102:042908. [PMID: 33212568 DOI: 10.1103/physreve.102.042908] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Accepted: 10/01/2020] [Indexed: 11/07/2022]
Abstract
We carry out numerical simulations to study the behavior of an athermal mixture of frictionless circular disks and elongated rods in two dimensions, under three different types of global linear deformation at a finite strain rate: (i) simple shearing, (ii) pure shearing, and (iii) isotropic compression. We find that the fluctuations induced by such deformations lead to depletion forces that cause rods to group in parallel oriented clusters for the cases of simple and pure shear, but not for isotropic compression. For simple shearing, we find that as the fraction of rods increases, this clustering increases, leading to an increase in the average rate of rotation of the rods, and a decrease in the magnitude of their nematic ordering.
Collapse
Affiliation(s)
- Theodore A Marschall
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - S Teitel
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
7
|
Binaree T, Azéma E, Estrada N, Renouf M, Preechawuttipong I. Combined effects of contact friction and particle shape on strength properties and microstructure of sheared granular media. Phys Rev E 2020; 102:022901. [PMID: 32942352 DOI: 10.1103/physreve.102.022901] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2020] [Accepted: 07/16/2020] [Indexed: 11/07/2022]
Abstract
We present a systematic numerical investigation concerning the combined effects of sliding friction and particle shape (i.e., angularity) parameters on the shear strength and microstructure of granular packings. Sliding friction at contacts varied from 0 (frictionless particles) to 0.7, and the particles were irregular polygons with an increasing number of sides, ranging from triangles to disks. We find that the effect of local friction on shear strength follows the same trend for all shapes. Strength first increases with local friction and then saturates at a shape-dependent value. In contrast, the effect of angularity varies, depending on the level of sliding friction. For low friction values (i.e., under 0.3), the strength first increases with angularity and then declines for the most angular shapes. For high friction values, strength systematically increases with angularity. At the microscale, we focus on the connectivity and texture of the contact and force networks. In general terms, increasing local friction causes these networks to be less connected and more anisotropic. In contrast, increasing particle angularity may change the network topology in different directions, directly affecting the macroscopic shear strength. These analyses and data constitute a first step toward understanding the joint effect of local variables such as friction and grain shape on the macroscopic rheology of granular systems.
Collapse
Affiliation(s)
- Theechalit Binaree
- Department of Mechanical Engineering, Faculty of Engineering, Chiang Mai University, 239 Huay Kaew Road, Chiang Mai 50200, Thailand
| | - Emilien Azéma
- LMGC, Université de Montpellier, CNRS, Montpellier, France
| | - Nicolas Estrada
- Departamento de Ingeniería Civil y Ambiental, Universidad de Los Andes, Bogotá, Colombia
| | - Mathieu Renouf
- LMGC, Université de Montpellier, CNRS, Montpellier, France
| | - Itthichai Preechawuttipong
- Department of Mechanical Engineering, Faculty of Engineering, Chiang Mai University, 239 Huay Kaew Road, Chiang Mai 50200, Thailand
| |
Collapse
|
8
|
Heussinger C. Packings of frictionless spherocylinders. Phys Rev E 2020; 102:022903. [PMID: 32942494 DOI: 10.1103/physreve.102.022903] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2020] [Accepted: 07/01/2020] [Indexed: 06/11/2023]
Abstract
We present simulation results on the properties of packings of frictionless spherocylindrical particles. Starting from a random distribution of particles in space, a packing is produced by minimizing the potential energy of interparticle contacts until a force-equilibrated state is reached. For different particle aspect ratios α=10⋯40, we calculate contacts z, pressure as well as bulk and shear modulus. Most important is the fraction f_{0}(α) of spherocylinders with contacts at both ends, as it governs the jamming threshold z_{c}(α)=8+2f_{0}(α). These results highlight the important role of the axial "sliding" degree of freedom of a spherocylinder, which is a zero-energy mode but only if no end contacts are present.
Collapse
Affiliation(s)
- Claus Heussinger
- Institute for Theoretical Physics, Georg August University Göttingen, 37077 Göttingen, Germany
| |
Collapse
|
9
|
Marschall TA, Teitel S. Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Spatial structure and correlations. Phys Rev E 2020; 101:032907. [PMID: 32289919 DOI: 10.1103/physreve.101.032907] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2020] [Accepted: 02/29/2020] [Indexed: 06/11/2023]
Abstract
We use numerical simulations to study the flow of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. We study the resulting spatial structure of the sheared system, and compute correlation functions of the velocity, the particle density, the nematic order parameter, and the particle angular velocity. Correlations of density, nematic order, and angular velocity are shown to be short ranged both below and above jamming. We compare a system of size-bidisperse particles with a system of size-monodisperse particles, and argue how differences in spatial order as the packing increases lead to differences in the global nematic order parameter. We consider the effect of shearing on initially well ordered configurations, and show that in many cases the shearing acts to destroy the order, leading to the same steady-state ensemble as found when starting from random initial configurations.
Collapse
Affiliation(s)
- Theodore A Marschall
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - S Teitel
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
10
|
Marschall TA, Van Hoesen D, Teitel S. Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Particle rotations and orientational ordering. Phys Rev E 2020; 101:032901. [PMID: 32290000 DOI: 10.1103/physreve.101.032901] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 02/15/2020] [Indexed: 06/11/2023]
Abstract
We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. Considering a range of packing fractions ϕ and particle asphericities α at small γ[over ̇], we study the angular rotation θ[over ̇]_{i} and the nematic orientational ordering S_{2} of the particles induced by the shear flow, finding a nonmonotonic behavior as the packing ϕ is varied. We interpret this nonmonotonic behavior as a crossover from dilute systems at small ϕ, where single-particle-like behavior occurs, to dense systems at large ϕ, where the geometry of the dense packing dominates and a random Poisson-like process for particle rotations results. We also argue that the finite nematic ordering S_{2} is a consequence of the shearing serving as an ordering field, rather than a result of long-range cooperative behavior among the particles. We arrive at these conclusions by consideration of (i) the distribution of waiting times for a particle to rotate by π, (ii) the behavior of the system under pure, as compared to simple, shearing, (iii) the relaxation of the nematic order parameter S_{2} when perturbed away from the steady state, and (iv) by construction, a numerical mean-field model for the rotational motion of a particle. Our results also help to explain the singular behavior observed when taking the α→0 limit approaching circular disks.
Collapse
Affiliation(s)
- Theodore A Marschall
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Daniel Van Hoesen
- Department of Physics, Washington University, St. Louis, Missouri 63130, USA
| | - S Teitel
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|