Search

Thermal fluctuations, geometric exclusion, and external driving all govern the mechanical response of dense particulate suspensions. Here, we measure the stress-strain response of quasi-two-dimensional flow-stabilized microsphere heaps in a regime in which all three effects are present using a microfluidic device. We observe that the elastic modulus and the mean interparticle separation of the heaps are tunable via the confining stress provided by the fluid flow. Furthermore, the measured stress-strain curves exhibit a universal nonlinear shape, which can be predicted from a thermal van der Waals equation of state with excluded volume. This analysis indicates that many-body interactions contribute a significant fraction of the stress supported by the heap. [ABSTRACT FROM AUTHOR]

A liquid droplet resting on a soft gel substrate can deform that substrate to the point of material failure, whereby fractures develop on the gel surface that propagate outwards from the contact line in a starburst pattern. In this paper, we characterize (i) the initiation process, in which the number of arms in the starburst is controlled by the ratio of the surface tension contrast to the gel's elastic modulus, and (ii) the propagation dynamics showing that once fractures are initiated they propagate with a universal power law L ∝ t3/4. We develop a model for crack initiation by treating the gel as a linear elastic solid and computing the deformations within the substrate from the liquid-solid wetting forces. The elastic solution shows that both the location and the magnitude of the wetting forces are critical in providing a quantitative prediction for the number of fractures and, hence, an interpretation of the initiation of capillary fractures. This solution also reveals that the depth of the gel is an important factor in the fracture process, as it can help mitigate large surface tractions; this finding is confirmed with experiments. We then develop a model for crack propagation by considering the transport of an inviscid fluid into the fracture tip of an incompressible material and find that a simple energy-conservation argument can explain the observed material-independent power law. We compare predictions for both linear elastic and neo-Hookean solids, finding that the latter better explains the observed exponent. [ABSTRACT FROM AUTHOR]

Although jammed granular systems are athermal, several thermodynamiclike descriptions have been proposed which make quantitative predictions about the distribution of volume and stress within a system and provide a corresponding temperaturelike variable. We perform experiments with an apparatus designed to generate a large number of independent, jammed, two-dimensional configurations. Each configuration consists of a single layer of photoelastic disks supported by a gentle layer of air. New configurations are generated by cyclically dilating, mixing, and then recompacting the system through a series of boundary displacements. Within each configuration, a bath of particles surrounds a smaller subsystem of particles with a different interparticle friction coefficient than the bath. The use of photoelastic particles permits us to find all particle positions as well as the vector forces at each interparticle contact. By comparing the temperaturelike quantities in both systems, we find compactivity (conjugate to the volume) does not equilibrate between the systems, while the angoricity (conjugate to the stress) does. Both independent components of the angoricity are linearly dependent on the hydrostatic pressure, in agreement with predictions of the stress ensemble. [ABSTRACT FROM AUTHOR]

Experiments quantifying the rotational and translational motion of particles in a dense, driven, 2D granular gas floating on an air table reveal that kinetic energy is divided equally between the two translational and one rotational degrees of freedom. This equipartition persists when the particle properties, confining pressure, packing density, or spatial ordering are changed. While the translational velocity distributions are the same for both large and small particles, the angular velocity distributions scale with the particle radius. The probability distributions of all particle velocities have approximately exponential tails. Additionally, we find that the system can be described with a granular Boyle's law with a van der Waals-like equation of state. These results demonstrate ways in which conventional statistical mechanics can unexpectedly apply to nonequilibrium systems. [ABSTRACT FROM AUTHOR]

Granular materials can fail through spontaneous events like earthquakes or brittle fracture. However, measurements and analytic models which forecast failure in this class of materials, while of both fundamental and practical interest, remain elusive. Materials including numerical packings of spheres, colloidal glasses, and granular materials have been known to develop an excess of low-frequency vibrational modes as the confining pressure is reduced. Here, we report experiments on sheared granular materials in which we monitor the evolving density of excited modes via passive monitoring of acoustic emissions. We observe a broadening of the distribution of excited modes coincident with both bulk and local plasticity, and evolution in the shape of the distribution before and after bulk failure. These results provide a new interpretation of the changing state of the material on its approach to stick-slip failure. [ABSTRACT FROM AUTHOR]

