Any solution of the Navier–Stokes equations in a three-dimensional axisymmetric domain admits a Fourier expansion with respect to the angular variable, and it can be noted that each Fourier coefficient satisfies a variational problem on the meridian domain, all problems being coupled due to the nonlinear convection term. We propose a discretization of these equations which combines Fourier truncation and finite element methods applied to each two-dimensional system. We perform the a priori and a posteriori analysis of this discretization.