We propose to develop a computer code for simulating Pulse Detonation Engines (PDEs). This code will be based on the space-time conservation element-solution element method. The basic idea of this method is global and local flux conservation in a space-time domain. In this method conservation variables and their derivatives are treated as unknowns which are defined on the solution element. A time marching scheme is constructed using local conservation within each conservation element. It has been demonstrated by that this method yields high resolution for shocks, rarefaction waves, acoustic waves, vortices, detonation waves, and shock/acoustic,shock/vortices interactions. During Phase I of this project, we propose to develop robust computer codes, based on the above methodology, for unsteady, implicit one dimensional and quasi one dimensional PDE problems, with multistep chemistry. These codes will be extended to solve general combustion problems, validated with benchmark problems and compared with MOZART code.At the end of Phase II, PDE and combustion software packages will be available for axisymmetric and three dimensional geometries using structured and unstructured grids with full multistep chemistry.