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 t...