Microfluidic rheology of soft colloids above and below jamming

Brittle and ductile yielding in soft materials

Many soft materials yield under mechanical loading, but how this transition from solid-like behavior to liquid-like behavior occurs can vary significantly. Understanding the physics of yielding is of great interest for the behavior of biological, environmental, and industrial materials, including those used as inks in additive manufacturing and muds and soils. For some materials, the yielding transition is gradual, while others yield abruptly. We refer to these behaviors as being ductile and brittle. The key rheological signatures of brittle yielding include a stress overshoot in steady-shear-startup tests and a steep increase in the loss modulus during oscillatory amplitude sweeps. In this work, we show how this spectrum of yielding behaviors may be accounted for in a continuum model for yield stress materials by introducing a parameter we call the brittility factor. Physically, an increased brittility decreases the contribution of recoverable deformation to plastic deformation, which impacts the rate at which yielding occurs. The model predictions are successfully compared to results of different rheological protocols from a number of real yield stress fluids with different microstructures, indicating the general applicability of the phenomenon of brittility. Our study shows that the brittility of soft materials plays a critical role in determining the rate of the yielding transition and provides a simple tool for understanding its effects under various loading conditions.

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.

Revealing the characteristic length of random close packing via critical-like random pinning

By randomly pinning particles in fluidized states and finding the local energy minima, we form static packings of mono-disperse disks that resemble random close packing, when only n_c = 2.6% of the particles are pinned. The packings are isostatic and exhibit typical critical scalings of the jamming transition. The non-triviality of n_c is manifested mainly in two aspects. First, n_c acts as a critical point, leading to bifurcated critical scalings in its vicinity. The criticality of n_c is also demonstrated in the packings of weakly polydisperse disks. Second, n_c sets a length scale in agreement with the characteristic length of random close packing. With robust evidence, we show that this agreement is generally true for both mono- and poly-disperse particles and in both two and three dimensions. The randomness inherited from fluidized states by random pinning thus interprets the randomness of random close packing from a unique perspective.

Bulk rheology of sticky DNA-functionalized emulsions

We measure by experiment and particle-based simulation the rheology of concentrated, non-Brownian droplet emulsions functionalized with surface-bound single-stranded (ss), "sticky," DNA. In the absence of ssDNA, the emulsion viscosity increases with the dispersed phase volume fraction ϕ, before passing through a liquid-solid transition at a critical ϕ_{c} related to random close packing. Introducing ssDNA leads to a liquid-solid transition at ϕ<ϕ_{c}, the onset being set by the droplet valency N and the ssDNA concentration (or simulated binding strength ε). Using insight from simulation, we identify three key behaviors: (i) jammed suspensions (ϕ>ϕ_{c}≈0.64) show weak effects of functionalization, with elastic rheology instead governed by droplet stiffness; (ii) suspensions with ϕ<ϕ_{c} and N=1, 2 always exhibit viscous rheology, regardless of functionalization; and (iii) for ϕ<ϕ_{c} and N>3, functionalization leads to a controllable viscous-elastic transition. We present state diagrams showing the range of rheological tuning attainable by these means.

The structural, vibrational, and mechanical properties of jammed packings of deformable particles in three dimensions

We investigate the structural, vibrational, and mechanical properties of jammed packings of deformable particles with shape degrees of freedom in three dimensions (3D). Each 3D deformable particle is modeled as a surface-triangulated polyhedron, with spherical vertices whose positions are determined by a shape-energy function with terms that constrain the particle surface area, volume, and curvature, and prevent interparticle overlap. We show that jammed packings of deformable particles without bending energy possess low-frequency, quartic vibrational modes, whose number decreases with increasing asphericity and matches the number of missing contacts relative to the isostatic value. In contrast, jammed packings of deformable particles with non-zero bending energy are isostatic in 3D, with no quartic modes. We find that the contributions to the eigenmodes of the dynamical matrix from the shape degrees of freedom are significant over the full range of frequency and shape parameters for particles with zero bending energy. We further show that the ensemble-averaged shear modulus 〈G〉 scales with pressure P as 〈G〉 ∼ P^β, with β ≈ 0.75 for jammed packings of deformable particles with zero bending energy. In contrast, β ≈ 0.5 for packings of deformable particles with non-zero bending energy, which matches the value for jammed packings of soft, spherical particles with fixed shape. These studies underscore the importance of incorporating particle deformability and shape change when modeling the properties of jammed soft materials.

