Abstract
The paper deals with the statement and the application of the entropy principle, through the Clausius–Duhem inequality, in continuum physics. The conceptual role is taken from the Coleman–Noll paper of 1963 thus leading to the physical admissibility of constitutive equations. The statement is generalized by letting the rate of entropy production be a constitutive property per se. This generalization proves essential in connection with the modelling of hysteretic phenomena. As to the application, the view of the Coleman–Noll procedure is maintained but a representation formula is shown to generalize the consequences of the entropy principle; as an example the modelling of heat conduction is investigated. Furthermore, while applying the entropy principle to magnetic materials, it is shown an interesting connection between the balance of angular momentum and the thermodynamic restrictions. Also, the modelling through rate-type equations shows the need of Lagrangian fields to obey objectivity and that of the entropy production as a constitutive function to account for the difference between loading and unloading processes.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
Continuum physics approaches are based on two essential tools namely the balance equations and the constitutive functions. Once these tools are settled, practical methods can be applied for solving mathematical problems involving mainly differential equations, boundary-value formulations, and integral equations. To fix ideas here we restrict attention to single-phase continua though the basic concepts apply also to mixtures and continua with internal structure.
The balance equations are established through axioms that are commonly based on global properties. As to the balance of linear momentum, angular momentum, and energy the pertinent equation is assumed as the analogue of the balance for systems of particles (see, e.g., [1], ch. I). Classical thermodynamics provides the analogous reference for the first and second law of continuous bodies; the balance of energy is based on both thermodynamic and mechanical analogues.
In local form, the second law leads to the entropy or Clausius–Duhem (CD) inequality. As we show in the next section, though we are rightly accustomed to the term inequality, the second law too is expressed by an equality.
While the balance of linear momentum and energy involve pertinent external supplies, upon accounting for the balance of energy the second law resulted (and results) in a supply-free inequality, for appropriate constitutive functions. The question about the conceptual role of the CD inequality was answered by Coleman and Noll in 1963 [2] by stating that the CD inequality has to be valid for any set of functions which are compatible with the balance equations. Hence if the balance equations are valid, the CD inequality is a test for the admissibility of the constitutive functions. Though with different methods, the paper by Coleman and Noll initiated much research into the restrictions imposed by thermodynamics on materials models.
The purpose of this paper is threefold. First, to revisit the formulations of the second law (or entropy principle) subsequent to the original one of Coleman and Noll. Secondly, to exhibit the methods of exploitation of the second law, through the CD inequality. Thirdly, to point out recent methods of exploitation and to apply them to the modelling of dissipative processes via rate equations. In particular, attention is addressed to hysteretic models and rate-type equations for heat conductors.
Notation. We let \(\Omega\) be the time-dependent region occupied by the continuous body in the three-dimensional space. The position vector of a point in \(\Omega\) is denoted by \(\textbf{x}\). Hence \(\rho (\textbf{x}, t)\) and \(\textbf{v}(\textbf{x}, t)\) are the mass density and the velocity fields at \(\textbf{x}\), at time \(t \in \mathbb R\). The symbol \(\nabla\) denotes the gradient, with respect to \(\textbf{x}\), while \(\nabla \cdot\) is the divergence operator. For any pair of vectors \(\textbf{u}, \textbf{w}\), or tensors \(\textbf{A}, \textbf{B}\), the notations \(\textbf{u}\cdot \textbf{w}\) and \(\textbf{A}\cdot \textbf{B}\) denote the inner product. Cartesian coordinates are used and then, in the suffix notation, \(\textbf{u}\cdot \textbf{w}= u_i w_i\), \(\textbf{A}\cdot \textbf{B}= A_{ij} B_{ij}\), the summation over repeated indices being understood. Also, \(\textrm{sym}\textbf{A}\) and \(\textrm{skw}\textbf{A}\) denote the symmetric and skew-symmetric parts of \(\textbf{A}\) while \(\textrm{Sym}\) is the space of symmetric tensors. A superposed dot denotes the total time derivative and hence, for any function \(f(\textbf{x}, t)\) on \(\Omega \times \mathbb R\) we have \(\dot{f} = \partial _t f + (\textbf{v}\cdot \nabla )f\). The symbol \(\textbf{L}\) denotes the velocity gradient, \(L_{ij} = \partial _{x_j} v_i\), while \(\textbf{D}= \textrm{sym}\textbf{L}\) and \(\textbf{W}= \textrm{skw}\textbf{L}\). Further, \(\textbf{T}\) is the Cauchy stress tensor, \(\textbf{b}\) is the specific body force, and \(\otimes\) denotes the dyadic product.
Let \(\textrm{R}\) be the region occupied by the body in a reference configuration. Any point in \(\textrm{R}\) is associated with the position vector \(\textbf{X}\) relative to a chosen origin. The motion of the body is a \(C^2\) function \(\varvec{\chi }(\textbf{X}, t): \textrm{R}\times \mathbb R\rightarrow \Omega = \varvec{\chi }(\textrm{R}, t)\). The gradient, with respect to \(\textbf{X}\), of \(\varvec{\chi }\) is the deformation gradient \(\textbf{F}\), \(F_{i{\scriptscriptstyle K}} = \partial _{X_{\scriptscriptstyle K}} \chi _i\).
2 Balance equations
Let \(\textrm{P}\subset \textrm{R}\) be any sub-region of \(\textrm{R}\) and denote by \(\mathcal{P}_t\) the image through the motion in that \(\mathcal{P}_t = \varvec{\chi }(\textrm{P}, t)\); we can say that \(\mathcal{P}_t\) convects with the body.
Let \(\Phi (\textbf{x}, t)\) be a density, per unit volume. If \(\Phi\) is differentiable then
Let \(\rho\) be the mass density and hence
is the mass of the convecting region \(\mathcal{P}_t\). The conservation of mass is expressed by
Applying Eq. (1) and using the continuity of the integrand and the arbitrariness of \(\mathcal{P}_t\) we find that
Consequently, for any specific density \(\varphi\) we have
this result being referred to as Reynolds’ transport theorem.
The global form of a balance law is stated by assuming that the rate of a volume integral consists of volume and surface integrals. If \(\rho \beta\) is the supply, per unit volume, and \(\lambda\) the supply per unit area then the balance law is given the form
If \(\varphi = \textbf{v}\) this is the balance of linear momentum and \(\beta = \textbf{b}\) is the body force while \(\lambda = \textbf{t}\) is the tension. Cauchy’s theorem implies the existence of the stress tensor so that \(\textbf{t}(\textbf{x}, \textbf{n}, t) = \textbf{T}(\textbf{x}, t) \textbf{n}\). Hence we find the local form
By the balance of angular momentum it follows that \(\textbf{T}= \textbf{T}^T\); for materials with couple stress tensors, the Cauchy stress \(\textbf{T}\) is no longer symmetric.
A purely mechanical balance of energy is not consistent except for particular cases (e.g. incompressible fluids). A general form is then considered with an energy density \(\varepsilon\), a supply \(\rho r\) and a surface power density h,
A Cauchy-like proof shows that there is a vector field, say \(- \textbf{q}\), such that
Using Eqs. (3) and (5) we then find the local form
the consistency with the classical first law of thermodynamics allows us to view \(\varepsilon\) as the specific internal energy, r the specific heat supply, and \(\textbf{q}\) the heat flux vector.
Let \(\eta\) be the specific entropy density. The balance of entropy too can be assumed in the form (4) with
Again a Cauchy-like theorem shows that there is a vector field, say \(\textbf{j}\), such that
\(\textbf{j}\) being viewed as the entropy flux. We then obtain the balance of entropy in the local form as
Simple models in classical thermodynamics (see, e.g., [3], ch. 9; [4], ch. 6) indicate that the rate of entropy is related to the rate of heat divided by the absolute temperature \(\theta\). Further, since entropy is not conserved then the balance of entropy is subject to the general requirement (e.g. [4], §6.5)
This suggests that
where \(\gamma\) is the specific (rate of) production of entropy. Hence the local form of the balance of entropy can be written as
The assumption
then results in the relation
let us call it the entropy inequality. The definition (7) denotes that both \(r/\theta\) and \(\gamma\) happen as volume densities; \(r/\theta\) arises from an external source [12], while \(\gamma\) is an internal production.
3 Entropy principle and entropy inequality
The non-negative character of the entropy production has been taken by Coleman and Noll as a criterion for the selection of admissible constitutive equations. To make this point precise let a thermodynamic process be the set of fields entering the model of the continuum (through balance and constitutive equations). Then we state the following
Entropy principle A thermodynamic process is admissible if the inequality
is valid.
It is natural to ask about the presumed argument, by Coleman and Noll, underlying the statement of the entropy principle. Also, we need operative methods for the exploitation of the entropy principle as a procedure for the selection of admissible constitutive equations. Before answering these questions we point out that different expressions have been established of the entropy production while preserving the entropy principle.
Consistent with the view that the rate of entropy is related to the rate of heat divided by the absolute temperature, Coleman and Noll [2] and Truesdell and Toupin ([5], §257), assume \(\textbf{j}= \textbf{q}/\theta\). Hence the entropy inequality (9) becomes
The requirement (10) is given in [5], §258, and is regarded as a postulate of irreversibility. Accordingly, sometimes (e.g. [6]) equation (9) is called the Clausius–Duhem–Truesdell–Toupin inequality; the name Clausius–Duhem inequality would be more appropriate for the special case \(r = 0\). Really the name Clausius–Duhem inequality was given to (10) by Truesdell himself [7], Eq. (28.5), for the case \(r =0\). Following a wide literature, we keep denoting (9), and the particular case (10), as the entropy or Clausius–Duhem inequality.
In 1967 Müller [8] generalized the entropy inequality by letting \(\textbf{j}\) as given by an unknown constitutive function and hence avoiding the assumption \(\textbf{j}= \textbf{q}/\theta\). In this connection it is standard to let
\(\textbf{k}\) being viewed as the extra-entropy flux. Accordingly the expression of \(\gamma\) is given by
Some years later, Green and Naghdi [9] while keeping \(\textbf{j}= \textbf{q}/\theta\) as in the original CD inequality, observed that Eq. (8) is an equation
and furthermore that the entropy production \(\gamma\), as well as \(\psi , \eta , \textbf{T}, \textbf{q}\), is specified by a constitutive equation (while \(\textbf{k}= \textbf{0}\) in Ref. [9]). If \(\textbf{j}= \textbf{q}/\theta + \textbf{k}\) then (12) is replaced by
a priori, the sign of \(\nabla \cdot \textbf{k}\) is undetermined whereas \(\gamma\) is required to be non-negative. If \(r = 0\) then Eq. (12) can be written in the form
still with \(\gamma \ge 0\) and sign of \(\nabla \cdot \textbf{k}\) a priori undetermined. Letting both \(\textbf{j}\) and \(\gamma\) be given by constitutive equations is a generalization of the entropy principle by Coleman and Noll [2] and, as shown by a number of examples, allows for a consistent framework of behaviours in matter.
To give a possible motivation underlying the statement of the entropy principle we observe that the balance of linear momentum and energy express the equilibrium between the pertinent fields and the supplies (\(\textbf{b}, r\)). The balance of mass instead is a constraint on the fields \(\rho\) and \(\textbf{v}\). Rather, once we substitute for r from (6) into (9), the CD inequality (11), possibly with \(\textbf{k}= \textbf{0}\), involves constitutive equations but not supplies. It is then a constraint on the constitutive functions subject to the condition \(\gamma \ge 0\). If \(r = 0\) then (12), with \(\gamma \ge 0\), is a constraint.
3.1 Remarks on the rate of entropy production
It is worth reporting the objection, by an anonymous reviewer, to viewing the (rate of) entropy production as independent of other conventional constitutive equations. The entropy production \(\gamma\) is just an abbreviation for the left-hand-side (LHS) of e.g. Eq. (9). Indeed, in examining some of the standard books in continuum mechanics this is made explicitly clear. For example, in Ch. 5 of [10], \(\gamma\) is not even defined; Liu just uses directly the inequality displayed in Eq. (9). Or one can examine § 4.6 of [11] or §§ 4.3 and 4.11 of [12] where it is explicitly noted that \(\gamma\) is defined by the LHS of Eq. (9).
Answering the objection, we say that for many models the entropy production \(\gamma\) is in fact given by the other constitutive functions. For instance, in Fourier-like models of heat conduction the CD inequality reduces to
and hence \(\gamma = - (\textbf{q}\cdot \nabla \theta )/\rho \theta ^2\) is determined by the function \(\textbf{q}\), provided \(\gamma \ge 0\). Instead there are cases where \(\gamma\) is independent of the other constitutive functions. In Sect. 7.2 we show that this happens for the modelling of magnetic hysteresis where, e.g., \(\gamma = \zeta |\dot{H}|\) is non-negative and independent of the other functions, see Eqs. (35) and (36). Further examples are given, e.g., in [13, 14].
4 Exploitation methods of the entropy inequality
In light of the entropy principle, and the relation (9) for \(\rho \theta \gamma\), the CD inequality can be written in the form
The entropy inequality and the entropy principle hold for any single-phase, non-polar continuum; things change according to the nature of the continuum as to the balance equations and hence for the analogue of (11). To fix ideas we let (13) be the form of the CD inequality; later on we will consider deformable magnetic continua.
For technical reasons it is often convenient to work within the Lagrangian description. In this connection we multiply Eq. (13) by \(J = \det \textbf{F}> 0\) and write
Let
denote the second Piola stress, the referential heat flux and the referential extra-entropy flux. Notice that, for any function \(f(\textbf{x}, t)\),
Further,
where \(\textbf{E}= (\textbf{F}^T \textbf{F}- \textbf{1})/2\). By (14), since \(J \rho = \rho _{\scriptscriptstyle R}\) is the referential density then we can write the referential entropy inequality in the form
where \(\psi _{\scriptscriptstyle R}= \rho _{\scriptscriptstyle R}\psi , \eta _{\scriptscriptstyle R}= \rho _{\scriptscriptstyle R}\eta\).
Granted the balance equations, the entropy principle, through the CD inequality, becomes a criterion of admissibility of the assumed constitutive functions. To derive the restrictions on the constitutive equations two procedures are used.
4.1 The Coleman–Noll procedure
To exploit the requirement of the entropy principle it is essential to envisage the freedom allowed by the balance equations for the admissible thermodynamic processes. In dynamic problems \(\textbf{b}\) and r are usually regarded as assigned a priori. Here we take the view that conceptually the functions \(\textbf{b}(\textbf{x}, t)\) and \(r(\textbf{x}, t)\) can be selected arbitrarily. Hence for any value of the right-hand sides of
the balance equations hold because \(\textbf{b}\) and r are just taken to be given by (16). Analogous claims hold for the space and time derivatives. Hence Eq. (16) provide the body force \(\textbf{b}\) and the heat supply r, and possibly their derivatives, that must be applied to support the process. Based on this view, we find that, given any spatial point \(\textbf{x}\) and any time t, it is possible to find a process such that \(\dot{\theta }, \ddot{\theta }, \nabla \dot{\theta }\) and \(\textbf{L}, \dot{\textbf{L}}, \ddot{\textbf{L}}, \nabla \textbf{L}\) or \(\dot{\textbf{F}}, \ddot{\textbf{F}}\) have arbitrarily prescribed values at that point and time (see, e.g., [16], §3.4).
Let \(\Sigma\) be a set of variables and assume the entropy inequality takes the form
If \(\alpha\) can take arbitrary (positive or negative) values, independent of \(\Sigma\) and \(\xi\), then (17) implies \(f = 0\) and \(p \ge 0\). Likewise, if vectors are involved then
implies \(\textbf{f}= \textbf{0}\). If tensors are concerned in the form
than the arbitrariness of \(\textbf{A}\) implies \(\varvec{\mathcal{F}}= \textbf{0}\). If \(\textbf{A}\in \textrm{Sym}\) (or \(\textrm{Skw}\)) then \(\varvec{\mathcal{F}}\in \textrm{Skw}\) (or \(\textrm{Sym}\)).
If (17) is changed to
then dividing by \(|\alpha |\) we have
Hence it follows that \(f(\Sigma ) = 0\). The analogue holds for \(\textbf{f}\) and \(\varvec{\mathcal{F}}\).
4.2 The Liu procedure
In a different procedure established by Liu [17] the entropy inequality is exploited for supply-free bodies. Hence all of the balance equations are constraints and the CD (or entropy) inequality is written in the (Eulerian) form
\(\Lambda _\rho , \varvec{\Lambda }_\textbf{v}, \Lambda _\varepsilon\) being the Lagrange multipliers, viewed as additional constitutive functions of the chosen independent functions. The linearity and arbitrariness of \(\partial _t \rho , \partial _t \theta , \partial _t \textbf{v}\) imply that
whence we obtain formally the Lagrange multipliers \(\Lambda _\rho , \Lambda _\varepsilon , \varvec{\Lambda }_\textbf{v}\). Further restrictions are derived via the reduced inequality and the appropriate constitutive assumptions.
The conceptual difference of the two methods is given by the use of the supplies \(\textbf{b}, r\). In the Coleman–Noll procedure we make use of the conceivable arbitrariness of the functions \(\textbf{b}(\textbf{x}, t), r(\textbf{x}, t)\) whereas in Liu’s procedure \(\textbf{b}= \textbf{0}, r =0\). Henceforth we apply the Coleman–Noll procedure and show how further methods are applicable to investigate the compatibility with the entropy principle.
5 Exploitation through a representation formula
From the mathematical side, the entropy principle requires that we determine the constitutive functions that satisfy the balance equations and inequality (13). About the inequality, and still into line with Coleman–Noll’s procedure, we need an improvement when the entropy inequality has the form
for the unknown constitutive tensor (or vector) function \(\dot{\textbf{A}}\). We then look for the solution \(\dot{\textbf{A}}\) to
5.1 Representation formula
Let \(\textbf{N}\) be a tensor, of any order n, subject to \(\textbf{N}\cdot \textbf{N}= 1\). Any tensor \(\textbf{Z}\), of the same order, can be decomposed in the form
Assume \(\textbf{Z}\cdot \textbf{N}\) is known, say \(\textbf{Z}\cdot \textbf{N}= g\). If \(\textbf{Z}_{\scriptscriptstyle {\perp }}\) is unknown we can represent it in the form
where \(\textbf{G}\) is an arbitrary tensor of order n and \(\varvec{{\mathsf I}}\) is the identity 2n-order tensor.
Let \(\textbf{Z}\) be a vector subject to
and select
Consequently
and then
If \(\Gamma\) is the set of variables then g and \(\textbf{G}\) are functions of \(\Gamma\).
5.2 Application to the modelling of heat conduction
We consider a thermoelastic body and follow a Lagrangian description. We have in mind that heat conduction is described by a rate equation and hence we let \(\theta , \textbf{E}, \textbf{q}_{\scriptscriptstyle R}, \nabla _{\scriptscriptstyle R}\theta\) be the set of variables for the constitutive functions
Computation of \(\dot{\psi }_{\scriptscriptstyle R}\) and substitution in Eq. (15) result in
By applying repeatidly the argument associated with (17) we conclude that the factors of \(\dot{\theta }, \dot{\textbf{E}}, \nabla _{\scriptscriptstyle R}\dot{\theta }\) should vanish, which implies that
Likewise we find that \(\textbf{k}_{\scriptscriptstyle R}\) is independent of \(\textbf{E}, \textbf{q}_{\scriptscriptstyle R}, \nabla _{\!\!{\scriptscriptstyle R}}\theta\) and hence
This term might be associated with \(\textbf{q}_{\scriptscriptstyle R}\cdot \nabla _{\scriptscriptstyle R}\theta\); yet the isotropy of the body implies the vanishing of a vector \(\textbf{k}_{\scriptscriptstyle R}\) produced by the scalar \(\theta\). Thus we let \(\textbf{k}_{\scriptscriptstyle R}= \textbf{0}\). Hence we are left with
Equation (19), along with the balance equations (2), (5), and (6), gives a system of differential equations for the unknowns \(\theta , \varvec{\chi }, \textbf{q}_{\scriptscriptstyle R}\) and \(\rho = \rho _{\scriptscriptstyle R}/J, J = [\det (2\textbf{E}+ \textbf{1})]^{1/2}\), with \(\rho _{\scriptscriptstyle R}\) a given constant. The form of (19) depends on the functions \(\psi _{\scriptscriptstyle R}(\theta , \textbf{E}, \textbf{q}_{\scriptscriptstyle R})\) and \(\gamma (\theta , \textbf{E}, \textbf{q}_{\scriptscriptstyle R}, \nabla _{\scriptscriptstyle R}\theta )\). Using the representation formula (18) we now derive explicit forms of \(\dot{\textbf{q}}_{\scriptscriptstyle R}\).
Assume \(\partial _{\textbf{q}_{\scriptscriptstyle R}}\psi _{\scriptscriptstyle R}\ne \textbf{0}\) and let \(\textbf{N}= \partial _{\textbf{q}_{\scriptscriptstyle R}}\psi _{\scriptscriptstyle R}/|\partial _{\textbf{q}_{\scriptscriptstyle R}}\psi _{\scriptscriptstyle R}|\). Using the representation formula (18), applied to vectors, and Eq. (19) we find
\(\textbf{g}\) being any vector function of \(\theta , \textbf{E}, \textbf{q}_{\scriptscriptstyle R}, \nabla _{\!\!{\scriptscriptstyle R}}\theta\). For definiteness, if \(\textbf{g}= \textbf{J}_{\scriptscriptstyle R}\nabla _{\!\!{\scriptscriptstyle R}}\theta\), where \(\textbf{J}_{\scriptscriptstyle R}\) is a second-order tensor possibly dependent on \(\theta , \textbf{E}, \textbf{q}_{\scriptscriptstyle R}\), then we have
where
We now go back to Eq. (20) and consider the particular case \(\textbf{g}= \textbf{0}\). Hence we have
Since \(\gamma > 0\) it follows that, in stationary conditions,
this is a necessary condition, subject only to \(\partial _{\textbf{q}_{\scriptscriptstyle R}} \psi _{\scriptscriptstyle R}\ne \textbf{0}\), and provides the classical result of heat conduction inequality.
An interesting particular case of (20) follows by letting \(\partial _{\textbf{q}_{\scriptscriptstyle R}} \psi _{\scriptscriptstyle R}= \alpha \textbf{q}_{\scriptscriptstyle R}\) so that
Hence selecting \(\textbf{g}= - \nabla _{\!\!{\scriptscriptstyle R}}\theta /\theta \alpha\) we have
Equation (21) has the standard form of the Maxwell-Cattaneo equation (see, e.g., [18, 19]) with relaxation time \(\tau\) and conductivity \(\kappa\) given by
If, rather, the conductivity \(\kappa\), e.g. as a function of \(\theta\), is given then we determine the entropy production,
as a function of the variables \(\textbf{q}_{\scriptscriptstyle R}\) and \(\theta\). The positive valuedness of the entropy production \(\gamma\) implies the positive valuedness of the heat conductivity \(\kappa\).
6 Entropy inequality and selection of variables
It is customary to assume that a thermodynamic process satisfies the balance equations and then constitutive restrictions are determined so that the entropy inequality holds. There are cases though where the entropy inequality can be employed to select appropriate independent variables. This feature is exemplified through the modelling of electromagnetic continua.
For simplicity consider a deformable magnetic body. Let \(\textbf{m}\) and \(\textbf{M}\) be the magnetization per unit mass and unit volume. The balance of angular momentum leads to ([16], §2.16.1)
Owing to the balance of energy,
we obtain the entropy inequality in the form
Since
we let
and write the entropy inequality in the form
If \(\textbf{F}, \textbf{H}\) and e.g. \(\theta , \nabla \theta\) are the variables then we have
the dots denoting the terms in \(\dot{\theta }, \dot{\textbf{H}}, (\nabla \theta )\dot{}\). Upon substitution of \(\dot{\phi }\) we can write Eq. (24) in the form
The linearity an arbitrariness of \((\nabla \theta )\dot{}, \dot{\theta }, \textbf{L}, \dot{\textbf{H}}\) imply that
We now investigate the compatibility of the symmetry requirement (22) with the thermodynamic restrictions (25). The function \(\phi\) might depend separately on \(\textbf{F}\) (e.g. through \(|\textbf{F}|\) or \(|\textbf{E}|\)) and \(\textbf{H}\) (through \(|\textbf{H}|\)). Hence we notice that
Consequently it follows
or
and
Hence
Thus the symmetry requirement (22) holds in the special case where \(\textbf{T}\in \textrm{Sym}\) and \(\textbf{H}\otimes \textbf{M}\in \textrm{Sym}\).
More realistic models should allow for a cross dependence on \(\textbf{F}\) and \(\textbf{H}\) so that a magnetomechanical coupling [29] is obtained. Both \(\textbf{F}^T \textbf{H}\) (i.e. \((\textbf{F}^T \textbf{H})_{\scriptscriptstyle P}= F_{i {\scriptscriptstyle P}} H_i\)) and \(\textbf{F}^{-1} \textbf{H}\) are invariant under Euclidean transformations (SO(3) invariant). Based on SO(3) invariance we might look for a dependence on \(\textbf{F}^T \textbf{H}\) or \(\textbf{F}^{-1} \textbf{H}\) possibly multiplied by any power of the scalar invariant \(J = \det \textbf{F}\). First we check the SO(3) invariance of \(J^n \textbf{F}^T \textbf{H}\) and \(J^n \textbf{F}^{-1} \textbf{H}\), with \(n \in \mathbb N\).
Under the Euclidean transformation
and the like for \(J^n \textbf{F}^{-1} \textbf{H}\).
Now let
and assume \(\phi = \phi (\theta , \textbf{E}, \textbf{A})\). We find that
Though \(\textbf{T}\) and \(\textbf{H}\otimes \textbf{M}\) are non-symmetric we have
and then \(\textbf{T}+ \mu _0 \textbf{H}\otimes \textbf{M}\in \textrm{Sym}\). Hence the dependence on \(J^n \textbf{F}^T \textbf{H}\) makes the symmetry requirement (22) valid for any \(n \in \mathbb N\). It is worth remarking that this conclusion is consistent with the literature where the field \(\textbf{A}= \textbf{F}^T \textbf{H}\) (\(n = 0\)) is adopted as the Lagrangian form of the magnetic field [20,21,22,23].
Now, \(J^n \textbf{F}^{-1} \textbf{H}\) is also objective and hence we may look for the possible function \(\phi (\theta , \textbf{E}, \varvec{\mathcal{A}})\), with
Since
then we find
Consequently
Hence \(\textbf{T}+ \mu _0 \textbf{H}\otimes \textbf{M}\in \textrm{Sym}\) and also the dependence on \(\varvec{\mathcal{A}}\) is objective and satisfies the requirement (22). Nevertheless the dependence on \(\varvec{\mathcal{A}}\) seems to be uncommon in the literature.
7 Dissipation through the entropy production or the dissipation potential
Dissipative phenomena are associated with a positive entropy production. Indeed, relative to the CD inequality, dissipative phenomena are characterized by the reduced inequality. Here we show how different approaches to the analysis of the reduced inequality lead to qualitatively different models of continua.
For definiteness we consider magnetic solids subject to the fields \(\textbf{H}, \textbf{M}, \textbf{B}\). For simplicity it is assumed that the material is thermally and electrically non-conducting. The possible electric field \(\textbf{E}\) is then viewed as a vector parameter.
7.1 Relation to some reversible-irreversible decompositions
We review Ref. [24] and start with the balance of energy in the form
where u is the internal energy density. The mechanical power \(\textbf{T}\cdot \dot{\varvec{\varepsilon }}\) means that we assume \(\textbf{T}\in \textrm{Sym}\) and \(\textbf{D}\) is replaced with the time derivative of the infinitesimal strain \(\varvec{\varepsilon }\). The magnetic power \(\textbf{H}\cdot \dot{\textbf{B}}\) is consistent with the expression \(\mu _0 \textbf{H}\cdot \dot{\textbf{M}}\) in that \(\rho\) is assumed to be constant and hence (26) is consistent with (23) by letting \(\rho \varepsilon = \rho u - \mu _0 |\textbf{H}|^2/2\). In light of (26) the CD inequality reads
We assume the additive decomposition of \(\textbf{B}\) and \(\varvec{\varepsilon }\) in the form
the irreversible parts are considered to describe dissipative effects, here magnetostriction or hysteresis. A possible internal variable \(\varvec{\xi }\) is here ignored. Let S be the entropy density, \(\Psi = \rho u - \theta S\) the free energy density. We let \(\Psi , S, \textbf{T}, \textbf{H}\) be functions of the variables
Inequality (27) becomes
The arbitrariness of \(\dot{\theta }, \dot{\textbf{B}}\), and \(\dot{\varvec{\varepsilon }}^\textrm{rev}\) implies
and hence
For technical convenience the free energy \(\Psi\) is then assumed to be a function with appropriate dependences, \(\Psi = \Psi (\theta , \textbf{B}- \textbf{B}^\textrm{irr}, \varvec{\varepsilon }^\textrm{rev})\). Hence it follows
and the reduced inequality (28) can be written in the form
To satisfy (29), appeal is made to a function \(\varphi (\textbf{H}, \textbf{T})\) and the rates \(\dot{\textbf{B}}^\textrm{irr}, \dot{\varvec{\varepsilon }}^\textrm{irr}\) are required to be given by
while
Granted these conditions, \(\varphi\) is called a pseudopotential [25]. If, further, \(\varphi = \hat{\varphi }(|\textbf{H}|, |\textbf{T}|)\) then the inequality reads
and holds if \(\hat{\varphi }\) is convex, relative to \(|\textbf{H}|\) and \(|\textbf{T}|\), with \(\hat{\varphi }(0,0) = 0\). More generally we might have
where \(\textbf{A}\) (\(\textbf{C}\)) is a second-order (fourth-order) positive definite tensor. In that case we have
if only the derivatives \(f_1, f_2\), with respect to the first and second variable, are non-negative.
Thus the model is characterized by two potentials, \(\Psi\) and \(\varphi\). Yet, as a general remark on the use of pseudopotentials, we observe that the functions \(\Psi\) and \(\varphi\) depend on different variables,
respectively.
A similar approach is developed in [26] where the variables are
where \(\textbf{C}\) is the Cauchy-Green tensor, \(\textbf{B}\) is the Lagrangian induction field; the additive decomposition is assumed for \(\textbf{B}\),
\(\textbf{B}_e, \textbf{B}_v\) being viewed as elastic and viscous parts. Instead, the deformation gradient \(\textbf{F}\) is factorized as
and \(\textbf{C}_e = \textbf{F}_e^T \textbf{F}_e\), \(\textbf{C}_v = \textbf{F}_v^T \textbf{F}_v\). A free energy density \(\Omega\) is considered and the CD inequality is written in the form
where \(\textbf{T}_{{\scriptscriptstyle R}{\scriptscriptstyle R}}\) is the second Piola stress to within a pressure term due to the possible incompressibility constraint. It follows
Next, for simplicity of the model it is assumed that the separate inequalities
are satisfied while
The second requirement of (30) is satisfied by letting
It is worth mentioning that this type of solution is adopted also by appealing to relaxation toward equilibrium [27]. Indeed, if e.g. \(\Omega _v\) depends on \(\textbf{B}_v\) in the form
then
if \(\textbf{B}\) is constant then \(\textbf{B}_v - \textbf{B}\) approaches zero as t approaches infinity.
Differently from the approach of [24], in [26] the set of variables is kept unchanged when establishing the equation for the rate \(\dot{\textbf{B}}_v\).
In both approaches the set of variables does not involve the time derivative \(\dot{\textbf{B}}\), or \(\dot{\textbf{B}}^\textrm{irr}\), and hence hysteretic models would seem out of reach. Now, in [28] the dissipation potential is assumed to be a function of \(\dot{\textbf{B}}^\textrm{irr}\) whereas the constitutive functions are not dependent on \(\dot{\textbf{B}}^\textrm{irr}\). Next, by Legendre transformations, a constitutive equation is established in the form
and experimental data are simulated through
where \(\mu\) is the permeability and \(\nu , Z\) are parameters. This system follows by suitable assumptions on the dissipation potential, the dependence on an additional internal variable, \(\varvec{\xi }\), and Legendre transformations. Things would be different by maintaining the same variables for all of the constitutive functions.
7.2 Hysteretic models through the entropy production
Consider the CD inequality (24) and, for formal simplicity, neglect heat conductivity so that
Hysteretic models necessarily involve the rates of the magnetic fields \(\textbf{H}, \textbf{B}, \textbf{M}\). Since the time rates \(\dot{\textbf{H}}, \dot{\textbf{B}}, \dot{\textbf{M}}\) are not objective then we consider an objective analogue of \(\textbf{H}\), namely
which is Euclidian invariant. Now,
and
Substitution in (32) results in
where
Since we have in mind hysteretic effects in deformable ferromagnets, the minimal set of variables is \(\theta , \textbf{F}, \pmb {{\mathfrak {H}}}, \pmb {{\mathfrak {M}}}, \dot{\pmb {{\mathfrak {H}}}}\) or possibly with \(\dot{\pmb {{\mathfrak {H}}}}\) replaced by \(\dot{\pmb {{\mathfrak {M}}}}\). The free energy \(\phi\) is Euclidean invariant and hence we let
Computation of \(\dot{\phi }\) and substitution in (33) yields
The linearity and arbitrariness of \(\dot{\theta }, \ddot{\pmb {{\mathfrak {H}}}}, \textbf{W}\) imply
If \(\dot{\textbf{E}}\) is independent of \(\dot{\pmb {{\mathfrak {H}}}}\) and \(\dot{\pmb {{\mathfrak {M}}}}\) then it follows
Consequently we obtain the reduced inequality
where \(\gamma (\theta , \textbf{F}, \pmb {{\mathfrak {H}}}, \pmb {{\mathfrak {M}}}, \dot{\pmb {{\mathfrak {H}}}})\).
Inequality (34) is the condition for the occurrence of dissipative processes. To describe hysteretic effects we need a distinction between increasing and decreasing of the source field. We may consider the magnetic field as the source field and then we model the entropy production in the form
Since \(\gamma \ge 0\) then (34) is satisfied and we can view
as the rate equation governing the hysteretic processes.
For simplicity, and also with reference to experimental setups [29, 30], we restrict attention to one-dimensional models. Let \(\textbf{e}\) be a fixed direction and let
We observe that \(\dot{H}\) and \(\dot{M}\) are Euclidean invariants. This is so because
and the like for M. Hence H, M are scalar Euclidean invariants and so are \(\dot{H}, \dot{M}\).
In a linear approximation, relative to \(H, M, \dot{F}\) we can write
Hence Eq. (36) is approximated by
Subject to \(\dot{H} \ne 0\), it is
that is the differential magnetic susceptibility. Thus dividing (37) by \(\dot{H}\) we have
The occurrence of \(\textrm{sgn}\,\dot{H}\) produces different dependences of M in terms of H, in hysteretic systems, depending on whether H increases or decreases. Furthermore Eq. (38) is determined by the free energy \(\phi\) and the entropy production \(\gamma\) through \(\zeta\). It is worth remarking that the crucial term \(\textrm{sgn}\,\dot{H}\) arises because of the constitutive property of the entropy production. The free energy \(\phi\) cannot depend on \(\dot{\pmb {{\mathfrak {H}}}}\) because \(\phi\) occurs through the time derivative \(\dot{\phi }\); instead \(\gamma\) occurs through the value possibly dependent on \(\dot{\pmb {{\mathfrak {H}}}}\). Examples of hysteresis cycles determined by equations like (38) are given in [16], ch. 15.
8 Conclusions
The paper deals with the application of the entropy principle, through the Clausius–Duhem inequality, for the analysis of constitutive properties in continuum physics. Though the application is meant in the sense of Coleman–Noll procedure [2], some generalizations are established. First, the use of the representation formula (18) allows us to solve equations of the form (19) which are customary when dealing with vectors and tensors.
Secondly, the Clausius–Duhem inequality may lead to the suitable selection of variables as is the case of the restrictions for magnetic solids,
subject to the balance equation
The Lagrangian variable \(\pmb {{\mathfrak {H}}}= \textbf{F}^T \textbf{H}\) proves to satisfy (39) and in addition is invariant under Euclidean transformations.
Thirdly, letting the entropy production \(\gamma\) be a constitutive function per se allows the modelling of hysteresis. As a remarkable example, in ferromagnetism we let the constitutive functions depend also on \(\dot{\pmb {{\mathfrak {H}}}} \simeq \dot{\textbf{H}}\) and the dependence of \(\gamma\) on \(\dot{\pmb {{\mathfrak {H}}}}\) allows us to select a non-negative function \(\zeta |\dot{\textbf{H}}|\) which, upon division by \(\dot{H}\), eventually distinguishes the susceptibility dM/dH in magnetizing and demagnetizing processes. In relation to other approaches in the literature, the analogue dependence on \(\dot{\textbf{B}}^\textrm{irr}\) is obtained by using a Legendre transformation thus establishing the model with two different sets of variables.
References
Truesdell C (1977) A first course in rational continuum mechanics, vol I. Academic Press, New York
Coleman BD, Noll W (1963) The thermodynamics of elastic materials with heat conduction and viscosity. Arch Rational Mech Anal 13:167–178
Zemansky MW (1968) Heat and thermodynamics. McGraw-Hill, New York
Moran MJ, Shapiro HN (2006) Fundamentals of engineering thermodynamics. Wiley, Chichester
Truesdell C, Toupin RA (1960) The classical field theories. In: Flügge S (ed) Encyclopedia of physics, vol III/1, pringer-Verlag, Berlin
Man C-S (1989) Thermodynamics based on a work inequality. Arch Rational Mech Anal 106:1–61
Truesdell C (1952) The mechanical foundations of elasticity and fluid dynamics. J Rational Mech Anal 1:125–300
Müller I (1967) On the entropy inequality. Arch Rational Mech Anal 26:118–141
Green AE, Naghdi PM (1977) On thermodynamics and the nature of the second law. Proc R Soc Lond A 357:253–270
Liu I-S (2002) Continuum mechanics. Springer, Berlin
Eringen AC (1980) Mechanics of continua. Krieger, Florida
Paolucci S (2016) Continuum mechanics and thermodynamics of matter. Cambridge University Press, Cambridge
Giorgi C, Morro A (2021) A thermodynamic approach to rate-type models of elastic-plastic materials. J Elast 147:113–148
Giorgi C, Morro A (2022) Materials with memory: viscoelasticity and hysteresis. In: Rega G (ed) 50+ years of AIMETA, Springer, pp 243–260
Gurtin ME, Fried E, Anand L (2011) The mechanics and thermodynamics of continua. Cambridge University Press, Cambridge
Morro A, Giorgi C (2023) Mathematical modelling of continuum physics. Birchkäuser, Cham
Liu I-S (1972) Method of Lagrange multipliers for exploitation of the entropy principle. Arch Rational Mech Anal 46:131–148
Straughan B (2011) Heat waves. Springer, Berlin
Sellitto A, Cimmelli VA, Jou D (2016) Mesoscopic theories of heat transport in nanosystems. Springer, Berlin
Ottenio M, Destrade M, Ogden RW (2008) Incremental magnetoelastic deformations, with application to surface instability. J Elast 90:19–42
Saxena P, Hossain M, Steinmann P (2014) Nonlinear magneto-viscoelasticity of transversally isotropic magneto-active polymers. Proc R Soc A 470:0082
Garcia-Gonzalez D (2019) Magneto-visco-hyperelasticity for hard-magnetic soft materials: theory and numerical applications. Smart Mater Struct 28:085020
Haldar K (2021) Constitutive modeling of magneto-viscoelastic polymers, demagnetization correction, and field-induced Poynting effect. Int J Eng Sci 165:103488
Ho K (2016) A constitutive model for magnetostriction based on thermodynamic framework. J Magn Magn Mater 412:250–254
Germain P, Nguyen QS, Suquet P (1983) Continuum thermodynamics. Trans ASME 50:1010–1020
Saxena P, Hossain M, Steinmann P (2013) A theory of finite deformation magneto-viscoelasticity. Int J Solids Struct 50:3886–3897
Gurtin ME (1980) Generalized Ginzburg–Landau and Cahn–Hilliard equations based on a microforce balance. Physica D 92:178–192
Ho K (2016) A thermodynamically consistent model for magnetic hysteresis. J Magn Magn Mater 357412:250–254
Evans PG, Dapino MJ (2013) Measurement and modeling of magnetic hysteresis under field and stress application in iron–gallium alloys. J Magn Magn Mater 330:37–48
Smith RC, Dapino MJ, Seelecke S (2003) Free energy model for hysteresis in magnetostrictive transducers. J Appl Phys 93:458–466
Acknowledgements
The research leading to this work has been developed under the auspices of INDAM-GNFM.
Funding
Open access funding provided by Università degli Studi di Genova within the CRUI-CARE Agreement.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The author declares that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Morro, A. On the entropy inequality and its exploitation in continuum physics. Meccanica 59, 1731–1743 (2024). https://doi.org/10.1007/s11012-024-01804-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11012-024-01804-3