\documentclass[12pt]{article}
\usepackage{amsmath,amssymb,amsfonts}
\begin{document}
Inspired by thermodynamical dissipative phenomena, we consider constant 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 to describe 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.700^{+0.021}_{-0.023} $, $\gamma=0.1406^{+0.0064}_{-0.0065}$ and $H_0=69.1^{+3.5}_{-2.8}$ 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.
\end{document}