Unit 6: THERMOACOUSTIC INSTABILITY

I. Orientation — Thermoacoustic Coupling in Combustion Systems

Thermoacoustic instability is a self-excited oscillation produced when unsteady heat release from combustion couples constructively with acoustic pressure fluctuations. In gas turbines and aeroengines, the feedback loop involves chamber acoustics, flame dynamics, fuel–air mixing, and boundary impedances.

A. Basic concepts of thermoacoustics

The central principle is that an acoustic disturbance grows when heat is added at a favorable phase of the pressure cycle.

• Feedback mechanism: A pressure or velocity fluctuation perturbs fuel–air mixing, flame position, or burning rate; the resulting heat-release fluctuation generates acoustic waves that return to perturb the flame.
  • Positive feedback produces instability.
  • Negative feedback damps the original disturbance.
• Acoustic variables: Flow quantities are decomposed into mean and fluctuating components:
  
  TEXT

  p(x,t) = p̄(x) + p′(x,t)
  u(x,t) = ū(x) + u′(x,t)
  q̇(x,t) = q̄(x) + q̇′(x,t)

  
  Here, (p) is pressure, (u) is velocity, (\dot q) is volumetric heat-release rate, overbars denote mean values, and primes denote fluctuations.
• Rayleigh criterion: Acoustic energy increases when pressure and heat-release fluctuations are positively correlated over one period:
  
  TEXT

  R = ∫₀ᵀ ∫ᵥ p′(x,t) q̇′(x,t) dV dt > 0

  
  (R) is the Rayleigh integral, (T) is the oscillation period, and (V) is the combustion volume. Growth occurs only when this driving exceeds acoustic and aerodynamic losses.
• Phase condition: For harmonic disturbances, maximum driving occurs when (p′) and (\dot q′) are nearly in phase; heat addition during pressure maxima raises acoustic energy.
• Resonance condition: Strong instability occurs when the combustion-response frequency lies near an acoustic eigenfrequency of the combustor, injector, plenum, or exhaust system.
• Mode quantities: Each mode is characterized by frequency (f) in hertz, angular frequency (\omega=2\pi f), mode shape, damping rate, and growth rate.
• Limit cycle: Linear growth cannot continue indefinitely. Nonlinear flame response, vortex saturation, changing losses, and mode interactions commonly produce a finite-amplitude periodic state called a limit-cycle oscillation.
• Core assumptions in linear analysis:
  • Perturbations are small compared with mean quantities.
  • Mean flow and mean flame are statistically steady.
  • Thermodynamic and flow equations may be linearized.
  • A mode is unstable when combustion-induced gain exceeds total damping.

II. Classification of Instabilities — Frequency, Mode Shape, and Physical Origin

Thermoacoustic instabilities can be classified by acoustic geometry, frequency range, combustion mechanism, and amplitude behavior; several classifications may apply to one oscillation simultaneously.

A. Types of thermoacoustic instabilities

The practical type of instability is identified from its mode shape, characteristic frequency, and dominant feedback path.

• Longitudinal modes: Pressure varies mainly along the combustor axis. For a simple uniform duct of length (L), approximate frequencies are:
  
  TEXT

  fₙ = n c/(2L)          closed–closed or open–open duct
  fₙ = (2n−1)c/(4L)      closed–open duct

  
  (f_n) is the (n)th natural frequency, (c) is sound speed, and (n=1,2,\ldots). Actual aeroengine boundaries require impedance corrections.
• Transverse modes: Pressure varies across the chamber diameter or annulus width. Their high frequencies can cause strong local pressure gradients and uneven loading of liners and injectors.
• Azimuthal modes: Waves propagate around annular combustors as clockwise or counter-clockwise spinning modes, or combine into standing modes. Circumferential order (m) gives an approximate annular relation (f_m\approx mc/(2\pi R)), where (R) is mean annulus radius.
• Radial modes: Pressure changes principally in the radial direction. They are less common at low order but can become relevant in wide annular chambers.
• Low-frequency bulk or Helmholtz modes: The chamber behaves as an acoustic compliance and the inlet or nozzle as an inertive neck:
  
  TEXT

  f_H = c/(2π) √(A/(V Lₑ))

  
  (A) is neck area, (V) is cavity volume, and (L_e) is effective neck length.
