Reacting Systems and Combustion Thermodynamics

Introduction and Context

A reacting system is a thermodynamic system in which the chemical identity of one or more constituents changes. Combustion is the most important engineering example: a fuel reacts with an oxidizer, usually air, and releases chemical energy as heat and work. Gas turbines, reciprocating engines, boilers, furnaces, rocket engines, industrial heaters, waste-to-energy plants, and many propulsion systems depend on the analysis of reacting systems.

Combustion thermodynamics connects the first and second laws of thermodynamics with chemical equilibrium, material balances, adiabatic flame temperature, exergy, emissions, and process efficiency. At postgraduate level, the subject must be treated not merely as a heat-release problem but as a coupled problem involving composition, temperature, pressure, irreversibility, transport, and environmental constraints.

Learning Objectives

After studying this chapter, the reader should be able to:

1. formulate mass and energy balances for reacting closed and open systems;
2. use stoichiometric equations, equivalence ratio, and air–fuel ratio correctly;
3. distinguish heating value, enthalpy of formation, sensible enthalpy, and chemical exergy;
4. calculate adiabatic flame temperature using appropriate assumptions;
5. explain chemical equilibrium through Gibbs free energy and equilibrium constants;
6. evaluate combustion efficiency, exergy destruction, and pollutant formation;
7. identify limitations of idealized combustion models and select suitable computational methods.

1. Chemical-Reaction Representation

A general reaction may be written as

where Aᵢ denotes chemical species, and νᵢ′ and νᵢ″ are stoichiometric coefficients on the reactant and product sides. The net stoichiometric coefficient is νᵢ = νᵢ″ − νᵢ′. A positive νᵢ denotes a product and a negative value denotes a reactant.

For complete combustion of methane with oxygen,

If water is treated as vapour, this equation represents the higher-temperature combustion state. If liquid water is formed after cooling, the thermodynamic energy released is greater because condensation releases latent heat.

Air is commonly approximated as 21 mol% oxygen and 79 mol% nitrogen. Thus,

The stoichiometric methane–air reaction is therefore

The nitrogen is approximately inert in a first idealization, although at high temperature it participates in dissociation and nitric-oxide formation.

Stoichiometric and Actual Mixtures

The stoichiometric air–fuel ratio is the mass of air required for complete reaction with a unit mass of fuel:

For methane, using molecular masses 16.04 kg/kmol for CH₄ and 28.97 kg/kmol for air,

The equivalence ratio is

A mixture is fuel-rich when φ > 1 and fuel-lean when φ < 1. The air ratio, often denoted λ, is the reciprocal:

These definitions must not be confused with excess-air percentage, which is

and applies principally to λ ≥ 1.

2. Conservation Laws for Reacting Systems

Species Mass Conservation

For a control volume at steady state, the balance on species i is

where ṅ is molar flow rate. Chemical reaction changes species internally, so the overall mass balance remains valid but individual species balances include reaction generation and consumption.

Atomic conservation is more fundamental than species conservation. For every chemical element k,

where aₖᵢ is the number of atoms of element k in species i and nᵢ is the number of moles. This equation is valid even when equilibrium produces many species.

First-Law Balance

For a closed reacting system,

is not the most convenient general form unless chemical internal energies are explicitly included. A safer formulation is to assign each species a total specific enthalpy that includes formation and sensible contributions. For a steady-flow control volume, neglecting kinetic and potential-energy changes,

where Ẇ_s is shaft-power output and hᵢ is the molar enthalpy of species i. This equation is valid for reacting and non-reacting flow when the enthalpy reference is used consistently.

For an adiabatic combustor with no shaft work,

The apparent temperature rise is therefore determined by the chemical energy of the reactants, heat capacities, dissociation, and inlet conditions.

3. Enthalpy of Formation and Heating Value

The standard enthalpy of formation, ∆h°f,i, is the enthalpy change required to form one mole of species i from its elements in their standard reference states, commonly at 298.15 K and 1 bar. The standard enthalpy change of a reaction is

The enthalpy of a species at temperature T can be represented as

