Shear Thickening and Jamming of Dense Suspensions: The "Roll" of Friction

Order-disorder transition during shear thickening in bidisperse dense suspensions

Shear thickening of multimodal suspensions has proven difficult to understand because the rheology depends largely on the microscopic details of stress-induced frictional contacts at different particle size distributions (PSDs). Our discrete particle simulations below a critical volume fraction ϕc over a broad range of shear rates and PSDs elucidate the basic mechanism of order-disorder transition. Around the theoretical optimal PSD (relative content of small particles ζ1= 0.26), particles order into a layered structure in the Newtonian regime. At the onset of shear thickening, this layered structure transforms to a disordered one, accompanied by an abrupt viscosity jump. Minor increase in large-large particle contacts after the order-disorder transition causes apparent increase in radial force along the compressional axis. Bidisperse suspensions with less regular but stable layered structure at ζ1= 0.50 show good fluidity in the shear thickening regime. This work shows that in inertial flows where particle collisions dominate, order-disorder transition could play an essential role in shear thickening for bidisperse suspensions.

Optimizing non-Newtonian fluids for impact protection of laminates

Non-Newtonian fluids can be used for the protection of flexible laminates. Understanding the coupling between the flow of the protecting fluid and the deformation of the protected solids is necessary in order to optimize this functionality. We present a scaling analysis of the problem based on a single coupling variable, the effective width of a squeeze flow between flat rigid plates, and predict that impact protection for laminates is optimized by using shear-thinning, and not shear-thickening, fluids. The prediction is verified experimentally by measuring the velocity and pressure in impact experiments. Our scaling analysis should be generically applicable for non-Newtonian fluid-solid interactions in diverse applications.

Viscosimetric squeeze flow of suspensions

The rheology of particle suspensions has been extensively explored in the case of a simple shear flow, but less in other flow configurations which are also important in practice. Here we investigate the behavior of a suspension in a squeeze flow, which we revisit using local pressure measurements to deduce the effective viscosity. The flow is generated by approaching a moving disk to a fixed wall at constant velocity in the low Reynolds number limit. We measure the evolution of the pressure field at the wall and deduce the effective viscosity from the radial pressure drop. After validation of our device using a Newtonian fluid, we measure the effective viscosity of a suspension for different squeezing speeds and volume fractions of particles. We find results in agreement with the Maron-Pierce law, an empirical expression for the viscosity of suspensions that was established for simple shear flows. We prove that this method to determine viscosity remains valid in the limit of large gap width. This makes it possible to study the rheology of suspensions within this limit and therefore suspensions composed of large particles, in contrast to Couette flow cells which require small gaps.

Network physics of attractive colloidal gels: Resilience, rigidity, and phase diagram

Colloidal gels exhibit solid-like behavior at vanishingly small fractions of solids, owing to ramified space-spanning networks that form due to particle-particle interactions. These networks give the gel its rigidity, and with stronger attractions the elasticity grows as well. The emergence of rigidity can be described through a mean field approach; nonetheless, fundamental understanding of how rigidity varies in gels of different attractions is lacking. Moreover, recovering an accurate gelation phase diagram based on the system's variables has been an extremely challenging task. Understanding the nature of colloidal clusters, and how rigidity emerges from their connections is key to controlling and designing gels with desirable properties. Here, we employ network analysis tools to interrogate and characterize the colloidal structures. We construct a particle-level network, having all the spatial coordinates of colloids with different attraction levels, and also identify polydisperse rigid fractal clusters using a Gaussian mixture model, to form a coarse-grained cluster network that distinctly shows main physical features of the colloidal gels. A simple mass-spring model then is used to recover quantitatively the elasticity of colloidal gels from these cluster networks. Interrogating the resilience of these gel networks shows that the elasticity of a gel (a dynamic property) is directly correlated to its cluster network's resilience (a static measure). Finally, we use the resilience investigations to devise [and experimentally validate] a fully resolved phase diagram for colloidal gelation, with a clear solid-liquid phase boundary using a single volume fraction of particles well beyond this phase boundary.

Stress-activated friction in sheared suspensions probed with piezoelectric nanoparticles

