The method of intrinsic low dimensional manifolds (ILDM) is applied to gas phase kinetic equations associated with combustion of solid energetic materials. Solution to common models is hindered by a necessity to track the evolution of dozens of species over times scales which vary by a factor as much as 1011. The ILDM method is able to systematically reduce the number of species which must be tracked, reduce the stiffness of the governing equations, and maintain fidelity to the full solution at user-specified coarse time scales. This comes at the expense of losing information at fine time scales and developing large tables which describe the evolution of species on slow time scales. The method is applied to a common model for the combustion of the explosive HMX.