at high temperatures [@karsch07] from recent lattice QCD calculations [@Cheng:2007jq], as $$\label{13}
P = \omega \rho,\hspace*{1cm}T = \beta \rho^r,\hspace*{1cm}\xi = \alpha \rho + \frac{9}{\omega_0} T_c^4,$$ with $\omega = (\gamma-1)$, $\gamma \simeq 1.183$, $r\simeq 0.213$, $\beta\simeq 0.718$, $$\alpha = \frac{1}{9\omega_0}  \frac{9\gamma^2-24\gamma+16}{\gamma-1},$$ and $\omega_0 \simeq 0.5-1.5$ GeV. In the following we assume that $\alpha \rho >> 9/\omega_0 T_c^4$, and therefore we take $\xi \simeq \alpha \rho$. In order to close the system of the cosmological equations, we have also to give the expression of the relaxation time $\tau $, for which we adopt the expression [@Ma95], $$\label{tau}
\tau=\xi\rho^{-1}\simeq\alpha .$$

Eqs. (\[13\]) are standard in the study of the viscous cosmological models, whereas the equation for $\tau$ is a simple procedure to ensure that the speed of viscous pulses does not exceed the speed of light. Eq. (\[tau\]) implies that the relaxation time in our treatment is constant but strongly depends on EoS. These equations are without sufficient thermodynamical motivation, but in the absence of better alternatives, we shall follow the practice of adopting them in the hope that they will at least provide some indication of the range of bulk viscous effects. The temperature law is the simplest law guaranteeing positive heat 