where cₚ,i is molar heat capacity and Tref is the reference temperature. If cₚ is treated as constant,

This constant-heat-capacity approximation is useful for preliminary calculations but may produce significant error at flame temperatures.

The lower heating value (LHV) assumes that product water remains vapour. The higher heating value (HHV) assumes that product water condenses to liquid. Hence HHV exceeds LHV for hydrogen-containing fuels. The appropriate value depends on whether condensation heat is recovered in the device.

4. Adiabatic Flame Temperature

The adiabatic flame temperature is the product temperature obtained when combustion occurs adiabatically, with specified reactant temperature, pressure, composition, and product model. It is not necessarily the actual temperature in an engine because heat transfer, finite-rate chemistry, dissociation, mixing, and radiation alter the result.

For reactants initially at T₀ and products at T_ad, the energy equation is

For complete combustion, fixed composition, constant heat capacities, and reactants at Tref,

where ∆H°ᵣₓₙ is negative for an exothermic reaction. This expression is a special case and should not be used when dissociation or strongly temperature-dependent heat capacities are important.

Worked Example: Approximate Methane Flame Temperature

Consider stoichiometric methane–air at 298 K and 1 bar. Take the LHV-based heat release as approximately 802,000 kJ/kmol CH₄. Assume complete products CO₂, 2H₂O(g), and 7.52N₂, with representative constant molar heat capacities: CO₂ = 55, H₂O = 44, and N₂ = 35 kJ/(kmol·K). The product heat-capacity rate per kmol of methane is

Thus,

This is an illustrative estimate. Actual equilibrium calculations generally predict a lower temperature because high-temperature dissociation absorbs energy and because heat capacities increase with temperature.

5. Entropy, Irreversibility, and Exergy

The entropy balance for a steady control volume is

where Ṡgen ≥ 0 by the second law. In combustion, entropy generation arises from chemical reaction, mixing, heat transfer across finite temperature differences, friction, and pressure loss.

The exergy destroyed is given by the Gouy–Stodola theorem:

where T₀ is the environmental temperature. Exergy destruction measures the lost potential to produce useful work. A combustor may have high thermal efficiency in the narrow sense while destroying substantial chemical exergy because fuel is converted irreversibly into hot products.

The specific flow exergy of a stream may be written as

The first two terms are physical exergy; eₓ,chemical accounts for departure from environmental chemical composition. Chemical exergy is essential when comparing fuels, carbon capture options, and fuel-conversion technologies.

6. Chemical Equilibrium

Complete combustion is an engineering approximation. At high temperature, products such as CO₂ and H₂O dissociate:

The equilibrium composition minimizes total Gibbs free energy at specified temperature, pressure, and elemental composition. For a reaction at equilibrium,

The standard Gibbs energy change is related to the equilibrium constant K by

where R is the universal gas constant and T is absolute temperature. For ideal gases, the reaction quotient is constructed from partial pressures:

where p° is the standard pressure. At equilibrium, Qp = Kp. For non-ideal gases, fugacity replaces partial pressure. Equilibrium constants depend strongly on temperature and, for reactions involving different total mole numbers, the equilibrium composition also depends on pressure.

A practical equilibrium calculation requires: elemental balances, an equation of state, thermodynamic property data, equilibrium relations, and non-negativity constraints on species amounts. Newton methods, Gibbs-minimization algorithms, and validated chemical-kinetics software are commonly used.

7. Combustion Efficiency and Emissions

Combustion efficiency may be defined on an energy basis as

Its numerical value depends on whether CO, unburned hydrocarbons, soot, and hydrogen are included. A high-temperature flame can still have poor efficiency if mixing is inadequate or residence time is insufficient.

Important pollutants include:

• CO and unburned hydrocarbons: indicators of incomplete oxidation, poor mixing, quenching, or inadequate residence time.
• NOₓ: thermal NOₓ increases strongly with high temperature and oxygen availability; prompt and fuel-bound mechanisms may also contribute.
• SOₓ: produced from sulfur in the fuel.
• Particulate matter and soot: associated especially with locally fuel-rich regions and heavy hydrocarbons.
• CO₂: not conventionally toxic at ambient concentrations but central to greenhouse-gas accounting.

