Yielding and flow in adhesive and nonadhesive concentrated emulsions

Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Detachment forces during parallel-plate gap separation mediated by a simple yield-stress fluid

In this work we have monitored the multiple stages of the normal traction force response of a yield-stress fluid confined between two circular parallel plates that are separated at constant velocity. At narrow initial gaps, the air-fluid interface suffers from the Saffman-Taylor instability, confirmed by visual inspection of fingering patterns imprinted on the fluid. At larger initial gaps, the fluid preserves the initially imposed circular symmetry of the confining plates, indicating the absence of instability. Due to the system characteristics and experimental environment, the multiple traction force contributions occurred in cascade, permitting us to isolate the adhesion responses associated with viscosity, capillarity, and yield stress. Employing a standard Herschel-Bulkley model, we assessed the scaling of the traction force in multiple regimes-specifically, evaluating the dependencies of the fingering to yield-stress transitions.

From shear bands to earthquakes in a model granular material with contact aging

We perform molecular dynamics simulations of homogeneous athermal systems of poly-disperse soft discs under shear. For purely repulsive interactions between particles, and under a confining external pressure, a monotonous flow curve (strain rate vs. stress) starting at a critical yield stress is obtained, with deformation distributing uniformly in the system, on average. Then we add a short range attractive contribution to the interaction potential that increases its intensity as particles remain in contact for a progressively longer time, mimicking an aging effect in the system. In this case the flow curve acquires a reentrant behavior, namely, a region where shear stress decreases with increasing strain. Within this region the deformation is seen to localize in a shear band with a well defined width that decreases as the global strain rate does. At very low strain rates the shear band becomes very thin and deformation acquires a prominent stick-slip behavior. This regime can be described as the system possessing a fault in which deformation occurs with an earthquake-resembling phenomenology. In this way the system we are analyzing connects a regime of uniform deformation at large strain rates, a localized deformation regime in the form of shear bands at intermediate stain rates, and seismic phenomena at very low strain rate. The unifying ingredient of this phenomenology is the existence of a reentrant flow curve, originating in the aging mechanisms present in the model.

Predicting yielding in attractive colloidal gels

One of the defining characteristics of soft glassy materials is their ability to exhibit a yield stress, which can result in an overall elasto-visco-plastic mechanics. To design soft materials with specific properties, it is essential to gain a comprehensive understanding of the topological and structural failure points that occur during yielding. However, predicting these failure points, which lead to yielding, is challenging due to the dynamic nature of structure development and its cooccurrence with other complicated processes, such as local rearrangements and anisotropy. In this study, we employ a series of tools from network science to investigate colloidal gels as a model for soft glassy materials during yielding. Our findings reveal that edge betweenness centrality can be utilized as a universal predictor for yielding across various state variables, including the volume fraction of solids, the strength, and the range of attraction between colloids.

Role of confinement in the shear banding and shear jamming in noncolloidal fiber suspensions

The jamming effect is critical in processing short fiber-reinforced thermoplastics (FRTs). Fiber jamming can induce discontinuous shear thickening (DST) in simple shear and result in fiber-matrix separation in more complex flows such as injection molding and compression molding of FRTs. The confinement effect commonly induces local jams and strongly enhances fiber jamming. However, the transient evolution of local fiber jams under confinement and its correlation with the tumbling of fibers are still elusive. In this study, we adopted rheo-PIV (particle image velocity) techniques to study this effect for glass fiber-reinforced thermoplastics (FRTs). The translational and tumbling motion of fiber were determined during rheological measurements, and the distribution of fiber orientation was determined by X-ray CT. Three shear banding regions appeared after the viscosity overshoot under high shear stress in suspensions with high fiber content, which was associated with the three regions of fiber orientation across the gap due to confinement. Shear banding was ascribed to the different tumbling speeds across the gap because of the different initial orientations and different wall confinements near and far from the wall. The local shear thickening and jamming behavior became most significant under intermediate confinement, and were affected by shear strain, shear stress, and fiber contents. 3D state diagrams were constructed to show the confinement effect on the evolution of shear banding and jamming.

Permanent shear localization in dense disordered materials due to microscopic inertia

In this work using computer simulations of 3D model of dense disordered solids we show, for the first time, the appearance of shear localization in the stationary flow under homogeneous driving conditions. To rationalize our simulation results we develop a continuum model, that couples the dynamics of the local flow to the evolution of a kinetic temperature field related to the local inertial dynamics. Our model predicts that the coupling of the flow field to this additional destabilizing field appears only as a necessary condition for shear localization, a minimum system size is necessary to accommodate the flow instability. Moreover we show that this size criterion resulting from our continuum description is in quantitative agreement with our particle-based simulation results.

Dynamic phase diagram of plastically deformed amorphous solids at finite temperature

The yielding transition that occurs in amorphous solids under athermal quasistatic deformation has been the subject of many theoretical and computational studies. Here, we extend this analysis to include thermal effects at finite shear rate, focusing on how temperature alters avalanches. We derive a nonequilibrium phase diagram capturing how temperature and strain rate effects compete, when avalanches overlap, and whether finite-size effects dominate over temperature effects. The predictions are tested through simulations of an elastoplastic model in two dimensions and in a mean-field approximation. We find a scaling for temperature-dependent softening in the low-strain rate regime when avalanches do not overlap, and a temperature-dependent Herschel-Bulkley exponent in the high-strain rate regime when avalanches do overlap.

Non-Invasive Rheo-MRI Study of Egg Yolk-Stabilized Emulsions: Yield Stress Decay and Protein Release

