Next Article in Journal
Packed Bed Thermal Energy Storage System: Parametric Study
Previous Article in Journal
An Optimized Artificial Neural Network Model of a Limaçon-to-Circular Gas Expander with an Inlet Valve
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Second Law of Thermodynamics in Continuum Physics

1
Dipartimento di Ingegneria Civile Ambiente Territorio Architettura e Matematica, Università di Brescia, Via D. Valotti 9, 25133 Brescia, Italy
2
Dipartimento di Informatica, Bioingegneria, Robotica e Ingegneria dei Sistemi, Università di Genova, Via All’Opera Pia 13, 16145 Genova, Italy
*
Author to whom correspondence should be addressed.
Thermo 2024, 4(2), 273-294; https://doi.org/10.3390/thermo4020015
Submission received: 16 May 2024 / Revised: 5 June 2024 / Accepted: 7 June 2024 / Published: 11 June 2024

Abstract

:
The paper revisits the formulation of the second law in continuum physics and investigates new methods of exploitation. Both the entropy flux and the entropy production are taken to be expressed by constitutive equations. In three-dimensional settings, vectors and tensors are in order and they occur through inner products in the inequality representing the second law; a representation formula, which is quite uncommon in the literature, produces the general solution whenever the sought equations are considered in rate-type forms. Next, the occurrence of the entropy production as a constitutive function is shown to produce a wider set of physically admissible models. Furthermore the constitutive property of the entropy production results in an additional, essential term in the evolution equation of rate-type materials, as is the case for Duhem-like hysteretic models. This feature of thermodynamically consistent hysteretic materials is exemplified for elastic–plastic materials. The representation formula is shown to allow more general non-local properties while the constitutive entropy production proves essential for the modeling of hysteresis.

1. Introduction

Understanding of the thermodynamics of continuous media has made decisive progress in the twentieth century where the general scheme has been established in terms of balance laws and constitutive relations. The list of balance laws identifies the theory of physics under consideration, e.g., mechanics, electrodynamics, theory of mixtures. The constitutive relations characterize the nature of the continuum, e.g., solid, fluid, gas, hysteretic material. The view that the balance of entropy eventually results in requirements on the physically admissible constitutive relations is due to a well-known paper by Coleman and Noll [1]. The associated postulate is the content of the corresponding second law of thermodynamics and it initiated far-reaching research on the exploitation of the second law for the constitutive relations. It is the purpose of this paper to show some new approaches to the exploitation of the second law. For this, we revisit the various formulations of the second law in Section 3.
It is a common feature of the various statements of the second law that the admissible constitutive relations are subject to the requirement that the entropy production be non-negative. The exploitation of this requirement depends on the form of the constitutive relations (functions, functionals, rate equations). Furthermore, we need to know the proper mathematical expression of the second law and, in particular, to know the expression of the entropy production. Indeed, we regard the entropy production as a constitutive property per se, in addition to being related to other constitutive properties.
The purpose of this paper is to emphasize new aspects associated with the formulation and the exploitation of the second law in continuum physics. Following Müller [2,3], we let the entropy flux, say j , be a constitutive function and not merely the heat flux q divided by the absolute temperature θ . Furthermore, we let the entropy production be a constitutive function.
Three main points have to emerge from this paper. First, the occurrence of a nonzero difference j q / θ proves essential whenever we look for non-local terms involving higher-order gradients of temperature and deformation. Second, in three-dimensional settings, vectors and tensors are in order and they occur through inner products in the inequality representing the second law. A representation formula, quite uncommon in the literature, produces the general solution whenever the sought equations are expressed in rate-type forms. Third, the occurrence of the entropy production as a constitutive function is essential in the thermodynamically consistent modeling of hysteretic materials.
The entropy production allows the completion of rate-type hysteretic equations, as with Duhem-like models. This feature is exemplified in this paper for elastic–plastic materials, though the analogue can be performed for magnetic or electric hysteresis [4]. As is shown in this paper, both the use of the representation formula and the entropy production as a constitutive function turn out to be decisive improvements in the elaboration of material modeling. The representation formula allows for more general non-local properties while the constitutive entropy production results in a direct method for the description of hysteretic materials.

2. Notation and Balance Equations

A body occupies a time-dependent region Ω in the three-dimensional space. The position vector of a point in Ω is denoted by x . Hence, ρ ( x , t ) and v ( x , t ) are the mass density and the velocity fields at x at time t R . The symbol ∇ denotes the gradient with respect to x , while · is the divergence operator. For any pair of vectors u , w , or tensors A , B , the notation u · w and A · B denotes the inner product. Cartesian coordinates are used, and then, in the suffix notation, u · w = u i w i , A · B = A i j B i j , the summation over repeated indices is understood. Also, sym A and skw A denote the symmetric and skew-symmetric parts of A , while Sym is the space of symmetric tensors. A superposed dot denotes the total time derivative, and hence, for any function f ( x , t ) on Ω × R we have f ˙ = t f + ( v · ) f . The symbol L denotes the velocity gradient, L i j = x j v i , while D = sym L and W = skw L . Further, T is the Cauchy stress tensor, b is the specific body force, and ⊗ denotes the dyadic product.
Let R be the region occupied by the body in a reference configuration. Any point in R is associated with the position vector X relative to a chosen origin. The motion of the body is a C 2 function χ ( X , t ) : R × R Ω = χ ( R , t ) . The gradient, with respect to X , of χ is the deformation gradient F , F i K = X K χ i .
The balance of mass is expressed by the continuity equation:
ρ ˙ + ρ · v = 0 .
The equation of motion is written in the form
ρ v ˙ = · T + ρ b .
We assume that there is no internal structure, and then, let T Sym .
Let ε be the specific internal energy density. The balance of energy leads to
ρ ε ˙ = T · D + ρ r · q ,
where r is the heat supply per unit mass, and q is the flux vector.

3. Statements of the Second Law

