Part 2.C: Introduction to Thermochemistry
[SB&VW-14.1-14.6]
Until now, we have specified the heat given to the devices analyzed, and not concerned
ourselves with how this heat might be produced. In this section, we examine the issue of how we
obtain the heat needed for work production. For the most part, this is from converting chemical
energy into heat, so the discussion will be on reacting mixtures of gas which are involved in
chemical combustion processes.
The topic addressed is “thermochemistry”, which is the combining of thermodynamics with
chemistry to predict such items as how much heat is released from a chemical reaction. This is the
“Q” or “q” that we have used in the cycle analysis. The principal components of the approach are
use of a chemical balance plus the steady flow energy equation (SFEE) which equates the sum of
shaft work (from) and heat transfer (to) a control volume to the difference in control volume inlet
and exit enthalpy fluxes.
2.C.1 Fuels
There are a wide variety of fuels used for aerospace power and propulsion. A primary one
is jet fuel (octane, essentially kerosene) which has the chemical formula CH
18
. Other fuels we
8
consider are hydrogen
(H
2
)
and methane ( CH
4
).
The chemical process in which a fuel, for example methane, is burned consists of (on a very basic
level—there are many intermediate reactions that need to be accounted for when computations of
the combustion process are carried out):
CH
4
+ 2O
2
→ CO
2
+ 2H O .
2
(Reactants) (Products)
The reactions we describe are carried out in air, which can be approximated as 21% O
2
and
79% N
2
. This composition is referred to as “theoretical air”. There are other components of air
(for example Argon, which is roughly 1%), but the results given using the theoretical air
approximation are more than adequate for our purposes. With this definition, for each mole of O
2
,
3.76 (79/21) moles of N
2
are involved:
CH
4
+ 2O
2
+ 23 76 )N → CO
2
+ 2H O + 7.52N
2
( .
2 2
Even if the nitrogen is not part of the combustion process, it leaves the combustion chamber at the
same temperature as the other products, and this change in state (change in enthalpy) needs to be
accounted for in the steady flow energy equation. At the high temperatures achieved in internal
combustion engines (aircraft and automobile) reaction does occur between the nitrogen and
oxygen, which gives rise to oxides of nitrogen, although we will not consider these reactions.
The condition at which the mixture of fuel and air is such that both completely participate
in the reaction is called stoichiometric. In gas turbines, excess air is often used so that the
2C-1
temperatures of the gas exiting the combustor is kept to within desired limits (see Figures A-8, A-
9, A-11 in Part 1 for data on these limits.)
Muddy points
Why is there 3.76 N
2
? (MP 2C.1)
What is the most effective way to solve for the number of moles in the reactions? (MP
2C.2)
2.C.2 Fuel-Air Ratio
The reaction for aeroengine fuel at stoichiometric conditions is:
CH +12.5 O
2
+12 5 (376 ) N →8 CO
2
+9 H O + 47.0 N . .
2818 2 2
On a molar basis, the ratio of fuel to air is [1/(12.5+47.0)] = 1/59.5 = 0.0167.
To find the ratio on a mass flow basis, which is the way in which the aeroengine industry discusses
it, we need to “weight” the molar proportions by the molecular weight of the components. The
fuel molecular weight is 114, the oxygen molecular weight is 32 and the nitrogen molecular weight
is (approximately) 28. The fuel/air ratio on a mass flow basis is thus
1 ×114
Fuel-air ratio = = 0 0664 .
12 5 ×32 +12 5 ×376 ×28 . . .
If we used the actual constituents of air we would get 0.0667, a value about 0.5% different.
Muddy points
Do we always assume 100% complete combustion? How good an approximation is this?
(MP 2C.3)
2.C.3 Enthalpy of formation
The systems we have worked with until now have been of fixed chemical composition.
Because of this, we could use thermodynamic properties relative to an arbitrary base, since all
comparisons could be made with respect to the chosen base. For example, the specific energy
u
f
(.
o
.001 C) = 0 0 for steam. If there are no changes in composition, and only changes in
properties of given substances, this is adequate. If there are changes in composition, however, we
need to have a reference state so there is consistency for different substances.
The convention used is that the reference state is a temperature of 25
o
C (298 K) and a pressure of
0.1 MPa. (These are roughly room conditions.) At these reference conditions, the enthalpy of the
elements (oxygen, hydrogen, nitrogen, carbon, etc.) is taken as zero.
The results of a combustion process can be diagrammed as below. The reactants enter at
standard conditions; the combustion (reaction) takes place in the volume indicated. Downstream
of the reaction zone there is an appropriate amount of heat transfer with the surroundings so that
the products leave at the standard conditions. For the reaction of carbon and oxygen to produce
CO
2
, the heat that has to be extracted is Q
CV
=?393 522 kJ/kmole ; this is heat that comes out of,
the control volume.
2C-2
C0
2
C + 0
2
1 kmole C
1 kmole C0
2
25
o
C, 0.1 MPa
25
o
C, 0.1 MPa Vo l ume
1 kmole C0
2
= -393,522 KJ, heat is out of control volumeQ
cv
Figure C-1: Constant pressure combustion
There is no shaft work done in the control volume and the first law for the control volume (SFEE)
reduces to:
mass flow of enthalpy in + rate of heat addition = mass flow of enthalpy out.
We can write this statement in the form
˙
∑mh +Q
˙
CV
=∑m h
e
(C.3.1)˙
ii e
R P
In Eq. (C.3.1) the subscripts “R” and “P” on the summations refer to the reactants (R) and products
(P) respectively. The subscripts on the mass flow rates and enthalpies refer to all of the
components at inlet and at exit.
The relation in terms of mass flows can be written in molar form, which is often more convenient
for reacting flow problems, by using the molecular weight, M
i
, to define the molar mass flow rate,
n˙
i
, and molar enthalpy, h
i
, for any individual i
th
(or e
th
) component as
n˙
i
= m˙
i
/ M ; mass flow rate in terms of kmoles/sec
i
h
i
= M h
i
; enthalpy per kmole
i
The SFEE is, in these terms,
˙
∑nh +Q
˙
CV
=∑n h
e
. (C.3.2)˙
ii e
R P
The statements that have been made do not necessarily need to be viewed in the context of flow
processes. Suppose we have one unit of C and one unit of O
2
at the initial conditions and we carry
out a constant pressure reaction at ambient pressure, P
amb
. If so,
2C-3
= Q ?WU
final
?U
initial
= Q
CV
? P
amb
(
V
final
?V
initial
)
,
since P = P
f
= P
amb
. Combining terms,
i
U
final
+ P
final
V
final
?
(
U
initial
+ P
initial
V
initial
)
= Q
CV
,
or,
= Q
CV
.H
final
? H
initial
In terms of the numbers of moles and the specific enthalpy this is
∑nh +Q
CV
=∑n h
e
(C.3.3)
ii e
R P
The enthalpy of CO
2
, at 25
o
C and 0.1 MPa, with reference to a base where the enthalpy of the
o
elements is zero, is called the enthalpy of formation and denoted by h
f
. Values of the heat of
formation for a number of substances are given in Table A.9 in SB&VW.
The enthalpies of the reactants and products for the formation of CO
2
are:
= h
C
= 0
h
h
O
2
o
f
,For one kmole: Q
CV
=∑n h
e
= H
P
=
()
=?393 522 kJ/kmole.
e
P
CO
2
The enthalpy of CO
2
in any other state (T,P)is given by
h
TP
=
()
+
(
?h
)
.
,
h
0
f
298K,. 0 1 Mpa→T P 0 1 MPa 298K,. ,
These descriptions can be applied to any compound. For elements or compounds that exist in
more than one state at the reference conditions (for example, carbon exists as diamond and as
graphite), we also need to specify the state.
Note that there is a minus sign for the heat of formation. The heat transfer is out of the control
,volume and is thus negative by our convention. This means that h
elements
> =? 393 522 kJ .h
CO
2
Muddy points
Is the enthalpy of formation equal to the heat transfer out of the combustion during the
formation reaction? (MP 2C.4)
Are the enthalpies of H
2
and H (monoatomic hydrogen) both zero at 298K? (MP 2C.5)
2C-4
2.C.4 First Law Analysis of Reacting Systems
The form of the first law for the control volume is (there is no shaft work):
˙
∑nh+Q
˙
CV
=∑nh .˙
ii e e
R P
This is given in terms of the moles of the different constituents, and it reduces to the more familiar
form for a single fluid (say air) with no reactions occurring. We need to specify one parameter as
the basis of the solution; 1 kmole of fuel, 1 kmole of air, 1 kmole total, etc. We use 1 kmole of
fuel as the basic unit and examine the burning of hydrogen.
H
2
+20 → 2HO
2 2
The reactants and the products are both taken to be at 0.1MPa and 25
o
C, so the inlet and
exit P and T are specified. The control volume is the combustion chamber. There is no shaft work
done and the SFEE is in the form of Equation (C.1.2). The enthalpy of the entering gas is zero for
both the hydrogen and the oxygen (elements have enthalpies defined as zero at the reference state).
If the exit products are in the gaseous state, the exit enthalpy is therefore related to the enthalpy of
formation of the product by:
o
h
f
n˙
e
HO
h
e
HO
= n˙
e
()
2
()HO HOg
2
2
2
= 2 x (- 241,827)kJ = - 483,654 kJ; gaseous state at exit.
If the water is in a liquid state at the exit of the process:
o
h
f
n˙
e
HO
h
e
HO
= n˙
e
()
2
()HO HOl
2
2
2
= 2 x (- 285,783) kJ = - 571, 676.
There is more heat given up if the products emerge as liquid. The difference between the two
values is the enthalpy needed to turn the liquid into gas at 25
o
C: h
fg
= 2442 kJ/kmole.
A more complex example is provided by the burning of methane (natural gas) in oxygen,
producing .
CH
4
+2O
2
→ CO
2
+2H O ()l
2
The components in this reaction equation are three ideal gases (methane, oxygen, and CO
2
) and
liquid water. We again specify that the inlet and exit states are at the reference conditions so that:
∑nh h
f
o
,
ii
=
()
=?74 873 kJ
R
CH
4
ee
=
()
+2
()
nh h∑
f
o
CO
2
h
f
o
HOl
P
2
()
Q
=-393,522 + 2(-285,838) = -965,198 kJ
CV
=?965 198 kJ – (-74,873) = -890,325 kJ.
,
2C-5
Suppose the substances which comprise the reactants and the products are not at 25
o
C and 0.1MPa.
If so, the expression that connects the reactants and products is;
????
?
?
?
?
?
h
o
f
+
h
P
i
and reference
conditions
?
?
?
?
?
?
{ e
?
?
?
?
?
h
o
f
+
h
P
e
and reference
conditions
?
?
?
?
?
?
{
.
(C.4.1)Q
CV
+∑ =∑
n n
i
R P
Between Between
T
i
T
e
, ,
i e
Equation (C.4.1) shows that we must compute the enthalpy difference ?h between the reference
conditions and the given state if the inlet or exit conditions are not the reference pressure and
temperature.
There are different levels of approximation for the computation: (a) assume the specific
heat is constant over the range at some average value, (b) use the polynomial expressions (Table
A.6) in the integral, and (c) use tabulated values. The first is the simplest and the crudest.
Combustion processes often involve changes of a thousand degrees or more and, as Figure C-2
shows, the specific heat for some gases can change by a factor of two or more over this range,
although the changes for air are more modest. This means that, depending on the accuracy desired,
one may need to consider the temperature dependence of the specific heat in computing ?h .
Figure C-2: Specific heat as a function of temperature [from SB&VW]
Muddy points
When doing cycle analysis, do we have to consider combustion products and
their effect on specific heat ratio (γ is not 1.4)? (MP 2C.6)
2C-6
2.C.5 Adiabatic Flame Temperature
For a combustion process that takes place adiabatically with no shaft work, the temperature
of the products is referred to as the adiabatic flame temperature. This is the maximum temperature
that can be achieved for given reactants. Heat transfer, incomplete combustion, and dissociation,
all result in lower temperature. The maximum adiabatic flame temperature for a given fuel and
oxidizer combination occurs with a stoichiometric mixture (correct proportions such that all fuel
and all oxidizer are consumed). The amount of excess air can be tailored as part of the design to
control the adiabatic flame temperature. The considerable distance between present temperatures
in a gas turbine engine and the maximum adiabatic flame temperature at stoichiometric conditions
is shown in Figure A-11 of Part 1, based on a compressor exit temperature of 1200
o
F (922 K).
An initial view of the concept of adiabatic flame temperature is provided by examining two
reacting gases, at a given pressure, and
asking what the end temperature is. The
process is shown schematically at the right,
1
2
?h
1
Constant P
a
Actual path
?h
f = Final state
where temperature is plotted versus the
percentage completion of the reaction.
T
?h
2
The initial state is i and the final state is f,
with the final state at a higher Constant P
temperature than the initial state. The
State i 2
solid line in the figure shows a
representation of the “actual” process.
0
Percentage 100%
To see how we would arrive at the final
completion
state the dashed lines break the state
of reaction
change
into two parts. Process (1) is reaction
at constant T and P. To carry out such a
process, we would need to extract heat.
Suppose the total amount of heat extracted
per unit mass is q
1
. The relation between Schematic of adiabatic flame temperature
the enthalpy changes in Process (1) is
h
2
?=? q
1
=
()unit
h
i
h
o
f
mass
where q
1
is the “heat of reaction”.
For Process (2), we put this amount back into the products to raise their temperature to the final
level. For this process, h ?h
2
= q
1
, or, if we can approximate the specific heat as constant (using
f
some appropriate average value) c
p
av.
(
T ?T
2
)
= q
1
. For the overall process there is no work done
f
and no heat exchanged so that the difference in enthalpy between initial and final states is zero:
?h
1
+?h = 0.
2
=?h
adiabatic
The temperature change during this second process is therefore given by (approximately)
2C-7
o
h
f
()unit
mass
1
(
T ?T
2
)
=
q
= . (C.5.1)
f
c c
p
av.
p
av.
The value of the adiabatic flame temperature given in Equation (C.5.1) is for 100%
completion of the reaction. In reality, as the temperature increases, the tendency is for the degree
of reaction to be less than 100%. For example, for the combustion of hydrogen and oxygen, at
high temperatures the combustion product (water) dissociates back into the simpler elemental
reactants. The degree of reaction is thus itself a function of temperature that needs to be computed.
We used this idea in discussing the stoichiometric ramjet, when we said that the maximum
temperature was independent of flight Mach number and hence of inlet stagnation temperature. It
is also to be emphasized that the idea of a constant (average) specific heat, c
p
av.
, is for illustration
and not inherently part of the definition of adiabatic flame temperature.
An example computation of adiabatic flame temperature is furnished by the combustion of liquid
octane at 25
o
C with 400% theoretical air. The reaction is
CH ()+12.5O +12.5 3.76N
2
)+3 12.5O +12.5 3.76N
2
)
]
→8CO +9H O g
818
l
2
(
[
2
(
2 2
()+37.5O
2
+188N
2
.
For an adiabatic process
o o
∑nh
f
+?h
)
=∑n
e
(
h
f
+?h
)
. (C.5.2)
i
(
i e
R P
At adiabatic flame temperature
We can again think of the general process in steps:
a) Bring reactants to 25
o
C
[
the term
(
?h
)
]
from the initial temperature, using whatever
i
heat transfer, q
a
, is needed. In this example we do not need step (i) because we are already
at the reference temperature.
o
b) Reaction at 25
o
C -
?
the term
()
?
. There will be some heat transfer in this step,h
?
f
reactants→products ?
??
q
b
, out of the combustor.
c) Put back heat q +q
b
into the products of combustion. The resulting temperature is the
a
adiabatic flame temperature.
In the present case Equation (C.1.6) is, explicitly:
o o o
h
f
()=8h
f
CO
2
+9h
f
+
{
?h
CO
2
+9?h
H
2
O
+37.5?h +188?h
N
2
}
C
8
H
18
l
H
2
O
O
2
We can examine the terms in the SFEE separately, starting with the heat of formation terms,
o o o o
h
f
: 8h
f
CO
2
+9h
f
?h
f
()= 8 (-393,522) + 9 (-241,827) – (-249,952)l
H
2
O C
8
H
18
2C-8
= -5.075 X 10
6
kJ/kmole.
The exit state at the adiabatic flame temperature is specified by:
nh∑ ? = 5.075 X 10
6
kJ/kmole
e e
P
We find the adiabatic flame temperature in three ways, approximate solution using an average
value of c
p
, a more accurate one using the tabulated evolution of c
p
with temperature and a more
precise solution using the tabulated values for gas enthalpy in Table A.8 of SB&VW.
a) Approximate solution using “average” values of specific heat:
From Figure C-2 we can use the values at 500K as representative. These are:
Gas c
p
(kJ/kmole)
CO
2
45
H
2
O 35
O
2
30
N
2
30.
Using ?hc
p
?T ,=
"ave"
? c c c c∑nh
e
=?T
{
8
()
+9
()
+37.5
()
+188
()
}
e p
CO
2
p
H
2
O
p
O
2
p
P
N
2
o
=where ?TT
final
?25 C = T
final
?298 K
?∑nh
e
= 7440?T kJ
e
P
?T = 682 ? T
final
= 980 K .
b) Solution for adiabatic flame temperature using evolutions of specific heats with temperature
Tables give the following evolutions of specific heats with temperature:
Gas Evolution of c
p
/ R with T (kJ/kmol)
CO
2
2.401+8.735.10
-3
xT-6.607.10
-6
xT
2
+2.002.10
-9
xT
3
H
2
O 4.070-1.108.10
-3
xT+4.152.10
-6
xT
2
-2.964.10
-9
xT
3
+0.807.10
-12
xT
4
O
2
3.626-1.878.10
-3
xT+7.055.10
-6
xT
2
-6.764.10
-9
xT
3
+2.156.10
-12
xT
4
N
2
3.675-1.208.10
-3
xT+2.324.10
-6
xT
2
-0.632.10
-9
xT
3
-0.226.10
-12
xT
4
T
f
in K
Using ?h =
∫
c ()T .dT and the same equation as above, we obtain:
p
298 K
T
f
= 899 K
c) Solution for adiabatic flame temperature using tabulated values for gas enthalpy:
?h
CO
2
?h
H
2
O
?h
O
2
?h
N
2
T= 900 K 28,041 21,924 19,246 18,221 kJ/kmole
T=1000K 33,405 25,978 22,707 21,460 kJ/kmole
2C-9
Plugging in the numbers shows the answer is between these two conditions. Linearly interpolating
gives a value of T
final
= 962 K .
Muddy points
Does "adiabatic flame temperature" assume 100% combustion? (MP 2C.7)
What part of the computation for adiabatic flame temperature involves iteration? (MP
2C.8)
2C-10
Muddiest points on part 2C
2C.1 Why is there 3.76 N
2
?
This is to represent the components other than oxygen that are in air. From SB&VW,
page 525, “The assumption that air is 21% oxygen and 79% nitrogen by volume leads to
the conclusion that for each mole of oxygen, 79/21 =3.76 moles of nitrogen are
involved.”
2C.2 What is the most effective way to solve for the number of moles in the reactions?
What we are doing is basically counting atoms on both sides of the reactants?products
statement. See Section 14.2 of SB&VW for a description of the combustion process and
for going from ratios in terms of moles to ratios in terms of mass.
2C.3 Do we always assume 100% complete combustion? How good an approximation is
this?
In the problems we do, we will only consider 100% combustion. It is a good
approximation for the range of problems that we address.
2C.4 Is the enthalpy of formation equal to the heat transfer out of the combustion during
the formation reaction?
As defined, the enthalpy of formation relates to a process in which the initial and final
states are at the same temperature. If there is combustion in between, heat will have to
be removed for this condition to occur. The enthalpy of formation is equal to the
negative of the magnitude of the heat outflow. If we consider the combustion as
occurring in a control volume, then per kmole h
f
o
=?Q
cv
, where Q
cv
is the heat transfer
out of the control volume per kmole. This is not in accord with our convention, and if
you please feel free to transform it back into the notation we have used before. (I find
that if I do this there are too many minus signs to keep track of easily.)
2C.5 Are the enthalpies of H
2
and H (monoatomic hydrogen) both zero at 298K?
The enthalpies of the elements are taken as zero at 298 and 0.1 MPa. In some cases there
are more than one form of the element. In that case the form chosen to have the value of
zero is that which is chemically stable at the reference state. The other forms then have
an enthalpy which is consistent with the reaction that produces this form of the element.
For hydrogen H
2
has zero enthalpy at the reference conditions and H has an enthalpy of
217,999 kJ/kmole (see Table A.8 in SB&VW), consistent with the idea that energy has to
be supplied to break the molecule apart.
2C.6 When doing a cycle analysis, do we have to consider combustion products and their
effect on specific heat ratio (γ is not 1.4)?
The specific heat ratio does depend on combustion products but the effect is not large
because the fuel air ratio is small. For example, for conditions of fuel air ratio ).034, the
specific heat ratio at room temperature is about 1.38. A larger variation encountered in
practice is with temperature; for a temperature of 1750K the specific heat ratio of pure air
is 1.3.
2C.7 Does "adiabatic flame temperature" assume 100% combustion?
Yes. This is the maximum temperature that could be produced. Incomplete combustion
will lower the temperature, as will heat transfer out of the combustion region.
2C.8 What part of the computation for adiabatic flame temperature involves iteration?
If the specific heat is not a simple analytic function of temperature (i.e., suppose it is
known only in tabular form), we cannot get a closed form solution for the adiabatic flame
temperature. We can, however, readily solve the enthalpy balance (SFEE) numerically
(this is where the iteration comes in) to find at what temperature the products have to
come out to have the same enthalpy (including the enthalpy of formation) as the
reactants. We did not do this calculation yet, but we will do it to show what the iteration
is all about.
Remember that the assumption of constant specific heat is just that, an assumption.
While this is an excellent assumption for many practical problems, if the precision of the
answer needed is very high, or if the range of temperatures is large (see Figure 5.11 in
SB&VW), then we cannot assume constant specific heat.