A comprehensive understanding of the time-dependent flow behavior of concentrated oil-in-water emulsions is of considerable industrial importance. Along with conventional rheology measurements, localized flow and structural information are key to gaining insight into the underlying mechanisms causing time variations upon constant shear. In this work, we study the time-dependent flow behavior of concentrated egg-yolk emulsions with (MEY) or without (EY) enzymatic modification and unravel the effects caused by viscous friction during shear. We observe that prolonged shear leads to irreversible and significant loss of apparent viscosity in both emulsion formulations at a mild shear rate. The latter effect is in fact related to a yield stress decay during constant shearing experiments, as indicated by the local flow curve measurements obtained by rheo-MRI. Concurrently, two-dimensional D-T₂ NMR measurements revealed a decrease in the T₂ NMR relaxation time of the aqueous phase, indicating the release of surface-active proteins from the droplet interface towards the continuous water phase. The combination of an increase in droplet diameter and the concomitant loss of proteins aggregates from the droplet interface leads to a slow decrease in yield stress.

Quantifying local rearrangements in three-dimensional granular materials: Rearrangement measures, correlations, and relationship to stresses

Quantifying the ways in which local particle rearrangements contribute to macroscopic plasticity is one of the fundamental pursuits of granular mechanics and soft matter physics. Here we examine local rearrangements that occur naturally during the deformation of three samples of 3D granular materials subjected to distinct boundary conditions by employing in situ x-ray measurements of particle-resolved structure and stress. We focus on five distinct rearrangement measures, their statistics, interrelationships, contributions to macroscopic deformation, repeatability, and dependence on local structure and stress. Our most significant findings are that local rearrangements (1) are correlated on a scale of three to four particle diameters, (2) exhibit volumetric strain-shear strain and nonaffine displacement-rotation coupling, (3) exhibit correlations that suggest either rearrangement repeatability or that rearrangements span multiple steps of incremental sample strain, and (4) show little dependence on local stress but correlate with quantities describing local structure, such as porosity. Our results are presented in the context of relevant plasticity theories and are consistent with recent findings suggesting that local structure may play at least as important of a role as local stress in determining the nature of local rearrangements.

Critical yielding rheology: from externally deformed glasses to active systems

We use extensive computer simulations to study the yielding transition under two different loading schemes: standard simple shear dynamics and self-propelled dense active systems. In the active systems, a yielding transition toward an out-of-equilibrium flowing state known as the liquid phase is observed when self-propulsion is increased. The range of self-propulsions in which this pure liquid regime exists appears to vanish upon approaching the so-called 'jamming point' at which the solidity of soft-sphere packings is lost. Such an 'active yielding' transition shares similarities with the generic yielding transition for shear flows. A Herschel-Bulkley law is observed along the liquid regime in both loading scenarios, with a clear difference in the critical scaling exponents between the two, suggesting the existence of different universality classes for the yielding transition under different driving conditions. In addition, we present the direct measurements of growing length and time scales for both driving scenarios. A comparison with theoretical predictions from the recent literature reveals poor agreement with our numerical results.

Continuum modeling of shear startup in soft glassy materials

Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a critical stress σ_{y} such that YSFs are solid for stresses σ imposed below σ_{y}, whereas they flow like liquids for σ>σ_{y}. Under an applied shear rate γ[over ̇], the solid-to-liquid transition is associated with a complex spatiotemporal scenario that depends on the microscopic details of the system, on the boundary conditions, and on the system size. Still, the general phenomenology reported in the literature boils down to a simple sequence that can be divided into a short-time response characterized by the so-called "stress overshoot," followed by stress relaxation towards a steady state. Such relaxation can be either (1) long-lasting, which usually involves the growth of a shear band that can be only transient or that may persist at steady state or (2) abrupt, in which case the solid-to-liquid transition resembles the failure of a brittle material, involving avalanches. In the present paper, we use a continuum model based on a spatially resolved fluidity approach to rationalize the complete scenario associated with the shear-induced yielding of YSFs. A key feature of our model is to provide a scaling for the coordinates of the stress overshoot, i.e., stress σ_{M} and strain γ_{M} as a function of γ[over ̇], which shows good agreement with experimental and numerical data extracted from the literature. Moreover, our approach shows that the power-law scaling σ_{M}(γ[over ̇]) is intimately linked to the growth dynamics of a fluidized boundary layer in the vicinity of the moving boundary. Yet such scaling is independent of the fate of that layer, and of the long-term behavior of the YSF, i.e., whether the steady-state flow profile is homogeneous or shear-banded. Finally, when including the presence of "long-range" correlations, we show that our model displays a ductile to brittle transition, i.e., the stress overshoot reduces into a sharp stress drop associated with avalanches, which impacts the scaling σ_{M}(γ[over ̇]). This generalized model nicely captures subtle avalanche-like features of the transient shear banding dynamics reported in experiments. Our work offers a unified picture of shear-induced yielding in YSFs, whose complex spatiotemporal dynamics are deeply connected to nonlocal effects.

Characterising shear-induced dynamics in flowing complex fluids using differential dynamic microscopy

Microscopic dynamics reveal the origin of the bulk rheological response in complex fluids. In model systems particle motion can be tracked, but for industrially relevant samples this is often impossible. Here we adapt differential dynamic microscopy (DDM) to study flowing highly-concentrated samples without particle resolution. By combining an investigation of oscillatory flow, using a novel "echo-DDM" analysis, and steady shear, through flow-DDM, we characterise the yielding of a silicone oil emulsion on both the microscopic and bulk level. Through measuring the rate of shear-induced droplet rearrangements and the flow velocity, the transition from a solid-like to liquid-like state is shown to occur in two steps: with droplet mobilisation marking the limit of linear visco-elasticity, followed by the development of shear localisation and macroscopic yielding. Using this suite of techniques, such insight could be developed for a wide variety of challenging complex fluids.

