Burt Rutan's White Knight and SpaceShip One, Photo Courtesy of Scaled Composites
Thermodynamics and Propulsion
Subsections


6.9 Effect of Departures from Ideal Behavior -- Real Cycle behavior

To conclude this chapter, we will now improve our estimates of cycle performance by including the effects of irreversibility. We will use the Brayton cycle as an example. What are the sources of non-ideal performance and departures from reversibility?

  • Losses (entropy production) in the compressor and the turbine.
  • Stagnation pressure decrease in the combustor.
  • Heat transfer.

We take into account here only irreversibility in the compressor and in the turbine. Because of these irreversibilities, we need more work, $ \Delta h_\textrm{comp}$ (the changes in kinetic energy from inlet to exit of the compressor are neglected), to drive the compressor than in the ideal situation. We also get less work, $ \Delta h_\textrm{turb}$ , back from the turbine. The consequence, as can be inferred from Figure 6.16, is that the net work from the engine is less than in the cycle with ideal components.

Figure 6.16: Gas turbine engine (Brayton) cycle showing effect of departure from ideal behavior in compressor and turbine
Image fig5BraytonRealCycle_web

To develop a quantitative description of the effect of these departures from reversible behavior, consider a perfect gas with constant specific heats and neglect kinetic energy at the inlet and exit of the turbine and compressor. We define the turbine adiabatic efficiency as

$\displaystyle \eta_\textrm{turb} = \frac{w_\textrm{turb}^\textrm{actual}}{w_\textrm{turb}^\textrm{ideal}} = \frac{h_4-h_5}{h_4-h_{5s}},$    

where $ w_\textrm{turb}^\textrm{actual}$ is specified to be at the same pressure ratio as $ w_\textrm{turb}^\textrm{ideal}$ . There is a similar metric for the compressor, the compressor adiabatic efficiency:

$\displaystyle \eta_\textrm{comp} = \frac{w_\textrm{comp}^\textrm{ideal}}{w_\textrm{comp}^\textrm{actual}} = \frac{h_{3s}-h_0}{h_3-h_0},$

again for the same pressure ratio. Note that for the turbine the ratio is the actual work delivered divided by the ideal work, whereas for the compressor the ratio is the ideal work needed divided by the actual work required. These are not thermal efficiencies, but rather measures of the degree to which the compression and expansion approach the ideal processes.

We now wish to find the net work done in the cycle and the efficiency. The net work is given either by the difference between the heat received and rejected or the work of the compressor and turbine, where the convention is that heat received is positive and heat rejected is negative and work done is positive and work absorbed is negative.

$\displaystyle \textrm{Net work}=\begin{cases}\underbrace{q_A}_\textrm{heat in} ...
...0) w_\textrm{turb} - w_\textrm{comp} &= (h_4 - h_5) - (h_3 -h_0). \end{cases}$    

The thermal efficiency is

$\displaystyle \eta_\textrm{thermal} = \frac{\textrm{Net work}}{\textrm{Heat input}}.$

We need to calculate $ T_3$ , $ T_5$ .

From the definition of $ \eta_\textrm{comp}$ ,

$\displaystyle T_3 - T_0 =\frac{(T_{3s} - T_0)}{\eta_\textrm{comp}}= T_0\frac{(
T_{3s}/T_0 - 1)}{\eta_\textrm{comp}}.$

With

