Triangular instability of a strained Batchelor vortex

This study combines theoretical analysis and numerical simulations to demonstrate that a stationary triangular strain field induces resonant instability in a Batchelor vortex, where the introduction of axial flow reduces critical layer damping to activate additional unstable mode pairs and ultimately shifts the dominant instability from a specific mode combination to one involving the first branches of both azimuthal wavenumbers.

The Gist

Imagine a giant, swirling tornado in the sky. In physics, we call this a vortex. Usually, these swirls are perfectly round, like a spinning top. But in the real world—like behind an airplane wing or a wind turbine—these swirls are often surrounded by other smaller swirls.

This paper is about what happens when a big, central swirl (the "Hub") is surrounded by three smaller swirls (the "Satellites") arranged in a perfect triangle.

Here is the story of the paper, broken down into simple concepts:

1. The Setup: The Dance of Three

Imagine a ballroom dance.

• The Hub: A big, strong dancer spinning in the center. In this study, this dancer is also walking forward (this is the axial flow).
• The Satellites: Three smaller dancers spinning around the big one, holding hands in a triangle.
• The Strain: As the three small dancers spin, they pull and push on the big dancer. This creates a "triangular squeeze" or strain on the big dancer's shape.

The scientists wanted to know: Does this triangular squeeze make the big dancer wobble and eventually fall apart?

2. The Problem: The "Critical Layer" Wall

In the past, scientists studied similar setups but without the big dancer walking forward (no axial flow). They found that the big dancer could only wobble in one specific way: a pair of waves with specific shapes (called wavenumbers -1 and 2).

Why only that one? Because of a "wall" called the Critical Layer.

• Think of the Critical Layer as a soundproof wall inside the vortex.
• If the big dancer tries to wobble in a different pattern (like shapes 0 and 3, or 1 and 4), the soundproof wall absorbs the energy, and the wobble dies out immediately. It's like trying to shout through a thick concrete wall; the sound never gets through.

3. The Discovery: The "Walking" Breaks the Wall

The big discovery in this paper is that when the big dancer starts walking forward (axial flow), the soundproof wall gets a hole in it.

• Before: Only one specific wobble pattern could survive.
• After: As the walking speed increases, the "wall" becomes thinner. Suddenly, new patterns (like 0 and 3, or 1 and 4) can get through!
• The Result: The big dancer can now wobble in many more dangerous ways, not just the one way it used to.

4. The "Switch" in Leadership

Here is the most interesting twist. The scientists found that the "most dangerous" wobble changes depending on how fast the dancer is walking.

• Slow Walk: The most dangerous wobble is a specific mix of the second branch of one wave and the first branch of another.
• Fast Walk: As the walking speed increases, that old champion gets weaker. A new champion takes over: a mix of the first branches of both waves.
• The Takeaway: Once the walking speed passes a certain point, this new "First Branch" wobble becomes the king of instability. It stays the most dangerous across a wide range of speeds and conditions.

5. Why Should We Care? (The Real World)

Why do we care about a dancing vortex?

• Airplanes: The wake behind a plane is full of these swirling vortices. If they become unstable too quickly, they might break up faster, which is actually good for safety (less turbulence for the next plane).
• Wind Turbines & Ships: The hubs of wind turbines and ship propellers create these exact triangular patterns. Understanding how they wobble helps engineers design better, safer, and more efficient machines.

Summary Analogy

Think of the vortex as a rubber band being stretched by three fingers (the satellites).

• Without walking: The rubber band can only snap in one specific direction because the material is too stiff in other directions.
• With walking: The rubber band gets heated up (by the axial flow). It becomes stretchy. Now, it can snap in many different directions, and the direction it snaps in changes depending on how hot (fast) it gets.

The Bottom Line:
This paper proves that movement changes the rules of stability. By adding a forward flow to a swirling vortex, we unlock new ways for it to become unstable, and we can predict exactly which "wobble" will win the race to break the vortex apart. This helps us understand the hidden chaos in the wakes of everything from airplanes to wind turbines.

Technical Summary

Here is a detailed technical summary of the paper "Triangular instability of a strained Batchelor vortex" by Ayapilla, Hattori, and Le Dizès.

1. Problem Statement

The study investigates the triangular instability of a Batchelor vortex (a vortex with a Gaussian vorticity profile and an axial flow component) subjected to a stationary, triangular strain field. This strain field is generated by three symmetrically positioned "satellite" vortices surrounding a central "hub" vortex.

• Context: While elliptic instability (driven by quadrupolar strain) and curvature instability are well-studied, multipolar instabilities (driven by higher-order strain) are less understood. Previous work by the authors (AHL25) established triangular instability in a Lamb-Oseen vortex (no axial flow), showing it was restricted to the azimuthal wavenumber pair $(m_A, m_B) = (-1, 2)$.
• Gap: Real-world vortices (e.g., aircraft trailing vortices, rotor wakes) often possess a non-zero axial flow component. The Batchelor vortex is the standard model for such flows. The paper aims to determine how axial flow influences the resonance conditions, mode selection, and growth rates of triangular instability.

2. Methodology

The authors employ a combined approach of theoretical linear stability analysis and numerical simulations.

A. Numerical Framework (DNS)

• Base Flow: The system consists of a central Batchelor vortex and three satellite Lamb-Oseen vortices. The satellite vortices induce a triangular strain on the hub. The base flow is computed by solving the incompressible Navier-Stokes equations in a cylindrical coordinate system until a quasi-steady state is reached.
• Linear Stability: The linearized Navier-Stokes equations are solved around this base flow using Direct Numerical Simulation (DNS).
• Parameters:
  • Reynolds number: $Re = 10^4$.
  • Straining strength: $\epsilon = (a/R)^3 = 0.008$ (where $a$ is core radius, $R$ is distance to satellites).
  • Axial flow parameter ($W_0$): Varied from $0$to$-0.4$.