Granular media, whose features range from the particle scale to the force-chain scale and the bulk scale, are usually modeled as either particulate or continuum materials. In contrast with each of these approaches, network representations are natural for the simultaneous examination of microscopic, mesoscopic, and macroscopic features. In this paper, we treat granular materials as spatially embedded networks in which the nodes (particles) are connected by weighted edges obtained from contact forces. We test a variety of network measures to determine their utility in helping to describe sound propagation in granular networks and find that network diagnostics can be used to probe particle-, curve-, domain-, and system-scale structures in granular media. In particular, diagnostics of mesoscale network structure are reproducible across experiments, are correlated with sound propagation in this medium, and can be used to identify potentially interesting size scales. We also demonstrate that the sensitivity of network diagnostics depends on the phase of sound propagation. In the injection phase, the signal propagates systemically, as indicated by correlations with the network diagnostic of global efficiency. In the scattering phase, however, the signal is better predicted by mesoscale community structure, suggesting that the acoustic signal scatters over local geographic neighborhoods. Collectively, our results demonstrate how the force network of a granular system is imprinted on transmitted waves [ABSTRACT FROM AUTHOR]

Fluctuations of the local volume fraction within granular materials have previously been observed to decrease as the system approaches jamming. We experimentally examine the role of boundary conditions and interparticle friction μ on this relationship for a dense granular material of bidisperse particles driven under either constant volume or constant pressure. Using a radical Voronoï tessellation, we find the variance of the local volume fraction ∅ monotonically decreases as the system becomes more dense, independent of boundary condition and μ. We examine the universality and origins of this trend using experiments and the recent granocentric model [M. Clusel, E. I. Corwin, A. O. N. Siemens, and J. Brujić, Nature (London) 460, 611(2009); E. I. Corwin, M. Clusel, A. O. N. Siemens, and J. Brujić, Soft Matter 6, 2949 (2010)], modified to draw particle locations from an arbitrary distribution P(s) of neighbor distances s. The mean and variance of the observed 2(s) are described by a single length scale controlled by∅-. Through the granocentric model, we observe that diverse functional forms of P(s) all produce the trend of decreasing fluctuations, but only the experimentally observed P(s) provides quantitative agreement with the measured ∅ fluctuations. Thus, we find that both P(s) and P(∅) encode similar information about the ensemble of observed packings and are connected to each other by the local granocentric model. [ABSTRACT FROM AUTHOR]

We show how rigidity emerges in experiments on sheared two-dimensional frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a critical contact number zc=2.4±0.1, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a spongelike morphology. While the pressure within rigid structures always exceeds the pressure outside the rigid structures, they are not identified with the force chains of shear jamming. These findings highlight the need to focus on mechanical stability arising through arch structures and hinges at the mesoscale. [ABSTRACT FROM AUTHOR]

We report on a series of experiments in which a grain-sized intruder is pushed by a spring through a two-dimensional granular material composed of photoelastic disks in a Couette geometry. We study the intruder dynamics as a function of packing fraction for two types of supporting substrates: A frictional glass plate and a layer of water for which basal friction forces are negligible. We observe two dynamical regimes: Intermittent flow, in which the intruder moves freely most of the time but occasionally gets stuck, and stick-slip dynamics, in which the intruder advances via a sequence of distinct, rapid events. When basal friction is present, we observe a smooth crossover between the two regimes as a function of packing fraction, and we find that reducing the interparticle friction coefficient causes the stick-slip regime to shift to higher packing fractions. When basal friction is eliminated, we observe intermittent flow at all accessible packing fractions. For all cases, we present results for the statistics of stick events, the intruder velocity, and the force exerted on the intruder by the grains. Our results indicate the qualitative importance of basal friction at high packing fractions and suggest a possible connection between intruder dynamics in a static material and clogging dynamics in granular flows. [ABSTRACT FROM AUTHOR]

