2.2 Heats of Reaction, Formation, and Adiabatic Flame Temperatures

Heats of Reaction, Formation, and Adiabatic Flame Temperatures

Chemical transformations involve energy exchanges that are orders of magnitude greater than simple sensible heat transfers. Whether designing a cooling jacket to prevent thermal runaway in a continuous stirred-tank reactor (CSTR) or calculating the fuel efficiency of a fired process heater, professional chemical engineers rely on precise thermochemical balances.

---

1. Thermochemistry Fundamentals: Hess's Law and Standard States

The standard heat of reaction ($\Delta H^\circ_{rxn}$) is the enthalpy change when stoichiometric quantities of reactants at standard state ($T^\circ = 298.15\text{ K} = 25^\circ\text{C}$ and $P^\circ = 1\text{ bar}$) react completely to form products at the same temperature and pressure.

Standard Heat of Formation Formulation

The standard heat of formation ($\Delta H^\circ_{f,i}$) is the enthalpy change associated with the formation of one mole of a compound from its constituent elements in their naturally occurring reference states at $25^\circ\text{C}$ and $1\text{ bar}$. By definition, $\Delta H^\circ_f = 0$ for all pure elemental species in their standard reference state (e.g., $\text{O}_2(g)$, $\text{N}_2(g)$, $\text{H}_2(g)$, and $\text{C}(graphite)$).

By Hess's Law, the reaction enthalpy is:

$\Delta H^\circ_{rxn} = \sum_{i} \nu_i \Delta H^\circ_{f,i} = \sum_{products} |\nu_i| \Delta H^\circ_{f,i} - \sum_{reactants} |\nu_i| \Delta H^\circ_{f,i}$

Where $\nu_i$ is the stoichiometric coefficient (conventionally positive for products, negative for reactants).

• If $\Delta H^\circ_{rxn} < 0$, the reaction is exothermic (releases heat to surroundings).
• If $\Delta H^\circ_{rxn} > 0$, the reaction is endothermic (absorbs heat from surroundings).

Standard Heat of Combustion Formulation

For fuels, hydrocarbons, and organic compounds, heats of combustion ($\Delta H^\circ_{c,i}$) are frequently tabulated. The heat of combustion is the enthalpy released during complete oxidation of a compound to form $\text{CO}_2(g)$, $\text{H}_2\text{O}(l)$, $\text{SO}_2(g)$, and $\text{N}_2(g)$ at $25^\circ\text{C}$.

When using combustion data, the signs reverse relative to formation data:

$\Delta H^\circ_{rxn} = -\sum_{i} \nu_i \Delta H^\circ_{c,i} = \sum_{reactants} |\nu_i| \Delta H^\circ_{c,i} - \sum_{products} |\nu_i| \Delta H^\circ_{c,i}$

---

2. Higher Heating Value (HHV) versus Lower Heating Value (LHV)

In combustion engineering, water is produced as a combustion product. The state of this water dictates the reported heating value of the fuel:

1. Higher Heating Value (Gross Heating Value, HHV): All water produced by combustion is assumed to condense into liquid water at $25^\circ\text{C}$. The latent heat of vaporization of water is recovered.
2. Lower Heating Value (Net Heating Value, LHV): All water produced remains in the vapor state at $25^\circ\text{C}$. The latent heat of vaporization is not credited.

The Mathematical Conversion

At $25^\circ\text{C}$ ($298.15\text{ K}$), the latent heat of vaporization of water is: $\Delta H_{vap,\text{H}_2\text{O}}(298.15\text{ K}) = 44.012\text{ kJ/mol} = 2,442.5\text{ kJ/kg} = 1,050\text{ Btu/lb}$

The two heating values are related by:

$\text{HHV} = \text{LHV} + n_{\text{H}_2\text{O}} \times \Delta H_{vap,\text{H}_2\text{O}}(25^\circ\text{C})$

Where $n_{\text{H}_2\text{O}}$ is the moles of water generated per mole of fuel burned. In real combustion processes, flue gases exit boilers and furnaces at temperatures well above $100^\circ\text{C}$ ($150^\circ\text{C}$ to $250^\circ\text{C}$) to avoid acid dew point condensation. Therefore, water leaves as vapor, and LHV represents the practical thermal energy available to the process.

---

