Combustion Reaction Calculator

Heats of combustion and stoichiometry for common fuels

Combustion Reaction Calculator: Complete Guide

The complete combustion of an organic compound CcHhOoNnSs in excess oxygen produces CO2, H2O, N2 and SO2. The balanced equation requires: moles of O2 = c + h/4 - o/2 + s - 3n/4 (where N oxidises to N2, not NOx). This calculator balances the equation, calculates the heat of combustion using Hess's law, and gives the stoichiometric masses of oxygen consumed and products formed from any mass of fuel.

Balancing combustion equations

For complete combustion of CcHhOo to CO2 and H2O: CcHhOo + (c + h/4 - o/2) O2 -> c CO2 + (h/2) H2O. Example 1: methane CH4. O2 = 1 + 4/4 = 2. CH4 + 2 O2 -> CO2 + 2 H2O. Example 2: octane C8H18. O2 = 8 + 18/4 = 12.5. C8H18 + 12.5 O2 -> 8 CO2 + 9 H2O. Multiply through by 2: 2 C8H18 + 25 O2 -> 16 CO2 + 18 H2O. Example 3: glucose C6H12O6. O2 = 6 + 3 - 3 = 6. C6H12O6 + 6 O2 -> 6 CO2 + 6 H2O. Example 4: sulfur-containing compound ethanethiol C2H5SH (C2H6S). O2 = 2 + 6/4 + 1 = 4.5. 2 C2H6S + 9 O2 -> 4 CO2 + 6 H2O + 2 SO2. The coefficient of O2 may be a half-integer -- multiply all coefficients by 2 to give integers if required by your course convention.

Heat of combustion by Hess's law

Standard enthalpy of combustion (delta-Hc) is calculated by Hess's law: delta-Hc = sum(products delta-Hf) - sum(reactants delta-Hf). For CcHhOo + O2 -> c CO2 + (h/2) H2O (liquid): delta-Hc = c x delta-Hf(CO2, g) + (h/2) x delta-Hf(H2O, l) - delta-Hf(fuel) - 0 (O2 elements = 0). Standard values at 298 K: delta-Hf(CO2, g) = -393.51 kJ/mol; delta-Hf(H2O, l) = -285.83 kJ/mol. Example: octane C8H18, delta-Hf = -250.1 kJ/mol. delta-Hc = 8 x (-393.51) + 9 x (-285.83) - (-250.1) = -3148.1 - 2572.5 + 250.1 = -5470.5 kJ/mol. Specific energy = 5470.5 / 114.23 (M_r C8H18) = 47.9 kJ/g. The heat of combustion for liquid water product (higher heating value) exceeds that for steam product (lower heating value) by the latent heat of vaporisation of water: 44.0 kJ/mol H2O x (h/2) moles.

Combustion in engines, power generation and safety

Air-fuel ratio (AFR): stoichiometric AFR = mass of air per mass of fuel for complete combustion. Air is 23.2% O2 by mass. For octane: theoretical O2 = 12.5 mol x 32 g/mol = 400 g per mol fuel (114.23 g). O2/fuel = 3.50 g/g. AFR = 3.50 / 0.232 = 15.1 (stoichiometric AFR for petrol engines). Engine management systems maintain lambda = actual AFR / stoichiometric AFR = 1.0 for three-way catalytic converter efficiency. Lean mixtures (lambda greater than 1) give lower emissions but higher NOx; rich mixtures (lambda less than 1) give CO and unburnt hydrocarbons. In safety engineering, the lower explosive limit (LEL) and upper explosive limit (UEL) define the flammable range for fuel-air mixtures: methane LEL = 5% v/v, UEL = 15%; hydrogen LEL = 4%, UEL = 75%. Below the LEL the mixture is too lean to ignite; above the UEL it is too rich.

Step-by-step worked example

A radiochemist measures the activity of a wood sample at 4.20 disintegrations per minute per gram of carbon (dpm/g C). Modern living wood gives 13.56 dpm/g C (the Libby standard). The half-life of C-14 is 5730 years (Cambridge standard). Step 1: calculate the activity ratio R = 4.20 / 13.56 = 0.3097. Step 2: apply the radiocarbon age formula: t = -t_half / ln(2) x ln(R) = -5730 / 0.6931 x ln(0.3097) = -8268 x (-1.1718) = 9688 years BP. Step 3: apply the reservoir correction if the sample is marine (typically -400 years) or from a known reservoir-effect region. Step 4: convert to calendar years using an established calibration curve (IntCal20 for terrestrial samples; Marine20 for marine samples). Calibration accounts for past variations in atmospheric C-14 production caused by solar activity and geomagnetic field changes. The calibrated date will carry an uncertainty of plus or minus 50 to 200 years at 2-sigma (95% confidence) depending on the sample age and position on the calibration curve. This calculated age of approximately 9700 years BP places the sample in the early Holocene -- consistent with early post-glacial woodland recolonisation in temperate regions. All calculations in this suite run locally in the browser with no data transmitted to any server.