A hallmark of concentrated suspensions is non-Newtonian behavior, whereby the viscosity increases dramatically once a characteristic shear rate or stress is exceeded. Such strong shear thickening is thought to originate from a network of frictional particle-particle contact forces, which forms under sufficiently large stress, evolves dynamically, and adapts to changing loads. While there is much evidence from simulations for the emergence of this network during shear thickening, experimental confirmation has been difficult. Here, we use suspensions of piezoelectric nanoparticles and exploit the strong local stress focusing within the network to activate charge generation. This charging can then be detected in the measured ac conductance and serve as a signature of frictional contact formation. The direct link between stress-activated frictional particle interactions and piezoelectric suspension response is further demonstrated by tracking the emergence of structural memory in the contact network under oscillatory shear and by showing how stress-activated friction can drive mechano-transduction of chemical reactions with nonlinear reaction kinetics. Taken together, this makes the ac conductance of piezoelectric suspensions a sensitive in-situ reporter of the micromechanics associated with frictional interactions.

Shear thickening in presence of adhesive contact forces: The singularity of cornstarch

HYPOTHESIS

A number of dense particle suspensions experience a dramatic increase in viscosity with the shear stress, up to a solid-like response. This shear-thickening process is understood as a transition under flow of the nature of the contacts - from lubricated to frictional - between initially repellent particles. Most systems are now assumed to fit in with this scenario, which is questionable.

EXPERIMENT

Using an in-house pressure sensor array, we provide a spatio-temporal map of the normal stresses in the flows of two shear-thickening fluids: a stabilized calcium carbonate suspension, known to fit in with the standard scenario, and a cornstarch suspension, which spectacular thickening behavior remains poorly understood.

FINDINGS

We evidence in cornstarch a unique, stable heterogeneous structure, which moves in the velocity direction and does not appear in calcium carbonate. Its nature changes from a stress wave to a rolling solid jammed aggregate at high solid fraction and small gap width. The modeling of these heterogenities points to an adhesive force between cornstarch particles at high stress, also evidenced in microscopic measurements. Cornstarch being also attractive at low stress, it stands out of the classical shear-thickening frame, and might be part of a larger family of adhesive and attractive shear-thickening fluids.

Strain-Controlled Critical Slowing Down in the Rheology of Disordered Networks

Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.

Scaling relationships between viscosity and diffusivity in shear-thickening suspensions

Dense suspensions often exhibit a dramatic response to large external deformation. The recent body of work has related this behavior to transition from an unconstrained lubricated state to a constrained frictional state. Here, we use numerical simulations to study the flow behavior and shear-induced diffusion of frictional non-Brownian spheres in two dimensions under simple shear flow. We first show that both viscosity η and diffusivity D/ of the particles increase under characteristic shear stress, which is associated with lubrication to frictional transition. Subsequently, we propose a one-to-one relationship between viscosity and diffusivity using the length scale ξ associated with the size of collective motions (rigid clusters) of the particles. We demonstrate that η and D/ are controlled by ξ in two distinct flow regimes, i.e. in the frictionless and frictional states, where the one-to-one relationship is described as a crossover from D/ ∼ η (frictionless) to η1/3 (frictional). We also confirm that the proposed power laws are insensitive to the interparticle friction and system size.

Dense bidisperse suspensions under non-homogeneous shear

We study the rheological behaviour of bidisperse suspensions in three dimensions under a non-uniform shear flow, made by the superimposition of a linear shear and a sinusoidal disturbance. Our results show that (i) only a streamwise disturbance in the shear-plane alters the suspension dynamics by substantially reducing the relative viscosity, (ii) with the amplitude of the disturbance determining a threshold value for the effect to kick-in and its wavenumber controlling the amount of reduction and which of the two phases is affected. We show that, (iii) the rheological changes are caused by the effective separation of the two phases, with the large or small particles layering in separate regions. We provide a physical explanation of the phase separation process and of the conditions necessary to trigger it. We test the results in the whole flow curve, and we show that the mechanism remains substantially unaltered, with the only difference being the nature of the interactions between particles modified by the phase separation.

Effective viscosity and elasticity in dense suspensions under impact: Toward a modeling of walking on suspensions

