Constitutive relations for colloidal gel

This paper challenges traditional continuum theories that assume a stress-free reference state for colloidal gels and proposes new, more accurate constitutive relations validated through large-scale numerical simulations of depletion and frictional gels.

The Gist

The Mystery of the "Stubborn Jelly": Why Old Science Couldn't Explain Gels

Imagine you are trying to predict how a giant block of Jell-O will behave when you squeeze it. For decades, scientists used a standard "recipe" (mathematical formulas) to predict this. This recipe assumed that if you squeeze the jelly, every tiny part of it moves in a perfectly predictable, uniform way—like a disciplined marching band where every soldier takes the exact same step at the exact same time.

But there’s a problem: Colloidal gels aren't marching bands; they are more like a chaotic crowd at a music festival.

This paper, written by researchers at IIT Ropar, explains why our old "marching band" math fails for gels and introduces a new way to understand them.

---

1. The Problem: The "Ghost" in the Machine

Traditional science assumes that before you touch a material, it is "relaxed" and "stress-free"—like a calm lake.

However, the authors point out that gels are born stressed. Because of how they form (particles clumping together quickly), they are full of "quenched stresses." Imagine building a Lego castle, but halfway through, you realize some bricks are being pushed inward and others are being pulled outward by invisible rubber bands. Even before you touch the castle, it is already under tension and compression.

Because of this internal "tug-of-war," when you squeeze the gel, it doesn't react uniformly. Some parts resist, some parts give way, and some parts push sideways in ways that old math says are impossible.

2. The Failed Prediction: The "Constant Ratio" Error

The old theory (by a scientist named Zaccone) predicted that if you squeeze a gel, the ratio of how it pushes down versus how it pushes outward should always be a constant number.

It’s like saying, "No matter how hard you squeeze this sponge, it will always bulge out exactly 2 inches."

But in real life, experiments showed this wasn't true. As the gel gets denser, the way it pushes sideways changes drastically. The old math was blind to the "personality" of the gel.

3. The New Solution: The "Force Skeleton"

Instead of looking at the gel as a uniform block, the researchers looked at the "Force Skeleton."

Think of a gel like a complex web of spiderwebs. Even if the web looks the same in every direction (isotropic), the strength of the silk threads isn't the same everywhere. When you squeeze the web, the force travels along specific "highways" of strength.

The researchers discovered something fascinating:

• The Structure is Boring: The actual physical arrangement of the particles stays mostly the same (it looks like a random, messy web).
• The Force is Exciting: Even though the web looks random, the forces traveling through it are highly organized. When you squeeze the gel, the particles create "force chains"—invisible highways of intense pressure that align with the squeeze.

4. The "New Recipe"

The authors proposed a new set of formulas. Instead of focusing on how much the material stretches (strain), they focused on how the internal force network organizes itself.

They used advanced math (spherical harmonics—think of these as describing shapes on a ball) to map out how forces are distributed. Their new model doesn't just look at the "marching band"; it looks at the "highways of force" running through the crowd.

Why does this matter?

Gels aren't just for dessert. They are used to make:

• Ceramics (for high-tech tools)
• Cosmetics (the texture of your lotion)
• Paints and Drilling Fluids (for oil and gas)

By understanding the "Force Skeleton," engineers can better predict how these materials will behave during manufacturing, preventing cracks, collapses, or unexpected textures.

In short: The researchers moved from treating gels like a simple, predictable block to treating them like a complex, living network of invisible force-highways.

Technical Summary

Technical Summary: Constitutive Relations for Colloidal Gel

Authors: Saikat Roy and Yezaz Ahmed Gadi Man (IIT Ropar)

1. Problem Statement

Traditional continuum theories for describing the mechanical response of colloidal gels (specifically depletion gels) rely on two fundamental assumptions:

1. Stress-free reference state: The theory assumes the material can be expanded around an unstressed state.
2. Affine deformation: The theory assumes that microscopic particle displacements follow the macroscopic strain field.

The authors argue that these assumptions fail for amorphous gel materials. Gels are inherently heterogeneous, containing large voids and a "quenched" internal stress field (local regions of tension and compression) that exists from the moment of formation. Consequently, the strain field is non-affine. Previous models (e.g., Zaccone et al.) fail to predict critical experimental observations, such as the Normal Stress Difference (NSD) during uniaxial consolidation and the specific power-law growth of stress with volume fraction. Specifically, existing theories predict a constant, non-evolving NSD, which contradicts experimental evidence showing that NSD evolves as the gel approaches a glassy state.

2. Methodology

The researchers employed large-scale Langevin dynamics simulations to test their theories across a vast parameter space, covering two distinct types of gels:

• Depletion Gels: Modeled using a Morse potential to simulate short-ranged entropic attractions (mimicking non-adsorbing polymer depletants).
• Frictional Gels: Modeled using the Discrete Element Method (DEM), incorporating non-hydrodynamic contact friction (sliding and rolling resistance) to simulate solid-solid interparticle bonds.

Simulation Setup:

• System: 20,000 bi-disperse spheres in a 3D box.
• Deformation: Quasistatic uniaxial compression (slow axial movement of walls) to reach various volume fractions ($\phi = 0.25–0.5$).
• Analysis: The authors used spherical harmonic expansions to represent the probability density functions of contact orientations and the angular modulation of interaction forces (normal and tangential).

3. Key Contributions

The paper proposes a new framework for constitutive relations that shifts the focus from macroscopic strain to microscopic internal state parameters.

• New Constitutive Relations: Instead of expanding free energy around an unstressed state, the authors relate the global stress tensor ($\sigma_{ij}$) to the volume averages of the dyadic products of pairwise interaction forces and branch vectors.
• Anisotropy Parameters: They introduce a method to quantify stress through three specific anisotropy parameters:
  • $a$: Contact orientation anisotropy.
  • $a_n$: Normal force anisotropy.
  • $a_t$: Tangential force anisotropy.
• Symmetry-Based Modeling: By using scalar and vector spherical harmonics, they developed equations (Eqs. 9–11) that correctly account for the azimuthal symmetry and $\pi$-periodicity of the force networks.

4. Results

• Failure of Old Models: The authors demonstrated that the traditional Cauchy-Born expansion predicts a constant $N_1 = N_2 = 2/3$, which is qualitatively incorrect compared to simulation and experimental data.
• Origin of Stress Anisotropy: A major finding is that while the contact network remains almost isotropic ($a \approx 0$) due to strong attractions, the force network is highly anisotropic. The stress anisotropy is driven by the "force skeletons" (anisotropic normal and tangential forces) rather than the arrangement of the particles themselves.
• Predictive Accuracy: The proposed relations (Eqs. 9–11) showed "remarkable" agreement with simulation data for both axial stress ($\sigma_{zz}$), transverse stress ($\sigma_{xx}$), and the evolution of the Normal Stress Difference (NSD) across varying volume fractions and attraction strengths.
• Universality: The angular modulation of forces follows a universal form regardless of whether the interaction is central (depletion) or non-central (frictional).

5. Significance

This work represents a paradigm shift in the rheology of soft matter. It moves the theoretical description of gels away from cluster-length-scale mechanisms (which the authors argue are inefficient at capturing elastic response) toward particle-length-scale phenomena.

By successfully encoding the non-linear effects of quenched stresses and force-chain anisotropy, this model provides a robust tool for predicting the mechanical behavior of complex colloidal systems during industrial processes like consolidation, drying, and fabrication.