A stream function implicit finite difference scheme for 2D incompressible flows of Newtonian fluids