A spatially non-dissipative, implicit numerical method to simulate turbulent reacting flows over a range of Mach numbers, is described. The compressible Navier-Stokes equations are rescaled so that the zero Mach number equations are discretely recovered in the limit of zero Mach number. The dependent variables are co-located in space, and thermodynamic variables are staggered from velocity in time. The algorithm discretely conserves kinetic energy in the incompressible, inviscid, non-reacting limit. The chemical source terms are implicit in time to allow for stiff chemical mechanisms. The algorithm is readily extended to complex chemical mechanisms. Numerical examples using both simple and complex chemical mechanisms are presented.
Bibliographical noteFunding Information:
This work was supported by the Air Force Office of Scientific Research under Grant FA9550-04-1-0341 and the Department of Energy under the Stanford ASC alliance. Computing resources were provided by the Minnesota Supercomputing Institute, the San Diego Supercomputing Center and the National Center for Supercomputing Applications.
- Turbulent reacting flows