Bulk viscosity in heavy ion collision

The effect of a temperature dependent bulk viscosity to entropy density ratio~($\zeta/s$) along with a constant shear viscosity to entropy density ratio~($\eta/s$) on the space time evolution of the fluid produced in high energy heavy ion collisions have been studied in a relativistic viscous hydrodynamics model. The boost invariant Israel-Stewart theory of causal relativistic viscous hydrodynamics is used to simulate the evolution of the fluid in 2 spatial and 1 temporal dimension. The dissipative correction to the freezeout distribution for bulk viscosity is calculated using Grad's fourteen moment method. From our simulation we show that the method is applicable only for $\zeta/s<0.004$.