Influence of shaped boundaries on propagating compaction bands in brittle porous media.

The compression of brittle porous media can lead to the propagation of compaction bands. Although such localization phenomena have been observed in different geometries, including cuboidal and axisymmetric uniaxial compression, the role of boundary geometry on compaction features has yet to be explored, despite its relevance in geological conditions and industrial processes. To this end, we investigate the influence of shaped boundaries and inhomogeneous inclusions in a model brittle material made of puffed rice cereal. Using a variety of geometries, we show that compaction bands assume the shape of nearby boundaries, but return to a default planar form a distance away from them. Remarkably, the band aligns parallel to characteristic lines of minor principal stress obtained from a simple linear elastic model. The compelling correlation between the rotation of the principal stress directions and compaction band orientation holds implications for the geological interpretation of localized patterns in rocks and for comprehending the formation of weak planes in pharmaceutical tablets.

Microstructural Origin of Propagating Compaction Patterns in Porous Media.

Porous rocks, foams, cereals, and snow display a diverse set of common compaction patterns, including propagating or stationary bands. Although this commonality across distinct media has been widely noted, the patterns' origin remains debated-current models employ empirical laws for material-specific processes. Here, using a generic model of inelastic structured porous geometries, we show that the previously observed patterns can be attributed to a universal process of pore collapse. Furthermore, the pattern diversity can be mapped in a phase space of only two dimensionless numbers describing material strength and loading rate.

The mechanics of brittle granular materials with coevolving grain size and shape.

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

Measuring grain size fractions of bidisperse granular materials using X-ray radiography.

The mechanical properties of granular materials such as sand, snow and rice are inherently tied to the size of the constituent particles. When a system is composed of particles of various sizes, it is common for these particles to segregate by size when disturbed. There is therefore a need to measure the particle size distribution within granular media as it evolves over time. However, there are very few experimental techniques available which can measure the particle sizes in situ without disturbing the medium. Here we present a technique to determine the volume fractions of the grain sizes in bidisperse granular materials with the aid of dynamic X-ray radiography. As a result of the penetration of the X-rays into the medium, radiography minimises the effect of walls and boundaries on experimental measurements, which typically dominate optical measurements. The technique proposed here is based on using Fourier transforms of X-ray radiographs to extract local measurements evolving over time that can be related to the particle size distribution. For the case of bidisperse granular media, with two distinct particle sizes, we show that this technique can measure the relative concentration of the two species, which we determine via a heuristic calibration parameter. We validate this technique by comparing discrete element simulations of mixtures of known concentration with experimental measurements derived from X-ray radiography of glass beads. In the future, this technique could be used to measure the grain size distribution in systems of bidisperse dense granular media where the concentration of particles is not known a priori. Additionally, the technique can be used to analyse granular segregation as it evolves dynamically.

The mechanism of delayed release in earthquake-induced avalanches.

Snow avalanches can be triggered by strong earthquakes. Most existing models assume that snow slab avalanches happen simultaneously during or immediately after their triggering. Therefore, they cannot explain the plausibility of delayed avalanches that are released minutes to hours after a quake. This paper establishes the basic mechanism of delays in earthquake-induced avalanche release using a novel analytical model that yields dynamics consistent with three documented cases, including two from Western Himalaya and one from central Italy. The mechanism arises from the interplay between creep, strain softening and strain-rate sensitivity of snow, which drive the growth of a basal shear fracture. Our model demonstrates that earthquake-triggered delayed avalanches are rare, yet possible, and could lead to significant damage, especially in long milder slopes. The generality of the model formulation opens a new approach for exploring many other problems related to natural slab avalanche release.

X-ray rheography uncovers planar granular flows despite non-planar walls.

Extremely useful techniques exist to observe the interior of deforming opaque materials, but these methods either require that the sample is replaced with a model material or that the motion is stopped intermittently. For example, X-ray computed tomography cannot measure the continuous flow of materials due to the significant scanning time required for density reconstruction. Here we resolve this technological gap with an alternative X-ray method that does not require such tomographs. Instead our approach uses correlation analysis of successive high-speed radiographs from just three directions to directly reconstruct three-dimensional velocities. When demonstrated on a steady granular system, we discover a compressible flow field that has planar streamlines despite curved confining boundaries, in surprising contrast to Newtonian fluids. More generally, our new X-ray technique can be applied using synchronous source/detector pairs to investigate transient phenomena in various soft matter such as biological tissues, geomaterials and foams.

Tracking time with ricequakes in partially soaked brittle porous media.

When brittle porous media interact with chemically active fluids, they may suddenly crumble. This has reportedly triggered the collapse of rockfill dams, sinkholes, and ice shelves. To study this problem, we use a surrogate experiment for the effect of fluid on rocks and ice involving a column of puffed rice partially soaked in a reservoir of liquid under constant pressure. We disclose localized crushing collapse in the unsaturated region that produces incremental global compaction and loud audible beats. These "ricequakes" repeat perpetually during the experiments and propagate upward through the material. The delay time between consecutive quakes grows linearly with time and is accompanied by creep motion. All those new observations can be explained using a simple chemomechanical model of capillary-driven crushing steps progressing through the micropores.

Evidence That Strain-Rate Softening Is Not Necessary for Material Instability Patterns.

