Transport processes in anisotropic gravitational collapse
  In this paper we introduce a new method to study the influence of thermal conduction and viscous processes in anisotropic gravitational collapse. To this end we employ the HJR method to solve the Einstein equations. The Maxwell-Cattaneo type transport equations are used to find the temperature and bulk and shear viscous pressures. Under some conditions Maxwell-Cattaneo transport equations comply with relativistic causality. Thus, it is advisable to use them instead of Eckart transport equations. In the inner layers of the star the temperature ceases to be sensitive to the boundary condition. This behavior, which can be explained in terms of the Eddington approximation, allows us to find the thickness of the neutrinosphere. The dynamics of collapsing dense stars is deeply influenced by the neutrino emission/absorption processes. These cool the star and drive it to a new equilibrium state. Therefore, the calculation of transport coefficients is based on these processes. 
Justino Martinez