Local Rectangular Refinement with Application to Nonreacting and Reacting Fluid Flow Problems

• Published in: Journal of computational physics 1999-05, Vol.151 (2), p.684-727
• Main Authors: Bennett, Beth Anne V, Smooke, Mitchell D
• Format: Article
• Language: English
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1006/jcph.1999.6214

A new solution-adaptive gridding method has been developed for the solution of discretized systems of coupled nonlinear elliptic partial differential equations on rectangular domains. Such a method is required for the numerical solution of realistic combustion problems, in which physical quantities may vary by orders of magnitude over one-tenth of a millimeter at atmospheric pressure, or over micrometers at higher pressures. The local rectangular refinement (LRR) method maintains orthogonality at grid-line intersections but lifts the tensor product restriction common to traditional grids, producing unstructured grids. Governing equations are discretized throughout the domain using newly derived forms, and Newton's method is used to solve the resulting system. On a simple test case with a known solution, the LRR method and its new discretizations are found to be more accurate than gridding methods representative of those appearing previously in the literature. For the more realistic problem of nonreacting driven square cavity flow, the LRR solution agrees very well with previously published data. When the LRR method is applied to a practical reacting flow (a rich axisymmetric laminar Bunsen flame with complex chemistry, multicomponent transport, and an optically thin radiation submodel), grid spacing highly influences the inner flame's position, which stabilizes only with adequate refinement. The vorticity–velocity formulation of the governing equations is shown to produce valid results when used in conjunction with the LRR gridding technique. Furthermore, each LRR grid is used to form a nonuniform equivalent tensor product (ETP) grid and also, in most cases, an equispaced fully refined (FR) grid; these additional grids are supersets of the LRR grids and thus contain refinement in exactly the same regions. Performance comparisons between the LRR, ETP, and FR grids indicate that the LRR method provides substantial savings in execution time and computer memory requirements, without compromising solution accuracy.