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


15.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 3.24(b), based on a compressor exit temperature of $ 1200^\circ\textrm{F}$ (922 K).

Figure 15.3: Schematic of adiabatic flame temperature
Image fig7AdiabaticFlameTemperature_web

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 in Figure 15.3, where temperature is plotted versus the percentage completion of the reaction. The initial state is $ i$ and the final state is $ f$ , with the final state at a higher temperature than the initial state. The solid line in the figure shows a representation of the ``actual'' process.

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 the enthalpy changes in Process (1) is

$\displaystyle h_2 -h_i=- q_1 =(h_f^\circ)_\textrm{unit 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,

$\displaystyle h_f -h_2 =q_1,$

or, if we can approximate the specific heat as constant (using some appropriate average value)

$\displaystyle c_\textrm{p,avg} (T_f -T_2)=q_1.$

For the overall process there is no work done and no heat exchanged so that the difference in enthalpy between initial and final states is zero:

$\displaystyle \Delta h_1 +\Delta h_2 = \Delta h_\textrm{adiabatic} = 0.$

The temperature change during this second process is therefore given by (approximately)

$\displaystyle (T_f -T_2)= \frac{q_1}{c_\textrm{p,avg}} = \frac{\vert(h_f^\circ)_\textrm{unit mass}\vert}{c_\textrm{p,avg}}.$ (15..5)

The value of the adiabatic flame temperature given in Equation (15.5) 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_\textrm{p,avg}$ , 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^\circ\textrm{C}$ with 400% theoretical air. The reaction is

\begin{multline*}
C_8H_{18}(l) +12.5O_2 +12.5(3.76N_2)+3[12.5O_2 +12.5(3.76N_2)]
\longrightarrow \\
8CO_2 +9H_2 O(g) +37.5O_2 +188N_2.
\end{multline*}

For an adiabatic process

$\displaystyle \sum_R n_i(\overline{h}_f^\circ+\Delta\overline{h})_i = \sum_P n_...
...rline{h}} _{\substack{\textrm{At adiabatic} \textrm{flame temperature}}})_e .$ (15..6)

We can again think of the general process in steps:
  1. Bring reactants to $ 25^\circ\textrm{C}$ [the term $ (\Delta \overline{h})_i$ ] from the initial temperature, using whatever heat transfer, $ q_a$ , is needed. In this example we do not need step (i) because we are already at the reference temperature.
  2. Reaction at $ 25^\circ C$ [the term $ (h_f^\circ)_
\textrm{reactants$\longrightarrow$products}$ ]. There will be some heat transfer in this step, $ q_b$ , out of the combustor.
  3. Put back heat $ q_a +q_b$ into the products of combustion. The resulting temperature is the adiabatic flame temperature.

In the present case Equation (15.6) is, explicitly:

$\displaystyle (\overline{h}_f^\circ)_{C_8H_{18}}(l) = 8(\overline{h}_f^\circ)_{...
...rline{h}_{H_2O} +37.5\Delta \overline{h}_{O_2}
+188\Delta \overline{h}_{N_2}\}.$

We can examine the terms in the SFEE separately, starting with the heat of formation terms, and keeping track of units:

$\displaystyle \overline{h}_f^\circ : \qquad$ $\displaystyle 8(\overline{h}_f^\circ)_{CO_2} + 9(\overline{h}_f^\circ)_{H_2O} -(\overline{h}_f^\circ)_{C_8H_{18}}(l)$    
  $\displaystyle = 8\textrm{ kmole} (-393,522\textrm{ kJ/kmole}) + 9\textrm{ kmole} (-241,827\textrm{ kJ/kmole})$    
  $\displaystyle \qquad- 1\textrm{ kmole}(-249,952\textrm{ kJ/kmole})$    
  $\displaystyle = -5.075 \times 10^6 \textrm{ kJ}.$    

The exit state at the adiabatic flame temperature is specified by:

$\displaystyle \sum_P n_e \Delta\overline{h}_e = 5.075 \times 10^6\textrm{ kJ}.$

We find the adiabatic flame temperature in three ways:
  • an approximate solution using an average value of $ c_p$ ,
  • a more accurate one using the tabulated evolution of $ c_p$ with temperature,
  • or a more accurate solution using the tabulated values for gas enthalpy in Table A.8 of SB&VW.

15.5.1 Approximate solution using ``average'' values of specific heat

From Figure 15.2 we can use the values at 500 K as representative. These are:

Gas $ c_p$ (kJ/kmole)
$ CO_2$ 45
$ H_2O$ 35
$ O_2$ 30
$ N_2$ 30
Using $ \Delta \overline{h} = \overline{c}_\textrm{p,avg} \Delta T$ ,

$\displaystyle \sum_P n_e \Delta \overline{h}_e =\Delta T\left\{8(c_p)_{CO_2} +9(c_p)_{H_2O}
+37.5(c_p)_{O_2} +188(c_p)_{N_2}\right\},$

where $ \Delta T =
T_\textrm{final} -25^\circ\textrm{C} =T_\textrm{final} -298\textrm{
K}.$

$\displaystyle \sum_P n_e \Delta \overline{h}_e = 7440 \Delta T \textrm{ kJ/K},$

and using the exit state calculated above, find that

$\displaystyle \Delta T = 682\textrm{ K} \quad \Rightarrow \quad T_\textrm{final} = 980\textrm{ K}.$

15.5.2 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/\overline{R}$ with $ T$ (kJ/kmol)
$ CO_2$ $ 2.401+(8.735\times10^{-3})T-(6.607\times10^{-6})T^2+(2.002\times10^{-9})T^3$
$ H_2O$ $ 4.070-(1.108\times10^{-3})T+(4.152\times10^{-6})T^2-(2.964\times10^{-9})T^3+(0.807\times10^{-12})T^4$
$ O_2$ $ 3.626-(1.878\times10^{-3})T+(7.055\times10^{-6})T^2-(6.764\times10^{-9})
T^3+(2.156\times10^{-12})T^4$
$ N_2$ $ 3.675-(1.208\times10^{-3})T+(2.324\times10^{-6})T^2-(0.632\times10^{-9})T^3-(0.226\times10^{-12})T^4$
Using

$\displaystyle \Delta h =\int_\textrm{298 K}^\textrm{$T_f$ in K} c_p(T) dT$

and the same equation as above, we obtain

$\displaystyle T_f= 899\textrm{ K}.$

15.5.3 Solution for adiabatic flame temperature using tabulated values for gas enthalpy

  $ \Delta \overline{h}_{CO_2}$ $ \Delta \overline{h}_{H_2O}$ $ \Delta \overline{h}_{O_2}$ $ \Delta
\overline{h}_{N_2}$
$ T= 900\textrm{ K}$ 28,041 21,924 19,246 18,221 kJ/kmole
$ T=1000\textrm{ K}$ 33,405 25,978 22,707 21,460 kJ/kmole

Plugging in the numbers shows the answer is between these two conditions. Linearly interpolating gives a value of

$\displaystyle T_\textrm{final} = 962\textrm{ K}.$





Muddy Points

Does ``adiabatic flame temperature'' assume 100% combustion? (MP 15.7)

What part of the computation for adiabatic flame temperature involves iteration? (MP 15.8)

UnifiedTP