Grain shape effects in bed load sediment transport.

Bed load sediment transport, in which wind or water flowing over a bed of sediment causes grains to roll or hop along the bed, is a critically important mechanism in contexts ranging from river restoration1 to planetary exploration2. Despite its widespread occurrence, predictions of bed load sediment flux are notoriously imprecise3,4. Many studies have focused on grain size variability5 as a source of uncertainty, but few have investigated the role of grain shape, even though shape has long been suspected to influence transport rates6. Here we show that grain shape can modify bed load transport rates by an amount comparable to the scatter in many sediment transport datasets4,7,8. We develop a theory that accounts for grain shape effects on fluid drag and granular friction and predicts that the onset and efficiency of transport depend on the coefficients of drag and bulk friction of the transported grains. Laboratory experiments confirm these predictions and reveal that the effect of grain shape on sediment transport can be difficult to intuit from the appearance of grains. We propose a shape-corrected sediment transport law that collapses our experimental measurements. Our results enable greater accuracy in predictions of sediment transport and help reconcile theories developed for spherical particles with the behaviour of natural sediment grains.

Mechanistic framework for reduced-order models in soft materials: Application to three-dimensional granular intrusion.

Soft materials often display complex behaviors that transition through apparent solid- and fluid-like regimes. While a growing number of microscale simulation methods exist for these materials, reduced-order models that encapsulate the macroscale physics are often desired to predict how external bodies interact with soft media. Such an approach could provide direct insights in diverse situations from impact and penetration problems to locomotion over natural terrains. This work proposes a systematic program to develop three-dimensional (3D) reduced-order models for soft materials from a fundamental basis using continuum symmetries and rheological principles. In particular, we derive a reduced-order, 3D resistive force theory (3D-RFT), which is capable of accurately and quickly predicting the resistive stress distribution on arbitrary-shaped bodies intruding through granular media. Aided by a continuum description of the granular medium, a comprehensive set of spatial symmetry constraints, and a limited amount of reference data, we develop a self-consistent and accurate 3D-RFT. We verify the model capabilities in a wide range of cases and show that it can be quickly recalibrated to different media and intruder surface types. The premises leading to 3D-RFT anticipate application to other soft materials with strongly hyperlocalized intrusion behavior.

A general constitutive model for dense, fine-particle suspensions validated in many geometries.

Fine-particle suspensions (such as cornstarch mixed with water) exhibit dramatic changes in viscosity when sheared, producing fascinating behaviors that captivate children and rheologists alike. Examination of these mixtures in simple flow geometries suggests intergranular repulsion and its influence on the frictional nature of granular contacts is central to this effect-for mixtures at rest or shearing slowly, repulsion prevents frictional contacts from forming between particles, whereas when sheared more forcefully, granular stresses overcome the repulsion allowing particles to interact frictionally and form microscopic structures that resist flow. Previous constitutive studies of these mixtures have focused on particular cases, typically limited to 2D, steady, simple shearing flows. In this work, we introduce a predictive and general, 3D continuum model for this material, using mixture theory to couple the fluid and particle phases. Playing a central role in the model, we introduce a microstructural state variable, whose evolution is deduced from small-scale physical arguments and checked with existing data. Our space- and time-dependent model is implemented numerically in a variety of unsteady, nonuniform flow configurations where it is shown to accurately capture a variety of key behaviors: 1) the continuous shear-thickening (CST) and discontinuous shear-thickening (DST) behavior observed in steady flows, 2) the time-dependent propagation of "shear jamming fronts," 3) the time-dependent propagation of "impact-activated jamming fronts," and 4) the non-Newtonian, "running on oobleck" effect, wherein fast locomotors stay afloat while slow ones sink.

Interplay between hysteresis and nonlocality during onset and arrest of flow in granular materials.

The jamming transition in granular materials is well-known for exhibiting hysteresis, wherein the level of shear stress required to trigger flow is larger than that below which flow stops. Although such behavior is typically modeled as a simple non-monotonic flow rule, the rheology of granular materials is also nonlocal due to cooperativity at the grain scale, leading for instance to increased strengthening of the flow threshold as system size is reduced. We investigate how these two effects - hysteresis and nonlocality - couple with each other by incorporating non-monotonicity of the flow rule into the nonlocal granular fluidity (NGF) model, a nonlocal constitutive model for granular flows. By artificially tuning the strength of nonlocal diffusion, we demonstrate that both ingredients are key to explaining certain features of the hysteretic transition between flow and arrest. Finally, we assess the ability of the NGF model to quantitatively predict material behavior both around the transition and in the flowing regime, through stress-driven discrete element method (DEM) simulations of flow onset and arrest in various geometries. Along the way, we develop a new methodology to compare deterministic model predictions with the stochastic behavior exhibited by the DEM simulations around the jamming transition.

