WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
Stanislaw Sieniutycz, Piotr Kuran
Variational Formulations and Extended Thermodynamics for Resting Irreversible Fluids with Heat Flow STANISLAW SIENIUTYCZ and PIOTR KURAN Faculty of Chemical and Process Engineering Warsaw University of Technology PL 00-645, 1 Waryńskiego Street, Warsaw POLAND
[email protected] [email protected] http://www. ichip.pw.edu.pl Abstract: - Nonequilibrium statistical mechanics helps to estimate corrections to the entropy and energy of the fluid with heat flux in terms of the nonequilibrium distribution function, f. This leads to the coefficients of wave model of heat: relaxation time, propagation speed and thermal inertia. With these data a quadratic Lagrangian and a variational principle of Hamilton’s type follows for the fluid in the field representation of fluid’s motion. We analyze canonical conservation laws and show the satisfaction of the second law under the constraint of these conservation laws. Key-Words: - Grad solution, variational calculus, wave equations, conservation laws, entropy. rather than mathematical criteria. An analysis shows that the approach is particularly useful in the field (Eulerian) description of transport phenomena, where equations of thermal field follow from variational principles based on the state adjoints rather than on the original physical variables. Exemplifying process is hyperbolic heat transfer, but the approach can also be applied to coupled parabolic transfer of heat, mass and electric charge. With various gradient or non-gradient representations of physical fields in terms of state adjoints useful actiontype criteria emerge. Symmetry principles are effective, and components of the formal energy-momentum tensor can be found. Focusing on heat flow, our work represents, in fact, an approach that shows the advantage of approaches borrowed from the optimal control theory in the variational setting of irreversible transport. The limiting reversible process provides a suitable reference frame for more involved irreversible evolutions.
1 Introduction Extended thermodynamics of fluids can be applied to set variational principles for the irreversible energy transfer. Significant help can be obtained from nonequilibrium statistical theories when evaluating kinetic or fluxdependent terms in energies and macroscopic Lagrangians. Especially, we can treat statistical aspects of nonequilibrium fluids with heat flow by applying an analysis that uses Grad’s results [1] to determine nonequilibrium corrections ∆s or ∆e to the energy e or entropy s in terms of the nonequilibrium density distribution function f. To find the required corrections to the energy e and kinetic potential L we exploit corrections ∆s and a relationship that links energy and entropy representations of thermodynamics. We may also evaluate coefficients of wave model of heat, such as: relaxation time, propagation speed and thermal inertia factors, g and θ. With these data we can formulate a variational principle of Hamilton’s or least action type for fluids with heat flux in the field or Eulerian representation of fluid motion. To find a variational formulation we apply here an approach that adjoints a given set of constraints to a kinetic potential L and transfers the original variational formulation to the space of the Lagrange multipliers (also called state adjoints). Considering limiting reversible process we evaluate canonical components of energy-momentum tensor and associated conservation laws. The approach works efficiently; it leads to exact imbedding of constraints in the potential space of Lagrange multipliers, implying that the appropriateness of the constraining set should be verified by physical
ISSN: 1991-8747
2 Optimization Type Approach Statistical theories are useful [1] to evaluate nonequilibrium corrections to the energy and other thermodynamic potentials in situations when a continuum is inhomogeneous because of the presence of irreversible fluxes. To illustrate benefits resulting from nonequilibrium statistical thermodynamics, heat transfer in locally non-equilibrium fluids is analyzed [2]. Quite essential therein is the connection between various representations of thermodynamics and a relationship (resembling the Gouy-Stodola law) that links energy and entropy pictures. With this relationship nonequilibrium corrections to the energy can be found from those known
759
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
For a given nonequilibrium state at a point A in Fig.1, two equilibrium reference states, at points B and C, say, correspond, respectively, to the energy and entropy representation. A researcher knowing entropy s (e. g., from distribution function f corresponding to A) formulates his description of state A in terms of equilibrium parameters at B, for a set of variables, here the entropy flux js. Yet, one who knows energy e can base his view on the heat flux q and equilibrium at C. When point A moves the equilibrium states (B and C) vary. The conventional picture of motion in terms of Hamilton's principle corresponds to following the behavior of B and the kinetic energy of entropy flux, whereas the kinetic theory view corresponds to tracking of C and the deviation of entropy from equilibrium. The transition from one view to the other is possible [2]. It is important to realize that for a single nonequilibrium state the use of the entropy representation and energy representation establishes two different equilibrium states located on the Gibbs surface. This, of course, is because of the difference in what is held constant. The distance between two discussed equilibrium states (B and C) increases with the distance of the state A from the Gibbs surface. This distance can also be measured in terms of the modulus of the flux js or in terms of the differences ∆e = AB or ∆s = AC. When the curvature of the Gibb's surface can be neglected, corresponding to the near-equilibrium situation, the two disequilibrium excesses are linked by an equality resembling the Gouy-Stodola law
for the entropy. These energy corrections can next be used to find a suitable kinetic potentials L. The present approach is optimization-type; it differs from the conventional variational ones in that the action functional is systematically constructed rather than assumed from the beginning. Equations of constraints (reversible or not) follow in the form of their counterparts in the space of Lagrange multipliers; they are extremum conditions for the action based on a composite (constraint involving) Lagrangian Λ or its gauge counterpart. As long as the representations of physical variables are explicit in terms of Lagrange multipliers, the whole variational formalism can be transferred to the adjoint space of multipliers, i.e. a variational principle can be formulated in the adjoint space. The Lagrangian can also be used to obtain the matter tensor and associated conservation laws. Finally we show that the use of the canonical conservation laws constructed for the reversible process and variational extremum conditions assures the satisfaction of the second law of thermodynamics, the property that renders the variational theory considered a candidate to be the physical one.
3 Energy and Entropy Representations Now our task is to recall some basic knowledge on the thermodynamics of heat flow without local equilibrium. This will help us to construct Lagrangians, variational principles and conservation laws. We work in the framework of extended thermodynamics of fluids [3] and restrict to incompressible, one-component fluid with heat flow. Consider a nonequilibrium state, say A, off but near the Gibbs surface, when the local equilibrium assumption is inapplicable, Fig. 1. The energy at the state A is the nonequilibrium internal energy. This internal energy depends not only on the usual state variables, but also on some nonequilibrium variables such as heat flux or diffusive entropy flux. Here we select the heat flux, q, as the nonequilibrium variable of choice. It is treated as an internal variable which may relax to an equilibrium. Nonequilibrium energy density of a continuum, ρe, or its specific energy e, is a function of density ρ, specific entropy s and diffusive entropy flux js or heat flux q. For a continuum as a stable macrosystem, its equilibrium internal energy density ρeeq is the minimum of ρe with respect to unconstrained js or q, at constant ρ and s. As ρ = v-1, the reciprocal of specific volume, the minimum of ρe (or e itself) with respect to js or q occurs at the constant entropy s and volume v which are proper variables at which the energy attains minimum at equilibrium. This is in agreement with basic thermodynamics [4]. Since js or q are a diffusive fluxes, the minimum occurs for js = 0 or q=0.
ISSN: 1991-8747
Stanislaw Sieniutycz, Piotr Kuran
(e − eeq ) s, ρ = −T ( s − seq )e, ρ
(1)
This formula states that the energy released during the isoentropic relaxation equals the product of the absolute temperature and the entropy deficiency in the system caused by the presence of the heat flux q or the entropy flux, js.
760
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
Fig. 1. Diverse reference equilibrium states (B, D, C, etc.) for a given nonequilibrium state A.
would correspond to the minimum energy. The second order correction to the entropy density (in entropy representation) is 1 ρ s( 2 ) = ρ∆s = − k B f eqϕ12 dc (7) 2 Hence, in view of the relation between ∆e and ∆s implied by Fig. 1 or Eq. (1) 1 ∆e = − k B Tρ −1 f eqϕ12 dc (8) 2 Since the state A is close to the equilibrium surface, the multiplicative coefficients involving usual thermostatic variables can always be evaluated at arbitrary equilibrium points (B or C in Fig. 1). However, in Eqs. (1), (9) and (10), they are evaluated (in the kinetic theory) for the case of the isoenergetic equilibrium (point C, Fig. 1). The function ϕ1, obtained in Grad's method when the system's disequilibrium is maintained by a heat flux q is 2 1 5 φ1 = (m / Pk B2 T 2 )( mC 2 − k B T )C.q (9) 5 2 2 where m is the mass of a molecule ([1], [3]). From Eqs. (7), (8) and (9) one obtains for the entropy deviation 1 ∆s = − (m / ρPk B T 2 )q 2 (10) 5 and for the energy deviation, Eq. (1), with entropy flux js = qT-1 1 1 ∆e = (m 2 / k B ρ 2 ) j 2s = ρ − 2 gj 2s (11) 5 2 Equations (10) and (11) hold to the accuracy of the thirteenth moment of the velocity [1]. When passing from Eq. (10) to (11) state equation P = ρkBTm-1 is used and a constant g is defined as 2mTρ 2m 2 g≡ = 2. (12)
4 Nonequilibrium Corrections to Energy and Entropy
∫
It is essential that the entropy representation is assumed in the Grad’s formalism of the kinetic theory [1]. Hence the specific energy of an ideal gas or fluid with heat at the point A is equal to the specific energy at equilibrium C in Fig. 1. The reference temperatures and pressures that appear in the expressions of kinetic theory are T(C) and P(C). From the formalism one finds disequilibrium corrections ∆s or ∆e in terms of the non-equilibrium density distribution function f. Here we recapitulate the results of several different works [3]-[6] all using Grad’s [1] solution of the Boltzmann equation in macroscopic predictions for dilute gas of rigid spheres. The molecular velocity distribution function, f, out of equilibrium but close to it is given as f (C) = f eq (C)(1 + ϕ1 )
∫
(2)
where feq is the local equilibrium (Maxwell-Boltzmann) distribution pertaining to the entropy representation equilibrium (point C, Fig. 1). f and feq are scalars, but functions of the peculiar velocity C = c - u, and ϕ1 is a function of the deviation from equilibrium. This deviation is expressed in terms of the gradT in the Chapman-Enskog method and in terms of the heat flux q in the Grad's method. Using Eq. (2) in the entropy definition, one integrates the expression flnf over all of the space of the molecular velocity c, ρ s = −k B
∫ f ln fdc
(3)
5Pk B
Proceeding with development of ρs up to second order in φ1, one obtains ρ s = ρ seq + ρ s(1) + ρ s( 2) .... with the local equilibrium entropy ρ seq = −k B
∫f
eq
ln f
eq
dc
(4)
∫f
eq
ϕ1 ln f eq dc = 0
(5)
Again, this proves that one deals with the entropy representation where the entropy is maximum at equilibrium. A counterpart of the above equation in the energy representation
∫
ρ e(1) = − f eq ϕ1 mc 2 dc = 0,
ISSN: 1991-8747
5k B
Here we abandoned the entropy representation. Pressure in Eqs. (9) and (12) is that of an ideal gas, given by the definition used in the kinetic theory (Grad 1958 [1]). Eq. (11) with constant g defined by Eq. (12) is the characteristic feature of the ideal monoatomic gas (dilute Boltzmann gas composed of hard spheres). For arbitrary fluids (polyatomic gases, dense monoatomic gases and liquids) one can retain the form of the last expression in Eq. (11) by using a general formula for by noting that 2 g ( ρ, s ) ≡ ρ eq (∂ 2 e / ∂j 2s ) eq (13)
and the nonequilibrium correction ρ s(1) = −k B
Stanislaw Sieniutycz, Piotr Kuran
In the ideal gas case the derivative ∂2e/ ∂j2s = (2/5)(m2/kB2ρ2) from Eq. (11) and the definition (12) is recovered form definition (13). quation (13) is consistent with a hypothesis about the equality of the kinetic and static nonequilibrium energy corrections in a thermal shock-wave front [5]. The hypothesis can be used to compute (∂2e/∂js2)eq for arbitrary fluids as T/(ρcpG) and
(6)
761
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
hence g as Tρ/(cpG), where G is the shear modulus. Equilibrium values of thermodynamic parameters can be applied in such expressions. For an ideal gas the shear modulus is just the pressure P (the result of Maxwell) and cp = 5kB/(2m). These results allow one to recover definition (12) from the expression g = Tρ/(cpG); they support the hypothesis mentioned above. Yet, for the purpose of general considerations the use of the implicit dependence of g on the basic variables (ρ, s) is often enough, i. e., function g(ρ, s) will be used when passing to arbitrary fluids. Some entropy flux adjoints, as and is, are useful. They are defined, respectively, by equations a s ( s, ρ, j s ) = ∂∆e( s, ρ , j s ) ρ , s / ∂j s = gρ −2 j s
(14)
i s ( s, ρ, j s ) = gρ −1 j s = gsv s = gs (u s − u)
(15)
1/ 2
a
1/ 2
τ
.
(17)
Substituting to this formula the ideal gas data, i.e. g of Eq. (12) and cp = 5kB/(2m), yields propagation speed in the ideal gas 1/ 2
T c0 = cp g
1/ 2
k T = B m
(18)
(thermal speed). Thus the results of nonequilibrium statistical mechanics help to estimate coefficients of the heat transfer model. Data ofτ and c0 are used below in a variational principle for heat transfer. One more coefficient that is quite useful in the wave theory of heat is that describing a thermal mass per unit of entropy θ=T c0−2 ,
The entropy diffusion velocity vs = us - u = js/ρs appears in Eq. (15). One may also introduce the product kBgs which has the dimension of mass. For the ideal gas this product yields ms = 2/5(m2skB-1), a measure of heat inertia. In the model of a constant g, nonequilibrium temperature T(B) is equal to the equilibrium temperature T(ρ, s) which is both the measure of mean kinetic energy of an equilibrium and the derivative of energy with respect to the entropy. This equality emerges because, above, we have chosen the entropy flux js, not the heat flux q, as the nonequilibrium variable in energy e. If one differentiates the nonequilibrium entropy s with respect to the energy holding q constant, then a quantity T(C) of Jou et al ([3]) follows, which differs from the reciprocal of the related temperature Teq by a term quadratic in q. In general, "nonequilibrium temperatures" (understood as the fifth moment of the nonequilibrium density functions) are not the measures of mean kinetic energy. The knowledge of inertial coefficients, such as g, from statistical considerations helps to calculate two basic quantities in the model of heat transfer with finite wave speed. They are: thermal relaxation time τ and and the propagation speed, c0. Of several formulae available that link τ and g, probably the following expression
(19)
[6]. For an ideal gas, Eq. (18) yields θ as θ = mk B−1
(20)
We can now set a variational model of the heat problem.
5 Approaches Ajoining Constraints to a Kinetic Potential For a heat conduction described in an Eulerian frame by the Cattaneo equation and conservation law for the internal energy, the constraints are q ∂q + 2 + ∇ρ e = 0 2 c0 ∂t c0τ
(21)
∂ρ e + ∇.q = 0 , ∂t
(22)
and
where the density of the thermal energy ρe satisfies dρe = ρcvdT, c0 is propagation speed for the thermal wave, τ is thermal relaxation time, and D=c02τ is the thermal diffusivity. Equation (22) assumes the conservation of thermal energy (rigid medium and neglect of the viscous dissipation). For simplicity we assume constant values of involved fields at the boundary. We ignore the vorticity properties of the heat flux. The energy-representation of the Cattaneo equation,
(16)
is most useful ([6], p. 199). It links thermal relaxation time τ with thermal conductivity k, inertia g and state parameters of the system. As, by definition, the propagation speed of the thermal wave c0=(a/τ)1/2, where a= k /(ρcp) is thermal diffusivity, the quantity c0 may be determined from the useful formula
ISSN: 1991-8747
T = cp g
c0 = ( )
and
τ = kg ( ρT ) −1
Stanislaw Sieniutycz, Piotr Kuran
∂js cs2∂t
+
js cs2τ
+ ∇T = 0
(23)
uses diffusive entropy flux js instead of heat flux q. The coefficient cs is defined as cs ≡ ( ρcvθ −1 )1 / 2 (24)
762
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
where θ=T c0−2 , and thermal diffusivity k ≡ ρcv c02τ . Equation (23) is Kaliski’s equation [6]. For an incompressible medium one may apply this equation in the form ∂js c02∂t
+
js c02τ
+ ∇ρ s = 0
q=
.
(28) ∂φ . ∂t
ρ e = −∇ψ −
(25)
(29)
These equations enable one to transfer variational formulation to the space of Lagrange multipliers. For the accepted structure of L, the action A, Eq. (26), in terms of the adjointsψ and φ is
)
2 1 ∂ψ ψ 2 2 1 ∂φ 1 2 A = ε 2 − + c0 ∇φ − ∇.ψ + − ε dVdt.. 2 ∂t 2 2c0 ∂t τ t1,V (30) Its Euler equations with respect to ψ and φ are t2
∫
t2
∫
∂ψ ψ − + c02∇φ ∂t τ
and
which uses the entropy density ρs as a field variable. Yet, we focus here on the action and extremum conditions in the entropy representation (Eqs. (21) and (22) in variables q and ρe). For Eqs. (23) - (25) another approach will be developed in a complementary paper. Action approaches should be distinguished from entropy-production approaches [6], [7]. Here an action is assumed that absorbs constraints (21) and (22) by Lagrange multipliers, the vector ψ and the scalar φ 1 q2 1 2 1 A = ε −1{ 2 − ρe − ε 2 2 c0 2 2 t ,V
Stanislaw Sieniutycz, Piotr Kuran
(26)
−1
1
∂ρ ∂q q + ψ.( 2 + 2 + ∇ρe ) + φ( e + ∇.q)}dVdt. ∂t c0 ∂t c0τ
∂ 1 ∂ψ ψ 2 1 ∂ψ ψ 2 ∂φ − +c0 ∇φ+ 2 − +c0 ∇φ −∇∇.ψ+ = 0 ∂t c02 ∂t τ ∂ t ∂t τ τc0
As kinetic potentials can be diverse, the conservation laws for energy and momentum substantiate the form (26). In Eq. (26), ε is the energy density at an equilibrium reference state, the constant which ensures action dimension for A, but otherwise is unimportant. Yet we assume that the actual energy densityρ e is close to ε, so that the variable ρe can be identified with the constant ε in suitable approximations. We call the multiplier-free term of the integrand of Eq. (26) L≡
1 −1 q 2 ε { 2 − ρe 2 − ε 2} 2 c0
and ∂ ∂φ ∂ψ ψ − ∇.ψ + + ∇. − + c02∇φ = 0. ∂t ∂t ∂t τ
(32)
It is easy to see that (31) and (32) are the original equations of the thermal field, eqs. (21) and (22), in terms of potentials ψ and φ. Their equivalent form below shows the damped wave nature of the transfer process. In fact, Lagrange multipliers ψ and φ of this problem satisfy certain inhomogeneous wave equations. In terms of the modified quantities Ψ and Φ satisfying 2 Ψ = ψ τc02 and Φ = − φτ c 0 these equations are
(27)
the kinetic potential of Hamilton type for heat transfer. It is based on the quadratic form of an indefinite sign, and it has usual units of the energy density. Not far from equilibrium, where ρe is close to ε, two static terms of L yield altogether the density of thermal energy, ρe. Indeed, in view of admissibility of the approximation ρe =ε in Eqs. (27), the kinetic potential (27) represents - in the framework of the linear heat theory - the Hamiltonian structure of a difference between “kinetic energy of heat”, and the nonequilibrium internal energy, ρe. To secure proper conservation laws, no better form of L was found in the entropy representation. The theory obtained in the present case is a linear one. Vanishing variations of action A with respect to multipliers ψ and φ recover constraints, whereas those with respect to state variables q and ρe yield representations of state variables in terms of ψ and φ . For the accepted Hamilton-like structure of L,
ISSN: 1991-8747
(31)
∇ 2Ψ −
∂ 2Ψ
(33)
∂ 2Φ ∂Φ + = ρe . c02∂t 2 τc02∂t
(34)
2
+
∂Ψ
= q.
c0 ∂t 2
τc0 2∂t
and ∇ 2Φ −
As both state variables (q, ρe) and adjoints (ψ ψ, φ) appear herein, they represent mixed formulations of the theory.
6 Canonical Conservation Laws The energy-momentum tensor is defined as G jk ≡ ∑ l
∂vl ∂χ j
∂Λ − δ jk Λ k χ ∂ ∂ v ∂ l
(
)
(35)
where δjk is the Kronecker delta and χ = (x, t) comprises the spatial coordinates and time. Our approach here follows those in [8] and [9], where components of Gjk
763
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
are calculated for a reversible Λ whose gauged form is obtained from the reversible limit of Eq. (26) at ψ=0 by use of the divergence theorem and the differentiation by parts. In our problem 1 ∂φ 1 1 Λ = ε −1{ ( ) 2 − (c0∇φ) 2 − ε 2 } 2 ∂t 2 2
∂S v q q2 q2 + ∇.( ) = ≡ = aj 2s , 2 2 ∂t T ετc 0 T kT
(36)
ρe α q ≅ c0− 2qα ε
(37)
Sv = ρ s −
Clearly, Γ vanishes in the Fourier’s case (c0 → ∞ ). The stress tensor Tab has the form T αβ = −ε −1c0−2qα qβ − δ αβ Λ
(39)
The energy density follows as the Legendre transform of the Lagrangian Λ with respect to rate change of φ 1 q2 1 2 1 2 1 q2 + ρ + ε } ≅ + ρe e 2 c02 2 2 2 εc02
(40)
Finally, density of energy flux follows as
− ∇.q =
G 4 β = Q β = ε −1 ρ e q β ≅ q β .
(41) The associated conservation laws for the energy and momentum show the role of thermal inertia effects
(
1 ∂ ( ε −1c0− 2q 2 + ρ e ) / ∂t = −∇. ε −1ρ eq 2
{
∂(c 0−2 ε −1 ρ e q α ) = ∇. ε −1 (−c0−2 q α q β + δ αβ ∂t
(
1 2 −2 q c0 2
)
− 12 ρ e2 + 12 ε 2
q = (1 − 2a ′φc 02 ) −1 (
(43)
A=
t2
∫ε
t1,V
∂ψ ψ − + c 02 ∇φ) ∂t τ
(47)
)
1 2 ψ 2 − 2 ∂ψ − + c 0 2 ∇φ dVdt 2 (1 − 2a ′φc 0 ) 2c 0 ∂t τ
−1
−
Calculating the four-divergence of the entropy flow ( ∇, ∂ / ∂t ) and using the global conservation law for the total energy E we obtain the entropy balance .
(46)
Substituting Eqs. (29) and (47) into action A of Eq. (26) (L of Eq. (27)) shows that in terms of the potentials the action acquires the form
7 Satisfaction of Second Law
∂Sv q 1 ∂q q ∂q + ∇.( ) = q. ∇T −1 − = . − k∇T − τ 2 kT 2 T t ∂t ∂ ∂t ε Tc 0
and
The energy representation (29) is unchanged, but the heat flux representation follows in a generalized form
These are canonical conservation laws which are obtained on the reversible paths. Now the important question is if these laws can be recovered for irreversible processes in which the entropy is produced and the Second Law holds.
2 ∂φ 1 2 −1 1 ε ∇ . ψ + + ε dVdt . ∫ 2 ∂ t 2 t1,V t2
(48)
However the Euler-Lagrange equations for this action are not the process constraints in terms of potentials, i.e. the method fails to provide a correct variational formulation for constraints with sources. The way to improve the situation is to substitute the obtained representations to a transformed augmented action in
(44)
Applying Cattaneo equation we arrive at the second law
ISSN: 1991-8747
∂ρ e − a' q 2 ∂t
(42)
))} .
energy
However, while simple and useful, the method of construction of a suitable action A in the space of potentials by the direct substitution of the representation equations to the kinetic potential L is limited to cases with linear constraints that do not contain sources. This may be exemplified when the internal energy balance contains a source term a’q2, where a’ is a positive constant. The augmented action integral, generalizing Eq. (26), should now contain the negative term - a’q2 in its φ term.
whence (after substituting the stationary Lagrangian)
E = G 44 = ε −1{
τ q2 2 2λ T
8 Sources of internal conservation laws
(38)
1 1 1 T αβ = −ε −1c0−2qα q β − δ αβ ε −1{ ρe2 − q2c0−2 − ε 2} 2 2 2
(45)
where a = k-1 is the thermal resistance. It should be kept in mind that the density Sv appearing in the above equations is not identical with the density of classical equilibrium entropy ρs. As shown by Eq. (61) in Sec.8, Sv is rather the density of the nonequilibrium entropy which contains the heat flux as an extra variable and appears in extended irreversible thermodynamics (EIT; [3]).
Below the results for tensor G = Gjk are discussed in physical variables. The momentum density of heat is Γα = −Gα 4 = c0− 2
Stanislaw Sieniutycz, Piotr Kuran
764
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
Stanislaw Sieniutycz, Piotr Kuran
Somewhat surprisingly, it follows that a source term in the internal energy balance, as in Eq. (46), should be the necessary property of the Cattaneo model, else the energy conservation will be violated. Indeed, aimed at the evaluation of energy conservation we multiply the non-truncated Cattaneo formula
which the only terms rejected are total time or space derivatives. The latter can be selected via partial differentiation within the integrand of the original action A. (As we know from the theory of the functional extrema the addition of subtraction of terms with total derivatives and divergences do not change extremum properties of a functional.) When this procedure is applied to the considered problem and total derivatives are rejected, a correct action follows in the form
∂q
c02 ∂t
+
q c02τ
+ ∇ρ e = 0
(21)
by heat flux q. The result is an equation A=
t2
∫ε
)
1 2 ψ 2 −1 ∂ψ − + c 0 2 ∇φ dVdt 2 (1 − 2a ′φc 0 ) τ 2c 0 ∂t
−1
t1,V
∂q 2 2εc 02 ∂t
2 ∂φ 1 − ∫ ε ∇.ψ + + ε 2 dVdt (49) ∂t 2 2 t1,V t2
−1 1
t2
−1 2 ∫ ε {(1 − 2a' φc0 )
t1 ,V
∂ρ e + ∇.q = 0 , ∂t
(51)
(22)
leads to a differential result ∂q 2 q2 ρe ∂ρe −1 + + ∇ .( q ) = − ε ρ e εc02τ 2εc02∂t ε ∂t
(52)
which – under the linear approximations of the present theory – yields not a conservation law for the total energy, but a balance formula with an energy source q 2 . (53) ∂ q 2 −1 ε ρ + + ρ ∇ .( q ) = − e e ∂t 2εc 02 εc 02τ
q2 1 2 1 2 − ρ e 0 + ε }dVdt (50) 2 2c02 2
This shows violation of the total energy conservation for the source-less internal energy of the model, and leads to the conclusion that the model composed of the Cattaneo equation and source-less balance of internal energy is physically admissible only in the reversible case of an infinite τ. Certainly, this is not a demanded property of the energy transfer model, thus a further analysis is required. The solution of the dilemma seems to admit a properly large, yet a non-vanishing, source in the internal energy balance. Admitting a source of the internal energy, as in equation (46), and using Eq. (46) in Eq. (51) we find
This form of A shows that, when original state space is involved, the state space required to accommodate the variational principle must be enlarged by inclusion of the Lagrange multiplier φ as an extra variable. In fact, Eq. (50) proves that original state space (“physical space”) is lacking sufficient symmetry (Vainberg’s theorem [6]). Yet, as Eq. (50) shows, the adjoint space of potentials (ψ ψ, φ), while also four-dimensional as space (q, ρe0), can accommodate the variational formulation. Why is this so? Because the representation equations do adjust themselves to the extremum requirement of A at given constraints, whereas the given constraints without controls cannot exhibit any flexibility.
ISSN: 1991-8747
ρe q2 ∇.q + ∇.(ε −1qρ e ) = − 2 ε εc 0 τ
which describes an energy balance We observe that the combination of this balance with the sourceless balance of internal energy
This form differs from that of Eq. (34) only by the power of the coefficient containing the constant a’, related to the source. With the related representation equations (29) and (47), action (49) yields the proper Cattaneo constraint (21) and the generalized balance of internal energy which extends equation (22) by the source term a’q2. Equation (49) proves that four-dimensional potential space (ψ ψ, φ) is sufficient to accommodate an exact variational formulation for the problem with a source. Yet, due to the presence of this source, the formulation does not exist in the original four-dimensional original space (q, ρe), and, if somebody insists to exploit this space plus possibly a necessary part of the potential space, the following action is obtained from Eqs. (29), (47) and (49) A=
−
ρ ∂ρ ρ ∂q2 q2 + e e − e a' q2 + ∇.(ε −1qρe ) = − 2 . 2 ε εc0τ 2εc0 ∂t ε ∂t
From this formula we observe that for a’ satisfying
765
Issue 8, Volume 3, August 2008
(54)
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
a' =
1 1 1 = = . 2 ρ e 0 c0 τ ρcv DhT kT
(55)
Stanislaw Sieniutycz, Piotr Kuran
The total energy is conserved in the linearized form (56) or in an exact form
and - under the approximation of the linear theory conservation of total energy is satisfied by the sourceless equation ∂ q 2 + ρ e + ∇.(ε −1qρ e ) = 0 2 ∂t 2εc 0
∂ Tτj2s + ∇.(Tj s ) = 0 + ρ e ∂t 2k
A similar scheme of reasoning can be applied in the energy representation where the entropy balance is the process constraint.
(56)
The positive value of a’ in Eq. (55) implies generation of the internal energy in the heat transfer process, as described by Eq. (46). Still, the internal energy equation with the positive source caused by the quadratic heat flux is a dubious structure. Nonetheless, it is a step forward in comparison with sourceless internal energy balance as we may observe that this generation awkwardly mimics the entropy production with the coefficient a=Ta’=1/k. This finally leads us to the conclusion that it is the entropy balance with the source j2/k than that should replace Eq. (46) in the entropy representation. Replacing Eq. (46) with a’ of Eq. (55) by its conserved counterpart i.e. total energy balance expressed in entropy terms
T
∂ρ s = −∇.q ∂t
9 Entropy balance as process constraint Assume that a variational functional can be constructed in the entropy representation for the Kaliski’s equation (25) [10]. This assumption means that after transforming the variational stationarity condition of an action with respect to ψ ∂js c02∂t
∂j s c02 ∂t
c s2 ∂t
and using in this result the Cattaneo equation (21) in the equivalent form
∂j 2s
we obtain
2c s2 ∂t
(60)
js c02 τ
= −ρcv T −1∇T
(63)
(64)
+
js c s2 τ
+ ∇T = 0
(23)
+ ∇.(Tj s ) − T∇.j s = −
j 2s
c s2 τ
.
(65)
Under the assumption of the entropy conservation Eq.(65) would yield an equation
and hence the second law balance for the entropy of extended thermodynamics
∂j 2s ∂ρ s j 2s + ∇ .( T j ) + T = − s ∂t 2c s2 ∂t c s2τ
2
ISSN: 1991-8747
+
After multiplying Kaliski’s equation (23) by js we find
(59)
τ ∂ q q 2 + ∇.(qT −1 ) = ρs − 2 ∂t 2λ T λT 2
(25)
we obtain Kaliski’s equation of entropy transfer with explicit temperature gradient
∂j s
∂ρ s q q τ ∂q = −∇.( ) + q. 2 + ∂t T λT ∂t λT 2
+ ∇ρ s = 0
c s2 = c02 (ρcvT −1 )
(58)
q τ ∂q + ∇T −1 − =0 2 λT ∂t λT 2
js c02τ
and using the relation between coefficients cs and c0
(where ρs is equilibrium entropy density) we expect to obtain a pronounced result. After rearranging
−
+
into the form
(57)
∂ρ s q = −∇.( ) + q.∇T −1 ∂t T
(62)
(61)
thus leading to the following energy balance
766
Issue 8, Volume 3, August 2008
(66)
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
j2s ∂ j 2s + ρ e 0 ) + ∇.(Tj s ) = − 2 . ( ∂t 2c s2 cs τ
thermodynamics. (67)
10 Conclusions The most important physical properties resulting from the obtained generalization are: a finite momentum of heat and the non-classical terms in the stress tensor caused by the heat flux. And with all this one still has both the first and the second law satisfied. Here the product of temperature T and positive entropy source of Eq. (68) (quadratic in js) emerges as the kinetic energy of heat in the energy balance (70) or (71). While this result is simple, its role is nontrivial for the existence of the variational formulation. Let us recall that in the variational dynamics of real fluids it is always extremely difficult to simultaneously satisfy both laws of thermodynamics (classical Hamilton's principle holds only for reversible processes i.e. those without entropy production). The optimizing approach overcomes familiar difficulties resulting from the presence of both odd and even derivatives with respect to time in the differential equations of heat, and, consequently, sets a variational formulation for the irreversible heat conduction. Extremum conditions, Eqs. (33) and (34), show that for given q and ρe heat transfer can be broken down to potentials. This is similar to the case of electromagnetic field theory or gravitation theory, where the knowledge of sources defines the field potentials. An important case is a “ballistic” transfer with τ → ∞ , when undamped thermal waves propagate with speed c0 and satisfy d'Alembert's equation. As shown by Eq. (45) the results are consistent with the second law in an identically satisfied form; this holds in both classical irreversible theory (CIT) and extended irreversible thermodynamics (EIT; [3]). Associated approaches with Lagrange multipliers and Dirac brackets are available [11,12,13]. It is interesting that not only variational calculus but also other optimization methods may also be fruitful in the context of problems considered here, such as, e.g., simulated annealing [14] and dynamic programming [15, 16].
This result implies a non-vanishing source of the energy and means that the assumption of the conserved entropy would result in the violation of the energy conservation law. Yet, for the ‘physical’ entropy satisfying the nonconserved balance
− ∇.j s =
∂ρ s − aj 2s ∂t
(68)
with the coefficient a equal to the reciprocal of thermal conductivity
a=
1 1 = . 2 Tc s τ k
(69)
total energy balance follows from Eqs. (65), (68) and (69) in the form
j2s ∂ j 2s 2 ( + ρ e ) + ∇.(Tj s ) − Taj s = − 2 ∂t 2c s2 cs τ
(70)
which is conservative because the two source terms mutually cancel. This is show in Eq. (71) below. Thus the energy conservation is preserved whenever the coefficient a satisfies Eq. (69). Importantly, even when the conservation laws are satisfied in irreversible processes in their canonical form the related extremum action and potential representations of physical variables do explicitly contain potentials not only their derivatives. We may thus claim that whenever a irreversible process occurs with the coefficient a satisfying Eq. (69), conservation laws in their canonical form (the same as for reversible processes) can be obtained from the Noether’s theorem. In particular the energy conservation law has the form
∂ j2s + ∇.(Tj s ) = 0 . + ρ e 2 ∂t 2c s
Acknowledgements This work was supported by TU Warsaw statute grant Simulation and optimization of thermal and mechanical separation processes in 2008.
(71)
References: [1] H. Grad, “Principles of the theory of gases”, in Handbook der Physik 12, S. Flugge, Ed., Berlin: Springer, 1958. [2] S. Sieniutycz and R. S. Berry, “Conservation laws from Hamilton's principle for nonlocal
As the coefficient cs2=k(Tτ)-1, the above equation is equivalent with Eq. (62). In conclusion, the mathematical scheme obtained preserves the conservation of energy and simultaneous production of the entropy, in accordance with laws of
ISSN: 1991-8747
Stanislaw Sieniutycz, Piotr Kuran
767
Issue 8, Volume 3, August 2008
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS
Stanislaw Sieniutycz, Piotr Kuran
thermodynamic equilibrium fluids with heat flow”, Phys. Rev. A. 40 (1989) 348-361. [3] D. Jou, J. Casas-Vazquez and G. Lebon, Extended Irreversible Thermodynamics, Heidelberg: Springer, 1993. [4] H. Callen, Thermodynamics and an Introduction to Thermostatistics, New York: Wiley, 1998. [5] S. Sieniutycz, “Thermodynamics of coupled heat, mass and momentum transport with finite wave speed. I.,II.Intern. J. Heat Mass Transfer 24 (1981) 1723-1732 and 24 (1981) 1759-1769. [6] S. Sieniutycz, Conservation Laws in Variational Thermo-Hydrodynamics, Dordrecht:, Kluwer, 1994. [7] S. Sieniutycz, “Variational thermomechanical processes and chemical reactions in distributed Systems”, Intern. J. Heat and Mass Transfer, 40 (1997) 3467-3485. [8] J.J. Stephens, “Alternate forms of the Herrivel-Lin variational principle”, Phys. Fluids 10 (1967) 76-77. [9] R. L. Seliger and G. B. Whitham, “Variational principles in continuum mechanics”, Proc. Roy. Soc. 302A (1968) 1-25. [10]. S. Sieniutycz and H. Farkas, eds. Variational and Extremum Principles in Macroscopic Systems, Oxford, Elsevier, 2005. (See especially Chap.7 of part II.) [11]G.L.Liu, The first exact variational formulation of 3D Navier-Stokes, Comp. Mech. VI, (2004) 438-443. [12]G.L.Liu, A unified variational formulation of aeroelasticity in 3D unsteady transonic flow, IUTAM Symp.Transsonicum IV (in Gottingen), H. Sobieczky (ed.), Kluwer Acad. Publ. (2003) 85-90. [13] S. Nguyen and Ł. A. Turski, Examples of the Dirac approach to dynamics of systems with constraints, Physica A 290 (3-4) (2001), 431-444. [15] H.G Arango and G.L. Tores, Spatial electric load distribution forecasting using simulated annealing, WSEAS Transactions on Systems, 3(1) 2004, 14-19. [16] M.T. Suzuki, A dynamic programming approach to search similar portions of 3D models, WSEAS Transactions on Systems, 3(1) 2004, 14-19. [17] S. Sieniutycz, Optimization of multistage electrochemical systems of fuel cell type by dynamic programming, WSEAS Transactions on Mathematics, 2(3) 2003, 228-232.
ISSN: 1991-8747
768
Issue 8, Volume 3, August 2008