Multilevel Adaptive Methods for Turbulent Combustion with Detailed Chemistry: Part I

C. Liu,C. Liao,Z. Liu,X. Zheng,S. McCormick

1995-01-01

Abstract:A highly accurate and efficient numercial method is developed for modeling reacting flows with detailed chemistry (e.g. 48 species and 22 finite- rate reation steps) in general curvilinear coordinates. The prediction of NOx formation is made by solving a fully-coupled system of flow and chemistry equations. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. We use a fully-implicit backward Euler technique and a third-order monotones upwind-biased scheme on a staggered grid for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line distributive relaxation os used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implictly. A preselected-direction least squares iterative method is developed for solving the species equations. A parallel version if the code (P4 form) is used to enhance efficiency on advanced architectures. To obtain good resolution for reacting flow around holes, a fine patch grid is embedded for local refinement, which is shown to be very efficient for improving the entire solution. Calculatons were performed for 3-D plane and axisymmetric co-flow diffusion flames and for strong swirling combustion in a 3-D gas turbine combusture. The current computational method is found to be two to three orders of magnitude faster than traditional iterative methods. For a 49x65x65x frid of 207,025 grid points, the total calculation takes about 200 time steps and 21.3 CRAY-YMP hours