Efficacy of simple continuum models for diverse granular intrusions.

Granular intrusion is commonly observed in natural and human-made settings. Unlike typical solids and fluids, granular media can simultaneously display fluid-like and solid-like characteristics in a variety of intrusion scenarios. This multi-phase behavior increases the difficulty of accurately modeling these and other yielding (or flowable) materials. Micro-scale modeling methods, such as DEM (Discrete Element Method), capture this behavior by modeling the media at the grain scale, but there is often interest in the macro-scale characterizations of such systems. We examine the efficacy of a macro-scale continuum approach in modeling and understanding the physics of various macroscopic phenomena in a variety of granular intrusion cases using two basic frictional yielding constitutive models. We compare predicted granular force response and material flow to experimental data in four quasi-2D intrusion cases: (1) depth-dependent force response in horizontal submerged-intruder motion; (2) separation-dependent drag variation in parallel-plate vertical-intrusion; (3) initial-density-dependent drag fluctuations in free surface plowing, and (4) flow zone development during vertical plate intrusions in under-compacted granular media. Our continuum modeling approach captures the flow process and drag forces while providing key meso- and macro-scopic insights. The modeling results are then compared to experimental data. Our study highlights how continuum modeling approaches provide an alternative for efficient modeling as well as a conceptual understanding of various granular intrusion phenomena.

Power-Law Scaling in Granular Rheology across Flow Geometries.