• High-frequency acoustic instabilities: Their time scales are comparable to flame and injector response times; transverse and azimuthal modes often belong to this class.
• Mixture-coupled instability: Velocity or pressure disturbances alter equivalence ratio, and the convected mixture fluctuation later changes heat release.
• Vortex-driven instability: Acoustic velocity modulates shear-layer vortex formation; coherent vortices wrinkle the flame and periodically increase burning area.
• Linear and nonlinear regimes:
  1. Linear instability: Infinitesimal disturbances grow exponentially at a rate predicted by eigenvalue analysis.
  2. Nonlinear instability: Finite-amplitude oscillations exhibit saturation, harmonics, hysteresis, intermittency, or switching between modes.

III. Aeroengine Consequences — Operability, Durability, and Efficiency

Combustion instability affects the engine by converting a normally steady combustion process into oscillatory thermal, mechanical, and aerodynamic loading.

A. Effect of thermoacoustic instabilities on aeroengine performance

The severity depends not only on pressure amplitude but also on frequency, mode shape, duration, and the locations of pressure antinodes.

• Combustion efficiency: Oscillatory flame displacement can create incomplete reaction zones, local extinction, and unburned fuel transport, reducing combustion efficiency (\eta_c).
• Pressure loss: Large velocity fluctuations increase mixing and dissipation, raising combustor total-pressure loss:
  
  TEXT

  Δpₜ/pₜ,in = (pₜ,in − pₜ,out)/pₜ,in

  
  (p{t,in}) and (p{t,out}) are combustor inlet and outlet total pressures.
• Thrust and specific fuel consumption: Reduced total pressure and nonuniform turbine-inlet conditions lower turbine work or require additional fuel to maintain thrust, increasing thrust-specific fuel consumption.
• Operability limits: Oscillations can promote lean blowout, flashback in premixed systems, autoignition upstream of the intended flame zone, or compressor–combustor interaction.
• Thermal loading: Periodic flame impingement increases liner-temperature fluctuations and local heat flux, accelerating oxidation, creep, and thermal fatigue.
• Mechanical fatigue: Dynamic pressure excites liners, fuel injectors, transition ducts, and turbine components. Damage is especially severe when acoustic and structural natural frequencies coincide.
• Turbine forcing: Spatially nonuniform exit temperature, entropy waves, and acoustic waves produce periodic forcing on nozzle guide vanes and rotor blades.
• Emissions: Instability changes local equivalence ratio and residence time. Hot excursions can increase nitrogen oxides, while extinction and incomplete combustion can increase carbon monoxide and unburned hydrocarbons.
• Monitoring measures: Engineers track root-mean-square pressure (p'{\mathrm{rms}}), peak-to-peak pressure, spectral peaks, mode growth rate, and normalized amplitude (p'{\mathrm{rms}}/\bar p).

IV. Reduced-Order Prediction — Acoustic-Network Modelling

A one-dimensional model represents the combustor as connected ducts and compact elements, allowing rapid prediction of frequencies, mode shapes, and stability trends.

A. One-dimensional calculation method

The method solves linear acoustic propagation along the mean-flow direction and introduces combustion through a flame-response relation.

• Governing equation: For a uniform, stationary, lossless gas outside the flame:
  
  TEXT

  ∂²p′/∂t² − c² ∂²p′/∂x² = 0

  
  (x) is axial position and (t) is time. Harmonic solutions contain downstream- and upstream-travelling waves.
• Wave representation:
  
  TEXT

  p′(x,t) = [P⁺e^(−ikx) + P⁻e^(ikx)]e^(iωt),   k = ω/c

  
  (P^+) and (P^-) are complex wave amplitudes and (k) is acoustic wavenumber.
• Network construction: Plenum, injector, flame region, combustor, and nozzle are represented by transfer matrices relating pressure and volume-velocity perturbations at component inlet and outlet.
• Boundary conditions: An open end may be approximated by (p′=0), a rigid end by (u′=0), while realistic injectors and nozzles use complex acoustic impedance:
  
  TEXT

  Z(ω) = p̂/û

  
  (Z) is impedance, and hats denote complex harmonic amplitudes.
• Flame transfer function: The normalized heat-release response is commonly related to an upstream velocity disturbance:
  
  TEXT

  q̂/q̄ = F(ω) ûᵣ/ūᵣ,     F(ω) = n(ω)e^(−iωτ)

  
  (n) is interaction gain, (\tau) is time delay, and subscript (r) denotes a reference location.
