Search

The flow of granular materials through constricted openings is important in many natural and industrial processes. These complex flows - featuring dense, dissipative flow in the bulk but low-dissipation, low density outpouring in the vicinity of the orifice - have long been characterized empirically by the Beverloo rule and, recently, modeled successfully using energy balance. The dependence of flow rate on the silo's angle with respect to gravity, however, is not captured by current models. We experimentally investigate the role of tilt angle in this work using a quasi-2D monolayer of grains in a silo. We measure mass flow rate, the average exit velocities of grains, and the packing fraction along the orifice with varying tilt angles. We propose a model that describes our results (and earlier findings with 3D systems [H. G. Sheldon and D. J. Durian, Granul. Matter 12, 579 (2010)]) by considering the dependence of outpouring speed and angle with respect to the orifice angle and, importantly, the angle of stagnant zones adjacent to the orifice. We conclude by posing questions about possible extensions of our model in order to describe spatial variations of exit velocity and density along the orifice cross section., 10 pages and 11 figures

We test the standard model for the length contraction of a bundle of strings under twist, and find deviation that is significantly greater than typically appreciated and that has a different nature at medium and large twist angles. By including volume conservation, we achieve better fits to data for single-, double-, and triple-stranded bundles of Nylon monofilament as an ideal test case. This gives a well-defined procedure for extracting an effective twist radius that characterizes contraction behavior. While our approach accounts for the observed faster-than-expected contraction up to medium twist angles, we also find that the contraction is nevertheless slower than expected at large twist angles for both Nylon monofilament bundles and several other string types. The size of this effect varies with the individual-string braid structure and with the number of strings in the bundle. We speculate that it may be related to elastic deformation within the material. However, our first modeling attempt does not fully capture the observed behavior., 8 pages, 8 figures

Physical networks, such as biological neural networks, can learn desired functions without a central processor, using local learning rules in space and time to learn in a fully distributed manner. Learning approaches such as equilibrium propagation, directed aging, and coupled learning similarly exploit local rules to accomplish learning in physical networks such as mechanical, flow, or electrical networks. In contrast to certain natural neural networks, however, such approaches have so far been restricted to the quasistatic limit, where they learn on time scales slow compared to their physical relaxation. This quasistatic constraint slows down learning, limiting the use of these methods as machine learning algorithms, and potentially restricting physical networks that could be used as learning platforms. Here we explore learning in an electrical resistor network that implements coupled learning, both in the lab and on the computer, at rates that range from slow to far above the quasistatic limit. We find that up to a critical threshold in the ratio of the learning rate to the physical rate of relaxation, learning speeds up without much change of behavior or error. Beyond the critical threshold, the error exhibits oscillatory dynamics but the networks still learn successfully., 26 pages, 5 figures

We present a solvable model inspired by dimensional analysis for the time-dependent spreading of droplets that partially wet a substrate, where the spreading eventually stops and the contact angle reaches a nonzero equilibrium value. We separately consider small droplets driven by capillarity and large droplets driven by gravity. To explore both regimes, we first measure the equilibrium radius versus a comprehensive range of droplet volumes for four household fluids, and we compare the results with predictions based on minimizing the sum of gravitational and interfacial energies. The agreement is good, and gives a reliable measurement of an equilibrium contact angle that is consistent in both small and large droplet regimes. Next we use energy considerations to develop equations of motion for the time dependence of the spreading, in both regimes, where the driving forces are balanced against viscous drag in the bulk of the droplet and by friction at the moving contact line. Our approach leads to explicit prediction of the functional form of the spreading dynamics. It successfully describes prior data for a small capillary-driven droplet, and it fits well to new data we obtain for large gravity-driven droplets with a wide range of volumes.

