This is devoted to numerical simulation of time domain inverse acoustic scattering problems with a point-like scatterer, multiple point-like scatterers or normal size scatterers. Using Green’s functions and time convolution, we propose direct sampling methods to reconstruct the location of the scatterer. The methods involve only integral calculus without solving any equations and are easy to implement. Numerical experiments show the effectiveness and robustness of the methods.