In our project we discuss a finite element time-relaxation method for high Reynolds number flows. The key idea consists of using local projections on polynomials defined on macro element of each pair of two elements sharing a face. We give the formulation for the scalar convection–diffusion equation and a numerical illustration.