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