After reading a lot of papers on beyond standard cosmology recently, I'm becoming confused and a bit nervous about my own knowledge on physics! By definition, a perfect fluid is described by an energy-momentum tensor that has the general following form (I'm using the $\eta = \mathrm{diag}(1, -1, -1, -1)$ convention, and $c = 1$): $$\tag{1} T_{ab} = (\rho + p) \, u_a \, u_b - p \, \eta_{ab}. $$ In principle, the isotropic pressure $p$ could be any non-linear function of the energy density $\rho$ (the perfect fluid may be self-interacting). A scalar field is a good example of this.
But then, fluids with bulk viscosity $\zeta$ (no shear viscosity, no heath flux) have the same energy-momentum expression as (1), but add a new term to the isotropic pressure: $$\tag{2} p \quad \Rightarrow \quad \tilde{p} = p - \zeta \Theta, $$ where $\Theta = \nabla_{\mu} u^{\mu}$ is the expansion scalar. In isotropic/homogeneous FLRW cosmology, we have $\Theta = 3 \dot{a} / a$.
Fluids with bulk viscosity are generally described as non-perfect, since they generate entropy, while perfect fluid have constant entropy. But then, this is where I'm getting confused! What is a "perfect fluid" in general relativity? Is it simply a fluid that has (1) as energy-momentum tensor (so this would include all the fluids with bulk viscosity as "perfect fluids"), or do we have to add a constraint to (1)? If so, how can we translate that constraint, mathematically, to get a "perfect fluid"?
And is it true that all perfect fluids have constant entropy?
This question is related (but not the same) as this old one from mine: