where \(q_k\)  is the heat flux and \(-I=-I_0=-\hat{\gamma }T\)  the energy
density dissipation rate in the absence of any additional forces, i.e. with
potential energy \(\phi =0\) .
For the explicit formulation of the coefficients see
[1]}, [2]}, [3]}, [4]}.
