Contact network topology in tapped granular media

Understanding slow compression of frictional granular particles by network analysis

We consider frictional granular packings exposed to quasi-static compression rates, with a focus on systems above the jamming transition. For frictionless packings, earlier work (S. Luding et al., Soft Matter, 2022, 18(9), 1868-1884) has uncovered that the system evolution/response involves smooth evolution phases, interrupted by fast transitions (events). The general finding is that the force networks' static quantities correlate closely with the pressure, while their evolution resembles the kinetic energy for both frictionless and frictional packings. The former represents reversible (elastic) particle deformations with affine and non-affine components, whereas the latter also involves much stronger, irreversible (plastic) rearrangements of the packings. Events are associated with jumps in the overall kinetic energy as well as dramatic changes in the force networks describing the particle micro-structure. The frictional nature of particle interactions affects both their frequency and the relevant time scale magnitude. For intermediate friction, events are often followed by an unexpected slow-down during which the kinetic energy drops below its average value. We find that these slow-downs are associated with a significant decrease in the non-affine dynamics of the particles, and are strongly influenced by friction. Friction modifies the structure of the networks, both through the typical number of contacts of a particle, and by influencing topological features of the resulting networks. Furthermore, friction modifies the dynamics of the networks, with larger values of friction leading to smaller changes of the more stable networks.

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.

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.

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.

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.

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

The force network of a granular assembly, defined by the contact network and the corresponding contact forces, carries valuable information about the state of the packing. Simple analysis of these networks based on the distribution of force strengths is rather insensitive to the changes in preparation protocols or to the types of particles. In this and the companion paper [Kondic et al., Phys. Rev. E 93, 062903 (2016)10.1103/PhysRevE.93.062903], we consider two-dimensional simulations of tapped systems built from frictional disks and pentagons, and study the structure of the force networks of granular packings by considering network's topology as force thresholds are varied. We show that the number of clusters and loops observed in the force networks as a function of the force threshold are markedly different for disks and pentagons if the tangential contact forces are considered, whereas they are surprisingly similar for the network defined by the normal forces. In particular, the results indicate that, overall, the force network is more heterogeneous for disks than for pentagons. Such differences in network properties are expected to lead to different macroscale response of the considered systems, despite the fact that averaged measures (such as force probability density function) do not show any obvious differences. Additionally, we show that the states obtained by tapping with different intensities that display similar packing fraction are difficult to distinguish based on simple topological invariants.

Scaling properties of force networks for compressed particulate systems

We consider, computationally and experimentally, the scaling properties of force networks in the systems of circular particles exposed to compression in two spatial dimensions. The simulations consider polydisperse and monodisperse particles, both frictional and frictionless, and in experiments we use monodisperse and bidisperse frictional particles. While for some of the considered systems we observe consistent scaling exponents describing the behavior of the force networks, we find that this behavior is not universal. In particular, we find that frictionless systems, independently of whether they partially crystallize under compression or not, show scaling properties that are significantly different compared to the frictional disordered ones. The findings of nonuniversality are confirmed by explicitly computing fractal dimension for the considered systems. The results of the physical experiments are consistent with the results obtained in simulations of frictional disordered systems.

Percolation and jamming transitions in particulate systems with and without cohesion

We consider percolation and jamming transitions for particulate systems exposed to compression. For the systems built of particles interacting by purely repulsive forces in addition to friction and viscous damping, it is found that these transitions are influenced by a number of effects, and in particular by the compression rate. In a quasistatic limit, we find that for the considered type of interaction between the particles, percolation and jamming transitions coincide. For cohesive systems, however, or for any system exposed to even slow dynamics, the differences between the considered transitions are found and quantified.

Evolution of force networks in dense particulate media

We discuss sets of measures with the goal of describing dynamical properties of force networks in dense particulate systems. The presented approach is based on persistent homology and allows for extracting precise, quantitative measures that describe the evolution of geometric features of the interparticle forces, without necessarily considering the details related to individual contacts between particles. The networks considered emerge from discrete element simulations of two-dimensional particulate systems consisting of compressible frictional circular disks. We quantify the evolution of the networks for slowly compressed systems undergoing jamming transition. The main findings include uncovering significant but localized changes of force networks for unjammed systems, global (systemwide) changes as the systems evolve through jamming, to be followed by significantly less dramatic evolution for the jammed states. We consider both connected components, related in a loose sense to force chains, and loops and find that both measures provide a significant insight into the evolution of force networks. In addition to normal, we consider also tangential forces between the particles and find that they evolve in the consistent manner. Consideration of both frictional and frictionless systems leads us to the conclusion that friction plays a significant role in determining the dynamical properties of the considered networks. We find that the proposed approach describes the considered networks in a precise yet tractable manner, making it possible to identify features which could be difficult or impossible to describe using other approaches.

Topological analysis of tapped granular media using persistent homology

We use the first Betti number of a complex to analyze the morphological structure of granular samples in mechanical equilibrium. We investigate two-dimensional granular packings after a tapping process by means of both simulations and experiments. States with equal packing fraction obtained with different tapping intensities are distinguished after the introduction of a filtration parameter which determines the particles (nodes in the network) that are joined by an edge. This is accomplished by just using the position of the particles obtained experimentally and no other information about the possible contacts, or magnitude of forces.