• Goal: Identify the most unstable modes, their growth rates, frequencies, and spatial structures.

B. Theoretical Framework

• Resonance Condition: The instability arises from the resonant coupling of two quasi-neutral Kelvin modes of the unstrained vortex with the background strain. The resonance condition requires:
  • $\omega_A = \omega_B$ (frequencies match)
  • $k_A = k_B$ (axial wavenumbers match)
  • $m_B - m_A = 3$ (azimuthal wavenumber difference matches the strain order).
• Asymptotic Analysis: Using a multiscale analysis (similar to Moore & Saffman, 1975), the authors derive an expression for the growth rate of the resonant pair.
• Kelvin Mode Classification: The study utilizes the WKB asymptotic framework (Le Dizès & Lacaze, 2005) to classify modes into:
  • Regular vs. Singular: Singular modes encounter a critical layer (where phase speed matches local flow speed), leading to damping.
  • Core vs. Ring: Core modes are localized near the vortex axis; ring modes are localized away from the center.
• Axial Flow Effect: The theory analyzes how axial flow shifts dispersion curves, potentially allowing modes previously damped by critical layers to become resonant.

3. Key Contributions

1. Extension to Batchelor Vortex: First comprehensive study of triangular instability in a Batchelor vortex, bridging the gap between idealized Rankine/Lamb-Oseen models and realistic wake flows.
2. Role of Axial Flow: Demonstrates that axial flow fundamentally alters the instability landscape by:
   • Reducing critical-layer damping for specific modes.
   • Enabling new resonant pairs beyond $(-1, 2)$, specifically $(0, 3)$, $(1, 4)$, and $(2, 5)$.
   • Shifting the dominant unstable mode from the second branch of $m=-1$ (in no-flow cases) to the first branch of both $m=-1$ and $m=2$ as axial flow increases.
3. Mode Transition: Identifies the transition of certain modes from core modes to ring modes as axial flow strength increases, characterized by a shift in eigenfunction localization away from the vortex center.
4. Validation: Provides rigorous validation between theoretical growth rate predictions and DNS results, showing excellent agreement for growth rates and frequencies.

4. Key Results

A. Resonant Pairs and Axial Flow

• No Axial Flow ($W_0=0$): Only the pair $(-1, 2)$ is unstable. The most unstable mode is $(-1, 2, [2, 1])$ (second branch of $m=-1$, first of $m=2$). Other pairs like $(0, 3)$ are heavily damped by critical layers.
• With Axial Flow ($W_0 < 0$):
  • As $|W_0|$ increases, critical-layer damping decreases for modes that were previously suppressed.
  • New resonant pairs become unstable: $(0, 3)$, $(1, 4)$, and $(2, 5)$.
  • Dominance Shift: The mode $(-1, 2, [2, 1])$ becomes less unstable as $|W_0|$ increases. It is eventually overtaken by the mode $(-1, 2, [1, 1])$ (first branch of both wavenumbers), which becomes the dominant unstable mode across a wide range of axial flow strengths.

B. Growth Rates and Parameters

• Growth Rate Trends: The maximum growth rate for the dominant mode $(-1, 2, [1, 1])$ peaks at $|W_0| \approx 0.2\text{--}0.3$.
• Reynolds Number: At $Re=10^4$, the triangular instability modes are significantly more unstable than viscous center modes (which only become unstable for $W_0 < -0.33$ for $m=1$).
• Straining Strength: Instability persists even for very weak straining ($\epsilon \approx 0.001$), though growth rates scale with $\epsilon^2$.
• Ring Modes: For pairs $(0, 3)$, $(1, 4)$, and $(2, 5)$, modes involving the first branch of $m_A$ transition to ring modes at higher $|W_0|$. These ring modes exhibit narrower growth-rate bands and lower magnitudes compared to core modes.

C. Structural Evolution

• Vorticity Structure: At low $|W_0|$, the unstable mode exhibits a spiral pattern due to critical-layer waves. As $|W_0|$ increases, the critical layer influence weakens, the spiral disappears, and the mode adopts a more regular structure.
• Amplitude Ratio: The ratio of the amplitude of the $m_B$ component to the $m_A$ component increases monotonically with $|W_0|$, eventually exceeding unity, indicating the $m_B$ mode becomes the dominant contributor.

5. Significance and Implications

• Rotor Wakes: The findings are directly applicable to the wakes of three-bladed rotors (e.g., wind turbines, ship propellers). The hub vortex in such systems is subjected to a third-order (triangular) strain from the tip/root vortices. This study suggests that triangular instability is a viable mechanism for transition to turbulence in these wakes, particularly when axial flow is present.
• Competition with Elliptic Instability: While elliptic instability (second-order strain) is generally dominant, triangular instability can become the primary mechanism if the flow is near a third-order symmetrical state or if the elliptic instability is suppressed.
• Future Directions: The paper highlights the need to study non-stationary strain fields (rotating satellites) and helical vortex structures, which are more representative of real rotor wakes. It also notes that nonlinear effects and interactions with global long-wavelength instabilities (e.g., Crow instability) require further investigation.

In summary, this work establishes that axial flow is a critical parameter that unlocks new resonant pathways for triangular instability in Batchelor vortices, shifting the dominant instability mechanism and altering the spatial structure of the perturbations, with significant implications for the stability analysis of rotating machinery wakes.