Undriven, incompressible Kolmogorov flow in two dimensional doubly periodic strongly coupled dusty plasma is modelled using generalised hydrodynamics, both in linear and nonlinear regime. A complete stability diagram is obtained for low Reynolds numbers R and for a range of viscoelastic relaxation time tau(m) [0 < tau(m) < 10]. For the system size considered, using a linear stability analysis, similar to Navier Stokes fluid (tau(m) = 0), it is found that for Reynolds number beyond a critical R, say R-c, the Kolmogorov flow becomes unstable. Importantly, it is found that R-c is strongly reduced for increasing values of tau(m). A critical tau(c)(m) is found above which Kolmogorov flow is unconditionally unstable and becomes independent of Reynolds number. For R < R-c, the neutral stability regime found in Navier Stokes fluid (tau(m) = 0) is now found to be a damped regime in viscoelastic fluids, thus changing the fundamental nature of transition of Kolmogorov flow as function of Reynolds number R. A new parallelized nonlinear pseudo spectral code has been developed and is benchmarked against eigen values for Kolmogorov flow obtained from linear analysis. Nonlinear states obtained from the pseudo spectral code exhibit cyclicity and pattern formation in vorticity and viscoelastic oscillations in energy. (C) 2014 AIP Publishing LLC.