How hydrodynamic interactions alter polymer stretching in turbulence

Using Brownian dynamics simulations in homogeneous isotropic turbulence, this study reveals that hydrodynamic interactions significantly alter polymer stretching by inducing a steeper coil-stretch transition, limiting the range of power-law extensions, and delaying state migration, thereby highlighting the necessity of incorporating extension-dependent drag forces in coarse-grained models of turbulent polymer solutions.

The Gist

Imagine a long, floppy noodle floating in a swirling, chaotic river. This noodle represents a polymer molecule, and the river represents a turbulent fluid. Scientists have long known that if you pull this noodle through calm water, the water itself pushes back on different parts of the noodle in a way that changes how it stretches. This is called Hydrodynamic Interaction (HI).

However, when the river is a raging storm (turbulence), nobody was sure if this "water pushing back" still mattered. This paper uses computer simulations to figure out exactly how these interactions change the noodle's behavior in the storm.

Here is the breakdown of their findings using simple analogies:

1. The "Two-Beetle" vs. The "Long Train"

To study this, the researchers modeled the polymer in two ways:

• The Dumbbell (Two Beads): Imagine the polymer as just two heavy beads connected by a single spring. It's like a dumbbell.
• The Chain (Many Beads): Imagine the polymer as a long train of many beads connected by springs.

The Big Surprise:
When they added the "water pushing back" (HI) to the Dumbbell, it barely changed anything. The two beads are so far apart that they don't really hide each other from the water's flow.

• Analogy: It's like two people standing far apart in the rain; neither one shields the other from getting wet.

But when they added the same "water pushing back" to the Long Chain, the results changed dramatically.

• Analogy: Now imagine a long line of people holding hands. The people in the middle are shielded from the rain by the people on the outside. The whole group gets wet much slower than if they were just two people standing apart.

The Lesson: You cannot understand how a long, complex polymer behaves in a storm just by looking at a simple two-bead model. The "shielding" effect only happens when you have enough beads to actually coil up.

2. The "Coil-Stretch" Dance

In a turbulent flow, these polymers are constantly being stretched out by the current and then snapping back into a ball (coiling) when the current relaxes.

• Without HI: The polymer stretches and snaps back relatively easily.
• With HI (The Long Chain): The "shielding" effect acts like a heavy anchor.
  • When the chain is coiled up (like a ball of yarn), the outer beads shield the inner ones, making the whole ball feel "heavier" and harder to pull apart. It stays coiled longer.
  • When the chain is stretched out, the beads are far apart, the shielding disappears, and the water drags on them more easily.

The Result: The transition between being a tight ball and a stretched string becomes much sharper. The polymer gets "stuck" in one state or the other for longer periods. It's like a door that is hard to open but hard to close; once it's open, it stays open, and once it's closed, it stays closed.

3. The "Traffic Jam" of Shapes

The researchers looked at how often the polymer is in a "coiled" state versus a "stretched" state.

• Without HI: The polymer spends a decent amount of time in the middle ground—somewhat stretched, somewhat coiled.
• With HI: The polymer avoids the middle ground. It is either very tightly coiled or very fully stretched. The "middle" range disappears.

The Analogy: Imagine a traffic light that usually cycles through Red, Yellow, and Green. With HI, the light seems to skip the Yellow phase entirely, flipping instantly between Red and Green. The polymer spends almost no time in the "in-between" state.

4. Why the "Dumbbell" Model Fails

Many computer simulations of turbulent fluids use the simple "dumbbell" model because it's easy to calculate. This paper argues that this is a mistake if you want to be accurate.

• Because a dumbbell can't actually coil up (it's just two beads), it can't experience the "shielding" effect.
• Therefore, adding HI to a dumbbell model doesn't fix the problem; it just gives you the wrong answer. To see the real physics, you need a model with enough "beads" to actually form a coil.

5. A Simpler Way to Simulate

Finally, the researchers tested if they could replace the complex, real-world turbulent river with a simpler, made-up "random flow" (a mathematical model that looks like turbulence but is easier to generate).

• The Finding: Surprisingly, the simple random model worked just as well as the complex real turbulence for predicting how these polymers stretch.
• Why it matters: This means scientists can use this simpler, faster computer model to test new theories about polymers without needing to run massive, expensive simulations of real turbulence.