Stress-based ensembles incorporating temperaturelike variables have been proposed as a route to an equation of state for granular materials. To test the efficacy of this approach, we perform experiments on a two-dimensional photoelastic granular system under three loading conditions: uniaxial compression, biaxial compression, and simple shear. From the interparticle forces, we find that the distributions of the normal component of the coarse-grained force-moment tensor are exponential tailed, while the deviatoric component is Gaussian distributed. This implies that the correct stress-based statistical mechanics conserves both the force-moment tensor and the Maxwell-Cremona force-tiling area. As such, two variables of state arise: the tensorial angoricity (α̂) and a new temperaturelike quantity associated with the force-tile area which we name keramicity (κ). Each quantity is observed to be inversely proportional to the global confining pressure; however, only κ exhibits the protocol independence expected of a state variable, while α̂ behaves as a variable of process. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We identify and characterize a new class of fingering instabilities in liquid metals; these instabilities are unexpected due to the large interfacial tension of metals. Electrochemical oxidation lowers the effective interfacial tension of a gallium-based liquid metal alloy to values approaching zero, thereby inducing drastic shape changes, including the formation of fractals. The measured fractal dimension (D=1.3±0.05) places the instability in a different universality class than other fingering instabilities. By characterizing changes in morphology and dynamics as a function of droplet volume and applied electric potential, we identify the three main forces involved in this process: interfacial tension, gravity, and oxidative stress. Importantly, we find that electrochemical oxidation can generate compressive interfacial forces that oppose the tensile forces at a liquid interface. The surface oxide layer ultimately provides a physical and electrochemical barrier that halts the instabilities at larger positive potentials. Controlling the competition between interfacial tension and oxidative (compressive) stresses at the interface is important for the development of reconfigurable electronic, electromagnetic, and optical devices that take advantage of the metallic properties of liquid metals. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

As a result of extreme weather conditions, such as heavy precipitation, natural hillslopes can fail dramatically; these slope failures can occur on a dry day, due to time lags between rainfall and pore-water pressure change at depth, or even after days to years of slow motion. While the prefailure deformation is sometimes apparent in retrospect, it remains challenging to predict the sudden transition from gradual deformation (creep) to runaway failure. We use a network science method-multilayer modularity optimization-to investigate the spatiotemporal patterns of deformation in a region near the 2017 Mud Creek, California landslide. We transform satellite radar data from the study site into a spatially embedded network in which the nodes are patches of ground and the edges connect the nearest neighbors, with a series of layers representing consecutive transits of the satellite. Each edge is weighted by the product of the local slope (susceptibility to failure) measured from a digital elevation model and ground surface deformation (current rheological state) from interferometric synthetic aperture radar (InSAR). We use multilayer modularity optimization to identify strongly connected clusters of nodes (communities) and are able to identify both the location of Mud Creek and nearby creeping landslides which have not yet failed. We develop a metric, i.e., community persistence, to quantify patterns of ground deformation leading up to failure, and find that this metric increased from a baseline value in the weeks leading up to Mud Creek's failure. These methods hold promise as a technique for highlighting regions at risk of catastrophic failure.

We report experimental evidence of a Gardner-like crossover from variable to persistent force contacts in a two-dimensional bidisperse granular crystal by analyzing the variability of both particle positions and force networks formed under uniaxial compression. Starting from densities just above the freezing transition and for variable amounts of additional compression, we compare configurations to both their own initial state and to an ensemble of equivalent reinitialized states. This protocol shows that force contacts are largely undetermined when the density is below a Gardner-like crossover, after which they gradually transition to being persistent, being fully so only above the jamming point. We associate the disorder that underlies this effect with the size of the microscopic asperities of the photoelastic disks used, by analogy to other mechanisms that have been previously predicted theoretically.

We show how rigidity emerges in experiments on sheared two-dimensional frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a critical contact number z_{c}=2.4±0.1, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a spongelike morphology. While the pressure within rigid structures always exceeds the pressure outside the rigid structures, they are not identified with the force chains of shear jamming. These findings highlight the need to focus on mechanical stability arising through arch structures and hinges at the mesoscale.

