Viscosity and thermalization

I solve the relativistic Navier stokes equations assuming boost invariance and rotational symmetry. I compare the resulting numerical solutions for two limiting models of the shear viscosity. In the first model the shear viscosity is made proportional to the temperature. Thus, η ∝ T/σ₀ where σ₀ is some fixed cross section (perhaps σ₀ ∼ λ_QCD^-2). This viscosity model is typical of the classical Boltzmann simulations of Gyulassy and Molnar. In the second model the shear viscosity is made proportional to T³. This model is typical of high temperature QCD. When the initial mean free path of the T³ model is four times larger than the T/σ₀ model, the two models of viscosity produce the same radial flow. This result can be understood with simple scaling arguments. Thus, the large transport opacity needed in classical Boltzmann simulations is in part an artefact of the fixed scale σ₀ in these models.