Spatial instability analysis and mode transition of a viscoelastic jet in a co-flowing gas stream

This study employs spatial linear instability analysis and energy budget methods to demonstrate that increasing Weber numbers and elasticity in a viscoelastic jet within a co-flowing gas stream drive a transition from axisymmetric to helical modes, revealing a distinct elasticity-enhanced shear-driven instability mechanism validated by experimental flow-focusing results.

The Gist

Imagine you are squeezing a tube of toothpaste. If the toothpaste is just water, it shoots out in a straight line and breaks into neat, round drops. But if you mix in some slime (polymer), the toothpaste becomes stretchy and "springy." It doesn't just break; it might wiggle, twist, or spiral like a corkscrew before it snaps.

This paper is a deep dive into why that happens. The researchers are studying a "viscoelastic jet"—a stream of stretchy liquid surrounded by a fast-moving stream of gas (like air blowing on the toothpaste). They want to know: When does the stream stay straight, and when does it start twisting into a spiral?

Here is the breakdown of their discovery, using simple analogies:

1. The Setup: The "Hair Dryer" Experiment

Think of the experiment as a high-tech version of a hair dryer blowing on a stream of honey.

• The Liquid: A stretchy fluid (like a polymer solution).
• The Gas: A fast stream of air blowing alongside it.
• The Goal: To see how the liquid stream breaks apart. Does it stay a straight, round rope (axisymmetric), or does it start wobbling and twisting like a snake (helical)?

2. The Two Main "Dancers" (Modes of Instability)

The researchers found that the liquid stream has two main ways of dancing before it breaks:

• The "Pulsing" Dance (Axisymmetric): The stream gets fat and thin in a rhythmic pattern, like a sausage. It stays perfectly round but changes thickness. This usually happens when the air isn't blowing too hard.
• The "Twisting" Dance (Helical): The stream starts to spiral. The center of the stream wiggles side-to-side, creating a corkscrew shape. This happens when the air blows harder or when the liquid is very stretchy.

3. The Secret Ingredient: "Elasticity"

The big surprise in this paper is the role of elasticity (the "springiness" of the liquid).

• Low Elasticity (Weak Spring): If the liquid is only slightly stretchy, the air pressure is the main boss. The air pushes the liquid, and the liquid breaks based on how hard the air blows.
• High Elasticity (Strong Spring): As the liquid gets more stretchy, something new happens. The liquid's own internal "memory" and tension take over. The researchers discovered a new mechanism they call "Elasticity-Enhanced Shear-Driven Instability."

The Analogy: Imagine a group of people running side-by-side.

• In a normal crowd (Newtonian fluid), if someone pushes from the side, the whole group wobbles together.
• In a stretchy crowd (viscoelastic fluid), the people are holding rubber bands between them. If the person on the outside runs faster, the rubber bands stretch. The energy from that stretching doesn't just sit there; it gets transferred into the people in the middle, making them twist and turn violently. The "spring" inside the liquid actually helps the twisting happen faster.

4. The Map: When to Expect a Twist

The authors created a "weather map" (a Phase Diagram) for these jets.

• X-axis: How hard the air is blowing (Weber number).
• Y-axis: How stretchy the liquid is (Elasticity number).

The Rules of the Map:

• Low Air + Low Stretch: The stream stays straight and round (Pulsing Dance).
• High Air: The stream starts twisting (Twisting Dance), even if the liquid isn't very stretchy.
• High Stretch: If the liquid is very stretchy, it starts twisting much earlier, even with just a little bit of air. The elasticity acts like a catalyst, speeding up the transition to the spiral shape.

5. Why Their Method Matters

Previous scientists tried to predict this by watching a single spot in time (like taking a photo of a runner). But jets are moving fast; the instability grows as it travels downstream.

• The Old Way: Like watching a runner stand still and guessing how fast they will run later. It often got the answer wrong.
• The New Way (This Paper): They used a "Spatial" analysis, which is like watching the runner move down the track. They tracked how the wobble grows as the liquid moves away from the nozzle. This method matched their real-world experiments perfectly, proving that you have to watch the whole journey, not just the starting line, to understand how the jet breaks.

The Bottom Line

This paper explains that stretchy liquids don't just break; they dance.
When you blow air on a stretchy liquid, the liquid's internal "springs" can grab the energy from the air and the speed difference, turning a simple wobble into a dramatic spiral. This knowledge helps engineers design better inkjet printers, create perfect micro-fibers for clothing, and make consistent medicine droplets, ensuring they don't accidentally twist and break in the wrong place.

Technical Summary

1. Problem Statement

The study investigates the linear instability and mode transition mechanisms of a viscoelastic liquid jet surrounded by a co-flowing gas stream, a configuration typical of flow-focusing devices used for generating monodisperse droplets and microfibers.

• Context: While Newtonian jet instability is well-understood (transitioning from capillary to Kelvin-Helmholtz modes), the behavior of viscoelastic jets (e.g., polymer solutions) is complex due to fluid elasticity.
• Gap: Previous studies on viscoelastic jets have largely relied on temporal instability analysis, which assumes disturbances grow in time at a fixed position. This approach is insufficient for open flows like jets, where instabilities are convective (growing as they travel downstream). Furthermore, the specific mechanisms driving the transition between axisymmetric ($n=0$) and helical ($n=1$) modes in viscoelastic co-flows, and the role of elasticity in this transition, remain unclear.

