A highly accurate numerical shock-fitting scheme composed of fifth order spatial and third order temporal discretizations was applied to the two-dimensional reactive Euler equations for several models with single-step irreversible kinetics in a slab geometry. The first case examined was for a calorically perfect mixture of ideal gases with a prescribed material interface motion. The steady detonation phase speed that results from this scheme was compared with the equivalent from shock-capturing. It was found that, the errors while using shock-fitting at the coarser solutions were significantly reduced in comparison with shock-capturing. Furthermore, for an error on the order of 10 m/s, which is similar to that observed in experiments, shock-fitting offers a computational savings on the order of 1000 Additionally, the behavior of the detonation phase speed was examined for several slab widths while utilizing the Wescott-Stewart-Davis (WSD) model, which is commonly used high explosive (HE) modeling, to evaluate the detonation performance of PBX 9501. It was found that the thickness effect curve resulting from this equation of state and reaction model using published values is dramatically more steep than observed in recent experiments. However, adjusting the reaction rate parameters improves the correlation of the model to experimental results. Moreover, similar results are predicted when the shock-fitting scheme is used to examine rate stick behavior of the steady detonation phase speed in PBX 9501.
展开▼