Combustion Science and Technology, Vol.141, No.1-6, 59-81, 1999
Study of alternative descriptions of methane oxidation for CFD modeling of turbulent combustors
Four alternative ways to describe methane oxidation that all can be used in computational fluid dynamics (CFD) have been investigated. The purpose has been to investigate how advanced a description of the methane oxidation is needed for an accurate description of the main components in practical CFD calculations. Especially the prediction of the intermediate combustion products CO and H-2 has been stressed, The four alternative approaches are: I) thermodynamic equilibrium, II) a three-step simplified mechanism including infinitely fast in;eversible global reaction steps for the oxidation of methane and hydrogen, and a global, kinetically controlled, reaction step for carbon monoxide oxidation, III) a four-step simplified mechanism including two global reaction steps for the break-down of the hydrocarbon, a global reaction step for the oxidation of hydrogen and a global reaction step describing the reversible water-shift reaction and IV) a comprehensive reaction mechanism including 32 chemical species and 156 reversible elementary reactions. The comparison has been made by studying the concentration predictions for perfectly stirred reactor (PSR) conditions as a function of nominal residence time and by a CFD study of a diffusion flame. The results from the PSR computations indicate that the three-step simplified mechanism is a good alternative for the fuel lean conditions and for flames far from extinction. The four-step simplified mechanism was found to be a good compromise between accuracy and computational speed. The four-step simplified mechanism also performed well at moderately fuel rich conditions. However, in the CFD study of the diffusion flame, the difference in the results obtained with the four-step simplified formulation and the comprehensive mechanism was surprisingly large. The assumption of instantaneous local thermodynamic equilibrium in the reacting fraction of the fluid proved to be an unsuccessful alternative in the CFD modeling.