A robust implicit compact scheme for two-dimensional unsteady flows with a biharmonic stream function formulation