Emission reduction methods include lean premixed combustion, exhaust-gas recirculation, staged combustion, water or steam injection, selective catalytic reduction, oxidation catalysts, improved atomization, and carbon capture. Each method creates trade-offs. For example, lowering flame temperature can reduce thermal NOₓ but may increase CO and unburned hydrocarbons.

8. Combustion in Practical Devices

A real combustor combines thermodynamics with fluid mechanics, heat transfer, and chemical kinetics. Key design variables include pressure, inlet temperature, residence time, turbulence, spray characteristics, wall temperature, cooling-air fraction, and fuel distribution.

In gas turbines, compressor pressure ratio raises reactant temperature and affects equilibrium, combustor stability, and NOₓ. In internal-combustion engines, ignition timing, spray penetration, charge stratification, and cyclic variability are decisive. In boilers and furnaces, excess air, flame shape, radiative transfer, and furnace-wall temperature strongly influence performance.

The Damköhler number compares flow time to chemical time:

Large Da generally indicates that chemistry is fast relative to flow, while small Da indicates that finite-rate chemistry may control conversion. The interpretation depends on the selected characteristic scales and reaction mechanism; it is not a universal single-value criterion.

9. Worked Application: Lean Methane Combustion

A combustor burns methane with 20% excess air. Determine the equivalence ratio and write the idealized product equation, assuming complete combustion.

Twenty percent excess air means λ = 1.20. Therefore,

The stoichiometric oxygen requirement is 2 kmol O₂ per kmol CH₄. The actual oxygen supplied is 2(1.20) = 2.40 kmol. Hence,

The residual oxygen is evidence of lean operation. The nitrogen quantity is 2.40 × 3.76 = 9.024 kmol. This equation is a complete-combustion model; it does not predict CO, NOₓ, dissociation, or soot.

10. Illustrative Case Study: Industrial Gas-Turbine Combustor

Background

An industrial gas turbine operates with natural gas at elevated pressure. Operators observe increasing NOₓ during periods of high load, while reducing fuel flow causes unstable combustion and increased CO.

Problem

The design team must reduce NOₓ without violating combustion-stability limits or sacrificing turbine output.

Analysis

A first-law analysis shows that the turbine inlet temperature is controlled primarily by fuel flow, compressor discharge conditions, and dilution air. An equilibrium calculation predicts that NO formation is favoured by high flame temperature and excess oxygen. However, equilibrium alone cannot predict actual emissions because NO formation is kinetically limited and depends on local temperature histories. Computational fluid dynamics with finite-rate chemistry is therefore required after preliminary thermodynamic screening.

Approach

The team evaluates lean premixed operation, improved fuel staging, and controlled exhaust-gas recirculation. The design is constrained by flashback, blow-off, pressure oscillations, liner cooling, and allowable CO. Temperature distribution is monitored rather than relying only on average flame temperature. A validated chemical mechanism and calibrated emissions model are used, with uncertainty bounds for fuel composition and sensor error.

Outcome

Fuel staging produces a more uniform lean mixture and reduces peak temperature. A modest increase in residence time limits CO, while acoustic monitoring prevents operation near instability boundaries. The final operating map reports NOₓ, CO, efficiency, pressure loss, and stability margin together rather than optimizing one metric.

Learning Points

1. Average adiabatic flame temperature is insufficient for emissions prediction.
2. Thermodynamic, kinetic, transport, and control models must be integrated.
3. Environmental performance is a constrained multi-objective design problem.
4. Experimental validation is necessary because reaction mechanisms and turbulence models contain uncertainty.

11. Limitations, Governance, and Contemporary Developments

Ideal-gas assumptions become questionable at high pressure, and constant heat capacities can create large flame-temperature errors. Complete-combustion models conceal local rich zones and pollutant precursors. Equilibrium predicts the ultimate composition, not the time required to reach it. Detailed kinetics may be computationally expensive, while reduced mechanisms may omit important pathways.

