large eddy simulation of azimuthal instability in a model ...zhi chen | ukctrf annual meeting,...
TRANSCRIPT
![Page 1: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/1.jpg)
Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019
Large Eddy Simulation of Azimuthal Instability in a Model Annular Gas Turbine Combustor
Zhi X. Chen & N. SwaminathanUniversity of Cambridge
Department of Engineering
longitudinal mode transition azimuthal mode
![Page 2: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/2.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �2
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 3: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/3.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �3
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 4: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/4.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �4
Thermoacoustic instabilities
• Fuel lean combustion is prone to thermoacoustic instabilities in gas turbines.
1. W. Polifke, 2nd Colloquium on Combustion Dynamics and Combustion Noise Menaggio, Sept. 20-22, 20162. Altunlu et al., J. Eng. Gas Turbines Power 136 (2014) 051510.3. N. Noiray, https://caps.ethz.ch/research/thermoacoustics.html.
Image sources:
1
1
![Page 5: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/5.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �4
Thermoacoustic instabilities
• Fuel lean combustion is prone to thermoacoustic instabilities in gas turbines.
1. W. Polifke, 2nd Colloquium on Combustion Dynamics and Combustion Noise Menaggio, Sept. 20-22, 20162. Altunlu et al., J. Eng. Gas Turbines Power 136 (2014) 051510.3. N. Noiray, https://caps.ethz.ch/research/thermoacoustics.html.
Image sources:
1
1
2
![Page 6: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/6.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �4
Thermoacoustic instabilities
• Fuel lean combustion is prone to thermoacoustic instabilities in gas turbines.
1. W. Polifke, 2nd Colloquium on Combustion Dynamics and Combustion Noise Menaggio, Sept. 20-22, 20162. Altunlu et al., J. Eng. Gas Turbines Power 136 (2014) 051510.3. N. Noiray, https://caps.ethz.ch/research/thermoacoustics.html.
Image sources:
1
1
2
3
![Page 7: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/7.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �5
Azimuthal mode
• Azimuthal mode instabilities are more complex, damaging and not yet well understood.
1. WSJ research; Rolls Royce2. M. Bauerheim et al., Phys. Fluids 28 (2016) 021303.
Image sources:
1
![Page 8: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/8.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �5
Azimuthal mode
• Azimuthal mode instabilities are more complex, damaging and not yet well understood.
1. WSJ research; Rolls Royce2. M. Bauerheim et al., Phys. Fluids 28 (2016) 021303.
Image sources:
1
![Page 9: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/9.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �5
Azimuthal mode
• Azimuthal mode instabilities are more complex, damaging and not yet well understood.
1. WSJ research; Rolls Royce2. M. Bauerheim et al., Phys. Fluids 28 (2016) 021303.
Image sources:
1
2
Smaller acoustic length compared to longitudinal modes
![Page 10: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/10.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �5
Azimuthal mode
• Azimuthal mode instabilities are more complex, damaging and not yet well understood.
1. WSJ research; Rolls Royce2. M. Bauerheim et al., Phys. Fluids 28 (2016) 021303.
Image sources:
1
2
Smaller acoustic length compared to longitudinal modes
Higher frequency and amplitude
![Page 11: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/11.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �5
Azimuthal mode
• Azimuthal mode instabilities are more complex, damaging and not yet well understood.
1. WSJ research; Rolls Royce2. M. Bauerheim et al., Phys. Fluids 28 (2016) 021303.
Image sources:
1
2
Smaller acoustic length compared to longitudinal modes
Higher frequency and amplitude
LES Onset of azimuthal instabilities
![Page 12: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/12.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �6
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 13: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/13.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �7
LES and combustion modelling
• 3D filtered N-S equations (mass, momentum and energy) are solved
• Compressibility effects are considered: pressure work heating, equation of state
1. Z.X. Chen et al., Combust. Flame 203 (2019) 279.2. Z.X. Chen et al., Proc. Combust. Inst. 37 (2019) 2325.
Refs:
Premixed
Non-premixed
![Page 14: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/14.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �8
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 15: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/15.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �9
Cambridge/NTNU gas turbine model annular combustor
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
![Page 16: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/16.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �9
Cambridge/NTNU gas turbine model annular combustor
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
![Page 17: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/17.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
![Page 18: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/18.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
![Page 19: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/19.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
![Page 20: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/20.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
U = A Ubsin(2πft)
Ub = 15 m /s
f = 1800 Hz(1st azim. mode)
![Page 21: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/21.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
U = A Ubsin(2πft)
Ub = 15 m /s
f = 1800 Hz(1st azim. mode)
FTF =·q′�/ ·q
U′�0/U0
= G(ω) eiΦω
![Page 22: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/22.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
U = A Ubsin(2πft)
Ub = 15 m /s
f = 1800 Hz(1st azim. mode)
FTF =·q′�/ ·q
U′�0/U0
= G(ω) eiΦω
![Page 23: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/23.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �10
Single sector – grid sensitivity and flame transfer function
1. Worth et al., Combust. Flame 160 (2013) 2476–24892. Mazur et al., Proc. Combust. Inst. 37 (2019) 5129–5136.
Refs:
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
U = A Ubsin(2πft)
Ub = 15 m /s
f = 1800 Hz(1st azim. mode)
FTF =·q′�/ ·q
U′�0/U0
= G(ω) eiΦω
·q = ∫VHRR(x, t) dV
![Page 24: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/24.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �11
Single sector – grid sensitivity for temperature and velocity fields
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Unforced
![Page 25: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/25.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �11
Single sector – grid sensitivity for temperature and velocity fields
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Unforced
x /D = 21
0.5
![Page 26: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/26.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �11
Single sector – grid sensitivity for temperature and velocity fields
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Unforced
x /D = 21
0.5
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
0
500
1000
1500
2000
2500
<! T>
[K]
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
-10
0
10
20
30
40
<! U>
[m/s]
1M2M
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
-10
0
10
20
30
40
<! U>
[m/s]
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
-10
0
10
20
30
40<
! U>
[m/s]
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
0
500
1000
1500
2000
2500
<! T>
[K]
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2r/D
0
500
1000
1500
2000
2500
<! T>
[K]
Mean axial velocity Mean temperature
![Page 27: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/27.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �12
Single sector – grid sensitivity for flame transfer function
U = A Ubsin(2π ft)
·q = ∫VHRR(x, t) dV
FTF =·q′�/ ·q
U′�0 /U0
= G(ω) eiΦω
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
![Page 28: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/28.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �12
Single sector – grid sensitivity for flame transfer function
U = A Ubsin(2π ft)
·q = ∫VHRR(x, t) dV
FTF =·q′�/ ·q
U′�0 /U0
= G(ω) eiΦω
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Two excitation amplitudes are tested:5% & 10%
p′�/p < 1 %
![Page 29: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/29.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �12
Single sector – grid sensitivity for flame transfer function
U = A Ubsin(2π ft)
·q = ∫VHRR(x, t) dV
FTF =·q′�/ ·q
U′�0 /U0
= G(ω) eiΦω
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Two excitation amplitudes are tested:5% & 10%
p′�/p < 1 %
![Page 30: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/30.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �12
Single sector – grid sensitivity for flame transfer function
U = A Ubsin(2π ft)
·q = ∫VHRR(x, t) dV
FTF =·q′�/ ·q
U′�0 /U0
= G(ω) eiΦω
Grid G1 G2
2 1.2
No. Cells 1.4M 2.1M
Δx /δ0th
Two excitation amplitudes are tested:5% & 10%
p′�/p < 1 %
![Page 31: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/31.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �13
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 32: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/32.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �14
Computed pressure signals in the plenum tubes
longitudinal mode transition azimuthal mode
![Page 33: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/33.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �14
Computed pressure signals in the plenum tubes
longitudinal mode transition azimuthal mode
![Page 34: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/34.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �14
Computed pressure signals in the plenum tubes
longitudinal mode transition azimuthal mode
![Page 35: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/35.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �14
Computed pressure signals in the plenum tubes
longitudinal mode transition azimuthal mode
Computational cost: 48 hrs on wall-clock for 0.1s using 1248 cores, 20M cells —> 900 kAUs
![Page 36: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/36.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �15
Longitudinal mode
![Page 37: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/37.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �16
Azimuthal mode
![Page 38: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/38.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �17
Longitudinal -> Azimuthal mode transition in frequency domain
![Page 39: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/39.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �18
Heat release rate fluctuations and acoustic mode structure
![Page 40: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/40.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �19
Outline
• Background – thermoacoustic instabilities & azimuthal mode
• LES formulation and combustion modelling
• Simulation of a model annular combustor
‣ Single sector – grid sensitivity and flame transfer function (ext. forcing)
‣ 12 burner full annular – self-excited oscillations
• Summary and future work
![Page 41: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/41.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �20
• Azimuthal instability in a model annular combustor
is studied using LES
• LES are performed in two steps:
‣ Single sector – grid sensitivity and flame transfer function (external forcing)
‣ 12 burner full annular – self-excited oscillations:
transition from longitudinal to azimuthal mode
• Projected future works:
‣ Analyse the LES data to study mode switching
‣ Non-adiabatic simulations
Summary and future work
![Page 42: Large Eddy Simulation of Azimuthal Instability in a Model ...Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 2 /21 Outline • Background – thermoacoustic](https://reader036.vdocument.in/reader036/viewer/2022081404/5f0609c77e708231d415fba0/html5/thumbnails/42.jpg)
/21Zhi Chen | UKCTRF Annual Meeting, Imperial College, London | 11.09.2019 �21
Thank you for listening!
Questions?
Acknowledgements