“…with the subscript V, E, L denoting the homogeneous viscous part, Ericksen part for the static state, and Leslie part for the non-equilibrium state, respectively. Now we introduce a strictly concave mathematical entropy function η = ρs(ν, u, d, ∇d, C, K, l, h), (67) where ν = 1/ρ, (C, K) are tensors with the same size of (σ, π), and (l, h) are vectors with the same size of (g, q). (C, K, l) are used to describe the viscous-elastic effects of nematic liquid crystal flows, and h characterizes the heat conduction induced by temperature gradients.…”