_{1}

^{*}

An extension of the linear irreversible thermodynamics is proposed through the inclusion of the first gradients of velocity and of the classical local state parameters as additional independent variables in the fundamental energy state equation of a fluid system. We show that consistency of this hypothesis with the energy balance equation leads to generalized nonlinear constitutive equations, which we discuss in terms of an isotropic non-Newtonian viscous fluid.

Non-Newtonian fluid mechanics is an area of active and intense research. Most of the materials that are of interest in a variety of manufacturing processes, exhibit non-Newtonian behaviors, which implies that the shear stress is not proportional to the shear rate [

Basic principles of continuum mechanics, namely, the balances of mass, momentum, energy, and entropy, lead to the following fundamental equations [

D ρ = 0 , (1)

D ( ρ v ) = ρ f + ∇ ⋅ T , (2)

D ( ρ ε ) = T T : G − ∇ ⋅ q , (3)

D ( ρ η ) = σ − ∇ ⋅ j , (4)

where, for any function F , operator D is defined by

D F ≡ d F / d t + F ∇ ⋅ v ≡ F ˙ + F ∇ ⋅ v . (5)

In total, Equations (1)-(4) constitute 6 independent equations for 22 unknown field variables, namely, mass density ρ, velocity v, specific internal energy ε, specific entropy η, stress tensor T, heat flux q, entropy flux j, and the entropy production σ (per unit volume and time). No distribution of heat sources is being considered, and the specific body force f is assumed to be known. Clearly, the foregoing basic equations are not adequate for the determination of these unknowns except for some trivial situations (for example, rigid body motions in the absence of heat conduction). Hence, 16 additional equations must be supplied in order to make the problem well-posed. Restriction to non-polar materials (i.e. when there are no assigned traction couples or body couples and no couple stresses), the moment of momentum principle establishes the symmetry of the stress tensor in general without any assumption of equilibrium or of uniformity of the stress distribution [1, pp. 215-216]. Then, in the case of non-polar materials, the number of unknown variables is reduced to 19.

In the derivation of the Equations (1) to (4), no specification has been required about any particular material. Therefore, the foregoing equations are not sufficient to explain fully the motions of materials having various types of physical properties. The character of the material is brought into the formulation through the constitutive equations, which specify the mechanical and thermal properties of particular materials based upon their internal constitution. These equations describe the relationships among the kinematic, mechanical, and thermal field variables, allowing the formulations of well-posed problems of continuum mechanics. Constitutive equations, from the physical standpoint, define various idealized materials that serve as models for the behavior of real materials.

In Gibbsian linear irreversible thermodynamics (LIT), the local equilibrium hypothesis (LEH) postulates the existence of a fundamental state equation for energy (or entropy), and the state equations of pressure and temperature are derived from it. However, the constitutive equations for the stress tensor and the heat flux remain as an open problem to be solved by the constitutive theory, which strives to find the form or at least restricts the generality of the constitutive functions [

For Newtonian fluids, the constitutive equation for the stress tensor is particularly simple: the deviatoric stress is linearly proportional to the rate of strain and the coefficient of proportionality is the viscosity. For non-Newtonian fluids, the constitutive equations can be much more complicated, implying, in general, a non-linear dependence of the stress tensor on the rate of strain. These constitutive relations must be built to reflect the macroscopic properties engendered by the fluid microstructure. One can go about this construction in several different ways.

In extended irreversible thermodynamics (EIT), for example, the LEH of LIT is replaced by an extension of the state variables space, accomplished by rising the thermodynamic fluxes (such as stress tensor and heat flux vector) to the status of additional independent variables. Then, the constitutive equations in this theory are intended to be recovered as solutions of the evolution equations of the new independent variables [

Our main objective in this paper is the development of a phenomenological gradient approach that consents to obtain generalized expressions for the constitutive equations of non-equilibrium fluid systems. To do it, we presumed that the local values of the ordinary state variables (mass density and specific entropy or energy) provide only a first approximation to the complete definition of the local state of a non-equilibrium fluid. In consequence, we also assumed that the specific internal energy of an isotropic, viscous and heat conducting fluid might depend not only on the ordinary state variables (mass density and specific entropy) but also on the first gradients of the mass density, specific entropy, and velocity fields. Our main hypothesis implies a fundamental state equation of energy that contains new information about the local behavior of the ordinary state variables and endows the fluid with a non-local character. Notwithstanding, the LEH is not relaxed completely in the theoretical framework of this approach, and the mass, momentum, energy and entropy balance equations, Equations (1)-(4), are still accepted as the complete set of fundamental physical laws that govern the fluid behavior.

We showed that consistence of our hypothesis with the energy balance equation brings out generalized non-linear expressions for the pressure and temperature state equations and for the constitutive equation of the stress tensor. We discussed this method for isotropic generalized viscous and heat conducting fluids. As simple examples, we applied this method to obtain the explicit forms of the constitutive equations for a second-order fluid [

The paper is organized as follows. In Section 2, the theoretical framework that supports the method we are proposing to derive constitutive equations for complex fluids is presented, and the method is developed. In particular, the generalized equations that consent the derivation of pressure, temperature, and stress tensor from the derivatives of the internal energy are drawn. In Section 3, we discussed the proposed method in terms of an isotropic fluid. Here, by considering the frame-indifferent and isotropic scalar character of the energy function, and with the base on the representation theorems for isotropic scalar functions of tensor arguments, the generalized expressions for pressure, temperature, and stress tensor are obtained. In Section 4, we applied the method to three cases of non-Newtonian viscous fluids for which an expression of the internal energy can be written explicitly as a function of the extended state variables such as the gradient of mass density and the shear rate, and the corresponding generalized constitutive equations were obtained. Finally, a brief paragraph of conclusions is presented in Section 5.

The LIT description of a system presupposes that one can find a partition of a non-uniform and non-equilibrium system with the property that each subsystem can be considered uniform and be locally described by exactly the same fundamental state equation of entropy that is valid under equilibrium conditions. This statement is the local equilibrium hypothesis and it implies that the Gibbs fundamental relation, which combines the first and second laws of thermodynamics, remains valid as a first approximation in the LIT description of a non-equilibrium system. The LEH restricts the applicability of LIT to systems close to equilibrium and presupposes that the gradients of the intensive properties are small enough that can be ignored locally. In case of fluid systems, LIT is only applicable to describe the Newtonian fluid behavior not far from equilibrium.

It is interesting to underline, however, that LIT definition of the entropy fundamental equation is the simplest among all possibilities because the second law requires merely that the entropy of an adiabatically isolated system should not decrease when a transition occurs between equilibrium states. Therefore, acceptable definitions may render an entropy function dependent on the values of the intensive variables at a point in the system, but also on their gradients, in such a way that one recovers the equilibrium entropy when the gradients go to zero [

Consequently, we postulate, within the framework of the energy representation, that the specific internal energy (ε) depends on the specific entropy (η) and mass density (ρ), but also on its first gradients ∇η, ∇ρ, and on the gradient of velocity, G (=∇v):

ε = ε ( ρ , ∇ ρ , η , ∇ η , G ) (6)

Under this hypothesis, we presume that a better description of the thermodynamic behavior of a non-equilibrium and non-uniform fluid system can be obtained.

The specific energy given by Equation (6), however, must be a frame-indifferent scalar function that reduces itself to the local equilibrium specific internal energy of the fluid, ε_{0}(ρ,η), when the gradients become negligible small:

ε ( ρ , ∇ ρ , η , ∇ η , G ) → ∇ ρ , ∇ η , G → 0 ε 0 ( ρ , η ) (7)

Then, our main hypothesis in this work comprises the assumptions expressed by the Equation (6) and Equation (7), and the acceptance of the Equations (1) to (4) as the complete set of fundamental physical laws that will govern the fluid behavior.

Equation (6) allows writing the following generalized Gibbs equation,

D ( ρ ε ) = θ D ( ρ η ) + ( p / ρ ) ρ ˙ − [ A + B + D ( ρ H ) ] : G + ∇ ⋅ h + D ( ρ H : G ) (8)

where

ρ θ = ρ ( ∂ ε / ∂ η ) − ∇ ⋅ [ ρ ( ∂ ε / ∂ ∇ η ) ] (9)

p / ρ = ρ ( ∂ ε / ∂ ρ ) − ∇ ⋅ [ ρ ( ∂ ε / ∂ ∇ ρ ) ] (10)

A = ρ [ ( ∂ ε / ∂ ∇ ρ ) ∇ ρ + ( ∂ ε / ∂ ∇ η ) ∇ η ] (11)

B = ρ ( G T ⋅ H − H ⋅ G T ) (12)

H = ∂ ε / ∂ G (13)

h = ρ [ ( ∂ ε / ∂ ∇ ρ ) ρ ˙ + ( ∂ ε / ∂ ∇ η ) η ˙ ] (14)

Equation (9) and Equation (10) are generalized state equations for temperature θ and pressure p, respectively. When the specific energy does not depend on the gradients of specific entropy and mass density, these equations give the local equilibrium temperature and pressure.

Now, with the help of the mass and entropy balance equations (Equation (1) and Equation (4), respectively), Equation (8) can be written as

D ( ρ ε ) = − [ p I + A + B + D ( ρ H ) ] : G − ∇ ⋅ ( θ j − h ) + [ θ σ + j ⋅ ∇ θ + D ( ρ H : G ) ] (15)

where I denotes the unit tensor.

Then, from the comparison of Equation (15) with Equation (3), we obtain the following expressions for the stress tensor T, the entropy flux j, and the entropy production σ:

T = − [ p I + A + B + D ( ρ H ) ] T (16)

j = ( q + h ) / θ (17)

σ = − [ j ⋅ ∇ θ + D ( ρ H : G ) ] / θ (18)

In principle, the Equations (9), (10) and (16) allow to calculate the temperature (θ), pressure (p), and the stress tensor (T) of the fluid from the derivatives of the generalized specific internal energy with respect the classical state variables, its gradients, and the gradient of velocity. Equation (17) and Equation (18), however, provide expressions for the entropy flux (j) and entropy production (σ) which also involve the derivatives of the generalized internal energy, but these properties are not completely determined by these energy derivatives because the heat flux (q) remains as an unknown within this formalism.

Within the framework of a phenomenological theory, the internal energy function must be determined by experiment for each particular fluid. However, it is possible to obtain important information about the constitutive equations for an isotropic fluid if we consider the frame-indifferent and isotropic scalar character of the energy function^{1}.

Therefore, in the case of an isotropic fluid, the specific internal energy can be a function of the velocity gradient, G, only through the strain tensor

D = 1 2 ( G + G T ) , (19)

because G itself is non-indifferent [

ε ( ρ , ∇ ρ , η , ∇ η , G ) = ε ( ρ , ∇ ρ , η , ∇ η , D ) . (20)

Now, according to the representation theorems for isotropic scalar functions of tensor arguments, the specific internal energy of an isotropic fluid can be expressed as a function of the scalar invariants I_{α} of ∇ρ, ∇η and D, that is,

ε ( ρ , ∇ ρ , η , ∇ η , D ) = ε ( ρ , η , I α ) . (21)

Therefore, the partial derivatives of the energy function with respect to the variables ∇ρ, ∇η and D can be expressed in the form

( ∂ ε / ∂ ∇ u ) = ∑ A α ( ∂ I α / ∂ ∇ u ) ≡ ∑ ( ∂ ε / ∂ I α ) ( ∂ I α / ∂ ∇ u ) (22)

where the sums run over the number of the independent scalar invariants. In general, the coefficients A_{α} depend on the same variables as the internal energy.

Hereafter, we will be considering a specific internal energy dependent only on the first nine independent scalar invariants of ∇ρ, ∇η and D. These scalar invariants involve the variables ∇ρ, ∇η and D up to third order:

I 1 = ∇ η ⋅ ∇ η , I 2 = ∇ η ⋅ ∇ ρ , I 3 = ∇ ρ ⋅ ∇ ρ , I 4 = T r ( D ) , I 5 = T r ( D 2 ) , I 6 = T r ( D 3 ) , I 7 = ∇ η ⋅ D ⋅ ∇ η , I 8 = ∇ η ⋅ D ⋅ ∇ ρ , I 9 = ∇ ρ ⋅ D ⋅ ∇ ρ . (23)

^{1}A scalar function F of a tensor argument X is said to be indifferent if it remains unchanged after a change of frame, so that F ( X * ) = F ( X ) , where X^{*} is the image of X under the change of frame. In addition, the function F is an isotropic scalar function of the tensor argument X, if it satisfies X * = Q ⋅ X ⋅ Q T , for any orthogonal tensor Q [

In this case, Equation (16) leads to the following expression for the generalized stress tensor:

T = − ( p + π ) I + 2 μ { D + α 1 [ D ˙ + ( D ⋅ Ω − Ω ⋅ D ) ] + α 2 D 2 + α 3 ( D ⋅ D ˙ + D ˙ ⋅ D ) } − 2 ρ ( A 1 ∇ η ∇ η + A 2 ∇ η ∇ ρ ¯ ¯ + A 3 ∇ ρ ∇ ρ ) − D [ ρ ( A 7 ∇ η ∇ η + A 8 ∇ η ∇ ρ ¯ ¯ + A 9 ∇ ρ ∇ ρ ) ] (24)

where only the terms up to second order in ∇ρ, ∇η and D have been retained. In this equation, we have introduced the vorticity tensor

Ω = 1 2 ( G − G T ) , (25)

and the following definitions

π = ρ A ˙ 4 , μ = − ρ A ˙ 5 , ∇ η ∇ ρ ¯ ¯ = 1 2 ( ∇ η ∇ ρ + ∇ ρ ∇ η ) , (26)

α 1 = A 5 / A ˙ 5 , α 2 = 3 2 ( A ˙ 6 / A ˙ 5 ) , α 3 = 3 2 ( A 6 / A ˙ 5 ) . (27)

The generalized state equations for temperature and pressure are obtained in a similar way starting from Equation (9) and Equation (10), respectively. The results are

θ = ( ∂ ε / ∂ η ) − ( B 1 | ∇ η | 2 + B 2 ∇ η ⋅ ∇ ρ + B 3 | ∇ ρ | 2 ) − [ 2 ( A 1 I + A 7 D ) : N + ( A 2 I + A 8 D ) : M ] − [ 2 A 7 ( ∇ ⋅ D ) ⋅ ∇ η + A 8 ( ∇ ⋅ D ) ⋅ ∇ ρ ] (28)

p = ρ 2 ( ∂ ε / ∂ ρ ) − ( C 1 | ∇ η | 2 + C 2 ∇ η ⋅ ∇ ρ + C 3 | ∇ ρ | 2 ) − ρ 2 [ ( A 2 I + A 8 D ) : N + 2 ( A 3 I + A 9 D ) : M ] − ρ 2 [ A 8 ( ∇ ⋅ D ) ⋅ ∇ η + 2 A 9 ( ∇ ⋅ D ) ⋅ ∇ ρ ] (29)

where

B 1 = 2 ( ∂ A 1 / ∂ η ) , B 2 = ( 1 / ρ ) [ 2 A 1 + 2 ( ∂ A 1 / ∂ ρ ) + ( ∂ A 2 / ∂ η ) ] , B 3 = ( 1 / ρ ) [ A 2 + ρ ( ∂ A 2 / ∂ ρ ) ] , (30)

C 1 = ρ 2 ( ∂ A 2 / ∂ η ) , C 2 = ρ [ A 2 + ρ ( ∂ A 2 / ∂ ρ ) + 2 ρ ( ∂ A 3 / ∂ η ) ] , C 3 = 2 ρ [ A 3 + ( ∂ A 3 / ∂ ρ ) ] , (31)

M = ∇ ( ∇ ρ ) , N = ∇ ( ∇ η ) . (32)

Equation (24), Equation (28) and Equation (29) reduce themselves to the corresponding linear expressions of LIT when the time rate terms are negligible and the gradients are small enough to neglect the quadratic and higher order terms involving them. When this is not the case, the additional terms in these equations represent nonlinear corrections that become relevant when gradients are not locally negligible. In particular, the expression for the stress tensor contains some additional terms that one can find in rheological constitutive equations [

T = − p I + 2 μ 0 D + 4 μ 1 D 2 + 2 μ 2 [ D ˙ + ( D ⋅ Ω − Ω ⋅ D ) ] . (33)

In this section, as simple applications of the proposed method, we used the generalized equations for pressure, temperature, and the stress tensor that we found in the previous section, in order to obtain explicit forms of these properties for three cases already reported in literature. In the first case, we started with the square gradient approximation for the specific internal energy, and the corresponding constitutive equations were obtained. In the second case, we proposed an expression for the internal energy of a power law viscous fluid, and then the known constitutive equations were derived. Finally, we carry out a comparison with results reported from non-equilibrium molecular dynamics simulations.

From the standpoint of the square gradient approximation [

ε ( η , ρ , ∇ ρ ) = ε 0 ( η , ρ ) + ( b / 2 ρ ) | ∇ ρ | 2 , (34)

where ε_{0} is the local equilibrium energy and b is a constant. Application of Equation (24) to the energy given by Equation (34), leads to the following stress tensor

T = − [ p 0 − b ( 1 2 | ∇ ρ | 2 + ρ ∇ 2 ρ ) ] I − b ∇ ρ ∇ ρ , (35)

with p_{0} being the local equilibrium pressure given by

p 0 = ρ 2 ( ∂ ε 0 / ∂ ρ ) . (36)

The pressure p is defined by

p = − 1 3 T r ( T ) = p 0 − 1 6 b | ∇ ρ | 2 − b ρ ∇ 2 ρ . (37)

Equation (35) is the Lovett expression [

A simple generalization of the Newtonian fluid is a phenomenological model that assumes that the flow regime changes the fluid viscosity (the rate of dissipation), but the tensorial structure of the Newtonian constitutive equation remains unchanged. The constitutive equation for this class of fluid models has the following form

T = − p I + 2 μ ( D ) D , (37)

where μ(D) is the viscosity, which is different from the Newtonian viscosity due to its possible dependence upon the fluid flow [

γ = D : D = T r ( D 2 ) ,

(i.e., on the square root of the second scalar invariant of the tensor D). Then

μ ( D ) = μ ( γ ) .

The material properties of a generalized Newtonian fluid are determined entirely by the behavior of the function μ(γ). Often, over an interval of shear rate, a straight line can approximate the relationship between the apparent viscosity and the shear rate, when plotted in log-log coordinates, i.e.,

μ ( γ ) = K γ n − 1 . (38)

This expression for the viscosity is called the Ostwald de Weale model or the power law model [

T ˜ = 2 K γ n − 1 D . (39)

Now, for the power law fluid, let us assume a specific internal energy with the form

ε ( ρ , η , γ ) = ε 0 ( ρ , η ) + ε 1 ( ρ , η ) γ n + 1 , (40)

where ε_{0} denote the local equilibrium value of energy, ε_{1} denote a non-equilibrium coefficient, and n is the power law index.

This expression for energy constitutes a particular case of the more general expression given by the Equation (21) because it depends on the gradients of the velocity field through the shear rate. Then, for the power law fluid with an energy given by Equation (40), the Equation (28), Equation (29) and Equation (24) can be used to obtain the temperature and pressure generalized state equations and the stress tensor. Respectively, the results were

θ ( ρ , η , γ ) = θ 0 ( ρ , η ) + θ 1 ( ρ , η ) γ n + 1 , (41)

p ( ρ , η , γ ) = p 0 ( ρ , η ) + p 1 ( ρ , η ) γ n + 1 , (42)

T = − p ( ρ , η , γ ) I + 2 [ β ( ρ , η ) γ n − 1 ] D . (43)

Here, the terms involving time derivatives of shear rate were ignored, and the following symbols were introduced

θ i = ( ∂ ε i / ∂ η ) , i = 0 , 1 , (44)

p i = ρ 2 ( ∂ ε i / ∂ ρ ) , i = 0 , 1 , (45)

β ( ρ , η ) = − 1 2 ( n + 1 ) ρ ε ˙ 1 . (46)

Equation (43) gives the deviatoric stress tensor with an apparent viscosity with the same dependence on shear rate as in Equation (38) of the power law model by Ostwald de Weale [

The properties of a fluid under a shear have been investigated with non-equilibrium molecular dynamics (NEMD) simulations [

μ = μ 0 − μ 1 γ 1 2 , (47)

p = p 0 + p 1 γ 3 2 , (48)

ε = ε 0 + ε 1 γ 3 2 , (49)

where the quantities with zero subscript are the zero strain rate quantities, so that, in particular, μ_{0} is the Newtonian viscosity. At temperatures at or near the triple point of a Lennard-Jones Fluid, the Equation (47) for viscosity is in good agreement with the simulation results. However, better overall agreement for other temperatures and densities has been obtained using other values of the exponent of the shear rate [

For n = 1⁄2, Equation (40) and Equation (42) give the specific energy and the pressure with the same dependence on shear rate γ as in Equation (49) and Equation (48). For n = 1⁄2, however, Equation (43) gives a deviatoric stress tensor with an apparent viscosity with dependence on shear rate to the power minus one-half (γ^{−1⁄2}) which is different from that predicted by mode-coupling theory in the NEMD simulations.

It is worth to mention, however, that the power law model, Equation (39), is perhaps the most widely used model in literature dealing with process engineering applications, and that there exist a variety of substances whose viscosity presents a dependence on γ^{−1⁄2} (power law index n = 1/2), such as chocolate, guava puree, papaya puree, and tomato paste, among others [

A slight generalization of the local equilibrium energy equation carried out with the inclusion of the first gradients of mass density, specific entropy, and velocity as additional state variables allowed developing a method to derive non-linear expressions for pressure, temperature, and the stress tensor for complex fluids. The proposed method was discussed in terms of an isotropic viscous fluid. As simple examples, it was shown that the generalized stress tensor could be reduced, under specific constraints, to the stress tensor of second-order fluids, to the Lovett expression for the stress tensor of a non-uniform fluid, and to stress tensor of a power law fluid. We also compared the results of the proposed method with results of NEMD simulations reported in the literature, and we could reproduce the pressure expression obtained for a Lennard-Jones fluid under shear.

Salcido, A. (2018) A Phenomenological Gradient Approach to Generalized Constitutive Equations for Isotropic Fluids. Journal of Applied Mathematics and Physics, 6, 1494-1506. https://doi.org/10.4236/jamp.2018.67126