This work is devoted to the superconvergence in space approximation of a fully discrete scheme for the incompressible time-dependent Navier-Stokes Equations in three-dimensional domains. We discrete by Inf-Sup-stable Finite Element in space and by a semi-implicit backward Euler (linear) scheme in time.
Using an extension of the duality argument in negative-norm for elliptic linear problems (see for instance [1]) to the mixed velocity-pressure formulation of the Stokes problem, we prove some superconvergence in space results for the velocity with respect to the energy-norm, and for a weaker norm of L2(0, T; L2(Ω)) type (this latter holds only for the case of Taylor-Hood approximation). On the other hand, we also obtain optimal error estimates for the pressure without imposing constraints on the time and spatial discrete parameters, arriving at superconvergence in the H1 (Ω)-norm again for Taylor-Hood approximations. These results are numerically verified by several computational experiments, where two splitting in time schemes are also considered.