where q = qn is the heat flux vector with the unit vector n and T0 is
the equilibrium temperature. s is the total reaction time which is
defined by the mean reaction time s0 multiplied by the particle
number of the system, and k is the total mean free path which is defined
by the mean free path k0 multiplied by the particle number of
the system