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".