Glass-liquid and glass-gel transitions of soft-shell particles

We study the structure and dynamics of colloidal particles with a spherical hard core and a thermo-responsive soft shell over the whole phase diagram by means of small-angle x-ray scattering and x-ray photon correlation spectroscopy. By changing the effective volume fraction by temperature and particle concentration, liquid, repulsive glass. and attractive gel phases are observed. The dynamics slow down with increasing volume fraction in the liquid phase and reflect a Vogel-Fulcher-Tamann behavior known for fragile glass formers. We find a liquid-glass transition above 50 vol.% that is independent of the particles' concentration and temperature. In an overpacked state at effective volume fractions above 1, the dispersion does not show a liquid phase but undergoes a gel-glass transition at an effective volume fraction of 34 vol.%. At the same concentration, extrema of subdiffusive dynamics are found in the liquid phase at lower weight fractions. We interpret this as dynamic precursors of the glass-gel transition.

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.

Carbomer microgels as model yield-stress fluids

The review presents current research results for Carbopol-based microgels as yield-stress materials, covering three aspects: chemical, physical and rheological. Such a joint three-aspect study has no analog in the literature. The chemical aspects of Carbopol polymers are presented in terms of a cross-linking polymerization of acrylic acid, their molecular structure, microgel formulation, polyacid dissociation and neutralization, osmotic pressure and associated immense microgel swelling. The physical characterization is focused on models of the shear-induced solid-to-liquid transition of microgels, which are formed of mesoscopic particles typical for soft matter materials. Models that describe interparticle effects are presented to explain the energy states of microgel particles at the mesoscale of scrutiny. Typical representatives of the models utilize attributes of jamming dispersions, micromechanical and polyelectrolyte reactions. Selected relationships that result from the models, such as scaling rules and nondimensional flow characteristics are also presented. The rheological part presents the discussion of problems of yield stress in 2D and 3D deformations, appearance and magnitude of the wall slip. The theory and characteristics of Carbopol microgel deformation in rotational rheometers are presented with graphs for the steady-state measurements, stress-controlled oscillation and two types of transient shear deformation. The review is concluded with suggestions for future research.

Glass and Jamming Rheology in Soft Particles Made of PNIPAM and Polyacrylic Acid

The phase behaviour of soft colloids has attracted great attention due to the large variety of new phenomenologies emerging from their ability to pack at very high volume fractions. Here we report rheological measurements on interpenetrated polymer network microgels composed of poly(N-isopropylacrylamide) (PNIPAM) and polyacrylic acid (PAAc) at fixed PAAc content as a function of weight concentration. We found three different rheological regimes characteristic of three different states: a Newtonian shear-thinning fluid, an attractive glass characterized by a yield stress, and a jamming state. We discuss the possible molecular mechanisms driving the formation of these states.

Effect of D-Mannitol on the Microstructure and Rheology of Non-Aqueous Carbopol Microgels

