High-order bicompact schemes for numerical modelling of multispecies multi-reaction gas flows

Euler equations for multidimensional inviscid gas flows with multispecies and multi-reaction are considered. Using the Marchuk–Strang splitting method, an implicit numerical scheme for this system is constructed. Its convection part is computed by the bicompact scheme SDIRK3B4 of fourth order in space and third order in time, while its chemical part is computed by the $L$-stable Runge-Kutta method of second order. The SDIRK3B4 scheme is compared to the WENO5/SR scheme in case of one- and two-dimensional flows with detonation waves. It is shown, that the SDIRK3B4 scheme has the same factual accuracy as the WENO5/SR, but the former needs 20–40 times less time steps and does not require any special algorithms to suppress non-physical breakdown of detonation waves on relatively coarse meshes.