We investigate the origins of the widely observed "en passant" crack pattern, which forms through interactions between two approaching cracks. A rectangular elastic plate is notched on each long side and then subjected to quasistatic uniaxial strain from the short side. The two cracks propagate along approximately straight paths until they pass each other, after which they curve and release a lenticular fragment. We find that, for materials with diverse mechanical properties, the shape of this fragment has an aspect ratio of 2:1, with the length scale set by the initial crack offset s and the time scale set by the ratio of s to the pulling velocity. The cracks have a universal square root shape, which we understand by using a simple geometric model of the crack-crack interaction.

Granular materials segregate by size under shear, and the ability to quantitatively predict the time required to achieve complete segregation is a key test of our understanding of the segregation process. In this paper, we apply the Gray-Thornton model of segregation (developed for linear shear profiles) to a granular flow with an exponential shear profile, and evaluate its ability to describe the observed segregation dynamics. Our experiment is conducted in an annular Couette cell with a moving lower boundary. The granular material is initially prepared in an unstable configuration with a layer of small particles above a layer of large particles. Under shear, the sample mixes and then resegregates so that the large particles are located in the top half of the system in the final state. During this segregation process, we measure the velocity profile and use the resulting exponential fit as input parameters to the model. To make a direct comparison between the continuum model and the observed segregation dynamics, we map the local concentration (from the model) to changes in packing fraction; this provides a way to make a semiquantitative comparison with the measured global dilation. We observe that the resulting model successfully captures the presence of a fast mixing process and relatively slower resegregation process, but the model predicts a finite resegregation time, while in the experiment resegregation occurs only exponentially in time.

The vertical size segregation of granular materials, a process commonly associated with the Brazil-nut effect, has generally been thought to proceed faster the greater the size difference of the particles. We experimentally investigate sheared dense bidisperse granular materials as a function of the size ratio of the two species and find that the mixing rate at low confining pressure behaves as expected from percolation-based arguments. However, we also observe an anomalous effect for the resegregation rates, wherein the segregation rate is a nonmonotonic function of the particle size ratio with a maximum for intermediate particle size ratio. Combined with the fact that increasing the confining pressure significantly suppresses both mixing and segregation rates of particles of sufficiently dissimilar size, we propose that the anomalous behavior may be attributed to a species-dependent distribution of forces within the system.

We report a novel surface-tension driven instability observed for droplets spreading on a compliant substrate. When a droplet is released on the surface of an agar gel, it forms arms or cracks when the ratio of surface-tension gradient to gel strength is sufficiently large. We explore a range of gel strengths and droplet surface tensions and find that the onset of the instability and the number of arms depend on the ratio of surface tension to gel strength. However, the arm length grows with an apparently universal law L proportional t(3/4).

Experiments on spherical particles in a 3D annular shear cell vibrated from below and sheared from above show a hysteretic freezing or melting transition. Under sufficient vibration a crystallized state is observed, which can be melted by sufficient shear. The critical line for this transition coincides with equal kinetic energies for vibration and shear. The force distribution is double peaked in the crystalline state and single peaked with an approximately exponential tail in the disordered state. Continuous relations between pressure and volume (with dP/dV>0) exist for a continuum of partially and/or intermittently melted states over a range of parameters.

We investigate a novel bursting state in inclined layer thermal convection in which convection rolls exhibit intermittent, localized, transverse bursts. With increasing temperature difference, the bursts increase in duration and number while exhibiting a characteristic wave number, magnitude, and size. We propose a mechanism which describes the duration of the observed bursting intervals and compare our results to bursting processes in other systems.

We report experimental results on the defect turbulent state of undulation chaos in inclined layer convection of a fluid with Prandtl number approximately 1. By measuring defect density and undulation wave number, we find that the onset of undulation chaos coincides with the theoretically predicted onset for stable, stationary undulations. At stronger driving, we observe a competition between ordered undulations and undulation chaos, suggesting bistability between a fixed-point attractor and spatiotemporal chaos. In the defect turbulent regime, we measured the defect creation, annihilation, entering, leaving, and rates. We derive a universal probability distribution function which agrees with the experimental findings.