The elastic response of dense suspensions under an impact is studied using coupled lattice Boltzmann method and discrete element method (LBM-DEM) and its reduced model. We succeed to extract the elastic force acting on the impactor in dense suspensions, which can exist even in the absence of percolating clusters of suspended particles. We then propose a reduced model to describe the motion of the impactor and demonstrate its relevancy through the comparison of the solution of the reduced model and that of LBM-DEM. Furthermore, we illustrate that the perturbation analysis of the reduced model captures the short-time behavior of the impactor motion quantitatively. We apply this reduced model to the impact of a foot-spring-body system on a dense suspension, which is the minimal model to realize walking on the suspension. Due to the spring force of the system and the stiffness of the suspension, the foot undergoes multiple bounces. We also study the parameter dependencies of the hopping motion and find that multiple bounces are suppressed as the spring stiffness increases.

3-D rotation tracking from 2-D images of spherical colloids with textured surfaces

Tracking the three-dimensional rotation of colloidal particles is essential to elucidate many open questions, e.g. concerning the contact interactions between particles under flow, or the way in which obstacles and neighboring particles affect self-propulsion in active suspensions. In order to achieve rotational tracking, optically anisotropic particles are required. We synthesise here rough spherical colloids that present randomly distributed fluorescent asperities and track their motion under different experimental conditions. Specifically, we propose a new algorithm based on a 3-D rotation registration, which enables us to track the 3-D rotation of our rough colloids at short time-scales, using time series of 2-D images acquired at high frame rates with a conventional wide-field microscope. The method is based on the image correlation between a reference image and rotated 3-D prospective images to identify the most likely angular displacements between frames. We first validate our approach against simulated data and then apply it to the cases of: particles flowing through a capillary, freely diffusing at solid-liquid and liquid-liquid interfaces, and self-propelling above a substrate. By demonstrating the applicability of our algorithm and sharing the code, we hope to encourage further investigations in the rotational dynamics of colloidal systems.

Leveraging the Polymer Glass Transition to Access Thermally Switchable Shear Jamming Suspensions

Suspensions of polymeric nano- and microparticles are fascinating stress-responsive material systems that, depending on their composition, can display a diverse range of flow properties under shear, such as drastic thinning, thickening, and even jamming (reversible solidification driven by shear). However, investigations to date have almost exclusively focused on nonresponsive particles, which do not allow in situ tuning of the flow properties. Polymeric materials possess rich phase transitions that can be directly tuned by their chemical structures, which has enabled researchers to engineer versatile adaptive materials that can respond to targeted external stimuli. Reported herein are suspensions of (readily prepared) micrometer-sized polymeric particles with accessible glass transition temperatures (T _g) designed to thermally control their non-Newtonian rheology. The underlying mechanical stiffness and interparticle friction between particles change dramatically near T _g. Capitalizing on these properties, it is shown that, in contrast to conventional systems, a dramatic and nonmonotonic change in shear thickening occurs as the suspensions transition through the particles' T _g. This straightforward strategy enables the in situ turning on (or off) of the system's ability to shear jam by varying the temperature relative to T _g and lays the groundwork for other types of stimuli-responsive jamming systems through polymer chemistry.

How surface roughness affects the interparticle interactions at a liquid interface

HYPOTHESIS

Colloidal particles can be trapped at a liquid interface, which reduces the energetically costly interfacial area. Once at an interface, colloids undergo various self-assemblies and structural transitions due to shape-dependent interparticle interactions. Particles with rough surfaces receive increasing attention and have been applied in material design, such as Pickering emulsions and shear-thickening materials. However, the roughness effects on the interactions at a liquid interface remain less understood.

EXPERIMENTS

Experimentally, particles with four surface roughnesses were designed and compared via isotherm measurements upon a uniaxial compression. At each stage of the compression, micrographic observations were conducted via the Blodgett method. Numerically, the compression of monolayer was simulated by using Langevin dynamics. Rough colloids were modelled as particles with capillary attraction and tangential constraints.

FINDINGS

Sufficiently rough systems exhibit a non-trivial intermediate state between a gas-like state and a close-packed jamming state. This state is understood as a gel state due to roughness-induced capillary attraction. Roughness-induced friction lowers the jamming point. Furthermore, the tangential contact force owing to surface asperities can cause a gradual off-plane collapse of the compressed monolayer.

Shear jamming and fragility in fractal suspensions under confinement