We investigate the local- and long-range structure of four different space-filling cellular patterns: bubbles in a quasi-2d foam plus Voronoi constructions made around points that are uncorrelated (Poisson patterns), low discrepancy (Halton patterns), and displaced from a lattice by Gaussian noise (Einstein patterns). We study distributions of local quantities including cell areas and topological features; the former is the widest for bubbles in a foam making them locally the most disordered but the latter show no significant differences between the cellular patterns. Long-range structure is probed by the spectral density and also by converting the real-space spectrum of number density or volume fraction fluctuations for windows of diameter $D$ to the effective distance $h(D)$ from the window boundary where these fluctuations occur. This real-space hyperuniformity disorder length spectroscopy is performed on various point patterns which are determined by the centroids of the bubbles in the foam, by the points patterns around which the Voronoi cells are created and by the centroids of the Voronoi cells. These patterns are either unweighted or weighted by the area of the cells they occupy. The unweighted bubble centroids have $h(D)$ that collapses for the different ages of of the foam with random Poissonian fluctuations at long distances. All patterns of area-weighted points have constant $h(D)=h_e$ for large $D$; $h_e=0.084 \sqrt{\left< a\right>}$ for the bubble centroids is the smallest value, meaning they are most uniform. All the weighted centroids collapse to the same constant $h_e=0.084 \sqrt{\left< a\right>}$ as for the foams. A similar analysis is performed on the edges of the cells where the spectra of $h(D)$ for the foam edges show $h(D) \sim D^{1-��}$ where $��=0.3$, The paper is 14 pages, 10 figures

The distribution of string tension on the contact line between an ideal string and a massive pulley is a frequently-discussed but incompletely-posed problem that confronts students in introductory mechanics. We highlight ambiguities in the usual presentation of this problem by the massive Atwood's machine and discuss two compact resolutions that treat situations where the pulley or the string elastically deform. We propose experiments that can be developed in an intermediate laboratory to determine the tension profile., Comment: American Journal of Physics (in press): pdfLaTeX: 8 pages, 6 figures

Quasi-static tensile experiments were performed for a model disordered solid consisting of a two-dimensional raft of polydisperse floating granular particles with capillary attractions. The ductility is tuned by controlling the capillary interaction range, which varies with the particle size. During the tensile tests, after an initial period of elastic deformation, strain localization occurs and leads to the formation of a shear band at which the pillar later fails. In this process, small particles with long-ranged interactions can endure large plastic deformations without forming significant voids, while large particles with short-range interactions fail dramatically by fracturing at small deformation. Particle-level structure was measured, and the strain-localized region was found to have higher structural anisotropy than the bulk. Local interactions between anisotropic sites and particle rearrangements were the main mechanisms driving strain localization and the subsequent failure, and significant differences of such interactions exist between ductile and brittle behaviors.

We investigate the coupling between the interstitial medium and granular particles by studying the hopper flow of dry and submerged systems experimentally and numerically. In accordance with earlier studies, we find that the dry hopper empties at a constant rate. However, in the submerged system we observe the surging of the flow rate. We model both systems using the discrete element method, which we couple with computational fluid dynamics in the case of a submerged hopper. We are able to match the simulations and the experiments with good accuracy by fitting the particle–particle contact friction for each system separately. Submerging the hopper changes the particle–particle contact friction from µvacuum = 0.15 to µsub = 0.13, while all the other simulation parameters remain the same.

Combining new and old data, we find that particle settling speeds are well-described by the empirical Richardson-Zaki function up to jamming, but surprisingly, there are separate branches for Brownian and non-Brownian particles with a crossover at an extraordinarily large Peclet number of ${10}^{8}$.

Structural defects within amorphous packings of symmetric particles can be characterized using a machine learning approach that incorporates structure functions of radial distances and angular arrangement. This yields a scalar field, \emph{softness}, that correlates with the probability that a particle is about to rearrange. However, when particle shapes are elongated, as in the case of dimers and ellipses, we find the standard structure functions produce imprecise softness measurements. Moreover, ellipses exhibit deformation profiles in stark contrast to circular particles. In order to account for effects of orientation and alignment, we introduce new structure functions to recover predictive performance of softness, as well as provide physical insight to local and extended dynamics. We study a model disordered solid, a bidisperse two-dimensional granular pillar, driven by uniaxial compression and composed entirely of monomers, dimers, or ellipses. We demonstrate how the computation of softness via support vector machine extends to dimers and ellipses with the introduction of new orientational structure functions. Then, we highlight the spatial extent of rearrangements and defects, as well as their cross-correlation, for each particle shape. Finally, we demonstrate how an additional machine learning algorithm, recursive feature elimination, provides an avenue to better understand how softness arises from particular structural aspects. We identify the most crucial structure functions in determining softness and discuss their physical implications., 16 pages, 12 figures

