On Entropy Flux of Anisotropic Elastic Bodies

The framework of irreversible thermodynamics is fundamental in development of constitutive models. One of the important aspects of the extended irreversible thermodynamics is the relationship between the entropy flux and the heat flux, especially for phenomena far from equilibrium. In this paper, we demonstrate that the assumption that Lagrange multiplier conjugated to the energy balance equation (in the expression for the second law of thermodynamics) is a function of temperature Λε=Λε(θ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varLambda^{\varepsilon } = \varLambda^{\varepsilon } (\theta ) $$\end{document} is a sufficient condition to derive the entropy flux–heat flux relation for all isotropic materials as well as for a number of crystal classes including transverse isotropy, orthotropy, triclinic systems and rhombic systems. For all considered crystal classes, the entropy flux–heat flux relation was derived explicitly. Further, we demonstrate that for some crystal classes heat flux is nonzero even when temperature gradient vanishes (as stated by Eringen). The anisotropic functions, with respect to the symmetry groups of the crystal classes, were expressed in terms of isotropic functions. The proposed procedure is very general in the sense that it can be used with nonlinear constitutive relations as demonstrated here. The presented results confirm that the all crystal elastic bodies considered are hyperelastic.


Introduction
One of the important problems in thermodynamics is the relationship between the entropy flux and the heat flux for phenomena far from equilibrium.

B R. Vignjevic
The entropy principle based on the Clausius-Duhem inequality where ρ is density, η is the specific entropy density, q the heat flux and r the heat supply, has been widely adopted in the development of modern rational thermodynamics after the fundamental work of Coleman and Noll [1]. The main assumptions, motivated by the result of classical thermostatics, are that the entropy flux and the entropy supply s are proportional to the heat flux and the heat supply, respectively. Moreover, both constants of proportionality are assumed to be the reciprocal of the absolute temperature, i.e., These main assumptions, while tacit in the classical theory of continuum mechanics, do not hold particularly well for materials in general. In fact, it is known that they are inconsistent with the kinetic theory of ideal gases and are also found to be inadequate to account for the thermodynamics of diffusion. Further, Sellitto et al. [2] demonstrated that the proportionality relations (2) do not hold for nonlocal heat transfer at nanoscale and provide an illustrative example of cylindrical nanodevice connected to a graphene layer.
There is an extended formulation of the second law of thermodynamics which has been applied to nonequilibrium thermodynamics by Serrin [3] and Silhavy [4] and summarized by Truesdell and Bharatha [5]. See also Muschik [6], Müller [7] and Domínguez-Cascante and Jou [8]. A comprehensive review of the related literature and detailed derivations can be found in Lebon et al. [9] and Jou et al. [10,11].
The extended form of the second law, usually called the entropy inequality, seems to be the most general formulation of the continuous second law of thermodynamics proposed so far. In this theory, the assumptions given by Eq. (2) were abandoned and the entropy flux and the heat flux vector q are treated as independent constitutive quantities and hence leaving the entropy inequality in its general form ρη + div − ρs ≥ 0. ( Liu [12] proposed a method, reminiscent of the classical method of Lagrange multipliers, for expanding the inequality (3). Instead of this inequality restricting the solution of field equations, he considered solutions of an extended inequality which should hold for all fields. This can be done if one considers the field equations as constraints on solutions of the energy inequality.
Further, Liu [13] analyzed the thermodynamic theory of viscoelastic bodies and proved that for isotropic viscoelastic materials the results are identical to the classical results given by Eq. (2). In the same paper, he also proved that the body is hyperelastic because of the Lagrange multiplier Λ ε Λ(θ ).
However, for anisotropic elastic materials in general, the validity of the classical entropy flux relation is yet to be demonstrated.
The first contribution in this direction has been given by Liu [14,15], who proved by considering transversely isotropic elastic bodies and transversely isotropic rigid heat conductors that the classical entropy flux relation (2) need not be valid in general.
Three years later, Bargmann et al. [16] considered the energy influx-entropy influx relation in the Green-Naghdi type III theory of heat conduction and showed that the entropy influx and the energy influx are proportional via the absolute temperature if heat conduction is isotropic. Further, they demonstrated that influx proportionality cannot be postulated in general by giving counterexample of transversely isotropic conduction. Their proof is based on a representation formula for isotropic vectorvalued mappings of two vector arguments where Lagrange multiplier is assumed to have the following general form: Λ Λ (|u|, |v|, u · v). Using Green-Naghdi, set of internal variables S{α,α, ∇α, ∇α} where α is thermal displacement they showed, through an algebraic analysis, that the multiplier Λ Λ (α,α) is at least proportional to coldness (the inverse of absolute temperature and does not depend on the thermal displacement gradient ∇α nor on the temperature gradient ∇α) for isotropic conduction.
Podio-Guidugli [17] considered energy and entropy inflows in the theory of heat conduction and demonstrated proportional via the absolute temperature for isotropic conduction using similar procedure.
We do not see the possibility for the Bargmann et al. procedure to be generalized and extended to conduction in continuum with other symmetry types.
In this paper, we investigate the functional dependence of the Lagrange multiplier conjugated to the energy balance equation Λ ε Λ(θ ), irrespective of whether the classical entropy flux relation is valid. This enables derivation of the entropy flux-heat flux relations for a number of crystal classes including transverse isotropy, orthotropy, triclinic systems, monoclinic systems and rhombic systems.
The paper is organized as follows: In Sect. 2, the basic ideas and formulas typically used in this field are given as a starting point of our investigation. In Sect. 3, the entropy flux relation for viscoelastic bodies and transverse isotropic elastic materials is reconsidered assuming the Lagrange multiplier dependency on temperature, i.e., Λ ε Λ(θ ). In Sect. 4, the procedure introduced in Sect. 3 is extended for the derivation of the entropy flux relation of anisotropic elastic materials defined by the crystal classes listed above. In Sect. 5, the entropy inequality for anisotropic bodies is examined further. Conclusions related to the outcome of this work are given in Sect. 6.

The Entropy Principle
In this section, the basic framework of the entropy principle for viscoelastic materials is presented. The balance laws of mass, linear momentum and energy can be stated in current configuration as:ρ where T is the Cauchy stress tensor, b is external body force and r is external heat supply. Note that for solid bodies, it is more convenient to use a referential description. Also, since constitutive relations do not depend on external supplies, it suffices to consider only supply-free bodies. Consequently, the balance laws can be rewritten as: and the entropy inequality Here, the first Piola-Kirchhoff stress tensor T κ , the material heat flux vector q κ and the material entropy flux vector κ are related to the Cauchy stress tensor T, the heat flux vector q and the entropy flux vector by where F is the deformation gradient in referential coordinates and J |det F|. "Div" is the divergence operator with respect to the referential coordinates. It is well known that the entropy principle imposes severe restrictions on constitutive functions and the exploitation of such restrictions based on the Clausius-Duhem inequality and is relatively easy. For elastic materials, in general, the thermodynamic restrictions can be easily obtained by the well-known Coleman-Noll procedure [1].
The derivation of the relation between the entropy flux and the heat flux based on the entropy principle, referred to as entropy flux relation, is a typical problem in this new theory.
Here, we outline the consideration for isotropic viscoelastic materials with isotropic elasticity as a special case within the framework of continuum mechanics. The thermodynamics of the continua sufficiently close to equilibrium so that the principle of equipresence can be used. Consequently, the local constitutive relations for viscoelastic materials can be written as functions of the state variables i.e., whereḞ is time derivative of the deformation gradient, ε is the specific internal energy, g κ ∇θ is temperature gradient, θ is an empirical temperature, which is some convenient measure of the hotness (or coldness) of the thermodynamic state. Note that the density field ρ(X, t) is completely determined by the motion x(X, t) and the density ρ κ (X) in the reference configuration. Therefore, the thermodynamic process is defined as the solution of the field equations (the balance laws of the linear momentum and energy) and integration of the constitutive relations for T κ , q κ and ε.
The determination of the restrictions imposed on the constitutive functions by the entropy principle is one of the major objectives in modern continuum thermodynamics.

Method of Lagrange Multipliers
According to the entropy principle, there exist Lagrange multipliers Λ v conjugated to the momentum balance equation and Λ ε conjugated to the energy balance equation which depend on the state variables, such that the inequality is valid under no additional constraints, i.e., valid for any field x(X, t), θ(X, t). Further, we invoke the condition of material objectivity, which implies the following reduced constitutive equations for viscoelastic materials where C F T F is the right Cauchy-Green tensor.
First, note that (10) is linear with respect toẍ. Consequently, ρ κ Λ v the coefficient ofẍ must be equal to zero, i.e., Thus, (10) becomes Next, we consider terms iṅ Note that the term The other terms in the Eq. (15) have equivalent component forms. After substituting (14) and (15) into (13), by inspection, we conclude that this inequality is also linear with respect to the following derivatives θ ,ẍ, . g κ ,C, ∇g κ , ∇C, ∇Ċ .
As the inequality must hold for arbitrary fields, we have eliminated the constraints imposed by the field equations. The coefficients of the above derivatives must vanish identically. Otherwise, we could choose the fields in such a way that one negative term would dominate all others and the inequality would be violated. Hence, we obtain the following equations Then, the entropy inequality (13) reduces to Making use of second Piola-Kirchhoff tensor S κ F −1 T κ J F −1 TF −T and the right Cauchy-Green tensorĊ Ḟ T F + F TḞ 2 Ḟ T F sym (which are symmetric tensors), the inequality (18) can be written in a more compact form as: where σ is entropy production density. Moreover, from (16) we obtain which has the form of the thermostatic Gibbs relation.

Entropy Flux Relation for Viscoelastic Materials
For further evaluation of the consequences of the entropy principle, particularly in connection with relations (17), we invoke the material symmetry condition that has to be satisfied by {ε, T κ , q κ , κ , η} for isotropic viscoelastic bodies. For instance, this condition for heat flux can be expressed as: where O is the full orthogonal group. Note thatq κ is an isotropic vector-valued function of C,Ċ, θ, g κ and that (21) imposes restriction on its form. After a lengthy calculation starting from (15), Liu [13] proved that the following entropy flux relation holds Further, based on (22) and (17), I-Shih Liu concluded that Λ ε must be independent of C,Ċ and g κ . Thus, Accordingly, (20) becomes By comparison with the classical Gibbs relation in thermostatics, the function Λ ε can be identified as the reciprocal of the absolute temperature θ , i.e., which leads to the classical entropy flux relation (2). Consequently, one can come to the conclusion that, following I-Shih Liu's procedure for viscoelastic bodies outlined above, that for isotropic elastic materials with state variables (F, θ, g κ ) the relation (25) holds.

Entropy Flux of Anisotropic Elastic Materials
In this subsection, the derivation of the relation between the entropy flux and the heat flux for anisotropic materials with state variables (F, θ, g κ ) is considered. In this case, the requirements (16) and (17) reduce to (26) and (27), respectively.
In addition, the material symmetry condition for anisotropic elastic bodies has to be satisfied by {ε, T κ , q κ , κ , η}.
The first paper considering entropy flux for transversely isotropic elastic bodies, was published by Liu [14], which relies on his paper [18]. In this, anisotropic materials properties in preferential directions were characterized by a number of unit vectors m 1 , . . . , m a and tensors M 1 , . . . , M b . Ig g is a group of transformations which preserve these characteristics, i.e., Here, v is a vector, A is a second-order tensor and f is either scalar-valued, vector-valued or tensor-valued function. Particularly, iff (v, A, m, M) is an isotropic function then: • for a scalar-valued function • for a vector-valued function • and for a tensor-valued function In the same paper, I-Shih Liu gives a list of 14 such groups g for some crystal classes. He uses the following notation: n i , i 1, 2, 3 where n i are orthonormal vectors, i.e., N i are skew-symmetric tensors defined by In order to write the exact form of the constitutive functions used in our further investigation, we need several functional relations particularly among n i and N i . These basic relations are given in "Appendix." Liu [14] considered only two different classes of transversally isotropic bodies g 2 (O; Qn 1 n 1 ), where n 1 is the preferred direction of transverse isotropy. In applying the isotropic representation of constitutive functions, instead of the Cauchy-Green strain tensor C, he used the Green-St. Venant strain tensor E, i.e., which vanishes when there is no deformation and considers constitutive functions q κ and κ of (E, θ, g κ ) up to bilinear terms in E and g κ , i.e., q κ (a 1 + a 2 trE + a 3 n · En)g κ + a 4 Eg κ where all the material coefficients are functions of the temperature θ only. For the class of transversally isotropic bodies defined by he was able to prove that Λ ε is a function of the temperature only, i.e., Further, he obtained the entropy flux and the heat flux relation as: Therefore, for this class of transversally isotropic bodies, the classical result does not hold in general.
For the class of transversally isotropic bodies defined by which is identical to the classical result (2). For these cases, a functional form of Λ ε had to be found first to determine the relation between entropy flux and heat flux. It appears that Λ ε Λ ε (θ ) holds in all these cases, i.e., Λ ε Λ ε (θ ) is a necessary condition for the determination of a relation between the entropy flux and the heat flux. This is where we pose the question whether this is sufficient condition? Having this in mind, we proceed first by re-examining the above cases.
This assumption differs substantially from Green and Laws [19], as well as from Hutter [20] and Bargmann and Steinmann [21]. In Green and Laws [19], the entropy flux and heat flux relationships are defined as κ 1 ϕ q κ , where ϕ is a constitutive function which reduces to the absolute temperature θ in equilibrium. Hutter [20] postulated the classical entropy flux heat flux relation. In their contribution, Bargmann and Steinmann [21] adopted the Green and Naghdi approach for non-classical theory of thermos-elasticity for isotropic materials and to obtain the entropy flux-heat flux relation (4).The consequence of the assumption that Lagrange multiplier conjugated to the energy balance equation is function of temperature only Λ ε Λ ε (θ ).
Eringen stated [22] that "it is always possible to express the entropy change as a sum of entropy flux and entropy source as" where 1 is the entropy change due to all other effects except heat input. In our consideration, we do not use Eringen's postulate and do not make any assumption about the entropy flux relation. The only assumption used in our derivation is that Lagrange multiplier Λ ε is a function of temperature θ only, i.e., Λ ε Λ ε (θ ). The implication of this assumption for isotropic viscoelastic bodies is considered first. In this, the starting point is the set of Eq. (17) which are restated below Recall that κ f κ C,Ċ, θ, g κ and q κ q κ C,Ċ, θ, g κ , consequently k k C,Ċ, θ, g κ . Using k, Eq. (35) can be written as: Thus, k k (θ, g κ ) and the set of Eq. (36) reduces to ∂k ∂g κ sym 0 (37) For clarity reasons, the following few steps are written in index notation or component notation. Equation (37) when written in component form is ∂k i ∂g j + ∂k j ∂g i 0, where k i and g i are components of vectors k and g κ , respectively. Differentiation of (37) with respect to g κ yields The remaining two relations obtained by cyclic index permutation are From them, we have The solution of this simple set of differential equations is This solution can be rewritten, using symbolic notation, as: where A(θ ) is skew symmetric. Since we are dealing with isotropic viscoelastic bodies, k(θ, g κ ) must be a vector-valued isotropic function. Thus must hold for all Q ∈ O and an arbitrary g κ . Equivalently, (39) can be restated as: Particularly, for the case of Q −I, it follows that a(θ ) 0 and consequently (40) reduces to which can be rewritten as: Since (42) must hold for all Q ∈ O and an arbitrary g κ , we have that since A(θ ) is skew symmetric, only A(θ ) 0 satisfies (43), and This result is identical to Liu's [13] for isotropic viscoelastic bodies, which validates the new procedure and led to its application to the anisotropic materials considered below. It is important to observe that the assumption Λ ε Λ ε (θ ) significantly simplifies and shortens the procedure.

Entropy Flux Relation for Anisotropic Elastic Materials in General
This section considers anisotropic bodies characterized with This must hold for all Q ∈ O, and hence i.e., A(θ, m, M) is skew-symmetric tensor-valued isotropic function of its arguments. For vector-valued and skew-symmetric tensor-valued isotropic functions in R 3 see Smith [23] and Spencer [24]. All groups of crystal classes given by Liu [18] are considered below.

Transversally Isotropic Material Bodies
In this section, we consider transversely isotropic material bodies divided into four cases. The cases (a) characterized by g 2 (O; n 1 ) and the case (b) characterized by g 5 (O; n 1 ⊗ n 1 ) were selected in order to validate the proposed procedure against the results obtained by Liu [14] for the same materials. The results for the cases (c) characterized by g 1 (O; n 1 , N 1 ) and d) characterized g 1 (O; n 1 , N 1 ) are new and, to the best of our knowledge, not available in the published literature.
In this case, k(θ, n 1 ⊗ n 1 , g κ ) A(θ, n 1 ⊗ n 1 )g κ + a(θ, n 1 ⊗ n 1 ) where a(θ, n 1 ⊗ n 1 ) and A(θ, n 1 ⊗ n 1 ), as isotropic functions, must vanish, resulting in The result (59) agrees with I-Shih Liu's result for transversally isotropic bodies with group symmetry group g 5 [14]. Therefore, for the above cases, we demonstrated that Λ ε Λ ε (θ ) is not only a necessary but also a sufficient condition to determine the entropy flux-heat flux relation.
To demonstrate the generality of the proposed procedure, we consider the other crystal classes for which representation of anisotropic invariants is given by Liu [18]. Note, in all these cases, representations of anisotropic invariant function are obtained using the tables for isotropic functions. This more general entropy principle delivers relations which are the same as those obtained by the Clausius-Duhem inequality and hold for all materials in classical thermodynamics. Heuristically, we do not see any physically based objection which would contradict the possibility that Λ ε could be different to Λ ε Λ ε (θ ) for crystal classes from G3 to G14. It is also in accordance with assumption given by Hutter [25], p. 211, who demonstrated that Λ ε (θ ) is independent of the material properties for the heat conducting compressible fluids. Further, this could be related to the statement by Ingo Muller and Tommaso Ruggeri in relation to extended thermodynamics: "Physicists firmly believe that the differential equations of nature should be hyperbolic so as to exclude action at a distance. This incompatibility between the expectation of the physicist and the classical laws of thermodynamics has prompted the formulation of extended thermodynamics". In addition, we do not see any physical reason that Λ ε would have a different form to Λ ε Λ ε (θ ) for these crystal classes. The result for the following cases is new (to the best of our knowledge not available in the published literature).
In this case,
isotropic elastic materials as well as for some crystal classes including transverse isotropy, orthotropy, triclinic systems and rhombic systems. First, we re-examined the entropy flux-heat flux relation for viscoelastic materials, isotropic elastic materials and transversely isotropic elastic bodies and demonstrated that our results agree with the results obtained by Liu [12][13][14][15]27]. All these cases confirm that Λ ε Λ(θ ) is a necessary and sufficient condition for the determination of the entropy flux-heat flux relation. Furthermore, we derived the entropy flux-heat flux relations, based on the assumption that Λ ε Λ(θ ), for all the following crystal classes: transverse isotropy, orthotropy, triclinic systems, monoclinic systems and rhombic systems for which representations of anisotropic functions with respect to their symmetry groups can be expressed in terms of isotropic functions. Our derivation is very general in the sense that the constitutive relations are nonlinear. One of our main results is the proof that all crystal elastic bodies, we considered, are hyperelastic. This represents a generalization of I-Shih Liu's finding for transversely isotropic bodies, the only case he analyzed.
We would like to draw attention to the following three points: (i) The vector function a and skew-symmetric function A are isotropic functions depending only on the set (θ, m, M) which simplifies the procedure. (ii) Generally, the classical entropy flux-heat flux relation does not hold; it is true for all crystal classes investigated here except for g 2 , g 6 g 14 . (iii) The heat flux in the absence of a temperature gradient is not zero for all crystal classes. This confirms Eringen's statement that 1 θ q + 1 where 1 is the entropy change due to all other effects except heat input.
We repeat that heuristically, we do not see any physically based objection which would contradict the possibility that Λ ε could be different to Λ ε Λ ε (θ ) for crystal classes we considered. Of course, the assumption Λ ε Λ(θ ) and all our predictions have to be verified experimentally. This is a task for future research in collaboration with some experimentalists.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
where v and w are arbitrary vectors.