Continuum modeling of Soft Glassy Materials under shear

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a jar of thick, sticky toothpaste, a glob of hair gel, or a bucket of dense paint. These are Soft Glassy Materials (SGMs). At rest, they act like solids; you can poke them, and they hold their shape. But if you stir them hard enough, they suddenly turn into liquids and flow.

This paper is about understanding exactly how and when that switch happens, and why the process is often messy, uneven, and full of surprises.

Here is the story of the paper, broken down into simple concepts and analogies.

1. The Problem: The "Stress Overshoot"

Imagine you are trying to push a heavy, stuck sofa across a carpet.

The Push: You start pushing gently. The sofa doesn't move, but you feel the resistance building up in your arms. This is the elastic response (like a spring stretching).
The Snap: Suddenly, the sofa jerks forward. For a split second, the resistance you feel spikes to a maximum before dropping as the sofa starts sliding easily.
The Analogy: In the world of these materials, this spike is called a Stress Overshoot. It's the moment the material "breaks" its own internal structure to start flowing.

The researchers wanted to know: How does the size of this spike change if you push faster or slower? And how long does it take for the whole jar to turn from a solid to a liquid?

2. The Solution: A New "Traffic Map"

To solve this, the authors created a mathematical model. Instead of tracking every single tiny particle (which would be like counting every grain of sand on a beach), they used a concept called Fluidity.

The Analogy: Think of the material as a crowded dance floor.
- Solid State: Everyone is frozen in place, holding hands.
- Fluid State: Everyone is dancing wildly.
- Fluidity: This is a measure of how "wiggly" the crowd is at any specific spot. High fluidity = dancing; Low fluidity = frozen.

The model treats the material like a fluid map where "wiggly-ness" spreads from the edges inward, like a ripple in a pond.

3. The "Shear Band": The Unfair Breakup

When you start stirring the material, it doesn't always melt evenly.

The Phenomenon: Often, the material near the stirring spoon (the wall) melts first, creating a river of liquid. But the stuff in the middle stays frozen. This is called Shear Banding.
The Analogy: Imagine a crowd of people trying to exit a stadium. The people right next to the exit door start running (fluid), but the people in the middle of the stands are still sitting down (solid). The "river" of runners slowly eats its way through the crowd until everyone is moving.

The paper's model successfully predicts how fast this "river" grows and how long it takes for the whole crowd to start moving (the Fluidization Time).

4. The Secret Ingredient: "Cooperativity"

Why doesn't the whole thing melt instantly? Why does it happen in a wave?
The model introduces a concept called Cooperativity.

The Analogy: Think of a game of "Dominoes" or a "Wave" in a stadium. One person standing up (a plastic event) makes their neighbors more likely to stand up too.
The model says that when one part of the material breaks, it "shouts" to its neighbors, helping them break too. This "shout" travels a short distance. The model uses this to explain why the melting front is sharp and moves at a specific speed.

5. The Twist: The "Slippery Wall" Effect

The authors realized that in real life, the walls of the container matter. If the container is smooth, the soft particles might slide against the wall like ice on ice. If it's rough, they get stuck.

The Analogy: Imagine pushing a block of jelly on a table.
- Rough Table: The jelly sticks to the table, and you have to push the whole block to move it.
- Wet Table: The jelly slides easily. You might push the top, but the bottom just slides away without the whole thing moving together.

The paper updated their model to include this "slippery" effect (called Elasto-Hydrodynamic interactions). They found that when the walls are very smooth, the rules for how the material breaks change completely. The "Stress Overshoot" behaves differently, following a new set of math rules that depend on how slippery the wall is.

6. The Big Picture: Why This Matters

This paper is a "pedagogical introduction," meaning it's a guidebook for other scientists. It shows that:

We can predict the chaos: Even though these materials look messy, their behavior follows strict mathematical laws.
History matters: How you push the material (fast vs. slow) changes how it breaks.
Boundaries are key: What happens at the edge of the container (the wall) controls what happens in the middle.

In Summary:
The authors built a smart "traffic simulator" for soft materials. It explains why toothpaste sometimes shoots out of the tube in a sudden burst, why paint might flow unevenly, and how to predict exactly when a solid will turn into a liquid. By understanding these rules, engineers can design better 3D printers, adhesives, and food products that behave exactly the way we want them to.

1. Problem Statement

Soft Glassy Materials (SGMs), such as colloidal glasses, foams, and emulsions, exhibit complex rheological behaviors when subjected to shear. Key phenomena include:

Shear-induced yielding: A transition from a solid-like state to a fluid-like state.
Stress overshoot: A non-monotonic stress response where stress rises to a maximum ( $\sigma_M$ ) before relaxing to a steady state.
Transient shear banding: The formation of spatially heterogeneous flows where a fluidized region (shear band) coexists with a solid-like arrested region.
Fluidization time: The time required for the entire sample to transition to a homogeneous flowing state.

Existing models often struggle to quantitatively capture the interplay between these transient phenomena, the dependence on shear rate, and the influence of boundary conditions (e.g., wall slip) within a unified continuum framework. The authors aim to provide a pedagogical introduction to a non-local fluidity model that rationalizes these complex flows and extends it to include elasto-hydrodynamic (EHD) interactions.

2. Methodology

The authors employ a continuum modeling approach based on the concept of fluidity ( $f$ ), defined as the local rate of plastic events. The model is built upon the following pillars:

