Structure of force networks in tapped particulate systems of disks and pentagons. I. Clusters and loops

Jammed solids with pins: Thresholds, force networks, and elasticity

The role of fixed degrees of freedom in soft or granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our two-dimensional simulated system consists of soft particles and fixed "pins," both of which harmonically repel overlaps. On the one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology, we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.

On intermittency in sheared granular systems

We consider a system of granular particles, modeled by two dimensional frictional soft elastic disks, that is exposed to externally applied time-dependent shear stress in a planar Couette geometry. We concentrate on the external forcing that produces intermittent dynamics of stick-slip type. In this regime, the top wall remains almost at rest until the applied stress becomes sufficiently large, and then it slips. We focus on the evolution of the system as it approaches a slip event. Our main finding is that there are two distinct groups of measures describing system behavior before a slip event. The first group consists of global measures defined as system-wide averages at a fixed time. Typical examples of measures in this group are averages of the normal or tangent forces acting between the particles, system size and number of contacts between the particles. These measures do not seem to be sensitive to an approaching slip event. On average, they tend to increase linearly with the force pulling the spring. The second group consists of the time-dependent measures that quantify the evolution of the system on a micro (particle) or mesoscale. Measures in this group first quantify the temporal differences between two states and only then aggregate them to a single number. For example, Wasserstein distance quantitatively measures the changes of the force network as it evolves in time while the number of broken contacts quantifies the evolution of the contact network. The behavior of the measures in the second group changes dramatically before a slip event starts. They increase rapidly as a slip event approaches, indicating a significant increase in fluctuations of the system before a slip event is triggered.

Understanding slow compression and decompression of frictionless soft granular matter by network analysis

We consider dense granular systems in three spatial dimensions exposed to slow compression and decompression, below, during, above and well above jamming. The evolution of granular systems under slow deformation is non-trivial and involves smooth, continuous, reversible (de)compression periods, interrupted by fast, discontinuous, irreversible transition events. These events are often, but not always, associated with rearrangements of particles and of the contact network. How many particles are involved in these transitions between two states can range from few to almost all in the system. An analysis of the force network that is built on top of the contact network is carried out using the tools of persistent homology. Results involve the observation that kinetic energy is correlated with the intensity of rearrangements, while the evolution of global mechanical measures, such as pressure, is strongly correlated with the evolution of the topological measures quantifying loops in the force network. Surprisingly, some transitions are clearly detected by persistent homology even though motion/rearrangement of particles is much weaker, i.e., much harder to detect or, in some cases, not observed at all.

Stress propagation in locally loaded packings of disks and pentagons

The mechanical strength and flow of granular materials can depend strongly on the shapes of individual grains. We report quantitative results obtained from photoelasticimetry experiments on locally loaded, quasi-two-dimensional granular packings of either disks or pentagons exhibiting stick-slip dynamics. Packings of pentagons resist the intruder at significantly lower packing fractions than packings of disks, transmitting stresses from the intruder to the boundaries over a smaller spatial extent. Moreover, packings of pentagons feature significantly fewer back-bending force chains than packings of disks. Data obtained on the forward spatial extent of stresses and back-bending force chains collapse when the packing fraction is rescaled according to the packing fraction of steady state open channel formation, though data on intruder forces and dynamics do not collapse. We comment on the influence of system size on these findings and highlight connections with the dynamics of the disks and pentagons during slip events.

Micromechanical description of the compaction of soft pentagon assemblies

We analyze the isotropic compaction of assemblies composed of soft pentagons interacting through classical Coulomb friction via numerical simulations. The effect of the initial particle shape is discussed by comparing packings of pentagons with packings of soft circular particles. We characterize the evolution of the packing fraction, the elastic modulus, and the microstructure (particle rearrangement, connectivity, contact force, and particle stress distributions) as a function of the applied stresses. Both systems behave similarly: the packing fraction increases and tends asymptotically to a maximum value ϕ_{max}, where the bulk modulus diverges. At the microscopic scale we show that particle rearrangements occur even beyond the jammed state, the mean coordination increases as a square root of the packing fraction, and the force and stress distributions become more homogeneous as the packing fraction increases. Soft pentagons experience larger particle rearrangements than circular particles, and such behavior decreases proportionally to the friction. Interestingly, the friction between particles also contributes to a better homogenization of the contact force network in both systems. From the expression of the granular stress tensor we develop a model that describes the compaction behavior as a function of the applied pressure, the Young modulus, and the initial shape of the particles. This model, settled on the joint evolution of the particle connectivity and the contact stress, provides outstanding predictions from the jamming point up to very high densities.

Enhancement of the flow of vibrated grains through narrow apertures by addition of small particles

