  The renormalization group method is applied to the study of homogeneous and
flat Friedmann-Robertson-Walker type Universes, filled with a causal bulk
viscous cosmological fluid. The starting point of the study is the
consideration of the scaling properties of the gravitational field equations,
of the causal evolution equation of the bulk viscous pressure and of the
equations of state. The requirement of scale invariance imposes strong
constraints on the temporal evolution of the bulk viscosity coefficient,
temperature and relaxation time, thus leading to the possibility of obtaining
the bulk viscosity coefficient-energy density dependence. For a cosmological
model with bulk viscosity coefficient proportional to the Hubble parameter, we
perform the analysis of the renormalization group flow around the scale
invariant fixed point, therefore obtaining the long time behavior of the scale
factor.