Order Parameter: The fluidity $f(y,t)$ serves as the order parameter, where $f=0$ represents a solid-like (arrested) state and $f>0$ represents a fluidized state.
Non-local Effects: The model incorporates a cooperativity length scale ( $\xi$ ), which accounts for the spatial propagation of plastic rearrangements. This is derived from the continuum limit of microscopic equations (Hébraud-Lequeux model).
Governing Equations:
1. Fluidity Dynamics: The evolution of fluidity is governed by a gradient flow equation derived from a free energy functional $F[f]$ :
  $\frac{\partial f}{\partial \tilde{t}} = f \left( \xi^2 \Delta f + m f - f^{3/2} \right)$
  Here, $\tilde{t}$ is a rescaled time, and $m$ depends on the stress $\Sigma$ . The term $\kappa[f] \sim f$ (mobility) is crucial; it ensures that a solid-like region ( $f \approx 0$ ) is unstable and will eventually fluidize, allowing for transient shear banding.
2. Stress Evolution (Maxwell Model): The stress $\Sigma$ evolves based on the total strain rate, decomposed into elastic and plastic parts:
  $\frac{d\Sigma}{d\tilde{t}} = \frac{1}{\tau} \left[ \dot{\Gamma} - \langle f \rangle \Sigma \right]$
  where $\langle f \rangle$ is the spatial average of fluidity.
3. Boundary Conditions: The model imposes specific boundary conditions at the moving wall ( $y=0$ ). A "wall fluidity" $f_w = m^2(\Sigma)$ is imposed when stress exceeds the yield stress, driving the nucleation of the shear band.
Extension to EHD Interactions: To address dense assemblies of deformable particles (e.g., microgels), the authors modify the Maxwell equation to include an EHD strain contribution ( $\Gamma_{EHD}$ ). This introduces a term proportional to $\dot{\Gamma}_0 \Sigma^2$ , representing lubrication flows between particles.

3. Key Contributions

Unified Framework for Transient Flows: The paper demonstrates that a non-local fluidity model can quantitatively capture the entire shear start-up process, from the initial elastic response to the stress overshoot, transient shear banding, and final homogeneous flow.
Scaling Laws for Stress Overshoot: The authors derive analytical scaling laws for the stress overshoot magnitude ( $\Sigma_M$ $Σ_{M}$ ) as a function of the applied shear rate ( $\dot{\Gamma}$ $\dot{Γ}$ ). They identify two distinct regimes:
- Diffusive regime (low $\dot{\Gamma}$ ): $\Sigma_M - 1 \sim \dot{\Gamma}^{\beta}$ with $\beta = 2n/3$ .
- Asymptotic regime (high $\dot{\Gamma}$ ): $\Sigma_M - 1 \sim \dot{\Gamma}^{\alpha}$ with $\alpha = 4n/(9-n)$ .
  These predictions match experimental data on Carbopol microgels remarkably well.
Incorporation of EHD Effects: The model is successfully extended to include EHD interactions. This modification explains the "kink" observed in steady-state flow curves for smooth walls (indicating wall slip) and alters the stress overshoot scaling at low shear rates to $\Sigma_M \sim \dot{\Gamma}^{1/2}$ , independent of the Herschel-Bulkley exponent $n$ .
Fluidization Time Prediction: The model predicts the fluidization time ( $T_f$ ) scaling as $T_f \sim \dot{\Gamma}^{-9/4}$ for imposed shear rate, a result consistent with experimental observations and independent of the shear-thinning index $n$ .

4. Results

Stress Overshoot: Numerical simulations of the model reproduce the characteristic stress overshoot seen in experiments. The model successfully predicts the transition between the "unsteady" regime (where data deviates from the steady-state Herschel-Bulkley curve) and the "quasi-steady" regime (where data collapses onto the curve).
Shear Banding: The model captures the growth of a shear band ( $\ell_b$ ) from the moving wall. The interface between the fluidized band and the solid core remains sharp, and the band growth dynamics are governed by the cooperativity length $\xi$ .
EHD Impact: When EHD interactions are included, the model reproduces the deviation from the standard Herschel-Bulkley law at low shear rates for smooth surfaces. It also correctly predicts that the stress overshoot scaling exponent changes from the HB-dependent value to $1/2$ when EHD effects dominate.
Boundary Condition Sensitivity: The study highlights that the nature of fluidization (ductile vs. brittle) is highly sensitive to boundary conditions. Fixing the fluidity gradient to zero at the wall (instead of fixing the fluidity value) can suppress shear band nucleation, leading to brittle-like failure characterized by an abrupt stress drop.

5. Significance

This work provides a versatile and quantitative theoretical framework for understanding the yielding transition in Soft Glassy Materials.

Predictive Power: It offers specific, testable scaling laws for rheological variables (stress overshoot, fluidization time) that depend on material parameters (Herschel-Bulkley exponent $n$ , cooperativity length $\xi$ ).
Bridging Scales: By linking microscopic plastic event rates to macroscopic continuum fields via non-local terms, it bridges the gap between particle-level simulations and macroscopic rheology.
Practical Relevance: The inclusion of EHD interactions and the analysis of boundary conditions are critical for engineering applications involving SGMs (e.g., additive manufacturing, coating), where wall slip and particle deformability significantly alter flow behavior.
Open Questions: The paper identifies future directions, including the application of the model to time-dependent protocols (shear ramps) to study hysteresis, and the need for deeper physical insight into the microscopic mechanisms governing boundary conditions.