🔬 condensed matter

Numerical analysis of capillarity-driven thinning rheometry for polydisperse polymer solutions

This study demonstrates that the characteristic thinning timescale ( $\tau_{EC}$ ) measured in capillarity-driven rheometry of polydisperse polymer solutions does not represent the full molecular weight distribution but rather reflects a specific sub-ensemble of high-molecular-weight chains actively stretched by the flow, making $\tau_{EC}$ an experiment-specific quantity dependent on concentration and setup parameters rather than an intrinsic fluid property.

Original authors: Isaac Pincus, Vincenzo Calabrese, Simon J. Haward, Gareth H. McKinley

Published 2026-07-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Isaac Pincus, Vincenzo Calabrese, Simon J. Haward, Gareth H. McKinley

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ 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 drop of sticky, stretchy liquid—like a super-charged honey mixed with tiny, invisible rubber bands—suspended between two plates. When you pull the plates apart, the liquid forms a thin thread that tries to snap back into a ball. This is the "CaBER" experiment, a way scientists study how fluids stretch and break.

For a long time, scientists thought this stretching thread behaved like a simple, single-speed rubber band. They believed that by measuring how fast the thread thinned down, they could find one single "relaxation time" (let's call it the snap-back speed) that described the entire fluid. It was like assuming a whole orchestra played at the exact same tempo.

But this new paper, by Isaac Pincus and his team, pulls the curtain back on that idea. They show that for real-world polymer solutions (which are messy mixes of rubber bands of all different lengths), that single "snap-back speed" is a bit of a trick.

The Great Rubber Band Filter

The researchers used a computer model to simulate what happens inside that stretching thread. They treated the fluid not as one big blob, but as a crowd of individual polymer chains, some short and some very long.

Here is the magic trick they discovered: The flow acts like a filter.

When the thread starts stretching, the fluid moves faster and faster. The short, snappy rubber bands (low molecular weight) are like little kids trying to run a marathon; they get tired and stop stretching almost immediately. They relax back to their curly shape and stop contributing to the tension.

The long, lazy rubber bands (high molecular weight), however, are like marathon runners. Because they are so long, they take a long time to relax. When the fluid stretches, these long chains stay stretched out for much longer. They are the only ones strong enough to hold the thread together against the surface tension trying to snap it.

The paper finds that only the chains that are stretched enough to be "awake" (specifically, those with a "Weissenberg number" greater than 0.5) actually do the work of holding the thread. The short chains just sit there, relaxed and useless.

The "Fake" Relaxation Time

Because only the long chains are doing the heavy lifting, the speed at which the thread thins down (the elastocapillary timescale, or $\tau_{EC}$ ) doesn't tell you about the average speed of all the chains. Instead, it tells you about the speed of just the longest, strongest chains that are still stretched.

The paper explicitly argues against the idea that $\tau_{EC}$ is a fixed, intrinsic property of the fluid (like a fingerprint that never changes). Instead, their simulations show that $\tau_{EC}$ is a "chameleon" that changes based on:

The mix: If you add a tiny bit of long chains to a sea of short chains, the long chains take over completely. The short chains become invisible to the measurement.
The concentration: How much polymer is in the water changes which chains get stretched.
The experiment setup: Even how fast you pull the plates apart (pre-stretch) or how wide the thread starts (initial diameter) changes the result.

In their simulations, they found that for a mixture of two polystyrene samples (labeled PS7 and PS16), when the concentration of the high-molecular-weight PS16 reached just 50 ppm, there was very little difference in the thinning speed whether the solution contained 1000 ppm of the short PS7 chains or none at all. However, once the PS16 concentration went up, it completely dominated the behavior, making the thread stretch much longer. The short chains (PS7) were effectively ignored.

The "Stress-Weighted" Truth

The authors ran detailed computer simulations (using a model called FENE-PM) to prove this. They didn't just guess; they calculated the stress contribution of every single "sub-species" of chain in the mix.

They found that the "average" molecular weight you might calculate on paper (like a simple average of all the chains) is wrong. The real "average" that matters is a stress-weighted average. This means you only count the chains that are actually pulling their weight.

For example, in their mix of 1000 ppm PS7 and 500 ppm PS16, the short chains (PS7) were actually more numerous, but the long chains (PS16) were the ones holding the stress. The measured time scale ( $\tau_{EC}$ ) shifted to match the long chains, ignoring the short ones entirely.

What This Means (and What It Doesn't)

The paper confirms that the "snap-back speed" we measure in the lab is not the relaxation time of the whole fluid. It is a specific, experiment-dependent number that reflects only the subset of chains that are stretched enough to matter at that exact moment.

They also ruled out the idea that this happens because the chains are bumping into each other (interacting) in a simple way. Instead, it's about the flow itself deciding which chains get to play. The short chains simply don't get stretched enough to stay in the game.

The authors are careful to say this is based on their simulations and validated against specific experiments with polystyrene in a solvent called dioctyl phthalate (DOP). They suggest that for other fluids, the same rules likely apply, but they haven't tested every possible liquid.

So, the next time you see a sticky thread stretching, remember: it's not a team effort. It's a solo act by the longest, strongest chains, while the rest of the crowd has already gone home. The "speed" you measure is really just the speed of the stars of the show, not the whole cast.

Technical Summary: Numerical analysis of capillarity-driven thinning rheometry for polydisperse polymer solutions

Problem Statement
Capillary Breakup Extensional Rheometry (CaBER) is a standard technique for characterizing the extensional viscosity of polymer solutions. In idealized monodisperse systems described by the Oldroyd-B model, the filament thinning rate in the elastocapillary (EC) regime is inversely proportional to the fluid's relaxation time ( $\tau$ ), allowing the measured decay rate ( $\dot{\epsilon}_{EC}$ ) to define a characteristic relaxation time $\tau_{EC} = 2/(3\dot{\epsilon}_{EC})$ . However, for real, polydisperse polymer solutions, the interpretation of $\tau_{EC}$ is ambiguous. Previous studies have noted discrepancies between $\tau_{EC}$ and relaxation times measured via shear or linear viscoelasticity. Specifically, recent experiments on bidisperse polystyrene blends (Calabrese et al., 2025) revealed that adding low-molecular-weight (MW) species to high-MW solutions at high concentrations had negligible impact on the measured $\tau_{EC}$ , suggesting that the thinning dynamics are governed by a specific subset of the molecular weight distribution rather than the ensemble average. A comprehensive framework linking the molecular weight distribution, flow history, and the resulting $\tau_{EC}$ was lacking.

Methodology
The authors developed a numerical model to simulate the thinning of a cylindrical polymer filament, explicitly incorporating polydispersity. The approach combines:

Force Balance: A one-dimensional force balance equation for the filament radius $a(t)$ , accounting for inertial, viscous, capillary, and polymeric stresses.
Constitutive Modeling: A multi-mode FENE-PM (Finitely Extensible Nonlinear Elastic - Peterlin Mean) model. The polymer solution is treated as a mixture of distinct "species" (representing different MW distributions, specifically PS7 and PS16 polystyrene samples).
- Linear Regime: The model utilizes the Zimm model with pre-averaged hydrodynamic interactions to describe linear viscoelasticity, generating a spectrum of relaxation times $\tau_{is}$ for each species.
- Nonlinear Regime: The model incorporates finite extensibility via FENE springs and a conformation-dependent drag term. This drag term modifies the relaxation time of the longest mode ( $\hat{\tau}_{i1}$ ) as a function of chain extension, accounting for the change in hydrodynamic friction as chains stretch from coils to ellipsoids.
Simulation: The system of ordinary differential equations governing the filament radius and the conformation tensors of each species/mode is solved simultaneously. The model is validated against experimental data from Calabrese et al. (2025) involving bidisperse blends of PS7 and PS16 in dioctyl phthalate (DOP).

Key Contributions and Results

Selective Stretching Mechanism: The simulations demonstrate that during the EC regime, the elongational strain rate $\dot{\epsilon}$ decreases over time. Consequently, only polymer chains with an effective Weissenberg number $Wi = \dot{\epsilon}\tau > 0.5$ (the coil-stretch transition threshold) remain stretched and contribute significantly to the elastic stress. Chains with shorter relaxation times (lower MW) fall below this threshold, relax, and cease to support the filament.
Dominance of High-MW Species: Because relaxation time scales with molecular weight ( $\tau \sim M^{1.5}$ ), high-MW species maintain $Wi > 0.5$ longer than low-MW species. The measured $\tau_{EC}$ is therefore determined by the "stress-contributing sub-ensemble" of chains, which is heavily biased toward the high-MW tail of the distribution. This explains why adding low-MW polymer to a high-MW solution does not alter $\tau_{EC}$ once the high-MW concentration is sufficient to carry the total stress.
Conformation-Dependent Drag: The inclusion of conformation-dependent drag was necessary to match experimental data at higher concentrations. Without this term, the model predicted $\tau_{EC}$ values that were too low compared to experiments, failing to capture the increase in relaxation time observed as chains stretch and interact.
Experimental Parameter Dependence: The study confirms that $\tau_{EC}$ is not an intrinsic fluid property but depends on experimental conditions. Specifically, the initial filament diameter ( $a_0$ ) and the pre-stretch history (plate separation rate) influence the time available for short modes to relax before the EC regime is established. Larger initial diameters and higher pre-stretch rates lead to longer measured $\tau_{EC}$ values.
Re-evaluation of Averages: The paper contrasts the "extensibility-weighted molecular weight" ( $M_L$ ) with a "stress-weighted" molecular weight. While $M_L$ captures the potential of high-MW chains to contribute stress, the stress-weighted average (derived from the simulation at the minimum strain rate) more accurately reflects the actual subset of chains active during the EC regime.

Significance and Claims
The paper argues that the characteristic thinning timescale $\tau_{EC}$ extracted from CaBER experiments on polydisperse fluids should not be interpreted as a single, intrinsic relaxation time of the fluid. Instead, it is an experiment-specific quantity reflecting the relaxation time of the specific subset of chains that are actively stretched by the flow at the moment the elastocapillary balance is achieved.

The authors conclude that for polydisperse solutions, $\tau_{EC}$ is a function of:

The molecular weight distribution (specifically the high-MW tail).
The total polymer concentration.
Experimental parameters such as initial filament diameter and pre-stretch history.

This finding resolves the discrepancy between shear/linear viscoelastic relaxation times and extensional thinning times in complex fluids. It suggests that the "spinnability" and stability of polymer filaments are disproportionately influenced by the high-MW fraction, a critical insight for understanding phenomena ranging from industrial fiber spinning to droplet breakup in biological systems. The work provides a computational framework for analyzing CaBER data on polydisperse fluids, moving beyond the limitations of single-mode Oldroyd-B interpretations.

The Great Rubber Band Filter

The "Fake" Relaxation Time

The "Stress-Weighted" Truth

What This Means (and What It Doesn't)

More like this