This work is licensed under the Creative Commons Attribution 4.0 International License.
Introduction
The classical irreversible thermodynamic with internal variables (CIT-IV) means a general approach to the study of different interactions among irreversible thermodynamic processes using classical exstensive variables. After the pioneering works of Lars Onsager [1, 2], many excellent scientists: Eckart [3,4,5,6], Prigogine [7], Groot et al. [8], Meixner et al. [9], Colemann et al. [10], Truesdell [11], Kluitenberg [12,13,14,15,16], Maugin et al. [17, 18] developed the basics of the theory CIT-IV characterized by the introduction of thermodynamic variables, called internal (or hidden) variables.
The flexibility of the methodologies used in CIT-IV consists in the fact that “a priori” the physical meaning of the internal thermodynamical variables is not specified but only their influence on particular types of occurring phenomena in the considered medium is assumed. For this reason, in the following discussion, we will call these thermodynamical variables: dynamical variables since they have also been used, successfully, for the study of problems in dielectric and magnetic relaxation phenomena [19,20,21,22], diffusion and anelastic deformation [23,24,25,26,27,28,29,30,31,32]. It’s important to note that the many theoretical results have been confirmed by the experimental data [33,34,35,36,37,38,39,40].
Parallel approaches were worked out in the frame of Extended Irreversible Thermodynamics (EIT) [41] and in Rational Extended Thermodynamics (RET) [42, 43] which, in consideration of suggestions from the kinetic theory, considers an entropy depending on the fluxes as well as on the classical variables. It was shown [29,30,31] that by using the usual procedure of CIV-IV it is possible to describe the relaxation of thermal phenomena thus obtaining results which are generally justified by additional hypothesis suggested by kinetic theory.
In this paper, we show that thermoconduction phenomena in visco-anelastic media may be studied by using the systematic procedure of CIT-IV. In order to study the behavior of rheological media, in the presence of viscous and anelastic phenomena, besides determining the stress-strain relationships it is important to obtain the heat propagation equation. For this purpose the tensor
\varepsilon _{\alpha \beta }^{(1)}
(inelastic strain) and the vector ξ are introduced as internal variables which charaterize both stress-strain relationships and which generalize the Fourier heat equation (parabolic type equation) [44] and Maxwell-Cattaneo-Vernotte (hyperbolic type equation) [45,46,47].
In Section 2, introducing a vectorial dynamical variable j and a stress field
\tau _{\alpha \beta }^{(eq)}
which is of a thermoelastic nature, an explicit form for the entropy production is derived.
In Sections 3–7, the phenomenological equations, which generalize the Kelvin, Jeffreys and Poynting-Thomson bodies are obtained.
In Section 8, by virtue of the dynamical variable, influencing the thermal transport phenomena, the heat equation is obtained. In particular, in the isotropic case, when the medium has symmetry properties that, under orthogonal transformations, are invariant with respect to all rotations and inversions of the frame of axes, it is obtained that the heat flux can be split in two parts: a first contribution J(0), governed by Fourier law and a second contribution J(1), obeying the Maxwell-Cattaneo-Vernotte equation (MCV) [45,46,47,48] in which a relaxation time is present.
Finally, in Section 9, we obtain a general temperature equation which generalizes the analogous equations of Fourier and MCV.
The balance equation of entropy
In the contest of irreversible processes an important role is played by the flow of heat which, classically, is not considered to be a state variable. Therefore, we will suppose that the specific entropy s, depends not only on the specific internal energy u and the total strain ɛαβ, the tensor
\varepsilon _{\alpha \beta }^{(1)}
describing the inelastic strain and also on a vectorial dynamic variable, ξ that is an odd function of microscopic particles velocities that have influence on the propagation phenomena which occur in the medium.
s = s(u,{\varepsilon _{\alpha \beta }},\varepsilon _{\alpha \beta }^{(1)},{\xi _\alpha }){\kern 1pt} ,
where ξα (α = 1,2,3) is the α-component of the vector ξ.
Theorem 1 (Gibbs relation)
By using(1)we can obtain the following Gibbs relationTds = du - \nu {\kern 1pt} \tau _{\alpha \beta }^{(eq)}d{\varepsilon _{\alpha \beta }} + \nu {\kern 1pt} \tau _{\alpha \beta }^{(1)}d\varepsilon _{\alpha \beta }^{(1)} + \nu {\kern 1pt} {j_\alpha }{\kern 1pt} d{\xi _\alpha }{\kern 1pt} ,where ν is the specific volume,\tau _{\alpha \beta }^{(eq)}is the equilibrium stress tensor,\tau _{\alpha \beta }^{(1)}is the affinity-stress conjugate to\varepsilon _{\alpha \beta }^{(1)}
, T is the absolute temperature and jα is the α-component of the vector j conjugate of dξα /dt.
Proof
We define the absolute temperature{T^{ - 1}}{\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} \frac{\partial }{{\partial u}}s(u,{\varepsilon _{\alpha \beta }},\varepsilon _{\alpha \beta }^{(1)},{\xi _\alpha }){\kern 1pt} ,
the equilibrium-stress tensor\tau _{\alpha \beta }^{(eq)}{\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} - \rho T\frac{\partial }{{\partial {\varepsilon _{\alpha \beta }}}}s(u,{\varepsilon _{\alpha \beta }},\varepsilon _{\alpha \beta }^{(1)},{\xi _\alpha }){\kern 1pt} ,
the affinity-stress conjugate to
\varepsilon _{\alpha \beta }^{(1)}\tau _{\alpha \beta }^{(1)}{\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} \rho T\frac{\partial }{{\partial \varepsilon _{\alpha \beta }^{(1)}}}s(u,{\varepsilon _{\alpha \beta }},\varepsilon _{\alpha \beta }^{(1)},{\xi _\alpha }){\kern 1pt} ,
and the vector j conjugate to the internal vector variable ξ{j_\alpha }{\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} \rho T\frac{\partial }{{\partial {\xi _\alpha }}}s(u,{\varepsilon _{\alpha \beta }},\varepsilon _{\alpha \beta }^{(1)},{\xi _\alpha }){\kern 1pt} ,
where
\rho {\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} {\nu ^{ - 1}}
is the mass density. By using equations (3)–(6) from (2), we obtain the differential ds of s:
Tds = du - \nu {\kern 1pt} \tau _{\alpha \beta }^{(eq)}d{\varepsilon _{\alpha \beta }} + \nu {\kern 1pt} \tau _{\alpha \beta }^{(1)}{\kern 1pt} d\varepsilon _{\alpha \beta }^{(1)} + \nu {\kern 1pt} {j_\alpha }{\kern 1pt} d{\xi _\alpha }{\kern 1pt} .
The relation (7) is usually called Gibbs relation, in which the usual summation convection for dummy is used.
From (7), we have:
\rho T\frac{{ds}}{{dt}} = \rho \frac{{du}}{{dt}} - {\kern 1pt} \tau _{\alpha \beta }^{(eq)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}} + \tau _{\alpha \beta }^{(1)}\frac{{d\varepsilon _{\alpha \beta }^{(1)}}}{{dt}} + {j_\alpha }{\kern 1pt} \frac{{d{\xi _\alpha }}}{{dt}}{\kern 1pt} ,
where
\frac{d}{{dt}} = \frac{\partial }{{\partial t}} + {\bf v} \cdot \nabla ,
is the substantial derivative respect to time and v is the velocity field and ∇ is the gradient.
To analyze phenomena due to viscous flows (analogous to those which occur during flows in ordinary viscous liquids and gases) we introduce the following viscous stress tensor
\tau _{\alpha \beta }^{(vi)}{\kern 1pt} \mathop = \limits^{\rm def} {\kern 1pt} {\tau _{\alpha \beta }} - \tau _{\alpha \beta }^{(eq)}{\kern 1pt} ,
where ταβ is the mechanical stress tensor which occurs in the equation of motion
\rho \frac{{{d{\rm v}_\alpha }}}{{dt}} = \rho {\kern 1pt} {F_\alpha } + \frac{{\partial {\tau _{\alpha \beta }}}}{{\partial {x^\beta }}},
and in the first law of thermodynamics
\rho \frac{{du}}{{dt}} = - \nabla \cdot {{\bf J}^{(q)}} + {\tau _{\alpha \beta }}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}}.
In (11) the force Fα is the volume force per unit of mass and in (12) the vector J(q) is the heat flux.
Theorem 2 (Entropy production)
By using the first law of thermodynamics(12)the balance equation of the entropy can be obtained\rho \frac{{ds}}{{dt}} = - \nabla \cdot \left( {\frac{{{{\boldsymbol J}^{(q)}}}}{T}} \right) + {\sigma ^{(s)}},and the entropy production{\sigma ^{(s)}} = {T^{ - 1}}\left[ {{{\boldsymbol J}^{(q)}} \cdot \left( { - {T^{ - 1}}\nabla {\kern 1pt} T} \right) + \tau _{\alpha \beta }^{(vi)}{\kern 1pt} \frac{{d{\kern 1pt} {\varepsilon _{\alpha \beta }}}}{{dt}}{\kern 1pt} + {\kern 1pt} \tau _{\alpha \beta }^{(1)}{\kern 1pt} \frac{{d\varepsilon _{\alpha \beta }^{(1)}}}{{dt}}{\kern 1pt} + {\kern 1pt} {j_\alpha }{\kern 1pt} \frac{{d{\xi _\alpha }}}{{dt}}} \right] \ge 0.
Proof
The equation (12), by virtue of (10) becomes
\rho \frac{{du}}{{dt}} = - \nabla \cdot {{\bf J}^{(q)}} + \left( {\tau _{\alpha \beta }^{(eq)} + \tau _{\alpha \beta }^{(vi)}} \right)\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}},
which substituted in the equation (8) gives
\rho {\kern 1pt} \frac{{ds}}{{dt}} = {T^{ - 1}}{\kern 1pt} \left( { - \nabla \cdot {{\bf J}^{(q)}} + \tau _{\alpha \beta }^{(vi)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}} + \tau _{\alpha \beta }^{(1)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}} + {j_\alpha }\frac{{d{\xi _\alpha }}}{{dt}}} \right){\kern 1pt} .
Using the following identity
{T^{ - 1}}\nabla \cdot {{\bf J}^{(q)}} = \nabla \cdot \left( {\frac{{{{\bf J}^{(q)}}}}{T}} \right) + {T^{ - 2}}{{\bf J}^{(q)}} \cdot \nabla T,
the equation (16) gives (13) (balance equation of the entropy) and (14) (entropy production of entropy).
It is observed that the entropy production is due to three types of phenomena: the first term on the right-hand of (14) gives the contribution of the heat conduction phenomena, the second sum is the contribution of viscous phenomena, the last sum is the contribution of the variation of the dynamical variable. Each term is a production of flux (J(q),
\tau _{\alpha \beta }^{(vi)}
,
d\varepsilon _{\alpha \beta }^{(1)}/dt
, jα of the process and affinities conjugated to them: T−1∇T, d ɛαβ /dt,\tau _{\alpha \beta }^{(1)}
, dξα /dt, respectively.
Phenomenological equations
According to the usual procedure of non-equilibrium thermodynamics, by virtue of the form (14) for the entropy production, as anisotropic media, we have the following phenomenological equations:
{J}_\alpha ^{(q)} = L_{\alpha \beta }^{(q)(q)}{\kern 1pt} \left( { - {T^{ - 1}}\frac{{\partial T}}{{\partial {x^\beta }}}} \right) + L_{\alpha (\mu \nu )}^{(q)(0)}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}}{\kern 1pt} + L_{\alpha (\mu \nu )}^{(q)(1)}\tau _{\mu \nu }^{(1)}{\kern 1pt} + {\kern 1pt} L_{\alpha \mu }^{(q)(\xi )}{\kern 1pt} \frac{{d{\xi _\mu }}}{{dt}}{\kern 1pt} ,\tau _{\alpha \beta }^{(vi)} = L_{(\alpha \beta )\mu }^{(0)(q)}{\kern 1pt} \left( { - {T^{ - 1}}\frac{{\partial T}}{{\partial {x^\mu }}}} \right) + L_{(\alpha \beta )(\mu \nu )}^{(0)(0)}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}}{\kern 1pt} + L_{(\alpha \beta )(\mu \nu )}^{(0)(1)}\tau _{\mu \nu }^{(1)} + {\kern 1pt} L_{(\alpha \beta )\mu }^{(0)(\xi )}{\kern 1pt} \frac{{d{\xi _\mu }}}{{dt}}{\kern 1pt} ,\frac{{d\varepsilon _{\alpha \beta }^{(1)}}}{{dt}} = L_{(\alpha \beta )\mu }^{(1)(q)}{\kern 1pt} \left( { - {T^{ - 1}}\frac{{\partial T}}{{\partial {x^\mu }}}} \right) + L_{(\alpha \beta )(\mu \nu )}^{(1)(0)}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}}{\kern 1pt} + L_{(\alpha \beta )(\mu \nu )}^{(1)(1)}\tau _{\mu \nu }^{(1)} + {\kern 1pt} L_{(\alpha \beta )\mu }^{(1)(\xi )}{\kern 1pt} \frac{{d{\xi _\mu }}}{{dt}}{\kern 1pt} ,{j_\alpha } = L_{\alpha \beta }^{(\xi )(q)}{\kern 1pt} \left( { - {T^{ - 1}}\frac{{\partial T}}{{\partial {x^\beta }}}} \right) + L_{\alpha (\mu \nu )}^{(\xi )(0)}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}}{\kern 1pt} + L_{\alpha (\mu \nu )}^{(\xi )(1)}{\kern 1pt} \tau _{\mu \nu }^{(\xi )(1)}{\kern 1pt} + L_{\alpha \beta }^{(\xi )(\xi )}{\kern 1pt} \frac{{d{\xi _\beta }}}{{dt}}{\kern 1pt} .
The tensors L are called phenomenological tensors and the indices of these tensors enclosed in round brackets mean that they are symmetrical because the tensors ɛμν,
\varepsilon _{\alpha \beta }^{(1)}
,
\tau _{\alpha \beta }^{(1)}
and
\tau _{\alpha \beta }^{(vi)}
are symmetric. The first of these equations may be regarded as a generalization of Fourier’s law. Equation (19) describes the viscous flow phenomenon and it may be considered to be a generalization of Stockes-Navier’s law. Finally, the equations (20) and (21) are the phenomenological equations for the irreversible process of the dynamic degrees of freedom.
Substituting (18)–(21) into (14) one has
\begin{array}{*{20}{l}}{T{\sigma ^{(s)}} = }&{L_{\alpha \beta }^{(q)(q)}\left( {{T^{ - 2}}\frac{{\partial T}}{{\partial {x^\alpha }}}\frac{{\partial T}}{{\partial {x^\beta }}}} \right) + \left[ {L_{\alpha \beta }^{(\xi )(q)} + L_{\alpha \beta }^{(q)(\xi )}} \right]\frac{{d{\xi _\alpha }}}{{dt}}\left( { - {T^{ - 1}}\frac{{\partial T}}{{\partial {x^\beta }}}} \right) + }\\{}&{ + L_{(\alpha \beta )(\mu \nu )}^{(0)(0)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}} + }\\{}&{ + \left[ {L_{(\alpha \beta )(\mu \nu )}^{(1)(0)} + L_{(\mu \nu )(\alpha \beta )}^{(0)(1)}} \right]{\kern 1pt} \tau _{\alpha \beta }^{(1)}{\kern 1pt} \frac{{d{\varepsilon _{\mu \nu }}}}{{dt}} + }\\{}&{ + L_{(\alpha \beta )(\mu \nu )}^{(1)(1)}\tau _{\alpha \beta }^{(1)}{\kern 1pt} \tau _{\mu \nu }^{(1)} + }\\{}&{ + L_{\alpha \beta }^{(\xi )(\xi )}\frac{{d{\xi _\alpha }}}{{dt}}{\kern 1pt} \frac{{d{\xi _\beta }}}{{dt}}{\kern 1pt} .}\end{array}
This form for the entropy production (s(s) ≥ 0) is a useful detail if isotropic media are considered in Section 5.
Symmetric relations
Since the time derivatives dɛαβ /dt,
d\varepsilon _{\alpha \beta }^{(1)}/dt
, the heat flow J(q) and j vector conjugate to the internal variable ξα are odd functions of the microscopic particle velocities and the stress
\tau _{\alpha \beta }^{(vi)}
,
\tau _{\alpha \beta }^{(1)}
, dξα /dt and the temperature gradient ∂T /∂t are even functions of these velocities, the Onsager-Casimir reciprocity relation reads [8]:
L_{\alpha \beta }^{(q)(q)} = L_{\beta \alpha }^{(q)(q)}{\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} L_{\alpha \beta }^{(\xi )(q)} = L_{\beta \alpha }^{(\xi )(q)}{\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} L_{\alpha \beta }^{(q)(\xi )} = L_{\beta \alpha }^{(q)(\xi )}{\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} L_{\alpha \beta }^{(\xi )(\xi )} = L_{\beta \alpha }^{(\xi )(\xi )}{\kern 1pt} ,L_{(\alpha \beta )(\mu \nu )}^{(0)(0)} = L_{(\mu \nu )(\alpha \beta )}^{(0)(0)}{\kern 1pt} {\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} {\kern 1pt} L_{(\alpha \beta )(\mu \nu )}^{(1)(1)} = L_{(\mu \nu )(\alpha \beta )}^{(1)(1)}{\kern 1pt} ,L_{(\alpha \beta )(\mu \nu )}^{(1)(0)} = L_{(\mu \nu )(\alpha \beta )}^{(1)(0)}{\kern 1pt} {\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} {\kern 1pt} L_{(\alpha \beta )(\mu \nu )}^{(0)(1)} = L_{(\mu \nu )(\alpha \beta )}^{(0)(1)}{\kern 1pt} ,L_{(\alpha \beta )(\mu \nu )}^{(1)(0)} = - L_{(\mu \nu )(\alpha \beta )}^{(0)(1)}{\kern 1pt} .
The equations (23)–(26) reduce the number of indipendent components of the phenomenological tensors.
Phenomenological equations for isotropic media
In this section, we consider perfect isotropic media for which the symmetry properties are invariant under orthogonal transformations with respect to all rotations and to inversion of the frame of axes.
In this case, it can be shown [8, 49] that the phenomenological tensors of odd order vanish while even tensors can have the following forms:
L_{\alpha \beta }^{(q)(q)} = T{\lambda ^{(q,q)}}{\delta _{\alpha \beta }}{\kern 1pt} \quad ,\quad L_{\alpha \beta }^{(q)(\xi )} = {\lambda ^{(q,\xi )}}{\kern 1pt} {\delta _{\alpha \beta }}{\kern 1pt} ,L_{\alpha \beta }^{(\xi )(q)} = T{\lambda ^{(\xi ,q)}}{\delta _{\alpha \beta }}{\kern 1pt} \quad ,\quad {L^{(\xi )(\xi )}} = {\lambda ^{(\xi ,\xi )}}{\delta _{\alpha \beta }}{\kern 1pt} ,L_{(\alpha \beta )(\mu \nu )}^{(i)(k)} = \frac{1}{2}\eta _s^{(i,k)}({\delta _{\alpha \mu }}{\delta _{\beta \nu }} + {\delta _{\alpha \nu }}{\delta _{\beta \mu }}) + \frac{1}{3}(\eta _v^{(i,k)} - \eta _s^{(i.k)}){\delta _{\alpha \beta }}{\delta _{\mu \nu }}{\kern 1pt} \quad (i,k = 0,1).
By using (27)–(29), the equations (18)–(21) become
{J}_\alpha ^{(q)} = - {\lambda ^{(q,q)}}\frac{{\partial T}}{{\partial {x^\alpha }}} + {\lambda ^{(q,\xi )}}{\kern 1pt} \frac{{d{\xi _\alpha }}}{{dt}}{\kern 1pt} ,\tau _{\alpha \beta }^{(vi)} = \eta _s^{(0,0)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}} + (\eta _v^{(0,0)} - \eta _s^{(0,0)}){\kern 1pt} \frac{{d\varepsilon }}{{dt}}{\kern 1pt} {\delta _{\alpha \beta }} + \eta _s^{(0,1)}\tau _{\alpha \beta }^{(1)} + (\eta _v^{(0,1)} - \eta _s^{(0,1)}){\kern 1pt} {\tau ^{(1)}}{\kern 1pt} {\delta _{\alpha \beta }}{\kern 1pt} ,\frac{{d\varepsilon _{\alpha \beta }^{(1)}}}{{dt}} = \eta _s^{(1,0)}\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}} + (\eta _v^{(1,0)} - \eta _s^{(1,0)}){\kern 1pt} \frac{{d\varepsilon }}{{dt}}{\kern 1pt} {\delta _{\alpha \beta }} + \eta _s^{(1,1)}\tau _{\alpha \beta }^{(1)} + (\eta _v^{(1,1)} - \eta _s^{(1,1)}){\kern 1pt} {\tau ^{(1)}}{\kern 1pt} {\delta _{\alpha \beta }}{\kern 1pt} ,{j_\alpha } = - {\lambda ^{(\xi ,q)}}{\kern 1pt} \frac{{\partial T}}{{\partial {x^\alpha }}} + {\lambda ^{(\xi ,\xi )}}{\kern 1pt} \frac{{d{\xi _\alpha }}}{{dt}}{\kern 1pt} .
The scalar quantities λ,
\eta _s^{(0,0)}
and
\eta _v^{(0,0)}
are called phenomenological coefficients. In particular
\eta _s^{(0,0)}
and
\eta _v^{(0,0)}
may be called the shear viscosity and the volume viscosity, respectively. These coefficients also occur in the theory of ordinary (Stokes-Navier) viscous fluids. Equation (22), by virtue of (30)–(33), becomes
\begin{array}{*{20}{l}}{T{\kern 1pt} {\sigma ^{(s)}} = }&{{\kern 1pt} {\lambda ^{(q,q)}}{\kern 1pt} {{(\frac{{\partial T}}{{\partial {x^\alpha }}})}^2} + \eta _s^{(0,0)}{{(\frac{{d{\varepsilon _{\alpha \beta }}}}{{dt}})}^2}{\kern 1pt} + 3(\eta _v^{(0,0)} - \eta _s^{(0,0)}){{(\frac{{d\varepsilon }}{{dt}})}^2} + }\\{}&{ + \eta _s^{(1,1)}{{({\tau ^{(1)}})}^2}{\kern 1pt} + 3(\eta _v^{(1,1)} - \eta _s^{(1,1)}){{({\tau ^{(1)}})}^2} + {\lambda ^{(\xi ,\xi )}}{{(\frac{{d{\xi _\alpha }}}{{dt}})}^2}{\kern 1pt} .}\end{array}
By virtue of the positive definite character of the entropy production, several inequalities for the phenomenological coefficients may be derived.
For example, we have
{\lambda ^{(q,q)}} \ge 0\quad ,\quad \eta _s^{(0,0)} \ge 0\quad ,\quad \eta _v^{(0,0)} \ge \eta _s^{(0,0)}{\kern 1pt} ,\eta _s^{(1,1)} \ge 0\quad ,\quad \eta _v^{(1,1)} \ge \eta _s^{(1,1)}\quad ,\quad {\lambda ^{(\xi ,\xi )}} \ge 0{\kern 1pt} .
The coefficients λ(q,ξ) and λ(ξ,q) represent possible cross effects between the heat flux and the irreversible process of the dynamic degrees of freedom.
Linear equations of state
Let f the specific free energy of the medium
f = u - T{\kern 1pt} s.
With the aid of Gibbs relation (7), we have:
df = \nu {\kern 1pt} \tau _{\alpha \beta }^{(eq)}d{\varepsilon _{\alpha \beta }} - \nu \tau _{\alpha \beta }^{(1)}d\varepsilon _{\alpha \beta }^{(1)} - \nu {\kern 1pt} {j_\alpha }\;d{\xi _\alpha } - sdT{\kern 1pt} ,
where
\nu {\kern 1pt} \tau _{\alpha \beta }^{(eq)} = \frac{{\partial f}}{{\partial {\varepsilon _{\alpha \beta }}}},\quad {\kern 1pt} \nu {\kern 1pt} \tau _{\alpha \beta }^{(1)} = - \frac{{\partial f}}{{\partial \varepsilon _{\alpha \beta }^{(1)}}},\quad \nu {\kern 1pt} {j_\alpha } = - {\kern 1pt} \frac{{\partial f}}{{\partial {\xi _\alpha }}},\quad s = - {\kern 1pt} \frac{{\partial f}}{{\partial T}}.
Let us choose a configuration Σ(0) with uniform temperature T0 and in which s0 is the specific entropy, ν0 is the specific volume and
{(\tau _{\alpha \beta }^{(eq)})_0}
,
{(\tau _{\alpha \beta }^{(1)})_0}
,
{({\tau _{\alpha \beta }})_0}
,
{({\varepsilon _{\alpha \beta }})_0}
,
{(\varepsilon _{\alpha \beta }^{(1)})_0}
and (jα)0 are zero.
Let us suppose that the deviation from Σ(0) is sufficiently small (with ρ ≈ ρ0) and the free energy for isotropic medium can be written in the following form
\begin{array}{*{20}{l}}{f = {\kern 1pt} {\nu _0}{\kern 1pt} }&{\left\{ {\frac{1}{2}\left[ {a_{\alpha \beta \mu \nu }^{(0)(0)}{\kern 1pt} {\varepsilon _{\alpha \beta }}{\kern 1pt} {\varepsilon _{\mu \nu }}{\kern 1pt} + {\kern 1pt} a_{\alpha \beta \mu \nu }^{(1)(1)}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} \varepsilon _{\mu \nu }^{(1)} + a_{\alpha \beta }^{(\xi )(\xi )}{\xi _\alpha }{\kern 1pt} {\xi _\beta }} \right]} \right.}\\{}&{ + a_{\alpha \beta \mu \nu }^{(0)(1)}{\kern 1pt} {\varepsilon _{\alpha \beta }}{\kern 1pt} \varepsilon _{\mu \nu }^{(1)}{\kern 1pt} + {\kern 1pt} a_{\alpha \beta \mu }^{(0)(\xi )}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} {\xi _\mu }{\kern 1pt} + {\kern 1pt} a_{\alpha \beta \mu }^{(1)(\xi )}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} {\xi _\mu }}\\{}&{\left. { + (T - {T_0})\left[ {a_{\alpha \beta }^{(0)(T)}{\varepsilon _{\alpha \beta }}{\kern 1pt} + {\kern 1pt} a_{\alpha \beta }^{(1)(T)}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} + {\kern 1pt} a_\alpha ^{(\xi )(T)}{\xi _\alpha }} \right]} \right\} - \Psi (T){\kern 1pt} .}\end{array}
where Ψ is some function of the temperature. In (40), the tensors “a” and the vector
a_\alpha ^{(\xi )(T)}
, are constants (i.e. they do not depend on temperature and the strains) and are determined by the physical properties of the medium in the reference state.
From (39)4 we have
s{\kern 1pt} = {\kern 1pt} - {\kern 1pt} {\nu _0}\left( {a_{\alpha \beta }^{(0)(T)}{\varepsilon _{\alpha \beta }} + a_{\alpha \beta }^{(1)(T)}\varepsilon _{\alpha \beta }^{(1)} + a_\alpha ^{(\xi )(T)}{\kern 1pt} {\xi _\alpha }} \right){\kern 1pt} + {\kern 1pt} \frac{{d{\kern 1pt} \Psi }}{{dT}}{\kern 1pt} ,
and as u = f + T s we obtain
\begin{array}{*{20}{l}}{u{\kern 1pt} = {\kern 1pt} }&{{\nu _0}\left\{ {\frac{1}{2}\left[ {a_{\alpha \beta \mu \nu }^{(0)(0)}{\kern 1pt} {\varepsilon _{\alpha \beta }}{\kern 1pt} {\varepsilon _{\mu \nu }}{\kern 1pt} + {\kern 1pt} a_{\alpha \beta \mu \nu }^{(1)(1)}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} \varepsilon _{\mu \nu }^{(1)}{\kern 1pt} + {\kern 1pt} a_{\alpha \beta }^{(\xi )(\xi )}{\kern 1pt} {\xi _\alpha }{\kern 1pt} {\xi _\beta }} \right]} \right.}\\{\kern 1pt} &{\left. { - {T_0}\left[ {a_{\alpha \beta }^{(0)(T)}{\varepsilon _{\alpha \beta }} + a_{\alpha \beta }^{(1)(T)}{\kern 1pt} \varepsilon _{\alpha \beta }^{(1)}{\kern 1pt} + {\kern 1pt} a_\alpha ^{(\xi )(T)}{\kern 1pt} {\xi _\alpha }} \right]} \right\}{\kern 1pt} + {\kern 1pt} T{\kern 1pt} \frac{{d\Psi }}{{dT}}{\kern 1pt} - {\kern 1pt} \Psi (T){\kern 1pt} .}\end{array}
The specific heat at constant deformation, c(ɛ), may be defined by
{c_{(\varepsilon )}}{\kern 1pt} = {\kern 1pt} \frac{{\partial {\kern 1pt} u}}{{\partial T}}{\kern 1pt} ,
and from (42) it obtains
{c_{(\varepsilon )}}{\kern 1pt} = {\kern 1pt} T{\kern 1pt} \frac{{{d^2}{\kern 1pt} \Psi }}{{d{T^2}}}{\kern 1pt} .
By integrating the equation (44), one has
\Psi {\kern 1pt} = {\kern 1pt} {c_{(\varepsilon )}}T{\kern 1pt} \log {\kern 1pt} {\kern 1pt} \left( {\frac{T}{{{T_0}}}} \right){\kern 1pt} + {\kern 1pt} {s_0},T - {c_{(\varepsilon )}}(T - {T_0}) - {u_0}{\kern 1pt} .
If c(ɛ) is constant where s0 and u0 are integration constants that are the specific entropy and the specific energy in the reference state, respectively.
Using the phenomenological equations (Section 5) and the linear equations of state (Section 6) we obtain the following results:
Theorem 3 (Decomposition of the heat flow)
By using equation(56)and(33)the heat flow J(q)can be split in two parts{{\boldsymbol J}^{(q)}} = {{\boldsymbol J}^{(0)}}{\kern 1pt} + {\kern 1pt} {{\boldsymbol J}^{(1)}}{\kern 1pt} ,where J(0)satisfies the Fourier’s law and J(1)obeys the Maxwell-Cattaneo-Vernotte law.
Proof
Assuming λ(ξ)(ξ) ≠ 0, by virtue of equation (56), from the equation (33) we obtain
\frac{{d{\xi _\alpha }}}{{dt}} = - \frac{{{a^{(\xi ,\xi )}}}}{{{\lambda ^{(\xi ,\xi )}}}}{\kern 1pt} {\xi _\alpha }{\kern 1pt} + {\kern 1pt} \frac{{{\lambda ^{(\xi ,q)}}}}{{{\lambda ^{(\xi ,\xi )}}}}\frac{{\partial T}}{{\partial {x^\alpha }}}{\kern 1pt} .
Substituting the equation (89) into equation (30), one obtains
{{\bf J}^{(q)}} = {{\bf J}^{(0)}}{\kern 1pt} + {\kern 1pt} {{\bf J}^{(1)}}{\kern 1pt} ,
where
{{\bf J}^{(0)}} = {\kern 1pt} - {\lambda ^{(0)}}{\kern 1pt} \nabla T{\kern 1pt} ,
with
{\lambda ^{(0)}}{\kern 1pt} = {\kern 1pt} {\lambda ^{(q,q)}}{\kern 1pt} - {\kern 1pt} \frac{{{\lambda ^{(q,\xi )}}{\kern 1pt} {\lambda ^{(\xi ,q)}}}}{{{\lambda ^{(\xi ,\xi )}}}}{\kern 1pt} ,
and
{{\bf J}^{(1)}} = - \frac{{{\lambda ^{(q,\xi )}}{\kern 1pt} {a^{(\xi ,\xi )}}}}{{{\lambda ^{(\xi ,\xi )}}}}{\kern 1pt} {\boldsymbol \xi} {\kern 1pt} .
By virtue of equation (89), we see that J(1) satisfies the following equation:
{t_r}{\kern 1pt} \frac{{d{{\bf J}^{(1)}}}}{{dt}} + {{\bf J}^{(1)}}{\kern 1pt} = {\kern 1pt} - \frac{{{\lambda ^{(\xi ,q)}}{\kern 1pt} {\lambda ^{(q,\xi )}}}}{{{\lambda ^{(\xi ,\xi )}}}}{\kern 1pt} \nabla T{\kern 1pt} ,
where
{t_r}{\kern 1pt} = {\kern 1pt} \frac{{{\lambda ^{(\xi ,\xi )}}}}{{{a^{(\xi ,\xi )}}}}{\kern 1pt} ,
is the relaxation time.
From (90) we can see that the heat current is split into two parts: the first one, J(0), is governed by Fourier’s law (91) while the second part, J(1), is governed by the MCV type constitutive equation (94).
Temperature equation
We consider an undeformable medium at rest, i.e.
\frac{d}{{dt}} = \frac{\partial }{{\partial t}},
and we assume that the specific internal energy u is related to the temperature by
du = {c_{(v)}}dT{\kern 1pt} ,
being c(v) the heat capacity per unit of mass at constant volume.
Theorem 4 (Generalization of temperature equation)
By using(91)and the first law of thermodynamics we obtain the following temperature equation{t_r}\frac{{{\partial ^2}T}}{{\partial {t^2}}} + \frac{{\partial T}}{{\partial t}} = \alpha '\Delta T{\kern 1pt} + {\kern 1pt} {t_r}{\kern 1pt} \eta '\frac{\partial }{{\partial t}}\Delta T,where α′ and η′ are parameters expressed through the phenomenological and state coefficients.
Proof
By virtue of (96) and (97) the first law of thermostatic (12) becomes
\rho {\kern 1pt} {c_{(v)}}{\kern 1pt} \frac{{\partial T}}{{\partial t}}{\kern 1pt} + {\kern 1pt} \nabla \cdot ({{\bf J}^{(0)}} + {{\bf J}^{(1)}}) = 0,
and from (99) we obtain
\rho {\kern 1pt} {c_{(v)}}{\kern 1pt} {t_r}{\kern 1pt} \frac{{{\partial ^2}T}}{{\partial {t^2}}}{\kern 1pt} + {\kern 1pt} {t_r}{\kern 1pt} \nabla \cdot \left[ {\frac{{\partial {{\bf J}^{(0)}}}}{{\partial t}} + \frac{{\partial {{\bf J}^{(1)}}}}{{\partial t}}} \right] = 0{\kern 1pt} .
Summing (99) and (100), by using (91) and (94), we have the following temperature equation
{t_r}{\kern 1pt} \frac{{{\partial ^2}T}}{{\partial {t^2}}}{\kern 1pt} + {\kern 1pt} \frac{{\partial T}}{{\partial t}} = {\kern 1pt} \alpha '{\kern 1pt} \Delta T{\kern 1pt} + {\kern 1pt} {t_r}{\kern 1pt} \eta '{\kern 1pt} \frac{{\partial \Delta T}}{{\partial t}},
where
\left\{ {\begin{array}{*{20}{l}}{\alpha '{\kern 1pt} = {\kern 1pt} \frac{{{\lambda ^{(0)}}{\lambda ^{(\xi ,\xi )}} + {\lambda ^{(\xi ,q)}}{\lambda ^{(q,\xi )}}}}{{\rho {\kern 1pt} {c_{(v)}}{\kern 1pt} {\lambda ^{(\xi ,\xi )}}}}}\\{\eta ' = \frac{{{\lambda ^{(0)}}}}{{\rho {\kern 1pt} {c_{(v)}}}}{\kern 1pt} .}\end{array}} \right.
The equation (101) generalizes the temperature equation with analogous equations of Fourier and Maxwell-Cattaneo-Vernotte.
Special cases
If the relaxation time tr(95) can be neglected with respect to the heat propagation time t (tr ≪ t) the equation (101) becomes
\frac{{\partial T}}{{\partial t}} = \alpha '\Delta T{\kern 1pt} .
i.e. the Fourier equation.
• If the phenomenological coefficients satisfy the following condition
{\lambda ^{(q,q)}} = \frac{{{\lambda ^{(q,\xi )}},{\lambda ^{(\xi ,q)}}}}{{{\lambda ^{(\xi ,\xi )}}}}{\kern 1pt} ,
from (92) one has: λ(0) = 0 and the equations (102) become
\left\{ {\begin{array}{*{20}{l}}{\alpha '{\kern 1pt} = {\kern 1pt} \frac{{{\lambda ^{(\xi ,q)}}{\lambda ^{(q,\xi )}}}}{{\rho {\kern 1pt} {c_{(v)}}{\kern 1pt} {\lambda ^{(\xi ,\xi )}}}}}\\{\eta ' = 0}\end{array}} \right.{\kern 1pt} ,
and the equation (101) becomes
{t_r}{\kern 1pt} \frac{{{\partial ^2}T}}{{\partial {t^2}}}{\kern 1pt} + {\kern 1pt} \frac{{\partial T}}{{\partial t}}{\kern 1pt} - {\kern 1pt} \alpha '{\kern 1pt} \Delta T{\kern 1pt} = {\kern 1pt} 0{\kern 1pt} ,
being the Maxwell-Cattaneo-Vernotte.
Remarks and conclusion
The Fourier equation (103) is a parabolic equation and therefore the heat flux can be considered of infinite speed. Several scientists view this result as a paradox. In 1992, Fichera [51] states that this claim is unfounded as Fourier’s theory is not correctly interpreted. In our opinion the Fourier result is justified by the value of the phenomenological parameters of the medium and therefore it cannot be said that Fourier’s theory leads to a paradox.
The Maxwell-Cattaneo-Vernotte law leads to a hyperbolic equation obtained with intuitive hypotheses but without justifying it in the context of a general theory. Various authors, for instance [41,42,43], have proposed several theories suggested by kinetic theory in order to justify the Maxwell-Cattaneo-Vernotte equation by assuming that entropy depends on fluxes instead of state variables.
In this paper, without a priori hypotheses and the important role of phenomenological and state coefficients, we obtained the stress-strain relations for viscoanelastic media and an equation for the temperature which generalizes the Fourier and MCV temperature equations.
Declarations
Conflict of interest
There is no conflict of interests regarding the publication of this paper.
Funding
There is no funding regarding the publication of this paper.
Author’s contribution
V.C.-Writing original draft, Writing review editing, Methodology.
Acknowledgement
The author deeply appreciates the anonymous reviewers for their helpful and constructive suggestions, which can help further improve this paper.
Data availability statement
All data that support the findings of this study are included within the article.
Using of AI tools
The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.