Connections to related tools and broader chemistry suite

The LazyTools organic chemistry suite covers the major quantitative calculations in organic and physical organic chemistry. The Degree of Unsaturation Calculator connects structural formula to the number of rings and pi bonds. The Double Bond Equivalent Calculator derives molecular formulas from elemental analysis data. The Combustion Analysis Calculator converts CO2 and H2O masses from CHN analysis into empirical and molecular formulas. The Chemical Oxygen Demand Calculator links molecular structure to oxygen requirement, connecting synthetic chemistry to environmental engineering. The Crude Protein Calculator connects elemental nitrogen analysis to nutritional biochemistry. The Combustion Reaction Calculator balances combustion equations and calculates heat release from standard enthalpies of formation. Together these tools cover the quantitative skills required from A-level through undergraduate organic chemistry, and the specific calculations used daily in pharmaceutical, environmental, food and materials analysis laboratories. Results copy with one click for direct transcription into lab records, electronic laboratory notebooks or regulatory submissions.

Laboratory precision, significant figures and error propagation

In quantitative organic chemistry, attention to significant figures prevents false precision in reported results. The Kjeldahl nitrogen determination is typically accurate to plus or minus 0.1 to 0.3% absolute nitrogen content; this propagates to approximately plus or minus 1% crude protein at a 6.25 factor. Radiocarbon dating has an inherent analytical uncertainty of plus or minus 0.3 to 0.5% in the C-14/C-12 ratio (plus or minus 25 to 50 years), with additional systematic uncertainty from isotopic fractionation (corrected by measuring delta-C-13) and reservoir effects. Combustion analysis (CHN) is typically accurate to plus or minus 0.3% for C and H, permitting unambiguous formula determination only if the molecular formula mass is independently known. COD by dichromate titration has a precision of plus or minus 5 mg/L O2 or plus or minus 2% relative, whichever is greater, under standard laboratory conditions with calibrated glassware and standardised reagents.

Worked calculation and practical application

A process engineer needs to verify the inventory of a cryogenic ethylene storage tank. The insulated horizontal cylindrical tank has a working volume of 2500 m3 and is maintained at -120 deg C. Step 1: determine liquid ethylene density at -120 deg C using the NIST saturation data polynomial: approximately 457 kg/m3. Step 2: calculate the mass of ethylene: 2500 m3 x 457 kg/m3 = 1,142,500 kg = 1142.5 tonnes. Step 3: verify against custody transfer meter readings (typical accuracy plus or minus 0.1%). Step 4: calculate the vapour equivalent using the ideal gas law: 1,142,500 kg / 0.028054 kg/mol x 22.414 L/mol = 913 million litres of ethylene gas at STP (0 deg C, 1 bar). This inventory calculation is performed daily at ethylene terminals and cracker units for safety management, custody transfer and production scheduling. The density-temperature relationship is also essential for correcting custody transfer measurements from flow meters that measure volumetric flow at line temperature to mass flow for billing purposes.

Accuracy, data sources and range of validity

The density values in this calculator are fitted to data from the NIST Chemistry WebBook (https://webbook.nist.gov) for liquid ethylene (ethene) at the saturation pressure corresponding to each temperature. The polynomial fit reproduces the NIST data to within approximately 0.5 kg/m3 across the full liquid range from the melting point (-169.2 deg C) to the normal boiling point (-103.7 deg C). For temperatures outside this range -- supercooled liquid (below melting point) or gas phase (above boiling point at 1 bar) -- the calculation is not valid and an error is returned. For high-pressure liquid ethylene (above the saturation pressure at a given temperature), the compressed liquid density is slightly higher than the saturation value and can be calculated using an equation of state (Peng-Robinson, NIST REFPROP). For process engineering purposes at pressures below 10 MPa, the saturation density is typically used as a conservative estimate with less than 2% error. Temperature measurement should use calibrated platinum resistance thermometers (PRTs) with accuracy of plus or minus 0.1 deg C for inventory calculations at the precision required for custody transfer. All calculations in this tool are for engineering estimation purposes -- consult the current NIST WebBook or REFPROP software for safety-critical applications.

Frequently asked questions