Streamer front propagation involves steep gradients in charge density and electric field. Since the front has to be meshed with a sufficient number of points, adaptive meshing is essential for fast and accurate numerical simulations. In this paper a one-dimensional (1D) moving mesh method recently developed by Tang and Tang (2003 SIAM J. Numer. Anal. 41 487–515) is successfully applied to the simulation of streamer discharge. One-and-half (1.5D) and two-dimensional (2D) simulations of streamer discharges in nitrogen at atmospheric pressure are presented. The moving mesh method is combined with the third order ULTIMATE QUICKEST scheme (Leonard 1991 Comput. Math. Appl. Mech. Eng. 88 17–74) to solve the advection part of the plasma continuity ...