I. Introduction
The importance of granular matter for human activities is discreetly striking. It not only covers 90% of the Earth’s crust, but also accounts for about 70% of the materials used in industry, represents the second-largest volume of matter carried by human beings - just behind water -, and consumes more energy in its production than that consumed in human transportation. For all these reasons, a thorough understanding of the behavior of granular matter, made of all sorts of materials, is of tremendous importance for human endeavours ranging from environmental preservation and control of geohazards to space exploration. Yet, this is an extremely complex task since these materials span several scales, from colloids to screes to asteroids. These materials can also be quite different - made of virtually any sort of matter - they can behave like a solid, a liquid, or a gas depending on loading conditions and local interactions. Over the last three decades, a framework and state equations to describe the flow of granular systems made of rigid, virtually undeformable, particles have been set up (4-6). Similarly, the jamming transition of these rigid grains has been fruitfully explored and these findings have been extended to the case of soft grains (7-10). Nevertheless, it is unsettling that in much of the literature involving deformable granular media the material these particles are made of is assumed to stay in the linear elastic regime so that the Hertz contact assumption remains valid (11, 12). This means that the level of this material deformation does not exceed 10%. However, we can naturally wonder about the consequences of such an assumption in granular systems undergoing more intense loading, well above 10% deformation (13). The existing framework mostly neglects the grains’ change in shape, limiting the number of real-life scenarios that can be described: This approach is good at describing the loading of corn grains but not of pop-corn. As shown in Fig. 1, many systems, ranging from biology to industrial applications, are made of particles brought to large deformation regimes. To differentiate these systems from those composed of soft grains, we propose to call squishy granular materials all discrete matter whose individual bodies undergo large nonlinear deformations. This definition implies that the global stresses undergone by the system are substantially higher than the Young modulus of the particles, i.e., Ematerial ≪ σloading.
In addition, stating that a granular system is squishy also means that the particles’ change of shape must be considered to describe their local interactions, displacement, and, more generally, to understand the global behavior of the granular medium. In other words, the exact nature of a squishy material cannot be hidden behind a linear elastic law. The intrinsic rheology of the particles becomes paramount and must be considered through the nonlinear theory of elasticity (18). Going further, the exact geometry of the particles must be thoroughly considered since it can, for example, induce discontinuous and sudden phenomena such as buckling. This is the case for shell particles (19,20) or, more generally, for particles made of any metamaterial (21,22). Naturally, a large diversity of geometries and material behaviors must be considered to properly handle squishy media. This can encompass thermal and anisotropic effects, as well as magnetic interactions (23, 24). As a limitation in the broad category of squishy particles, we assume that the number of entities in the system stays stable all along the loading process. The squeezable nature of the particles allows them to support any loading with no merging or fragmentation. Thus, the edge of each particle is continuously described by a closed function.
Fortunately, the granular matter physics field is not starting from zero as far as it concerns squishy particles. The most advanced studies in this domain come from the rheology of liquid foams (25), materials composed of gas bubbles within liquids. Numerous experimental and numerical works have significantly contributed to the understanding of these materials under different loadings (14,26-28) and direct links have been drawn to granular matter physics (29). However, emulsions and colloids wildly differ from other discrete materials when considering their material rheology and interparticle frictional properties (30-33). Consequently, it is difficult to generalize the results found for foams, colloids, and emulsions to the case of squishy granular matter. Motivated by technological applications, very soft grains have also been studied under compression, as in the case of sintered material manufacturing. The compaction of ceramic, metal, gels, and drug powders has been central in the research of compacted materials (34-36) or selfformed colloids (37). In all these cases, however, studies are still limited to the plastic behavior of particles under compression. Other examples of squishy grains can be found in living organisms. Cells forming a tissue or clogging in a vein can certainly be described using squishy particle assemblies (16,31,38-40).
From a numerical point of view, only recently a number of models have emerged for squishy granular matter. Some of them, based on the discrete element method (DEM), have allowed researchers to catch a glimpse of the behavior of multiparticle, highly deformable systems (20, 41, 42), but remain unsatisfactory in terms of material mimicking. Other methods based on finite-elements (FEM) are better at mimicking particles’ materials, but have difficulties to properly take contact interactions into account (43-46). In both cases, major improvements are necessary to simulate any squishy granular system.
The rest of this article is dedicated to summarizing the different loading conditions that squishy granular matter can undergo and various approaches to dealing with their characterization, simulation, and analysis. More particularly, we focus on the experimental and numerical challenges to understanding the behavior of squishy grains under compression and shearing.
II. Browsing deep in the jammed state i Compressing squishy particle assemblies
The mechanical behavior of a granular system is ruled by the individual behavior of each particle within it. For example, when hard particles - i.e., those made of a material with high stiffness relatively to the applied stress - are compressed inside a box, they tend to rearrange compactly. In this case, the rearrangement ends when a stable force network is reached. This state, called the jammed state, depends mainly on the morphology of the particles, their size and shape (47), and their interaction laws (47-49). If one considers assemblies of soft particles with low stiffness relatively to the applied stress, they are able to undergo large deformations, and the compression can continue beyond the jammed state. Exploring the compaction of granular matter and, especially, systems very deep in the jammed state is challenging, mainly due to the nature of the particle rearrangements and the evolution of the force-chain network. For assemblies of squishy particles, the compactness - i.e., the ratio between the volume of the particles and the volume of the box - surpasses the compactness of the jammed state, the random close packing, and tends to a maximal value close to 1. Naturally, this high compactness comes with an increasing number of contacts and contact surfaces. Therefore, it is seen as an evolution of the probability distribution of the forces and an evident reorganization of the force-chain network that drastically changes the mechanical properties of the assembly.
Innovative experiments and advanced numerical methods have allowed one to take a small step forward in understanding the microstructural evolution beyond the jamming point (see section III). However, theoretical modeling of the compaction process remains a complicated task. Although many equations have been proposed in the last decades attempting to link the confining pressure to the compactness (50-53), most of them are based on empirical strategies with no apparent physical bases, requiring several fitting parameters. Only recently, a systematic approach, based entirely on micromechanical parameters and free of any ad hoc parameters, has been proposed and numerically tested in the case of hyperelastic circular, polygonal, and spherical shapes under different frictional conditions (45, 46). This method was derived from the micromechanical expression of the granular stress tensor together with the approximation of the Hertz contact law between two particles at small strain and the theoretical solution of the hydro-static compaction of a single particle.
Although the above micromechanical approach constitutes a substantial progress in the challenge of understanding and modeling the compaction of deformable particles, it only describes convex particles under isotropic compaction with continuous laws of elasticity. Many unknowns remain to be resolved concerning particles with a more complex shape or composition characteristics, such as nonconvex particles or particles that do not follow elastic or continuous deformation behavior (for example, particles suffering buckling). Also, we still need to better understand the behavior of granular systems under different loading conditions, such as uniaxial compaction or shear loading.
ii. Shearing squishy particle assemblies
The second and more complex loading that a squishy granular system can undergo is shearing. This loading type induces not only large particle deformations but also large and often erratic particle rearrangements. When this latter regime dominates, the granular system yields and flows. Below the jammed state, for rigid grains, the flowing regime is described in many configurations by well-known state equations (4-6). Although squishy grains could present state equations reminiscent of those for rigid grains, this has not been shown yet. So far, for a dense flow of soft grains, it has been
shown that grains rearrange through a succession of local mechanisms called T1-events (14, 54) inducing specific glassy rheology. These T1-events, presented in Fig. 2(e), mainly happen in a relatively narrow area of the system called the shear band. This has been observed for liquid foams (14) and living cells (54), but one may wonder how the shear-band evolves depending on the material rheology. Do these T1-events still exist for other material types, such as plastic ones?
The shearing of granular systems often exhibits erratic dynamic mechanisms at the global scale - also known as avalanches - even under constant loading (55,56). For soft grains, depending on the particle stiffness, strain rate, and confining pressure, the global dynamics can display different statistical signatures (57). In densely packed squishy systems, the space for particles to rearrange can be extremely reduced and the material rheology has a strong influence on the possible rearrangement mechanisms. As such, it is expected that the avalanche dynamics are strongly modified by the exact particle nature and frictional properties. But we still do not know how or up to which extent this is the case.
In the case of rigid or soft particles, plasticity, seen as an irreversible deformation of the material, happens at the meso-scale: grains rearrange so that they cannot recover their previous position when the loading is released. In the case of squishy granular systems not only this plastic mechanism is at play, but also the intrinsic plasticity of the particles. Hence, the plasticity of these amorphous materials is by definition of a different nature than that of soft grains.
Concerning the shearing of squishy particles, the main challenge is to evidence the local rearrangement mechanisms, their dynamics and spatial distribution, and to understand how they induce the flowing of the sheared system. The second challenge consists in understanding the coupling between meso and micro-scale plasticity processes in squishy amorphous materials.
III. Experimental & numerical challenges i Physical experiments
i. Physical experiments
Experimentally studying granular matter is challenging due to its multi-scale nature, so tests have to be able to catch main features across scales. One needs to follow the behavior of individual grains that end up inducing the global material behavior from specific particle properties. In the case of squishy granular matter, one also needs to go one scale smaller since it is necessary to consider the material behavior inside each grain at the subgranular scale (17). Although several studies are possible without analyzing this sub-granular level or measuring the local stresses - as has been done in the early stages of granular matter physics - a comprehensive study of squishy granular matter requires such a multiscale approach.
So far, regardless of whether it be a 2D or 3D configuration, the existing ways to measure grain displacements, contacts, and stresses are all based on an inverse problem method. Usually, the Hertz contact law (11, 12), or a certain surface tension in the case of liquid foams (14), is assumed, and some local field measurements permit detection of contacts and guesses at interaction forces. These local measurements can be either the grain boundary deformation (14,30,58), a photoelastic (59,60) or thermoelastic (61) signal, or digital image correlation (DIC) rough measurements (62). By definition, these measurements are limited by the method’s sensitivity and the domain of validity of the assumed contact law. For example, in the case of photoelasticity, when deformation is too large, interference fringes overlap and the measured signal saturates. Assuming that there is no saturation of the measured signal for large deformations, it would be possible to extend these methods by finding contact laws that include complex material rheologies. For example, if the Hertz contact law could be extended beyond the case of a linear elastic behavior with a finite friction coefficient (13), this would permit extending these methods to squishy grains. Still, one can imagine that the number of guessedat parameters in the problem inversion would be much larger, making it complicated to implement.
Recently it has been proposed not to make any assumption and to directly measure displacement fields at the grain scale, even in the large deformation regime. This process is based on a large field and an accurate imaging method (13,17,63), which allows one to measure all the tensors of the non-linear elasticity theory and to guess the stress tensor knowing the loading history and the material rheology via large deformation DIC. However, this only works for small systems (∼ 100 grains) under quasistatic conditions. This method - and the fact that it provides strain fields inside the grains - brings a new outlook to this topic based upon the notion of force chain and its relevance or validity in the study of densely packed configuration.
Although recent experimental approaches have permitted to follow the evolution of squashy granular systems, the development of an experimental method to characterize squishy particle assemblies under dynamic loading remains a challenging task.
ii. Numerical simulations
The modeling of granular media is historically associated with discrete-element approaches (DEM) to describe the dynamical evolution of a collection of hard grains. DEM can be split into two main approaches: the (i) smooth and (ii) nonsmooth methods. This classification is related to the way a contact between two particles is considered. In the case (i), regularised contact laws link the contact force to a slight overlap between particles in contact. In the case (ii), the unilaterality of contacts (non-regularized laws) traduces the nonpenetration of the particles. Both approaches have proven to produce equivalent results as far as the basic hypotheses of each method are respected - i.e., binary contact and small grain deformations compared to the grain size.
Extending the numerical approaches, first to soft, and further to squishy grains, implies new assumptions and new numerical strategies. In the framework of the smooth approach, one alternative is to increase the allowed overlapping between bodies that can be later interpreted as a virtual grain deformation. Although this approach is debatable in terms of realistic representation of the grain shapes, it has shown to be valid in the small deformation domain, as in the case of the compaction of hard grains assemblies just after reaching the jammed state (65-67). Nevertheless, going deeper into the jammed state leads to the breaking of some of the model assumptions and the classical methods need to be extended to meet this new range of deformation.
As a first approximation, a squishy grain can be considered as an ‘aggregate’ of elements (in the numerical sense) bonded together by bilateral or cohesive bonds. The bonding law can be adjusted to mimic the macroscopic behaviour of a given grain material. This aggregate representation can be performed throughout the volume of the grain. This approach is known as the Bonded-Particle Method (BPM) (68-71). Another alternative is to consider sub-elements only on the particle’s surface, as proposed in the Deformable Particle Method (DPM) (20, 72) (see Fig. 2(a)). The main advantages of these methods are their low computational cost and ease of implementation in the DEM framework.
However, in the case of ‘aggregate’ methods, very large deformations can lead to non-physical situations. Another option consists in using constitutive models to describe the bulk behavior of the particles. Then, it is possible to use standard numerical methods such as a finite elements (FEM) or meshless approaches.
In a meshed approach, a volumetric rheological model associated with a state equation is considered for each grain. In this case, FEM-like methods coupled with DEM are used to describe the particle deformations, keeping - or not - its integrity (73-76). The advantage of these strategies is to closely mimic the geometry and bulk properties of the particles. However, they end up being computationally expensive. As with the classical DEM approaches, it is possible to distinguish two classes of methods: the Multi-Particle Finite Element Method (MPFEM) (73,77-79), based on regularized contact interactions, and the NonSmooth Contact Dynamic Method (NSCD) (80,81), based on non-regularized contact laws. MPFEM simulations have been performed by Mesarovic et al. (82) to study the deformation mechanisms in a periodic configuration of dissimilar elastoplastic spheres. For squishy granular packings, the first MPFEM studies were conducted by Gethin et al. in 2002 (77) and Chen et al. in 2006 (83) in 2D and 3D, in the context of metal compaction. The NSCD was used for the first time in 2000 by Acary et al. (84) with a work on masonry. They were followed by Cao et al. in 2011 (85), who tested the numerical approach for the compaction and shearing of soft packings. The latest works related to compaction of squishy granular packings using NSCD corresponds to Vu et al. (44,86) in 2D, followed by Cardenas et al. in 2D and 3D (46) (see Fig. 2(d)).
In the framework of mesh-free models, squishy grains have been simulated via the Material Point Method (MPM) (41) (see Fig. 2(b)), or the Soft Discrete Element Method (SDEM) (64) (see Fig. 2(c)), where simple kinematics are postulated in order to represent the ovalisation of the grains and their local deformations around interparticle contacts. In the last case, a continuous deformation field is interpolated between the nodes belonging to a grain by means of a moving least square algorithm. These approaches require the calibration of some fundamental computational parameters from experimental data, such as the stiffness and viscosity (64, 87). Despite the development of these alternative methods, the different results obtained via FEM-DEM compare better to experiments and show that the behavior of squishy grain collections can only be understood using these advanced numerical approaches, as time-consuming as they may be.
Most of today’s simulations involving squishy grains are performed in the quasi-static regime and under uniaxial deformation. These conditions remain relatively simple and involve few dynamical effects. Other loading conditions such as shearing would rapidly increase the complexity of the simulations, making it difficult to conduct numerical experiments at higher shear rates or compression levels. Moreover, while several thousands of squishy grains can be modeled in 2D, only a few hundred can be simulated in 3D in a reasonable computing time. Therefore, additional challenges concern (i) the development of scalable and highly parallel algorithms, (ii) the implementation of fully periodic conditions, allowing for the volume changes to account for the expansion/contraction of the particles, and (iii) considering realistic mixtures that include grains with different size, shape, and bulk behavior within the same system.
IV. Summary & perspectives
Deep in the jammed state, when grain deformation is well above 10%, the soft granular matter framework is no longer relevant when describing the evolution of granular systems. As the packing fraction grows far beyond the jamming point, it is necessary to consider the large change in particle shape, which effects: (i) creation of new contacts, (ii) interaction stresses depending on the particle material, and (iii) possible large particle rearrangements. To make up for the lack of a correct framework for these systems, we propose the concept of squishy granular materials in which the linear elastic behavior is overpassed to consider hyperelastic, plastic or other complex material rheologies.
As regards squishy materials, rapid increases in the complexity of the experimental, numerical and theoretical approaches are needed for their study. This involves several technical and scientific challenges including the understanding of the local deformation mechanisms that lead to a global compaction law. Such a law is expected to include the effect of variability in the particles’ rheology. Concerning the shear process of squishy systems, the main challenge consists in getting insights into the local rearrangement mechanisms, generalizing the T1-event concept. This should work as a building block for the description of the deformation of amorphous squishy materials. Studying the dynamics of this phenomenon and the spatial distribution of these events will also constitute a major step forward in understanding the plastic processes in these systems.
These challenges cannot be overcome without tackling serious numerical and experimental difficulties. The first consists in finding an experimental method to measure the evolution of dynamically loaded 2D and 3D systems down to the sub-granular scale. The numerical counterpart of this task is equally challenging since simulation methods can currently deal primarily with the quasistatic compression of a limited number of grains. The complexity of the simulations can rapidly increase as soon as shearing tests are considered, or the implementation of realistic and strongly deformed particle geometries.
Setting up a new framework to analyse the behavior of squishy particle systems deep in the jammed state is an open and complex task yet to be developed in granular matter research. Some aspects concerning the experimental and numerical study of squishy matter are currently maturing, but major challenges are still ahead. Fortunately, these systems are so ubiquitous that scientists brave enough to take up these challenges can count on the support of a broad and diverse community gathering physicists, geophysicists, biologists and manufacturers, to name a few.