A numerical technique for solving unsteady non-Newtonian free surface flows