We examine linear perturbation theory to evaluate the contribution of viscosity coefficient in the growing of dark matter perturbations in the context of the bulk viscous dark energy model inspired by thermodynamical dissipative phenomena proposed by [B. Mostaghel et al., 2017, The European Physical Journal C, 77, 541]. As the cosmological implementations, we investigate the {\it Integrated Sachs-Wolfe} (ISW) auto-power spectrum, the ISW-galaxy cross-power spectrum and derive limits on $f\sigma_8$. The dimensionless bulk viscosity coefficient ($\gamma$) in the absence of interaction between dark sectors, modifies the Hubble parameter and the growth function, while the Poisson equation remains unchanged. Increasing $\gamma$ reduces the dark matter growing mode at the early epoch while a considerable enhancement will be achieved at the late time. This behavior imposes non-monotonic variation in the time evolution of gravitational potential generating a unique signature on the CMB photons. The bulk viscous dark energy model leads to almost a decreasing in ISW source function at the late time. Implementation of the Redshift Space Distortion (RSD) observations based on "Gold-2017" catalogue, shows $\Omega^0_{\rm m} =0.303^{+0.044}_{-0.038}$, $\gamma=0.033^{+0.098}_{-0.033}$ and $\sigma_8= 0.769^{+0.080}_{-0.089}$ at $1\sigma$ level of confidence. Finally, tension in the $\sigma_8$ is alleviated in our viscous dark energy model.

Inspired by thermodynamical dissipative phenomena, we consider bulk viscosity for dark fluid in a spatially flat two-component Universe. Our viscous dark energy model represents Phantom crossing avoiding Big-Rip singularity. We propose a non-minimal derivative coupling scalar field with zero potential leading to accelerated expansion of Universe in the framework of bulk viscous dark energy model. In this approach, coupling constant ($\kappa$) is related to viscosity coefficient ($\gamma$) and energy density of dark energy at the present time ($\Omega_{\rm DE}^0$). This coupling is bounded as $\kappa\in [-1/9H_0^2(1-\Omega_{\rm DE}^0), 0]$ and for $\gamma=0$ leads to $\kappa=0$. To perform robust analysis, we implement recent observational data sets including Joint Light-curve Analysis (JLA) for SNIa, Gamma Ray Bursts (GRBs) for most luminous astrophysical objects at high redshifts, Baryon Acoustic Oscillations (BAO) from different surveys, Hubble parameter from HST project, {\it Planck} data for CMB power spectrum and CMB Lensing. Joint analysis of JLA$+$GRBs$+$BAO$+$HST shows that $\Omega_{\rm DE}^0=0.696\pm 0.010$, $\gamma=0.1404\pm0.0014$ and $H_0=68.1\pm1.3$ at $1\sigma$ confidence interval. {\it Planck} TT observation provides $\gamma=0.32^{+0.31}_{-0.26}$ at $68\%$ confidence limit for viscosity coefficient. Tension in Hubble parameter is alleviated in this model. Cosmographic distance ratio indicates that current observed data prefer to increase bulk viscosity. Finally, the competition between Phantom and Quintessence behavior of viscous dark energy model can accommodate cosmological old objects reported as a sign of age crisis in $\Lambda$CDM model.