A new experimental facility has been designed and constructed to study driven granular media in a low-gravity environment. This versatile instrument, fully automatized, with a modular design based on several interchangeable experimental cells, allows us to investigate research topics ranging from dilute to dense regimes of granular media such as granular gas, segregation, convection, sound propagation, jamming, and rheology—all without the disturbance by gravitational stresses active on Earth. Here, we present the main parameters, protocols, and performance characteristics of the instrument. The current scientific objectives are then briefly described and, as a proof of concept, some first selected results obtained in low gravity during parabolic flight campaigns are presented.A new experimental facility has been designed and constructed to study driven granular media in a low-gravity environment. This versatile instrument, fully automatized, with a modular design based on several interchangeable experimental cells, allows us to investigate research topics ranging from dilute to dense regimes of granular media such as granular gas, segregation, convection, sound propagation, jamming, and rheology—all without the disturbance by gravitational stresses active on Earth. Here, we present the main parameters, protocols, and performance characteristics of the instrument. The current scientific objectives are then briefly described and, as a proof of concept, some first selected results obtained in low gravity during parabolic flight campaigns are presented.

When a large blunt object penetrates into a granular material, the force response exhibits an initial nonlinear relationship with depth that is widely attributed to the formation of a shear jammed stagnant zone. We present an experimental study of a model 2D granular system that preserves common aspects of slow granular penetration, such as a linear force law for large depths, while also allowing for direct and precise measurement of a stagnant zone. We also probe structural anisotropies throughout the bed and examine their characteristics inside and outside the stagnant zone. Finally, we simultaneously consider effects of concave elongated particle shapes (dimers) and find that while most results are not affected, the relationship between local structure and strain suggests a distinct fragility in the stagnant zone.

Diffusing wave spectroscopy (DWS), which makes use of multiple scattered light in concentrated and highly turbid colloidal systems,has emerged as a very reliable and promising technique to probe their dynamics and rheology. However, the application of DWS to the multiple emulsion systems, particularly double emulsions (DEs),has been relatively limited so far despite their significant potential in the diverse fields of applications including food, pharmaceuticals,and cosmetics. The potential of DSW could facilitate gaining a more detailed understanding and quantitative measurement of the factors that govern the dynamics, structure, and properties of DEs,and most importantly their stability. This review briefly describes fundamental principles of DWS and its potential application to DEs to study the evolution of dynamics and structure at various time and length scales, as well as the rheological properties. Refereed/Peer-reviewed

We investigate the short-, medium-, and long-range structure of soft-disk configurations for a wide range of area fractions and simulation protocols by converting the real-space spectrum of volume fraction fluctuations for windows of width $L$ to the distance $h(L)$ from the window boundary over which fluctuations occur. Rapidly quenched unjammed configurations exhibit size-dependent super-Poissonian long-range features that surprisingly approach the totally random limit even close to jamming. Above and just below jamming, the spectra exhibit a plateau $h(L)={h}_{e}$ for $L$ larger than particle size and smaller than a cutoff ${L}_{c}$ beyond which there are long-range fluctuations. The value of ${h}_{e}$ is independent of protocol and characterizes the putative hyperuniform limit. This behavior is compared with that for Einstein solids, with and without hyperuniformity-destroying defects. We find that key structural features of the particle configurations are more evident, as well as easier and more intuitive to quantify, using the real-space spectrum of hyperuniformity lengths rather than the spectral density.

Deformation of a fluid interface caused by the presence of objects at the interface can lead to large lateral forces between objects. We explore these fluid-mediated attractive force between partially submerged vertical cylinders. Forces are experimentally measured by slowly separating cylinder pairs and cylinder triplets after capillary rise is initially established for cylinders in contact. For cylinder pairs, numerical computations and a theoretical model are found to be in good agreement with measurements. The model provides insight into the relative importance of the contributions to the total force. For small separations, the lateral force is dominated by the fluid pressure acting over the wetted cylinder surfaces. At large separations, the surface tension acting along the contact line dominates the lateral force. A crossover between the two regimes occurs at a separation of around half of a capillary length. The experimentally measured forces between cylinder triplets are also in good agreement with numerical computations, and we show that pairwise contributions account for nearly all of the attractive force between triplets. For cylinders with an equilibrium capillary rise height greater than the height of the cylinder, we find that the attractive force depends on the height of the cylinders above the submersion level, which provides a means to create precisely controlled tunable cohesive forces between objects deforming a fluid interface.

