Strong interlocking of nonconvex particles in random packings

Geometric cohesion in two-dimensional systems composed of star-shaped particles

Using a discrete element method, we investigate the phenomenon of geometric cohesion in granular systems composed of star-shaped particles with 3 to 13 arms. This was done by analyzing the stability of columns built with these particles and by studying the microstructure of these columns in terms of density and connectivity. We find that systems composed of star-shaped particles can exhibit geometric cohesion (i.e., a solidlike behavior, in the absence of adhesive forces between the grains), depending on the shape of the particles and the friction between them. This phenomenon is observed up to a given critical size of the system, from which a transition to a metastable behavior takes place. We also have evidence that geometric cohesion is closely linked to the systems' connectivity and especially to the capability of forming interlocked interactions (i.e., multicontact interactions that hinder the relative rotation of the grains). Our results contribute to the understanding of the interesting and potentially useful phenomenon of geometric cohesion. In addition, our work supplements an important set of experimental observations and sheds light on the complex behavior of real, three-dimensional, granular systems.

Investigation on the Influencing Factors of K0 of Granular Materials Using Discrete Element Modelling

Earth pressure coefficient at rest K0 is commonly estimated by empirical equations, which to date has had insufficient accuracy and universality. For better prediction, the investigation on the factors influencing K0 is required. A series of discrete element method (DEM) simulations of oedometer tests are conducted to verify the key factors influencing K0 of granular materials. The influences of initial fabric anisotropy, particle shape, initial void ratio, inter-particle friction angle is investigated. The evolution of microstructure is monitored during the tests to reveal the relationship between the microstructure evolution and K0 values. The results show that the effect of fabric anisotropy exists but is limited. Particle shape, initial void ratio, and inter-particle friction angle all significantly affect the K0 values alone. According to the DEM results, an attempt is made to propose a more reasonable empirical equation in which K0 is a function of relative density, critical state friction angle, and “shape factor”. This new empirical equation has higher accuracy and can consider the effect of particle shape, inspiring the determination of K0 values in practical engineering.

Link between packing morphology and the distribution of contact forces and stresses in packings of highly nonconvex particles

An external load on a particle packing is distributed internally through a heterogeneous network of particle contacts. This contact force distribution determines the stability of the particle packing and the resulting structure. Here, we investigate the homogeneity of the contact force distribution in packings of highly nonconvex particles both in two-dimensional (2D) and three-dimensional (3D) packings. A recently developed discrete element method is used to model packings of nonconvex particles of varying sphericity. Our results establish that in 3D packings the distribution of the contact forces in the normal direction becomes increasingly heterogeneous with decreasing particle sphericity. However, in 2D packings the contact force distribution is independent of particle sphericity, indicating that results obtained in 2D packings cannot be extrapolated readily to 3D packings. Radial distribution functions show that the crystallinity in 3D packings decreases with decreasing particle sphericity. We link the decreasing homogeneity of the contact force distributions to the decreasing crystallinity of 3D packings. These findings are complementary to the previously observed link between the heterogeneity of the contact force distribution and a decreasing packing crystallinity due to an increasing polydispersity of spherical particles.

Two-stage athermal solidification of semiflexible polymers and fibers

We study how solidification of model freely rotating polymers under athermal quasistatic compression varies with their bond angle θ0. All systems undergo two discrete, first-order-like transitions: entanglement at φ = φE(θ0) followed by jamming at φ = φJ(θ0) ≃ (4/3 ± 1/12)φE(θ0). For φ < φE(θ0), systems are in a "gas" phase wherein all chains remain free to translate and reorient. For φE(θ0) ≤ φ ≤ φJ(θ0), systems are in a liquid-like phase wherein chains are entangled. In this phase, chains' rigid-body-like motion is blocked, yet they can still locally relax via dihedral rotations, and hence energy and pressure remain extremely small. The ability of dihedral relaxation mechanisms to accommodate further compression becomes exhausted, and systems rigidify, at φJ(θ0). At and slightly above φJ, the bulk moduli increase linearly with the pressure P rather than jumping discontinuously, indicating these systems solidify via rigidity percolation. The character of the energy and pressure increases above φJ(θ0) can be characterized via chains' effective aspect ratio αeff. Large-αeff (small-θ0) systems' jamming is bending-dominated and is similar to that observed in systems composed of straight fibers. Small-αeff (large-θ0) systems' jamming is dominated by the degree to which individual chains' dihedrals can collapse into compact, tetrahedron-like structures. For intermediate θ0, chains remain in highly disordered globule-like configurations throughout the compression process; jamming occurs when entangled globules can no longer even locally relax away from one another.