Under applied stress, the viscosity of many dense particulate suspensions increases drastically, a response known as discontinuous shear-thickening (DST). In some cases, the applied stress can even transform the suspension into a solid-like shear jammed state. Although shear jamming (SJ) has been probed for dense suspensions with particles having well-defined shapes, such a phenomenon for fractal objects has not been explored. Here, using rheology and in situ optical imaging, we study the flow behaviour of ultra-dilute fractal suspensions of multi-walled carbon nanotubes (MWCNT) under confinement. We show a direct transition from flowing to SJ state without a precursory DST in fractal suspensions at an onset volume fraction, ϕ ∼ 0.5%, significantly lower than that of conventional dense suspensions (ϕ ∼ 55%). The ultra-low concentration enables us to demonstrate the fragility and associated contact dynamics of the SJ state, which remain experimentally unexplored in suspensions. Furthermore, using a generalized Wyart-Cates model, we propose a generic phase diagram for fractal suspensions that captures the possibility of SJ without prior DST over a wide range of shear stress and volume fractions.

Chiral Liquid Crystalline Properties of Cellulose Nanocrystals: Fundamentals and Applications

By using an independent self-assembly process that is occasionally controlled by evaporation, cellulose nanocrystals (CNCs) may create films (pure or in conjunction with other materials) that have iridescent structural colors. The self-forming chiral nematic structures and environmental safety of a new class of photonic liquid crystals (LCs), referred to as CNCs and CNC-embedded materials, make them simple to make and treat. The structure of the matrix interacts with light to give structural coloring, as opposed to other dye pigments, which interact with light by adsorption and reflection. Understanding how CNC self-assembly constructs structures is vital in several fields, including physics, science, and engineering. To constructure this review, the colloidal characteristics of CNC particles and their behavior during the formation of liquid crystals and gelling were studied. Then, some of the recognized applications for these naturally occurring nanoparticles were summarized. Different factors were considered, including the CNC aspect ratio, surface chemistry, concentration, the amount of time needed to produce an anisotropic phase, and the addition of additional substances to the suspension medium. The effects of alignment and the drying process conditions on structural changes are also covered. The focus of this study however is on the optical properties of the films as well as the impact of the aforementioned factors on the final transparency, iridescent colors, and versus the overall response of these bioinspired photonic materials. Control of the examined factors was found to be necessary to produce reliable materials for optoelectronics, intelligent inks and papers, transparent flexible support for electronics, and decorative coatings and films.

Structure and Dynamics of Force Clusters and Networks in Shear Thickening Suspensions

Dense suspensions can exhibit shear thickening in response to large deformation. A consensus has emerged over the past few years on the formation of force networks, that span the entire system size, that lead to increased resistance to motion. Nonetheless, the characteristics of these networks are to a large extent poorly understood. Here, force networks formed in continuous and discontinuous shear thickening dense suspensions (CST and DST, respectively) are studied. We first show the evolution of the network formation and its topological heterogeneities as the applied stress increases. Subsequently, we identify force communities and coarse grain the suspension into a cluster network, and show that cluster-level dynamics are responsible for stark differences between the CST and DST behavior. Our results suggest that the force clusters formed in the DST regime are considerably more constrained in their motion, while CST clusters are loosely connected to their surrounding clusters.

In Situ Observation of Shear-Induced Jamming Front Propagation during Low-Velocity Impact in Polypropylene Glycol/Fumed Silica Shear Thickening Fluids

Shear jamming, a relatively new type of phase transition from discontinuous shear thickening into a solid-like state driven by shear in dense suspensions, has been shown to originate from frictional interactions between particles. However, not all dense suspensions shear jam. Dense fumed silica colloidal systems have wide applications in the industry of smart materials from body armor to dynamic dampers due to extremely low bulk density and high colloid stability. In this paper, we provide new evidence of shear jamming in polypropylene glycol/fumed silica suspensions using optical in situ speed recording during low-velocity impact and explain how it contributes to impact absorption. Flow rheology confirmed the presence of discontinuous shear thickening at all studied concentrations. Calculations of the flow during impact reveal that front propagation speed is 3–5 times higher than the speed of the impactor rod, which rules out jamming by densification, showing that the cause of the drastic impact absorption is the shear jamming. The main impact absorption begins when the jamming front reaches the boundary, creating a solid-like plug under the rod that confronts its movement. These results provide important insights into the impact absorption mechanism in fumed silica suspensions with a focus on shear jamming.