3. Kirchhoff's Law: Enthalpy of Reaction at Operating Temperature

Chemical reactors rarely operate at $25^\circ\text{C}$. To evaluate the actual heat of reaction $\Delta H_{rxn}(T)$ at operating temperature $T$, we construct a state-function cycle:

1. Cool reactants from $T$ to $T_{ref}$ ($25^\circ\text{C}$): $\Delta H_1 = -\int_{T_{ref}}^{T} \sum_{react} |\nu_i| C_{p,i} dT$
2. React at standard state $T_{ref}$: $\Delta H_2 = \Delta H^\circ_{rxn}(T_{ref})$
3. Heat products from $T_{ref}$ to $T$: $\Delta H_3 = \int_{T_{ref}}^{T} \sum_{prod} |\nu_i| C_{p,i} dT$

Summing these terms yields Kirchhoff's Law:

$\Delta H_{rxn}(T) = \Delta H^\circ_{rxn}(T_{ref}) + \int_{T_{ref}}^{T} \Delta C_p(T) dT$

Where the difference in heat capacities between products and reactants is: $\Delta C_p(T) = \sum_{i} \nu_i C_{p,i}(T) = \sum_{products} |\nu_i| C_{p,i}(T) - \sum_{reactants} |\nu_i| C_{p,i}(T)$

• If $\Delta C_p > 0$: $\Delta H_{rxn}(T)$ becomes algebraically more positive (less exothermic / more endothermic) as temperature rises.
• If $\Delta C_p < 0$: $\Delta H_{rxn}(T)$ becomes algebraically more negative (more exothermic) as temperature rises.

---

4. Reactor Energy Balance Formulations

There are two mathematically equivalent frameworks for solving continuous, steady-state non-isothermal reactor energy balances:

Method A: Heat of Reaction Method

This method uses standard heats of reaction evaluated at $T_{ref} = 25^\circ\text{C}$ combined with sensible heat paths:

$\dot{Q} - \dot{W}_s = \dot{\xi} \Delta H^\circ_{rxn}(T_{ref}) + \sum_{out} \dot{n}_{i,out} \int_{T_{ref}}^{T_{out}} C_{p,i} dT - \sum_{in} \dot{n}_{i,in} \int_{T_{ref}}^{T_{in}} C_{p,i} dT$

Where $\dot{\xi}$ is the extent of reaction ($\text{mol/s}$ or $\text{kmol/h}$), defined for any participating species $i$ as: $\dot{\xi} = \frac{\dot{n}_{i,out} - \dot{n}_{i,in}}{\nu_i}$ This method is ideal when a single primary reaction occurs and its conversion or extent is readily known.

Method B: Heat of Formation Method

This method computes absolute total stream enthalpies referenced to elements at $25^\circ\text{C}$:

$\dot{Q} - \dot{W}_s = \sum_{out} \dot{n}_{i,out} \hat{H}_{i}(T_{out}) - \sum_{in} \dot{n}_{i,in} \hat{H}_{i}(T_{in})$

Where the specific molar enthalpy of compound $i$ at temperature $T$ is: $\hat{H}_i(T) = \Delta H^\circ_{f,i}(298.15\text{ K}) + \int_{298.15}^{T} C_{p,i}(T) dT$ Method B is immensely powerful for complex multi-reaction networks, reforming furnaces, and combustion chambers because you do not need to determine extents of reaction or track multiple heats of reaction; the energy balance is satisfied directly from inlet and outlet atom/species balances.

---

5. Adiabatic Flame Temperature (AFT)

The adiabatic flame temperature ($T_{ad}$) is the maximum theoretical temperature achieved by a combustion system when fuel and oxidizer react completely in a thermally insulated chamber with no shaft work ($\dot{Q} = 0$, $\dot{W}_s = 0$):

$\Delta \dot{H} = 0 \implies \dot{H}_{products}(T_{ad}) = \dot{H}_{reactants}(T_{in})$

If reactants enter at $T_{ref} = 25^\circ\text{C}$ and combustion is complete:

$-\Delta H^\circ_{rxn}(298.15\text{ K}) = \sum_{i,products} n_{i,out} \int_{298.15}^{T_{ad}} C_{p,i}(T) dT$

Using mean product heat capacities $\bar{C}_{p,i}$ over the temperature span:

