Thermodynamically Consistent Evolution Equations in Continuum Mechanics

1. Introduction

The constitutive properties of materials are often described through evolution equations whereby the time derivative of suitable fields is eventually related to appropriate fields and their derivatives of various order. The dependence on the space and time derivatives is a way of accounting for non-local properties in space and time. Relative to the description of response of the material in terms of space and time integrals, differential evolution equations show the advantage of greater flexibility.

Starting from rheological models [1,2], examples of differential equations are settled in the mechanical domain. For example, in connection with viscoelastic materials a simple equation is that associated with the Maxwell model. In the three-dimensional setting, we can write the stress–strain relation in the form,

$\dot{T} = E \dot{ϵ} - λ T,$

(1)

where $T$ is the stress tensor and $ϵ$ is the (infinitesimal) strain tensor. Likewise, the Kelvin–Voigt model of fluids is expressed by

$\dot{D} = a_{1} D + a_{2} T,$

(2)

where $T$ is the viscous stress and $D$ is the stretching tensor. Here, $E, λ$ and $a_{1}, a_{2}$ are constants; we might view these coefficients as functions of temperature and mass density. Other models may involve higher-order time differentiation, as seen in the Burgers material, which involves second-order time derivatives (see, e.g., [3], ch. 4).

As for thermal models, it is worth mentioning the Maxwell–Cattaneo equation [4,5,6,7]

$\dot{q} + λ \nabla θ + ν q = 0$

(3)

for the temperature $θ$ and the heat flux vector $q$ . In nanoscale systems continuum models might be applicable provided non-local and/or rate properties are properly described. The evolution is often modelled by the Guyer–Krumhansl [8,9] equation, namely

$\dot{q} + λ \nabla θ + ν q - ζ (Δ q + \nabla (\nabla \cdot q)) = 0,$

(4)

where $q$ is the heat flux vector, ∇ is the gradient operator, $Δ$ is the Laplacian, and the superposed dot denotes the material time derivative. Equation (4) is then a generalization of (18) with the dependence of $\dot{q}$ on the higher-order spatial derivatives $Δ q$ and $\nabla (\nabla \cdot q)$ .

Within continuum mechanics, the constitutive properties are required to be objective, meaning they remain form-invariant under the set of Euclidean transformations, SO(3). However, Equations (1)–(4) are not objective, merely because the time derivative is not objective. Hence, these equations deserve an objective form. Furthermore, the physical admissibility of any constitutive equation, like, e.g., (1)–(4), is required to be consistent with the second law of thermodynamics. Consequently, we need a procedure to ascertain the physical admissibility of constitutive equations.

The physical admissibility is by now a well-known topic; however, it is performed in various ways. A formal statement about the thermodynamic consistency traces back to Coleman and Noll [10]. Next, Müller [11] generalized the statement by letting the entropy flux be a constitutive quantity to be determined. Lately, a further generalization has been accomplished by regarding also the entropy production as a constitutive quantity per se ([2], ch. 3).

The purpose of this paper is to establish a systematic procedure to ascertain whether given constitutive equations are thermodynamically consistent. While this is classical in the literature, albeit with technical variants, this paper shows how the outcome is in general a set of constitutive equations that possibly includes a selected equation, such as Equations (1)–(4). Mathematically, the procedure is made operative by using a representation formula that determines a wide class of equations once a set of variables is chosen depending on the given particular equation.

This paper revisits the thermodynamic approach and emphasizes the operative character of the second law in connection with physically admissible models. Next, some significant examples of evolution equations are investigated, both in the Lagrangian and in the Eulerian description.

2. Notation and Balance Equations