Instantaneous Normal Modes Reveal Structural Signatures for the Herschel-Bulkley Rheology in Sheared Glasses

The Herschel-Bulkley law, a universal constitutive relation, has been empirically known to be applicable to a vast range of soft materials, including sheared glasses. Although the Herschel-Bulkley law has attracted public attention, its structural origin has remained an open question. In this Letter, by means of atomistic simulation of binary Lennard-Jones glasses, we report that the instantaneous normal modes with negative eigenvalues, or so-called imaginary modes, serve as the structural signatures for the Herschel-Bulkley rheology in sheared glasses.

Stringiness of hyaluronic acid emulsions

Objective

Cosmetic emulsions containing hyaluronic acid are ubiquitous in the cosmetic industry. However, the addition of (different molecular weight) hyaluronic acid can affect the filament stretching properties of concentrated emulsions. This property is often related to the “stringiness” of an emulsion, which can affect the consumer's choice for a product. It is thus very important to investigate and predict the effect of hyaluronic acid on the filament stretching properties of cosmetic emulsions.

Methods

Model emulsions and emulsions with low and high molecular weights are prepared and their filament stretching properties are studied by the use of an extensional rheometer. Two different stretching speeds are employed during the stretching of the emulsions, a low speed at 10 µm/s and a high speed at 10 mm/s. The shear rheology of the samples is measured by rotational rheology.

Results

We find that filament formation only occurs at high stretching speeds when the emulsion contains high molecular weight hyaluronic acid. The formation of this filament, which happens at intermediate states of the break‐up, coincides with an exponential decay in the break‐up dynamics. The beginning and end of the break‐up of high molecular weight hyaluronic acid emulsions show a power law behaviour, where the exponent depends on the initial stretching rate. At a lower stretching speed, no filament is observed for both high molecular weight and low molecular weight hyaluronic acid emulsions and the model emulsion. The emulsions show a power law behaviour over the whole break‐up range, where the exponent also depends on the stretching rate. No significant difference is observed between the shear flow properties of the emulsions containing different molecular weights hyaluronic acid.

Conclusion

In this work, we underline the importance of the molecular weight of hyaluronic acid on the elongational properties of concentrated emulsions. The filament formation properties, for example the stringiness, of an emulsion is a key determinant of a product liking and repeat purchase. Here, we find that high molecular weight hyaluronic acid and a high stretching speed are the control parameters affecting the filament formation of an emulsion.

Unification of the Rheological Physics of Yield Stress Fluids

The physics above and below the yield stress is unified by a simple model for viscoplasticity that accounts for the nonlinear rheology of multiple yield stress fluids. The model has a rate-dependent relaxation time, allows for plastic deformation below the yield stress, and indicates that rapid elastic deformation aids yielding. A range of commonly observed rheological behaviors are predicted, including the smooth overshoot in the loss modulus and the recently discovered contributions from recoverable and unrecoverable strains in amplitude sweeps.

Interparticle attraction controls flow heterogeneity in calcite gels

We couple rheometry and ultrasonic velocimetry to study experimentally the flow behavior of gels of colloidal calcite particles dispersed in water, while tuning the strength of the interparticle attraction through physico-chemistry. We unveil, for the first time in a colloidal gel, a direct connection between attractive interactions and the occurrence of shear bands, as well as stress fluctuations.

Elucidating the G″ overshoot in soft materials with a yield transition via a time-resolved experimental strain decomposition

Materials that exhibit yielding behavior are used in many applications, from spreadable foods and cosmetics to direct write three-dimensional printing inks and filled rubbers. Their key design feature is the ability to transition behaviorally from solid to fluid under sufficient load or deformation. Despite its widespread applications, little is known about the dynamics of yielding in real processes, as the nonequilibrium nature of the transition impedes understanding. We demonstrate an iteratively punctuated rheological protocol that combines strain-controlled oscillatory shear with stress-controlled recovery tests. This technique provides an experimental decomposition of recoverable and unrecoverable strains, allowing for solid-like and fluid-like contributions to a yield stress material's behavior to be separated in a time-resolved manner. Using this protocol, we investigate the overshoot in loss modulus seen in materials that yield. We show that this phenomenon is caused by the transition from primarily solid-like, viscoelastic dissipation in the linear regime to primarily fluid-like, plastic flow at larger amplitudes. We compare and contrast this with a viscoelastic liquid with no yielding behavior, where the contribution to energy dissipation from viscous flow dominates over the entire range of amplitudes tested.

Computational study of transient shear banding in soft jammed solids

We have designed three-dimensional numerical simulations of a soft spheres model, with size polidispersity and in athermal conditions, to study the transient shear banding that occurs during yielding of jammed soft solids. We analyze the effects of different types of drag coefficients used in the simulations and compare the results obtained using Lees-Edwards periodic boundary conditions with the case in which the same model solid is confined between two walls. The specific damping mechanism and the different boundary conditions indeed modify the load curves and the velocity profiles in the transient regime. Nevertheless, we find that the presence of a stress overshoot and of a related transient banding phenomenon, for large enough samples, is a robust feature for overdamped systems, where their presence do not depend on the specific drag used and on the different boundary conditions.

Two step yielding in soft materials

A review is presented on the topic of two-step yielding observed in complex fluids that cover a broad variety of materials ranging from colloidal gels, attractive glasses, emulsions, suspensions, and several commercial paste-like materials. The common features in various systems displaying two-step yielding behavior are the presence of two characteristic forces between the interacting particles or two varying representative length or time scales. This focused review aims to provide physical insights, mechanistic understanding of the two-step yielding and other associated rheological consequences of this nonlinear behavior. A discussion is provided on the microstructural details with an overview of different experimental systems exhibiting double-yielding studied so far highlighting the similarities and differences among them. Particularly, the effects of continuous phase properties, dispersed particle phase factors (size, shape, softness and surface charge) and external force field (electric, magnetic, thermal and shear flows) on two-step yielding are considered.

