In this paper we recall a stabilization technique for finite element methods for convection-diffusion-reaction equations, originally proposed by Douglas and Dupont. The method uses least square stabilization of the gradient jumps a across element boundaries. We prove that the method is stable in the hyperbolic limit and prove optimal a priori error estimates. We address the question of monotonicity of discrete solutions and present some numerical examples illustrating the theoretical results