High order finite difference methods with subcell resolution for stiff multispecies detonation capturing

Authors

Wei Wang, Florida International UniversityFollow
Chi-Wang Shu, Brown UniversityFollow
Helen C. Yee, NASA Ames Research CenterFollow
Dmitry V. Kotov, Center for Turbulence ResearchFollow
Bjorn Sjögreen, Lawrence Livermore National LaboratoriesFollow

Date of this Version

2-2015

Citation

Published in Communications in Computational Physics, Volume 17, Issue 2 (February 2015), pp. 317-336. doi: 10.4208/cicp.250214.130814a

Comments

U.S. government work.

Abstract

In this paper, we extend the high order finite-difference method with subcell resolution (SR) in [33] for two-species stiff one-reaction models to multispecies and multireaction chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.