A simple and efficient numerical model is presented for the simulation of pulse combustors. It is based on the numerical solution of the quasi-1D unsteady flow equations and on phenomenological sub-models of turbulence and combustion. The gas dynamics equations are solved by using the Flux Difference Splitting (FDS) technique, a finite-volume upwind numerical scheme, and ENO reconstructions to obtain secondorder accurate non-oscillatory solutions. The numerical fluxes computed at the cell interfaces are used to transport also the reacting species, their formation energy and the turbulent kinetic energy. The combustion progress in each cell is evaluated explicitly at the end of each time step according to a second-order overall reaction kinetics. In this way, the computations of gas dynamic evolution and heat release are decoupled, which makes the model particularly simple and efficient. A comprehensive set of measurements has been performed on a small Helmholtz type pulse-jet in order to validate the model. Air and fuel consumptions, wall temperatures, pressure cycles in both combustion chamber and tail-pipe, and instantaneous thrust have been recorded in different operating conditions of the device. The comparison between numerical and experimental results turns out to be satisfactory in all the working conditions of the pulse-jet. In particular, accurate predictions are obtained of the device operating frequency and of shape, amplitude and phase of the pressure waves in both combustion chamber and tail-pipe.