Thermodynamic property data must be traceable to recognized databases. Engineers should document reference states, fuel composition, pressure, temperature, product assumptions, and uncertainty. Emissions calculations must comply with applicable air-quality regulations and reporting standards. Safety analysis must address explosion limits, ignition sources, toxic exposure, hydrogen embrittlement where relevant, and carbon-monoxide hazards.

Hydrogen and ammonia combustion are important contemporary topics. Hydrogen has no carbon in the fuel but can produce substantial NOₓ because of high flame temperature and rapid chemistry. Ammonia offers carbon-free fuel potential but has low flame speed, difficult ignition, ammonia slip, and NOₓ challenges. Sustainable aviation fuels, oxy-fuel combustion, chemical looping, supercritical CO₂ cycles, electrified heating, and carbon capture are active research areas.

Machine learning is increasingly used for surrogate chemistry, reduced-order combustion models, sensor-based control, and anomaly detection. It should supplement, not replace, conservation laws, thermodynamic consistency, uncertainty quantification, and experimental validation. A model that predicts emissions accurately within its training range may fail under fuel switching, pressure changes, or unusual transient operation.

Best-Practice Workflow

1. Define the system boundary, steady or transient operation, and desired outputs.
2. Specify fuel composition, oxidizer composition, pressure, temperature, and reference state.
3. Balance elements before calculating energy or equilibrium.
4. Select complete combustion, equilibrium, or finite-rate chemistry according to the question.
5. Use temperature-dependent properties for high-temperature systems.
6. Check energy, atom, charge, and unit balances.
7. Report assumptions and conduct sensitivity analysis.
8. Validate important predictions against measurements or certified databases.
9. Evaluate efficiency, emissions, safety, and exergy together.

Chapter Summary

Reacting-system analysis combines elemental conservation, the first law, the second law, chemical thermodynamics, and combustion chemistry. Stoichiometric equations establish composition; air–fuel ratio, equivalence ratio, and excess air describe mixture strength. Enthalpies of formation and sensible enthalpies provide a consistent energy framework, while the adiabatic flame temperature is obtained from an energy balance subject to clearly stated assumptions. At high temperatures, dissociation and chemical equilibrium must be considered through Gibbs free energy and equilibrium constants. Entropy generation and exergy destruction quantify irreversibility. Practical combustion requires attention to kinetics, mixing, transport, heat transfer, emissions, safety, regulation, and uncertainty. Modern low-carbon combustion systems therefore require integrated thermodynamic, computational, experimental, and control approaches.

Key Terms

Reacting system: A system in which chemical species change through reaction.

Stoichiometric reaction: Reaction with exact proportions for complete chemical conversion.

Air–fuel ratio: Mass of air divided by mass of fuel.

Equivalence ratio: Actual fuel–air ratio divided by the stoichiometric fuel–air ratio.

Excess air: Air supplied beyond the stoichiometric requirement.

Enthalpy of formation: Enthalpy required to form a species from reference-state elements.

Heating value: Energy released by complete combustion of a unit quantity of fuel.

LHV: Heating value with product water remaining vapour.

HHV: Heating value including condensation of product water.

Adiabatic flame temperature: Product temperature for an adiabatic combustion process.

Chemical equilibrium: State of minimum Gibbs free energy at specified conditions.

Dissociation: Breakdown of a compound into simpler species at high temperature.

Reaction quotient: Composition-based expression used to assess reaction direction.

Entropy generation: Measure of irreversibility in a process.

Exergy: Maximum theoretical useful work relative to an environment.

Exergy destruction: Lost work potential equal to T₀Ṡgen.

Finite-rate chemistry: Description in which reaction rates are explicitly considered.

Thermal NOₓ: Nitrogen oxides formed primarily through high-temperature oxidation of atmospheric nitrogen.

Damköhler number: Ratio of flow time scale to chemical time scale.

Combustion efficiency: Fraction of supplied chemical energy released as intended.

Review Questions

5 Short Answer Questions

