We test the physical relevance of the full and the truncated versions of the Israel-Stewart theory of irreversible thermodynamics in a cosmological setting. Using a dynamical systems method, we determine the asymptotic future of plane symmetric Bianchi type I spacetimes with a viscous mathematical fluid, keeping track of the magnitude of the relative dissipative fluxes, which determines the applicability of the Israel-Stewart theory. We consider the situations where the dissipative mechanisms of shear and bulk viscosity are involved separately and simultaneously. It is demonstrated that the only case in the given model when the fluid asymptotically approaches local thermal equilibrium, and the underlying assumptions of the Israel-Stewart theory are therefore not violated, is that of a dissipative fluid with vanishing bulk viscosity. The truncated Israel-Stewart equations for shear viscosity are found to produce solutions which manifest pathological dynamical features and, in addition, to be strongly sensitive to the choice of initial conditions. Since these features are observed already in the case of an oversimplified mathematical fluid model, we have no reason to assume that the truncation of the Israel-Stewart transport equations will produce relevant results for physically more realistic fluids. The possible role of bulk and shear viscosity in cosmological evolution is also discussed.

We apply the causal Israel-Stewart theory of irreversible thermodynamics to model the matter content of the universe as a dissipative fluid with bulk and shear viscosity. Along with the full transport equations we consider their widely used truncated version. By implementing a dynamical systems approach to Bianchi type IV and V cosmological models with and without cosmological constant, we determine the future asymptotic states of such universes and show that the truncated Israel-Stewart theory leads to solutions essentially different from the full theory. The solutions of the truncated theory may also manifest unphysical properties. Finally, we find that in the full theory shear viscosity can give a substantial rise to dissipative fluxes, driving the fluid extremely far from equilibrium, where the linear Israel-Stewart theory ceases to be valid.

We investigate the three types of class B Bianchi cosmologies filled with a tilted perfect fluid undergoing velocity diffusion in a scalar field background. We consider the two most importantcases: dust and radiation. A complete numerical integration of the Einstein field equations coupled with the diffusion equations is done to demonstrate how the presence of diffusion can affect the dynamics of cosmological evolution, where the most attention is paid to changes to the late-time behaviour. We show that aside from quantitative effects, diffusion can result in significant qualitative differences. For example, the cosmologies may recollapse if diffusion is sufficiently strong, or evolve towards the de Sitter state otherwise. In constrast to the diffusionless case, radiation isotropizes in presence of diffusion, and the tilt decreases exponentially at later times: $V\sim e^{-0.25\tau}$; also, we determine the decay rates of energy density, which become slower when the diffusion term is non-zero.

We apply the dynamical systems approach to ever-expanding Bianchi type VIII cosmologies filled with a tilted $\gamma$-fluid undergoing velocity diffusion on a scalar field. We determine the future attractors and investigate the late-time behaviour of the models. We find that at late times the normalized energy density $\Omega$ tends to zero, while the scalar potential $\Phi$ approaches 1 and dominates the evolution. Moreover, we demonstrate that in presence of diffusion fluids with $\gamma<3/2$, which includes physically important cases of dust $(\gamma=1)$ and radiation $(\gamma=4/3)$, are asymptotically non-tilted; the velocity of the fluid with $\gamma=3/2$ tends to a constant value $0<\bar{V}<1$; and stiffer fluids evolve towards a state of extreme tilt. Finally, we show that diffusion significantly reduces the decay rates of energy density for dust and fluids stiffer than dust $(\gamma \geq 1)$; for example, at $\gamma=4/3$ (radiation) we obtain $\rho/H^2 \propto e^{-3H_0 t}$ at late times, while $\rho/H^2 \propto e^{-4H_0 t}$ when diffusion is absent.

We investigate a simple inhomogeneous anisotropic cosmology (plane symmetric $G_2$ model) filled with a tilted perfect fluid undergoing velocity diffusion on a scalar field. Considered are two types of fluid: dust and radiation. We solve the system of Einstein field equations and diffusion equations numerically and demonstrate how the universe evolves towards its future asymptotic state. Also, typical time scales of characteristic processes are determined. The obtained results for dust- and radiation-filled cosmologies are compared to each other and to those in the diffusionless case, giving a hint on which effects can be the result of including diffusion in more complicated models. For example, diffusion causes the accelerated expansion stage to arrive at later times.