  We study the evolution of a flat Friedmann-Robertson-Walker Universe, filled
with a bulk viscous cosmological fluid, in the presence of time varying
``constants''. The dimensional analysis of the model suggests a proportionality
between the bulk viscous pressure of the dissipative fluid and the energy
density. On using this assumption and with the choice of the standard equations
of state for the bulk viscosity coefficient, temperature and relaxation time,
the general solution of the field equations can be obtained, with all physical
parameters having a power-law time dependence. The symmetry analysis of this
model, performed by using Lie group techniques, confirms the unicity of the
solution for this functional form of the bulk viscous pressure. In order to
find another possible solution we relax the hypotheses assuming a concrete
functional dependence for the ``constants''.