Let θ be the absolute temperature and η the specific entropy density. We denote as a thermodynamic process the set of fields describing the evolution of the body, namely, ρ , v , T , b , ε , r , q , η , θ . We now revisit the statement of the second law in continuum physics and point out the various formulations that have appeared in the literature.
Let P t Ω be any sub-region that is convected by the motion. As with any balance equation we may express the balance of entropy by letting the rate consist of a volume integral and a surface integral,
d d t P t ρ η d v = P t ρ s d v P t j · n d a ,
and correspondingly viewing s as the entropy supply and j as the entropy flux. The arbitrariness of the region P t , the transport theorem, and the smoothness of the functions ρ , η , s , j imply that
ρ η ˙ = ρ s · j .
Borrowing from classical thermodynamics (e.g., [5]), Coleman and Noll [1] considered r / θ as the external volume supply of entropy, and likewise, assumed that j = q / θ . Hence, they considered the difference
γ = η ˙ r θ + 1 ρ · q θ
as the internal specific (rate of) production of entropy. Accordingly, they stated the following postulate:
For every process admissible in a body the inequality
γ 0
is valid.
This postulate, based on Definition (4), amounts to assuming that Inequality (5), and hence,
ρ η ˙ + · q θ ρ r θ 0 ,
selects the admissible processes. Inequality (5), or (6), is called the Clausius–Duhem (CD) inequality or entropy inequality, while the postulate is viewed as the second law of thermodynamics or entropy principle.
In light of (6), it follows that
ρ η ˙ + 1 θ ( · q ρ r ) 1 θ 2 q · θ 0 .
Replacing · q ρ r from (1) and multiplying by θ we have
ρ θ η ˙ ρ ε ˙ + T · D 1 θ q · θ 0 .
In terms of the Helmholtz free energy ψ = ε θ η , we find
ρ ( ψ ˙ + η θ ˙ ) + T · D 1 θ q · θ 0 .
In 1967, Müller [2] postulated the entropy balance in the form
ρ η ˙ + · j ρ r θ 0 ,
where the entropy flux j need not be equal to q / θ , and furthermore, j has to be determined as any constitutive function.
Next, Green and Laws [6] assumed a modified form of the entropy inequality by replacing the absolute temperature θ in (6) with a non-equilibrium temperature ϕ , which requires a constitutive function and, in equilibrium, reduces to θ .
In 1977, Green and Naghdi [7] wrote the balance of entropy in the form of an equality,
ρ η ˙ = ρ ( r θ + ξ ) · q ,
which, in the previous scheme, amounts to viewing ξ as the entropy production. Yet, they introduced two novelties. Firstly, the entropy production ξ is given by a constitutive relation. Secondly, ξ need not be a non-negative while, as for the postulate about the second law, they assumed that
ρ ( ψ ˙ + η θ ˙ ) + T · D 0
for all thermo-mechanical processes. Note that Equation (10) is recovered from (7) when q = 0 .
Some further comments and statements of the second law appeared later on. In 1990, Maugin [8] (see also [9]) wrote the second law in the form
S ˙ + · S 0 ,
where S = ρ η and S is the entropy flux taken in Müller’s form S = q / θ + k . Yet, the energy supply r is missing and the subsequent procedure leads to the requirement ( S · ) θ 0 , which is quite unusual.
Lately, “non-conventional” statements have been given and corresponding approaches have been developed in Refs. [10,11] by distinguishing equilibrium and non-equilibrium quantities. The stress power w and the heat flux q are considered in the forms w = w ¯ + w and q = q ¯ + q , where w ¯ , q ¯ are the values associated with the local equilibrium state. Hence, the entropy inequality is stated in the form
η ˙ + 1 ρ · q ¯ θ 0 ,
while the balance of energy is written in the form ρ ε ˙ = w ¯ · q ¯ . Again, the energy supply r is missing.
A further approach is due to Dunn and Serrin [12], who posited the existence of a rate of supply of mechanical energy, u, through the boundary of each sub-region, and hence, via a corresponding divergence term · u . So, they assumed the balance of energy and entropy in the form
ρ ε ˙ = T · L · q + · u + ρ r , ρ ( ε ˙ θ η ˙ ) T · L · u + 1 θ q · θ 0 ,
as though j = q / θ . If T Sym , then T · L = T · D .

Second Law and Thermodynamic Processes

Back to the general balance of entropy (3), we let
s = r θ + γ ,
where r / θ and γ denote the external and internal volume supply of entropy. As any flux, j may be viewed as an external entropy contribution to the pertinent sub-region P t . Accordingly, we view γ as a term of internal character, and then, we refer to γ as the (rate of) specific entropy production. Therefore, consistent with Postulate (5), we assume that
ρ η ˙ + · j ρ r θ = ρ θ γ 0
and regard both j and γ as expressed by constitutive relations. Hence, a process is the set P = ( ρ , v , T , ε , η , θ , q , j , γ ) expressed by constitutive relations, while b and r are arbitrary given time-dependent fields on Ω × R . If further fields are involved, such as, e.g., electromagnetic fields, the set P is completed accordingly. The Coleman–Noll postulate is then generalized as follows.
Second law of thermodynamics. For every process P admissible in a body, the inequality (11) is valid at any internal point.
As to boundary points and the required boundary condition, we recall the following.
Principle of the increase in entropy. The entropy of an isolated system cannot decrease in time.
Now, let
j = q θ + k ;
the vector field k is referred to as the extra-entropy flux [3]. Hence, the balance of entropy reads
ρ η ˙ ρ r θ + · q θ = ρ γ · k .
By the principle of the increase in entropy, when r = 0 on Ω and q · n = 0 on Ω , we have
d d t Ω ρ η d v = Ω ( ρ γ · k ) d v 0
or
Ω k · n d a Ω ρ γ d v .
The flow through the boundary Ω of the extra-entropy flux k is bounded by the entropy production in the body.
We append two comments on the properties of k . Firstly, keeping the inequality (6) as valid also when k 0 or letting (13) hold if Ω is replaced with any sub-region P t leads to
ρ γ · k 0 .
Next, we show the consequences of (15) and compare them with those of (11). Secondly, sometimes the boundary condition is taken in the form
Ω k · n d a = 0 .
This condition, which is consistent with (14) and Postulate (5), may be suggested by the mathematical modeling [6,13].
Since j = q / θ + k then Equation (11) can be written in the form
ρ η ˙ + 1 θ ( · q ρ r ) 1 θ 2 q · θ + · k = ρ γ 0 .
Upon replacing · q ρ r from (1), using the Helmholtz free energy,
ψ = ε θ η ,
and multiplying by θ we obtain
ρ ( ψ ˙ + η θ ˙ ) + T · D 1 θ q · θ + θ · k = ρ θ γ 0 .
As we show in the next section, the role of the extra-entropy flux k is crucial in the modeling of materials with higher-order gradients [14].
For later use we now derive the Lagrangian version of (16). Let J = det F > 0 and notice that
J ρ = ρ R ( X )
is the mass density in the reference configuration R . Next, let
T R R = J F 1 T F T , q R = J F 1 q , k R = J F 1 k ,
the referential stress T R R and vectors q R , k R ; T R R is referred to as the second Piola (or Piola–Kirchhoff) stress. The Green–Lagrange strain tensor
E = 1 2 ( F T F 1 )
is related to the stretching D by
E ˙ = F T D F .
Hence, it follows that
J T · D = T R R · E ˙ .
Furthermore, we have
R · k R = J · k , q R · R θ = J q · θ .
Hence, J times Equation (16) yields
ρ R ( ψ ˙ + η θ ˙ ) + T R R · E ˙ 1 θ q R · R θ + θ R · k R = ρ R θ γ 0 .

4. The Extra-Entropy Flux and Materials with Higher-Order Gradients

