Subgrid Stabilized Projection Method for 2D Unsteady Flows at High Reynolds Numbers