Origin of Two Distinct Stress Relaxation Regimes in Shear Jammed Dense Suspensions

Many dense particulate suspensions show a stress induced transformation from a liquidlike state to a solidlike shear jammed (SJ) state. However, the underlying particle-scale dynamics leading to such striking, reversible transition of the bulk remains unknown. Here, we study transient stress relaxation behaviour of SJ states formed by a well-characterized dense suspension under a step strain perturbation. We observe a strongly nonexponential relaxation that develops a sharp discontinuous stress drop at short time for high enough peak-stress values. High resolution boundary imaging and normal stress measurements confirm that such stress discontinuity originates from the localized plastic events, whereas system spanning dilation controls the slower relaxation process. We also find an intriguing correlation between the nature of transient relaxation and the steady-state shear jamming phase diagram obtained from the Wyart-Cates model.

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology

Despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow, which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

Stress propagation in locally loaded packings of disks and pentagons

The mechanical strength and flow of granular materials can depend strongly on the shapes of individual grains. We report quantitative results obtained from photoelasticimetry experiments on locally loaded, quasi-two-dimensional granular packings of either disks or pentagons exhibiting stick-slip dynamics. Packings of pentagons resist the intruder at significantly lower packing fractions than packings of disks, transmitting stresses from the intruder to the boundaries over a smaller spatial extent. Moreover, packings of pentagons feature significantly fewer back-bending force chains than packings of disks. Data obtained on the forward spatial extent of stresses and back-bending force chains collapse when the packing fraction is rescaled according to the packing fraction of steady state open channel formation, though data on intruder forces and dynamics do not collapse. We comment on the influence of system size on these findings and highlight connections with the dynamics of the disks and pentagons during slip events.

Friction and Adhesion Govern Yielding of Disordered Nanoparticle Packings: A Multiscale Adhesive Discrete Element Method Study

Recent studies have demonstrated that amorphous materials, from granular packings to atomic glasses, share multiple striking similarities, including a universal onset strain level for yield. This is despite vast differences in length scales and in the constituent particles' interactions. However, the nature of localized particle rearrangements is not well understood, and how local interactions affect overall performance remains unknown. Here, we introduce a multiscale adhesive discrete element method to simulate recent novel experiments of disordered nanoparticle packings indented and imaged with single nanoparticle resolution. The simulations exhibit multiple behaviors matching the experiments. By directly monitoring spatial rearrangements and interparticle bonding/debonding under the packing's surface, we uncover the mechanisms of the yielding and hardening phenomena observed in experiments. Interparticle friction and adhesion synergistically toughen the packings and retard plastic deformation. Moreover, plasticity can result from bond switching without particle rearrangements. These results furnish insights for understanding yielding in amorphous materials generally.

Jamming Distance Dictates Colloidal Shear Thickening

We report experimental and computational observations of dynamic contact networks for colloidal suspensions undergoing shear thickening. The dense suspensions are comprised of sterically stabilized poly(methyl methacrylate) colloids that are spherically symmetric and have varied surface roughness. Confocal rheometry and dissipative particle dynamics simulations show that the shear thickening strength β scales exponentially with the scaled deficit contact number and the scaled jamming distance. Rough colloids, which experience additional rotational constraints, require an average of 1.5-2 fewer particle contacts as compared to smooth colloids, in order to generate the same β. This is because the surface roughness enhances geometric friction in such a way that the rough colloids do not experience a large change in the free volume near the jamming point. The available free volume for colloids of different roughness is related to the deficiency from the maximum number of nearest neighbors at jamming under shear. Our results further suggest that the force per contact is different for particles with different morphologies.

Roughness induced rotational slowdown near the colloidal glass transition

HYPOTHESIS

In concentrated suspensions, the dynamics of colloids are strongly influenced by the shape and topographical surface characteristics of the particles. As the particles get into close proximity, surface roughness alters the translational and rotational Brownian motions in different ways. Eventually, the rotations will get frustrated due to geometric hindrance from interacting asperities.

EXPERIMENTS