Emergence and persistence of flow inhomogeneities in the yielding and fluidization of dense soft solids

In three-dimensional computer simulations of model non-Brownian jammed suspensions, we compute the time required to reach homogeneous flow upon yielding, by analyzing stresses and particle packing at different shear rates, with and without confinement. We show that the stress overshoot and persistent shear banding preceding the complete fluidization are controlled by the presence of overconstrained microscopic domains in the initial solids. Such domains, identifiable with icosahedrally packed regions in the model used, allow for stress accumulation during the shear startup. Their structural reorganization under deformation controls the emergence and the persistence of the shear banding.

Origin of the extremely high elasticity of bulk emulsions, stabilized by Yucca Schidigera saponins

We found experimentally that the elasticity of sunflower oil-in-water emulsions (SFO-in-W) stabilized by Yucca Schidigera Roezl saponin extract, is by >50 times higher as compared to the elasticity of common emulsions. We revealed that strong specific interactions between the phytosterols from the non-purified oil and the saponins from the Yucca extract lead to the formation of nanostructured adsorption layers which are responsible for the very high elasticity of the oil-water interface and of the respective bulk emulsions. Remarkably, this extra high emulsion elasticity inhibits the emulsion syneresis even at 65 vol% of the oil drops - these emulsions remain homogeneous and stable even after 30 days of shelf-storage. These results demonstrate that the combination of saponin and phytosterols is a powerful new approach to structure oil-in-water emulsions with potential applications for formulating healthier functional food.

Depletion attraction impairs the plasticity of emulsions flowing in a constriction

We study the elasto-plastic behavior of dense attractive emulsions under a mechanical perturbation. The attraction is introduced through non-specific depletion interactions between the droplets and is controlled by changing the concentration of surfactant micelles in the continuous phase. We find that such attractive forces are not sufficient to induce any measurable modification on the scalings between the local packing fraction and the deformation of the droplets. However, when the emulsions are flowed through 2D microfluidic constrictions, we uncover a measurable effect of attraction on their elasto-plastic response. Indeed, we measure higher levels of deformation inside the constriction for attractive droplets. In addition, we show that these measurements correlate with droplet rearrangements that are spatially delayed in the constriction for higher attraction forces.

Emulsions Stabilized by Inorganic Nanoclays and Surfactants: Stability, Viscosity, and Implications for Applications