Summary

In short, this paper tells us that complexity matters. If you want to know how a long polymer behaves in a storm, you can't just look at a simple two-part model. You need to account for how the different parts of the chain hide each other from the water. This "hiding" makes the polymer act more stubborn, staying coiled or stretched for longer times, and skipping the middle ground entirely.

Technical Summary

1. Problem Statement

The study addresses a critical gap in the understanding of polymer dynamics in turbulent flows: the role of Hydrodynamic Interactions (HI).

• Context: In laminar, extension-dominated flows, it is well-established that HI (solvent-mediated interactions between polymer segments) significantly alter polymer stretching, leading to phenomena like coil-stretch hysteresis and conformation-dependent drag.
• The Gap: In turbulent flows, where strain rates fluctuate rapidly, the role of HI has remained unclear. Most existing studies on Lagrangian polymer dynamics in turbulence ignore HI, often relying on "free-draining" models (like the Rouse model) or simplified dumbbell models (two beads connected by a spring).
• The Core Question: Does adding HI to a simple dumbbell model accurately capture the physics of HI in a more realistic, multi-bead polymer chain within a turbulent environment? Or does the coarse-graining level (number of beads) fundamentally change the outcome?

2. Methodology

The authors employed Brownian dynamics (BD) simulations to track the stretching of bead-spring chains in a fluctuating velocity field.

• Polymer Model:

  • Bead-Spring Chains: Modeled as freely jointed chains of $N_b$ beads connected by FENE (Finitely Extensible Nonlinear Elastic) springs.
  • Hydrodynamic Interactions: Implemented using the Rotne-Prager-Yamakawa (RPY) mobility tensor, which accounts for the shielding of drag forces between beads as the chain coils and stretches.
  • Coarse-Graining Hierarchy: Simulations were performed for chains with varying numbers of beads: $N_b = 2$ (dumbbell), $4$, $10$, and $20$.
  • Parameter Mapping: To ensure chains of different $N_b$ represented the same physical polymer, the authors used the Jin-Collins mapping. This ensures that the ratio of contour length to equilibrium radius of gyration ($R_m/R_{eq}$) remains constant across different levels of coarse-graining.
  • HI Strength: The strength of HI was controlled by the parameter $h^*$ (related to bead radius), tested at values $0$ (free-draining), $0.2$ (intermediate), and $0.49$ (maximum physical limit).

• Flow Field:

  • Turbulent Flow: The chains were transported in a homogeneous isotropic turbulence field generated by Direct Numerical Simulation (DNS) with a Taylor-microscale Reynolds number $Re_\lambda \approx 111$.
  • Random Flow Model: To test the robustness of the findings, a time-correlated Gaussian random flow (Brunk-Koch-Lion model) was also used. This model mimics the Lyapunov exponent and temporal correlations of the turbulent flow but lacks non-Gaussian extreme events.
  • Regime: The study focused on the ultra-dilute, one-way coupling regime, meaning the polymers are passive tracers and do not feedback to alter the flow (though the paper discusses the implications of this for future work).

• Key Metrics:

  • Mean extension $\langle R \rangle$.
  • Probability Distribution Functions (PDF) of the end-to-end extension $R$.
  • Persistence times in coiled vs. stretched states.
  • Weissenberg number ($Wi = \lambda \tau^D$), where $\lambda$ is the Lyapunov exponent and $\tau^D$ is the relaxation time.

3. Key Contributions

1. Validation of Coarse-Graining Sensitivity: The paper demonstrates that the effect of HI is not captured by simply adding HI to a dumbbell model. The qualitative and quantitative behavior of HI depends heavily on the number of beads ($N_b$).
2. Mechanism of Hydrodynamic Shielding: It elucidates that the suppression of stretching in multi-bead chains is due to hydrodynamic shielding in coiled states, a physical phenomenon a dumbbell ($N_b=2$) cannot replicate because it cannot form a physical coil.
3. Random Flow Surrogate: It validates that a time-correlated Gaussian random flow can accurately reproduce the stretching statistics of polymers with HI in turbulence, offering a computationally cheaper alternative to DNS for developing reduced-order models.

4. Key Results

