A comparative study of difference schemes for the solution of $$\frac{{\partial u_i }}{{\partial t}} + u_j \frac{{\partial u_i }}{{\partial x_j }} = \nu \frac{{\partial ^2 u_i }}{{\partial x_j \partial x_j }} + \frac{{\partial F}}{{\partial x_i }}$$