B.A.V. Bennett and M.D. Smooke,

Combustion Theory & Modelling, 2, 221–258, 1998

Abstract

Within realistic combustion devices, physical quantities may change by an order of magnitude over an extremely thin flame front, while remaining nearly unchanged throughout large areas nearby. Such behavior dictates the use of adaptive numerical methods. The recently developed local rectangular refinement (LRR) solution-adaptive gridding method produces robust unstructured rectangular grids, utilizes novel multiple-scale finite difference discretizations, and incorporates a damped modified Newton's method for simultaneously solving systems of governing elliptic PDEs. Here, the LRR method is applied to two axisymmetric laminar flames: a premixed Bunsen flame with one-step chemistry; and a diffusion flame employing various complex chemical mechanisms. The Bunsen flame's position depends highly upon grid spacing, especially on coarse grids; it stabilizes only with adequate refinement. The diffusion flame results show excellent agreement with experimental data for flame structure, temperature, and major species. For both flames, the LRR results on intermediate grids are comparable to those obtained on equivalently-refined conventional grids. Solution accuracy on the final LRR grids could not be achieved using conventional grids because the latter exceeded available computer memory. In general, the LRR method required about half the grid points, half the memory, and half the computation time of the solution process on conventional grids.