Abstract: This paper proposes a numerical method for handling the inviscid reactive Euler equations, which can capture the correct detonation wave structure with coarse grids and achieve numerical simulation of multi-component and multi chemical reactions problems. To calculate velocity of the precise shock wave, the level set method is applied to track the detonation front and capture a sharp shock discontinuity. By using an operator splitting technique, we can separate the stiff reaction source terms from the governing equations. A modified finite volume method is used to discretize the conservation law equations, and solve the problem of a loss of interface conservation caused by the level set method. Through numerical tests for multi-component and multi chemical reactions stiff detonation waves, our algorithm shows good reliability and robustness.