D-mannitol is a common polyol that is used as additive in pharmaceutical and personal care product formulations. We investigated its effect on the microstructure and rheology of novel non-aqueous Carbopol dispersions employing traditional and time-resolved rheological analysis. We considered two types of sample, (i) fresh (i.e., mannitol completely dissolved in solution) and aged (i.e., visible in crystalline form). The analysis of the intracycle rheological transitions that were observed for different samples revealed that, when completely dissolved in solution, mannitol does not alter the rheological behaviour of the Carbopol dispersions. This highlights that the chemical similarity of the additive with the molecules of the surrounding solvent allows preserving the swollen dimension and interparticle interactions of the Carbopol molecules. Conversely, when crystals are present, a hierarchical structure forms, consisting of a small dispersed phase (Carbopol) agglomerated around a big dispersed phase (crystals). In keeping with this microstructural picture, as the concentration of Carbopol reduces, the local dynamics of the crystals gradually start to control the integrity of the microstructure. Rheologically, this results in a higher elasticity of the suspensions at infinitesimal deformations, but a fragile yielding process at intermediate strains.

Critical scaling of compression-driven jamming of athermal frictionless spheres in suspension

We study numerically a system of athermal, overdamped, frictionless spheres, as in a non-Brownian suspension, in two and three dimensions. Compressing the system isotropically at a fixed rate ε[over ̇], we investigate the critical behavior at the jamming transition. The finite compression rate introduces a control timescale, which allows one to probe the critical timescale associated with jamming. As was found previously for steady-state shear-driven jamming, we find for compression-driven jamming that pressure obeys a critical scaling relation as a function of packing fraction ϕ and compression rate ε[over ̇], and that the bulk viscosity p/ε[over ̇] diverges upon jamming. A scaling analysis determines the critical exponents associated with the compression-driven jamming transition. Our results suggest that stress-isotropic, compression-driven jamming may be in the same universality class as stress-anisotropic, shear-driven jamming.

Rheology finds distinct glass and jamming transitions in emulsions

We study the rheology of monodisperse and bidisperse emulsions with various droplet sizes (1-2 μm diameter). Above a critical volume fraction φc, these systems exhibit solid-like behavior and a yield stress can be detected. Previous experiments suggest that for small thermal particles, rheology will see a glass transition at φc = φg ≈ 0.58; for large athermal systems, rheology will see a jamming transition at φc = φJ ≈ 0.64. However, simulations point out that at the crossover of thermal and athermal regimes, the glass and jamming transitions may both be observed in the same sample. Here we conduct an experiment by shearing four oil-in-water emulsions with a rheometer. We observe both a glass and a jamming transition for our smaller diameter droplets, and only a jamming transition for our larger diameter droplets. The bidisperse sample behaves similarly to the small droplet sample, with two transitions observed. Our rheology data are well-fit by both the Herschel-Bulkley model and the three component model. Based on the fitting parameters, our raw rheological data would not collapse onto a master curve. Our results show that liquid-solid transitions in dispersions are not universal, but depend on particle size.

Investigation of the swollen state of Carbopol molecules in non-aqueous solvents through rheological characterization

We explore how different types of solvent influence the rheological properties of non-aqueous Carbopol dispersions from the dilute to the jammed state. In novel non-aqueous formulations, polar solvents are used more and more frequently, because they can form Carbopol microgels without the need of any neutralizing agents. However, the swelling behaviour of Carbopol molecules in the absence of water, when ionic forces are weak, is still poorly understood. To this end, we study the swelling behaviour of Carbopol 974P NF in different polar solvents, i.e. glycerol, PEG400 and mixtures of the two solvents, by mapping the rheological behaviour of Carbopol suspensions from very dilute to highly concentrated conditions. The rheological study reveals that the onset of the jamming transition occurs at different critical polymer concentrations depending on the solvents used. Nevertheless, once the jammed state is reached, both elastic and yielding behaviours are scalable with the particle volume fraction. These results suggest that the type of solvent influences the final volume of the single Carbopol particles but does not alter the interactions between the particles. The final radius of the swollen particles is estimated from shear rheology measurements in dilute conditions, showing a decrease of the final swelling ratio of Carbopol molecules of almost 50% for PEG400 solutions, a result that confirms the shift to higher values of the critical jamming concentration obtained from linear viscoelasticity for the same solutions.

