Left panel shows the evolution of temperature in the one-dimensional thermal conductivity test computed with the scheme #1 (explicit evolution of the parabolic term) and with the scheme #2 (hyperbolic equation for heat flux) with three different time steps. Right panel shows a comparison between simulations using scheme #2, one including the saturation factor of Eq. (11) and another one without it.