Thermodynamics of Gas Turbine Cycles
Thermodynamics of Gas Turbine Cycles
1.1 Introduction
1
1.2. Thermodynamic Systems
1.2.1. System, Boundary and Surroundings
System. A system is a finite quantity of matter or a prescribed region of
space, see figure (1-1).
2
Fig.1-3: Closed system . Fig.1-4: Open system.
3
1.4 Pure Substance
A pure substance is one that has a homogeneous and invariable chemical
composition even though there is a change of phase. In other words, it is a
system which is (a) homogeneous in composition, (b) homogeneous in
chemical aggregation. Examples: Liquid, water, mixture of liquid water
and steam, mixture of ice and water. The mixture of liquid air and gaseous
air is not a pure substance.
4
1.6 State
State is the condition of the system at an instant of time as described or
measured by its properties. Or each unique condition of a system is called
a state.
It follows from the definition of state that each property has a single value
at each state.
Stated differently, all properties are state or point functions. Therefore, all
properties are identical for identical states.
— A variable is a property, if and only if, the change in its value between
any two pre-scribed equilibrium states is single-valued.
5
Note that the process path indicates a series of equilibrium states through
which the system passes during a process and has significance for quasi-
equilibrium processes only.
The Process Path defines the type of process undergone. Typical process
paths are:
6
internal energy which can be felt due to increase in pressure and
temperature.
1.8 Entropy
The ‘thermodynamic flow’ quantity to be considered in heat engines has
𝑄̇
the value of . This quantity is conserved in a reversible engine; and is
𝑇
𝛿𝑄
∮ ≤0 (1-1)
𝑇
𝛿𝑄
for reversible processes ∮ =0 (1-2)
𝑇
7
𝛿𝑄
While; for irreversible processes ∮ <0 (1-3)
𝑇
- Power cycles Thermodynamic power cycles are the basis for the
operation of heat engines, which supply most of the world's electric
power and run almost all motor vehicles. Power cycles can be
divided according to the type of heat engine they seek to model. The
most common cycles that model internal combustion engines are the
Otto cycle, which models gasoline engines and the Diesel cycle,
which models diesel engines. Cycles that model external combustion
engines include the Brayton cycle, which models gas turbines, and
the Rankine cycle, which models steam turbines.
- Heat pump cycles
Thermodynamic heat pump cycles are the models for heat pumps
and refrigerators. The difference between the two is that heat pumps
8
are intended to keep a place warm while refrigerators are designed
to cool it. The most common refrigeration cycle is the vapor
compression cycle, which models systems using refrigerants that
change phase. The absorption refrigeration cycle is an alternative
that absorbs the refrigerant in a liquid solution rather than
evaporating it.
Fig .1-6: Brayton cycle on P-V diagram Fig .1-7: Brayton cycle on T-φ diagram
The cold air at 3 is fed to the inlet of the compressor where it is compressed
along 3-4 and then fed to the combustion chamber where it is heated at
constant pressure along 4-1. The hot air enters the turbine at 1 and expands
adiabatically along 1-2 and is then cooled at constant pressure along 2-3.
9
Work done = Heat supplied – Heat rejected
=𝐾𝑝 (𝑇1 − 𝑇4 ) − 𝐾𝑝 (𝑇2 − 𝑇3 ) (1-7)
Thermal efficiency (η) of Brayton Cycle :
𝑊𝑜𝑟𝑘 𝑑𝑜𝑛𝑒 [𝐾𝑝 {(𝑇1 −𝑇4 )−𝐾𝑝 (𝑇2 −𝑇3 )}]
𝜂= =
𝐻𝑒𝑎𝑡 𝑆𝑢𝑝𝑝𝑙𝑖𝑒𝑑 𝐾𝑝 (𝑇1 −𝑇4 )
(𝑇 −𝑇3 )
𝜂 = 1 − (𝑇2 (1-8)
1 −𝑇4 )
(𝑇2 −𝑇3 )
𝜂 =1− (𝛾−1)
⁄𝛾
(1-11)
𝑃1
(𝑃 ) (𝑇2 −𝑇3 )
2
11
Fig.1-8: Rankine cycle in steam turbine
Fig.1-9: Rankine cycle Pressure Volume P-V Fig .1-10: Rankine cycle on S-T diagram
The Rankine cycle is the ideal cycle for vapor power plants; it includes the
following four reversible processes:
11
The thermal efficiency of the cycle is determined from:
(𝐻1 − 𝐻2 )
⁄(𝐻 − 𝐻 ) (1-12)
1 𝑊2
Where,
(𝑇1 −𝑇2 )
𝜂= (𝑇1 )
(1-13)
Where,
T2 = Temperature of receiver
12
1.10 Thermodynamic lows
“If two systems are both in thermal equilibrium with a third then they
are in thermal equilibrium with each other“.
13
1.10.2 First law of thermodynamics
The boundary and the interactions that are present at the boundary play
roles in the analysis devoted to solving a problem. One feature is whether
the boundary is crossed by the flow of mass. A system defined by a
boundary impermeable to mass flow is a closed system . Conversely,
systems whose defining boundaries can be crossed by the flow of mass are
open systems, or flow systems.
14
1.10.2`.1 First law of thermodynamics for closed system (Control Mass):
Energy interactions
(Non properties)
Fig.1-12: Graphic statements of the first law of thermodynamics for closed systems.
The difference between the net heat input Q1–2 and the net work output
W1–2 represents the change in the thermodynamic property called energy.
The first law proclaims the existence of energy as a thermodynamic
property. It defines the property “energy change.”
The energy change E2−E1 depends only on the end states, whereas the
energy interactions Q1–2 and W1–2 depend on the end states and the path of
15
the process that links the end states. This important distinction is stressed
under each term appearing in equation (1-14). Another way to stress this
difference is to use a different notation for the infinitesimal increments in
work and heat transfer relative to the exact differential notation that applies
to the infinitesimal change in E.
For this reason, the first law for a process between two states situated
infinitely close to one another is written as:
𝛿𝑄 − 𝛿𝑊 = 𝑑𝐸 (1-15)
In the same notation, the energy interactions on the left side of equation
(1-14) are:
2 2
𝑄1−2 = ∫1 𝛿𝑄 and 𝑊1−2 = ∫1 𝛿𝑊 (1-16)
The peculiar notation 𝛿 may not be the ideal way to emphasize the
difference between energy interactions and energy change. The alternative
is to introduce the concept of time in the description of the process, see the
figure (1-12) .
In this new description, state 1 is the condition of the system at time t1,
state 2 is the condition at time t2 , and the net energy interactions Ql–2 and
W1–2 are the time integrals :
𝑡 𝑡
𝑄1−2 = ∫𝑡 2 𝑄̇ 𝑑𝑡 , 𝑊1−2 = ∫𝑡 2 𝑊̇ 𝑑𝑡 (1-17)
1 1
Quantities 𝑄̇ and 𝑊̇ are the instantaneous heat transfer rate and the
mechanical power output, respectively. The first law of thermodynamics
for a closed system can be written on a per-unit-time basis as:
𝑑𝐸
𝑄̇ − 𝑤̇ = (1-18)
𝑑𝑡
16
Another way of stressing the path dependence of Q1–2 and W1–2 is
presented graphically in figure (1-13). The system can proceed from state
1 to state 2 along an infinity of paths.
Fig .1-13: Path dependence of the energy interactions Q1–2 and W1–2
Although the difference Q1–2 , W1–2 matches E2−E1 along both paths, the
sizes of Q1–2 and W1–2 vary from one path to the next. If the process
executed by the closed system is a cycle, the first law (1-14) reduces to:
∮ 𝛿𝑄 − ∮ 𝛿𝑊 = 0 (1-19)
The major difference between a Control Mass and Control Volume is that
mass crosses the system boundary of a control volume.
Figure (1-14) shows the main features of an open system: heat interactions
per unit time, 𝑄̇; work interactions per unit time 𝑊̇ ; and portions of the
boundary that are crossed by the flow of mass. For simplicity, the figure
shows only one of each type of boundary crossing, one inlet port labeled
“ in,” and one outlet port labeled “out.”
17
The open system, or the control volume, is the region contained between
the inlet and outlet ports; in other words, the stationary dashed lines labeled
in and out belong to the boundary of the open system. The work transfer
ratė 𝑊̇ refers to any mode or combination of work modes, P (dV/dt) ,̇ 𝑊𝑠ℎ
̇ ;
̇
𝑊𝑒𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑎𝑙 ̇
; 𝑊𝑚𝑎𝑔𝑛𝑒𝑡𝑖𝑐 , and so on .
Fig .1-14: Flow of a closed system (shaded area) through the space occupied by an
open system and the conversion of the first law for closed systems into a statement
valid for open systems.
For the process from state 1 (time t) to state 2 (time t + Δt) executed by the
18
The last two terms on the right side account for the P dV type of work
transfer associated with the deformation of the closed system from time t
to time t+Δt.
Note that in each term P is the local pressure: that is, the pressure in the
immediate vicinity of the port. Relations similar to eq. (1-20) express the
relative size of the energy inventories of the closed and open systems:
The ΔE’s and ΔV’s can be rewritten in terms of their per-unit mass as e
and v:
Like the port pressure P, the specific energy and volume (e and v,
respectively) are properties of the intensive state of the fluid that crosses
the boundary at time t.
Combining equations (1-21) and (1-23) for the purpose of eliminating the
terms that refer to the energy inventory of the closed system (Eclosed), we
obtain:
1 ∆𝑀 ∆𝑀
(𝐸𝑜𝑝𝑒𝑛,(𝑡+∆𝑡) − 𝐸𝑜𝑝𝑒𝑛 ,𝑡 = 𝑄̇ − 𝑊̇ + [(𝑒 + 𝑃𝑣 ) ] − [(𝑒 + 𝑃𝑣 ) ]
∆𝑡 ∆𝑡 𝑖𝑛 ∆𝑡 𝑜𝑢𝑡
(1-25)
Invoking the limit Δt →0, writing 𝑚̇ for the mass flow rate ΔM∕Δt ,
dropping the subscript “open” from the energy inventory of the control
volume, and assuming that more than one inlet port and outlet port exist,
we arrive at the most general statement of the first law of thermodynamics
for an open system :
𝑑𝐸
= 𝑄̇ − 𝑊̇ + ∑𝑖𝑛 𝑚̇ (𝑒 + 𝑃𝑣) − ∑𝑜𝑢𝑡 𝑚̇(𝑒 + 𝑃𝑣) (1-26)
𝑑𝑡
19
What makes this statement more general than the per-unit-time first law
for closed systems equation (1-18) are the terms 𝑚̇(𝑒 + 𝑃𝑣) : These terms
represent the energy transfer associated with the flow of mass across the
system boundary.
enthalpy:
ℎ = 𝑢 + 𝑃𝑣 (1-27)
Shows up explicitly in the terms accounting for energy transfer via mass
flow:
𝑑𝐸 1 1
= 𝑄̇ − 𝑊̇ + ∑𝑖𝑛 𝑚̇ (𝑢 + 𝑉 2 + 𝑔𝑧) − ∑𝑜𝑢𝑡 𝑚̇(𝑢 + 𝑉 2 + 𝑔𝑧) (1-28)
𝑑𝑡 2 2
In the fields of gas dynamics and compressible fluid mechanics, the group
1
(𝑢 + 𝑉 2 + 𝑔𝑧) is recognized as the local stagnation enthalpy of the
2
flowing fluid.
This equation spells out the difference between open systems and closed
systems (in the latter, the 𝑚̇ values are all zero and the mass inventory M
is a constant).
21
1.10.3 Second law of thermodynamics
21
1.10.3.1 Second law of thermodynamics for closed system
The essence of the second law is that in the Rumford and Joule experiments
the apparatuses received work and rejected heat. It was never the other way
around; in fact, all the attempts to construct a heat engine that would
operate cyclically as a closed system while in possible contact with a single
temperature reservoir have failed.
∮ 𝛿𝑊 ≤ 𝑜 (1-30)
or, in view of the first law for a cycle executed by a closed system equation
(1-19) :
∮ 𝛿𝑄 ≤ 𝑜 (1-31)
The next step in the direction of more complex cycles is made using the
concept of are eversible cycle, introduced by Sadi Carnot in 1824. This
concept became known as the Carnot cycle.
22
Fig.1-15: Closed system executing a reversible cycle while in communication with two
temperature reservoirs T1 and T2. No assumption is made regarding the relative
magnitude of temperaturesT1 andT2 and the sense of the cycle on the P−v and T−s
planes.
It is important that each state visited by the system during the cycle is one
of uniform pressure, temperature, and specific volume. This means that the
system expands and contracts at such a slow rate that at each point in time
the state of the system is one of equilibrium.
An isothermal volume change is “slow enough” when the system is in
thermal equilibrium with the temperature reservoir. Under such conditions,
the system can execute the same cycle in the reverse sense; that is, it can
pass through the same sequence of equilibrium states in reverse order.
Therefore, if WC is the net work produced by the Carnot cycle and Q1C and
Q2 Care the respective heat interactions with reservoirs T1 and T2 during
the same cycle, the energy interactions experienced by the reversed Carnot
cycle are:
( 𝑊𝐶 , 𝑄1𝐶 , 𝑄2𝐶 )𝑟𝑒𝑣𝑒𝑟𝑠𝑒𝑑 = − ( 𝑊𝐶 , 𝑄1𝐶 , 𝑄2𝐶 ) (1-32)
To extend the second law represented by eq. (1-30) to the cyclical operation
of closed systems in contact with two temperature reservoirs. Figure (1-16)
shows a closed system (A) executing an unspecified cycle while in contact
23
with reservoirs T1 and T2, with T1≠T2. The net energy interactions
experienced by this unspecified cycle, Q1,and Q2, are related through the
first law:
𝑄1 + 𝑄2 = 𝑊 (1-33)
Fig 1-16 closed system that executes a cycle while in communication with two
temperature reservoirs
𝑊>0 (1-34)
24
Consider now the arrangement in figure (1-17), the unspecified cycle (A)
and the Carnot cycle (C) of figure (1-16) share reservoirs T1 and T2.
𝑄1 + 𝑄1𝐶 = 0 (1-36)
Which means that the reservoir T1 also executes a cycle. Because the
composite system ( A+C+T1 ) completes a cycle while making contact
with only one reservoir (T2), the second law (1-31) requires :
𝑄2 + 𝑄2𝐶 ≤ 0 (1-37)
−𝑄2 𝑄2𝐶
≥ (1-38)
𝑄1 𝑄1𝐶
At this stage, equation (1-38) allows us to conclude that the ratio obtained
by dividing the absolute value of the negative heat inter-action by the value
of the positive heat interaction cannot be smaller than a certain value. The
limiting case corresponds to the equal sign in equation (1-38):
−𝑄2 𝑄2𝐶
= (1-39)
𝑄1 𝑄1𝐶
The equal sign in the second law (1-39) is associated with any reversible
cycle that is executed by the unspecified system (A). Therefore, the second
law can be stated by making reference only to the system of interest (A):
−𝑄2 −𝑄2
≥( ) (1-41)
𝑄1 𝑄1 𝑟𝑒𝑣
25
Where the subscript rev stands for “reversible”: in this case, for a reversible
cycle executed by an arbitrary closed system in contact with two
temperature reservoirs.
The limiting value (−Q2 ∕ Q1)rev can only be a function of parameters T1 and
T2, because the existence of two different temperature reservoirs was
assumed in the derivation of equation (1-41), Writing:
−𝑄2
( ) = 𝑓(𝑇1 , 𝑇2 ) (1-42)
𝑄1 𝑟𝑒𝑣
The unknown function f (T1,T2) has a special property that becomes visible
while invoking the definition (1-42) for two additional closed systems that
execute cycles while in contact with pairs of temperature reservoirs as
figure (1-17) .
Fig.1-17: Two closed systems that execute cycles while in contact with pairs of
temperature reservoirs
−𝑄3
( ) = 𝑓(𝑇1 , 𝑇3 ) (1-43)
𝑄1 𝑟𝑒𝑣
−𝑄3
( ) = 𝑓(𝑇2 , 𝑇3 ) (1-44)
𝑄2 𝑟𝑒𝑣
26
Dividing equation (1-43) by equation (1-44), along with using equation
(1-42) to eliminate (−Q2∕Q1)rev , yields :
𝑓(𝑇1 ,𝑇3 )
𝑓(𝑇1 , 𝑇2 ) = (1-45)
𝑓(𝑇2 ,𝑇3 )
The left side of eq. (1-45) does not depend on the constant T3 ; therefore,
the analytical form of f (T1,T2) is :
𝜓 (𝑇1 )
𝑓(𝑇1 , 𝑇2 ) =
𝜓(𝑇2 )
(1-46)
−𝑄3 𝜙𝑇2
( ) = (1-47)
𝑄1 𝑟𝑒𝑣 𝜙𝑇1
Equation (1-47) can be rewritten to express the ratio between the heat
interactions of a reversible cycle that absorbs one unit of energy (Q0) from
a reference reservoir of empirical temperature 𝜃0 and rejects the amount
–Q to an arbitrary reservoir (𝜃):
−𝑄
𝜙(𝜃) = 𝜙(𝜃0 ) ( ) (1-48)
𝑄0 𝑟𝑒𝑣
The numerical values obtained in this manner for function 𝜙 constitute the
thermodynamic temperature scale (symbol T); in other words, ≡ T or :
−𝑄
𝑇 = 𝑇0 ( ) (1-49)
𝑄0
One of the properties of the ideal gas Carnot cycle sketched in the P−v
plane is :
−𝑄2𝐶 𝜃2,𝐼𝑑𝑒𝑎𝑙 𝑔𝑎𝑠
= (1-50)
𝑄1𝐶 𝜃1,𝐼𝑑𝑒𝑎𝑙 𝑔𝑎𝑠
27
Comparing equation (1-47) with equation (1-50), we see that 𝜙 = 𝜃ideal gas :
𝜃𝑖𝑑𝑒𝑎𝑙 𝑔𝑎𝑠 = 𝑇 (1-52)
The reward for introducing the concept of a thermodynamic temperature
scale is that (−Q2∕Q1)rev =T2∕T1 and the second law (1-41) becomes :
𝑄1 𝑄2
+ ≤0 (1-53)
𝑇1 𝑇2
This new statement is general despite the assumption W > 0 made early in
the derivation of equation (1-38).
We see a sequence that leads to the second law for any cycle executed in
contact with any number of temperature reservoirs:
𝑄1 𝑄2 𝑄𝑛
+ + ⋯+ ≤0 (1-54)
𝑇1 𝑇2 𝑇𝑛
𝑄𝑖
∑𝑛+1
𝑖=1 ≤0 (1-55)
𝑇𝑖
Then the assumed equation (1-55) is valid. Consider, for this purpose, a
system A that executes a cycle while in thermal communication withn+1
temperature reservoirs ,T1,T2,...,Tn,Tn+1 . No assumption is made about
the direction (sign) of each of the heat transfer interactions,
Q1,Q2, . .., Qn,Qn+1.Letus return reservoirTn+1 to its original state by
placing it in contact with a reversible cycle (C) such that :
28
The second law for the composite system A+Tn+1+C, which
completes a cycle while in contact with temperature reservoirs,
we obtain :
𝑄𝑖 𝑄𝑛,𝐶
∑𝑛𝑖=1 + ≤0 (1-57)
𝑇𝑖 𝑇𝑛
𝑄𝑖 𝑄(𝑛+1),𝐶
∑𝑛𝑖=1 − ≤0 (1-58)
𝑇𝑖 𝑇(𝑛+1)
𝑄𝑖 𝑄(𝑛+1),𝐶
∑𝑛𝑖=1 + ≤0 (1-59)
𝑇𝑖 𝑇(𝑛+1)
𝛿𝑄
∮ ≤0 (1-60)
𝑇
𝛿𝑄𝑟𝑒𝑣
∮ ≤0 (1-61)
𝑇
29
reversible cycle, 𝛿Qrev ∕ T represents the change in a thermodynamic
property S:
δQrev
dS = (1-62)
T
This new property was named entropy by Clausius in 1865; however, the
same property was discovered and used earlier by Rankine , who called it
thermodynamic function, labeled it 𝜙, instead of S, and regarded equation
(1-62) as the “general equation of thermodynamics.”
Thinking of only that family of paths that are reversible (a path is reversible
if it can be a part of a reversible cycle), we integrate equation (1-62) and
obtain the definition of entropy change:
2 𝛿𝑄𝑟𝑒𝑣
𝑆1 − 𝑆2 = ∫1 (1-63)
𝑇
The arbitrary process 1→2 can be a part of a cycle 1→2→1, where the
return process 2→1 takes place along a reversible path; rewriting
equation(1-60) for this cycle and using the definition (1-63), we have, in
order,
2 𝛿𝑄 1 𝛿𝑄𝑟𝑒𝑣
∫1 + ∫2 ≤0 (1-64)
𝑇 𝑇
2 𝛿𝑄
∫1 𝑇
≤ 𝑆2 − 𝑆 1 (1-65)
Entropy Entropy
transfer change
(Nonproperty) (Property)
31
(i.e., a new definition) is the entropy generation or entropy production Sgen
, a quantity that is never negative:
2 𝛿𝑄
𝑆𝑔𝑒𝑛 = 𝑆2 − 𝑆1 − ∫1 ≥0 (1-66)
𝑇
and
𝑄̇𝑖
∆𝑆𝑔𝑒𝑛 = 𝑆𝑜𝑝𝑒𝑛 ,(𝑡+∆𝑡) − 𝑆𝑜𝑝𝑒𝑛 ,𝑡 − ∆𝑡 + (𝑚𝑠
̇ )𝑜𝑢𝑡 ∆𝑡 − (𝑚𝑠
̇ )𝑖𝑛 ∆𝑡 ≥ 𝑜
𝑇𝑖
(1-70)
Invoking the limit Δt →0, dropping the subscript “open,” and assuming the
existence of any number of spots with heat transfer i and mass flow (in,
out) on the control surface, we obtain :
𝑑𝑠 𝑄̇𝑖
̇
𝑆𝑔𝑒𝑛 = − ∑𝑖 + ∑𝑜𝑢𝑡 𝑚𝑠 ̇ 0
̇ − ∑𝑖𝑛 𝑚𝑠 ≥ (1-71)
𝑑𝑡 𝑇𝑖
31
1.10.4 Third law of thermodynamics
Figure (1- 18) shows the shape of the calculated S (T , ℋ) curves for a
paramagnetic salt such as gadolinium sulfate or gadolinium gallium garnet
at temperatures below 1 K . An extensive comparison of the various
properties of paramagnetic materials contemplated for magnetic
refrigeration systems has been provided by Barclay and Steyert.
Fig.1-18: Constant ℋ lines on the (T–S) plane of a paramagnetic salt and the process
of adiabatic demagnetization.
32
Fig.1-19: Unattainability of absolute zero by means of a finite sequence of isothermal
magnetization and adiabatic demagnetization processes of the type sown in Fig.1-18
From the construction of figure (1-19) is that the temperature T=0 cannot
be reached through a finite sequence of operations. Because it stems from
in the fact that the constant-ℋ lines intersect at T=0. The same geometric
feature can be expressed as:
33