Packing and dynamics of a protein solution approaching the jammed state

The packing of proteins and their collective behavior in crowded media is crucial for the understanding of biological processes. Here we study the structural and dynamical evolution of solutions of the globular protein bovine serum albumin with increasing concentration via drying using small angle X-ray scattering and dynamic light scattering. We probe an evolving correlation peak on the scattering profile, corresponding to the inter-protein distance, ξ, which decreases following a power law of the protein volume fraction, φ. The rate of decrease in ξ becomes faster above a protein concentration of ∼200 mg ml-1 (φ = 0.15). The power law exponent changes from 0.33, which is typical of colloidal or protein solutions, to 0.41. During the entire drying process, we observe the development and the growth of two-step relaxation dynamics with increasing φ as revealed by dynamic light scattering. We find three different regimes of the dependence of ξ as a function of φ. In the dilute regime (φ < 0.22), protein molecules are far apart from each other compared to their size. In this case, the dynamics mainly corresponds to Brownian motion. At an intermediate concentration (0.22 < φ < 0.47), inter-protein distances become comparable to the size of protein molecules, leading to a preferential orientation of the ellipsoidal protein molecules along with a possible deformation. In this regime, the dynamics shows two distinct relaxation times. At a very high concentration (φ > 0.47), the system reaches a jammed state. Subsequently, the secondary relaxation time in this state becomes extremely slow. In this state, the protein molecules have approximately one hydration layer. This study contributes to the understanding of protein molecular packing in crowded environments and the phenomenon of density-driven jamming for soft matter systems.

Soft granular particles sheared at a controlled volume: rate-dependent dynamics and the solid-fluid transition

