disturbance propagation and generation in reacting flows · comments on length scales • a coustic...
TRANSCRIPT
1
Disturbance Propagation and Generation in Reacting Flows
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
2
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
3
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
4
Motivation – Combustion Instabilities
• Large amplitude acoustic oscillations driven by heat release oscillations
• Oscillations occur at specific frequencies, associated with resonant modes of combustor
Heat release Acoustics
Video courtesy of S. Menon
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
5
Resonant Modes – You can try this at home
• Helmholtz Mode
• 190 Hz
• Longitudinal Modes
• 1,225 Hz
• 1,775 Hz
• Transverse Modes
• 3,719 Hz
• 10,661 Hz
Slide courtesy of R. Mihata, Alta Solutions
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
6
Key Problem: Flame Sensitive to Acoustic Waves
Video from Ecole Centrale – 75 Hz, Courtesy of S. Candel
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
7
Key Problem: Combustor System Sensitive to Acoustic Waves
Rubens Tube
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
8
Thermo-acoustics
• Rijke Tube (heated gauze in tube)
• Self-excited oscillations in cryogenic tubes
• Thermo-acoustic refrigerators/heat pumps
Purdue’s Thermoacoustic Refrigerator
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
9
Liquid Rockets
• F-1 Engine – Used on Saturn V – Largest thrust engine
developed by U.S – Problem overcome
with over 2000 (out of 3200) full scale tests
From Liquid Propellant Rocket Combustion Instability, Ed. Harrje and Reardon, NASA Publication SP-194
Injector face destroyed by combustion instability, Source: D. Talley
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
10
Solid Rockets
From Blomshield, AIAA Paper #2001-3875
• Examples: – Space shuttle booster- 1-3 psi
oscillations (1 psi = 33,000 pounds of thrust)
Adverse effects –thrust oscillations, mean pressure changes, changes in burning rates
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
11
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
12
Small Amplitude Propagation in Uniform, Inviscid Flows
• Assume infinitesimal perturbations superposed upon a spatially homogenous background flow.
• We will also introduce two additional assumptions: (1) The gas is non-reacting and calorically perfect, implying that the
specific heats are constants. (2) Neglect viscous and thermal transport.
Decomposition Approach – Decompose variables into the sum of a base and fluctuating
component; e.g., (2.9)
0 1
0 1
0 1
( , ) ( , )( , ) ( , )( , ) ( , )
p x t p p x tx t x t
u x t u u x tρ ρ ρ
= +
= += +
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
13
Perturbing the Continuity Equation
• Consider the continuity equation:
(2.10)
• Ds/Dt is identically zero because of the neglect of molecular transport terms and chemical reaction. Therefore …
• Expanding each variable into base and fluctuating components yields:
(2.12)
1 1 1
p
D D Dp uDt c Dt p Dtρ
ρ γ− = − = ∇⋅
s
1 Dp up Dtγ
− = ∇ ⋅
( ) ( )10 1 1 1
0 1
1 p u u p up p tγ
∂ − + + ⋅∇ = ∇⋅ + ∂
(2.11)
14
Linearizing Continuity
• Continuity Equation Expanded: (2.13)
• Linearize by neglecting products of perturbations. If perturbations are small in amplitude then products of perturbations are negligible.
• Definition of the substantial derivative, D0/Dt.
(2.15)
10 1 1 1 0 1 1 1
1 p u p u p p u p utγ
∂ − + ⋅∇ + ⋅∇ = ∇⋅ + ∇ ⋅ ∂
0 11 1
0
1 D p up Dtγ
− = ∇ ⋅ ≡ Λ
( ) ( ) ( )00
D uDt t
∂= + ⋅∇∂
(2.14)
15
Summary of Disturbance Equations
• Vorticity: (2.18)
• Acoustic: (2.20)
• Entropy:
0 1 0DDtΩ
=
22 20 10 12 0D p c p
Dt− ∇ =
0 1 0DDt
=s (2.27)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
16
Canonical Decomposition
• Further decompose perturbations by their origin
(2.28)
• Examples: – vorticity fluctuations induced by entropy fluctuations. – pressure fluctuations induced by vorticty fluctuations.
• Dynamical equations are linear and can be decomposed into subsystems
1 1 1 1
1 1 1 1
1 1 1 1p p p p
Λ Ω
Λ Ω
Λ Ω
Ω = Ω +Ω +Ω= + += + +
s
s
s
s s s s
1Ω
s1p Ω
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
17
Oscillations from Vorticity
• Oscillations associated with vorticity mode:
(2.29)
(2.30)
(2.31)
– Vorticity, and induced velocity, fluctuations are convected by the mean
flow. – Vorticity fluctuations induce no fluctuations in pressure, entropy,
temperature, density, or dilatation.
0 1 0DDt
ΩΩ=
1 1 1 1 0p T ρΩ Ω Ω Ω= = = =s
0 1 0D uDt
Ω =
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
18
Oscillations from Acoustics
• Oscillations associated with acoustic mode:
(2.32)
• The density, temperature, and pressure fluctuations are locally and algebraically related through their isentropic relations
22 20 102 01
D p c pDt
ΛΛ− ∇ =
1 1 0Λ ΛΩ = =
s01
1 12 20 0
pcp Tc c
ρρ Λ
Λ Λ= =
10 1
Du pDt
ρ ΛΛ= −∇
(2.33)
(2.34)
(2.35)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
19
Oscillations from Entropy
• Oscillations associated with entropy mode:
(2.37)
(2.38)
• Entropy oscillations do not excite vorticity, velocity, pressure, or dilatational disturbances
• They do excite density and temperature perturbations
0 11 1 1 0D p u
Dt= = Ω = =
ss s s
s
0 01 1 1
0p
Tc Tρ ρρ = − = −s s ss
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
20
Comments on the Decomposition
• In a homogeneous, uniform flow, these three disturbance modes propagate independently in the linear approximation.
• three modes are decoupled within the approximations of this analysis - vortical, entropy, and acoustically induced fluctuations are completely independent of each other.
• For example, velocity fluctuations induced by vorticity and acoustic disturbances, and are independent of each other and each propagates as if the other were not there.
• Moreover, there are no sources or sinks of any of these disturbance modes. Once created, they propagate with constant amplitude.
1u Ω
1u Λ
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
21
Comments on Length Scales
• Acoustic disturbances propagate at the speed of sound.
• Vorticity and entropy disturbances convect at bulk flow velocity, .
• Acoustic properties vary over an acoustic length scale, given by
• Entropy and vorticity modes vary over a convective length scale, given by .
• Entropy and vortical mode “wavelength” is shorter than the acoustic wavelength by a factor equal to the mean flow Mach number.
0u
0c fλΛ =
0c u fλ =
0 0c u c Mλ λΛ = =
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
22
Comments on Acoustic Disturbances
• Acoustic disturbances, being true waves, reflect off boundaries, are refracted at property changes, and diffract around obstacles.
FlameAcoustic wave
Image of instantaneous pressure field and flame front of a sound wave incident upon a turbulent flame from the left at three successive times. Courtesy of D. Thévenin.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
23
Reactants
Transverse plane wave
Refracted acoustic
wave
Comments on Acoustic Disturbances
• PIV data shows example of – refraction in a combustor environment – Simultaneous presence of acoustic and vortical disturbances
Illustration of acoustic refraction effects. Data courtesy of J. O'Connor
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
24
Summary
Acoustic Entropy Vorticity Pressure X Velocity X X Density X X
Temperature X X Vorticity X Entropy X
25
Example: Effects of Simultaneous Acoustic and Vortical Disturbances
• Consider superposition of two disturbances with different phase speeds:
• For simplicity, assume
• Velocity field oscillates harmonically at each point as cos(ωt)
• Amplitude of these oscillations varies spatially as due to interference
( )1 1, 1,0 ,0
, cos cosx
x xu x t u u t tc u
ω ωΛ Ω Λ Ω
= + − + −
= A A
Λ Ω= =A A A
( ) 0 0 0 01
0 0 0 0
, 2 cos cos2 2c u c uu x t x t x
c u c uω
− −= −
A
0 0
0 0
cos2c u x
c u −
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
26
Some Data Illustrating This Effect
• Fit parameter:
0 0.5 1 1.50
0.02
0.04
0.06
0.08
0.1
x/λc
|un′
|/U0
MeasurementsFit-u1,Λ+u1,Ω
u1,Λonly
u1, Ω onlyu
u′
0 0.5 1 1.5-5
-4
-3
-2
-1
0
1
x/λc∠
u n′/ π
Measurements
Fit-u1,Λ+u1,Ω
u1,Λonly
u1, Ω only
0.6Λ Ω =A A
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
27
Modal Coupling Processes
• We just showed that the small amplitude canonical disturbance modes propagate independently within the fluid domain in a homogeneous, inviscid flow.
• These modes couple with each other from: – Boundaries
• e.g., acoustic wave impinging obliquely on wall excites vorticity and entropy
– Regions of flow inhomogeneity • e.g., acoustic wave propagating through shear flow generates vorticity • Accelerating an entropy disturbance generates an acoustic wave
– Nonlinearities • e.g., large amplitude vortical disturbances generate acoustic waves (jet
noise)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
28
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
29
Energy Density and Energy Flux Associated with Disturbance Fields
• Although the time average of the disturbance fields may be zero,
they nonetheless contain non-zero time averaged energy and lead to energy flux whose time average is also non-zero. – Ex:
• Consider the energy equation.
(2.50)
(2.51)
E It
∂+∇ ⋅ = Φ
∂
VV A V
d EdV I ndA ddt
+ ⋅ = Φ∫∫∫ ∫∫ ∫∫∫
( )1 112kinE u uρ= ⋅
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
30
Acoustic Energy Equation
• Consider the acoustic energy equation, assuming: – Entropy and vorticity fluctuations are zero – Zero mean velocity, homogenous flow. – Combustion process is isomolar
(2.52)
(2.53)
(2.54)
(2.55)
( ) ( ) ( )21
0 1 1 1 1 1 120 0 0
11 12 2
pu u p u p qt c p
γρ
ρ γ− ∂
⋅ + +∇ ⋅ = ∂
( )1 1
0
1p q
pγγΛ
−Φ =
1 1I p uΛ =
( )21
0 1 1 20 0
1 12 2
pE u uc
ρρΛ = ⋅ +
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
31
Time Average of Products
• Time average of product of two fluctuating quantities depends on their relative phasing
• Example ( ) ( ) 1sin sin cos2
t tω ω θ θ+ =
( ) ( )sin sint tω ω θ+
t
t
t
θ
o0θ =
o90θ =
o180θ =Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
32
Comments on Acoustic Energy
• The energy density, , is a linear superposition of the kinetic energy associated with unsteady motions, and potential energy associated with the isentropic compression of an elastic gas.
• The flux term, , reflects the familiar “pumping work” done by pressure forces on a system.
• The source term, , shows that unsteady heat release can add or remove energy from the acoustic field, depending upon its phasing with the acoustic pressure.
EΛ
1 1p u
ΛΦ
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
33
Rayleigh’s Criterion
• Rayleigh’s Criterion states that unsteady heat addition locally adds energy to the acoustic field when the phases between the pressure and heat release oscillations is within ninety degrees of each other.
• Conversely, when these oscillations are out of phase, the heat addition oscillations damp the acoustic field.
• Figure illustrates that the highest pressure amplitudes are observed at conditions where the pressure and heat release are in phase.
/p p′-180
-135
-90
-45
0
45
90
135
180
0.00 0.04 0.08 0.12
pqθData courtesy of K. Kim and D. Santavicca.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
34
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
35
Linear and Nonlinear Stability of Disturbances
• The disturbances which have been analyzed arise because of underlying instabilities, either in the local flow profile or to the coupled flame-combustor acoustic systems (such as thermoacoustic instabilities).
• Consider a more general study of stability concepts by considering the time evolution of a disturbance
FA and FD denote processes responsible for amplification and damping of the disturbance, respectively.
( )DA
d t F Fdt
= −A (2.63)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
36
Linear Stable/Unstable Systems
• In a linearly stable/unstable system, infinitesimally small disturbances decay/grow, respectively.
• To illustrate these points, we can expand the functions FA and FD around their A=0 values in a Taylor’s series:
(2.64)
(2.65)
(2.66) A
mpl
ifica
tion/
Dam
ping
FA
A ALC
1
1εD
εA
FD
,A A A NLF Fε= +A
,D D D NLF Fε= +A
0A
AFε =∂
=∂ AA
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
37
Comments on Linear Behavior
• The amplification and damping curves intersect at the origin, indicating that a zero amplitude oscillation is a potential equilibrium point.
• However, this equilibrium point is unstable, since and any small disturbance makes FA larger than FD resulting in further growth of the disturbance.
• If the A=0 point is an example of an “attractor” in that disturbances are drawn toward it
• Linearized solution:
(2.67)
A Dε ε>
1( ) ( 0)exp(( ) )DAt t tε ε= = −A A
A Dε ε<
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
38
Comments on Nonlinear Behavior
• This linearized solution may be a reasonable approximation to the system dynamics if the system is linearly stable (unless the initial excitation, A(t=0) is large).
• However, it is only valid for small time intervals when the system is unstable, as disturbance amplitudes cannot increase indefinitely.
• In this situation, the amplitude dependence of system amplification/damping is needed to describe the system dynamics.
• The steady state amplitude is stable at ALC=0 because amplitude perturbations to the left (right) causes FA to become larger (smaller) than FD, thus causing the amplitude to increase (decrease).
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
39
Example of Stable Orbits: Limit Cycles
• In many other problems, is used to describe the amplitude of a fluctuating disturbance; for example,
• Limit cycle is example of orbit that encircles an unstable fixed point
– A stable limit cycle will be pulled back into this attracting, periodic orbit even when it is slightly perturbed.
(2.69) 0( ) ( )cos( )p t p t tω− =A
( )tA pt
∂∂
( )p t
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
40
Variation in System Behavior with Change in Parameter
• Consider a situation where some combustor parameter is systematically varied in such a way that εA increases while εD remains constant.
Am
plifi
catio
n/D
ampi
ng
A Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
41
Supercritical Bifurcations
• The εA=εD condition separates two regions of fundamentally different dynamics and is referred to as a supercritical bifurcation point. – Note smooth monotonic variation of limit cycle amplitude with parameter
StableUnstable
Am
plitu
de,A
εA- εD
00
Pressure amplitude
Nozzle velocity, m/s
pp′
0
0.005
0.01
0.015
0.02
18 21 24 27 30
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
42
Nonlinearly Unstable Systems
• Small amplitude disturbances decay in a nonlinearly unstable system, but disturbances with amplitudes exceeding a critical value, Ac, will grow.
• This type of instability is sometimes referred to as subcritical.
• Other examples: – hydrodynamic stability of shear flows
without inflection points – Certain kinds of thermoacoustic
instabilities in combustors. • historically referred to as “triggering”
in rocket instabilities
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
43
Subcritical Instability
• Figure provides example of amplitude dependence of FA and FD that produces subcritical instability.
• In this case, the system has three equilibrium points where the amplification and dissipation curves intersect.
• All disturbances with amplitudes A<AT return to the stable solution A=0 and disturbances with amplitudes A>AT grow until their amplitude attains the value A=ALC.
• Consequently, two stable solutions exist at this operating condition. The one observed at any point in time will depend upon the history of the system.
Am
plifi
catio
n/D
ampi
ng
AALC
1εD
FD
1εA
FA
AT
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
44
Variation in System Behavior with Change in Parameter
Am
plifi
catio
n/D
ampi
ng
A
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
45
Subcritical Bifurcation Diagram A
mpl
itude
, A
εΑ-εD
StableUnstable
AT
00
0.005
0.01
0.015
13 13.5 14 14.5 15 15.5Nozzle velocity, m/s
pp′
Pressure amplitude
• If a system parameter is monotonically increased to change the sign of from a negative to a positive value, the system’s amplitude jumps discontinuously from A=0 to A= ALC at .
• Note hysteresis as well
A Dε ε−
0A Dε ε− =
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
46
• Introduction
• Flashback and Flameholding
• Flame Stabilization and Blowoff
• Combustion Instabilities – Motivation – Disturbance propagation, amplification, and stability – Acoustic Wave Propagation Primer – Unsteady Heat Release Effects and Thermoacoustic Instability
• Flame Dynamics
Course Outline
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
47
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
48
Overview 1 of 2
• Acoustic Disturbances: – Propagate energy and information through the medium without requiring
bulk advection of the actual flow particles. – Details of the time averaged flow has relatively minor influences on the
acoustic wave field (except in higher Mach number flows). – Acoustic field largely controlled by the boundaries and sound speed
field.
• Vortical disturbances – Propagate with the local flow field. – Highly sensitive to the flow details. – No analogue in the acoustic problem to the hydrodynamic stability
problem.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
49
Overview 2 of 2
• Some distinctives of the acoustics problem: – Sound waves reflect off of boundaries and refract around bends or other
obstacles. • Vortical and entropy disturbances advect out of the domain where they are
excited. – An acoustic disturbance in any region of the system will make itself felt
in every other region of the flow. – Wave reflections cause the system to have natural acoustic modes;
oscillations at a multiplicity of discrete frequencies.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
50
Traveling and Standing Waves
• The acoustic wave equation for a homogeneous medium with no unsteady heat release is a linearized equation describing the propagation of small amplitude disturbances.
• We will first assume a one-dimensional domain and neglect mean flow, and therefore consider the wave equation:
(5.1)
2 221 102 2 0p pc
t x∂ ∂
− =∂ ∂
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
51
Harmonic Oscillations
• We will use complex notation for harmonic disturbances. In this case, we write the unsteady pressure and velocity as
(5.13)
(5.14)
• As such, the one-dimensional acoustic field is given by:
(5.15)
( )( )1 1ˆReal ( , , ) expp p x y z i tω= −
( )( )1 1Real ( , , ) expu u x y z i tω= −
( ) ( )( ) ( )( )1 Real exp exp expp ikx ikx i tω= + − −A B
( ) ( )( ) ( )( ),10 0
1 Real exp exp expxu ikx ikx i tc
ωρ
= − − −A B (5.16)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
52
Time Response of Harmonically Oscillating Traveling Waves
• The disturbance field consists of space-time harmonic disturbances propagating with unchanged shape at the sound speed.
• An alternative way to visualize these results is to write the pressure in terms of amplitude and phase as:
(5.17) ( ) ( ) ( )( )1 1ˆ ˆ expp x p x i xθ= −
p1
x
tTemporal variation of harmonically
varying pressure for rightward moving wave
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
53
Amplitude and Phase of Traveling Wave Disturbances
• The magnitude of the disturbance stays constant but the phase decreases/increases linearly with axial distance.
• The slope of these lines are .
• Harmonic disturbances propagating with a constant phase speed have a linearly varying axial phase dependence, whose slope is inversely proportional to the disturbance phase speed.
BA
leftward movingrightward moving
1p
x
x
0/ cω
0/ cω1
1θ
Spatial amplitude/phase variation of harmonically varying
acoustic disturbances.
0/d cdxθ ω=
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
54
Standing Waves
• Consider next the superposition of a left and rightward traveling wave of equal amplitudes, A=B, assuming without loss of generality that A and B are real.
(5.18)
• Such a disturbance field is referred to as a “standing wave”.
• Observations: – amplitude of the oscillations is not spatially constant, as it was for a single
traveling wave. – phase does not vary linearly with x, but has a constant phase, except across the
nodes where it jumps 180 degrees. – pressure and velocity have a 90 degree phase difference, as opposed to being
in-phase for a single plane wave.
( ) ( )1( , ) 2 cos cosp x t kx tω= A
( ) ( ),10 0
1 2 sin sinxu kx tc
ωρ
= A (5.19)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
55
Behavior of Standing Waves
Spatial dependence of pressure (solid) and velocity (dashed) in a standing wave.
p1ρ0c0ux,1
kx/2π
1
1,5
2,8 2,4
3,7
3
4,6
5
6,8
7
(a)
(degrees) 180o
kx/2π
270o
90o
0o
(c)
kx/2π
2A(b)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
56
One Dimensional Natural Modes
• Consider a duct of length L with rigid boundaries at both ends, . Applying the boundary condition at x=0 implies that A=B, leading to:
(5.62)
(5.64)
• These natural frequencies are integer multiples of each other, i.e., f2=2f1, f3=3f1, etc.
(5.65)
(5.66)
( ) ( ),1 ,10, , 0x xu x t u x L t= = = =
( ) ( )( ),10 0
1 Real 2 sin expxu i kx i tc
ωρ
= −A
nkLπ
=
0 0
2 2 2nkc ncf
Lωπ π
= = =
( )1( , ) 2 cos cosn xp x t tLπ ω =
A
( ),10 0
2( , ) sin sinxn xu x t t
c Lπ ω
ρ =
A
(5.63)
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
57
Course Outline
A) Introduction and Outlook
B) Flame Aerodynamics and Flashback
C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts
D) Disturbance Propagation and Generation in Reacting Flows
E) Flame Response to Harmonic Excitation
• Introduction
• Decomposition of Disturbances into Fundamental Disturbance Modes
• Disturbance Energy
• Nonlinear Behavior
• Acoustic Wave Propagation Primer
• Unsteady Heat Release Effects and Thermoacoustic Instability
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
58
Unsteady Heat Release Effects
• Oscillations in heat release generate acoustic waves. – For unconfined flames, this is manifested as broadband noise emitted
by turbulent flames. – For confined flames, these oscillations generally manifest themselves
as discrete tones at the natural acoustic modes of the system.
• The fundamental mechanism for sound generation is the unsteady gas expansion as the mixture reacts.
• To illustrate, consider the wave equation with unsteady heat release
(57)
• The two acoustic source terms describe sound wave production associated with unsteady gas expansion, due to either heat release (first term) or non isomolar combustion (second term).
22 21 10 1 02
1( 1)p q nc p p
t ntγ γ
∂ ∂− ∇ = − +∂∂
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
59
Unsteady Heat Release Effects – Confined Flames
• Unsteady heat release from an acoustically compact flame induces a jump in unsteady acoustic velocity, but no change in pressure:
(61)
(62)
• Unsteady heat release causes both amplification/damping and shifts in phase of sound waves traversing the flame zone. – The relative significance of these two effects depends upon the
relative phase of the unsteady pressure and heat release. • The first effect is typically the most important and can cause systems with
unsteady heat release to exhibit self-excited oscillations. • The second effect causes shifts in natural frequencies of the system.
1 1b ap p=
,1 ,1 10
1 1x b x au u Q
A pγγ−
− =
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
60
Combustor Stability Model Problem
• Assume rigid and pressure release boundary conditions at x=0 and x=L, respectively, and that the flame is acoustically compact and located at x=LF .
• Unsteady pressure and velocity in the two regions are given by:
Region I: (63)
Region II: (64) ( ) ( ) ( )( )( ) ( ) ( )( )
0, 0,
,
1,
II F II F
II F II F
ik x L ik x L i tII II II
ik x L ik x L i tII II II
II II
p x t e e e
u x t e e ec
ω
ω
ρ
− − − −
− − − −
= +
= −
A B
A B
( ) ( ) ( )( )( ) ( ) ( )( ),1
0 0
,
1,
I F I F
I F I F
ik x L ik x L i tI I I
ik x L ik x L i tx I I I
I I
p x t e e e
u x t e e ec
ω
ω
ρ
− − − −
− − − −
= +
= −
A B
A B
ux,1=0 p1=0
L
I II
Heat releasezone
a b
LF
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
61
Matching Conditions
• Applying the boundary/matching conditions leads to the following algebraic equations:
(Left BC)
(Right BC)
(Pressure Matching)
(Velocity Matching) 6)
0I F I Fik L ik LI Ie e− =A -B
I I II II+ = +A B A B
( ) ( ) ( ) 12
0, 0,
, , 1II F I FI I
Qu L t u L tc
γρ
+ −− = −
( ) ( ) 0II F II Fik L L ik L LII IIe e− − −+ =A B
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
62
Unsteady Heat Release Model
• Model for unsteady heat release is the most challenging aspect of combustion instability prediction
• We will used a “velocity coupled” flame response model
• Assumes that the unsteady heat release is proportional to the unsteady flow velocity, multiplied by the gain factor, n, and delayed in time by the time delay, τ.
• This time delay could originate from, for example, the convection time associated with a vortex that is excited by the sound waves.
( )2
0, 0,1 1 ,
1I I
F
cQ A u x L t
ρτ
γ= = −
− n
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
63
Heat Release Modeling Simplifications
• Solving the four boundary/matching conditions leads to:
• In order to obtain analytic solutions, we will next assume that the flame is located at the midpoint of the duct, i.e., L=2LF and that the temperature jump across the flame is negligible.
(72)
( )2cos sin 2ikL e kLωτ= n
( ) ( )( )
( ) ( ) ( )( )
0, 0,
0, 0,
cos cos0
1 sin sin
II III F II F
I I
iI F II F
ck L k L L
c
e k L k L Lωτ
ρρ
−
= − + − n
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
64
Approximate Solution for Small Gain Values
• We can expand the solution around the solution by looking at
a Taylor series in wavenumber k in the limit .
0(2 1)
2nk L π
=
−=n
0=n
( ) ( )0 20
0 0
exp1
2 sini
k k Ok L k L
ω τ==
= =
= − +
nn
n n
n n
<< 1n
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
65
Complex Wavenumbers
• Wavenumber and frequency have an imaginary component,. Considering the time component and expanding ω:
(75)
• Real component is the frequency of combustor response.
• Imaginary component corresponds to exponential growth or decay in time and space. Thus, ωi>0 corresponds to a linearly unstable situation, referred to as combustion instability.
(77)
exp( ) exp( )exp( )r ii t i t tω ω ω− = −
00
cos( )1 ( 1)(2 1)
nr n
ω τω ωπ=
=
= + − −
nn
n
0 0sin( )( 1)2
ni
cLω τω == − nn
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
66
Rayleigh’s Criterion
• Rayleigh's criterion: Energy is added to the acoustic field when the
product of the time averaged unsteady pressure and heat release is
greater than zero.
( ) ( ) ( ) ( )1 1 0 0
2 12, ( ) sin cos sin ( )
2 2n
p x L t Q t t tπ
ω ω τ= =
− = ∝ −
n nn
( ) ( )11 1 0( / 2, ) ( ) 0 1 sin 0np x L t Q t ω τ+
== > ⇒ − >
n
( ) ( )1 1 0
2 1( / 2, ) ( ) sin sin
4 2n
p x L t Q tπ
ω τ=
− = ∝
nn
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
67
Stability of Combustor Modes
• Unstable 1/4 wave mode (n=1)
• Unstable 3/4 wave mode (n=2)
• Sign of alternates with time delay
– Important implications on why instability prediction is so difficult- no monotonic dependence upon underlying parameters
S U S1/4 wave
mode
1/2 1 3/2 0
3/4 wavemode
Frequencyωn=0
S S S S UUUUU
1/4
τT
(2 1)n π−n
1/4
1/ 2m mτ− < <
T
1/4
1/ 23 3m mτ +< <
T
1 1p Q
• Largest frequency shifts occur at the values where oscillations are not amplified and that the center of instability bands coincides with points of no frequency shift.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
68
Thermo-acoustic Instability Trends
• Two parameters, the heat release time delay, τ, and acoustic period, T, control instability conditions.
• Data clearly illustrate the non-monotonic variation of instability amplitude with .
Air
Fuel
LFIShorter flame
LFL
Longerflame
Flame stabilizedat centerbody
Air
FuelLSO
Flame stabilizeddownstream
0
0.1
0.2
0.3
0.4
0.5
0.6
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
τ /T
( )ppsi
′
Data illustrating variation of instability amplitude with normalized time delay. Image courtesy of D. Santavicca
τ T
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
69
Understanding Time Delays
• Data from variable length combustor
Measured instability amplitude (in psi) of
combustor as a function of fuel/air ratio and combustor length. Data courtesy of D.
Santavicca.
0.60 0.65 0.70 0.75
2075200019251850
25
30
35
40
45
50
Adiabatic flame temperature, K
Equivalence ratio
Com
bust
or le
ngth
, in
4
3
32
1
12
DecreasingPeriod, T
Decreasing time delay, τ
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
70
More Instability Trends 1 of 2
• Note non-monotonic variation of instability amplitude with axial injector location, due to the more fundamental variation of fuel convection time delay, τ.
• Biggest change in frequency is observed near the stability boundary.
Measured dependence of instability amplitude and frequency upon axial
location of fuel injector. Data obtained from Lovett and Uznanski.
0
1
2
3
4
5
6
6 8 10 12 14 16 18 20 22 24 26
Fuel injector location, LFI
p′
6 8 10 12 14 16 18 20 22 24 26
Fuel injector location, LFI
0
30
60
90
120
150
180
Freq
uenc
y, H
z
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
71
More Instability Trends 2 of 2
Measured dependence of the excited instability mode amplitude upon the mean velocity in the combustor inlet.
Obtained from measurements by the author.
0 500 1000Frequency (Hz)
=18m/s
p′
u
Frequency (Hz)
=25m/s
0 500 1000
u
0 500 1000Frequency (Hz)
=36m/su
0
0.02
0.04
0.06
0.08
15 20 25 30 35 40Premixer Velocity (m/s)
430 Hz630 Hz
p′
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
72 Normalized time, f0t
pp′
2500 50000-0.015
-0.01
-0.005
0.0
0.005
0.01
0.015
Limit Cycles and Nonlinear Behavior
• As amplitudes grow nonlinear effects grow in significance and the system is attracted to a new orbit in phase space, typically a limit cycle.
• This limit cycle oscillation can consist of relatively simple oscillations at some nearly constant amplitude, but in real combustors the amplitude more commonly "breathes" up and down in a somewhat random or quasi-periodic fashion.
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
73
Combustion Process and Gas Dynamic Nonlinearities
• Gas dynamic nonlinearities introduced by nonlinearities present in Navier-Stokes equations
• Combustion process nonlinearities are introduced by the nonlinear dependence of the heat release oscillations upon the acoustic disturbance amplitude.
Dependence of unsteady heat release magnitude and phase upon
velocity disturbance amplitude. Graph generated from data obtained by Bellows et al.
0
30
60
90
120
150
0
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 1
QQ′
x xu u′
Phas
e, d
egre
es
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited
74
Boundary Induced Nonlinearities
• The nonlinearities in processes that occur at or near the combustor boundaries also affect the combustor dynamics as they are introduced into the analysis of the problem through nonlinear boundary conditions.
• Such nonlinearities are caused by, e.g., flow separation at sharp edges or rapid expansions, which cause stagnation pressure losses and a corresponding transfer of acoustic energy into vorticity.
• These nonlinearities become significant when .
• Also, wave reflection and transmission processes through choked and unchoked nozzles become amplitude dependent at large amplitudes.
( )1u u′ Ο
Copyright © T. Lieuwen 2013 Unauthorized Reproduction Prohibited