Combustion and Flame, Vol.123, No.3, 402-411, 2000
Dependence of CO emissions on the rate of product cooling
Detailed chemical kinetics simulations of the cooling of hydrocarbon combustion products show that the mole fraction of CO freezes well above equilibrium levels at a value that depends on the cooling rate. The frozen mole fraction of CO can be expressed as a function of the system pressure and the mole fractions of CO2, O-2, and H2O X-CO,X-f = 2.0 x 10(-9)[X-CO2/(XO2XH2O)](P tau)(-1.6) where tau is the temperature relaxation time, inversely proportional to the cooling rate. The ratio of mole fractions in this result can also be derived by partial equilibrium analysis. However, this derivation yields X-CO,X-f proportional to X-OH(2) rather than (P tau)(-1.6). The dependence on cooling rate is shown to be a kinetic effect due to the relatively slow forward rate of the principal CO oxidation reaction. This rate is proportional to PXOH, and X-CO freezes when X-OH is too small to maintain the consumption rate for CO that is needed to keep X-CO near equilibrium at a given cooling rate. Partial equilibrium analysis shows that X-OH,X-f, the value of X-OH when X-CO freezes, is proportional to (P tau)(-1). the temperature, T-f where X-OH equals X-OH,X-f is found by setting X-OH,X-f = X-OH,X-0 exp(-gamma(OH)/T-f), where gamma(OH) is a fitting parameter related to equilibrium constants. Using T-f and the equilibrium value of gamma(OH) with the temperature-dependent result of the partial equilibrium analysis gives X-CO proportional to (P tau)(-1.56). This gives values of frozen X-CO in reasonable agreement with the detailed chemical kinetics results. At high cooling rates, gamma(OH) is a function of tau that can be found by a curve fit to X-OH vs. T data from detailed kinetics calculations. Using this function to calculate T-f improves the agreement with the detailed chemical kinetics results.