In this study, proper orthogonal decomposition based reduced-order modelling and variational multiscale stabilization method are utilized for the incompressible Navier-Stokes equations. Also, the difficulties resulting from nonlinearity are eliminated by using the extrapolation. Theoretical analysis of the method is presented. To check the efficiency of the proposed method, we utilize the test problem of 2D channel flow past a cylinder.