We use model raspberry-like colloids to study the effect of roughness on the translational and rotational dynamics. Using Confocal Scanning Laser Microscopy and particle tracking, we simultaneously resolve the two types of Brownian motion and obtain the corresponding Mean Squared Displacements for varying concentrations up to the maximum packing fraction.

FINDINGS

Roughness not only lowers the concentration of the translational colloidal glass transition, but also generates a broad concentration range in which the rotational Brownian motion changes signature from high-amplitude diffusive to low-amplitude rattling. This hitherto not reported second glass transition for rough spherical colloids emerges when the particle intersurface distance becomes comparable to the roughness length scale. Our work provides a unifying understanding of the surface characteristics' effect on the rotational dynamics during glass formation and provides a microscopic foundation for many roughness-related macroscale phenomena in nature and technology.

Shear Jamming and Fragility of Suspensions in a Continuum Model with Elastic Constraints

Under an applied traction, highly concentrated suspensions of solid particles in fluids can turn from a state in which they flow to a state in which they counteract the traction as an elastic solid: a shear-jammed state. Remarkably, the suspension can turn back to the flowing state simply by inverting the traction. A tensorial model is presented and tested in paradigmatic cases. We show that, to reproduce the phenomenology of shear jamming in generic geometries, it is necessary to link this effect to the elastic response supported by the suspension microstructure rather than to a divergence of the viscosity.

Rotational diffusion and rotational correlations in frictional amorphous disk packings under shear

We show here that rotations of round particles in amorphous disk packing reveal various nontrivial microscopic features when the packing is close to rigidification. We analyze experimental measurements on disk packing subjected to simple shear deformation with various inter-particle friction coefficients and across a range of volume fractions where the system is known to stiffen. The analysis of measurements indicates that shear induces diffusive microrotation, that can be both enhanced and suppressed depending upon the volume fraction as well as the inter-particle friction. Rotations also display persistent anticorrelated motion. Spatial correlations in microrotation are observed to be directly correlated with system pressure. These observations point towards the broader mechanical relevance of collective dynamics in the rotational degree of freedom of particles.

Unifying disparate rate-dependent rheological regimes in non-Brownian suspensions

A typical dense non-Brownian particulate suspension exhibits shear thinning (decreasing viscosity) at low shear rate or stress followed by a Newtonian plateau (constant viscosity) at intermediate shear rate or stress values which transitions to shear thickening (increasing viscosity) beyond a critical shear rate or stress value and finally undergoes a second shear thinning transition at extremely high shear rate or stress values. In this study, we unify and quantitatively reproduce all the disparate rate-dependent regimes and the corresponding transitions for a dense non-Brownian suspension with increasing shear rate or stress. We employ discrete particle dynamics simulations based on the proposed mechanism to elucidate its accuracy. We find that a competition between interparticle interactions of hydrodynamic and nonhydrodynamic origins and the switching in the dominant stress scale with increasing the shear rate or stress lead to each of the above transitions. Inclusion of traditional hydrodynamic interactions, attractive or repulsive Derjaguin-Landau-Verwey-Overbeek (DLVO) interactions the interparticle contact interactions, and a constant friction (or other constraint mechanism) reproduces the initial thinning as well as the shear thickening transition. However, to quantitatively capture the intermediate Newtonian plateau and the second shear thinning, an additional nonhydrodynamic interaction of non-DLVO origin and a decreasing coefficient of friction, respectively, are essential, thus providing an explanation for the presence of the intermediate Newtonian plateau along with reproducing the second shear thinning in a single framework. Expressions utilized for various interactions and friction are determined from experimental measurements and hence result in excellent quantitative agreement between the simulations and previous experiments.

The role of friction in statistics and scaling laws of avalanches

We investigate statistics and scaling laws of avalanches in two-dimensional frictional particles by numerical simulations. We find that the critical exponent for avalanche size distributions is governed by microscopic friction between the particles in contact, where the exponent is larger and closer to mean-field predictions if the friction coefficient is finite. We reveal that microscopic "slips" between frictional particles induce numerous small avalanches which increase the slope, as well as the power-law exponent, of avalanche size distributions. We also analyze statistics and scaling laws of the avalanche duration and maximum stress drop rates, and examine power spectra of stress drop rates. Our numerical results suggest that the microscopic friction is a key ingredient of mean-field descriptions and plays a crucial role in avalanches observed in real materials.