Pickering emulsions, or emulsions with solid particles at the interface, have attracted significant interest in Enhanced Oil Recovery (EOR) processes, cosmetics, and drug delivery systems due to their ability to resist coalescence. Here, a synthetic clay nanoparticle, laponite®, is utilized to create oil-in-water (o/w) emulsions, and the addition of small-molecule surfactants induces a more stable emulsion. In this study, the stability of laponite® Pickering emulsions with and without the surfactants (dodecyltrimethylammonium bromide (DTAB), Pluronic F68 (F68), and sodium dodecyl sulfate (SDS) is investigated using dynamic light scattering (DLS), ζ-potential, optical microscopy, and rheology. With laponite® and no added surfactants, the DLS and ζ-potential results show formation of emulsion droplets with a diameter of 3 μm and a ζ-potential of -90 mV. With the addition of surfactants, both the droplet diameter and ζ-potential increase, suggesting adsorption of surfactants on the surface of laponite® particle. Optical microscopy suggests that the Pickering emulsion without surfactant undergoes flocculation, while the emulsion becomes stable to coalescence and creaming with addition of surfactants due to formation of a network structure. Regardless of the formation of network structure, the laponite®-F68 emulsion rheologically behaves as a Newtonian fluid, while the laponite®-SDS and laponite®-DTAB emulsions display shear thinning behavior. The difference in the rheological behavior can be attributed to the weak adsorption of F68 on laponite® and electrostatic interactions between laponite® and charged surfactants at oil-water interface.

Discrete fluidization of dense monodisperse emulsions in neutral wetting microchannels

The rheology of pressure-driven flows of two-dimensional dense monodisperse emulsions in neutral wetting microchannels is investigated by means of mesoscopic lattice Boltzmann simulations, capable of handling large collections of droplets, in the order of several hundreds. The simulations reveal that the fluidization of the emulsion proceeds through a sequence of discrete steps, characterized by yielding events whereby layers of droplets start rolling over each other, thus leading to sudden drops of the relative effective viscosity. It is shown that such discrete fluidization is robust against loss of confinement, namely it persists also in the regime of small ratios of the droplet diameter over the microchannel width. We also develop a simple phenomenological model which predicts a linear relation between the relative effective viscosity of the emulsion and the product of the confinement parameter (global size of the device over droplet radius) and the viscosity ratio between the disperse and continuous phases. The model shows excellent agreement with the numerical simulations. The present work offers new insights to enable the design of microfluidic scaffolds for tissue engineering applications and paves the way to detailed rheological studies of soft-glassy materials in complex geometries.

Role of interfacial elasticity for the rheological properties of saponin-stabilized emulsions

HYPOTHESIS

Saponins are natural surfactants which can provide highly viscoelastic interfaces. This property can be used to quantify precisely the effect of interfacial dilatational elasticity on the various rheological properties of bulk emulsions.

EXPERIMENTS

We measured the interfacial dilatational elasticity of adsorption layers from four saponins (Quillaja, Escin, Berry, Tea) adsorbed on hexadecane-water and sunflower oil-water interfaces. In parallel, the rheological properties under steady and oscillatory shear deformations were measured for bulk emulsions, stabilized by the same saponins (oil volume fraction between 75 and 85%).

FINDINGS

Quillaja saponin and Berry saponin formed solid adsorption layers (shells) on the SFO-water interface. As a consequence, the respective emulsions contained non-spherical drops. For the other systems, the interfacial elasticities varied between 2 mN/m and 500 mN/m. We found that this interfacial elasticity has very significant impact on the emulsion shear elasticity, moderate effect on the dynamic yield stress, and no effect on the viscous stress of the respective steadily sheared emulsions. The last conclusion is not trivial, because the dilatational surface viscoelasticity is known to have strong impact on the viscous stress of steadily sheared foams. Mechanistic explanations of all observed effects are described.

Criticality at a Finite Strain Rate in Fluidized Soft Glassy Materials

We study the emergence of critical dynamics in the steady shear rheology of fluidized soft glassy materials. Within a mesoscale elastoplastic model accounting for a shear band instability, we show how additional noise can induce a transition from a phase separated to homogeneous flow, accompanied by critical-like fluctuations of the macroscopic shear rate. Both macroscopic quantities and fluctuations exhibit power law behaviors in the vicinity of this transition, consistent with previous experimental findings on vibrated granular media. Altogether, our results suggest a generic scenario for the emergence of criticality when shear weakening mechanisms compete with a fluidizing noise.

Surface flow profiles for dry and wet granular materials by Particle Tracking Velocimetry; the effect of wall roughness

Two-dimensional Particle Tracking Velocimetry (PTV) is a promising technique to study the behaviour of granular flows. The aim is to experimentally determine the free surface width and position of the shear band from the velocity profile to validate simulations in a split-bottom shear cell geometry. The position and velocities of scattered tracer particles are tracked as they move with the bulk flow by analyzing images. We then use a new technique to extract the continuum velocity field, applying coarse-graining with the postprocessing toolbox MercuryCG on the discrete experimental PTV data. For intermediate filling heights, the dependence of the shear (or angular) velocity on the radial coordinate at the free surface is well fitted by an error function. From the error function, we get the width and the centre position of the shear band. We investigate the dependence of these shear band properties on filling height and rotation frequencies of the shear cell for dry glass beads for rough and smooth wall surfaces. For rough surfaces, the data agrees with the existing experimental results and theoretical scaling predictions. For smooth surfaces, particle-wall slippage is significant and the data deviates from the predictions. We further study the effect of cohesion on the shear band properties by using small amount of silicon oil and glycerol as interstitial liquids with the glass beads. While silicon oil does not lead to big changes, glycerol changes the shear band properties considerably. The shear band gets wider and is situated further inward with increasing liquid saturation, due to the correspondingly increasing trend of particles to stick together.

Sticky Matters: Jamming and Rigid Cluster Statistics with Attractive Particle Interactions

While the large majority of theoretical and numerical studies of the jamming transition consider athermal packings of purely repulsive spheres, real complex fluids and soft solids generically display attraction between particles. By studying the statistics of rigid clusters in simulations of soft particles with an attractive shell, we present evidence for two distinct jamming scenarios. Strongly attractive systems undergo a continuous transition in which rigid clusters grow and ultimately diverge in size at a critical packing fraction. Purely repulsive and weakly attractive systems jam via a first-order transition, with no growing cluster size. We further show that the weakly attractive scenario is a finite size effect, so that for any nonzero attraction strength, a sufficiently large system will fall in the strongly attractive universality class. We therefore expect attractive jamming to be generic in the laboratory and in nature.

On the universality of the flow properties of soft-particle glasses

We identify the minimal interparticle interactions necessary for a particle dynamics simulation to predict the structure and flow behaviour of soft particle glasses (SPGs). Generally, two kinds of forces between the particles must be accounted for in simulations of SPGs: viscous or frictional drag forces and elastic contact forces. Far field drag forces are required to dissipate energy in the simulations and capture the effect of the rheology of the suspending fluid. Elastic forces are found to be dominant compared to near-field drag or other forms of friction forces and are the most important component to compute the rheology. The shear stress, the first and second normal stress differences for different interparticle force laws collapse onto universal master curves of the Herschel-Bulkley form by non-dimensionalizing the stress with the yield stress and the shear rate with the viscosity of the suspending fluid divided by the low-frequency shear modulus. The Herschel-Bulkley exponents are close to 0.5 with a slight dependence on the repulsive pairwise elastic forces.

Analysis of network formation and long-term stability in silica nanoparticle stabilized emulsions

Emulsions are widely used in industrial applications, including in food sciences, cosmetics, and enhanced oil recovery. For these industries, an in depth understanding of the stability and rheological properties of emulsions under both static and dynamic conditions is vital to their successful application. Presented here is a thorough assessment of a model nanoparticle (NP) stabilized dodecane-in-water emulsion as a route to improved understanding of the relationship between NP properties, microstructure and droplet-droplet interactions on the stability and rheological properties of emulsions. Emulsions are obtained here with low NP loadings without the need for added electrolyte through the use of an optimized silica NP (SNP) surface modification procedure. The prepared emulsions were characterized via optical microscopy, cryo-scanning electron microscopy (cryo-SEM), zeta potential analysis and laser scanning confocal microscopy (LSCM), enabling quantification of the emulsion droplet size, SNP interfacial coverage/morphology and surface charge. The correlation of these properties with the rheology of the emulsions is investigated through small amplitude oscillatory shear experiments which provide significant insight into the origins of the emulsions' rheological behavior and their stability. In addition, long-term stability, droplet-droplet network formation and microstructural evolution are found to be readily detectable shortly after preparation through measured progression of the emulsion's rheological properties.

Microscopic processes controlling the Herschel-Bulkley exponent

The flow curve of various yield stress materials is singular as the strain rate vanishes and can be characterized by the so-called Herschel-Bulkley exponent n=1/β. A mean-field approximation due to Hebraud and Lequeux (HL) assumes mechanical noise to be Gaussian and leads to β=2 in rather good agreement with observations. Here we prove that the improved mean-field model where the mechanical noise has fat tails instead leads to β=1 with logarithmic correction. This result supports that HL is not a suitable explanation for the value of β, which is instead significantly affected by finite-dimensional effects. From considerations on elastoplastic models and on the limitation of speed at which avalanches of plasticity can propagate, we argue that β=1+1/(d-d_{f}), where d_{f} is the fractal dimension of avalanches and d the spatial dimension. Measurements of d_{f} then supports that β≈2.1 and β≈1.7 in two and three dimensions, respectively. We discuss theoretical arguments leading to approximations of β in finite dimensions.

Forming Sticky Droplets from Slippery Polymer Zwitterions

Polymer zwitterions are generally regarded as hydrophilic and repellant or "slippery" materials. Here, a case is described in which the polymer zwitterion structure is tailored to decrease water solubility, stabilize emulsion droplets, and promote interdroplet adhesion. Harnessing the upper critical solution temperature of sulfonium- and ammonium-based polymer zwitterions in water, adhesive droplets are prepared by adding organic solvent to an aqueous polymer solution at elevated temperature, followed by agitation to induce emulsification. Droplet aggregation is observed as the mixture cools. Variation of salt concentration, temperature, polymer concentration, and polymer structure modulates these interdroplet interactions, resulting in distinct changes in emulsion stability and fluidity. Under attractive conditions, emulsions encapsulating 50-75% oil undergo gelation. By contrast, emulsions prepared under conditions where droplets are nonadhesive remain fluid and, for oil fractions exceeding 0.6, coalescence is observed. The uniquely reactive nature of the selected zwitterions allows their in situ modification and affords a route to chemically trigger deaggregation and droplet dispersion. Finally, experiments performed in a microfluidic device, in which droplets are formed under conditions that either promote or suppress adhesion, confirm the salt-responsive character of these emulsions and the persistence of adhesive interdroplet interactions under flow.

Advances and challenges in the rheology of concentrated emulsions and nanoemulsions

We review advances that have been made in the rheology of concentrated emulsions and nanoemulsions, which can serve as model soft materials that have highly tunable viscoelastic properties at droplet volume fractions near and above the glass transition and jamming point. As revealed by experiments, simulations, and theoretical models, interfacial and positional structures of droplets can depend on the applied flow history and osmotic pressure that an emulsion has experienced, thereby influencing its key rheological properties such as viscoelastic moduli, yield stress and strain, and flow behavior. We emphasize studies of monodisperse droplets, since these have led to breakthroughs in the fundamental understanding of dispersed soft matter. This review also covers the rheological properties of attractive emulsions, which can exhibit a dominant elasticity even at droplet volume fractions far below maximal random jamming of hard spheres.

Flow-induced nanostructuring of gelled emulsions

Although the phase behavior of emulsions has been thoroughly investigated, the effect of flow on emulsion morphology, which is relevant for many applications, is far from being fully elucidated. Here, we investigate an emulsion based on two common nonionic surfactants in a range of water concentration where complex and diverse microstructures are found at rest, such as multilamellar and bicontinuous phases. In spite of such complexity, once subjected to shear flow, all the emulsions investigated are characterized by thinning filaments which eventually break up into a concentrated suspension of micro-sized water-based droplets dispersed in a continuous oil phase. The so-formed droplets tend to align in string-like structures. The emulsions exhibit a yield stress, whose value can be estimated by the plug-core velocity profiles in pressure-driven capillary flow, thus providing evidence of weakly attractive interdroplet interactions. The latter are consistent with droplet clustering and percolation observed at rest. These results can also be relevant to the flow behavior of other liquid-liquid systems, such as polymer blends, where the flow-induced microstructure is under debate as well.

Shear-induced slab-like domains in a directed percolated colloidal gel

We explore the structural changes of a gel-forming colloid polymer mixture under shear by employing Brownian dynamics simulations of a colloidal system with short-ranged attractive depletion interaction in a linear flow profile. While the structure of unpercolated systems changes only slightly under shearing, we discover the formation of slab-like clusters in sheared directed percolated gel networks that are confined between two walls. These gel-slabs are stable over a long time and seem to be related to the syneresis phenomena that can be observed in directed percolated colloidal gels. Only at large shear strength the slabs are destroyed and a homogeneous state with many unbounded particles can be observed. We also quantitatively analyze our results by determining void volumes.

Different universality classes at the yielding transition of amorphous systems

We study the yielding transition of a two-dimensional amorphous system under shear by using a mesoscopic elasto-plastic model. The model combines a full (tensorial) description of the elastic interactions in the system and the possibility of structural reaccommodations that are responsible for the plastic behavior. The possible structural reaccommodations are encoded in the form of a "plastic disorder" potential, which is chosen independently at each position of the sample to account for local heterogeneities. We observe that the stress must exceed a critical value σ_{c} in order for the system to yield. In addition, when the system yields a flow curve (relating stress σ and strain rate γ[over ̇]) of the form γ[over ̇]∼(σ-σ_{c})^{β} is obtained. Remarkably, we observe the value of β to depend on some details of the plastic disorder potential. For smooth potentials a value of β≃2.0 is obtained, whereas for potentials obtained as a concatenation of smooth pieces a value β≃1.5 is observed in the simulations. This indicates a dependence of critical behavior on details of the plastic behavior. In addition, by integrating out nonessential, harmonic degrees of freedom, we derive a simplified scalar version of the model that represents a collection of interacting Prandtl-Tomlinson particles. A mean-field treatment of this interaction reproduces the difference of β exponents for the two classes of plastic disorder potentials and provides values of β that compare favorably with those found in the full simulations.

Non-trivial rheological exponents in sheared yield stress fluids

In this work we discuss possible physical origins of non-trivial exponents in the athermal rheology of soft materials at low but finite driving rates. A key ingredient in our scenario is the presence of a self-consistent mechanical noise that stems from the spatial superposition of long-range elastic responses to localized plastically deforming regions. We study analytically a mean-field model, in which this mechanical noise is accounted for by a stress diffusion term coupled to the plastic activity. Within this description we show how a dependence of the shear modulus and/or the local relaxation time on the shear rate introduces corrections to the usual mean-field prediction, concerning the Herschel-Bulkley-type rheological response of exponent 1/2. This feature of the mean-field picture is then shown to be robust with respect to structural disorder and partial relaxation of the local stress. We test this prediction numerically on a mesoscopic lattice model that implements explicitly the long-range elastic response to localized shear transformations, and we conclude on how our scenario might be tested in rheological experiments.

Predicting and assessing rupture in protein gels under oscillatory shear

Soft materials may break irreversibly upon applying sufficiently large shear oscillations, a process whose physical mechanism remains largely elusive. In this work, the rupture of protein gels made of sodium caseinate under an oscillatory stress is shown to occur in an abrupt, brittle-like manner. Upon increasing the stress amplitude, the build-up of harmonic modes in the strain response can be rescaled for all gel concentrations. This rescaling yields an empirical criterion to predict the rupture point way before the samples are significantly damaged. "Fatigue" experiments under stress oscillations of constant amplitude can be mapped onto the former results, which indicates that rupture is independent of the temporal pathway in which strain and damage accumulate. Finally, using ultrasonic imaging, we measure the local mechanical properties of the gels before, during and after breakdown, showing that the strain field remains perfectly homogeneous up to rupture but suddenly gives way to a solid-fluid phase separation upon breakdown.

Understanding rheological hysteresis in soft glassy materials

Motivated by recent experimental studies of rheological hysteresis in soft glassy materials, we study numerically strain rate sweeps in simple yield stress fluids and viscosity bifurcating yield stress fluids. Our simulations of downward followed by upward strain rate sweeps, performed within fluidity models and the soft glassy rheology model, successfully capture the experimentally observed monotonic decrease of the area of the rheological hysteresis loop with sweep time in simple yield stress fluids, and the bell shaped dependence of hysteresis loop area on sweep time in viscosity bifurcating fluids. We provide arguments explaining these two different functional forms in terms of differing tendencies of simple and viscosity bifurcating fluids to form shear bands during the sweeps, and show that the banding behaviour captured by our simulations indeed agrees with that reported experimentally. We also discuss the difference in hysteresis behaviour between inelastic and viscoelastic fluids. Our simulations qualitatively agree with the experimental data discussed here for four different soft glassy materials.

Athermal rheology of weakly attractive soft particles

We study the rheology of a soft particulate system where the interparticle interactions are weakly attractive. Using extensive molecular dynamics simulations, we scan across a wide range of packing fractions (ϕ), attraction strengths (u), and imposed shear rates (γ[over ̇]). In striking contrast to repulsive systems, we find that at small shear rates generically a fragile isostatic solid is formed even if we go to ϕ≪ϕ_{J}. Further, with increasing shear rates, even at these low ϕ, nonmonotonic flow curves occur which lead to the formation of persistent shear bands in large enough systems. By tuning the damping parameter, we also show that inertia plays an important role in this process. Furthermore, we observe enhanced particle dynamics in the attraction-dominated regime as well as a pronounced anisotropy of velocity and diffusion constant, which we take as precursors to the formation of shear bands. At low enough ϕ, we also observe structural changes via the interplay of low shear rates and attraction with the formation of microclusters and voids. Finally, we characterize the properties of the emergent shear bands, and thereby, we find surprisingly small mobility of these bands, leading to prohibitively long time scales and extensive history effects in ramping experiments.

Shear Banding of Soft Glassy Materials in Large Amplitude Oscillatory Shear

We study shear banding in soft glassy materials subject to a large amplitude oscillatory shear flow (LAOS). By numerical simulations of the widely used soft glassy rheology model, supplemented by more general physical arguments, we demonstrate strong banding over an extensive range of amplitudes and frequencies of the imposed shear rate γ[over ˙](t)=γ[over ˙]_{0}cos(ωt), even in materials that do not permit banding as their steady state response to a steadily imposed shear flow γ[over ˙]=γ[over ˙]_{0}=const. Highly counterintuitively, banding persists in LAOS even in the limit of zero frequency ω→0, where one might a priori have expected a homogeneous flow response in a material that does not display banding under conditions of steadily imposed shear. We explain this finding in terms of an alternating competition within each cycle between glassy aging and flow rejuvenation. Our predictions have far-reaching implications for the flow behavior of aging yield stress fluids, suggesting a generic expectation of shear banding in flows of even arbitrarily slow time variation.

Attraction-induced jamming in the flow of foam through a channel

We study the flow of a pressure-driven foam through a straight channel using numerical simulations, and examine the effects of a tuneable attractive potential between bubbles. We show that the effect of an attractive potential is to introduce a regime of jamming and stick-slip flow in a channel, and report on the behaviour resulting from varying the strength of the attraction. We find that there is a force threshold below which the flow jams, and upon further increasing the driving force, a crossover from intermittent (stick-slip) to smooth flow is observed. This threshold force below which the foam jams increases linearly with the strength of the attractive potential. By examining the spectra of energy fluctuations, we show that stick-slip flow is characterized by low frequency rearrangements and strongly local behaviour, whereas steady flow shows a broad spectrum of energy drop events and collective behaviour. Our work suggests that the stick-slip and the jamming regimes occur due to the increased stabilization of contact networks by the attractive potential - as the strength of attraction is increased, bubbles are increasingly trapped within networks, and there is a decrease in the number of contact changes.

A flow visualization and superposition rheology study of shear-banding wormlike micelle solutions

In this paper, we use rheometry and flow visualization to study the dynamics of the interface between shear bands in a wormlike micellar solution sheared between concentric cylinders, i.e., in a Taylor-Couette (TC) cell, and to evaluate the stress diffusion coefficient and the stress correlation length in the Johnson-Segalman model. Two wormlike micellar solutions are studied: an aqueous solution of CTAB-NaNO3 and a solution of CPCl-NaSal in brine. These systems are highly elastic, exhibit Maxwellian behavior in linear viscoelasticity experiments, and shear banding in nonlinear experiments [S. Lerouge, et al., Soft Matter, 2008, 4, 1808-1819, M. A. Fardin, et al., Soft Matter, 2012, 8(39), 10072-10089, P. Ballesta, et al., J. Rheol., 2007, 51, 1047]. A large, custom-built, computer controlled TC cell allows us to rotate both cylinders independently and to visualize the flow in the r-z plane using a CCD camera. At low shear rates, the flow is stable and the fluid appears homogeneous throughout the gap between the cylinders. Above a critical shear rate, a shear banding transition occurs. This manifests itself in the formation of two distinct bands in the r-z plane, with an interface between the two bands. For sufficiently high ramp speeds, multiple steps of interface evolution are identified, as noted by Radulescu, Lerouge, and others [O. Redulescu, et al., Europhys. Lett., 2003, 62, 230, S. Lerouge, et al., Soft Matter, 2008, 4, 1808-1819]. We quantify the interface travel using direct visualization and use this measure, as well as superposition rheometry [P. Ballesta, et al., J. Rheol., 2007, 51, 1047], to determine the stress diffusion coefficient D and the stress correlation length ζ in the Johnson-Segalman model. These parameters are evaluated at different temperatures, shear rates, and gap sizes. We find that the stress diffusion coefficient and the stress correlation length exhibit a strong dependence on the gap of the Taylor-Couette cell for both shear-banding systems. For the CTAB-NaNO3 system, we report a linear dependence of the stress diffusion coefficient on temperature for the parameter range considered. In addition, we find that for this system, the stress diffusion coefficient is independent of shear rate. For the CPCl-NaSal system, we observe the same color changes in the sample reported by others on extended light exposure; however, we find that different histories of light exposure do not affect the measured stress diffusion coefficient.

Cooperativity flows and shear-bandings: a statistical field theory approach

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.

Probing relevant ingredients in mean-field approaches for the athermal rheology of yield stress materials

Although the notion of mechanical noise is expected to play a key role in the non-linear rheology of athermally sheared amorphous systems, its characterization has so far remained elusive. Here, we show using molecular dynamic simulations that in spite of the presence of strong spatio-temporal correlations in the system, the local stress exhibits normal diffusion under the effect of the mechanical noise in the finite driving regime. The diffusion constant appears to be proportional to the mean plastic activity. Our data suggests that the corresponding proportionality constant is density independent, and can be directly related to the specific form of the rheological flow curve, pointing the way to a generic way of modeling mechanical noise in mean-field equations.

Avalanche-size distributions in mean-field plastic yielding models

We discuss the size distribution N(S) of avalanches occurring at the yielding transition of mean-field (i.e., Hebraud-Lequeux) models of amorphous solids. The size distribution follows a power law dependence of the form N(S)∼S(-τ). However (contrary to what is found in its depinning counterpart), the value of τ depends on details of the dynamic protocol used. For random triggering of avalanches we recover the τ=3/2 exponent typical of mean-field models, which, in particular, is valid for the depinning case. However, for the physically relevant case of external loading through a quasistatic increase of applied strain, a smaller exponent (close to 1) is obtained. This result is rationalized by mapping the problem to an effective random walk in the presence of a moving absorbing boundary.

Physical chemistry of highly concentrated emulsions

This review explores the physics underlying the rheology of highly concentrated emulsions (HCEs) to determine the relationship between elasticity and HCE stability, and to consider whether it is possible to describe all physicochemical properties of HCEs on the basis of a unique physical approach. We define HCEs as emulsions with a volume fraction above the maximum closest packing fraction of monodisperse spheres, φm=0.74, even if droplets are not of polyhedron shape. The solid-like rheological behavior of HCEs is characterized by yield stress and elasticity, properties which depend on droplet polydispersity and which are affected by caging at volume fractions about the jamming concentration, φj. A bimodal size distribution in HCEs diminishes caging and facilitates droplet movement, resulting in HCEs with negligible yield stress and no plateau in storage modulus. Thermodynamic forces automatically move HCEs toward the lowest free energy state, but since interdroplet forces create local minimums - points beyond which free energy temporarily increases before it reaches the global minimum of the system - the free energy of HCEs will settle at a local minimum unless additional energy is added. Several attempts have been undertaken to predict the elasticity of HCEs. In many cases, the elastic modulus of HCEs is higher than the one predicted from classical models, which only take into account spatial repulsion (or simply interfacial energy). Improved models based on free energy calculation should be developed to consider the disjoining pressure and interfacial rheology in addition to spatial repulsion. The disjoining pressure and interfacial viscoelasticity, which result in the deviation of elasticity from the classical model, can be regarded as parameters for quantifying the stability of HCEs.