Densest versus jammed packings of bent-core trimers

We identify putatively maximally dense packings of tangent-sphere trimers with fixed bond angles (θ=θ_{0}), and contrast them to the disordered jammed states they form under quasistatic and dynamic athermal compression. Incommensurability of θ_{0} with three-dimensional (3D) close packing does not by itself inhibit formation of dense 3D crystals; all θ_{0} allow formation of crystals with ϕ_{max}(θ_{0})>0.97ϕ_{cp}. Trimers are always able to arrange into periodic structures composed of close-packed bilayers or trilayers of triangular-lattice planes, separated by "gap layers" that accommodate the incommensurability. All systems have ϕ_{J} significantly below the monomeric value, indicating that trimers' quenched bond-length and bond-angle constraints always act to promote jamming. ϕ_{J} varies strongly with θ_{0}; straight (θ_{0}=0) trimers minimize ϕ_{J} while closed (θ_{0}=120^{∘}) trimers maximize it. Marginally jammed states of trimers with lower ϕ_{J}(θ_{0}) exhibit quantifiably greater disorder, and the lower ϕ_{J} for small θ_{0} is apparently caused by trimers' decreasing effective configurational freedom as they approach linearity.

Cluster growth in driven granular gases

We investigate numerically and theoretically the internal structures of a driven granular gas in cuboidal cell geometries. Clustering is reported and particles are classified as gaseous or clustered via a local packing fraction criterion based on a Voronoi tessellation. We observe that small clusters arise in the corners of the box, elucidating early reports of partial clustering. These aggregates have a condensation-like surface growth. When a critical size is reached, a structural transition occurs and all clusters merge together, leaving a hole in the center of the cell. This hole then becomes the new center of particle capture. Taking into account all structural modifications and defining a saturation packing fraction, we propose an empirical model for the cluster growth.

Optimizing packing fraction in granular media composed of overlapping spheres

What particle shape will generate the highest packing fraction when randomly poured into a container? In order to explore and navigate the enormous search space efficiently, we pair molecular dynamics simulations with artificial evolution. Arbitrary particle shape is represented by a set of overlapping spheres of varying diameter, enabling us to approximate smooth surfaces with a resolution proportional to the number of spheres included. We discover a family of planar triangular particles, whose packing fraction of ϕ ∼ 0.73 is among the highest experimental results for disordered packings of frictionless particles. We investigate how ϕ depends on the arrangement of spheres comprising an individual particle and on the smoothness of the surface. We validate the simulations with experiments using 3D-printed copies of the simplest member of the family, a planar particle consisting of three overlapping spheres with identical radius. Direct experimental comparison with 3D-printed aspherical ellipsoids demonstrates that the triangular particles pack exceedingly well not only in the limit of large system size but also when confined to small containers.

Granular transport in driven granular gas

We numerically and theoretically investigate the behavior of a granular gas driven by asymmetric plates. The injection of energy in the dissipative system differs from one side to the opposite one. We prove that the dynamical clustering which is expected for such a system is affected by the asymmetry. As a consequence, the cluster position can be fully controlled. This property could lead to various applications in the handling of granular materials in low-gravity environment. Moreover, the dynamical cluster is characterized by natural oscillations which are also captured by a model. These oscillations are mainly related to the cluster size, thus providing an original way to probe the clustering behavior.

Celebrating Soft Matter's 10th Anniversary: toward jamming by design

