Flame Dynamics, Hydrodynamics, and Acoustics: Tim Lieuwen
Flame Dynamics, Hydrodynamics, and Acoustics: Tim Lieuwen
Tim Lieuwen
Georgia Institute of Technology
1
School of Aerospace Engineering
• Fluid Mechanics
– Hydrodynamic stability
of inhomogeneous flows
• Acoustics
– Wave propagation in
complex flows
• Combustion
– Response of reaction
fronts to disturbances
2
School of Aerospace Engineering
Combustion Instability
Fourier Transform
combustor 300
Frequency (Hz)
3
School of Aerospace Engineering
4
School of Aerospace Engineering
5
School of Aerospace Engineering
• Key mechanisms:
– Kelvin-Helmholtz Instability
– Rayleigh-Taylor Instability
6
School of Aerospace Engineering
• Dispersive
7
School of Aerospace Engineering
• Dispersive
8
School of Aerospace Engineering
Absolutely unstable
Emerson, B., O’Connor, J., Juniper, M. and Lieuwen, T., JFM, 2012.
10
School of Aerospace Engineering
• Dispersive
0.5
0
-4 -2 0 2 4 6
12
School of Aerospace Engineering
0.5
0
-4 -2 0 2 4 6
13
School of Aerospace Engineering
14
School of Aerospace Engineering
Q '(ω) / Qo
0.4 G=
u '(ω) / uo
Why does this saturate?
Q’/Qoo
0.3
CH*′ / CH*
15
School of Aerospace Engineering
Increased
Amplitude of Forcing
16
School of Aerospace Engineering
17
School of Aerospace Engineering
18
School of Aerospace Engineering
Power
Spectrum
L’(x, f0)
19
School of Aerospace Engineering
Convective wavelength:
λc= U0/f0at forcing
• Strong response
frequency
- ‒distance a disturbance
Non-monotonic spatial
propagates at mean flow
dependence
speed in one excitation
period
20
School of Aerospace Engineering
∂L ∂ G
: ∂L+ u ⋅∇ ∂L |
2 F. Williams,
G-equation =
G s | ∇
G
+ u − v = S L L1 + Combustion
∂t ∂t ∂x ∂x Theory
22
School of Aerospace Engineering
23
School of Aerospace Engineering
24
School of Aerospace Engineering
∂L′ ( x, t ) 1 x ∂un′ x − x′ 1 x
∂x
=
ut ∫ 0 ∂x
( x′, t −
ut
) dx′ + ⋅ un′ ( x =0, t =−
ut
t
ut
)
25
School of Aerospace Engineering
u0 uc,f
ut ut
sL
27
School of Aerospace Engineering
28
School of Aerospace Engineering
29
School of Aerospace Engineering
un′
= ε n cos(2π f (t − x / uc ,v ))
ut
• Linearized solution:
L' − i ⋅ ε n sin θ i 2π f ( y / ( uc ,v tan θ )− t ) i 2π f ( y / ( ut sin θ )− t )
= Real × e −e
ut f 2π ut cos θ / u c ,v − 1 ( )
30
School of Aerospace Engineering
Solution Characteristics
• Note interference pattern on
flame wrinkling λint
1
λint ( λ t sin θ ) =
| ut uc ,v − 1|
L’
y / ( λ t sin θ )
31
School of Aerospace Engineering
L’
L’
1
λint ( λ t sin θ ) = Acharya et al., ASME Turbo Expo, 2011
| ut uc ,v − 1|
• High phase acoustic speed disturbances (with long wavelengths)
where uc,v>>ut lead to short length scale disturbances
• Vortical disturbances where ut~uc,v lead to long length scale
32 wrinkling
School of Aerospace Engineering
Lre
disturbances key to
f
L’2
interference patterns Single frequency
excitation
• Example: convecting or
random disturbances to
simulate turbulent flow
disturbances
33
School of Aerospace Engineering
34
School of Aerospace Engineering
36
School of Aerospace Engineering
Q '(ω) / Qo
0.4 G=
u '(ω) / uo
Why does this saturate?
Q’/Qoo
0.3
CH*′ / CH*
37
School of Aerospace Engineering
38
School of Aerospace Engineering
40
School of Aerospace Engineering
• Combustion
instabilities emerged
during full engine
testing of GenX engine
41
School of Aerospace Engineering
42
School of Aerospace Engineering
Financial Times
23 July 99, Issue 49
Siemens, together with Italy’s Ansaldo, took the turnkey contract for the 350 MW plant in 1996 and
should have had it in service by May of last year. The startup was delayed till January. Since then
matters have worsened. There have been two major breakdowns and, says Colbun, there have been no
satisfactory explanations.
The trouble could not have come worse for Colbun. The manly hydroelectric generator, which is
controlled by a consortium made up of Belgium’s Tractebel, Spain’s Iberdrola and the local Matte and
Yaconi-Santa Cruz groups, has been crippled by severe drought in Chile, which has slashed its output
and thrown it back – without Nehuenco – onto a prohibitively expensive spot market.”
43
School of Aerospace Engineering
Wrinkle Convection
Flamelet Step Response problem: ua t<0
u=
Step increase axial velocity from ua to ub, : ub t≥0
? .
44
School of Aerospace Engineering
Wrinkle Convection
ua t<0
u=
ub t≥0
45
School of Aerospace Engineering
46
School of Aerospace Engineering
θ
u0 u
t
sL
47
School of Aerospace Engineering
48
School of Aerospace Engineering
x/λ c
Numerical Calculation Experimental Result
50