The reduction of the truncation error including numerical diffusion, has been one of the most important tasks in the development of numerical methods. The stream line method is used to cancel cross numerical diffusion and some of the non-diffusion type truncation error. The two-step stream line method which is the combination of the stream line method and finite difference methods is developed in this work for the solution of the governing equations of incompressible buoyant turbulent flow. This method is compared with the finite difference method. The predictions of both classes of numerical methods are compared with experimental findings. Truncation error analysis also has been performed in order to the compare truncation error of the stream line method with that of finite difference methods.