$T_{ad} = 298.15\text{ K} + \frac{-\Delta H^\circ_{rxn}(298.15\text{ K})}{\sum_{i} n_{i,out} \bar{C}_{p,i}}$

Engineering Factors Influencing AFT

1. Excess Air: Atmospheric air consists of approximately $21\text{ mol}% \text{ O}_2$ and $79\text{ mol}% \text{ N}_2$ ($3.762\text{ mol N}_2 / \text{mol O}_2$). Any air fed beyond stoichiometric requirements introduces excess $\text{O}_2$ and a large volume of non-reactive $\text{N}2$ ballast. This ballast absorbs combustion energy as sensible heat, drastically suppressing $T{ad}$.
2. Air Preheating: Preheating the incoming air above ambient temperature delivers additional sensible enthalpy $\dot{H}{reactants}$, elevating $T{ad}$ and increasing radiative heat transfer in furnaces.
3. Thermal Dissociation: At temperatures above $1,700^\circ\text{C}$ ($2,000\text{ K}$), combustion products dissociate endothermically ($\text{CO}_2 \rightleftharpoons \text{CO} + \frac{1}{2}\text{O}_2$ and $\text{H}_2\text{O} \rightleftharpoons \text{H}_2 + \frac{1}{2}\text{O}_2$). This dissociation absorbs heat, capping real flame temperatures $100^\circ\text{C}$ to $250^\circ\text{C}$ below unreacted theoretical models.

---

6. Summary Comparison: Thermochemical Energy Formulations

---

7. Worked Numerical Example: Natural Gas Combustor with Excess Air

Problem Statement

Methane ($\text{CH}4$) fuel gas at $25.0^\circ\text{C}$ is burned completely with $20.0%$ excess dry air in an adiabatic, continuous combustion chamber. The combustion air is preheated to $T{air,in} = 200.0^\circ\text{C}$ ($473.15\text{ K}$) using hot flue gas heat recovery.

Reaction: $\text{CH}_4(g) + 2\text{O}_2(g) \to \text{CO}_2(g) + 2\text{H}_2\text{O}(g)$

Thermochemical Data:

• $\Delta H^\circ_{f,\text{CH}_4(g)} = -74.85\text{ kJ/mol}$
• $\Delta H^\circ_{f,\text{CO}_2(g)} = -393.51\text{ kJ/mol}$
• $\Delta H^\circ_{f,\text{H}_2\text{O}(g)} = -241.82\text{ kJ/mol}$
• Mean heat capacity of air between $25^\circ\text{C}$ and $200^\circ\text{C}$: $\bar{C}_{p,air} = 29.50\text{ J/(mol}\cdot\text{K)}$

Mean Flue Gas Heat Capacities between $25.0^\circ\text{C}$ and $T_{ad}$:

• $\bar{C}_{p,\text{CO}_2} = 51.50\text{ J/(mol}\cdot\text{K)}$
• $\bar{C}_{p,\text{H}_2\text{O}(g)} = 41.20\text{ J/(mol}\cdot\text{K)}$
• $\bar{C}_{p,\text{O}_2} = 33.80\text{ J/(mol}\cdot\text{K)}$
• $\bar{C}_{p,\text{N}_2} = 32.10\text{ J/(mol}\cdot\text{K)}$

Calculate:

1. The standard heat of reaction at $25.0^\circ\text{C}$ (LHV basis).
2. The molar composition of the product flue gas per mole of $\text{CH}_4$ burned.
3. The sensible enthalpy added by the preheated combustion air.
4. The resulting adiabatic flame temperature ($T_{ad}$).

Step-by-Step Solution

Step 1: Standard Reaction Enthalpy (LHV Basis)

$\Delta H^\circ_{rxn} = [1\times (-393.51) + 2\times (-241.82)] - [1\times (-74.85) + 2\times (0)]$ $\Delta H^\circ_{rxn} = [-393.51 - 483.64] - [-74.85] = -877.15 + 74.85 = -802.30\text{ kJ/mol CH}_4$

Step 2: Material Balance on Combustor (Basis: $1.00\text{ mol CH}_4$)