We report two techniques to mitigate stripe artifacts in light-sheet fluorescence imaging. The first uses an image processing algorithm called the multidirectional stripe remover (MDSR) method to filter stripes from an existing image. The second uses an elliptical holographic diffuser (EHD) with strong scattering anisotropy to prevent stripe formation during image acquisition. These techniques facilitate accurate interpretation of image data, especially in denser samples. They are also facile and cost-effective., Supplementary Material available at the journal page

We probe the effects of particle shape on the global and local behavior of a two-dimensional granular pillar, acting as a proxy for a disordered solid, under uniaxial compression. This geometry allows for direct measurement of global material response, as well as tracking of all individual particle trajectories. In general, drawing connections between local structure and local dynamics can be challenging in amorphous materials due to lower precision of atomic positions, so this study aims to elucidate such connections. We vary local interactions by using three different particle shapes: discrete circular grains (monomers), pairs of grains bonded together (dimers), and groups of three bonded in a triangle (trimers). We find that dimers substantially strengthen the pillar and the degree of this effect is determined by orientational order in the initial condition. In addition, while the three particle shapes form void regions at distinct rates, we find that anisotropies in the local amorphous structure remain robust through the definition of a metric that quantifies packing anisotropy. Finally, we highlight connections between local deformation rates and local structure., 15 pages, 15 figures

We introduce the concept of a ``hyperuniformity disorder length'' $h$ that controls the variance of volume fraction fluctuations for randomly placed windows of fixed size. In particular, fluctuations are determined by the average number of particles within a distance $h$ from the boundary of the window. We first compute special expectations and bounds in $d$ dimensions, and then illustrate the range of behavior of $h$ versus window size $L$ by analyzing several different types of simulated two-dimensional pixel patterns---where particle positions are stored as a binary digital image in which pixels have value zero if empty and one if they contain a particle. The first are random binomial patterns, where pixels are randomly flipped from zero to one with probability equal to area fraction. These have long-ranged density fluctuations, and simulations confirm the exact result $h=L/2$. Next we consider vacancy patterns, where a fraction $f$ of particles on a lattice are randomly removed. These also display long-range density fluctuations, but with $h=(L/2)(f/d)$ for small $f$, and $h=L/2$ for $f\ensuremath{\rightarrow}1$. And finally, for a hyperuniform system with no long-range density fluctuations, we consider ``Einstein patterns,'' where each particle is independently displaced from a lattice site by a Gaussian-distributed amount. For these, at large $L,\phantom{\rule{0.16em}{0ex}}h$ approaches a constant equal to about half the root-mean-square displacement in each dimension. Then we turn to gray-scale pixel patterns that represent simulated arrangements of polydisperse particles, where the volume of a particle is encoded in the value of its central pixel. And we discuss the continuum limit of point patterns, where pixel size vanishes. In general, we thus propose to quantify particle configurations not just by the scaling of the density fluctuation spectrum but rather by the real-space spectrum of $h(L)$ versus $L$. We call this approach ``hyperuniformity disorder length spectroscopy''.

For dry foams, the transport of gas from small high-pressure bubbles to large low-pressure bubbles is dominated by diffusion across the thin soap films separating neighboring bubbles. For wetter foams, the film areas become smaller as the Plateau borders and vertices inflate with liquid. So-called "border-blocking" models can explain some features of wet-foam coarsening based on the presumption that the inflated borders totally block the gas flux; however, this approximation dramatically fails in the wet or unjamming limit where the bubbles become close-packed spheres and coarsening proceeds even though there are no films. Here, we account for the ever-present border-crossing flux by a new length scale defined by the average gradient of gas concentration inside the borders. We compute that it is proportional to the geometric average of film and border thicknesses, and we verify this scaling by numerical solution of the diffusion equation. We similarly consider transport across inflated vertices and surface Plateau borders in quasi-two-dimensional foams. And we show how the dA/dt=K_{0}(n-6) von Neumann law is modified by the appearance of terms that depend on bubble size and shape as well as the concentration gradient length scales. Finally, we use the modified von Neumann law to compute the growth rate of the average bubble area, which is not constant.