We analyze the flow and clogging of circular grains passing through a small aperture under vibration in two dimensions. Via discrete element method simulations, we show that when grains smaller than the original ones are introduced in the system as an additive, the net flow of the original species can be significantly increased. Moreover, there is an optimal radius of the additive particles that maximizes the effect. This finding may constitute the basis for technological applications not only concerning the flow of granular materials but also regarding active matter, including pedestrian evacuation.

Self-assembly of self-propelled magnetic grains

In this work, we study bidisperse mixtures of self-propelled magnetic particles of different shapes via discrete element method simulations. We show how these particles self-assemble into clusters and how these clusters depend on the ratio of the mixture, the magnetic interaction, and the shape of the grains. It is found that the mix ratio of the system controls the cluster size. Besides, the intensity of the magnetic dipoles and the shape of the grains in the mixture rule the average number of neighbors in contact and the shape of the clusters. By varying the intensity of the interactions, globular, linear and branched clusters were obtained.

Failure of confined granular media due to pullout of an intruder: from force networks to a system wide response

We investigate computationally the pullout of a spherical intruder initially buried at the bottom of a granular column. The intruder starts to move out of the granular bed once the pulling force reaches a critical value, leading to material failure. The failure point is found to depend on the diameter of the granular column, pointing out the importance of particle-wall interactions in determining the material response. Discrete element simulations show that prior to failure, the contact network is essentially static, with only minor rearrangements of the particles. However, the force network, which includes not only the contact information, but also the information about the interaction strength, undergoes nontrivial evolution. An initial insight is obtained by considering the relative magnitudes of normal and tangential forces between the particles, and in particular the proportion of contacts that reach Coulomb threshold. More detailed understanding of the processes leading to failure is reached by the analysis of both spatial and temporal properties of the force network using the tools of persistent homology. We find that the forces between the particles undergo intermittent temporal variations ahead of the failure. In addition to this temporal intermittency, the response of the force network is found to be spatially dependent and influenced by proximity to the intruder. Furthermore, the response is modified significantly by the interaction strength, with the relevant measures describing the response showing differing behaviors for the contacts characterized by large interaction forces.

Betweenness centrality as predictor for forces in granular packings

A load applied to a jammed frictional granular system will be localized into a network of force chains making inter-particle connections throughout the system. Because such systems are typically under-constrained, the observed force network is not unique to a given particle configuration, but instead varies upon repeated formation. In this paper, we examine the ensemble of force chain configurations created under repeated assembly in order to develop tools to statistically forecast the observed force network. In experiments on a gently suspended 2D layer of photoelastic particles, we subject the assembly to hundreds of repeated cyclic compressions. As expected, we observe the non-unique nature of the force network, which differs for each compression cycle, by measuring all vector inter-particle contact forces using our open source PeGS software. We find that total pressure on each particle in the system correlates to its betweenness centrality value extracted from the geometric contact network. Thus, the mesoscale network structure is a key control on individual particle pressures.

Characterizing granular networks using topological metrics

We carry out a direct comparison of experimental and numerical realizations of the exact same granular system as it undergoes shear jamming. We adjust the numerical methods used to optimally represent the experimental settings and outcomes up to microscopic contact force dynamics. Measures presented here range from microscopic through mesoscopic to systemwide characteristics of the system. Topological properties of the mesoscopic force networks provide a key link between microscales and macroscales. We report two main findings: (1) The number of particles in the packing that have at least two contacts is a good predictor for the mechanical state of the system, regardless of strain history and packing density. All measures explored in both experiments and numerics, including stress-tensor-derived measures and contact numbers depend in a universal manner on the fraction of nonrattler particles, f_{NR}. (2) The force network topology also tends to show this universality, yet the shape of the master curve depends much more on the details of the numerical simulations. In particular we show that adding force noise to the numerical data set can significantly alter the topological features in the data. We conclude that both f_{NR} and topological metrics are useful measures to consider when quantifying the state of a granular system.

Granular response to impact: Topology of the force networks

The impact of an intruder on granular matter leads to the formation of mesoscopic force networks, which were seen particularly clearly in the recent experiments carried out with photoelastic particles [Clark et al., Phys. Rev. Lett. 114, 144502 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.144502]. These force networks are characterized by complex structure and evolve on fast time scales. While it is known that total photoelastic activity in the granular system is correlated with the acceleration of the intruder, it is not known how the structure of the force network evolves during impact, and if there are dominant features in the networks that can be used to describe the intruder's dynamics. Here, we use topological tools, in particular persistent homology, to describe these features. Persistent homology allows quantification of both structure and time evolution of the resulting force networks. We find that there is a clear correlation of the intruder's dynamics and some of the topological measures implemented. This finding allows us to discuss which properties of the force networks are most important when attempting to describe the intruder's dynamics. In particular, we find that the presence of loops in the force network, quantified by persistent homology, is strongly correlated to the deceleration of the intruder. In some cases, particularly for the impact on soft particles, the measures derived from the persistence analysis describe the deceleration of the intruder even better than the total photoelastic activity. We are also able to define an upper bound on the relevant time scale over which the force networks evolve.