We denote by $Ω$ , the time-dependent, three-dimensional, region occupied by the body under consideration. 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 \in R$ . The symbol ∇ denotes the gradient, with respect to $x$ , while $\nabla \cdot$ is the divergence operator. For any pair of vectors $u, w$ , or tensors $A, B$ , the notations $u \cdot w$ and $A \cdot B$ denote the inner product. Cartesian coordinates are used and then, in the suffix notation, $u \cdot w = u_{i} w_{i}$ , $A \cdot B = A_{i j} B_{i j}$ , the summation over repeated indices being 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 $Ω \times R$ , we have $\dot{f} = \partial_{t} f + (v \cdot \nabla) f$ . The symbol $L$ denotes the velocity gradient, $L_{i j} = \partial_{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.

2.1. Balance Equations

The conservation of mass is expressed locally by the continuity equation

$\dot{ρ} + ρ \nabla \cdot v = 0 .$

The balance of linear momentum leads to the equation of motion

$ρ \dot{v} = \nabla \cdot T + ρ b,$

(5)

where ${(\nabla \cdot T)}_{i} = \partial_{x_{j}} T_{i j}$ . The balance of angular momentum implies that $T = T^{T}$ . To obtain the balance of energy, we consider the specific energy density $\frac{1}{2} v^{2} + ε$ (per unit mass) and the power $(T n) \cdot v - q \cdot n$ per unit area, with unit external normal $n$ , and $b \cdot v + r$ per unit mass; $ε, q, r$ are quantities of a non-mechanical character. In light of the equation of motion (5), it follows that the internal energy $ε$ is subject to the local equation

$ρ \dot{ε} = T \cdot D - \nabla \cdot q + ρ r .$

(6)

Let $η$ be the specific entropy density, $j$ the entropy flux, and $r / θ$ the entropy supply. The balance of entropy leads to

$ρ \dot{η} + \nabla \cdot j - \frac{ρ r}{θ} = ρ γ,$

where $γ$ is the (rate of) entropy production. It is assumed that $γ \geq 0$ , which implies the increase of entropy for isolated systems ( $j = 0, r = 0$ ).

We assume $b$ and r are arbitrarily given time-dependent fields on $Ω \times R$ . A process is the set $P = (ρ, v, T, ε, η, θ, q, j, γ)$ entering the balance equations. Since the number of components in a process is larger than the number of balance equations, then a corresponding set of constitutive equations are needed to make the two numbers equal. Hence, the balance equations provide $\dot{ρ}, \dot{v}, \dot{ε}$ ; the assumed arbitrariness of $b (x, t)$ and $r (x, t)$ allows us to view $\dot{v}$ , $\dot{ε}$ , and their derivatives as arbitrary.

2.2. Second Law of Thermodynamics

The second law of thermodynamics is assumed as follows.

Postulate.For every process $P$ admissible in a body, the inequality
$ρ \dot{η} + \nabla \cdot j - \frac{ρ r}{θ} = ρ γ \geq 0$
(7)
is valid at any internal point.

Henceforth, we refer to Equation (7) as the Clausius–Duhem (CD) inequality. It is convenient to let

$j = \frac{q}{θ} + k, ψ = ε - θ η;$

$k$ is said to be the extra-entropy flux and $ψ$ the Helmholtz free energy. Substitution of $\nabla \cdot q - ρ r$ from (6) makes (7) in the form

$- ρ (\dot{ψ} + η \dot{θ}) + T \cdot D + θ \nabla \cdot k - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ \geq 0 .$

(8)

This is the thermodynamic scheme for fields in the Eulerian description, where $Ω \times R$ is the space–time domain of the process. For later use, we now establish the corresponding relations in the space–time domain $R \times R$ , where $R$ is the region $Ω$ at some convenient time. To fix ideas, let $R = Ω_{t = 0}$ .

2.3. Lagrangian Description

Let

X

be the position vector of a point in

R

and hence,

x (X, t)

is the motion. We denote by

\nabla R = \partial_{X}

the gradient operator in

R

. We let

$F = \nabla_{R} x, F_{i K} = \partial_{X_{K}} x_{i},$

be the deformation gradient and hence, for any differentiable function $f (x, t) = f (x (X, t), t)$ , we have

Define

and let

be the Green–Lagrange strain tensor. Hence,

$\dot{F} = L F, \dot{E} = F^{T} D F .$

Let

$T_{R R} = J F^{- 1} T F^{- T}, q_{R} = J q F^{- T};$

$T_{R R}$ is the second Piola (or Piola–Kirchhoff) stress while $q_{R}$ is the heat flux in $R$ . We can prove that ([2], ch. 1; [12], §25)

$J T \cdot D = T_{R R} \cdot \dot{E}, J q \cdot \nabla θ = q_{R} \cdot \nabla R θ,$

Furthermore, letting $q_{R} = J q F^{- T}$ m, we have

$J \nabla \cdot k = \nabla R \cdot k_{R} .$

Since $J ρ = ρ_{R}$ is the mass density in the reference configuration $R$ , then J times the CD inequality (8) yields

$- ρ_{R} (\dot{ψ} + η \dot{θ}) + T_{R R} \cdot \dot{E} + θ \nabla R \cdot k_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ = ρ_{R} θ γ \geq 0 .$

(9)

Two points characterize the present approach. Firstly, the entropy production $γ$ is given by a constitutive equation per se as is the case for the flux $k_{R}$ . Secondly, the constitutive equation for ${\dot{q}}_{R}$ is not assumed a priori but is derived as a consequence of the other constitutive equations and the CD inequality. In particular, this shows how $k_{R}$ and $γ$ influence the equation for $q_{R}$ .

2.4. Representation Formula

Depending on the set of variables, the CD inequality (8) or (9) has the form

in the unknown $z$ with $w, z$ as vectors or tensors. For example, if $ψ$ depends on $q_{R}$ , then (9) eventually has the form (10) for the unknown $z = {\dot{q}}_{R}$ , with

$w = ρ_{R} \partial_{q_{R}} ψ, f = - ρ_{R} θ γ + θ \nabla \cdot k_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ .$

Notice that, for any unit vector

n

, a vector

z

can be represented in the form

$z = (z \cdot n) n + z_{⊥},$

where $(z \cdot n) n$ is the longitudinal part (in the direction $n$ ) while $z_{⊥}$ is the transverse part, $z_{⊥} \cdot n = 0$ . If only $z \cdot n$ is known, say $z \cdot n = g$ , then

Any vector $z_{⊥}$ perpendicular to $n$ can be represented in the form

$z_{⊥} = (1 - n \otimes n) ξ,$

for any vector $ξ$ , where $1$ is the identity tensor. This is so in that

$[(1 - n \otimes n) ξ] \cdot n = ξ \cdot n - (ξ \cdot n) n \cdot n = 0 .$

Back to (10), we have

$z = \frac{z \cdot w}{{| w |}^{2}} w + (1 - \frac{w \otimes w}{{| w |}^{2}}) ξ .$

(11)

w

and

z

are tensors, say

then we repeat the formal steps to conclude that

$Z = \frac{Z \cdot U}{{| U |}^{2}} U + (I - \frac{U \otimes U}{{| U |}^{2}}) Ξ,$

(12)

where $I$ is the fourth-order identity tensor and $Ξ$ is any second-order tensor.

3. Thermodynamic Restrictions

Our purpose is to determine evolution equations of

q_{R}

consistent with thermodynamics and hence, we look for a constitutive equation of

{\dot{q}}_{R}

. Possible nonlocal effects are described by a suitable dependence on the spatial derivatives. For definiteness, we direct our attention to a dependence on

\nabla R q_{R}

and

\nabla R \nabla R q_{R}

. A dependence on the strain

E

is allowed to account for the elastic properties of the solid. Hence, we let

$Γ = (θ, E, q_{R}, \nabla R θ, \nabla R q_{R}, \nabla R \nabla R q_{R})$

be the set of independent variables and assume $η, T_{R R}, k_{R}, γ$ , and ${\dot{q}}_{R}$ are continuous functions of $Γ$ , while $ψ$ is continuously differentiable.

Compute

\dot{ψ}

and substitute in (9) to obtain

$\begin{matrix} - ρ_{R} (\partial_{θ} ψ + η) \dot{θ} + (T_{R R} - ρ_{R} \partial_{E} ψ) \cdot \dot{E} - ρ_{R} \partial_{\nabla R θ} ψ \cdot \nabla R \dot{θ} - ρ_{R} \partial_{\nabla R q_{R}} ψ \cdot \nabla R {\dot{q}}_{R} \\ - ρ_{R} \partial_{\nabla R \nabla R q_{R}} ψ \cdot \nabla R \nabla R {\dot{q}}_{R} - ρ_{R} \partial_{q_{R}} ψ \cdot {\dot{q}}_{R} + θ \nabla_{R} \cdot k_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ = ρ_{R} θ γ \geq 0 . \end{matrix}$

(13)

The linearity and arbitrariness of $\nabla R \dot{θ}, \nabla R {\dot{q}}_{R}, \nabla R \nabla R {\dot{q}}_{R}, \dot{θ}, \dot{E}$ , imply that

$\partial_{\nabla R θ} ψ = 0, \partial_{\nabla R q_{R}} ψ = 0, \partial_{\nabla R \nabla R q_{R}} ψ = 0,$

$η = - \partial_{θ} ψ, T_{R R} = \partial_{E} ψ .$

(14)

Hence, (13) simplifies to

$- ρ_{R} \partial_{q_{R}} ψ \cdot {\dot{q}}_{R} + θ \nabla_{R} \cdot k_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ = ρ_{R} θ γ \geq 0 .$

(15)

Inequality (15) requires the compatibility among the constitutive functions $ψ, K_{R}, γ$ along with the non-negative valuedness of $γ$ . For the present purpose, we now use Equation (15) as an equation in the unknown ${\dot{q}}_{R}$ .

The representation formula (11) is now applied to obtain

{\dot{q}}_{R}

from (10). Assume that

\partial_{q_{R}} ψ \neq 0

and let

n = \partial_{q_{R}} ψ / | \partial_{q_{R}} ψ |

. By (11), it follows that

${\dot{q}}_{R} = \frac{- ρ_{R} θ γ + θ \nabla R \cdot k_{R} - (1 / θ) q_{R} \cdot \nabla R θ}{ρ_{R} {| \partial_{q_{R}} ψ |}^{2}} \partial_{q_{R}} ψ + (1 - \frac{\partial_{q_{R}} ψ \otimes \partial_{q_{R}} ψ}{| \partial_{q_{R}} {ψ |}^{2}}) ξ .$

(16)

Equation (16) is the general thermodynamic requirement on ${\dot{q}}_{R}$ under the assumption $\partial_{q_{R}} ψ \neq 0$ with variables $Γ$ .

The representation (16) and the arbitrariness of the field

ξ

show that any transverse term relative to

\partial_{q_{R}} ψ

is allowed, irrespective of its connection with the entropy production

γ

or the extra-entropy flux

k_{R}

. This property is apparent from the CD inequality (9), where

{\dot{q}}_{R}

occurs through the component along

\partial_{q_{R}} ψ

. Physical examples are given in the next section.

4. Remarkable Models of Heat Conduction

Models are now derived by taking different assumptions on $γ, k_{R}, ξ$ .

4.1. Non-Dissipative Heat Conduction

Let

γ = 0, k_{R} = 0, ξ = 0

. Hence, Equation (16) simplifies to

${\dot{q}}_{R} = - \frac{q_{R} \cdot \nabla R θ}{ρ_{R} θ {| \partial_{q_{R}} ψ |}^{2}} \partial_{q_{R}} ψ .$

The non-dissipative character is characterized by the assumption $γ = 0$ . Now, as $θ$ and $E$ are constants, we have

$\dot{ψ} = \partial_{q_{R}} ψ \cdot {\dot{q}}_{R} = - \frac{1}{θ ρ_{R}} q_{R} \cdot \nabla R θ .$

The classical entropy production $- 1 / ρ_{R} θ^{2}$ times $q_{R} \cdot \nabla R θ$ , which enters the heat conduction inequality, produces a corresponding rate of the free energy in this model.

4.2. Fourier Model

Let

k_{R} = 0

and

ξ = 0

. If, further, we let

ρ_{R} θ^{2} γ = - q_{R} \cdot \nabla R θ

, then we find

{\dot{q}}_{R} = 0

, which means that

q_{R}

is stationary. Now, if we let

$γ = \nabla R θ \cdot K (θ, E) \nabla R θ,$

where $K$ is positive definite, then it follows that

$q_{R} = - ρ_{R} θ^{2} K \nabla R θ .$

Thus, the Fourier law follows with heat conductivity tensor $ρ_{R} θ^{2} K$ .

4.3. Some Evolutionary Models

We now look for some equations involving

{\dot{q}}_{R}

and relate them to models that appear in the literature. Hence, for definiteness, hereafter, we let

$ψ = ψ_{0} (θ, E) + \frac{1}{2} α {| q_{R} |}^{2},$

where $α$ is possibly dependent on $θ$ and $E$ . Equation (16) then is written in the form

${\dot{q}}_{R} = \frac{- ρ_{R} θ γ + θ \nabla \cdot k_{R} - (1 / θ) q_{R} \cdot \nabla R θ}{ρ_{R} α {| q_{R} |}^{2}} q_{R} + (1 - \frac{\partial_{q_{R}} ψ \otimes \partial_{q_{R}} ψ}{| \partial_{q_{R}} {ψ |}^{2}}) ξ .$

(17)

Maxwell–Cattaneo Equation

Let

k_{R} = 0

. Hence, Equation (17) simplifies to

${\dot{q}}_{R} = \frac{- ρ_{R} θ γ - (1 / θ) q_{R} \cdot \nabla R θ}{ρ_{R} α {| q_{R} |}^{2}} q_{R} + (1 - \frac{q_{R} \otimes q_{R}}{| q_{R} |^{2}}) ξ .$

Assume $ξ = β \nabla R θ$ . If

it follows

${\dot{q}}_{R} + \frac{θ γ}{α | q_{R} |^{2}} q_{R} = - \frac{1}{ρ_{R} θ α} \nabla R θ .$

(18)

Equation (18) is just the Maxwell–Cattaneo (MC for short) equation with relaxation time $τ$ and heat conductivity $κ$ , given by

$τ = \frac{α | q_{R} |^{2}}{θ γ}, κ = \frac{| q_{R} |^{2}}{ρ_{R} θ^{2} γ} .$

The positive value of the entropy dissipation $γ$ implies the positive value of $κ$ . Instead, the sign of $α$ is left undetermined by thermodynamics; it is customary to assume $τ > 0$ , which follows from the assumption $α > 0$ .

4.4. Cross-Couplings between Deformation and Heat Conduction

Cross-coupling terms are modelled by the joint dependence on

E

and

q_{R}

. For simplicity, let

where $ζ$ is any scalar function of $Γ$ . Hence, Equation (16) takes the form

${\dot{q}}_{R} = \frac{- ρ_{R} θ γ + θ \nabla \cdot k_{R} - (1 / θ) q_{R} \cdot \nabla R θ - ζ ρ_{R} α q_{R} \cdot E q_{R}}{ρ_{R} α {| q_{R} |}^{2}} q_{R} + ζ E q_{R} .$

(19)

The selection of appropriate functions $γ \geq 0$ and $k_{R}$ completes the model.

4.5. Second-Order Non-Local Models

Non-local models are established in terms of the dependence on the spatial derivatives of

q_{R}

. A model often considered in the literature is named after Guyer and Krumhansl and involves second-order derivatives [8,9]. Indeed, we consider the generalized (Lagrangian) form

${\dot{q}}_{R} + λ \nabla R θ + ν q_{R} = c_{1} Δ_{R} q_{R} + c_{2} \nabla R (\nabla R \cdot q_{R});$

(20)

the coefficients $λ, ν, c_{1}, c_{2}$ are expected to be positive-valued functions of temperature.

With the purpose of obtaining an evolution equation like (20) and recalling the assumption

ξ = β \nabla R θ

made in connection with the MC Equation (18), we let

$ξ = β \nabla R θ + c_{1} Δ_{R} q_{R} + c_{2} \nabla R (\nabla R \cdot q_{R}),$

where $β, c_{1}, c_{2}$ are functions of $θ$ and $E$ . Letting $β = - 1 / ρ_{R} θ α$ , we can write (17) in the form

$\begin{matrix} {\dot{q}}_{R} = \frac{- ρ_{R} θ γ + θ \nabla \cdot k_{R}}{ρ_{R} α {| q_{R} |}^{2}} q_{R} - \frac{1}{ρ_{R} θ α} \nabla R θ + c_{1} Δ_{R} q_{R} + c_{2} \nabla R (\nabla R \cdot q_{R}) \\ - \frac{c_{1}}{| q_{R} |^{2}} (q_{R} \cdot Δ_{R} q_{R}) q_{R} - \frac{c_{2}}{| q_{R} |^{2}} [(q_{R} \cdot \nabla R) \nabla_{R} q_{R}] q_{R} . \end{matrix}$

We can get Equation (20) by the identifications $λ = 1 / ρ_{R} θ α$ and

$ν = \frac{1}{ρ_{R} α {| q_{R} |}^{2}} ρ_{R} θ γ - θ \nabla \cdot k_{R} + ρ_{R} α c_{1} q \cdot Δ_{R} q_{R} + ρ_{R} α c_{2} (q_{R} \cdot \nabla R) \nabla R \cdot q_{R} .$

(21)

The form (20) is then obtained, albeit with $ν$ as a function of the second-order derivatives $\nabla R \nabla R q_{R}$ . Owing to the occurrence of $γ$ and $k_{R}$ , we ask whether we can eliminate the dependence of $ν$ on the derivatives of $θ$ and $q_{R}$ .

Let

$h = \frac{ρ_{R} α c_{1}}{θ}, l = \frac{ρ_{R} α c_{2}}{θ} .$

We can then write Equation (21) in the form

$\frac{ρ_{R} α ν}{θ} {| q_{R} |}^{2} = ρ_{R} γ - \nabla \cdot k_{R} + h q \cdot Δ_{R} q_{R} + l (q_{R} \cdot \nabla R) \nabla R \cdot q_{R} .$

Notice that

$h q_{R} \cdot Δ_{R} q_{R} = \nabla R \cdot (h \nabla R \frac{1}{2} | q_{R} |^{2}) - [\nabla R (h q_{R})] \cdot [\nabla R q_{R}],$

$l (q_{R} \cdot \nabla R) \nabla R \cdot q_{R} = \nabla R \cdot (l q_{R} \nabla R \cdot q_{R}) - [\nabla R \cdot (l q_{R})] \nabla R \cdot q_{R} .$

Hence, we select

$k_{R} = h \nabla R \frac{1}{2} {| q_{R} |}^{2} - l q_{R} \nabla R \cdot q_{R}$

(22)

to find

$\frac{ρ_{R} α {| q_{R} |}^{2} ν}{θ} + [\nabla R (h q_{R})] \cdot [\nabla R q_{R}] + [\nabla R \cdot (l q_{R})] \nabla R \cdot q_{R} = ρ_{R} γ .$

If h and l are constants, then

$[\nabla R (h q_{R})] \cdot [\nabla R q_{R}] + [\nabla R \cdot (l q_{R})] \nabla R \cdot q_{R} = h | \nabla R q_{R} |^{2} + l {(\nabla R \cdot q_{R})}^{2} \geq 0,$

provided only that $h, l > 0$ and hence, $α c_{1}, α c_{2} > 0$ .

The generalized form (20) is thermodynamically consistent if h and l are positive constants. In that case the extra-entropy flux has the form (22), while the entropy production is

$γ = \frac{α | q_{R} |^{2} ν}{θ} + \frac{α c_{1}}{θ} (\nabla R q_{R}) \cdot (\nabla R q_{R}) + \frac{α c_{2}}{θ} {(\nabla R \cdot q_{R})}^{2} .$

5. Higher-Order Differential Equations

Experiments on heat propagation in graphene [13] show that the results are best fitted with the differential equation

$\frac{1}{α} \partial_{t} θ + \frac{τ_{q}}{α} \partial_{t}^{2} θ = Δ θ + τ_{θ} \partial_{t} Δ θ,$

(23)

for the temperature $θ (x, t)$ . There are various physical models leading to (23). The simplest one is to assume the conductor is rigid and let $c = \partial_{θ} ε$ be the specific heat so that the balance of energy simplifies to

$ρ c \partial_{t} θ = - \nabla \cdot q,$

with a zero energy supply. The evolution of the heat flux is assumed to be governed by the MC equation with a higher-order correction so that

$τ_{q} \dot{q} + q = - λ \nabla θ - λ τ_{θ} \partial_{t} \nabla θ .$

(24)

Assuming c is constant, upon some rearrangements, we find Equation (23) with $α = 1 / ρ c$ . Owing to the formal occurrence of two times $τ_{q}, τ_{θ}$ , Equation (23) is also referred to as a heat equation with two phase lags. Other physical motivations of (23) are given in [14].

We now examine the thermodynamic consistency of heat equations of the form (24). For generality, we let the body be deformable and hence, we follow again the Lagrangian description. We then consider equations of the form

$τ_{q} {\dot{q}}_{R} + q_{R} = - λ \nabla R θ - λ τ_{θ} \nabla R \dot{θ} .$

(25)

Owing to the term $\nabla R \dot{θ}$ in Equation (25), we let both $\nabla R θ$ and $\dot{θ}$ be among the variables. Hence, we let

$Γ = (θ, E, q_{R}, \nabla R θ, \dot{θ})$

be the set of variables and $ψ, η, T_{R R}, k_{R}, γ$ and ${\dot{q}}_{R}$ be given by constitutive functions. Computation of $\dot{ψ}$ and substitution in the CD inequality yields

$- ρ_{R} (\partial_{θ} ψ + η) \dot{θ} - ρ_{R} \partial_{q_{R}} ψ \cdot {\dot{q}}_{R} - ρ_{R} \partial_{\nabla R θ} ψ \cdot \nabla R \dot{θ} - ρ_{R} \partial_{\dot{θ}} ψ \ddot{θ} + θ \nabla R \cdot k_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ = ρ_{R} θ γ .$

(26)

The linearity and arbitrariness of $\ddot{θ}$ implies that

The occurrence of $\nabla R \dot{θ}$ affects the properties of the entropy $η$ . Hence, we define

${\hat{q}}_{R} = {\dot{q}}_{R} + d (θ, E) \nabla R \dot{θ} .$

In (25), it is $d = λ τ_{θ} / τ_{q}$ ; consistent with (25), we assume ${\dot{q}}_{R}$ , and hence, ${\hat{q}}_{R}$ , are independent of $\dot{θ}$ . Substitution of ${\dot{q}}_{R}$ in (26) and division by $θ$ result in

$- \frac{ρ_{R}}{θ} (\partial_{θ} ψ + η) \dot{θ} - \frac{ρ_{R}}{θ} \partial_{q_{R}} ψ \cdot {\hat{q}}_{R} - \frac{ρ_{R}}{θ} y \cdot \nabla R \dot{θ} + \nabla R \cdot k_{R} - \frac{1}{θ^{2}} q_{R} \cdot \nabla R θ = ρ_{R} γ,$

(27)

where

$y = d \partial_{q_{R}} ψ + \partial_{\nabla R θ} ψ .$

Now, observe that

$- \frac{ρ_{R}}{θ} y \cdot \nabla R \dot{θ} + \nabla R \cdot k_{R} = \nabla R \cdot (k_{R} - \frac{ρ_{R}}{θ} y \dot{θ}) + [\nabla R \cdot (\frac{ρ_{R}}{θ} y] \dot{θ} .$

Hence, we let

Consequently, we can write Equation (27) in the form

$- \frac{ρ_{R}}{θ} (δ_{θ} ψ + η) \dot{θ} - \frac{ρ_{R}}{θ} \partial_{q_{R}} ψ \cdot {\hat{q}}_{R} - \frac{1}{θ^{2}} q_{R} \cdot \nabla R θ = ρ_{R} γ,$

(28)

where

$δ_{θ} ψ : = \partial_{θ} ψ - \frac{θ}{ρ_{R}} \nabla R \cdot [\frac{ρ_{R}}{θ} (d \partial_{q_{R}} ψ + \partial_{\nabla R θ} ψ)] .$

If $ρ_{R}$ and $θ$ are uniform, that is $\nabla R ρ_{R} = 0, \nabla R θ = 0$ , then

$δ_{θ} ψ : = \partial_{θ} ψ - \nabla R \cdot (d \partial_{q_{R}} ψ + \partial_{\nabla R θ} ψ) .$

If, further, $d = 0$ , then $δ_{θ} ψ$ equals the variational derivative of $ψ$ with respect to $θ$ . Accordingly, we can say that $δ_{θ} ψ$ is a generalized variational derivative of $ψ$ .

While

\hat{q}

is independent of

\dot{θ}

, both

η

and

γ

can depend on

\dot{θ}

. If we let

$η = η_{0} (θ, E, q_{R}, \nabla R θ) + η_{1} (Γ) \dot{θ},$

then we find that

$η_{0} = - δ_{θ} ψ, η_{1} \leq 0, ρ_{R} γ = ρ_{R} γ_{0} + ρ_{R} γ_{1}, γ_{1} = - \frac{1}{θ} η_{1} {\dot{θ}}^{2},$

where $γ_{0} \geq 0$ is the value of the entropy production $γ$ at $\dot{θ} = 0$ .

Equation (28) reduces to

$- ρ_{R} \partial_{q_{R}} ψ \cdot {\hat{q}}_{R} - \frac{1}{θ} q_{R} \cdot \nabla R θ = ρ_{R} θ γ_{0} .$

(29)

Hence, using the representation (11) with $n = \partial_{q_{R}} ψ / | \partial_{q_{R}} ψ |$ , we can determine ${\hat{q}}_{R}$ in the form

${\hat{q}}_{R} = - \frac{(1 / θ) q_{R} \cdot \nabla R θ + ρ_{R} θ γ_{0}}{ρ_{R} {| \partial_{q_{R}} ψ |}^{2}} \partial_{q_{R}} ψ + (1 - \frac{\partial_{q_{R}} ψ \otimes \partial_{q_{R}} ψ}{| \partial_{q_{R}} {ψ |}^{2}}) ξ .$

We assume

$\partial_{q_{R}} ψ = α (θ, E, \nabla R θ) q_{R},$

whence

$ψ = ψ_{0} (θ, E, \nabla R θ) + \frac{1}{2} α (θ, E, \nabla R θ) {| q_{R} |}^{2} .$

The selection

leads to

${\hat{q}}_{R} = - \frac{ρ_{R} θ γ_{0}}{ρ_{R} α {| \partial_{q_{R}} |}^{2}} q_{R} - \frac{q_{R} \cdot \nabla R θ}{θ ρ_{R} α {| q_{R} |}^{2}} q_{R} + β \nabla R θ - β \frac{q_{R} \cdot \nabla R θ}{| q_{R} |^{2}} q_{R} .$

The particular case

results in

${\hat{q}}_{R} = - \frac{θ γ_{0}}{α | q_{R} |^{2}} q_{R} - \frac{1}{θ ρ_{R} α} \nabla R θ .$

This in turn leads to the rate equation for $q_{R}$ in the form

${\dot{q}}_{R} + \frac{θ γ_{0}}{α | q_{R} |^{2}} q_{R} = - \frac{1}{θ ρ_{R} α} \nabla R θ - d \nabla R \dot{θ},$

whence

$λ = \frac{| q_{R} |^{2}}{θ^{2} γ_{0} ρ_{R}} > 0, τ_{q} = \frac{α | q_{R} |^{2}}{θ γ_{0}}, τ_{θ} = θ ρ_{R} α d .$

In stationary conditions, we have

$q_{R} = - \frac{| q_{R} |^{2}}{θ^{2} ρ_{R} γ_{0}} \nabla R θ,$

which ascribes to

$K = \frac{| q_{R} |^{2}}{θ^{2} ρ_{R} γ_{0}}$

the meaning of heat conductivity. Since $γ_{0} > 0$ the conductivity K proves to be positive, as expected.

6. Rate Equations in the Eulerian Description

In accordance with the objectivity principle ([15], p. 35; [16], §41), the constitutive equations must be form-invariant under Euclidean transformations, namely under the map

x \mapsto x^{*}

such that

where the vector $c$ and the tensor $Q$ are time-dependent and $Q$ is orthogonal, $Q^{T} Q = 1$ . Consequently, if a constitutive equation is expressed by a rate equation then the time derivatives must be objective in that the objective time derivative of a vector transforms as a vector, and the same is true for tensors.

A natural question arises since there are infinitely many objective derivatives and the non-uniqueness remains, even though we require the thermodynamic consistency to hold. We now examine rate-type models. For formal simplicity, we consider separately the occurrence of rate effects for the stress and the heat flux.

6.1. The Navier–Stokes–Voigt Fluid

Roughly, the NSV (Navier–Stokes–Voigt) fluid model is the counterpart of the Kelvin–Voigt solid, whereby the stress is a linear combination of the strain and the strain-rate. Consistently, the NSV system of equations modifies the Navier–Stokes equations by addition of a term of form $ν Δ \partial_{t} v$ to the acceleration $\dot{v}$ . Really, this is the approximation induced by replacing the Navier–Stokes stress $μ D + λ (tr D) 1$ with the linear combination $μ D + λ (tr D) 1 + \hat{μ} \dot{D} + \hat{λ} (tr \dot{D}) 1$ . However, the total time derivative, along with the partial time derivative, is not objective and hence, the idea of the linear combination has to be revisited.

Here, for generality, we let the fluid be compressible, though usually in the literature the fluid is incompressible (see, e.g., [17,18,19]). Hence, we let

$T = - p (θ, ρ) 1 + T,$

where p is the thermodynamic pressure, assumed to be a function of the temperature $θ$ and the mass density $ρ$ . The stress $T$ is expected to be related to $D$ and the rate of $D$ . Furthermore, we keep modelling heat conduction through a rate equation. Hence, we let

$Λ = (θ, ρ, T, D, q, \dot{θ}, \nabla θ)$

be the set of variables for the constitutive equations of $ψ, η$ , and the rates of $D$ and $q$ . Indeed, the rate equations for $D$ and $q$ have to be established via objective time derivatives.

Compute the time derivative of

ψ (Λ)

and substitute in the CD inequality (8) to obtain

$\begin{matrix} - ρ (\partial_{θ} ψ + η) \dot{θ} + (ρ^{2} \partial_{ρ} ψ - p) \nabla \cdot v - ρ \partial_{D} ψ \cdot \dot{D} - ρ \partial_{q} ψ \cdot \dot{q} - ρ \partial_{\nabla θ} ψ \cdot \dot{(\nabla θ)} \\ - ρ \partial_{\dot{θ}} ψ \ddot{θ} + T \cdot D - \frac{1}{θ} q \cdot \nabla θ + θ \nabla \cdot k = ρ θ γ . \end{matrix}$

(30)

The linearity and arbitrariness of $\ddot{θ}$ imply that

By the thermodynamic character of the pressure, we conclude that

Thus, Equation (30) simplifies to

$\begin{matrix} - \frac{ρ}{θ} (\partial_{θ} ψ + η) \dot{θ} - \frac{ρ}{θ} \partial_{D} ψ \cdot \dot{D} - \frac{ρ}{θ} \partial_{q} ψ \cdot \dot{q} - \frac{ρ}{θ} \partial_{\nabla θ} ψ \cdot \dot{(\nabla θ)} \\ + \frac{1}{θ} T \cdot D - \frac{1}{θ^{2}} q \cdot \nabla θ + \nabla \cdot k = ρ γ, \end{matrix}$

(31)

where, for technical purposes, we have divided throughout by $θ$ .

Notice that

$\dot{(\nabla θ)} = \nabla \dot{θ} - L^{T} \nabla θ$

and then

$- \frac{ρ}{θ} \partial_{\nabla θ} ψ \cdot \dot{(\nabla θ)} + \nabla \cdot k = \nabla \cdot (k - \frac{ρ}{θ} \partial_{\nabla θ} ψ \dot{θ}) + \nabla \cdot (\frac{ρ}{θ} \partial_{\nabla θ} ψ) \dot{θ} .$

Consequently, we let

$k = \frac{ρ}{θ} \partial_{\nabla θ} ψ \dot{θ} .$

Hence, $θ$ times the remaining part of (31) results in

$- ρ (δ_{θ} ψ + η) \dot{θ} - ρ \partial_{D} ψ \cdot \dot{D} - ρ \partial_{q} ψ \cdot \dot{q} + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ,$

(32)

where $δ_{θ} ψ$ denotes the variational derivative,

$δ_{θ} ψ = \partial_{θ} ψ - \frac{θ}{ρ} \nabla \cdot (\frac{ρ}{θ} \partial_{\nabla θ} ψ) .$

By Equation (32), the set

Λ

of variables allows for a rich variety of properties. Formally, one of the rates, e.g.,

\dot{θ}

, is a function of

\dot{D}, \dot{q}, γ

and the remaining quantities. For definiteness, we look for separate effects or, mathematically, sufficient conditions to satisfy Equation (32).

Firstly, let

η

, as well as

\dot{D}, \dot{q}, γ

be independent of

\dot{θ}

. Hence, Equation (32) holds only if

and

$- ρ \partial_{D} ψ \cdot \dot{D} - ρ \partial_{q} ψ \cdot \dot{q} + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ .$

(33)

To proceed with the analysis of (33), it is worth observing that a rate equation of the form

$\dot{D} = F (θ, ρ, T, D, q, \nabla θ)$

does not comply with the objectivity principle in that $\dot{D}$ is not form invariant. We then look for a recourse to objective derivatives. To fix ideas, consider the objective derivative

$\overset{⋄}{D} = \dot{D} - W D - D W^{T} - 2 c_{1} D D + c_{2} (\nabla \cdot v) D, \overset{⋄}{q} = \dot{q} - W Q - c_{1} D q + c_{2} (\nabla \cdot v) q .$

If $c_{1} = 1, c_{2} = 1$ , then $\overset{⋄}{D}, \overset{⋄}{q}$ are the Truesdell derivatives of $D$ and $q$ . If, instead, $c_{1} = 0, c_{2} = 0$ , then $\overset{⋄}{D}$ and $\overset{⋄}{q}$ coincide with the corotational derivatives,

$\overset{\circ}{D} = \dot{D} - W D - D W^{T}, \overset{\circ}{q} = \dot{q} - W q .$

The other objective derivatives have as common part the corotational derivative [2]. Replacing

$\dot{D} = \overset{⋄}{D} + W D + D W^{T} + 2 D D - (\nabla \cdot v) D, \dot{q} = \overset{⋄}{q} + W Q + D q - (\nabla \cdot v) q,$

we can write Equation (33) in the form

$\begin{matrix} - ρ \partial_{D} ψ \cdot [\overset{⋄}{D} + W D + D W^{T} + 2 c_{1} D D - c_{2} (\nabla \cdot v) D] \\ - ρ \partial_{q} ψ \cdot [\overset{⋄}{q} + W q + c_{1} D q - c_{2} (\nabla \cdot v) q] + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ . \end{matrix}$

(34)

The linearity and arbitrariness of the spin $W$ in (34) imply that

$\partial_{D} ψ D \in Sym, \partial_{q} ψ \otimes q \in Sym,$

whence

$\partial_{D} ψ = α_{1} D, \partial_{q} ψ = α_{2} q .$

(35)

Hence, (34) simplifies to

$\begin{matrix} - ρ α_{1} D \cdot [\overset{⋄}{D} + 2 c_{1} D D - c_{2} (\nabla \cdot v) D] + T \cdot D \\ - ρ α_{2} q \cdot [\overset{⋄}{q} + c_{1} D q - c_{2} (\nabla \cdot v) q] - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ . \end{matrix}$

(36)

Furthermore, in light of (35), we have

$ψ = ψ_{0} (θ, ρ, T, \nabla θ) + \frac{1}{2} α_{1} {| D |}^{2} + \frac{1}{2} α_{2} {| q |}^{2} .$

(37)

Two significant instances are taken from the literature.

6.2. Non-Heat Conducting NSV Fluid

This approximation traces back to Oskolkov, though in the particular case of incompressible fluids [20]. We let

q = 0

and notice that

$D \cdot D = {| D |}^{2}, D \cdot (D D) = tr (D D D) .$

Consequently, Equation (36) simplifies to

$- ρ α D \cdot \overset{⋄}{D} - 2 ρ α c_{1} tr (D D D) + ρ α c_{2} {| D |}^{2} \nabla \cdot v + T \cdot D = ρ θ γ .$

(38)

We now apply the representation formula (11) for the unknown $\overset{⋄}{D}$ with

Hence,

$\begin{matrix} \overset{⋄}{D} & = \frac{ρ α D \cdot \overset{⋄}{D}}{{ρ α | D |}^{2}} D + [I - \frac{D \otimes D}{{| D |}^{2}}] Ξ \\ = - \frac{θ γ}{{α | D |}^{2}} D + \frac{T \cdot D}{{ρ α | D |}^{2}} D - \frac{2 c_{1} tr (D D D)}{{| D |}^{2}} D + c_{2} (\nabla \cdot v) D + [I - \frac{D \otimes D}{{| D |}^{2}}] Ξ, \end{matrix}$

where $Ξ$ is any tensor function of $θ, ρ, T, D, \nabla θ$ . If $c_{1} = 0, c_{2} = 0$ , then $\overset{⋄}{D}$ reduces to the corotational derivative; thus, we have

$\overset{\circ}{D} = - \frac{θ γ}{{α | D |}^{2}} D + \frac{T \cdot D}{{ρ α | D |}^{2}} D + [I - \frac{D \otimes D}{{| D |}^{2}}] Ξ .$

If $Ξ = β T$ , then

$\overset{\circ}{D} = - \frac{θ γ}{{α | D |}^{2}} D + \frac{T \cdot D}{{ρ α | D |}^{2}} D + β T - β \frac{T \cdot D}{{| D |}^{2}} D .$

The choice $β = 1 / ρ α$ results in

$T = \frac{ρ θ γ}{{| D |}^{2}} D + ρ α \overset{\circ}{D} .$

Borrowing from the Navier–Stokes theory we identify the coefficient of $D$ with the shear viscosity $μ$ and hence, we can write

which gives

$T = μ (D + τ \overset{\circ}{D}),$

where $τ = ρ α / μ$ represents the relaxation time.

In light of (35), the free energy has the form

$ψ = ψ_{0} (θ, ρ, T, \nabla θ) + \frac{1}{2} α {| D |}^{2} .$

Hence, also $η$ and $ε = ψ + η θ$ have additive terms proportional to ${| D |}^{2}$ .

Remark 1.

Two theories of NSV fluids have been developed with different schemes; however, both of them use a kinetic energy density proportional to

{| \nabla v |}^{2}

[21,22]. We then see a similarity with the free energy ψ, comprising a term proportional to

{| D |}^{2}

6.3. NSV Fluid with Non-Fourier Conduction

Back to Equation (34), we now focus our attention on the corotational derivative as the objective derivative (and hence, we let

c_{1} = 0, c_{2} = 0

). Furthermore, the model of heat conduction is sought by looking for a rate-type equation; thus, we obtain, in particular, a Maxwell-Cattaneo equation. In this way, a generalization is established of the model applied by Straughan [23].

By (36) and (37), it follows

$- ρ α_{1} D \cdot \overset{\circ}{D} - ρ α_{2} q \cdot \overset{\circ}{q} + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ .$

(39)

For simplicity, cross couplings (e.g., $D q$ terms) are set aside and hence, the entropy production is assumed in the separate form

$γ = γ_{T} (θ, ρ, T, D) + γ_{q} (θ, ρ, q, \nabla θ) .$

Both $γ_{T}$ and $γ_{q}$ are non-negative, in that $γ_{T}$ is the value of $γ$ when $T = 0, D = 0$ , and $γ_{q}$ is the value of $γ$ when $q = 0$ . Hence, (39) splits into two equations,

$- ρ α_{1} D \cdot \overset{\circ}{D} + T \cdot D = ρ θ γ,$

(40)

$- ρ α_{2} q \cdot \overset{\circ}{q} - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ_{q} .$

(41)

Equation (40) has just been considered in the previous subsection. As for (41), we apply the representation formula with $n = q / | q |$ to obtain

$\overset{\circ}{q} = \frac{ρ θ γ_{q} + (1 / θ) q \cdot \nabla θ}{- ρ α_{2} {| q |}^{2}} q + [1 - \frac{q \otimes q}{{| q |}^{2}}] ξ .$

(42)

If we assume

ξ = 0

, then it follows

$\overset{\circ}{q} = \frac{ρ θ γ_{q} + (1 / θ) q \cdot \nabla θ}{- ρ α_{2} {| q |}^{2}} q .$

If no action on the system happens ( $ξ = 0$ ), then $\overset{\circ}{q}$ is directed along $q$ . Furthermore, in stationary conditions, $\overset{\circ}{q} = 0$ , it follows

$- \frac{1}{θ} q \cdot \nabla θ = ρ θ γ_{q} \geq 0;$

that is, the standard thermodynamic condition on the heat flux when only the temperature gradient is considered to produce the heat flux.

If we assume

ξ = β \nabla θ

, then it follows

$\overset{\circ}{q} = \frac{ρ θ γ_{q} + (1 / θ) q \cdot \nabla θ}{- ρ α_{2} {| q |}^{2}} q + β \nabla θ - β \frac{q \cdot \nabla θ}{{| q |}^{2}} q .$

Hence, letting $β = - 1 / ρ θ α_{2}$ , we have

$\overset{\circ}{q} + \frac{θ γ_{q}}{α_{2} {| q |}^{2}} q = - \frac{1}{ρ θ α_{2}} \nabla θ .$

(43)

Equation (43) has the form of the MC equation where the time derivative is in the corotational (and hence objective) form. Indeed, we can view

$τ = \frac{α_{2} {| q |}^{2}}{θ γ_{q}}, κ = \frac{{| q |}^{2}}{ρ θ^{2} γ_{q}},$

as the relaxation time and the heat conductivity. The positive value of the entropy production $γ_{q}$ yields $κ > 0$ . Any assumption $γ_{q} = g (θ, ρ) {| q |}^{2}$ makes $κ$ a function of $θ$ and $ρ$ . The assumption $α_{2} > 0$ implies $τ > 0$ and hence, guarantees the boundedness of the solution $q$ in time. This procedure in turn shows how the MC equation is merely a particular case of the family of thermodynamically consistent evolution equations $\overset{\circ}{q} = \hat{q} (θ, ρ, q, \nabla θ)$ with $ξ = β \nabla θ$ .

7. Evolution Equations, Second Law, and Objective Derivatives

By analogy with the MC equation for the heat flux, in fluid-dynamics, evolution equations for the stress tensor are considered in the Eulerian description. In this context, we contemplate the selection of objective derivatives and the restrictions placed by the second law, as well as the objectivity of thermodynamically consistent evolution equations.

The stress

T

is represented in the form

T = - p 1 + T

, with p representing the thermodynamic pressure. Describe a thermo-viscous fluid through the variables

$θ, ρ, T, q, D, \nabla θ,$

for the constitutive functions $ψ, η, p, \dot{T}, \dot{q}, γ$ . The Clausius–Duhem inequality (with $k = 0$ ) reads

$\begin{matrix} - ρ (\partial_{θ} ψ + η) \dot{θ} - ρ \partial_{ρ} ψ \dot{ρ} - ρ \partial_{T} ψ \cdot \dot{T} - ρ \partial_{q} ψ \cdot \dot{q} - ρ \partial_{D} ψ \cdot \dot{D} - ρ \partial_{\nabla θ} ψ \cdot \dot{(\nabla θ)} \\ - p \nabla \cdot v + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ . \end{matrix}$

The linearity and arbitrariness of $\dot{D}, \dot{(\nabla θ)},$ and $\dot{θ}$ imply

$\partial_{D} ψ = 0, \partial_{\nabla θ} = 0, η = - \partial_{θ} ψ .$

Substitute $\dot{ρ}$ with $- ρ \nabla \cdot v$ and define the thermodynamic pressure by

The remaining condition is

$- ρ \partial_{T} ψ \cdot \dot{T} - ρ \partial_{q} ψ \cdot \dot{q} + T \cdot D - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ .$

For simplicity, we neglect cross-coupling terms (i.e.,

\dot{q}

is independent of

D

and

T

, while

\dot{T}

is independent of

q

and

\nabla θ

) and

$γ = γ_{T} (θ, ρ, T, D) + γ_{q} (θ, ρ, q, \nabla θ);$

both $γ_{T}$ and $γ_{q}$ are non-negative. Hence, we have

$- ρ \partial_{T} ψ \cdot \dot{T} + T \cdot D = ρ θ γ_{T} .$

(44)

$- ρ \partial_{q} ψ \cdot \dot{q} - \frac{1}{θ} q \cdot \nabla θ = ρ θ γ_{q} .$

(45)

By applying the representation formula to (44), it follows

$\dot{T} = \frac{T \cdot D - ρ θ γ_{T}}{ρ | \partial_{T} {ψ |}^{2}} \partial_{T} ψ + [I - \frac{\partial_{T} ψ \otimes \partial_{T} ψ}{| \partial_{T} {ψ |}^{2}}] Ξ,$

(46)

where $Ξ$ is any tensor function of $θ, ρ, T, D$ . Equation (46) is not objective because $\dot{T}$ is not. Thus, the second law per se does not provide objective evolution equations. Yet, we show that thermodynamic consistency and objectivity are compatible.

Consider the objective time derivative

$\overset{⋄}{T} = \dot{T} - W T - T W^{T} + c_{1} (D T + T D) + c_{2} (\nabla \cdot v) T,$

with generic coefficients $c_{1}, c_{2}$ , and examine the thermodynamic consistency of the sought constitutive function $\overset{⋄}{T} (θ, ρ, T, D)$ . Substitution of $\dot{T}$ in (44) results in

$- ρ \partial_{T} ψ \cdot [\overset{⋄}{T} + W T + T W^{T} - c_{1} (D T + T D) - c_{2} (\nabla \cdot v) T] + T \cdot D = ρ θ γ_{T} .$

(47)

Since

$\partial_{T} ψ \cdot [W T + T W^{T}] = W \cdot [\partial_{T} ψ T^{T} + {(\partial_{T} ψ)}^{T} T],$

then the linearity and arbitrariness of $W$ in (47) imply

$\partial_{T} ψ T^{T} + {(\partial_{T} ψ)}^{T} T \in Sym;$

(48)

the condition (48) holds if $ψ$ depends on $T$ through the invariants $tr T$ or $| T |$ . Hence, Equation (47) reduces to

$- ρ \partial_{T} ψ \cdot [\overset{⋄}{T} - c_{1} (D T + T D) - c_{2} (\nabla \cdot v) T] + T \cdot D = ρ θ γ_{T} .$

For definiteness, let

$ψ = ψ_{0} (θ, ρ) + \frac{1}{2} λ {| T |}^{2} .$

The representation formula for $\overset{⋄}{T}$ yields

$\overset{⋄}{T} = - \frac{θ γ_{T}}{{λ | T |}^{2}} T - \frac{c_{1} (D T + T D) \cdot T + c_{2} (\nabla \cdot v) {| T |}^{2}}{{| T |}^{2}} T + \frac{T \cdot D}{{ρ λ | T |}^{2}} T + [I - \frac{T \otimes T}{{| T |}^{2}}] Ξ .$

(49)

Equation (49) shows various terms for the objective derivative $\overset{⋄}{T}$ , namely the term induced by the entropy production, the terms due to the chosen derivative, the term induced by the power $T \cdot D$ and the term associated with the transverse tensor $[I - T \otimes T / | T |^{2}] Ξ$ . Among the particular cases, we may consider the tensor $Ξ = β D$ . Letting $β = 1 / ρ λ$ , we find that

$\overset{⋄}{T} + \frac{θ γ_{T}}{{λ | T |}^{2}} T = - \frac{c_{1} (D T + T D) \cdot T + c_{2} (\nabla \cdot v) {| T |}^{2}}{{| T |}^{2}} T + \frac{1}{λ ρ} D .$

If, further, $c_{1} = 0, c_{2} = 0$ , then we have

$\overset{⋄}{T} + \frac{θ γ_{T}}{{λ | T |}^{2}} T = \frac{1}{λ ρ} D .$

(50)

Equation (50) can be viewed as the tensor form of the Maxwell–Wiechert fluid model ([2], sec. 6.1.4). Accordingly,

$τ = \frac{{λ | T |}^{2}}{θ γ_{T}},$

can be viewed as the relaxation time. In stationary conditions ( $\overset{⋄}{T} = 0$ ), we have

$T = 2 μ D, μ = \frac{{| T |}^{2}}{2 θ γ_{T} ρ},$

thus ascribing to $μ$ the role of viscosity coefficient. The assumption $γ_{T} > 0$ implies that $μ > 0$ .

Remarks about Rate Equations and Thermodynamic Restrictions

Depending on the set of variables, it may happen that the CD inequality has the form (as, e.g., in (44))

$f (θ, ρ, T, D, \dot{T}) = 0 .$

If $\partial_{\dot{T}} f \neq 0$ , then by the representation formula, we can derive $\dot{T}$ , as in (46), but this rate equation would be non-objective. Hence, objectivity is an additional requirement. In fact, this requirement is realized by replacing the time derivative with an objective derivative, say $\overset{⋄}{T}$ , and then deriving the (objective) equation for $\overset{⋄}{T}$ . Owing to the non-uniqueness of objective derivatives, this shows that thermodynamics and objectivity lead to a class of rate-type models rather than to a single equation.

The need of objective constitutive equations and the non-uniqueness of the objective variables are well established in the literature. For instance, in connection with fluids of differential type [24,25], the dependence on the variables

L, \dot{L}, \ddot{L}, . . .

is replaced from the start with that on the Rivlin–Ericksen tensors

A_{1}, A_{2}, A_{3}, . . .

, e.g.,

$\hat{w} (L, \dot{L}, \ddot{L}, . . .) = \tilde{w} (A_{1}, A_{2}, A_{3}, . . .) .$

The literature also shows another way of modelling through objective derivatives. The starting assertion in [26] is that a scalar subjected to an objective rate is equal to its conventional time rate. Hence, one might expect the validity of the equation

$\dot{G} = \dot{a} \cdot b + a \cdot \dot{b} = \overset{⋄}{a} \cdot b + a \cdot \overset{⋄}{b} = \overset{⋄}{G}$

for the two vector functions $a, b$ and $G = a \cdot b$ ; in [26] $\overset{⋄}{G}$ is the logarithmic derivative. Yet, this equation in general is not true. For simplicity, let

$\overset{⋄}{a} = \dot{a} - W a + D a .$

We have

$\overset{⋄}{a} \cdot b + a \cdot \overset{⋄}{b} = \dot{a} \cdot b + a \cdot \dot{b} - (b \otimes a + a \otimes b) \cdot W + (b \otimes a + a \otimes b) \cdot D .$

Since $b \otimes a + a \otimes b \in Sym$ , then $(b \otimes a + a \otimes b) \cdot W = 0$ , but $(b \otimes a + a \otimes b) \cdot D$ need not be zero. Accordingly, in general, the most convenient procedure is to start from the time derivatives and then to replace them through the selected objective derivatives.

8. Conclusions

This paper addresses the modelling of material behaviour in continuous bodies. In general, non-locality effects, in space and time, have to be described by involving time and space dependence of the pertinent fields. In this connection, the paper looks at the modelling of constitutive properties in terms of differential equations as relations among space and time derivatives of various orders.

The occurrence of time derivatives has to comply with the objectivity principle. Objective equations are obtained in two ways. In the Lagrangian description some pertinent fields (e.g., $T_{R R}, q_{R}, E$ , as well as the densities $ψ, ε$ ) are invariant under the set of Euclidean transformations. Hence, their time derivatives are invariant and objective. Furthermore, the recourse to the Lagrangian description naturally allows models with finite deformations. Instead, in the Eulerian description, the standard fields (e.g., $T, q, D$ ) are not invariant and their derivatives are not objective. Objectivity is obtained by using objective time derivatives; the corotational derivative is the simplest one and seemingly the most appropriate.

Physically admissible models are required to comply with the second law of thermodynamics. Here, the postulate of the second law is eventually expressed in terms of the Clausius–Duhem inequality. Two features have to be emphasized with respect to the tradition. Firstly, the (non-negative) entropy production is viewed as a constitutive function per se and not merely deduced from the other constitutive properties. Secondly, the Clausius–Duhem inequality is viewed as a constraint on the pertinent fields and it is solved with respect to the interested vector or tensor by using a representation formula (Section 2; see also [2]). This allows for qualitatively different constitutive equations, as is emphasized in Section 4, Section 5 and Section 6. Rather than confining to the admissibility of a single equation, the procedure leads to the admissibility of a class of models characterized by the entropy production

γ

and the transverse vector or tensor (determined by

ξ

Ξ

Section 3, Section 4 and Section 5 are devoted to models of heat conduction, within Lagrangian descriptions, while Section 6 develops some models of the Navier–Stokes–Voigt fluid, within Eulerian descriptions. In connection with thermo-viscous fluids, Section 7 gives a systematic approach to the modelling of evolution equations in the Eulerian description. It is shown that the thermodynamic consistency is compatible with both objective and non-objective evolution equations, which emphasizes that objectivity is an independent principle of continuum mechanics.

Source link

Angelo Morro www.mdpi.com

Greenberg News

Thermodynamically Consistent Evolution Equations in Continuum Mechanics

1. Introduction

2. Notation and Balance Equations

2.1. Balance Equations

2.2. Second Law of Thermodynamics

2.3. Lagrangian Description

2.4. Representation Formula

3. Thermodynamic Restrictions

4. Remarkable Models of Heat Conduction

4.1. Non-Dissipative Heat Conduction

4.2. Fourier Model

4.3. Some Evolutionary Models

Maxwell–Cattaneo Equation

4.4. Cross-Couplings between Deformation and Heat Conduction

4.5. Second-Order Non-Local Models

5. Higher-Order Differential Equations

6. Rate Equations in the Eulerian Description

6.1. The Navier–Stokes–Voigt Fluid

6.2. Non-Heat Conducting NSV Fluid

6.3. NSV Fluid with Non-Fourier Conduction

7. Evolution Equations, Second Law, and Objective Derivatives

Remarks about Rate Equations and Thermodynamic Restrictions

8. Conclusions

Greenberg

1. Introduction

2. Notation and Balance Equations

2.1. Balance Equations

2.2. Second Law of Thermodynamics

2.3. Lagrangian Description

2.4. Representation Formula

3. Thermodynamic Restrictions

4. Remarkable Models of Heat Conduction

4.1. Non-Dissipative Heat Conduction

4.2. Fourier Model

4.3. Some Evolutionary Models

Maxwell–Cattaneo Equation

4.4. Cross-Couplings between Deformation and Heat Conduction

4.5. Second-Order Non-Local Models

5. Higher-Order Differential Equations

6. Rate Equations in the Eulerian Description

6.1. The Navier–Stokes–Voigt Fluid

6.2. Non-Heat Conducting NSV Fluid

6.3. NSV Fluid with Non-Fourier Conduction

7. Evolution Equations, Second Law, and Objective Derivatives

Remarks about Rate Equations and Thermodynamic Restrictions

8. Conclusions

Related Posts

Empowering Young Entrepreneurs to Drive Renewable Energy Solutions

U.S. gasoline prices this Thanksgiving are at their lowest since 2020

A Q&A about the ‘conspiracy’ of overconsumption

Greenberg