This paper deals with the construction of a class of high order accurate Residual Dis- tribution schemes for the Navier Stokes equations using conformal meshes. The approx- imation of the solution is obtained using standard Lagrangian finite elements, and the total residual of the problem is constructed taking into account both the advective and the diffusive terms in order to discretize within the same scheme both parts of the gov- erning equation. To cope with the fact that the normal component of the gradients of the numerical solution is discontinuous across the faces of the elements, the gradient of the numerical solution is recovered at each degree of freedom of the grid and then interpolated with the same shape functions used for the solution. The procedure is fully described for the scalar case, and formaly extended to the system case. Linear and non-linear schemes are constructed and their accuracy is first tested with the help of manufactured solutions, and then applied to several (2D and 3D) test cases.