Based on discrete element method simulations, we propose a new form of the constitutive equation for granular flows independent of packing fraction. Rescaling the stress ratio µ by a power of dimensionless temperature Θ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number I. The basic power-law structure appears robust to varying particle properties (e.g., surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the nonlocal granular fluidity model.

Intrusion rheology in grains and other flowable materials.

The interaction of intruding objects with deformable materials arises in many contexts, including locomotion in fluids and loose media, impact and penetration problems, and geospace applications. Despite the complex constitutive behaviour of granular media, forces on arbitrarily shaped granular intruders are observed to obey surprisingly simple, yet empirical 'resistive force hypotheses'. The physics of this macroscale reduction, and how it might play out in other media, has however remained elusive. Here, we show that all resistive force hypotheses in grains arise from local frictional yielding, revealing a novel invariance within a class of plasticity models. This mechanical foundation, supported by numerical and experimental validations, leads to a general analytical criterion to determine which rheologies can obey resistive force hypotheses. We use it to explain why viscous fluids are observed to perform worse than grains, and to predict a new family of resistive-force-obeying materials: cohesive media such as pastes, gels and muds.

Microscopic Description of the Granular Fluidity Field in Nonlocal Flow Modeling.

A recent granular rheology based on an implicit "granular fluidity" field has been shown to quantitatively predict many nonlocal phenomena. However, the physical nature of the field has not been identified. Here, the granular fluidity is found to be a kinematic variable given by the velocity fluctuation and packing fraction. This is verified with many discrete element simulations, which show that the operational fluidity definition, solutions of the fluidity model, and the proposed microscopic formula all agree. Kinetic theoretical and Eyring-like explanations shed insight into the obtained form.

Coupled dynamics of flow, microstructure, and conductivity in sheared suspensions.

We propose a model for the evolution of the conductivity tensor for a flowing suspension of electrically conductive particles. We use discrete particle numerical simulations together with a continuum physical framework to construct an evolution law for the suspension microstructure during flow. This model is then coupled with a relationship between the microstructure and the electrical conductivity tensor. Certain parameters of the joint model are fit experimentally using rheo-electrical conductivity measurements of carbon black suspensions under flow over a range of shear rates. The model is applied to the case of steady shearing as well as time-varying conductivity of unsteady flow experiments. We find that the model prediction agrees closely with the measured experimental data in all cases.

A predictive, size-dependent continuum model for dense granular flows.

Dense granular materials display a complicated set of flow properties, which differentiate them from ordinary fluids. Despite their ubiquity, no model has been developed that captures or predicts the complexities of granular flow, posing an obstacle in industrial and geophysical applications. Here we propose a 3D constitutive model for well-developed, dense granular flows aimed at filling this need. The key ingredient of the theory is a grain-size-dependent nonlocal rheology--inspired by efforts for emulsions--in which flow at a point is affected by the local stress as well as the flow in neighboring material. The microscopic physical basis for this approach borrows from recent principles in soft glassy rheology. The size-dependence is captured using a single material parameter, and the resulting model is able to quantitatively describe dense granular flows in an array of different geometries. Of particular importance, it passes the stringent test of capturing all aspects of the highly nontrivial flows observed in split-bottom cells--a geometry that has resisted modeling efforts for nearly a decade. A key benefit of the model is its simple-to-implement and highly predictive final form, as needed for many real-world applications.

Modeling growth paths of interacting crack pairs in elastic media.

The problem of predicting the growth of a system of cracks, each crack influencing the growth of the others, arises in multiple fields. We develop an analytical framework toward this aim, which we apply to the 'En-Passant' family of crack growth problems, in which a pair of initially parallel, offset cracks propagate nontrivially toward each other under far-field opening stress. We utilize boundary integral and perturbation methods of linear elasticity, linear elastic fracture mechanics, and common crack opening criteria to calculate the first analytical model for curved En-Passant crack paths. The integral system is reduced under a hierarchy of approximations, producing three methods of increasing simplicity for computing crack paths. The last such method is a major highlight of this work, using an asymptotic matching argument to predict crack paths based on superposition of simple, single-crack fields. Within the corresponding limits of the three methods, all three are shown to agree with each other. We provide comparisons to exact results and existing experimental data to verify certain approximation steps.

Modeling tensorial conductivity of particle suspension networks.

Significant microstructural anisotropy is known to develop during shearing flow of attractive particle suspensions. These suspensions, and their capacity to form conductive networks, play a key role in flow-battery technology, among other applications. Herein, we present and test an analytical model for the tensorial conductivity of attractive particle suspensions. The model utilizes the mean fabric of the network to characterize the structure, and the relationship to the conductivity is inspired by a lattice argument. We test the accuracy of our model against a large number of computer-generated suspension networks, based on multiple in-house generation protocols, giving rise to particle networks that emulate the physical system. The model is shown to adequately capture the tensorial conductivity, both in terms of its invariants and its mean directionality.

Nonlocal modeling of granular flows down inclines.

Flows of granular media down a rough inclined plane demonstrate a number of nonlocal phenomena. We apply the recently proposed nonlocal granular fluidity model to this geometry and find that the model captures many of these effects. Utilizing the model's dynamical form, we obtain a formula for the critical stopping height of a layer of grains on an inclined surface. Using an existing parameter calibration for glass beads, the theoretical result compares quantitatively to existing experimental data for glass beads. This provides a stringent test of the model, whose previous validations focused on driven steady-flow problems. For layers thicker than the stopping height, the theoretical flow profiles display a thickness-dependent shape whose features are in agreement with previous discrete particle simulations. We also address the issue of the Froude number of the flows, which has been shown experimentally to collapse as a function of the ratio of layer thickness to stopping height. While the collapse is not obvious, two explanations emerge leading to a revisiting of the history of inertial rheology, which the nonlocal model references for its homogeneous flow response.

Continuum modeling of secondary rheology in dense granular materials.

Recent dense granular flow experiments have shown that shear deformation in one region of a granular medium fluidizes its entirety, including regions far from the sheared zone, effectively erasing the yield condition everywhere. This enables slow creep deformation to occur when an external force is applied to a probe in the nominally static regions of the material. The apparent change in rheology induced by far-away motion is termed the "secondary rheology," and a theoretical rationalization of this phenomenon is needed. Recently, a new nonlocal granular rheology was successfully used to predict steady granular flow fields, including grain-size-dependent shear-band widths in a wide variety of flow configurations. We show that the nonlocal fluidity model is also capable of capturing secondary rheology. Specifically, we explore creep of a circular intruder in a two-dimensional annular Couette cell and show that the model captures all salient features observed in experiments, including both the rate-independent nature of creep for sufficiently slow driving rates and the faster-than-linear increase in the creep speed with the force applied to the intruder.

Bioinspired, ingestible electroceutical capsules for hunger-regulating hormone modulation.

The gut-brain axis, which is mediated via enteric and central neurohormonal signaling, is known to regulate a broad set of physiological functions from feeding to emotional behavior. Various pharmaceuticals and surgical interventions, such as motility agents and bariatric surgery, are used to modulate this axis. Such approaches, however, are associated with off-target effects or post-procedure recovery time and expose patients to substantial risks. Electrical stimulation has also been used to attempt to modulate the gut-brain axis with greater spatial and temporal resolution. Electrical stimulation of the gastrointestinal (GI) tract, however, has generally required invasive intervention for electrode placement on serosal tissue. Stimulating mucosal tissue remains challenging because of the presence of gastric and intestinal fluid, which can influence the effectiveness of local luminal stimulation. Here, we report the development of a bioinspired ingestible fluid-wicking capsule for active stimulation and hormone modulation (FLASH) capable of rapidly wicking fluid and locally stimulating mucosal tissue, resulting in systemic modulation of an orexigenic GI hormone. Drawing inspiration from Moloch horridus, the "thorny devil" lizard with water-wicking skin, we developed a capsule surface capable of displacing fluid. We characterized the stimulation parameters for modulation of various GI hormones in a porcine model and applied these parameters to an ingestible capsule system. FLASH can be orally administered to modulate GI hormones and is safely excreted with no adverse effects in porcine models. We anticipate that this device could be used to treat metabolic, GI, and neuropsychiatric disorders noninvasively with minimal off-target effects.

Nonlocal constitutive relation for steady granular flow.

Extending recent modeling efforts for emulsions, we propose a nonlocal fluidity relation for flowing granular materials, capturing several known finite-size effects observed in steady flow. We express the local Bagnold-type granular flow law in terms of a fluidity ratio and then extend it with a particular Laplacian term that is scaled by the grain size. The resulting model is calibrated against a sequence of existing discrete element method data sets for two-dimensional annular shear, where it is shown that the model correctly describes the divergence from a local rheology due to the grain size as well as the rate-independence phenomenon commonly observed in slowly flowing zones. The same law is then applied in two additional inhomogeneous flow geometries, and the predicted velocity profiles are compared against corresponding discrete element method simulations utilizing the same grain composition as before, yielding favorable agreement in each case.

Surprising simplicity in the modeling of dynamic granular intrusion.

Granular intrusions, such as dynamic impact or wheel locomotion, are complex multiphase phenomena where the grains exhibit solid-like and fluid-like characteristics together with an ejected gas-like phase. Despite decades of modeling efforts, a unified description of the physics in such intrusions is as yet unknown. Here, we show that a continuum model based on the simple notions of frictional flow and tension-free separation describes complex granular intrusions near free surfaces. This model captures dynamics in a variety of experiments including wheel locomotion, plate intrusions, and running legged robots. The model reveals that one static and two dynamic effects primarily give rise to intrusion forces in such scenarios. We merge these effects into a further reduced-order technique (dynamic resistive force theory) for rapid modeling of granular locomotion of arbitrarily shaped intruders. The continuum-motivated strategy we propose for identifying physical mechanisms and corresponding reduced-order relations has potential use for a variety of other materials.

Mitigating memory effects during undulatory locomotion on hysteretic materials.

While terrestrial locomotors often contend with permanently deformable substrates like sand, soil, and mud, principles of motion on such materials are lacking. We study the desert-specialist shovel-nosed snake traversing a model sand and find body inertia is negligible despite rapid transit and speed dependent granular reaction forces. New surface resistive force theory (RFT) calculation reveals how wave shape in these snakes minimizes material memory effects and optimizes escape performance given physiological power limitations. RFT explains the morphology and waveform-dependent performance of a diversity of non-sand-specialist snakes but overestimates the capability of those snakes which suffer high lateral slipping of the body. Robophysical experiments recapitulate aspects of these failure-prone snakes and elucidate how re-encountering previously deformed material hinders performance. This study reveals how memory effects stymied the locomotion of a diversity of snakes in our previous studies (Marvi et al., 2014) and indicates avenues to improve all-terrain robots.