We report on the nature of flow events for the gravity-driven discharge of glass beads through a hole that is small enough that the hopper is susceptible to clogging. In particular, we measure the average and standard deviation of the distribution of discharged masses as a function of both hole and grain sizes. We do so in air, which is usual, but also with the system entirely submerged under water. This damps the grain dynamics and could be expected to dramatically affect the distribution of the flow events, which are described in prior work as avalanche-like. Though the flow is slower and the events last longer, we find that the average discharge mass is only slightly reduced for submerged grains. Furthermore, we find that the shape of the distribution remains exponential, implying that clogging is still a Poisson process even for immersed grains. Per Thomas and Durian [Phys. Rev. Lett. 114, 178001 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.178001], this allows for an interpretation of the average discharge mass in terms of the fraction of flow microstates that precede, i.e., that effectively cause, a stable clog to form. Since this fraction is barely altered by water, we conclude that the crucial microscopic variables are the grain positions; grain momenta play only a secondary role in destabilizing weak incipient arches. These insights should aid ongoing efforts to understand the susceptibility of granular hoppers to clogging.

The concept of a hyperuniformity disorder length h was recently introduced for analyzing volume fraction fluctuations for a set of measuring windows [Chieco et al., Phys. Rev. E 96, 032909 (2017).PLEEE81539-375510.1103/PhysRevE.96.032909]. This length permits a direct connection to the nature of disorder in the spatial configuration of the particles and provides a way to diagnose the degree of hyperuniformity in terms of the scaling of h and its value in comparison with established bounds. Here, this approach is generalized for extended particles, which are larger than the image resolution and can lie partially inside and partially outside the measuring windows. The starting point is an expression for the relative volume fraction variance in terms of four distinct volumes: that of the particle, the measuring window, the mean-squared overlap between particle and region, and the region over which particles have nonzero overlap with the measuring window. After establishing limiting behaviors for the relative variance, computational methods are developed for both continuum and pixelated particles. Exact results are presented for particles of special shape and for measuring windows of special shape, for which the equations are tractable. Comparison is made for other particle shapes, using simulated Poisson patterns. And the effects of polydispersity and image errors are discussed. For small measuring windows, both particle shape and spatial arrangement affect the form of the variance. For large regions, the variance scaling depends only on arrangement but particle shape sets the numerical proportionality. The combined understanding permit the measured variance to be translated to the spectrum of hyperuniformity lengths versus region size, as the quantifier of spatial arrangement. This program is demonstrated for a system of nonoverlapping particles at a series of increasing packing fractions as well as for an Einstein pattern of particles with several different extended shapes.

When deformed beyond their elastic limits, crystalline solids flow plastically via particle rearrangements localized around structural defects. Disordered solids also flow, but without obvious structural defects. We link structure to plasticity in disordered solids via a microscopic structural quantity, "softness," designed by machine learning to be maximally predictive of rearrangements. Experimental results and computations enabled us to measure the spatial correlations and strain response of softness, as well as two measures of plasticity: the size of rearrangements and the yield strain. All four quantities maintained remarkable commonality in their values for disordered packings of objects ranging from atoms to grains, spanning seven orders of magnitude in diameter and 13 orders of magnitude in elastic modulus. These commonalities link the spatial correlations and strain response of softness to rearrangement size and yield strain, respectively.

