Numerical resolution of pulsating detonation waves

The canonical problem of the one-dimensional, pulsating, overdriven detonation wave has been studied for over 30 years, not only for its phenomenological relation to the evolution of multidimensional detonation instabilities, but also to provide a robust, reactive, high-speed flowfield with which to test numerical schemes. The present study examines this flowfield using high-order, essentially non-oscillatory schemes, systematically varying the level of resolution of the reaction zone, the size and retention of information in the computational domain, the initial conditions, and the order of the scheme. It is found that there can be profound differences in peak pressures as well as in the period of oscillation, not only for cases in which the reaction front is under-resolved, but for cases in which the computation is corrupted due to a too-small computational domain. Methods for estimating the required size of the computational domain to reduce costs while avoiding erroneous solutions are proposed and tested.