Spatially varying primordial magnetic fields may be efficiently dissipated prior to the epoch of recombination due to the large viscosity of the baryon-photon fluid. We show that this dissipation may result in observable chemical potential mu and Compton y distortions in the cosmic microwave background radiation (CMB) spectrum. Current upper limits on mu and y from FIRAS constrain magnetic fields to have strength B_0 < 3\times 10^{-8}Gauss (scaled to the present) between comoving coherence length \approx 400 pc and \approx 0.6 Mpc. This represents the strongest upper limit on small-scale primordial magnetic fields to date.

We examine the evolution of magnetic fields in an expanding fluid composed of matter and radiation with particular interest in the evolution of cosmic magnetic fields. We derive the propagation velocities and damping rates for relativistic and non-relativistic fast and slow magnetosonic, and Alfv\'en waves in the presence of viscous and heat conducting processes. The analysis covers all MHD modes in the radiation diffusion and the free-streaming regimes. When our results are applied to the evolution of magnetic fields in the early universe, we find that cosmic magnetic fields are damped from prior to the epoch of neutrino decoupling up to recombination. Our findings have multifold implications for cosmology. The dissipation of magnetic field energy into heat during the epoch of neutrino decoupling ensures that most magnetic field configurations generated in the very early universe satisfy big bang nucleosynthesis constraints. Further dissipation before recombination constrains models in which primordial magnetic fields give rise to galactic magnetic fields or density perturbations. Finally, the survival of Alfv\'en and slow magnetosonic modes on scales well below the Silk mass may be of significance for the formation of structure on small scales (abridged).