We study experimentally the dynamics of granular media in a discharging hopper. In such flows, there often appears to be a critical outlet size D_{c} such that the flow never clogs for DD_{c}. We report on the time-averaged velocity distributions, as well as temporal intermittency in the ensemble-averaged velocity of grains in a viewing window, for both DD_{c} and DD_{c}, near and far from the outlet. We characterize the velocity distributions by the standard deviation and the skewness of the distribution of vertical velocities. We propose a measure for intermittency based on the two-sample Kolmogorov-Smirnov D_{KS} statistic for the velocity distributions as a function of time. We find that there is no discontinuity or kink in these various measures as a function of hole size. This result supports the proposition that there is no well-defined D_{c} and that clogging is always possible. Furthermore, the intermittency time scale of the flow is set by the speed of the grains at the hopper exit. This latter finding is consistent with a model of clogging as the independent sampling for stable configurations at the exit with a rate set by the exiting grain speed [C. C. Thomas and D. J. Durian, Phys. Rev. Lett. 114, 178001 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.178001].

Characterizing structural inhomogeneity is an essential step in understanding the mechanical response of amorphous materials. We introduce a threshold-free measure based on the field of vectors pointing from the center of each particle to the centroid of the Voronoi cell in which the particle resides. These vectors tend to point in toward regions of high free volume and away from regions of low free volume, reminiscent of sinks and sources in a vector field. We compute the local divergence of these vectors, for which positive values correspond to overpacked regions and negative values identify underpacked regions within the material. Distributions of this divergence are nearly Gaussian with zero mean, allowing for structural characterization using only the moments of the distribution. We explore how the standard deviation and skewness vary with packing fraction for simulations of bidisperse systems and find a kink in these moments that coincides with the jamming transition.

Submerged granular hoppers exhibit an unexpected surge in discharge rate as they empty [Wilson et al. 2015]. With a more sensitive apparatus, we find that this surge depends on hopper diameter and also happens in air --- though the effect is smaller and previously unnoticed. We also find that the surge may be turned off by fixing the rate of fluid flow through the granular packing. With no flow control, dye injected on top of the packing gets drawn into the grains, at a rate that increases as the hopper empties. Thus we conclude that the surge is caused by a self-generated pumping of fluid through the packing. We successfully model this effect via a driving pressure set by the dilation of grains as they exit. This highlights a surprising and unrecognized role that interstitial fluid plays in setting the discharge rate, and likely also in controlling clog formation, for granular hoppers whether in air or under water., Contains main article (5 pages, 4 figures) and supplemental material (8 pages, 7 figures). The supplemental video can be found at https://youtu.be/gRdiKEZbhA4

We report on length statistics measured for bubbles in the equilibrium bamboo state, achieved by the coarsening of aqueous foam in long cylindrical tubes, such that the soap films are all flat and perpendicular to the axis of the tube. The average bubble length is found to be 0.88 times the tube diameter, independent of variation of the liquid filling fraction by a factor of nearly three. The actual distribution is well-approximated by a shifted Rayleigh form, with a minimum bubble size of 0.28 tube diameters. And, perhaps surprisingly, no correlations are found in the lengths of neighboring bubbles. The observed length distribution agrees with that of Fortes et al. for short bubbles, but not for long bubbles.

Experiments on the low-speed impact of solid objects into granular media have been used both to mimic geophysical events1,2,3,4,5 and to probe the unusual nature of the granular state of matter6,7,8,9,10. Observations have been interpreted in terms of conflicting stopping forces: product of powers of projectile depth and speed6; linear in speed7; constant, proportional to the initial impact speed8; and proportional to depth9,10. This is reminiscent of high-speed ballistics impact in the nineteenth and twentieth centuries, when a plethora of empirical rules were proposed11,12. To make progress, we developed a means to measure projectile dynamics with 100 nm and 20 μs precision. For a 1-inch-diameter steel sphere dropped from a wide range of heights into non-cohesive glass beads, we reproduce previous observations6,7,8,9,10 either as reasonable approximations or as limiting behaviours. Furthermore, we demonstrate that the interaction between the projectile and the medium can be decomposed into the sum of velocity-dependent inertial drag plus depth-dependent friction. Thus, we achieve a unified description of low-speed impact phenomena and show that the complex response of granular materials to impact, although fundamentally different from that of liquids and solids, can be simply understood.

Measurement of growing dynamical length scales and prediction of the jamming transition in a granular material