2. Methodology

The authors employ a comprehensive theoretical framework validated by controlled experiments:

• Theoretical Model:

  • Governing Equations: Incompressible Navier-Stokes equations for the liquid and gas phases.
  • Constitutive Model: The Oldroyd-B model is used to describe the viscoelastic liquid, capturing both solvent viscosity and polymeric stress.
  • Base Flow: Unlike previous studies using uniform profiles, this work adopts a non-uniform hyperbolic-tangent axial velocity profile to accurately represent the shear layer at the liquid-gas interface.
  • Analysis Type: Spatial linear instability analysis is used. Here, the frequency ($\omega$) is real, and the wavenumber ($k$) is complex. This directly resolves the downstream amplification of disturbances, making it suitable for convective flows.
  • Numerical Solution: The eigenvalue problem is solved using a Chebyshev spectral collocation method. The quadratic eigenvalue problem is linearized via an auxiliary variable transformation.
  • Mechanism Analysis: An energy budget analysis is performed to decompose the growth rate of disturbance kinetic energy into physical contributions (e.g., surface tension, gas pressure, shear stress, elastic stress).

• Experimental Validation:

  • Experiments were conducted using a flow-focusing device with polyethylene oxide (PEO) solutions in glycerol-water mixtures.
  • Parameters (Reynolds number $Re$, Weber number $We$, and Elasticity number $El$) were systematically varied by adjusting polymer molecular weight, concentration, and gas pressure.
  • Jet morphologies (axisymmetric vs. helical) were observed and compared against theoretical predictions.

3. Key Contributions

1. Spatial vs. Temporal Framework: The study demonstrates the superiority of spatial instability analysis over temporal analysis for viscoelastic jets. It shows that temporal analysis fails to predict the correct transition boundaries between modes, often overestimating the stability of the axisymmetric mode.
2. Discovery of a New Instability Mechanism: The authors identify a distinct instability mechanism termed Elasticity-Enhanced Shear-Driven Instability (ESI). Unlike traditional capillary or Kelvin-Helmholtz instabilities, ESI arises from enhanced energy exchange between the base-flow shear and velocity perturbations within the jet interior, driven by elasticity.
3. Mode Transition Map: A comprehensive $We$–$El$ phase diagram is constructed, delineating the boundaries between predominant axisymmetric and helical modes, which is validated against experimental data.

4. Key Results

A. Influence of Elasticity ($El$)

• Low Elasticity: The instability is primarily governed by gas-pressure fluctuations and interfacial shear (Kelvin-Helmholtz mechanism). The axisymmetric mode dominates.
• High Elasticity: As $El$ increases, the predominant mode transitions from axisymmetric ($n=0$) to helical ($n=1$).
• Mechanism Shift: In the high-elasticity regime, the instability mechanism shifts to ESI. The disturbance structures migrate from the interface toward the jet interior (core). The elastic stresses enhance the coupling between velocity perturbations and the basic flow shear, allowing the helical mode to extract energy more efficiently from the liquid core than the axisymmetric mode.
• Growth Rate: Elasticity initially enhances instability at low $El$ but suppresses it at very high $El$ due to nonlinear viscoelastic effects.

B. Influence of Weber Number ($We$)

• Low $We$: Surface tension dominates, leading to capillary instability (axisymmetric/varicose mode).
• High $We$: Inertia dominates. As $We$ increases, the system transitions to a Kelvin-Helmholtz instability driven by gas-pressure fluctuations and interfacial shear.
• Mode Transition: Increasing $We$ drives a transition from axisymmetric to helical modes. At high $We$, the phase speed of disturbances converges to the interfacial velocity, indicating an "interface mode."

C. Energy Budget Insights

• Axisymmetric Mode: Driven by surface tension (low $We$) or gas pressure/shear (high $We$, low $El$). Energy is extracted near the interface.
• Helical Mode (High $El$): The Reynolds stress term (REY) becomes the dominant energy source. This term represents the work done by the base flow shear on the disturbance. The analysis reveals that elasticity facilitates a redistribution of disturbance energy from the interface to the jet core, enabling the helical mode to thrive.

D. Validation

• The theoretical $We$–$El$ phase diagram shows excellent quantitative agreement with experimental observations.
• The spatial analysis correctly predicts the emergence of the helical mode at specific thresholds of $We$ and $El$, whereas temporal analysis predicts a much wider range of axisymmetric dominance, contradicting experimental reality.

5. Significance

• Theoretical Advancement: This work resolves the discrepancy between previous theoretical models and experimental observations in viscoelastic jetting by establishing the necessity of spatial analysis for convective flows.
• Mechanistic Insight: It uncovers the Elasticity-Enhanced Shear-Driven Instability (ESI), a mechanism unique to viscoelastic fluids that fundamentally alters how energy is transferred from the mean flow to perturbations.
• Practical Application: The derived phase diagrams and mechanistic understanding provide critical guidance for controlling jet breakup in industrial applications such as inkjet printing, electrospinning, and microfiber fabrication. By tuning elasticity and flow rates, operators can deliberately select between axisymmetric (monodisperse droplets) and helical (fibers or complex patterns) breakup regimes.

In summary, the paper provides a rigorous, experimentally validated framework for understanding viscoelastic jet instability, highlighting the critical role of spatial analysis and identifying a new elasticity-driven instability mechanism that governs the transition to helical flow.