With the attempts of extending the hydrodynamic framework of heavy-ion collision to proton-proton and other small and low-energy systems, we are confronted with the question of how small the system can get and still be safely modeled as a fluid. One of the transport coefficients required in the second-order relativistic viscous hydrodynamics is the shear relaxation time, inclusion of which solves the causality violation problem in the Navier-Stokes equation. In phenomenological studies this coefficient has been taken as a constant and much attention has gone into finding and fixing the shear viscosity to entropy density ratio, . This transport coefficient also happens to control the nonhydrodynamic mode of the out-of-equilibrium hydrodynamics theory. It has been predicted that for decreasing system size, observables become sensitive to variation in shear relaxation time as a result of increasing dominance of the nonhydrodynamic mode, which could potentially indicate the breakdown of hydrodynamics. In this study, we try to test this prediction in peripheral Pb-Pb collisions at 2.76 TeV and Au-Au collisions at 200 GeV, with IPGlasma initial condition and -dimensional viscous hydrodynamics. We find that elliptic flow does show adequate sensitivity to variation in relaxation time for decreasing system size. The multiplicity rapidity density limit for applicability of hydrodynamics is found to be around , with the possibility of refinement in this value given a way to improve the centrality resolution in experimental data for referencing in peripheral collisions.