• Solution sequence:
  1. Divide the geometry into acoustically uniform elements.
  2. Form propagation, area-change, flame, and boundary matrices.
  3. Assemble the global characteristic equation (D(\omega)=0).
  4. Find complex roots and reconstruct mode shapes.
• Interpretation: With the convention (e^{-i\omega t}), (\operatorname{Im}(\omega)>0) indicates growth; frequency is (\operatorname{Re}(\omega)/(2\pi)).
• Limitations: One-dimensional models cannot resolve transverse flame structure, annular mode splitting, local injector coupling, or strongly three-dimensional mean flow.

V. High-Fidelity Prediction — Multidimensional Eigenmode Modelling

Three-dimensional linear analysis resolves realistic geometry and spatially distributed combustion–acoustic coupling while retaining the efficiency of a small-disturbance formulation.

A. Three-dimensional linear combustion instability analysis method

The method computes complex eigenmodes of the linearized flow or acoustic equations coupled to a spatial flame-response model.

• Mean state: Pressure, density, temperature, velocity, and heat release are obtained from measurement, Reynolds-averaged simulation, or time-averaged large-eddy simulation.
• Linearization: Conservation equations are expanded about the mean state and products of perturbations are neglected. Approaches include the inhomogeneous Helmholtz equation and linearized Euler equations.
• Combustion source: Local heat-release perturbation may be prescribed through a distributed flame transfer function:
  
  TEXT

  q̂(x) = G(x,ω) û(xᵣ)

  
  (G) is the complex spatial response and (x_r) is the reference point.
• Discretization: Finite-element, finite-volume, or spectral methods convert the governing equations and impedance boundaries into:
  
  TEXT

  A(ω)φ = 0

  
  (A) is the frequency-dependent system matrix and (\phi) is the eigenmode vector.
• Outputs: Each eigenpair supplies oscillation frequency, growth or decay rate, pressure mode shape, velocity field, and local Rayleigh-index distribution.
• Mode diagnosis: Regions where the cycle-averaged product (\operatorname{Re}(\hat p^\hat q)) is positive drive the mode; negative regions damp it, with (^) denoting complex conjugation.
• Advantages and limitations: Complex annular modes and local coupling are resolved, but predictions depend strongly on mean-state accuracy, flame-response data, boundary impedance, mesh resolution, and linear-response assumptions.

VI. Mitigation Strategies — Breaking the Feedback Loop

Thermoacoustic control aims to reduce flame gain, shift the phase relation, detune acoustic modes, or increase damping without unacceptable penalties in emissions and performance.

A. Control of thermoacoustic instability

Effective control combines design-stage avoidance, passive damping, active actuation, and operational monitoring.

• Combustion modification: Injector staging, altered swirl number, fuel distribution, pilot fraction, or flame position changes the gain (n) and delay (\tau) in the flame response.
• Acoustic detuning: Changing combustor length, plenum volume, nozzle area, or injector impedance shifts eigenfrequencies away from strongly responsive flame frequencies.
• Passive damping: Helmholtz resonators, quarter-wave tubes, perforated liners, and acoustic cavities absorb energy near selected frequencies.
  • Resonators are robust but generally narrow-band.
  • Perforated liners can cover wider bands but may add pressure loss and cooling complexity.
• Mode symmetry breaking: In annular combustors, deliberately nonuniform burner spacing or impedance can prevent coherent coupling of all burners to one azimuthal mode.
• Active control: Sensors detect pressure oscillations, a controller determines phase and amplitude, and fuel valves, secondary-air actuators, or loudspeakers apply an opposing disturbance.
  
  TEXT

  actuator command = C(ω) × measured pressure

  
  (C(\omega)) is the controller transfer function, designed to provide stabilizing gain and phase.
• Operational control: Fuel staging schedules and rapid transit through unstable operating regions reduce exposure, although they do not remove the underlying mode.
• Design verification: Stable operation should be checked across altitude, inlet temperature, fuel composition, deterioration state, and manufacturing tolerances because each changes acoustic and flame-response properties.
• Control trade-offs: A successful measure must preserve relight capability, lean-blowout margin, emissions compliance, durability, pressure loss, weight, and actuator reliability.