We consider a pressure-stabilized, finite element approximation of incompressible flow problems in primitive velocity–pressure variables, which is based on a projection of the gradient of the discrete pressure onto the space of discrete functions. Equal order interpolation for the velocity and the pressure can be employed with this formulation. The method introduced here is specially developed to be used on anisotropic finite element meshes with large element aspect ratios.