Quantifying disorder in colloidal films spin-coated onto patterned substrates

Polycrystals of thin colloidal deposits, with thickness controlled by spin-coating speed, exhibit axial symmetry with local 4-fold and 6-fold symmetric structures, termed orientationally correlated polycrystals (OCPs). While spin-coating is a very facile technique for producing large-area colloidal deposits, the axial symmetry prevents us from achieving true long-range order. To obtain true long-range order, we break this axial symmetry by introducing a patterned surface topography and thus eliminate the OCP character. We then examine symmetry-independent methods to quantify order in these disordered colloidal deposits. We find that all the information in the bond-orientational order parameters is well captured by persistent homology analysis methods that only use the centers of the particles as input data. It is expected that these methods will prove useful in characterizing other disordered structures.

Evolution of network architecture in a granular material under compression

As a granular material is compressed, the particles and forces within the system arrange to form complex and heterogeneous collective structures. Force chains are a prime example of such structures, and are thought to constrain bulk properties such as mechanical stability and acoustic transmission. However, capturing and characterizing the evolving nature of the intrinsic inhomogeneity and mesoscale architecture of granular systems can be challenging. A growing body of work has shown that graph theoretic approaches may provide a useful foundation for tackling these problems. Here, we extend the current approaches by utilizing multilayer networks as a framework for directly quantifying the progression of mesoscale architecture in a compressed granular system. We examine a quasi-two-dimensional aggregate of photoelastic disks, subject to biaxial compressions through a series of small, quasistatic steps. Treating particles as network nodes and interparticle forces as network edges, we construct a multilayer network for the system by linking together the series of static force networks that exist at each strain step. We then extract the inherent mesoscale structure from the system by using a generalization of community detection methods to multilayer networks, and we define quantitative measures to characterize the changes in this structure throughout the compression process. We separately consider the network of normal and tangential forces, and find that they display a different progression throughout compression. To test the sensitivity of the network model to particle properties, we examine whether the method can distinguish a subsystem of low-friction particles within a bath of higher-friction particles. We find that this can be achieved by considering the network of tangential forces, and that the community structure is better able to separate the subsystem than a purely local measure of interparticle forces alone. The results discussed throughout this study suggest that these network science techniques may provide a direct way to compare and classify data from systems under different external conditions or with different physical makeup.

Topological and geometric measurements of force-chain structure

Developing quantitative methods for characterizing structural properties of force chains in densely packed granular media is an important step toward understanding or predicting large-scale physical properties of a packing. A promising framework in which to develop such methods is network science, which can be used to translate particle locations and force contacts into a graph in which particles are represented by nodes and forces between particles are represented by weighted edges. Recent work applying network-based community-detection techniques to extract force chains opens the door to developing statistics of force-chain structure, with the goal of identifying geometric and topological differences across packings, and providing a foundation on which to build predictions of bulk material properties from mesoscale network features. Here we discuss a trio of related but fundamentally distinct measurements of the mesoscale structure of force chains in two-dimensional (2D) packings, including a statistic derived using tools from algebraic topology, which together provide a tool set for the analysis of force chain architecture. We demonstrate the utility of this tool set by detecting variations in force-chain architecture with pressure. Collectively, these techniques can be generalized to 3D packings, and to the assessment of continuous deformations of packings under stress or strain.

Structure of force networks in tapped particulate systems of disks and pentagons. II. Persistence analysis

In the companion paper [Pugnaloni et al., Phys. Rev. E 93, 062902 (2016)10.1103/PhysRevE.93.062902], we use classical measures based on force probability density functions (PDFs), as well as Betti numbers (quantifying the number of components, related to force chains, and loops), to describe the force networks in tapped systems of disks and pentagons. In the present work, we focus on the use of persistence analysis, which allows us to describe these networks in much more detail. This approach allows us not only to describe but also to quantify the differences between the force networks in different realizations of a system, in different parts of the considered domain, or in different systems. We show that persistence analysis clearly distinguishes the systems that are very difficult or impossible to differentiate using other means. One important finding is that the differences in force networks between disks and pentagons are most apparent when loops are considered: the quantities describing properties of the loops may differ significantly even if other measures (properties of components, Betti numbers, force PDFs, or the stress tensor) do not distinguish clearly or at all the investigated systems.