• Stoichiometric $\text{O}_2$ required: $2.00\text{ mol}$
• With $20.0%$ excess air: $\text{O}_2\text{ fed} = 2.00 \times 1.20 = 2.40\text{ mol}$
• $\text{N}_2\text{ fed} = 2.40 \times (79.0 / 21.0) = 2.40 \times 3.7619 = 9.029\text{ mol}$
• Total air fed = $2.40 + 9.029 = 11.429\text{ mol air / mol CH}_4$

Exhaust Gas Moles:

• $n_{\text{CO}_2} = 1.00\text{ mol}$
• $n_{\text{H}_2\text{O}} = 2.00\text{ mol}$
• $n_{\text{O}_2,excess} = 2.40 - 2.00 = 0.40\text{ mol}$
• $n_{\text{N}_2} = 9.029\text{ mol}$
• Total dry + wet flue gas = $1.00 + 2.00 + 0.40 + 9.029 = 12.429\text{ mol}$

Step 3: Sensible Enthalpy of Preheated Air

Air enters at $200.0^\circ\text{C}$ ($\Delta T_{air} = 200.0 - 25.0 = 175.0\text{ K}$): $\Delta H_{air,in} = n_{air} \bar{C}_{p,air} \Delta T_{air} = (11.429\text{ mol}) \times (0.02950\text{ kJ/(mol}\cdot\text{K)}) \times (175.0\text{ K}) = 58.995\text{ kJ}$

Step 4: Product Mixture Total Heat Capacity

$\sum n_i \bar{C}_{p,i} = (1.00)(51.50) + (2.00)(41.20) + (0.40)(33.80) + (9.029)(32.10)$ $\sum n_i \bar{C}_{p,i} = 51.50 + 82.40 + 13.52 + 289.83 = 437.25\text{ J/K} = 0.43725\text{ kJ/K}$ Notice that $\text{N}_2$ accounts for $289.83 / 437.25 = 66.3%$ of the total heat absorption capacity of the flue gas!

Step 5: Adiabatic Energy Balance for $T_{ad}$

Because the combustor is adiabatic ($\dot{Q} = 0$, $\dot{W}s = 0$): $\Delta H_{total} = 0 \implies -\Delta H^\circ_{rxn} + \Delta H_{air,in} = \sum n_i \bar{C}_{p,i} (T_{ad} - 298.15)$ $802.30\text{ kJ} + 58.995\text{ kJ} = (0.43725\text{ kJ/K}) (T_{ad} - 298.15)$ $861.295\text{ kJ} = (0.43725\text{ kJ/K}) (T_{ad} - 298.15)$ $\Delta T_{ad} = \frac{861.295}{0.43725} = 1969.8\text{ K}$ $T_{ad} = 298.15 + 1969.8 = 2268.0\text{ K} = 1994.8^\circ\text{C} \approx 1995^\circ\text{C}$ If air had not been preheated, $\Delta T{ad} = 802.30 / 0.43725 = 1834.8\text{ K}$, giving $T_{ad} = 2133\text{ K} = 1860^\circ\text{C}$. Preheating the combustion air by $175^\circ\text{C}$ boosted flame temperature by $135^\circ\text{C}$.

---

8. Common PE Exam Traps & Calculation Pitfalls

Trap 1: Forgetting Nitrogen in Combustion Flue Gas
Fuel burns with air, not pure oxygen (unless explicitly specified as oxy-combustion). For every mole of $\text{O}_2$ consumed, $3.762\text{ moles of N}_2$ accompany it into the firebox. Omitting $\text{N}_2$ drops the product gas heat capacity by two-thirds, resulting in preposterously high calculated flame temperatures ($>4,000^\circ\text{C}$).

Trap 2: Mixing Up Higher and Lower Heating Values
When computing boiler efficiency ($\eta = \dot{Q}{steam} / (\dot{m}{fuel} \times \text{Heating Value})$), in the United States, utility power plant boilers are traditionally rated on an HHV basis, whereas European installations and gas turbine systems use LHV. Using HHV when LHV is specified (or vice versa) produces a $5%$ to $11%$ error on fuel consumption calculations.

Trap 3: Inverting Signs in Hess's Law When Using Heats of Combustion
Formations follow: $\Delta H^\circ_{rxn} = \sum \text{Products} - \sum \text{Reactants}$. Combustions follow: $\Delta H^\circ_{rxn} = \sum \text{Reactants} - \sum \text{Products}$. Reversing this sign convention makes exothermic reactions appear endothermic.