Beyond binary scission: a generalized three-species cascade breakage model for wormlike micellar solutions

Technical Summary: A Generalized Three-Species Cascade Breakage Model for Wormlike Micellar Solutions

Problem Statement
Wormlike micellar (WLM) solutions exhibit complex rheological behaviors, including viscoelasticity, shear thinning, and shear banding, driven by the continuous breakage and recombination of self-assembled micellar networks. Existing theoretical frameworks, particularly the widely used two-species models (e.g., the Vasquez-Cook-McKinley or VCM model), rely on a coarse binary representation of the micellar population (typically distinguishing only between long entangled chains and short fragments). While successful in capturing certain macroscopic features, these models possess a fundamental limitation: they fail to explicitly resolve the multiscale cascade breakage and the resulting continuous relaxation regimes inherent to real micellar networks. This binary approximation filters out intermediate structural states, leading to an inability to accurately describe broad relaxation spectra, particularly the high-frequency upturn in dynamic moduli observed in experiments, and limits the physical interpretation of transient responses and shear banding mechanisms.

Methodology
To address these limitations, the authors develop a three-species cascade breakage model that introduces an intermediate micellar state into the kinetic framework. The model assumes the micellar network consists of three discrete structural states:

1. Species A: A gel-network (primary, highly entangled).
2. Species B: Long chains (intermediate state).
3. Species C: Short chains (disentangled fragments).

The mathematical formulation is based on a kinetic theory approach using continuous phase-space distributions. The authors derive coupled Smoluchowski equations for the probability density functions of the three species, incorporating convective stretching, diffusion, and reaction terms. Key methodological features include:

• Hierarchical Kinetics: The model defines a cascade pathway where the gel-network (A) breaks into long chains (B), which subsequently break into short chains (C). Recombination pathways are also defined to ensure mass conservation.
• Flow-Dependent Scission: The breakage rates for the gel-network and long chains are coupled to the local macroscopic flow via a partially extending strand convection (PEC) mechanism. This creates a stress-dependent breakage mechanism where deformation accelerates scission.
• Constitutive Equations: The model derives dimensionless evolution equations for number densities and conformation tensors, coupled to the macroscopic Cauchy momentum equation. The total stress is the superposition of elastic stresses from all three species.
• Numerical Implementation: The authors analyze the model under both homogeneous flow conditions (to study linear viscoelasticity and steady shear response) and inhomogeneous cylindrical Couette flow (to investigate shear banding, transient startup, and relaxation).

Key Contributions and Results

1. Three-Mode Viscoelastic Response:
   Under homogeneous flow, the model analytically yields a discrete three-mode Maxwell structure. Unlike two-species models, this framework successfully captures the high-frequency upturn in the dynamic storage ($G'$) and loss ($G''$) moduli. The authors demonstrate that the additional relaxation mode provided by the intermediate species (B) and the short chains (C) significantly improves the fit to experimental master curves for CPyCl/NaSal solutions, particularly in the high-frequency regime.

2. Non-Monotonic Flow Curves and Shear Banding:
   The model predicts a non-monotonic constitutive curve with a region of strictly negative slope, a hallmark of constitutive instability. In inhomogeneous Couette flow simulations, this instability leads to the spontaneous formation of coexisting shear bands (high and low shear rate regions). The three-species model reveals that the intermediate species (B) forms a distinct transition layer at the band interface, linking the gel-network-dominated low shear band and the short-chain-dominated high shear band. This provides a more resolved microstructural description of the band interface than binary models allow.

3. Transient Dynamics and Hysteresis:
   The model captures complex transient behaviors:

   • Stress Overshoot: During startup, the competition between elastic loading and accelerated cascade breakage leads to pronounced stress overshoots.
   • Flow Hysteresis: The model exhibits hysteresis, where the final steady state depends on the flow history (startup from rest vs. step-down from high shear).
   • Multistep Relaxation: Upon shear cessation, the model predicts a multistep relaxation process. Instead of a single exponential decay, the stress relaxes through sequential stages: rapid release from short chains (C), intermediate adjustment by long chains (B), and slow recovery of the gel-network (A). This resolves the "prolonged recovery timescale" often observed in experiments but difficult to capture with reduced models.

4. Continuum Mechanical Consistency:
   The model satisfies the Lodge-Meissner relation ($N_1/\sigma_{r\theta} = \gamma_0$) during large amplitude step strains, confirming its tensorial consistency and physical validity under finite deformation.

Significance and Claims
The paper claims that the proposed three-species framework bridges the gap between discrete phenomenological network models and continuous kinetic theories (such as Cates' original reptation-reaction theory). By explicitly resolving the cascade kinetics governing intermediate structural states, the model provides a physically interpretable framework for wormlike micellar shear banding.

The authors assert that this approach captures the intrinsic connection between:

• Cascade microstructural evolution (A $\to$ B $\to$ C).
• Broad relaxation dynamics (three distinct time scales).
• Macroscopic flow localization (shear banding).

The significance lies in the model's ability to maintain macroscopic computational tractability while recovering missing intermediate relaxation modes. This offers a more comprehensive physical basis for interpreting the rich transient and steady-state dynamics of complex micellar fluids, potentially aiding in the interpretation of future microfluidic experiments and industrial applications involving WLMs. The paper does not propose new experimental protocols but rather offers a theoretical tool to better interpret existing rheological phenomena.