Non-locality properties in the modeling of materials are often described by a dependence on higher-order gradients. The corresponding thermodynamic consistency is crucially related to the occurrence of a nonzero extra-entropy flux k and to the way the flux k is applied.
For definiteness, here we examine materials where the non-locality is modeled by second-order gradients of temperature and mass density, and then, we let
Γ = ( θ , ρ , θ , ρ , θ ˙ , ρ ˙ , D , θ , ρ )
be the set of variables. The stress T is assumed to be in the form
T = p 1 + T .
We then apply the second law of thermodynamics to determine the class of thermodynamically consistent models based on the set Γ of variables.
Compute the time derivative ψ ˙ and replace it in (16) to obtain
ρ ( θ ψ + η ) θ ˙ ρ ρ ψ ρ ˙ p · v ρ θ ψ · ( θ ) ˙ ρ ρ ψ · ( ρ ) ˙ ρ θ ˙ ψ θ ¨ ρ ρ ˙ ψ ρ ¨ ρ D ψ · D ˙ ρ θ ψ · ( θ ) ˙ ρ ρ ψ · ( ρ ) ˙ + T · D 1 θ q · θ + θ · k = ρ θ γ 0 ,
where
· k = θ k · θ + ρ k · ρ + θ k · θ + ρ k · ρ + θ ˙ k · θ ˙ + ρ ˙ k · ρ ˙ + D k · D + θ k · θ + ρ k · ρ .
Two identities are convenient in the analysis of the inequality. They are
( θ ) ˙ = θ ˙ L T θ ,
[ ( θ ) ˙ ] j k = x j x k θ ˙ x i θ x j x k v i ( L i k x j + L i j x k ) x i θ ,
and similar with ρ in place of θ . Note that ( θ ) ˙ , ( ρ ) ˙ , θ ¨ , ρ ¨ , and D ˙ can take arbitrary (tensor or scalar) values at the point x and time t under consideration. The linearity (and arbitrariness) of these quantities imply
θ ψ = 0 , ρ ψ = 0 , θ ˙ ψ = 0 , ρ ˙ ψ = 0 , D ψ = 0 .
Now, observe that
ρ θ ψ · ( θ ) ˙ = ρ θ ψ · ( θ ˙ L T θ ) = θ · ρ θ θ ψ θ ˙ θ ˙ · ρ θ θ ψ ρ θ ψ · L T θ
and the like for ρ . Hence, the remaining inequality can be written in the form
ρ ( δ θ ψ + η ) θ ˙ ρ ( δ ρ ψ p ρ 2 ) ρ ˙ + θ { · k · ( ρ θ θ ψ θ ˙ ) · ( ρ θ ρ ψ ρ ˙ ) } + ρ ( θ θ ψ + θ ρ ψ ) · W + { T + ρ ( θ θ ψ + θ ρ ψ ) } · D 1 θ q · θ = ρ θ γ 0 ,
where δ θ and δ ρ denote generalized variational derivatives,
δ θ ψ = θ ψ θ ρ · ( ρ θ θ ψ ) , δ ρ ψ = ρ ψ θ ρ · ( ρ θ ρ ψ ) .
The linearity and arbitrariness of W in (20) imply that
θ θ ψ + ρ ρ ψ Sym .
Condition (22) holds if ψ depends on θ and ρ through | θ | , θ · ρ , | ρ | .
To within inessential divergence-free terms we can take the extra-entropy flux k in the form
k = ρ θ { θ ψ θ ˙ + ρ ψ ρ ˙ } .
Inequality (20) allows for a dependence of η on θ ˙ and p on ρ ˙ , e.g., by letting η = η ˜ c θ ˙ , p = p ˜ + d ρ ˙ with c , d > 0 [4,15]. Yet, for simplicity we neglect these dependencies for η and p, and then, it follows that
η = δ θ ψ , p = ρ 2 δ ρ ψ .
Consequently, Equation (20) reduces to
T d · D 1 θ q · θ = ρ θ γ 0 ,
where
T d = T + ρ ( θ θ ψ + θ ρ ψ ) .
Hence, the entropy production γ is amenable to the dissipative stress T d and the heat flux q . The classical Navier–Stokes–Fourier model for T d and q is just the simplest non-trivial model to account for the entropy production.
To summarize, a free energy ψ ( θ , ρ , | θ | , θ · ρ , | ρ | ) and the constitutive functions k , η , p , T d , q satisfying (23)–(25) make a non-local model thermodynamically consistent. Though the model might be more general (e.g., η , p dependent on θ ˙ , ρ ˙ ), the previous scheme allows for higher-order gradients. Indeed, we can say that the scheme is characterized by the free energy ψ and the entropy production γ .

4.1. Some Features of the Free Negentropy

It is of interest to examine some consequences of the dependence of constitutive properties on the gradients θ and ρ . The occurrence of ρ ψ / θ in the variational derivatives (21) suggests that we determine η and p in terms of the function
ζ = ρ θ ψ ,
which is the opposite of the Massieu potential [16,17]; borrowing from the terminology in [18] we can say that ζ is the Helmholtz free negentropy. We find that
η = 1 ρ [ ζ + θ d θ ζ ] , p = θ [ ζ ρ d ρ ζ ] ,
where d θ , d ρ stand for the classical variational derivatives
d θ = θ · θ , d ρ = ρ · ρ .

4.1.1. Convexity Relative to the Mass Density