Understanding the dynamics of fluid-driven sediment transport remains challenging, as it occurs at the interface between a granular material and a fluid flow. Boyer, Guazzelli, and Pouliquen [Phys. Rev. Lett. 107, 188301 (2011)]PRLTAO0031-900710.1103/PhysRevLett.107.188301 proposed a local rheology unifying dense dry-granular and viscous-suspension flows, but it has been validated only for neutrally buoyant particles in a confined and homogeneous system. Here we generalize the Boyer, Guazzelli, and Pouliquen model to account for the weight of a particle by addition of a pressure P_{0} and test the ability of this model to describe sediment transport in an idealized laboratory river. We subject a bed of settling plastic particles to a laminar-shear flow from above, and use refractive-index-matching to track particles' motion and determine local rheology-from the fluid-granular interface to deep in the granular bed. Data from all experiments collapse onto a single curve of friction μ as a function of the viscous number I_{v} over the range 3×10^{-5}≤I_{v}≤2, validating the local rheology model. For I_{v}

We use machine learning methods on local structure to identify flow defects - or regions susceptible to rearrangement - in jammed and glassy systems. We apply this method successfully to two disparate systems: a two dimensional experimental realization of a granular pillar under compression, and a Lennard-Jones glass in both two and three dimensions above and below its glass transition temperature. We also identify characteristics of flow defects that differentiate them from the rest of the sample. Our results show it is possible to discern subtle structural features responsible for heterogeneous dynamics observed across a broad range of disordered materials., 4 pages, 4 figures

We report a combined experimental and simulation study of deformation-induced diffusion in compacted quasi-two-dimensional amorphous granular pillars, in which thermal fluctuations play a negligible role. The pillars, consisting of bidisperse cylindrical acetal plastic particles standing upright on a substrate, are deformed uniaxially and quasistatically by a rigid bar moving at a constant speed. The plastic flow and particle rearrangements in the pillars are characterized by computing the best-fit affine transformation strain and nonaffine displacement associated with each particle between two stages of deformation. The nonaffine displacement exhibits exponential crossover from ballistic to diffusive behavior with respect to the cumulative deviatoric strain, indicating that in athermal granular packings, the cumulative deviatoric strain plays the role of time in thermal systems and drives effective particle diffusion. We further study the size-dependent deformation of the granular pillars by simulation, and find that different-sized pillars follow self-similar shape evolution during deformation. In addition, the yield stress of the pillars increases linearly with pillar size. Formation of transient shear lines in the pillars during deformation becomes more evident as pillar size increases. The width of these elementary shear bands is about twice the diameter of a particle, and does not vary with pillar size.

We present new data for the electrical conductivity of foams in which the liquid fraction ranges from two to eighty percent. We compare with a comprehensive collection of prior data, and we model all results with simple empirical formul\ae. We achieve a unified description that applies equally to dry foams and emulsions, where the droplets are highly compressed, as well as to dilute suspensions of spherical particles, where the particle separation is large. In the former limit, Lemlich's result is recovered; in the latter limit, Maxwell's result is recovered.

We report on the absorption of diffuse photons in aqueous foams by a dye added to the continuous liquid phase. For very wet and for dry foams, the absorption of the diffuse photons equals the absorption length of the liquid divided by the liquid volume fraction. This indicates that the diffuse photons propagate by a random walk, sampling each phase in proportion to its volume. Foams of intermediate wetness, by contrast, absorb photons more strongly than expected. A 2D computer simulation, modeling photons scattering in a foam crystal, also shows enhanced absorption. This encourages us to consider novel transport effects, such as the total internal reflection of photons inside the Plateau borders.

Fluid-sheared granular transport sculpts landscapes and undermines infrastructure, yet predicting the onset of sediment transport remains notoriously unreliable. For almost a century, this onset has been treated as a discontinuous transition at which hydrodynamic forces overcome gravity-loaded grain–grain friction. Using a custom laminar-shear flume to image slow granular dynamics deep into the bed, here we find that the onset is instead a continuous transition from creeping to granular flow. This transition occurs inside the dense granular bed at a critical viscous number, similar to granular flows and colloidal suspensions and inconsistent with hydrodynamic frameworks. We propose a new phase diagram for sediment transport, where ‘bed load’ is a dense granular flow bounded by creep below and suspension above. Creep is characteristic of disordered solids and reminiscent of soil diffusion on hillslopes. Results provide new predictions for the onset and dynamics of sediment transport that challenge existing models., The onset of sediment transport is commonly believed to be an abrupt process, which influences our understanding of river erosion and landscape dynamics. Houssais et al. show instead a continuous transition from a creeping to a granular-flow regime by monitoring internal motion in a sediment bed.

