In this article, we are proposing an updated form of Krylov subspace-based interpolatory projection techniques for the stabilization of incompressible Navier-Stokes flows. In the proposed techniques, we utilize the reduced-order modelling approach implicitly, where reduced-order matrices need not be gathered explicitly. To estimate the aimed optimal feedback matrix, only the factored solution of the desired continuous-time algebraic Riccati equation (CARE) needs to be stored through the classical eigenvalue decomposition. The sparse structure of the target systems will remain invariant within the matrix-vector operations for generating the bases of the projector matrices through Krylov subspace techniques, where a cohesive projection scheme will be incorporated to ensure the potency of the projector matrices. So, the proposed techniques will be feasible for memory allocation and enhance the rapid convergence of the simulation. Analysis of the target systems’ transient characteristics, such as eigenvalues and step-responses, will be used to ascertain the competence and reliability of the proposed techniques. Necessary computation will be done numerically through MATLAB. Stabilization of the transient behaviors and minimization of the simulation time are the prime concerns in this work. Eventually, by comparing with the contemporary techniques the advancement of the proposed techniques will be confirmed.