Subject to the approximation of a constant temperature, the propagation of linear acoustic waves is governed by the equation
t 2 ρ = Δ p ,
where Δ denotes the Laplacian. If p = p ( ρ , ρ ) then, neglecting the nonlinear terms in ρ , we have
p = p 0 ( ρ ) + p 1 ( ρ ) Δ ρ .
The governing equation becomes
t 2 ρ = ρ p 0 Δ ρ + p 1 Δ Δ ρ .
Harmonic plane waves exp ( i ( ω t k x ) ) occur with
ω 2 = ρ p 0 k 2 + p 1 k 4
only if ρ p 0 > 0 . This insight, along with the thermodynamic interest in the dependence of p on ρ and ρ , suggests that we look for the effect of non-locality (via ρ ). Now, by (26) we have
ρ p = θ ρ ρ 2 ζ θ [ · ρ ζ + ρ ρ · ρ ζ ] .
For definiteness suppose that ζ has the form
ζ = ζ 0 ( θ , ρ ) + 1 2 f ( ρ ) | ρ | 2 .
Thus, ρ ζ = f ( ρ ) ρ , and then,
· ρ ζ + ρ ρ · ρ ζ = [ f + ρ f ] | ρ | 2 + [ f + ρ f ] Δ ρ .
so that
ρ p = θ ρ ρ 2 ζ θ { [ f + ρ f ] | ρ | 2 + [ f + ρ f ] Δ ρ } .
A further simplification arises if ζ is independent of Δ ρ , which is the case if f + ρ f = 0 . This happens if
f ( ρ ) = λ ρ ,
with λ being a constant. With this function f, it follows that
f + ρ f = λ ρ 2
and then,
ρ p = θ ρ ρ 2 ζ θ { [ f + ρ f ] | ρ | 2 = θ ρ ρ 2 ζ λ θ ρ 2 | ρ | 2 .
In light of (29), it follows that
ρ p = ρ θ ρ 2 ζ 0 ( θ , ρ ) .
Hence, if the negentropy has the form (29), then the convexity of ζ 0 , relative to the mass density ρ , implies the positive value of ρ p . This in turn occurs if the free energy ψ has the form
ψ = 1 ρ ψ 0 ( θ , ρ ) + 1 2 θ λ ρ 2 | ρ | 2 ,
where ψ 0 is convex relative to ρ .
Incidentally, in view of (27), the function p 0 ( ρ ) = ρ 2 ρ ψ , at ρ = 0 , yields
ρ p 0 = ρ [ ρ 2 ρ ( θ ρ ζ 0 ) ] = ρ θ ρ 2 ζ 0 .
Hence, in the event (29), the requirement of the positiveness of ρ p 0 ( ρ ) coincides with that of ρ p ( ρ , ρ ) .
The convexity of ζ 0 ( θ , ρ ) , relative to ρ , is connected with the convexity of the free energy ψ 0 = θ ζ / ρ . Indeed,
p = ρ 2 ρ ψ = ρ 2 ρ ( θ ζ 0 / ρ ) = θ ( ζ + ρ ρ ζ ) ,
ρ p = ρ θ ρ 2 ζ 0 .
Also, let v = 1 / ρ the specific volume and define ψ ˜ 0 ( θ , v ) = ψ 0 ( θ , ρ ) . Hence,
ρ p ρ ( ρ 2 ρ ψ 0 ) = v ( v ψ ˜ ) ( v 2 ) = v 2 ψ ˜ 0 .
Thus,
ρ θ ρ 2 ζ 0 = v 2 ψ ˜ 0 .
and the convexity of ζ 0 ( θ , ρ ) amounts to the convexity of ψ ˜ 0 ( θ , v ) .

4.1.2. Convexity Relative to the Temperature

It is worth checking the influence of the temperature gradient θ on the specific heat θ ε = θ θ η . Since
η = δ θ ψ = θ ψ + θ ρ · ( ρ θ θ ψ )
then in terms of ζ we can write
η = θ ψ + θ ρ · ( θ ζ ) .
It follows that
θ η = θ 2 ψ + 1 ρ [ · ( θ ζ ) + θ θ · ( θ ζ ) .
For definiteness let
ζ = ρ θ ψ 0 ( θ , ρ ) + 1 2 g ( θ ) | θ | 2 .
Hence, we have
· ( θ ζ ) + ρ θ · ( θ ζ ) = [ g + θ g ] | θ | 2 + [ g + θ g ] Δ θ .
Consider the particular case g + θ g = 0 , where
g ( θ ) = ν θ ,
with ν being a constant, and
g + θ g = ν θ 2 .
Consequently, Equation (31) yields
θ η = θ 2 ψ + ν ρ θ 2 | θ | 2 .
We then notice that the definition ψ = ( θ / ρ ) ζ and Function (32) result in
θ θ η = θ θ 2 ψ 0 ( η , ρ ) + ν ρ θ | θ | 2 .
The specific heat θ θ η is positive for any values of θ , ρ , and θ provided ν > 0 and ψ 0 is convex, relative to θ .

4.2. Restrictions Placed by Inequality (15)

As a comment on inequality (15), which is not assumed to be valid, we point out that the consequences of (15) on the modeling of non-local materials would be different from those of γ 0 .
For formal simplicity we restrict attention to non-local effects of temperature, and hence, let
Γ = ( θ , ρ , θ , θ ˙ , D , θ )
be the set of variables. Inequality (15) implies that
ρ ( ψ ˙ + η θ ˙ ) p · v + T · D 1 θ q · θ = ρ θ γ θ · k 0 ,
which means that ρ ( ψ ˙ + η θ ˙ ) p · v + T · D ( 1 / θ ) q · θ , and ρ θ γ θ · k have to be non-negative in addition to being equal to each other. Now, ρ θ γ θ · k 0 results in
ρ γ ( Γ ) θ k · θ ρ k · ρ θ k · θ θ ˙ k · θ ˙ D k · D θ k · θ 0 .
The linearity and arbitrariness of ρ , θ ˙ , D , θ imply that
k = k ( θ , θ ) .
Likewise, from
ρ ( ψ ˙ + η θ ˙ ) p · v + T · D ( 1 / θ ) q · θ 0 ,
namely,
ρ ( θ ψ + η ) θ ˙ ( ρ ρ ψ p / ρ ) ρ ˙ ρ θ ψ · ( θ ) ˙ ρ θ ˙ ψ θ ¨ ρ D ψ · D ˙ θ ψ · ( θ ) ˙ + T · D ( 1 / θ ) q · θ 0
it follows that
ψ = ψ ( θ , ρ ) , η = θ ψ , p = ρ 2 ρ ψ ,
along with the reduced inequality
T · D ( 1 / θ ) q · θ 0 .
Different to what follows from the CD inequality (16), here ψ is required to be independent of θ , and so is for η and p. Furthermore, T does not involve the dyadic product θ θ (and this would be the same for ρ ρ ) as happens in the previous scheme.
This example shows that the assumption (15) on the entropy inequality would be unduly restrictive relative to the correct assumption (11). Having ψ ˙ and · k in distinct inequalities is more restrictive than a single condition on ρ ψ ˙ + θ · k .

5. Entropy Production as a Constitutive Function

Back to the CD inequality (16), we now show how the entropy production γ affects, or is affected by, the constitutive equations. This is exemplified by considering the temperature-rate dependence or by models of aging materials.

5.1. Models of Rigid Heat Conductors

For simplicity consider a rigid heat conductor with
θ , θ ˙ , θ ,
as the set of variables. The CD inequality becomes
ρ ( θ ψ + η ) θ ˙ ρ θ ψ · ( θ ˙ ) 1 θ q · θ + θ · k = ρ θ γ .
Since
ρ θ ψ · ( θ ˙ ) = θ { · ( ρ θ θ ψ θ ˙ ) [ · ( ρ θ θ ψ ) ] θ ˙ }
then we have
ρ ( δ θ ψ + η ) θ ˙ + θ · ( k ρ θ θ ψ θ ˙ ) 1 θ q · θ = ρ θ γ .
Hence, it follows that
k = ρ θ θ ψ θ ˙ .
No further dependence of k is allowed, otherwise · k would include terms with an undetermined sign. A sufficient pair of relations for the validity of the remaining requirement
ρ ( δ θ ψ + η ) θ ˙ 1 θ q · θ = ρ θ γ
is
η = δ θ ψ , q = κ θ , ρ θ γ = κ θ | θ | 2 , κ > 0 .
This is what follows if γ is only assumed to be non-negative; once γ 0 is satisfied, then γ is given by 1 / ρ θ times the left-hand side.
Things are different if γ 0 is defined per se; in this event, a family of relations follow depending on the form of γ . For definiteness, if γ = c 1 ( θ ) | θ ˙ | n + c 2 ( θ ) | θ | m , then we have the relations
η = δ θ ψ c 1 | θ ˙ | n 2 θ ˙ , q = ρ θ 2 c 2 | θ | m 2 e , c 1 , c 2 > 0 , e = θ / | θ | .

5.2. Models of Aging Thermoelastic Materials

Aging properties are described by letting the constitutive parameters depend explicitly on time. This feature is now developed in connection with thermoelastic solids.
Classically (linear) thermoelastic solids are modeled by letting the second Piola stress T R R be determined by strain and temperature in the form (see [19], ch. 59)
T R R = C E + ( θ θ 0 ) M ,
where θ 0 is an equilibrium reference temperature such that T R R = 0 when E = 0 and θ = θ 0 . Furthermore, the heat flux is assumed to be given by a Fourier-type law,
q R = K R θ .
The tensors C , M , and K are the classical thermoelastic tensors. Aging thermoelastic solids are characterized by letting C , M , and K depend on time.
This suggests that we consider a thermoelastic framework where the variables are
E , θ , R θ , t ,
with the occurrence of t accounting for the aging effects. Hence,
ψ = ψ ( E , θ , R θ , t )
and similar for η , T R R , q R . The Clausius–Duhem inequality is considered in form (17), with the formal change due to the partial dependence on t. Upon computation and substitution of ψ ˙ we have
ρ R t ψ ρ R ( θ ψ + η ) θ ˙ ρ R θ ψ · ( R θ ) ˙ + ( T R R ρ R E ψ ) · E ˙ 1 θ q R · R θ = ρ R θ γ 0 ;
without any loss of generality, for formal simplicity we have assumed k R = 0 from the start. The linearity and arbitrariness of ( R θ ) ˙ , θ ˙ , E ˙ imply that
R θ ψ = 0 , η = θ ψ , T R R = ρ R E ψ
and
1 θ R θ · K R θ ρ R t ψ = ρ R θ γ .
We now restrict attention to the constitutive Equations (34) and (35). By (37) we have
ρ R ψ = Ψ ( θ ) + 1 2 E C ( t ) E + ( θ θ 0 ) M ( t ) · E .
Hence, we have
t ψ = 1 2 E · C ˙ E + ( θ θ 0 ) M ˙ · E .
Likewise, we let K depend on time, and then, the reduced inequality (38) reads
1 θ R θ · K ( t ) R θ 1 2 E · C ˙ ( t ) E ( θ θ 0 ) M ˙ ( t ) · E = ρ R θ γ 0 .
The requirement (39) can be applied by following two views. Firstly, we let γ 0 be a reminder that the left-hand side has to be non-negative and the left-hand side is just the expression of ρ R θ γ . Secondly, the left-hand side is defined in terms of γ , of course subject to γ 0 . To illustrate the two views we simplify the model by letting the solid be isotropic so that
M = m 1 , C E = 2 μ E + λ ( tr E ) 1 = 2 μ E 0 + κ ( tr E ) 1 ,
where E 0 is the deviator of E , κ = λ + 2 μ / 3 , and μ and λ are the Lamé moduli. Hence,
T R R = 2 μ E 0 + κ ( tr E ) 1 + m ( θ θ 0 ) 1 .
In stress-free conditions, we have
E 0 = 0 , κ ( tr E ) + m ( θ θ 0 ) = 0 .
Since tr E R · u and R · u is the relative variation in the volume, then m / κ is the coefficient of thermal expansion (in R ). We assume that m < 0 , so that, since κ > 0 , the body expands when the temperature increases. For isotropic solids the free energy has the form
ρ R ψ = Ψ ( θ ) + μ E 0 · E 0 + κ ( tr E ) 2 + m ( θ θ 0 ) tr E
and hence,
ρ R t ψ = μ ˙ E 0 · E 0 + κ ˙ ( tr E ) 2 + m ˙ ( θ θ 0 ) tr E .
By (40) we have ( θ θ 0 ) tr E = ( κ / m ) ( tr E ) 2 . Consequently, it follows that
ρ R t ψ = μ ˙ E 0 · E 0 + ( κ ˙ κ m ˙ / m ) ( tr E ) 2 .
Hence, the reduced inequality
K θ | R θ | 2 μ ˙ | E 0 | 2 ( κ ˙ κ m ˙ / m ) ( tr E ) 2 = ρ R θ γ 0
implies
K 0 , μ ˙ 0 , κ ˙ κ m ˙ m 0 .
In the second view, we might fix the constitutive equation for γ . For example, let
ρ R θ γ = 1 θ K | R θ | 2 + c 1 ( θ ) μ n | E 0 | 2 + κ c 2 ( θ ) | R θ | n ( tr E ) 2 ,
where K , α , and β are positive parameters, while c 1 ( θ ) 0 , c 2 ( θ ) 0 . Hence, Equation (41) implies that
μ ˙ = c 1 ( θ ) μ n , κ ˙ κ m ˙ m = κ c 2 ( θ ) | R θ | n .
Accordingly, given the constitutive function of the entropy production the entropy inequality results in the aging rate of the thermoelastic parameters. A larger set of variables might allow a more realistic evolution equation for the parameters μ , κ , and m.
In these models, we can view γ as determined by the constitutive equations, but also, the constitutive equations as determined by γ . The next section shows that for hysteretic rate-type materials the complete form of the constitutive equation is given by the assumption on the constitutive property of the entropy production.
Some comments are in order about the inequalities (42). The requirement K 0 merely shows that K can increase or decrease because of aging but anyway K remains non-negative. Instead, aging produces a decrease in μ . The bulk modulus κ is positive, and then, we can write
κ ˙ κ + m ˙ | m | 0 .
In a thermoelastic material, aging results in a decrease in μ . A joint decrease in μ , κ , and m is consistent with thermodynamics. Yet, since m < 0 , then an increase in m looks more realistic, m ˙ = | m | ˙ 0 . In this event, the consistency is expressed by
m ˙ | m | κ ˙ κ , κ ˙ 0 .
The coefficient of thermal expansion α = m / 3 κ satisfies
α ˙ = m 3 κ ( m ˙ m κ ˙ κ )
and hence,
α ˙ 0 .
Accordingly, aging results in an increase in the ratio
tr E θ θ 0
so that the solid expands more and more per increment of temperature.

6. Hysteretic Models and Entropy Production

To show the essential role of the entropy production we now consider constitutive relations for elastic–plastic bodies. We let the strain E , the Piola stress T R R , and the derivatives E ˙ , T ˙ R R be among the independent variables. The common dependence on stress and strain is connected with the hysteretic behavior; otherwise we should allow T R R to depend on E through a multi-valued function or to add an internal variable (as in [19], ch. 76). Thermal properties are also modeled, and then, we let
Ξ = ( θ , E , T R R , R θ , E ˙ , T ˙ R R )
be the set of variables. Hence, we let ψ , η , q R be functions of Ξ and assume η and q R are continuous while ψ is continuously differentiable.
Upon computation of ψ ˙ and substitution into (17) we obtain
ρ R ( θ ψ + η ) θ ˙ + ( ρ R E ψ T R R ) · E ˙ + ρ R T R R ψ · T ˙ R R + ρ R R θ ψ · R θ ˙ + E ˙ ψ · E ¨ + T ˙ R R ψ · T ¨ R R + 1 θ q R · R θ + θ R · k R = ρ R θ γ 0 .
The linearity and arbitrariness of θ ˙ , R θ ˙ , E ¨ , T ¨ R R imply that ψ is independent of R θ , E ˙ , T ˙ R R , and hence,
ψ = ψ ( θ , E , T R R ) , η = θ ψ .
Likewise, we find that
k R = k R ( θ , R θ ) ,
subject to R θ k R Skw . No skew tensor is available in the model, and hence, R θ k R = 0 . Furthermore, the isotropic character of the solid implies that k R ( θ ) has to be zero. The remaining inequality is
( ρ R E ψ T R R ) · E ˙ + ρ R T R R ψ · T ˙ R R + 1 θ q R · R θ = ρ R θ γ 0 .
If E ˙ and T ˙ R R are independent, then it follows that
T R R ψ = 0 , T R R = ρ R E ψ ,
as happens for hyperelastic materials. Yet, here we consider hysteretic materials, and hence, E ˙ and T ˙ R R are not independent. A reasonable assumption is to assume q R is independent of E ˙ and T ˙ R R . In this event, Equation (43) splits into
( ρ R E ψ T R R ) · E ˙ + ρ R T R R ψ · T ˙ R R = ρ R θ γ T 0 ,
q R · R θ = ρ R θ 2 γ q 0 .
where γ T is the value of γ when R θ = 0 , while γ q is the value of γ when E ˙ = 0 and T ˙ R R = 0 . If, instead, q R depends on E ˙ and T ˙ R R , then (44) holds along with (25), whereas (45) no longer holds.
As to (45), a Fourier-like equation for q R is allowed in the form
q R = κ ( θ , E , T R R ) R θ , γ q = κ J θ 2 | R θ | 2 , κ 0 .
Since q R = J F 1 q and R θ = θ F , then in the corresponding Eulerian description, we have
q = 1 J F F T κ θ .
Equation (44) can be solved by finding, e.g., T ˙ R R , on the assumption that T R R ψ 0 . This problem is solved by using a representation formula for tensors ([4], §A.1.3). Given any tensor A and N = A / | A | , we can represent a tensor Z in the form
Z = ( Z · N ) N + Z ,
where Z · N = 0 . If Z · N is known, say Z · N = g , while Z is unknown, then we can write
Z = g N + ( I N N ) G ,
where I is the unit fourth-order tensor and G is any second-order tensor. As a check, N · ( I N N ) G = 0 while ( I N N ) Z = Z .
Let
N = T R R ψ | T R R ψ | .
By applying (46) to (44) we obtain
T ˙ R R = ρ R θ γ T + ( ρ R E ψ T R R ) · E ˙ ρ R | T R R ψ | 2 T R R ψ + ( I T R R ψ T R R ψ | T R R ψ | 2 ) G .
Depending on the choice of G we can find various models of rate-type materials. The simplest example is obtained by letting
E ψ = 0 , ψ = 1 2 β | T R R | 2 , G = 1 β ρ R E ˙ .
In this event, Equation (47) takes the form
T ˙ R R + θ γ T β | T R R | 2 T R R = 1 β ρ R E ˙ .
This is the referential version of the Maxwell–Wiechert fluid. Indeed, the quantity β | T R R | 2 / θ γ T plays the role of relaxation time.

6.1. One-Dimensional Models

Also, with a view to experimental settings, we observe that it is worth investigating the continuum in a one-dimensional geometry. This has the advantage of simplifying the model because we can apply the Eulerian description.
Let e 1 be the longitudinal direction of the one-dimensional domain and let T 11 = σ be the only nonzero stress component. Positive values of σ denote traction, negative values denote compression. The mechanical power T · D simplifies to
T · D = σ D 11 = σ F 1 F ˙ ,
where F is the longitudinal strain, F = F 11 > 0 . Consistent with the one-dimensional model, we assume J = 1 , and hence, ρ is constant while F 22 = F 33 = 1 / F 1 / 2 . For formal simplicity we neglect heat conduction. Hence, we write the counterpart of (44) in the form
ρ ( ψ ˙ + η θ ˙ ) + σ F 1 F ˙ = ρ θ γ 0 .
Since
F 1 F ˙ = ( ln F ) ˙
then, letting λ = ln F we can write the CD inequality (48) in the form
ρ ( ψ ˙ + η θ ˙ ) + σ λ ˙ = ρ θ γ 0 .
The scalars σ and λ are Euclidean invariants. Consider the Euclidean transformation ([19], ch. 20, 21; [4], §1.9)
x * = c ( t ) + Q ( t ) x t * = t .
where Q is a rotation tensor, Q T Q = 1 . Since σ = T 11 = e 1 · T e 1 , then, under a Euclidean transformation, we have
σ * = e 1 * · T * e 1 * = ( Q e 1 ) · ( Q T Q T ) Q e 1 = e 1 · ( Q T Q T Q T Q e 1 ) = e 1 · T e 1 = σ .
Likewise, letting e ^ 1 be the first referential unit vector we have
F 11 * = e 1 * · F * e ^ 1 = Q e 1 · Q F e ^ 1 = e 1 · F e ^ 1 = F 11 .
Consequently, λ , σ , λ ˙ , and σ ˙ are Euclidean invariants and can be used as constitutive variables.
Let
θ , λ , σ , λ ˙ , σ ˙
be the set of variables. It is standard to prove that ψ has to be independent of λ ˙ and σ ˙ , and that
η = θ ψ .
Hence, it follows from (49) that
( σ ρ λ ψ ) λ ˙ ρ σ ψ σ ˙ = ρ θ γ 0 .
Since ψ = ψ ( θ , λ , σ ) , then at constant temperature ψ ˙ = λ ψ λ ˙ + σ ψ σ ˙ . Hence, along any cyclic process on [ t 1 , t 2 ] we have
t 1 t 2 σ λ ˙ d t = ρ θ t 1 t 2 γ d t 0 .
The positiveness of this integral denotes that the area within the oriented loop is positive. Thus, in a cyclic process in the λ σ plane, the curve ( λ ( t ) , σ ( t ) ) is run in the clockwise sense.
In a hysteretic process, the rate σ ˙ is associated with a
d σ d λ
that depends on the sign of λ ˙ . This would not be the case if γ = 0 or even if γ λ ˙ and γ σ ˙ . Hence, necessarily the entropy production γ has to be a constitutive function qualitatively different from the left-hand side, say a constitutive function per se. The simplest attempt is to look for a function γ proportional to | λ ˙ | . Hence, we let
ρ θ γ = γ 0 ( θ , λ , σ ) | λ ˙ | .
Thus, Equation (50) takes the form
( σ ρ λ ψ ) λ ˙ ρ σ ψ σ ˙ = γ 0 | λ ˙ | 0 .
and becomes an operative model of hysteresis once ψ and γ 0 are determined.
Analogous models are obtained by letting γ = γ σ ( θ , λ , σ ) | σ ˙ | ; here, though, we restrict our attention to inequality (51).

6.2. A Thermoelastic Hysteretic Model

For formal convenience we let Ψ = ρ ψ . Assume σ ψ 0 . Except for times where λ ˙ = 0 , we can divide (51) by σ ψ λ ˙ to obtain
σ ˙ λ ˙ = σ λ Ψ σ Ψ γ 0 σ Ψ sgn λ ˙ .
Both σ and λ are functions of X, in the referential domain, and t R . At a fixed point X in the referential domain [ 0 , L ] , σ and λ are functions of t only. Hence,
σ ˙ λ ˙ = d σ d λ .
For formal convenience we put
χ 1 = σ λ Ψ σ Ψ , χ 2 = γ 0 σ Ψ .
Both χ 1 and χ 2 are functions of λ and σ , parameterized by the temperature θ . The uniaxial stress–strain slope is then expressed in the form
d σ d λ = χ 1 + χ 2 sgn λ ˙ .
If γ 0 = 0 , then χ 2 = 0 and
d σ d λ = χ 1 ( θ , λ , σ ) ;
the slope of the curve depends also on σ and we assume that χ 1 > 0 . Since the slope d σ / d λ is anyway supposed to be non-negative, we assume
χ 1 + χ 2 sgn λ ˙ 0 .
To determine the free energy Ψ we look for a function in the form
Ψ ( λ , σ ) = L ( σ G ( λ ) ) + F ( σ ) + H ( λ ) ,
where L , G , F , H are differentiable functions parameterized by θ . Substitution of σ Ψ and λ Φ yields χ 1 and χ 2 in the forms
χ 1 = σ + L ( σ G ( λ ) ) G ( λ ) H ( λ ) L ( σ G ( λ ) ) + F ( σ ) , χ 2 = γ 0 L ( σ G ( λ ) ) + F ( σ ) .
The function χ 1 is the elastic differential stiffness. Hence, we let
χ 1 = g ( λ ) > 0 .
Accordingly, we obtain the requirement
σ F ( σ ) g ( λ ) L ( σ G ( λ ) ) [ g ( λ ) G ( λ ) ] = H ( λ ) .
This condition is satisfied by letting F ( σ ) = 0 and
g ( λ ) G ( λ ) = α , H ( λ ) = G ( λ ) , L ( σ G ( λ ) ) = 1 α ( σ G ( λ ) ) ,
where α 0 is a suitable parameter for the model. Hence, we have
χ 2 = γ 0 α 1 σ G ( λ ) .
Furthermore,
Ψ ( λ , σ ) = 1 2 α [ σ G ( λ ) ] 2 + H ( λ ) , H ( λ ) = G ( λ ) .
To sum up, the whole model is determined by
g ( λ ) , α , γ 0 ( λ , σ )
and
G ( λ ) = 0 λ [ g ( y ) α ] d y , H ( λ ) = G ( λ ) .
For definiteness we now establish some examples of hysteretic solids. The corresponding loops are obtained by letting λ = Λ sin ( ω t ) , and then, solving the system
{ λ ˙ = ω Λ cos ω t , σ ˙ = ( d σ / d λ ) λ ˙ .
Since the model is rate-independent, the loops are not affected by the value of the angular frequency ω .
  • Plastic flow with asymptotic strength.
We start with a model based on a constant elastic differential stiffness α . Let G ( λ ) = 0 so that Ψ = σ 2 / 2 α and χ 1 = α . The hysteretic function γ 0 is taken in the form
γ 0 ( σ ) = σ 2 σ u , σ u > 0 .
Hence, the whole differential stiffness is
d σ d λ = α σ u ( σ u σ sgn λ ˙ ) .
In this event it follows that
d d λ ( σ σ u ) 2 = 2 α σ u ( σ σ u ) ( σ u σ sgn λ ˙ ) ,
where
d d λ ( σ σ u ) 2 = ( 2 α / σ u ) { ( σ σ u ) 2 if λ ˙ > 0 , σ 2 σ u 2 if λ ˙ < 0 .
Thus,
d d λ ( σ σ u ) 2 0
and the hysteresis loops are confined to the strip | σ | σ u .
The hysteresis loops in Figure 1 are obtained by solving the system (54) and using (55) with α = 1 and σ u = 1.5 .
  • Plastic flow with a nonlinear elastic function.
Let
G ( λ ) = tanh ( κ λ ) , κ > 0
and
γ 0 = 1 β { σ tanh ( κ λ ) } 2 , β > 0 .
Hence, by (53) and (52) it follows that
Ψ ( σ , λ ) = 1 2 α [ σ tanh ( κ λ ) ] 2 + 1 κ ln [ cosh ( κ λ ) ] ,
χ 1 = κ [ 1 tanh 2 ( κ λ ) ] , χ 2 = α [ σ tanh ( κ λ ) ] sgn λ ˙ .
The differential stiffness d σ / d λ can be given the form
d σ d λ = κ [ 1 tanh 2 ( κ λ ) ] + α β { β [ σ tanh ( κ λ ) ] sgn λ ˙ } .
The hysteresis loops in Figure 2 are obtained by solving the system (54) and using (56) with α = 1 , β = 1.5 , and κ = 4 . They well describe the hysteretic responses of lateral loads with respect to lateral displacements in a typical medium-rise building model (see, e.g., [20]).
  • Elastic–plastic model.
We now describe a solid undergoing linear behavior in the elastic regime. Hence, we let G ( λ ) = 0 and obtain
Ψ = σ 2 2 α , χ 1 = α .
To characterize γ we consider two stress levels, σ u > σ y > 0 , and assume hysteretic effects are confined to the region | σ | [ σ y , σ u ] in the form
γ 0 ( σ , sgn λ ˙ ) = { ( | σ | σ y ) / ( σ u σ y ) | σ | if | σ | σ y and σ λ ˙ > 0 , 0 if | σ | < σ y or | σ | σ y , σ λ ˙ < 0 .
The differential stiffness takes the form
d σ d λ = { α ( σ u | σ | ) / ( σ u σ y ) if | σ | σ y and σ λ ˙ > 0 , 0 if | σ | < σ y or | σ | σ y , σ λ ˙ < 0 .
Figure 3 shows the hysteresis loop obtained by solving the system (54) and using (57) with the parameters α = 1 , σ u = 2.5 , and σ y = 1.5 . Within the region σ y | σ | σ u , the material behaves elastically during unloading and plastically during loading.
As exemplified by the previous models of plastic materials, the hysteretic properties are represented in simple and direct ways by an appropriate form of the entropy production γ as a constitutive function. Analogous properties hold in the modeling of ferroics. In addition to the conceptual character of γ as a constitutive function, these examples prove the experimental evidence of the reported method on the exploitation of the entropy inequality.

7. Conclusions

This paper deals with the mathematical formulation and the use of the second law of thermodynamics in continuum physics. Conceptually the second law states that the rate of entropy in any sub-region of the continuum is greater than the external entropy supply. This amounts to the assumption that the (rate of) entropy production is non-negative. Mathematically, this leads to a procedure for the selection of physically admissible constitutive properties [21]. In the Coleman–Noll formulation, the entropy flux, say j , is q / θ , while j is a constitutive function in the Müller formulation. It is an important point of the present formulation that the entropy production γ is also a constitutive function (Section 3).
The constitutive property of γ is shown to have remarkable consequences on the whole thermodynamic scheme. Though quite uncommon in the literature, Section 5 shows that the aging properties of non-dissipative materials, that is, terms related to energy functions, result in positive entropy productions. Next, as is shown in Section 6, the occurrence of γ as a constitutive function is essential in the modeling of thermodynamically consistent hysteretic materials. In particular, this is shown for elastic–plastic materials, though the analogue can be performed for magnetic or electric hysteresis [22,23,24].
It is a further result, shown in Section 6, that a representation formula allows a complete description of the consequences of the second law inequality. This greater generality is apparent when the constitutive equations involve vectors or tensors in rate-type equations.
From the standpoint of the mathematical modeling, the role played by the entropy production as a constitutive function is decisive, at least in the case of hysteretic materials. This is so because hysteresis exhibits a different behavior depending on the sign of a time derivative (namely, in loading and unloading). If, e.g., the variables are λ , λ ˙ , then the term λ ˙ ψ λ ¨ in the entropy inequality leads to λ ˙ ψ = 0 . Hence, the dependence on λ ˙ happens through λ ψ λ ˙ . The dependence on the sign of λ ˙ is then allowed by letting γ depend on λ ˙ through | λ ˙ | . By this approach, the thermodynamic requirement results in a hysteretic Duhem-like model [25].
The procedure of Section 6, based on the constitutive function of entropy production, is likely to apply to hysteresis processes [26,27]. To our mind, a rate equation similar to (50) might well describe the time evolution of entropy–temperature loops.

Author Contributions

Investigation: C.G. and A.M. All authors have contributed substantially and equally to the work reported. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Acknowledgments

The research leading to this work has been developed under the auspices of INDAM-GNFM.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Coleman, B.D.; Noll, W. The thermodynamics of elastic materials with heat conduction and viscosity. Arch. Ration. Mech. Anal. 1963, 13, 167–178. [Google Scholar] [CrossRef]
  2. Müller, I. On the entropy inequality. Arch. Ration. Mech. Anal. 1967, 26, 118–141. [Google Scholar] [CrossRef]
  3. Müller, I. Thermodynamics; Pitman: London, UK, 1985. [Google Scholar]
  4. Morro, A.; Giorgi, C. Mathematical Modelling of Continuum Physics; Birkhäuser: Cham, Switzerland, 2023. [Google Scholar]
  5. Zemansky, M.W.; Dittman, R.H. Heat and Thermodynamics; McGraw-Hill: New York, NY, USA, 1997; Chapter 8. [Google Scholar]
  6. Green, A.E.; Laws, N. On the entropy production inequality. Arch. Ration. Mech. Anal. 1972, 45, 47–53. [Google Scholar] [CrossRef]
  7. Green, A.E.; Naghdi, P.M. On thermodynamics and the nature of the second law. Proc. R. Soc. Lond. A 1977, 357, 253–270. [Google Scholar]
  8. Maugin, G.A. Internal variables and dissipative structures. J. Non-Equilib. Thermodyn. 1990, 15, 173–192. [Google Scholar] [CrossRef]
  9. Maugin, G.A. On the thermomechanics of continuous media with diffusion and/or weak nonlinearity. Arch. Appl. Mech. 2006, 75, 217–249. [Google Scholar] [CrossRef]
  10. Alber, H.-D.; Hutter, K.; Tsakmakis, C. Nonconventional thermodynamics, indeterminate couple stress elasticity and heat conduction. Cont. Mech. Thermodyn. 2016, 28, 699–719. [Google Scholar] [CrossRef]
  11. Alber, H.-D.; Broese, C.; Tsakmakis, C.; Beskos, D.E. Non-convenctional thermodynamics and models of gradient elasticity. Entropy 2018, 20, 179–192. [Google Scholar] [CrossRef]
  12. Dunn, J.E.; Serrin, J. On the thermodynamics of interstitial working. Arch. Ration. Mech. Anal. 1985, 88, 95–133. [Google Scholar] [CrossRef]
  13. Fabrizio, M.; Giorgi, C.; Morro, A. A thermodynamic approach to non-isothermal phase-field evolution in continuum physics. Physica D 2006, 214, 144–156. [Google Scholar] [CrossRef]
  14. Bertram, A. Compendium on Gradient Materials, 1st ed.; Springer: Cham, Switzerland, 2023. [Google Scholar]
  15. Bogy, D.B.; Naghdi, P.M. On heat conduction and wave propagation in rigid solids. J. Math. Phys. 1970, 11, 917–923. [Google Scholar] [CrossRef]
  16. Massieu, M.F. Sur les fonctions caractéristiques des divers fluides. C.R. Acad. Sci. 1869, 69, 858–862. [Google Scholar]
  17. Massieu, M.F. Addition au precedent memoire sur les fonctions caractéristiques. C.R. Acad. Sci. 1869, 69, 1057–1061. [Google Scholar]
  18. Brillouin, L. Negentropy principle of information. J. Appl. Phys. 1953, 24, 1152–1163. [Google Scholar] [CrossRef]
  19. Gurtin, M.E.; Fried, E.; Anand, L. The Mechanics and Thermodynamics of Continua; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  20. Grubišić, M.; Ivošević, J.; Grubišić, A. Reliability analysis of reinforced concrete frame by finite element method with implicit limit state functions. Buildings 2019, 9, 119. [Google Scholar] [CrossRef]
  21. Morro, A.; Giorgi, C. Techniques for the thermodynamic consistency of constitutive equations. Thermo 2023, 3, 260–276. [Google Scholar] [CrossRef]
  22. Giorgi, C.; Morro, A. A thermodynamic approach to rate-type models of elastic-plastic materials. J. Elast. 2021, 147, 113–148. [Google Scholar] [CrossRef]
  23. Giorgi, C.; Morro, A. A thermodynamic approach to rate-type models in deformable ferroelectrics. Cont. Mech. Thermodyn. 2021, 33, 727–747. [Google Scholar] [CrossRef]
  24. Giorgi, C.; Morro, A. Thermodynamically-consistent modeling of ferromagnetic hysteresis. Materials 2023, 16, 2882. [Google Scholar] [CrossRef] [PubMed]
  25. Ikhouane, F. A survey of the hysteretic Duhem model. Arch. Computat. Methods Eng. 2018, 25, 965–1002. [Google Scholar] [CrossRef]
  26. Hess, T.; Maier, L.M.; Bachmann, N.; Corhan, P.; Schäfer-Welsen, O.; Wöllenstein, J.; Bartholomé, K. Thermal hysteresis and its impact on the efficiency of first-order caloric materials. J. Appl. Phys. 2020, 127, 075103. [Google Scholar] [CrossRef]
  27. Gottschall, T.; Skokov, K.P.; Burriel, R.; Gutfleisch, O. On the S(T) diagram of magnetocaloric materials with first-order transition: Kinetic and cyclic effects of Heusler alloys. Acta Mater. 2016, 107, 1–8. [Google Scholar] [CrossRef]
Figure 1. Plastic model with asymptotic bounds σ = ± σ u (dashed): hysteresis loops (solid) with Λ = 1 , 2 , 3 and starting values ( λ 0 , σ 0 ) = ( 0 , 0 ) , ( 0 , 0.1 ) .
Figure 1. Plastic model with asymptotic bounds σ = ± σ u (dashed): hysteresis loops (solid) with Λ = 1 , 2 , 3 and starting values ( λ 0 , σ 0 ) = ( 0 , 0 ) , ( 0 , 0.1 ) .
Thermo 04 00015 g001
Figure 2. Plastic model with nonlinear elastic function: hysteresis loops (solid) with Λ = 2 , 5 and starting values ( λ 0 , σ 0 ) = ( 0 , 0 ) , ( 0 , 0.1 ) ; asymptotic bounds σ = ± ( 1 + β ) (dashed).
Figure 2. Plastic model with nonlinear elastic function: hysteresis loops (solid) with Λ = 2 , 5 and starting values ( λ 0 , σ 0 ) = ( 0 , 0 ) , ( 0 , 0.1 ) ; asymptotic bounds σ = ± ( 1 + β ) (dashed).
Thermo 04 00015 g002
Figure 3. Elastic–plastic model with asymptotic bounds σ = ± σ u (dashed) and yielding thresholds σ = ± σ y (short dashed); hysteresis loops (solid) with amplitude Λ = 4 and starting value ( λ 0 , σ 0 ) = ( 0 , 0 ) .
Figure 3. Elastic–plastic model with asymptotic bounds σ = ± σ u (dashed) and yielding thresholds σ = ± σ y (short dashed); hysteresis loops (solid) with amplitude Λ = 4 and starting value ( λ 0 , σ 0 ) = ( 0 , 0 ) .
Thermo 04 00015 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Giorgi, C.; Morro, A. On the Second Law of Thermodynamics in Continuum Physics. Thermo 2024, 4, 273-294. https://doi.org/10.3390/thermo4020015

AMA Style

Giorgi C, Morro A. On the Second Law of Thermodynamics in Continuum Physics. Thermo. 2024; 4(2):273-294. https://doi.org/10.3390/thermo4020015

Chicago/Turabian Style

Giorgi, Claudio, and Angelo Morro. 2024. "On the Second Law of Thermodynamics in Continuum Physics" Thermo 4, no. 2: 273-294. https://doi.org/10.3390/thermo4020015

APA Style

Giorgi, C., & Morro, A. (2024). On the Second Law of Thermodynamics in Continuum Physics. Thermo, 4(2), 273-294. https://doi.org/10.3390/thermo4020015

Article Metrics

Back to TopTop