In materials science, high performance is typically associated with regularity and order, while disorder and the presence of defects are assumed to lead to sub-optimal outcomes. This holds for traditional solids such as crystals as well as for many types of nanoscale devices. However, there are circumstances where disorder can be harnessed to achieve performance not possible with approaches based on regularity. Recent research has shown opportunities specifically for soft matter. There, the phenomenon of jamming leads to unique emergent behavior that enables disordered, amorphous systems to switch reversibly between solid-like rigidity and fluid-like plasticity. This makes it possible to envision materials that can change stiffness or even shape adaptively. We review some of the progress in this direction, discussing examples where jamming has been explored from micro to macro scales in colloidal systems, suspensions, granular-materials-enabled soft robotics, and architecture. We focus in particular on how the jammed aggregate state can be tailored by controlling particle level properties and discuss very recent ideas that provide an important first step toward actual design of specifically targeted jamming behavior.

Particle shape effects on the stress response of granular packings

We present measurements of the stress response of packings formed from a wide range of particle shapes. Besides spheres these include convex shapes such as the Platonic solids, truncated tetrahedra, and triangular bipyramids, as well as more complex, non-convex geometries such as hexapods with various arm lengths, dolos, and tetrahedral frames. All particles were 3D-printed in hard resin. Well-defined initial packing states were established through preconditioning by cyclic loading under given confinement pressure. Starting from such initial states, stress-strain relationships for axial compression were obtained at four different confining pressures for each particle type. While confining pressure has the largest overall effect on the mechanical response, we find that particle shape controls the details of the stress-strain curves and can be used to tune packing stiffness and yielding. By correlating the experimentally measured values for the effective Young's modulus under compression, yield stress and energy loss during cyclic loading, we identify trends among the various shapes that allow for designing a packing's aggregate behavior.

How dynamical clustering triggers Maxwell's demon in microgravity

In microgravity, the gathering of granular material can be achieved by a dynamical clustering whose existence depends on the geometry of the cell that contains the particles and the energy that is injected into the system. By compartmentalizing the cell in several subcells of smaller volume, local clustering is triggered and the so formed dense regions act as stable traps. In this paper, molecular dynamics simulations were performed in order to reproduce the phenomenon and to analyze the formation and the stability of such traps. Depending on the total number N of particles present in the whole system, several clustering modes are encountered and a corresponding bifurcation diagram is presented. Moreover, an iterative model based on the measured particle flux F as well as a theoretical model giving the asymptotical steady states are used to validate our results. The obtained results are promising and can provide ways to manipulate grains in microgravity.

Cohesive granular materials composed of nonconvex particles

The macroscopic cohesion of granular materials made up of sticky particles depends on the particle shapes. We address this issue by performing contact dynamics simulations of 2D packings of nonconvex aggregates. We find that the macroscopic cohesion is strongly dependent on the strain and stress inhomogeneities developing inside the material. The largest cohesion is obtained for nearly homogeneous deformation at the beginning of unconfined axial compression and it evolves linearly with nonconvexity. Interestingly, the aggregates in a sheared packing tend to form more contacts with fewer neighboring aggregates as the degree of nonconvexity increases. We also find that shearing leads either to an isotropic distribution of tensile contacts or to the same privileged direction as that of compressive contacts.

Rheology of three-dimensional packings of aggregates: microstructure and effects of nonconvexity

We use three-dimensional contact dynamics simulations to analyze the rheological properties of granular materials composed of rigid aggregates. The aggregates are made from four overlapping spheres and described by a nonconvexity parameter depending on the relative positions of the spheres. The macroscopic and microstructural properties of several sheared packings are analyzed as a function of the degree of nonconvexity of the aggregates. We find that the internal angle of friction increases with the nonconvexity. In contrast, the packing fraction first increases to a maximum value but declines as the nonconvexity increases further. At a high level of nonconvexity, the packings are looser but show a higher shear strength. At the microscopic scale, the fabric and force anisotropy, as well as the friction mobilization, are enhanced by multiple contacts between aggregates and interlocking, thus revealings the mechanical and geometrical origins of shear strength.