Particle dynamics in two-dimensional point-loaded granular media composed of circular or pentagonal grains

Granular packings exhibit significant changes in rheological and structural properties when the rotational symmetry of spherical or circular particles is broken. Here, we report on experiments exploring the differences in dynamics of a grain-scale intruder driven through a packing of either disks or pentagons, where the presence of edges and vertices on grains introduces the possibility of rotational constraints at edge-edge contacts. We observe that the intruder’s stick-slip dynamics are comparable between the disk packing near the frictional jamming fraction and the pentagonal packing at significantly lower packing fractions. We connect this stark contrast in packing fraction with the average speed and rotation fields of grains during slip events, finding that rotation of pentagons is limited and the flow of pentagonal grains is largely confined in front of the intruder, whereas disks rotate more on average and circulate around the intruder to fill the open channel behind it. Our results indicate that grain-scale rotation constraints significantly modify collective motion of grains on mesoscopic scales and correspondingly enhance resistance to penetration of a local intruder.

The role of solvent molecular weight in shear thickening and shear jamming

The application of stress can drive a dense suspension into a regime of highly non-Newtonian response, characterized by discontinuous shear thickening (DST) and potentially shear jamming (SJ), due to the formation of a frictionally stabilized contact network. Investigating how the molecular weight of the suspending solvent affects the frictional particle-particle interactions, we report on experiments with suspensions of fumed silica particles in polyethylene glycol (PEG). Focusing on the monomer-to-oligomer limit, with n = 1 to 8 ethylene oxide repeat units, we find that increasing n enhances shear thickening under steady-state shear and even elicits rapidly propagating shear jamming fronts, as assessed by high-speed ultrasound imaging of impact experiments. We associate this behavior with a weakening of the solvation layers surrounding the particles as n is increased, which thereby facilitates the formation of frictional contacts. We argue that for n larger than the monomer-to-oligomer limit the trend reverses and frictional interactions are diminished, as observed in prior experiments. This reversal occurs because the polymeric solvent transitions from being enthalpically bound to entropically bound to the particle surfaces, which strengthens solvation layers.

Exploring the roles of roughness, friction and adhesion in discontinuous shear thickening by means of thermo-responsive particles

Dense suspensions of colloidal or granular particles can display pronounced non-Newtonian behaviour, such as discontinuous shear thickening and shear jamming. The essential contribution of particle surface roughness and adhesive forces confirms that stress-activated frictional contacts can play a key role in these phenomena. Here, by employing a system of microparticles coated by responsive polymers, we report experimental evidence that the relative contributions of friction, adhesion, and surface roughness can be tuned in situ as a function of temperature. Modifying temperature during shear therefore allows contact conditions to be regulated, and discontinuous shear thickening to be switched on and off on demand. The macroscopic rheological response follows the dictates of independent single-particle characterization of adhesive and tribological properties, obtained by colloidal-probe atomic force microscopy. Our findings identify additional routes for the design of smart non-Newtonian fluids and open a way to more directly connect experiments to computational models of sheared suspensions.

Autonomously Probing Viscoelasticity in Disordered Suspensions

Recent experiments show a strong rotational diffusion enhancement for self-propelled microrheological probes in colloidal glasses. Here, we provide microscopic understanding using simulations with a frictional probe-medium coupling that converts active translation into rotation. Diffusive enhancement emerges from the medium's disordered structure and peaks at a second-order transition in the number of contacts. Our results reproduce the salient features of the colloidal glass experiment and support an effective description that is applicable to a broader class of viscoelastic suspensions.

Frictional active Brownian particles

Frictional forces affect the rheology of hard-sphere colloids, at high shear rate. Here we demonstrate, via numerical simulations, that they also affect the dynamics of active Brownian particles and their motility-induced phase separation. Frictional forces increase the angular diffusivity of the particles, in the dilute phase, and prevent colliding particles from resolving their collision by sliding one past to the other. This leads to qualitatively changes of motility-induced phase diagram in the volume-fraction motility plane. While frictionless systems become unstable towards phase separation as the motility increases only if their volume fraction overcomes a threshold, frictional systems become unstable regardless of their volume fraction. These results suggest the possibility of controlling the motility-induced phase diagram by tuning the roughness of the particles.

Granular packings with sliding, rolling, and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.