Module-5: Constitutive Relations: Lecture-38: Thermomechanics of Fluids
Module-5: Constitutive Relations: Lecture-38: Thermomechanics of Fluids
Module-5: Constitutive Relations: Lecture-38: Thermomechanics of Fluids
τ = τ̂ (X, F , θ, g, Ḟ ),
q = q̂(X, F , θ, g, Ḟ ),
e = ê(X, F , θ, g, Ḟ ),
η = η̂(X, F , θ, g, Ḟ ),
for every unimodular tensor H, i.e., for every H such that det(H) = 1. Let L and ρ be
velocity gradient and density field, respectively. Then, choosing H = J 1/3 F −1 , we get
since ρ0 = ρJ.
Proof. The given constitutive relation τ = τ̃ (X, ρ, θ, g, L). Let F ∗ and F be fixed and
moving frames as discussed in the principle of material frame-indifference (see Lecture-34).
Then
τ̃ (X ∗ , ρ∗ , θ∗ , g ∗ , L∗ ) = Qτ̃ (X, ρ, θ, g, L)QT ,
where Q(t) is relative rotation between frames. Since X ∗ = X, ρ∗ = ρ, θ∗ = θ, g ∗ = g and
L∗ = QLQT + Q̇QT (see Lectures 31, 32 and 35) and the velocity gradient L = D + W
where D is symmetric part and W is skew symmetric part, above equation can be written
as
τ̃ (X, ρ, θ, g, QDQT + QW QT + Q̇QT ) = Qτ̃ (X, ρ, θ, g, L)QT , (2)
where Q is arbitrary rotation tensor. For any given skewsymmetric tensor W , we can
choose an orthogonal tensor such that
Q(0) = I and Q̇(0) = −W .
Substituting these relations in Eq. (2), we get
τ̃ (X, ρ, θ, g, L) = τ̃ (X, ρ, θ, g, D).
Thus, the principle of material frame-indifference implies above relation. Since the given
function is isotropic, the converse is trivial.
Applying similar steps that are discussed in Problem-1 to other constitutive relations,
we can show that they are also function of rate of deformation. In summary, we have
τ = τ̃ (X, ρ, θ, g, D)
q = q̃(X, ρ, θ, g, D)
e = ẽ(X, ρ, θ, g, D)
η = η̃(X, ρ, θ, g, D), (3)
The following Clausius-Duhem inequality (second law of thermodynamics) in terms of
internal energy e and entropy η can be obtained from Eqs. (10) and (13) of Lecture-30.
De Dη 1 q·g
−θ − τ :D+ ≤ 0. (4)
Dt Dt ρ ρθ
Since the constitutive relations are not direct function of time, we have
De ∂e ∂e ∂e ∂e Dη ∂η ∂η ∂η ∂η
= ρ̇ + θ̇ + · ġ + : Ḋ and = ρ̇ + θ̇ + · ġ + : Ḋ.
Dt ∂ρ ∂θ ∂g ∂D Dt ∂ρ ∂θ ∂g ∂D
Substituting these relations in above Clausius-Duhem inequality, we get
! ! !
∂e ∂η ∂e ∂η ∂e ∂η
−θ ρ̇ + −θ θ̇ + −θ · ġ
∂ρ ∂ρ ∂θ ∂θ ∂g ∂g
!
∂e ∂η 1 q·g
+ −θ : Ḋ − τ :D+ ≤ 0.
∂D ∂D ρ ρθ
Let v be a spatial velocity field. Then the conservation of mass can be written as ρ̇ =
−ρ ∇x · v = −ρ tr(L) = −ρ I : L. Substitution of this relation yields
! ! !
∂e ∂η ∂e ∂η ∂e ∂η
−θ θ̇ + −θ · ġ + −θ : Ḋ
∂θ ∂θ ∂g ∂g ∂D ∂D
! !
∂e ∂η 1 q·g
+ − ρ − ρθ I− τ :D+ ≤ 0. (5)
∂ρ ∂ρ ρ ρθ
Clearly, the constitutive relations are not function of θ̇, ġ and Ḋ as shown in Eq. (3).
Since θ̇, ġ and Ḋ are arbitrary, we get
∂e ∂η ∂e ∂η ∂e ∂η
−θ = 0, −θ = 0, and −θ =O (6)
∂θ ∂θ ∂g ∂g ∂D ∂D
! !
∂e ∂η 1 q·g
ρ − ρθ I+ τ :D− ≥0 (7)
∂ρ ∂ρ ρ ρθ
Differentiating Eq. (6)1 with g and Eq. (6)2 with θ, we can write
∂ 2e ∂ 2η ∂ 2e ∂ 2η ∂η
−θ = 0 and −θ − = 0.
∂g∂θ ∂g∂θ ∂θ∂g ∂θ∂g ∂g
Combining these two equations yield
∂η
= 0. (8)
∂g
∂e
= 0. (9)
∂g
Thus, the internal energy and entropy are independent of temperature gradient g. Similar
steps applied to Eq. (6)1 and Eq. (6)3 , we get
∂η ∂e
= O and =O (10)
∂D ∂D
In conclusion, constitutive relation to internal energy and entropy reduces to
δ(ρ, θ, g, D) ≥ 0 (12)
The thermodynamic equilibrium is defined as the process that is independent of time and
also uniform. Physically, the thermodynamic equilibrium refers to mechanical equilibrium
(force balance), thermal equilibrium (no heat transfer between system and surroundings)
and chemical equilibrium (no change in reactants and products with time). For the
thermodynamic equilibrium, we have the following velocity and temperature fields.
D = O and g = 0. (15)
δ(ρ, θ, 0, O) = 0. (16)
It is clear from Eqs. (12) and (16) that the dissipation function δ is minimum at thermody-
namic equilibrium. In order to have minimum dissipation at thermodynamic equilibrium,
we need to have the necessary condition
∂δ ∂δ
= O and =0 (17)
∂D ∂g
and positive definite Hessian (second derivative of dissipation) for sufficiency condition.
The Hessian of dissipation function is given by
∂ 2δ ∂ 2δ
∂D∂D ∂D∂g
H=
. (18)
∂ 2δ ∂ 2δ
∂D∂g ∂g∂g
Taking derivatives to the dissipation function (see Eq. (13)), we get
!
∂δ ∂e ∂η ∂τ 1 ∂q
= τ + ρ2 −θ I+ :D− · g, (19)
∂D ∂ρ ∂ρ ∂D θ ∂D
∂δ ∂τ q 1 ∂q
= :D− − · g, (20)
∂g ∂g θ θ ∂g
where
! !
∂τ ∂τkl ∂q ∂qk
:D = Dkl , ·g = gk ,
∂D ij
∂Dij ∂D ij
∂Dij
! !
∂τ ∂τjk ∂q ∂qj
:D = Djk , and ·g = gj .
∂g i
∂gi ∂g i
∂gi
In above relations the subscript T .E. represents thermodynamic equilibrium. Since the
derivatives should vanish at thermodynamic equilibrium (see Eq. (17)), we have
τ = −pI and q = 0, (22)
T .E. T .E.
Since the free energy ψ = e − ηθ (see Lecture-30), the thermodynamic pressure can be
written as
∂ψ
p = ρ2 . (24)
∂ρ
In other words, the stress tensor can be decomposed into two parts: (i) thermodynamic
equilibrium stress (ii) viscous stress. Therefore, the stress tensor
τ = −pI + σ. (27)
• The constitutive relations for internal energy and entropy are given by
where !
2 ∂e ∂η
p=ρ −θ .
∂ρ ∂ρ
• Let σ = σ̃(X, ρ, θ, g, D) be viscous stress then the total stress is given by
τ = −pI + σ.
σ = σ̃(ρ, θ, g, D),
q = q̃(ρ, θ, g, D),
e = ẽ(ρ, θ),
η = η̃(ρ, θ).
The isotropic functions of viscous stress and heat flux have the following representation1
It is clear that both the viscous stress σ and also the heat flux q are dependent on rate
of deformation D and temperature gradient g.
Uncoupled viscous heat-conducting fluids:
In case of uncoupled system, we assume that viscous stress is independent of temperature
gradient and the heat flux is independent of rate of deformation. The constitutive relations
for uncoupled system can be written as
where ID represents set of three invariants of D. We can have two special cases: (i) rate
of deformation is zero (D = O) and temperature gradient is non-zero (g 6= 0) (ii) rate
of deformation is non-zero (D 6= O) and temperature gradient is zero (g = 0). We note
that later case is also known as Reiner-Rivlin fluid.
Heat conducting fluid:
Since the rate of deformation is zero (D = O), we only have heat conduction inequality
q · g ≤ 0 =⇒ β0 (ρ, θ, g · g, O) g · g ≤ 0. (32)
k = −β0 . (33)
1
Liu, I-S., Continuum Mechanics, Springer-Verlag, Berlin, 2002.
q = −kg. (34)
This relation is known as Fourier law of heat conduction for fluids. Furthermore, the
Fourier law is similar to both solids and also for fluids. Through this law is important for
fluids, there are many instances where convective heat transfer is more than conduction.
Thus, in many practical problems involving fluids the conduction can be ignored. There
are few problems where the conduction is important. For example if the fluid does not have
sufficient space (thin fluid layer between walls) to convect then conduction is important.
Reiner-Rivlin fluids:
In this case, we assume that the temperature gradient is zero (g = 0). The constitutive
relation for the viscous stress can be written as
σ : D ≥ 0. (36)
References