We consider a generalization of the Lemaitre-Tolman-Bondi (LTB) solutions by keeping the LTB metric but replacing its dust matter source by an imperfect fluid with anisotropic pressure Πab. Assuming that total matter-energy density ρ is the sum of a rest mass term, ρ(m), plus a radiation ρ(r) = 3p density where p is the isotropic pressure, Einstein's equations are fully integrated without having to place any previous assumption on the form of Πab. Three particular cases of interest are contained: the usual LTB dust solutions (the dust limit), a class of FLRW cosmologies (the homogeneous limit) and of the Vaydia solution (the vacuum limit). Initial conditions are provided in terms of suitable averages and contrast functions of the initial densities of ρ(m), ρ(r) and the 3-dimensional Ricci scalar along an arbitrary initial surface t = ti. We consider the source of the models as an interactive radiation-matter mixture in local thermal equilibrium that must be consistent with causal Extended Irreversible Thermodynamics (hence Πab is shear viscosity). Assuming near equilibrium conditions associated with small initial density and curvature contrasts, the evolution of the models is qualitatively similar to that of adiabatic perturbations on a matter plus radiation FLRW background. We show that initial conditions exist that lead to thermodynamically consistent models, but only for the full transport equation of Extended Irreversible Thermodynamics. These interactive mixtures provide a reasonable approximation to a dissipative ‘tight coupling’ characteristic of radiation-matter mixtures in the radiative pre-decoupling era.