$\displaystyle (T_{3s}/T_0)= \textrm{isentropic
temperature ratio} =
\left(\frac...
...{comp}^{\frac{\gamma-1}{\gamma}} =
\Pi_\textrm{comp}^{\frac{\gamma-1}{\gamma}},$

$\displaystyle T_3=T_0\frac{\left(\Pi_\textrm{comp}^{\frac{\gamma-1}{\gamma}}-1\right)}{\eta_\textrm{comp}}+T_0.$

Similarly, by the definition

$\displaystyle \eta_\textrm{turb} = \frac{\textrm{actual work received}}
{\textrm{ideal work for same
$\frac{P_\textrm{exit}}{P_\textrm{inlet}}$}},$

we can find $ T_5$ :

$\displaystyle T_4-T_5$ $\displaystyle = \eta_\textrm{turb}(T_4-T_{5s}) = \eta_\textrm{turb}T_4\left(1-\...
... =\eta_\textrm{turb}T_4\left(1-\Pi_\textrm{turb}^\frac{\gamma-1}{\gamma}\right)$    
$\displaystyle T_5$ $\displaystyle =T_4-\eta_\textrm{turb}T_4\left(1-\Pi_\textrm{turb}^\frac{\gamma-1}{\gamma}\right).$    

The thermal efficiency can now be found:

$\displaystyle \eta_\textrm{thermal} = 1 + \frac{Q_L}{Q_H} = 1 - \frac{T_5-T_0}{T_4-T_3}.$

With

$\displaystyle \Pi_\textrm{comp}$ $\displaystyle = \frac{1}{\Pi_\textrm{turb}} = \Pi$    

and


$\displaystyle \tau_s$ $\displaystyle = \Pi^{\frac{\gamma-1}{\gamma}} = \textrm{the isentropic cycle temperature ratio},$    

$\displaystyle \eta_\textrm{thermal}$ $\displaystyle = 1-\cfrac{T_4\left[1-\eta_\textrm{turb}\left(1-\cfrac{1}{\tau_s}...
...ht)\right]-T_0} {T_4-T_0\left[\cfrac{1}{\eta_\textrm{comp}}(\tau_s-1)+1\right]}$    

or


$\displaystyle \eta_\textrm{thermal}$ $\displaystyle = \frac{\left[1-\cfrac{1}{\tau_s}\right]\left[\eta_\textrm{comp}\...
...0}-\tau_s\right]} {1+\eta_\textrm{comp}\left[\cfrac{T_4}{T_0}-1\right]-\tau_s}.$    

There are several non-dimensional parameters that appear in this expression for thermal efficiency. We list these in the two sections below and show their effects in accompanying figures.

6.9.1 Parameters reflecting design choices

$ (\Pi)$ : cycle pressure ratio    
$\displaystyle \frac{T_4}{T_0}$ $\displaystyle : \textrm{maximum turbine inlet temperature}$    

6.9.2 Parameters reflecting the ability to design and execute efficient components

$\displaystyle \eta_\textrm{comp}$ $\displaystyle :\textrm{compressor adiabatic efficiency}$    
$\displaystyle \eta_\textrm{turb}$ $\displaystyle :\textrm{turbine adiabatic efficiency}$    

In addition to efficiency, net rate of work is a quantity we need to examine,

$\displaystyle \dot{W}_\textrm{net} =\dot{W}_\textrm{turbine} -\dot{W}_\textrm{compressor}.$

Putting this in a non-dimensional form,

$\displaystyle \frac{\dot{W}_\textrm{net}}{\dot{m} c_p T_0} =
\underbrace{-\frac...
...\left(1-\frac{1}{\tau_s}\right)}
_\textrm{work extracted from flow by turbine}.$

$\displaystyle \frac{\dot{W}_\textrm{net}}{\dot{m}c_pT_0} =
(\tau_s-1)\left[\cfrac{\eta_\textrm{turb}\cfrac{T_4}{T_0}}{\tau_s}-\frac{1}{\eta_\textrm{comp}}\right]$

Figure 6.17: Non-dimensional power and efficiency for a non-ideal gas turbine engine [from Cumpsty, Jet Propulsion]
[Non-dimensional work as a function of cycle pressure ratio for different values of turbine entry temperature divided by compressor entry temperature] Image fig5GasTurbPowerEfficiencyWorkVSOPR_web [Overall cycle efficiency as a function of pressure ratio for different values of turbine entry temperature divided by compressor entry temperature] Image fig5GasTurbPowerEfficiencyEtaVSOPRTempRatios_web [Overall cycle efficiency as a function of cycle pressure ratio for different component efficiencies] Image fig5GasTurbPowerEfficiencyEtaVSOPRAdEff_web

Trends in net power and efficiency are shown in Figure 6.17 for parameters typical of advanced civil engines. Some points to note in the figure:

  • For any $ \eta_\textrm{comp}$ , $ \eta_\textrm{turb} \neq 1$ the optimum pressure ratio $ (\Pi)$ for maximum $ \eta_\textrm{th}$ is not the highest that can be achieved, as it is for the ideal Brayton cycle. The ideal analysis is too idealized in this regard. The highest efficiency also occurs closer to the pressure ratio for maximum power than in the case of an ideal cycle. Choosing this as a design criterion will therefore not lead to the efficiency penalty inferred from ideal cycle analysis.
  • There is a strong sensitivity to the component efficiencies. For example, for $ \eta_\textrm{turb}=\eta_\textrm{comp}= 0.85$ , the cycle efficiency is roughly two-thirds of the ideal value.
  • The maximum power occurs at a value of $ \tau_s$ or pressure ratio $ (\Pi)$ less than that for max $ \eta$ (this trend is captured by ideal analysis).
  • The maximum power and maximum $ \eta_\textrm{thermal}$ are strongly dependent on the maximum temperature, $ T_4/T_0$ .

A scale diagram of a Brayton cycle with non-ideal compressor and turbine behaviors, in terms of temperature-entropy ($ h$ -$ s$ ) and pressure-volume ($ P$ -$ v$ ) coordinates is given below as Figure 6.18.

Figure 6.18: Scale diagram of non-ideal gas turbine cycle. Nomenclature is shown in the figure. Pressure ratio $ 40$ , $ T_0 = 288\textrm{ K}$ , $ T_4 =1700\textrm{ K}$ , compressor and turbine efficiencies $ = 0.9$ [from Cumpsty, Jet Propulsion]
[] Image fig5BraytonScaleNonIdealTS_web [] Image fig5BraytonScaleNonIdealPV_web





Muddy Points

Isn't it possible for the mixing of two gases to go from the final state to the initial state? If you have two gases in a box, they should eventually separate by density, right? (MP 6.22)

How can $ T_4/T_0$ be the maximum turbine inlet temperature? (MP 6.23)

When there are losses in the turbine that shift the expansion in $ T$ -$ s$ diagram to the right, does this mean there is more work than ideal since the area is greater? (MP 6.24)

For an afterburning engine, why must the nozzle throat area increase if the temperature of the fluid is increased? (MP 6.25)

Why doesn't the pressure in the afterburner go up if heat is added? (MP 6.26)

Why is the flow in the nozzle choked? (MP 6.27)

What's the point of having a throat if it creates a retarding force? (MP 6.28)

Why isn't the stagnation temperature conserved in this steady flow? (MP 6.29)

Douglas Quattrochi 2006-08-06