Vortex breakdown and its topologies in turbulent flows within a typical swirl combustor geometry

This study utilizes large-eddy simulations to characterize vortex breakdown topologies and dynamics in turbulent, non-reacting swirl combustor flows across varying vane angles, establishing critical swirl-number thresholds, optimal measurement locations, and a comprehensive topology map for predicting vortex breakdown states.

The Gist

Imagine you are trying to keep a campfire burning steadily inside a windy, metal chimney. If the wind is just right, the fire stays put. If it's too weak, the fire blows out. If it's too strong, the fire gets chaotic and unstable.

This paper is about understanding exactly how to control that "wind" inside a special type of chimney called a swirl combustor. These are used in jet engines and power plants to burn fuel efficiently. The key to making them work is creating a swirling motion (like a tornado) that traps hot air in the center, acting as a permanent "hot spot" to keep the flame alive.

Here is the story of what the researchers discovered, explained simply:

1. The "Tornado" Problem (Vortex Breakdown)

When you spin water down a drain, it sometimes suddenly stops spinning and creates a bubble of still water in the middle. In aerodynamics, this is called Vortex Breakdown (VB).

• The Goal: The researchers wanted to know: "How hard do we need to spin the air to create this stable bubble?"
• The Tool: They used a super-powerful computer simulation (like a digital wind tunnel) to watch the air move. They tested different angles for the "fan blades" (vanes) that spin the air.

2. The "Speedometer" Dilemma

To measure how "swirly" the air is, scientists usually use a number called the "Swirl Number." But the researchers found that some old ways of measuring this were like using a broken speedometer that gave different readings depending on where you held it.

• The Fix: They proved that a specific, more complex formula (called the Generic Swirl Number) is the only reliable way to measure the spin.
• The Sweet Spot: They found that if you measure the spin just 40mm (about 1.5 inches) after the fan blades, you get the most accurate reading.
• The Magic Number: They discovered that the "tornado bubble" (which stabilizes the flame) only becomes stable when the swirl number hits a specific range between 0.21 and 0.35. Below that, the bubble is too wobbly to hold a flame. Above that, it's a solid, stable anchor.

3. The Shape of the Tornado (Topology)

Once the bubble forms, what does the "core" of the tornado look like?

• The Single Helix (The Main Act): In almost every case, the core of the swirl looked like a single corkscrew (a single helix). It was like a single, strong rope twisting through the center. This is the dominant shape that keeps the flame stable.
• The Double Helix (The Side Effect): In some cases, a second, weaker corkscrew appeared alongside the first one.
  • For moderate spins: This second rope was actually just a "shadow" or a reflection of the first one, created by the first rope bumping into itself. It wasn't a separate entity; it was just a side effect.
  • For the strongest spin (600° vane angle): The second rope became a real, independent partner. It wasn't just a shadow anymore; it was a distinct, second tornado twisting alongside the first.

4. The "Dancing" Core

The researchers also looked at how these swirls moved over time.

• The Steady Dancer: In some cases, the single corkscrew spun around the center like a figure skater on a perfect, steady rhythm. This is called a "limit cycle." It's predictable and stable.
• The Drunken Dancer: In other cases (very low or very high spin), the corkscrew would spin wildly, stop, start again, and wobble. It was being pushed around by random turbulence, like a leaf caught in a gusty wind. It wasn't a steady dance; it was a chaotic stumble.

5. The "Ghost" Strand

They also noticed a faint, weak strand of swirling air coming from the center of the fan.

• They realized this "ghost strand" wasn't a new tornado. It was actually just the main tornado's shadow, created because the main tornado was wobbling and hitting the air near the center. It proved that the main tornado was the boss, and everything else was just reacting to it.

The Big Takeaway

This paper is a "User Manual" for engineers designing jet engines and power plants.

1. Don't guess the spin: Use the specific "Generic Swirl Number" formula and measure it right after the fan blades.
2. Aim for the middle: You need a swirl number between 0.21 and 0.35 to get a stable flame anchor.
3. Know the shape: Usually, you get one strong corkscrew. If you push the spin too hard, you might accidentally create a second, independent corkscrew, which changes how the engine behaves.

In short, they figured out the exact recipe to create a stable, swirling "fire pit" in a jet engine, ensuring the flame stays lit without blowing out or getting too wild.

Technical Summary

1. Problem Statement

Swirl combustors are critical components in gas turbines and coal reactors, relying on Vortex Breakdown (VB) to create an Internal Recirculation Zone (IRZ) for flame stabilization. However, the onset of VB and the resulting flow topologies (e.g., single-helix vs. double-helix) are complex and depend on various parameters like swirl strength, Reynolds number ($Re$), and expansion ratio ($ER$).

Key gaps in existing literature identified by the authors include:

• Inconsistent Swirl Metrics: Different studies use varying definitions of swirl number (e.g., $SN_c$, $SN_s$), leading to conflicting critical values for VB onset.
• Ambiguity in Topologies: While single-helix (spiral) VB is common, stable double-helix VB has been observed in some recent studies but is rare in others. The mechanism behind this, particularly the role of bluff-body center-bodies (CB), is not fully understood.
• High Reynolds Number Gap: Most existing studies on VB topologies in swirl combustors are conducted at low Reynolds numbers ($Re < 4000$), whereas practical applications operate at $Re \ge 10^4$.
• Lack of Systematic Parametric Study: There is no comprehensive map linking specific swirl strengths to dominant VB topologies in isothermal, high-$Re$ swirl combustors.

2. Methodology

The study employs Large-Eddy Simulation (LES) to investigate non-reacting, isothermal turbulent flows in a canonical swirl combustor geometry.

