Wavemaker and endogeneity of gravitationally stretched weakly viscoelastic jets

This paper presents a unified one-dimensional model with full-curvature capillarity and Giesekus stress closure to analyze the global linear stability of gravitationally stretched viscoelastic jets, revealing how moderate elasticity shifts the critical jetting-dripping boundary and identifying the near-nozzle region as the dominant receptive location for instability onset through wavemaker and structural-sensitivity analysis.

The Gist

The Big Picture: The "Stretchy String" Problem

Imagine you are holding a bottle of honey or a thick polymer solution and letting it drip out of the nozzle. Sometimes, it falls as a steady, thin stream (like a faucet running smoothly). Other times, it breaks apart into droplets right at the tip (like a dripping tap).

Scientists have known for a long time how to predict when a normal liquid (like water or oil) will switch from a smooth stream to dripping. They use math to figure out the exact speed and thickness where this happens.

But what happens if the liquid is stretchy? Think of a liquid that acts a bit like a rubber band mixed with water (like a polymer solution or melted plastic). These liquids have "elasticity." When you pull them, they resist and try to snap back, just like a rubber band.

This paper asks: How does this "rubber band" behavior change the rules for when a liquid stream breaks into drops?

The Main Characters

1. The Jet: The thin stream of liquid falling down.
2. Gravity: The force pulling the liquid down, stretching it thinner and thinner as it falls.
3. Surface Tension: The "skin" of the liquid that wants to pull the stream into a ball (a drop).
4. Elasticity (The Rubber Band): The internal "memory" of the liquid. If you stretch it, it holds tension and tries to pull back.

The Discovery: It's Not Just About the Tip

The Old Way (Newtonian Fluids):
For normal liquids, scientists found that the "decision" to break into a drop happens almost entirely at the nozzle (the very top where the liquid comes out). It's like a traffic light at the start of a highway; if the light turns red, the traffic stops immediately. The instability is very localized.

The New Way (Viscoelastic Fluids):
The authors of this paper built a new computer model to simulate these stretchy liquids. They discovered that when you add elasticity, the rules change in two surprising ways:

1. The "Rubber Band" Delays the Break:
   The stretchy tension inside the liquid acts like a safety net. It fights against the surface tension trying to break the stream. This means you can actually run the liquid slower before it starts dripping. The "rubber band" holds the stream together longer, allowing for thinner, more stable threads.

2. The "Decision" Spreads Out:
   This is the most interesting part. In normal liquids, the instability starts at the top. But in stretchy liquids, the "decision" to break isn't just at the nozzle anymore. Because the liquid remembers being stretched, the instability spreads down the entire length of the stream.

   • Analogy: Imagine a line of people holding hands. If the person at the front stumbles, the whole line might wobble, but the reaction is immediate. If the people are holding elastic bungee cords, a stumble at the front sends a wave of tension down the line that takes time to travel and relax. The whole line becomes part of the "wobble," not just the front person.

The "Wavemaker" and the "Echo"

The researchers used a special mathematical tool called a "Wavemaker" analysis to see where the instability is born.

• In Normal Liquids: The "wavemaker" is a tiny, bright spot right at the nozzle. It's like a single drummer starting a beat.
• In Stretchy Liquids: The "wavemaker" gets fuzzy and spreads out downstream. The "beat" is now being made by the whole orchestra, not just the drummer. The elastic tension travels down the jet, creating a feedback loop that involves the liquid far away from the nozzle.

However, the sensitivity (where a tiny nudge can change the whole system) stays near the nozzle. It's like a microphone: even if the sound echoes through a large hall, the microphone is still right at the stage. You have to disturb the liquid at the start to change the rhythm, but the rhythm itself is influenced by the whole stretchy string.

Why Does This Matter?

This isn't just about dripping water. This physics is crucial for:

• 3D Printing: Making tiny, perfect fibers for medical implants.
• Medicine: Creating tiny droplets of medicine for inhalers or drug delivery.
• Manufacturing: Making thin films and fibers for electronics.

By understanding that elasticity spreads the instability down the stream, engineers can design better nozzles and control the flow to make thinner, more uniform threads without them breaking apart too early.

The Bottom Line

The paper shows that when you add "stretchiness" to a falling liquid stream:

1. It becomes harder to break (you can go slower).
2. The instability spreads out along the stream, involving the whole length of the thread, not just the tip.
3. The "rubber band" effect creates a complex dance between the pull of gravity, the squeeze of surface tension, and the snap-back of elasticity.

The authors have provided a new map for engineers to navigate this dance, helping them create better micro-threads and droplets for high-tech applications.

Technical Summary

1. Problem Statement

The paper addresses the stability of liquid jets stretched by gravity, a phenomenon central to drop generation, micro-thread formation, and extensional rheometry. While the transition from steady jetting to self-excited oscillations (dripping) is well-understood for Newtonian fluids, the behavior of weakly viscoelastic fluids (polymer solutions) under similar conditions remains less explored in a global stability context.

Key challenges include:

• Non-parallelism: In gravitationally stretched jets, the radius and velocity vary significantly along the axial direction. Local (parallel-flow) stability analyses often fail to capture the correct onset thresholds and frequencies because they neglect boundary conditions and the spatial development of the base flow.
• Viscoelasticity: Polymer solutions exhibit complex dynamics like "beads-on-a-string" and delayed pinch-off due to extensional stresses. Existing models often treat these locally or assume simplified constitutive laws.
• Mechanism Identification: There is a need to understand where and how viscoelasticity modifies the instability mechanism (specifically the feedback loop) compared to Newtonian jets.

