A numerical method for DNS/LES of turbulent reacting flows

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.