14.09.2014 Views

Impact of fuel supply impedance and fuel staging on gas turbine ...

Impact of fuel supply impedance and fuel staging on gas turbine ...

Impact of fuel supply impedance and fuel staging on gas turbine ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Numerical simulati<strong>on</strong><br />

can be calculated for a single jet injected into a circular duct. J represents the<br />

impulse ratio <str<strong>on</strong>g>of</str<strong>on</strong>g> the <str<strong>on</strong>g>fuel</str<strong>on</strong>g> <str<strong>on</strong>g>and</str<strong>on</strong>g> main stream <str<strong>on</strong>g>and</str<strong>on</strong>g> is defined as:<br />

J = ρ F u 2 F<br />

ρ A u 2 . (6.3)<br />

A<br />

The subscripts F <str<strong>on</strong>g>and</str<strong>on</strong>g> A represent again the <str<strong>on</strong>g>fuel</str<strong>on</strong>g> <str<strong>on</strong>g>and</str<strong>on</strong>g> main stream, respectively.<br />

In c<strong>on</strong>figurati<strong>on</strong> A, for example, the main flow velocity is equal to u A = 18.6<br />

m/s. The density has a value <str<strong>on</strong>g>of</str<strong>on</strong>g> ρ A = 1.167 kg/m 3 . With the properties for the<br />

<str<strong>on</strong>g>fuel</str<strong>on</strong>g> stream (u F = 85.6 m/s, ρ F = 0.66 kg/m 3 ) an impulse ratio <str<strong>on</strong>g>of</str<strong>on</strong>g> J = 12.0 is<br />

computed for the present case. The <str<strong>on</strong>g>fuel</str<strong>on</strong>g> injecti<strong>on</strong> including stream lines <str<strong>on</strong>g>and</str<strong>on</strong>g><br />

the <str<strong>on</strong>g>fuel</str<strong>on</strong>g> mass fracti<strong>on</strong> at a downstream locati<strong>on</strong> are presented exemplarily<br />

for the c<strong>on</strong>figurati<strong>on</strong> A in Fig. 6.4. The legend refers to the plane presented<br />

whereas the colors <str<strong>on</strong>g>of</str<strong>on</strong>g> the streamlines indicate the amplitude <str<strong>on</strong>g>of</str<strong>on</strong>g> the velocity.<br />

The <str<strong>on</strong>g>fuel</str<strong>on</strong>g> is distributed well over the radial directi<strong>on</strong>, with the main part <str<strong>on</strong>g>of</str<strong>on</strong>g> the<br />

<str<strong>on</strong>g>fuel</str<strong>on</strong>g> located in the middle between the outer wall <str<strong>on</strong>g>and</str<strong>on</strong>g> the center body. The<br />

maximum depth according to Eqn. (6.2) would be 0.005 m, which is slightly<br />

too low. The optimum number <str<strong>on</strong>g>of</str<strong>on</strong>g> circumferential holes was estimated according<br />

to Holdeman et al. [45] <str<strong>on</strong>g>and</str<strong>on</strong>g> lies in the range between 6.1 <str<strong>on</strong>g>and</str<strong>on</strong>g> 10.2 for the<br />

c<strong>on</strong>figurati<strong>on</strong> analyzed (see appendix A.2 for more details). CFD simulati<strong>on</strong>s<br />

with different numbers <str<strong>on</strong>g>of</str<strong>on</strong>g> holes showed that the c<strong>on</strong>figurati<strong>on</strong> with eight injecti<strong>on</strong><br />

holes exhibits adequate mixing characteristics.<br />

6.1.2 Numerical setup <str<strong>on</strong>g>of</str<strong>on</strong>g> the combusti<strong>on</strong> system<br />

To determine the flame transfer functi<strong>on</strong>s, a three-dimensi<strong>on</strong>al unsteady<br />

RANS computati<strong>on</strong> <str<strong>on</strong>g>of</str<strong>on</strong>g> the compressible turbulent reacting flow in the combustor<br />

has been performed using the s<str<strong>on</strong>g>of</str<strong>on</strong>g>tware package ANSYS-CFX. According<br />

to the discussi<strong>on</strong> in chapter 2.3.2, the computati<strong>on</strong>al domain can be reduced<br />

to the mixing secti<strong>on</strong>, the <str<strong>on</strong>g>fuel</str<strong>on</strong>g> injecti<strong>on</strong> holes, the axial swirler <str<strong>on</strong>g>and</str<strong>on</strong>g> the<br />

combusti<strong>on</strong> chamber. A sketch <str<strong>on</strong>g>of</str<strong>on</strong>g> the computati<strong>on</strong>al domain including the<br />

”measurement” planes for case A, at which the time series <str<strong>on</strong>g>of</str<strong>on</strong>g> the required velocity<br />

fluctuati<strong>on</strong>s are recorded, is shown in Fig. 6.5. To save computati<strong>on</strong>al<br />

time <strong>on</strong>ly a 90 ◦ segment is simulated using periodic boundary c<strong>on</strong>diti<strong>on</strong>s.<br />

A sec<strong>on</strong>d order backward Euler scheme <str<strong>on</strong>g>and</str<strong>on</strong>g> the high-resoluti<strong>on</strong> scheme were<br />

110

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!