1. Define equivalence ratio and explain the significance of φ = 1.
2. Why does HHV exceed LHV for hydrogen-containing fuels?
3. State the assumptions behind a constant-heat-capacity adiabatic flame-temperature calculation.
4. What is the physical meaning of exergy destruction?
5. Why can equilibrium calculations fail to predict actual pollutant emissions?

5 Long Answer Questions

1. Derive the steady-flow energy equation for an adiabatic reacting control volume and discuss the required enthalpy reference convention.
2. Explain the role of enthalpy of formation in calculating the heat of combustion and adiabatic flame temperature.
3. Develop the equilibrium formulation for a dissociating ideal-gas combustion mixture using elemental balances and Gibbs free energy.
4. Compare the thermodynamic and kinetic controls on CO, NOₓ, soot, and unburned hydrocarbons.
5. Discuss exergy destruction in a gas-turbine combustor and identify practical methods for reducing it.

3 Analytical / Application Questions

1. Propane burns with 10% excess air. Determine λ, φ, the idealized complete-combustion equation, and the residual oxygen coefficient.
2. A combustor operates at elevated pressure with preheated reactants. Explain how pressure and preheat affect flame temperature, dissociation, stability, and NOₓ, distinguishing equilibrium effects from kinetic effects.
3. Design a validation plan for a reduced-order model of a hydrogen-fuelled combustor. Specify measurements, conservation checks, uncertainty sources, and criteria for accepting the model.

Multiple Choice Questions with Answers

Question 1

For a mixture with 20% excess air, the equivalence ratio is:

A. 0.833

B. 1.20

C. 0.200

D. 1.833

Correct Answer: A. 0.833

Question 2

The stoichiometric oxygen requirement for one mole of methane is:

A. 1 mol O₂

B. 2 mol O₂

C. 3 mol O₂

D. 4 mol O₂

Correct Answer: B. 2 mol O₂

Question 3

The relation between exergy destruction and entropy generation is:

A. ĖD = Ṡgen/T₀

B. ĖD = T₀Ṡgen

C. ĖD = T₀/Ṡgen

D. ĖD = 0 for every real process

Correct Answer: B. ĖD = T₀Ṡgen

Question 4

For an adiabatic steady combustor with no shaft work and negligible kinetic and potential energy changes, which balance applies?

A. ∑in ṅᵢhᵢ = ∑out ṅᵢhᵢ

B. ∑in ṅᵢhᵢ = 0

C. Q̇ = Ẇs

D. ∑out ṅᵢ = 0

Correct Answer: A. ∑in ṅᵢhᵢ = ∑out ṅᵢhᵢ

Question 5

The equilibrium condition for a reacting mixture at fixed temperature and pressure is most directly associated with:

A. Maximum Gibbs free energy

B. Minimum Gibbs free energy

C. Zero enthalpy of formation for every species

D. Maximum entropy of each individual species

Correct Answer: B. Minimum Gibbs free energy

Question 6

If λ = 0.8, the mixture is:

A. Fuel-lean

B. Stoichiometric

C. Fuel-rich

D. Chemically inert

Correct Answer: C. Fuel-rich

Question 7

In the idealized methane–air equation, the approximate nitrogen coefficient for stoichiometric combustion is:

A. 2.00

B. 3.76

C. 7.52

D. 9.52

Correct Answer: C. 7.52

Question 8

LHV differs from HHV mainly because LHV assumes that product water:

A. Is dissociated into hydrogen and oxygen

B. Remains as vapour

C. Is absent from the products

D. Is converted to carbon dioxide

Correct Answer: B. Remains as vapour

Question 9

A Damköhler number much greater than unity generally indicates that:

A. The flow time is much greater than the chemical time

B. The chemical time is zero

C. No reaction can occur

D. The mixture must be stoichiometric

Correct Answer: A. The flow time is much greater than the chemical time

Question 10

Which change generally tends to reduce thermal NOₓ in a combustor?

A. Increasing peak flame temperature

B. Creating localized high-temperature zones

C. Lowering peak temperature through lean or staged combustion

D. Increasing residence time at maximum temperature without limit

Correct Answer: C. Lowering peak temperature through lean or staged combustion