A. Dumbbells vs. Multi-bead Chains (The "Dumbbell Failure")

• Dumbbells ($N_b=2$): HI has a negligible effect on mean stretching. In fact, HI slightly increases the extension of dumbbells (consistent with analytical predictions). This is because a dumbbell cannot coil; thus, it does not experience the hydrodynamic shielding that reduces drag in coiled configurations.
• Multi-bead Chains ($N_b \ge 4$): HI has a strong, suppressing effect on stretching. As $N_b$ increases, the shielding effect becomes more pronounced.
  • Stiff Polymers (Low $Wi$): HI causes chains to remain coiled longer, significantly reducing the mean extension compared to free-draining chains.
  • Elastic Polymers (High $Wi$): HI causes chains to stretch more than free-draining chains once the transition occurs, due to the delayed migration between states.
• Conclusion: A dumbbell with HI is not a valid surrogate for a realistic chain with HI. The effects of HI are qualitatively different (opposite direction of change in mean extension) and quantitatively distinct.

B. Steepening of the Coil-Stretch Transition

• For multi-bead chains, the presence of HI steepens the coil-stretch transition.
• The transition becomes sharper: chains stay coiled longer at lower $Wi$ and stretch more abruptly at higher $Wi$.
• This is attributed to an effective conformation-dependent drag: coiled chains experience high drag (due to shielding), while stretched chains experience lower drag (deshielding). This slows the migration between states, creating a "hysteresis-like" effect even in fluctuating flows where true hysteresis is usually washed out.

C. Modification of Extension Statistics (PDF)

• Power-Law Disruption: In free-draining chains, the PDF of extension $P(R)$ exhibits a power-law behavior ($R^{-\alpha}$) over a wide range. HI significantly limits the range of this power-law behavior.
• Bimodality: For chains with HI, the PDF develops a bimodal shape near the transition ($Wi \approx 0.5$), with distinct peaks near the coiled state ($R \approx R_{eq}$) and the stretched state ($R \approx R_m$).
• Skewness and Variance: The variance and skewness of the extension distribution increase significantly with $h^*$ for multi-bead chains, whereas they decrease slightly for dumbbells.

D. Persistence Times

• HI delays the migration between coiled and stretched states.
• The persistence time distributions for both coiled and stretched states show exponential tails that broaden significantly with increasing $h^*$. This confirms that HI acts as a "memory" mechanism, trapping chains in their current conformational state for longer durations.

E. Random Flow vs. DNS

• The Gaussian random flow model successfully reproduced the stretching statistics (mean extension, PDF shape, persistence times) of the DNS results for chains with HI.
• This implies that extreme strain-rate events (non-Gaussian tails in the velocity gradient distribution) are not the primary driver of HI effects; rather, the mean Lyapunov exponent and temporal correlations are sufficient to capture the physics.

5. Significance and Implications

• For Turbulent Drag Reduction (TDR): Since TDR is driven by highly stretched polymers, the finding that HI suppresses stretching in the transition regime (and alters the PDF) suggests that current continuum models (like Oldroyd-B or FENE-P) based on dumbbell assumptions may underestimate or misrepresent the efficiency of drag reduction for real polymers.
• Model Development: The study argues that to accurately simulate turbulent polymer solutions, one must move beyond simple dumbbells.
  • Recommendation: Develop reduced-order models (e.g., modified dumbbells) that incorporate extension-dependent drag or conformation-aware friction to mimic the shielding effects of multi-bead chains.
• Computational Efficiency: The validation of the Gaussian random flow as a surrogate allows researchers to test these new coarse-grained models without the immense computational cost of running full DNS, facilitating the development of data-driven models for viscoelastic turbulence.
• Future Directions: The authors highlight the need for "successive fine-graining" studies to determine the asymptotic behavior of HI as $N_b$ approaches the Kuhn step resolution ($N_k$), and to investigate these effects in anisotropic flows (e.g., pipe flow) where folding artifacts might differ.

In summary, the paper establishes that hydrodynamic interactions are a critical, coarse-graining-dependent phenomenon in turbulent polymer stretching. Ignoring the physical ability of chains to coil (as dumbbells do) leads to qualitatively incorrect predictions, necessitating more sophisticated modeling approaches for accurate prediction of viscoelastic turbulence.