We measure the fraction $F$ of flowing grain configurations that precede a clog, based on the average mass discharged between clogging events for various aperture geometries. By tilting the hopper, we demonstrate that $F$ is a function of the hole area projected in the direction of the exiting grain velocity. By varying the length of slits, we demonstrate that grains clog in the same manner as if they were flowing out of a set of smaller independent circular openings. The collapsed data for $F$ can be fit to a decay that is exponential in hole width raised to the power of the system dimensionality. This is consistent with a simple model in which individual grains near the hole have a large but constant probability to precede a clog. Such a picture implies that there is no sharp clogging transition, and that all hoppers have a nonzero probability to clog. See Supplemental Material for models of clogging as a discrete Poisson process, and for resulting alternative measures of l based on the standard deviation of the discharge mass distribution., Comment: Supplemental Material available at PRL and also with source files

We investigate the formation of fingered flow in dry granular media under simulated rainfall using a quasi-2D experimental set-up composed of a random close packing of mono-disperse glass beads. Using controlled experiments, we analyze the finger instabilities that develop from the wetting front as a function of fundamental granular (particle size) and fluid properties (rainfall, viscosity).These finger instabilities act as precursors for water channels, which serve as outlets for water drainage. We look into the characteristics of the homogeneous wetting front and channel size as well as estimate relevant time scales involved in the instability formation and the velocity of the channel finger tip. We compare our experimental results with that of the well-known prediction developed by Parlange and Hill [1976]. This model is based on linear stability analysis of the growth of perturbations arising at the interface between two immiscible fluids. Results show that in terms of morphology, experiments agree with the proposed model. However, in terms of kinetics we nevertheless account for another term that describes the homogenization of the wetting front. This result shows that the manner we introduce the fluid to a porous medium can also influence the formation of finger instabilities., 13 pages, 7 figures

We visualize the formation of fingered flow in dry model sandy soils under different raining conditions using a quasi-2d experimental set-up, and systematically determine the impact of soil grain diameter and surface wetting property on water channelization phenomenon. The model sandy soils we use are random closely-packed glass beads with varied diameters and surface treatments. For hydrophilic sandy soils, our experiments show that rain water infiltrates into a shallow top layer of soil and creates a horizontal water wetting front that grows downward homogeneously until instabilities occur to form fingered flows. For hydrophobic sandy soils, in contrast, we observe that rain water ponds on the top of soil surface until the hydraulic pressure is strong enough to overcome the capillary repellency of soil and create narrow water channels that penetrate the soil packing. Varying the raindrop impinging speed has little influence on water channel formation. However, varying the rain rate causes significant changes in water infiltration depth, water channel width, and water channel separation. At a fixed raining condition, we combine the effects of grain diameter and surface hydrophobicity into a single parameter and determine its influence on water infiltration depth, water channel width, and water channel separation. We also demonstrate the efficiency of several soil water improvement methods that relate to rain water channelization phenomenon, including pre-wetting sandy soils at different level before rainfall, modifying soil surface flatness, and applying superabsorbent hydrogel particles as soil modifiers.

We present experimental measurements of penetration depths for the impact of spheres into wetted granular media. We observe that the penetration depth in the liquid saturated case scales with projectile density, size, and drop height in a fashion consistent with the scaling observed in the dry case, but that penetration depths into saturated packings tend to be smaller. This result suggests that, for the range of impact energies observed, the stopping force is set by static contact forces between grains within the bed, and that the presence of liquid serves, primarily, to enhance these contact forces. The enhancement to the stopping force has a complicated dependence on liquid fraction, accompanied by a change in the drop-height dependence, that must be the consequence of accompanying changes in the conformation of the liquid phase in the interstices., Comment: 7 pages, 4 figures, to be published in Phys. Rev. E