We study the responses of fluid-immersed soft hydrogel spheres that are sheared under controlled volume fractions. Slippery, deformable particles along with the density-matched interstitial fluid are sandwiched between two opposing rough cones, allowing studies for a wide range of volume fraction φ both above and below the jamming of granular suspension. We utilize sudden cessations of shearing, accompanied by refraction-matched internal imaging, to supplement the conventional flow-curve measurements. At sufficiently high volume fractions, the settling of particles after the cessations exhibits a continuous yet distinct transition over the change of the shear rate. Such changes back out the qualitative difference in the state of flowing prior to the cessations: the quasi-static yielding of a tightly packed network, as opposed to the rapid sliding of particles mediated by the interstitial fluid whose dynamics depends on the driving rate. In addition, we determine the solid-fluid transition using two independent methods: the extrapolation of stress residues and the estimated yield stress from high values of φ, and the settling of particles upon shear cessations as φ goes across the transition. We also verify the power law on values of characteristic stress with respect to the distance from jamming φ - φc, with an exponent close to 2. These results demonstrate a multitude of relaxation timescales behind the dynamics of soft particles, and raise questions on how we extend the existing paradigms of the flow of a densely packed system when the softness is actively involved.

Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus G^{'} and the loss modulus G^{''} in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying G^{'}>0) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed G^{'} after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus G^{''} exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

Stress Relaxation above and below the Jamming Transition

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

Universal Relaxation Dynamics of Sphere Packings below Jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Viscosity Estimation of a Suspension with Rigid Spheres in Circular Microchannels Using Particle Tracking Velocimetry

Suspension flows are ubiquitous in industry and nature. Therefore, it is important to understand the rheological properties of a suspension. The key to understanding the mechanism of suspension rheology is considering changes in its microstructure. It is difficult to evaluate the influence of change in the microstructure on the rheological properties affected by the macroscopic flow field for non-colloidal particles. In this study, we propose a new method to evaluate the changes in both the microstructure and rheological properties of a suspension using particle tracking velocimetry (PTV) and a power-law fluid model. Dilute suspension (0.38%) flows with fluorescent particles in a microchannel with a circular cross section were measured under low Reynolds number conditions (Re ≈ 10^-4). Furthermore, the distribution of suspended particles in the radial direction was obtained from the measured images. Based on the power-law index and dependence of relative viscosity on the shear rate, we observed that the non-Newtonian properties of the suspension showed shear-thinning. This method will be useful in revealing the relationship between microstructural changes in a suspension and its rheology.

Nanogels and Microgels: From Model Colloids to Applications, Recent Developments, and Future Trends

Nanogels and microgels are soft, deformable, and penetrable objects with an internal gel-like structure that is swollen by the dispersing solvent. Their softness and the potential to respond to external stimuli like temperature, pressure, pH, ionic strength, and different analytes make them interesting as soft model systems in fundamental research as well as for a broad range of applications, in particular in the field of biological applications. Recent tremendous developments in their synthesis open access to systems with complex architectures and compositions allowing for tailoring microgels with specific properties. At the same time state-of-the-art theoretical and simulation approaches offer deeper understanding of the behavior and structure of nano- and microgels under external influences and confinement at interfaces or at high volume fractions. Developments in the experimental analysis of nano- and microgels have become particularly important for structural investigations covering a broad range of length scales relevant to the internal structure, the overall size and shape, and interparticle interactions in concentrated samples. Here we provide an overview of the state-of-the-art, recent developments as well as emerging trends in the field of nano- and microgels. The following aspects build the focus of our discussion: tailoring (multi)functionality through synthesis; the role in biological and biomedical applications; the structure and properties as a model system, e.g., for densely packed arrangements in bulk and at interfaces; as well as the theory and computer simulation.

From Yielding to Shear Jamming in a Cohesive Frictional Suspension

Simulations are used to study the steady shear rheology of dense suspensions of frictional particles exhibiting discontinuous shear thickening and shear jamming, in which finite-range cohesive interactions result in a yield stress. We develop a constitutive model that combines yielding behavior and shear thinning at low stress with the frictional shear thickening at high stresses, in good agreement with the simulation results. This work shows that there is a distinct difference between solids below the yield stress and in the shear-jammed state, as the two occur at widely separated stress levels, with an intermediate region of stress in which the material is flowable.

Crossover between Athermal Jamming and the Thermal Glass Transition of Suspensions

The non-Newtonian flow behavior of thermal and athermal disordered systems of dispersed uniform particles at high densities have strikingly similar features. By investigating the flow curves of yield-stress fluids and colloidal glasses having different volume fractions, particle sizes, and interactions, we show that both thermal and athermal systems exhibit power-law scaling with respect to the glass and jamming point, respectively, with the same exponents. All yield-stress flow curves can be scaled onto a single universal curve using the Laplace pressure as the stress scale for athermal systems and the osmotic pressure for the thermal systems. Strikingly, the details of interparticle interactions do not matter for the rescaling, showing that they are akin to usual phase transitions of the same universality class. The rescaling allows us to predict the flow properties of these systems from the volume fraction and known material properties.

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.

Critical scaling near the yielding transition in granular media

We show that the yielding transition in granular media displays second-order critical-point scaling behavior. We carry out discrete element simulations in the low-inertial-number limit for frictionless, purely repulsive spherical grains undergoing simple shear at fixed nondimensional shear stress Σ in two and three spatial dimensions. To find a mechanically stable (MS) packing that can support the applied Σ, isotropically prepared states with size L must undergo a total strain γ_{ms}(Σ,L). The number density of MS packings (∝γ_{ms}^{-1}) vanishes for Σ>Σ_{c}≈0.11 according to a critical scaling form with a length scale ξ∝|Σ-Σ_{c}|^{-ν}, where ν≈1.7-1.8. Above the yield stress (Σ>Σ_{c}), no MS packings that can support Σ exist in the large-system limit L/ξ≫1. MS packings generated via shear possess anisotropic force and contact networks, suggesting that Σ_{c} is associated with an upper limit in the degree to which these networks can be deformed away from those for isotropic packings.

Different routes into the glass state for soft thermo-sensitive colloids

We report an experimental and theoretical investigation of glass formation in soft thermo-sensitive colloids following two different routes: a gradual increase of the particle number density at constant temperature and an increase of the radius in a fixed volume at constant particle number density. Confocal microscopy experiments and the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory consistently show that the two routes lead to a dynamically comparable state at sufficiently long aging times. However, experiments reveal the presence of moderate but persistent structural differences. Successive cycles of radius decrease and increase lead instead to a reproducible glass state, indicating a suitable route to obtain rejuvenation without using shear fields.

Viscous dynamics of drops and bubbles in Hele-Shaw cells: Drainage, drag friction, coalescence, and bursting

In this review article, we discuss recent studies on drops and bubbles in Hele-Shaw cells, focusing on how scaling laws exhibit crossovers from the three-dimensional counterparts and focusing on topics in which viscosity plays an important role. By virtue of progresses in analytical theory and high-speed imaging, dynamics of drops and bubbles have actively been studied with the aid of scaling arguments. However, compared with three-dimensional problems, studies on the corresponding problems in Hele-Shaw cells are still limited. This review demonstrates that the effect of confinement in the Hele-Shaw cell introduces new physics allowing different scaling regimes to appear. For this purpose, we discuss various examples that are potentially important for industrial applications handling drops and bubbles in confined spaces by showing agreement between experiments and scaling theories. As a result, this review provides a collection of problems in hydrodynamics that may be analytically solved or that may be worth studying numerically in the near future.

Structure-property relationships from universal signatures of plasticity in disordered solids

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

Softening and yielding of soft glassy materials

Solids deform and fluids flow, but soft glassy materials, such as emulsions, foams, suspensions, and pastes, exhibit an intricate mix of solid- and liquid-like behavior. While much progress has been made to understand their elastic (small strain) and flow (infinite strain) properties, such understanding is lacking for the softening and yielding phenomena that connect these asymptotic regimes. Here we present a comprehensive framework for softening and yielding of soft glassy materials, based on extensive numerical simulations of oscillatory rheological tests, and show that two distinct scenarios unfold depending on the material's packing density. For dense systems, there is a single, pressure-independent strain where the elastic modulus drops and the particle motion becomes diffusive. In contrast, for weakly jammed systems, a two-step process arises: at an intermediate softening strain, the elastic and loss moduli both drop down and then reach a new plateau value, whereas the particle motion becomes diffusive at the distinctly larger yield strain. We show that softening is associated with an extensive number of microscopic contact changes leading to a non-analytic rheological signature. Moreover, the scaling of the softening strain with pressure suggest the existence of a novel pressure scale above which softening and yielding coincide, and we verify the existence of this crossover scale numerically. Our findings thus evidence the existence of two distinct classes of soft glassy materials - jamming dominated and dense - and show how these can be distinguished by their rheological fingerprint.

Design of yield-stress fluids: a rheology-to-structure inverse problem

We present a paradigm for the design of yield-stress fluids, using six archetypal materials for demonstration. By applying concepts of engineering design, we outline a materials design paradigm that includes (i) morphological organization based on jammed versus networked microstructures, (ii) collected scaling laws for predictive design, (iii) low-dimensional descriptions of function-valued flow data, (iv) consideration of secondary properties including viscous behavior, and (v) a strategy for material concept synthesis based on the juxtaposition of microstructures. By explicitly specifying these design strategies, we seek to create an ontology and database for the engineering of yield-stress fluids. Our proposed design strategy increases the likelihood of finding an optimal material and prevents design fixation by considering multiple material classes to achieve a desired rheological performance. This flips the typical structure-to-rheology analysis to become the inverse: rheology-to-structure with multiple possible materials as solutions.

Microsecond Structural Rheology

The relationship between the local structure of complex liquids and their response to shear is generally not well understood. This concerns, in particular, the formation of particle strings in the flow direction or hydroclusters, both important for the understanding of shear thinning and thickening phenomena. Here, we present results of a microfocus X-ray scattering experiment on spherical silica colloids in a liquid jet at high shear rates. Along and across the jet, we observe direction-dependent modifications of the structure factor of the suspension, suggesting the formation of differently ordered clusters in compression lines and as particle strings. With increasing distance from the orifice, the structure relaxes to the unsheared case with a typical relaxation 10 times larger as the time scale of Brownian motion. These results provide the first experimental flow characterization of a complex fluid at high shear rates detecting cluster formation and relaxation with micrometer and microsecond resolution.

Discontinuous change of shear modulus for frictional jammed granular materials

The shear modulus of jammed frictional granular materials with harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude γ_{0} discontinuously emerges at the jamming transition point. The storage modulus for small γ_{0} differs from that of frictionless grains even in the zero friction limit, whereas they are almost identical with each other for sufficiently large γ_{0}, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for larger strain. We propose a scaling law to interpolate between the states of small and large γ_{0}.

Slow creep in soft granular packings

Transient creep mechanisms in soft granular packings are studied numerically using a constant pressure and constant stress simulation method. Rapid compression followed by slow dilation is predicted on the basis of a logarithmic creep phenomenon. Characteristic scales of creep strain and time exhibit a power-law dependence on jamming pressure, and they diverge at the jamming point. Microscopic analysis indicates the existence of a correlation between rheology and nonaffine fluctuations. Localized regions of large strain appear during creep and grow in magnitude and size at short times. At long times, the spatial structure of highly correlated local deformation becomes time-invariant. Finally, a microscale connection between local rheology and local fluctuations is demonstrated in the form of a linear scaling between granular fluidity and nonaffine velocity.

Scaling Theory for the Frictionless Unjamming Transition

We develop a scaling theory of the unjamming transition of soft frictionless disks in two dimensions by defining local areas, which can be uniquely assigned to each contact. These serve to define local order parameters, whose distribution exhibits divergences as the unjamming transition is approached. We derive scaling forms for these divergences from a mean-field approach that treats the local areas as noninteracting entities, and demonstrate that these results agree remarkably well with numerical simulations. We find that the asymptotic behavior of the scaling functions arises from the geometrical structure of the packing while the overall scaling with the compression energy depends on the force law. We use the scaling forms of the distributions to determine the scaling of the total grain area A_{G} and the total number of contacts N_{C}.

Structural and cooperative length scales in polymer gels

Understanding the relationship between the material structural details, the geometrical confining constraints, the local dynamical events and the global rheological response is at the core of present investigations on complex fluid properties. In the present article, this problem is addressed on a model yield stress fluid made of highly entangled polymer gels of Carbopol which follows at the macroscopic scale the well-known Herschel-Bulkley rheological law. First, performing local rheology measurements up to high shear rates ([Formula: see text] s^-1)and under confinement, we evidence unambiguously the breakdown of bulk rheology associated with cooperative processes under flow. Moreover, we show that these behaviors are fully captured with a unique cooperativity length [Formula: see text] over the whole range of experimental conditions. Second, we introduce an original optical microscopy method to access structural properties of the entangled polymer gel in the direct space. Performing image correlation spectroscopy of fluorophore-loaded gels, the characteristic size D of carbopol gels microstructure is determined as a function of preparation protocol. Combining both dynamical and structural information shows that the measured cooperative length [Formula: see text] corresponds to 2-5 times the underlying structural size D, thus providing a strong grounding to the "Shear Transformation Zones" modeling approach.

Yielding of glass under shear: A directed percolation transition precedes shear-band formation

Under external mechanical loading, glassy materials, ranging from soft matter systems to metallic alloys, often respond via formation of inhomogeneous flow patterns, during yielding. These inhomogeneities can be precursors to catastrophic failure, implying that a better understanding of their underlying mechanisms could lead to the design of smarter materials. Here, extensive molecular dynamics simulations are used to reveal the emergence of heterogeneous dynamics in a binary Lennard-Jones glass, subjected to a constant strain rate. At a critical strain, this system exhibits for all considered strain rates a transition towards the formation of a percolating cluster of mobile regions. We give evidence that this transition belongs to the universality class of directed percolation. Only at low shear rates, the percolating cluster subsequently evolves into a transient (but long-lived) shear band with a diffusive growth of its width. Finally, the steady state with a homogeneous flow pattern is reached. In the steady state, percolation transitions also do occur constantly, albeit over smaller strain intervals, to maintain the stationary plastic flow in the system.

Scaling crossover in thin-film drag dynamics of fluid drops in the Hele-Shaw cell

We study both experimentally and theoretically the descending motion due to gravity of a fluid drop surrounded by another immiscible fluid in a confined space between two parallel plates, i.e., in the Hele-Shaw cell. As a result, we show a new scaling regime of a nonlinear drag friction in viscous liquid that replaces the well-known Stokes’ drag friction through a clear collapse of experimental data thanks to the scaling law. In the novel regime, the dissipation in the liquid thin film formed between the drop and cell walls governs the dynamics. The crossover of this scaling regime to another scaling regime in which the dissipation inside the droplet is dominant is clearly demonstrated and a phase diagram separating these scaling regimes is presented.

Effect of particle collisions in dense suspension flows

We study nonlocal effects associated with particle collisions in dense suspension flows, in the context of the Affine Solvent Model, known to capture various aspects of the jamming transition. We show that an individual collision changes significantly the velocity field on a characteristic volume Ω_{c}∼1/δz that diverges as jamming is approached, where δz is the deficit in coordination number required to jam the system. Such an event also affects the contact forces between particles on that same volume Ω_{c}, but this change is modest in relative terms, of order f_{coll}∼f[over ¯]^{0.8}, where f[over ¯] is the typical contact force scale. We then show that the requirement that coordination is stationary (such that a collision has a finite probability to open one contact elsewhere in the system) yields the scaling of the viscosity (or equivalently the viscous number) with coordination deficit δz. The same scaling result was derived [E. DeGiuli, G. Düring, E. Lerner, and M. Wyart, Phys. Rev. E 91, 062206 (2015)PLEEE81539-375510.1103/PhysRevE.91.062206] via different arguments making an additional assumption. The present approach gives a mechanistic justification as to why the correct finite size scaling volume behaves as 1/δz and can be used to recover a marginality condition known to characterize the distributions of contact forces and gaps in jammed packings.

Microfluidic viscometers for shear rheology of complex fluids and biofluids

The rich diversity of man-made complex fluids and naturally occurring biofluids is opening up new opportunities for investigating their flow behavior and characterizing their rheological properties. Steady shear viscosity is undoubtedly the most widely characterized material property of these fluids. Although widely adopted, macroscale rheometers are limited by sample volumes, access to high shear rates, hydrodynamic instabilities, and interfacial artifacts. Currently, microfluidic devices are capable of handling low sample volumes, providing precision control of flow and channel geometry, enabling a high degree of multiplexing and automation, and integrating flow visualization and optical techniques. These intrinsic advantages of microfluidics have made it especially suitable for the steady shear rheology of complex fluids. In this paper, we review the use of microfluidics for conducting shear viscometry of complex fluids and biofluids with a focus on viscosity curves as a function of shear rate. We discuss the physical principles underlying different microfluidic viscometers, their unique features and limits of operation. This compilation of technological options will potentially serve in promoting the benefits of microfluidic viscometry along with evincing further interest and research in this area. We intend that this review will aid researchers handling and studying complex fluids in selecting and adopting microfluidic viscometers based on their needs. We conclude with challenges and future directions in microfluidic rheometry of complex fluids and biofluids.