Bulk viscosity and cavitation in boost-invariant hydrodynamic expansion

We solve second order relativistic hydrodynamics equations for a boost-invariant 1 + 1-dimensional expanding fluid with an equation of state taken from lattice calculations of the thermodynamics of strongly coupled quark-gluon plasma. We investigate the dependence of the energy density as a function of proper time on the values of the shear viscosity eta, the bulk viscosity zeta, and second order coefficients, confirming that large changes in the values of the latter have negligible effects. Varying the shear viscosity between zero and a few times s/4 pi, with s the entropy density, has significant effects, as expected based on other studies. Introducing a nonzero bulk viscosity also has significant effects. In fact, if the bulk viscosity peaks near the crossover temperature T-c to the degree indicated by recent lattice calculations in QCD without quarks, it can make the fluid cavitate - falling apart into droplets. It is interesting to see a hydrodynamic calculation predicting its own breakdown, via cavitation, at the temperatures where hadronization is thought to occur in ultrarelativistic heavy ion collisions.