This paper studies fully discrete approximations to the evolutionary Navier{
Stokes equations by means of inf-sup stable H1-conforming mixed nite elements
with a grad-div type stabilization and the Euler incremental projection method in
time. We get error bounds where the constants do not depend on negative powers
of the viscosity. We get the optimal rate of convergence in time of the projection
method. For the spatial error we get a bound O(hk) for the L2 error of the velocity,
k being the degree of the polynomials in the velocity approximation. We prove
numerically that this bound is sharp for this method.