Strain-rate softening has been associated with a wide variety of material instabilities, from the Portevin-Le Chatelier effect in metal alloys to stick-slip motion in crust faults. Dynamic instability patterns have been recently discovered in brittle porous media: diffused, oscillatory, and erratic compaction. Using model simulations inspired by experiments with puffed rice, we question the link between these dynamic patterns and strain-rate sensitivity in such media. An important feature of our model is that it can recover strain-rate softening as an emergent phenomenon, without imposing it a priori at its microstructural scale. More importantly, the model also demonstrates that the full range of dynamic patterns can develop without presenting macroscopic strain-rate softening. Based on this counterexample model, we therefore argue that strain-rate softening should not be taken as a necessary condition for the emergence of instability patterns. Our findings in brittle porous media have implications on models that require strain-rate softening to explain earthquake and metal alloy instabilities.

Dynamic X-ray radiography reveals particle size and shape orientation fields during granular flow.

When granular materials flow, the constituent particles segregate by size and align by shape. The impacts of these changes in fabric on the flow itself are not well understood, and thus novel non-invasive means are needed to observe the interior of the material. Here, we propose a new experimental technique using dynamic X-ray radiography to make such measurements possible. The technique is based on Fourier transformation to extract spatiotemporal fields of internal particle size and shape orientation distributions during flow, in addition to complementary measurements of velocity fields through image correlation. We show X-ray radiography captures the bulk flow properties, in contrast to optical methods which typically measure flow within boundary layers, as these are adjacent to any walls. Our results reveal the rich dynamic alignment of particles with respect to streamlines in the bulk during silo discharge, the understanding of which is critical to preventing destructive instabilities and undesirable clogging. The ideas developed in this paper are directly applicable to many other open questions in granular and soft matter systems, such as the evolution of size and shape distributions in foams and biological materials.

Information thermodynamics of near-equilibrium computation.

In studying fundamental physical limits and properties of computational processes, one is faced with the challenges of interpreting primitive information-processing functions through well-defined information-theoretic as well as thermodynamic quantities. In particular, transfer entropy, characterizing the function of computational transmission and its predictability, is known to peak near critical regimes. We focus on a thermodynamic interpretation of transfer entropy aiming to explain the underlying critical behavior by associating information flows intrinsic to computational transmission with particular physical fluxes. Specifically, in isothermal systems near thermodynamic equilibrium, the gradient of the average transfer entropy is shown to be dynamically related to Fisher information and the curvature of system's entropy. This relationship explicitly connects the predictability, sensitivity, and uncertainty of computational processes intrinsic to complex systems and allows us to consider thermodynamic interpretations of several important extreme cases and trade-offs.

Complex networks in confined comminution.

The physical process of confined comminution is investigated within the framework of complex networks. We first characterize the topology of the unweighted contact networks as generated by the confined comminution process. We find this process gives rise to an ultimate contact network which exhibits a scale-free degree distribution and small world properties. In particular, if viewed in the context of networks through which information travels along shortest paths, we find that the global average of the node vulnerability decreases as the comminution process continues, with individual node vulnerability correlating with grain size. A possible application to the design of synthetic networks (e.g., sensor networks) is highlighted. Next we turn our attention to the physics of the granular comminution process and examine force transmission with respect to the weighted contact networks, where each link is weighted by the inverse magnitude of the normal force acting at the associated contact. We find that the strong forces (i.e., force chains) are transmitted along pathways in the network which are mainly following shortest-path routing protocols, as typically found, for example, in communication systems. Motivated by our earlier studies of the building blocks for self-organization in dense granular systems, we also explore the properties of the minimal contact cycles. The distribution of the contact strain energy intensity of 4-cycle motifs in the ultimate state of the confined comminution process is shown to be consistent with a scale-free distribution with infinite variance, thereby suggesting that 4-cycle arrangements of grains are capable of storing vast amounts of energy in their contacts without breaking.

The role of self-organization during confined comminution of granular materials.

During confined comminution of granular materials a power-law grain size distribution (gsd) frequently evolves. We consider this power law as a hint for fractal topology if self-similar patterns appear across the scales. We demonstrate that this ultimate topology is mostly affected by the rules that define the self-organization of the fragment subunits, which agrees well with observations from simplistic models of cellular automata. There is, however, a major difference that highlights the novelty of the current work: here the conclusion is based on a comprehensive study using two-dimensional 'crushable' discrete-element simulations that do not neglect physical conservation laws. Motivated by the paradigm of self-organized criticality, we further demonstrate that in uniaxial compression the emerging ultimate fractal topology, as given by the fractal dimension, is generally insensitive to alteration of global index properties of initial porosity and initial gsd. Finally, we show that the fractal dimension in the confined crushing systems is approached irrespective of alteration of the criteria that define when particles crush.

Soil mechanics: breaking ground.

In soil mechanics, student's models are classified as simple models that teach us unexplained elements of behaviour; an example is the Cam clay constitutive models of critical state soil mechanics (CSSM). 'Engineer's models' are models that elaborate the theory to fit more behavioural trends; this is usually done by adding fitting parameters to the student's models. Can currently unexplained behavioural trends of soil be explained without adding fitting parameters to CSSM models, by developing alternative student's models based on modern theories?Here I apply an alternative theory to CSSM, called 'breakage mechanics', and develop a simple student's model for sand. Its unique and distinctive feature is the use of an energy balance equation that connects grain size reduction to consumption of energy, which enables us to predict how grain size distribution (gsd) evolves-an unprecedented capability in constitutive modelling. With only four parameters, the model is physically clarifying what CSSM cannot for sand: the dependency of yielding and critical state on the initial gsd and void ratio.