2. Methodology

The authors formulate and analyze a unified one-dimensional (1D) slender-jet model that combines:

• Full-Curvature Capillarity: Unlike many slender-jet approximations that use mean curvature, this model retains the exact interfacial curvature term, which is crucial for accurately predicting stability in strongly stretched jets.
• Viscoelastic Constitutive Law: The model incorporates the Giesekus constitutive equation (with the limit $\alpha=0$ recovering the Oldroyd-B model) to describe polymeric extra-stresses.
• Global Linear Stability Analysis: Instead of local analysis, the authors perform a global eigenvalue analysis on spatially developing base states. This involves:
  • Solving the steady base-flow equations.
  • Linearizing the governing equations (mass, momentum, and stress evolution) about the base state.
  • Discretizing the resulting generalized eigenvalue problem $(A - \omega B)\hat{q} = 0$ using Chebyshev collocation.
• Direct-Adjoint Diagnostics: To interpret the instability mechanism, the authors compute:
  • Wavemaker (Structural Sensitivity): The overlap of direct and adjoint eigenfunctions to identify the spatial region most sensitive to feedback perturbations (the "core" of the oscillator).
  • Endogeneity: A local budget analysis quantifying the contribution of specific physical terms (continuity, momentum, stress) to the eigenvalue (growth rate and frequency).

3. Key Contributions

1. Unified Framework: Development of a global stability framework for gravitationally stretched viscoelastic jets that couples full-curvature capillarity with Giesekus stress closure.
2. Benchmarking: Successful reproduction of Newtonian benchmark results (Rubio-Rubio et al., 2013), validating the numerical implementation of the full-curvature model.
3. Mechanistic Insight via Diagnostics: Application of wavemaker and endogeneity decompositions to viscoelastic jets, revealing how elasticity redistributes the instability feedback loop.
4. Parametric Mapping: Quantification of how the Deborah number ($De$) and viscosity ratio ($\beta$) shift the critical Weber number ($We_c$) and onset frequency ($\omega_i$).

4. Key Results

A. Base Flow and Spectrum

• Newtonian Limit: The model accurately reproduces the marginal spectra and base profiles of Newtonian jets, confirming that the jetting-dripping transition is a global instability selected by the inlet and non-parallel flow.
• Viscoelastic Effects: Introducing elasticity (increasing $De$) stabilizes the jet. For a fixed Weber number, increasing $De$ shifts the leading eigenvalue to more negative growth rates (stabilizing) and reduces the oscillation frequency.
• Polymeric Tension: The base state develops significant polymeric tensile stress ($T_0 = \tau_{zz} - \tau_{rr}$) downstream, which competes with capillary thinning.

B. Wavemaker and Endogeneity (The Instability Mechanism)

• Newtonian Case: The wavemaker is sharply localized near the nozzle inlet. The endogeneity is dominated by the continuity/capillary pathway.
• Viscoelastic Case:
  • Spatial Broadening: As $De$ increases, the wavemaker region broadens significantly downstream. The instability feedback loop is no longer confined to the nozzle but extends over a larger portion of the jet due to the advection and relaxation of polymeric stresses.
  • Endogeneity Decomposition: A distinct polymeric stress contribution ($E_{\tau_{zz}}$) emerges. The marginal state results from a balance between the capillary/kinematic pathway and an elastic-stress feedback pathway.
  • Receptivity: Despite the broadening of the wavemaker, the adjoint eigenfunction (receptivity) remains strongly localized near the inlet. This implies that while the feedback loop is distributed, the selection of the mode is still controlled by perturbations near the nozzle.

C. Marginal Conditions

• Critical Flow Rate: Increasing $De$ lowers the critical Weber number ($We_c$) required for instability. This means viscoelasticity delays the onset of dripping, allowing steady jetting at lower flow rates.
• Frequency: The onset frequency ($\omega_i$) decreases with increasing $De$.
• Viscosity Dependence: The viscoelastic shift is most pronounced at moderate Bond numbers and lower Kapitza numbers (higher viscosity). In very low-viscosity regimes (high Kapitza), the shift is mild, suggesting the instability is dominated by rapid inertio-capillary balances near the inlet where polymers have less time to accumulate stress.

5. Significance

• Theoretical Advancement: The study moves beyond local stability analysis to provide a rigorous global framework for viscoelastic jets, correctly capturing the interplay between gravity, capillarity, and elasticity.
• Physical Interpretation: It clarifies that viscoelasticity does not just shift eigenvalues but fundamentally alters the spatial structure of the instability. The "memory" of the polymer stretches the feedback loop downstream, creating a distributed capillary-elastic coupling.
• Practical Implications:
  • Process Control: Understanding that the near-nozzle region is the dominant receptive location suggests that controlling inlet conditions (flow rate, perturbations) is critical for managing jet stability in viscoelastic fluids.
  • Micro-manufacturing: The ability to predict how elasticity shifts the jetting-dripping boundary aids in optimizing the production of micro-droplets and sub-micron threads, allowing for the generation of smaller features at lower flow rates.
  • Modeling: The results highlight that reduced-order models for viscoelastic jets must account for the distributed nature of the feedback loop, not just local stress effects.

In summary, the paper demonstrates that weak elasticity in gravitationally stretched jets transforms the instability from a localized, capillary-driven nozzle phenomenon into a spatially distributed, capillary-elastic feedback system, while maintaining the nozzle as the primary control point for mode selection.