• Geometry: A lab-scale combustor with a sudden expansion ($ER = 2.0$). The inlet tube contains an axial swirler with eight vanes and a sharp-edged (bluff-body) center-body.
• Flow Conditions:
  • Inlet Reynolds number: $Re_{inlet} = 2.0 \times 10^4$.
  • Working fluid: Air at 300 K.
  • Parametric Variation: Six cases were simulated. A baseline case (45° vane angle) was used for validation. Five additional cases varied the swirler vane angle (17°, 25°, 40°, 50°, 60°) to alter swirl strength while keeping all other geometric and flow parameters constant.
• Numerical Setup:
  • Solver: OpenFOAM-v10 (rhoPimpleFoam).
  • Discretization: Finite Volume Method (FVM) with second-order accuracy.
  • Turbulence Modeling: Implicit LES with a $k$-equation sub-grid scale (SGS) model.
  • Grid Resolution: Hexahedral mesh (1.7 million elements) ensuring >80% of Turbulent Kinetic Energy (TKE) is resolved, satisfying Pope's criterion.
• Analysis Techniques:
  • Swirl Number Evaluation: Comparison of Generic Swirl Number ($SN_g$), Conventional Swirl Number ($SN_c$), and Simplified Swirl Number ($SN_s$).
  • Visualization: $Q$-criterion iso-surfaces to identify Vortex Core (VC) structures.
  • Spectral Analysis: Two-point cross-spectral analysis and four-probe azimuthal FFT to determine coherent modes (azimuthal mode numbers $m=0, 1, 2$) and precession frequencies.
  • Non-linear Interaction: Squared cross-bicoherence analysis to detect quadratic self-interaction between modes.

3. Key Contributions

1. Validation of $SN_g$: The study establishes that the Generic Swirl Number ($SN_g$) is the most robust metric for characterizing swirl strength in isothermal flows. Unlike simplified formulations ($SN_c$, $SN_s$), $SN_g$ remains nearly constant in the inlet tube and decays monotonically downstream, avoiding the spatial variations that lead to inconsistent critical values in literature.
2. Critical Swirl Threshold: The authors determine the critical $SN_g$ for the onset of stable VB (defined by a persistent IRZ in the mean flow) to be between 0.21 and 0.35. This onset occurs at a vane angle of 25°.
3. Topology Mapping: A comprehensive map of VC topologies is provided for $Re \ge 10^4$. The study clarifies that while a single-helix topology is dominant across all cases, double-helix signatures appear but have distinct origins depending on the swirl strength.
4. Mechanism of Double-Helix Formation: The paper distinguishes between two types of double-helix signatures:
   • Quadratic Self-Interaction: For vane angles $\le 50°$, the double-helix signature ($m=2$) is a harmonic generated by the quadratic self-interaction of the dominant single-helix mode ($m=1$).
   • Distinct Hydrodynamic Mode: For the highest vane angle (60°), the double-helix signature is quadratically independent of the single-helix mode, indicating a distinct helical hydrodynamic mode.

4. Key Results

A. Onset of Vortex Breakdown

• Stable VB: First appears at a vane angle of 25° ($SN_g \approx 0.35$). Below this (17°), the IRZ is highly intermittent and does not appear in the mean flow.
• Measurement Location: $SN_g$ measured within 40 mm downstream of the swirler (in the inlet tube) provides the most reliable correlation with flow behavior, minimizing the effects of expansion geometry and downstream decay.

B. Vortex Core (VC) Topologies and Dynamics

• Dominant Topology: Across all investigated vane angles (up to 60°, $SN_g = 0.98$), the single-helix ($m=1$) topology prevails as the primary coherent structure.
• Weakly Coherent Strand: A secondary, weakly coherent strand originates from the wake of the center-body. Its precession frequency matches the primary VC, confirming it is driven by VC precession dynamics rather than independent instability.
• Dynamics by Case:
  • 40°–50° Vane Angles: Exhibit stable limit-cycle oscillations. The single-helix precession is marginally stable, and the double-helix signature is a harmonic ($m=2 \approx 2f_1$) generated via quadratic self-coupling.
  • 25° Vane Angle (Onset): The single-helix and axisymmetric ($m=0$) modes are quadratically independent. The oscillations wax and wane strongly, driven by stochastic forcing on slightly stable modes.
  • 60° Vane Angle (High Swirl): The IRZ merges with the center-body separation bubble, disrupting the global wavemaker. Both single-helix ($m=1$) and double-helix ($m=2$) modes are quadratically independent and sustained by stochastic forcing. The double-helix here represents a distinct hydrodynamic mode, not a harmonic.

C. Precession Frequencies

The lowest coherent precession frequencies increase with vane angle:

• 25°: ~87.5 Hz
• 40°: ~118.75 Hz
• 45°: ~131.25 Hz
• 50°: ~143.75 Hz
• 60°: ~168.75 Hz

5. Significance

• Design Guidelines: The study provides a "topology map" and critical swirl thresholds for engineers designing swirl combustors, ensuring stable flame anchoring by avoiding intermittent VB states.
• Standardization: By validating $SN_g$ as the superior metric and specifying the optimal measurement location (40 mm downstream of the swirler), the paper offers a standardized approach for comparing swirl combustor flows across different geometries and operating conditions.
• Fundamental Insight: The research resolves the debate on double-helix VB, demonstrating that it can be either a non-linear artifact of a single-helix flow or a distinct mode depending on the interaction between the IRZ and the center-body wake.
• Scalability: The findings are applicable to industrial-scale reactors ($Re \ge 10^4$, $ER \ge 1.5$), bridging the gap between low-Re laboratory studies and practical high-Re applications.

In conclusion, this work establishes a rigorous framework for predicting and interpreting vortex breakdown states in isothermal swirl combustors, emphasizing the dominance of single-helix dynamics and the nuanced origins of double-helix signatures.