Comprehensive Criteria for the Extrema in Entropy Production Rate for Heat Transfer in the Linear Region of Extended Thermodynamics Framework
Comprehensive Criteria for the Extrema in Entropy Production Rate for Heat Transfer in the Linear...
Verros, George D.
2020-10-08 00:00:00
axioms Article Comprehensive Criteria for the Extrema in Entropy Production Rate for Heat Transfer in the Linear Region of Extended Thermodynamics Framework George D. Verros Department of Chemistry, Aristotle University of Thessaloniki, P.O. Box 454 Plagiari, 57500 Epanomi, Greece; gdverros@yahoo.gr Received: 18 September 2020; Accepted: 5 October 2020; Published: 8 October 2020 Abstract: In this work comprehensive criteria for detecting the extrema in entropy production rate for heat transfer by conduction in a uniform body under a constant volume in the linear region of Extended Thermodynamics Framework are developed. These criteria are based on calculating the time derivative of entropy production rate with the aid of well-established engineering principles, such as the local heat transfer coecients. By using these coecients, the temperature gradient is replaced by the dierence of this quantity. It is believed that the result of this work could be used to further elucidate irreversible processes. Keywords: entropy production; non equilibrium thermodynamics; criteria 1. Introduction Non equilibrium thermodynamics is a rapidly growing branch of chemical physics with many applications in science and engineering, from power engineering to environmental sciences, from chaos to complex systems, and from life sciences to nanosciences. The fundamentals of non-equilibrium thermodynamics are given in several text books [1–5] and reviews [6–8]. One of the most fundamental aspects presented in these works is the calculation of entropy production rate which is strongly related to irreversible processes such as heat transfer or diusion. In 1912, Ehrenfest introduced the fundamental question, in relation to the non-equilibrium stationary states, on the existence of a function that achieves its extreme value as does entropy for stationary states in equilibrium thermodynamics. Today, it is evident in the literature [1–4] that entropy production rate in non-equilibrium thermodynamics plays the same role as entropy does for stationary states in equilibrium thermodynamics. Therefore, it is crucial to establish criteria for detecting the extrema (minimum or maximum) of entropy production rate. This has been the subject of intensive research in recent years as is reviewed by many workers in the field. More specifically, Onsager [9–12], by assuming linear phenomenological relations between fluxes (J ) and forces (X ) (Linearity Axiom: X = S R J ) based on statistical mechanics arguments, proved that i i i i ij the R coecients are symmetric (Onsager Reciprocal Relations, ORR). Moreover, Onsager introduced ij the least dissipation principle. This principle reduces to maximization of entropy production at fixed thermodynamic forces. Alternatively, the Prigogine theorem [5,6] states, for systems in the linear regime (both the linear phenomenological relations and ORR hold true) that total internal entropy production reaches a minimum value at non equilibrium stationary states. The Ziegler principle [13,14] states if the thermodynamic forces (X ) are preset, then true thermodynamic fluxes (J ) satisfying the entropy production rate equation ( = S J X ) give the i i i maximum value of the entropy production rate ((J)). Axioms 2020, 9, 113; doi:10.3390/axioms9040113 www.mdpi.com/journal/axioms Axioms 2020, 9, 113 2 of 7 Many authors [15–22] have developed dierent approaches, but criticism [23–26] has been made of each of them. Please note that criteria for detecting the minimum or maximum [27] or alternative approaches [28] are available as a possible answer to the Ehrenfest question. On the other hand, Extended Thermodynamics is a powerful tool for investigating physicochemical processes outside the region of local equilibrium. The aim of this work is not only to establish simple criteria for entropy production rate extrema in the linear region of Extended Thermodynamics Framework, but also to introduce in the area the well-established engineering concept of heat transfer coecients. 2. Theoretical Part and Results Irreversible thermodynamics close to equilibrium is based on three independent axioms above and beyond those of equilibrium Thermodynamics [1–6]: (1) The equilibrium thermodynamic relations apply to systems that are not in equilibrium, provided that the gradients are not too large (quasi-equilibrium axiom). (2) All the fluxes (J ) in the system may be written as linear relations involving all the thermodynamic forces, X (linearity axiom, X = R J ; i = 1, 2, ..., K). i i ij j=1 (3) In the absence of magnetic fields and assuming linearly independent fluxes or thermodynamic forces the matrix of coecients in the flux–force relations are symmetrical. This axiom is known as the Onsager Reciprocal Relations (ORR): R = R . ij ji The starting point of this work is the entropy (s) per unit mass balance in a local form [1–6]: 0 1 B C B C B J J C B q i C B C ds B C i=1 B C B C +r = ; = XJ (1) B C B C dt B T C B C B C @ A where is the mass density, symbol stands for the entropy production rate per unit volume, T is the absolute temperature, J represents the heat flux, t is time, stands for the chemical potential of i-th substance, J is the i-th substance molar flux defined relative to the center of molar mass and K is the total number of substances participating in the diusion. In the Extended Thermodynamics framework, the axiom of local equilibrium is abandoned. The entropy production term per unit volume () for heat transfer by conduction in a uniform body under constant volume in the extended thermodynamics area is written as [5,6,23]: 0 1 @J B C rT q B C B C = B a CJ (2) @ A q @t The Cattaneo equation is directly derived from the above equation by further applying the linearity axiom [5,6,23]: 0 1 @J @J B C q rT q B C 2 2 B C J = L B a C ; = J + krT ; k = L /T ; = L a; a = /L = /kT (3) qq qq qq qq q @ A q @t @t All the symbols are explained in the nomenclature section. The relaxation time for heat transfer by conduction in a uniform body under constant volume is a positive constant [5,6,23]. A detailed review of the Cattaneo equation is given elsewhere [5,6,23]. In this work, the heat flux in the extended thermodynamics area is written in terms of the local heat transfer coecient close to equilibrium (h ) and the residual local heat transfer coecient (h ) res loc is written as: J = h h (T T ) (4) res q loc 0 Axioms 2020, 9, 113 3 of 7 The residual local heat transfer coecient may be parameterized as a positive quantity. The local heat transfer coecient (h ) close to equilibrium is a vector defined as: loc @T k = h (T T ) ; j = 1, 2, 3 (5) loc,j 0 @x If h = 1, then = 0 and the Fourier law is directly obtained from the Cattenao equation. res The reference temperature T is defined in such a way that h h (T T ) > 0 for positively defined 0 loc res 0 heat flux. In this way the temperature gradient is replaced by the temperature dierence. This is not a new idea; the origin of this idea could be found in many textbooks [1,29] as the definition of the local heat transfer coecients close to equilibrium (h ). These coecients in the most general case loc,j involving an industrial process are functions not only of the process conditions, but also of the reference temperature T . They are calculated in terms of dimensionless groups such as the Nusselt number by using experimental correlations of other dimensionless groups such as the Prandtl number [1,29]. The main advantage of using these coecients is that one could use these coecients at dierent scales based on the similarity principle [1] without resorting to analytical or numerical solutions to calculate the heat flux. Moreover, by assuming that the solution to the Cattaneo equation could be written as T T = C(t) C x and by replacing this solution to Equation (5), one could directly show that the local 1j j j=1 heat transfer coecient close to equilibrium (h ) is independent of time. Furthermore, by introducing loc the heat transfer coecients into the Cattaneo equation and by further taking into account that the local heat transfer coecient close to equilibrium is independent of time, it can be directly shown that the following equation holds true: dh (T T ) dh C(t) res 0 res = (1 h )(T T ) or = (1 h )C(t) res 0 res dt dt (6) d ln C(t) dh res or = h + (1 h )/ res res dt dt The above equation shows that the residual local heat transfer coecient (h ) is independent of res position and it is function of time. By neglecting viscous dissipation and by assuming an absence of external forces acting on the system, the time derivative of temperature is calculated as [1–4]: @T 1 1 h kh res res = r J = r h h (T T ) = r h (T T ) = r T (7) loc res 0 loc 0 @t c c c c v v v v where c represents the specific heat capacity per unit mass under a constant volume. The above equation in terms of the solution T T = C(t) C (x ) is written as: 0 1j j j=1 0 1 B C @C(t) C (x ) B 1j j C B C B C j=1 h C(t) res B C B C = r h C (x ) B C loc 1j j B @t C v B C j=1 @ A (8) @C(t) d ln(C(t)) or = h C(t) ; = h res res @t dt where is a constant having the inverse of time as units. Since h >0 the condition for a bounded res solution to Equation (8) requires > 0. Based on the above result Equation (6) is re-written as: dh res = h + (1 h )/ (9) res res dt Axioms 2020, 9, 113 4 of 7 The above equation is a constitutive equation for the residual local heat transfer coecient. The entropy production P under constant volume by using the definitions of heat transfer in coecients and the Cattaneo equation is given as: R R R 2 2 2 P = dV = J r1/T a@J /@t dV = h h (T T ) /kT dV = in q q res 0 loc V V V (10) 2 2 = kh (rT) /T dV res According to the above equation the entropy production per unit volume P equals to zero for in uniform temperature (rT = 0). This result for the linear region is in accordance with the literature (p. 41, ref [23]). However, the heat flux calculated by the Cattaneo equation can be non-zero in the case of uniform temperature (rT = 0). Therefore, in the most general case of the nonlinear region the entropy production per unit volume is not zero [30,31], even in the case of uniform temperature (rT = 0). One can directly show that if the derivative with respect to time of entropy production P inside in a non-equilibrium thermodynamic system with a constant volume is greater than zero or in other dP in words, = dV > 0, then entropy production monotonically increases with time. Therefore, dt dt the achieved value for entropy production is at a maximum and vice versa for dP /dt <0 [1–4]. in The derivative of entropy production P with respect to time for heat transfer under constant in volume is given as: Z Z Z d J . r1/T a@J /@t d J .(r1/T (1 h )r1/T) res q q q dP d in = dV = dV = dV (11) dt dt dt dt V V V In the derivation of the above equation the definitions of heat transfer coecients as well as the Cattaneo equation were used. By using the identity r (A f ) = fr A + A (r f ) and by taking into account that the residual heat transfer coecient is independent of position the above integral is written as: R R d J r(h /T) d J h r(1/T) ( res ) ( res ) dP q q in = dV = dV = dt dt dt V V ( ) (12) d(rJ h /T) d(h /TrJ ) q res res q = dV dt dt By using Gauss’s theorem, the first term of the above integral can be written as surface integral. Since boundary conditions are assumed to be time-independent on the boundary [4,5], this surface integral vanishes, and the above equation is written as: d h /Tr J res dP q in = dV (13) dt dt Moreover, the following equation is directly obtained by further using the above equation and Equation (7) as well as the chain rule for partial derivatives: ! ! d h /Tr J res q @T @T ( ) = c h /T + c dh /dt /T (14) v res v res dt @t @t If h = 1 then = 0 then the Fourier law is valid and the above equation is reduced to the res well-known criterion for entropy production extrema for heat transfer close to equilibrium found in Axioms 2020, 9, 113 5 of 7 many textbooks [1–4]. The second term on the right-hand side of the above equation by further using Equation (9) is analyzed as: @T 1 @T (dh /dt)/T = h + (1 h )/ (15) res res res @t T @t Finally, the criterion for minimum entropy production is formulated by combining Equations (14) and (15) with Equation (11) as: c c dP v @T v @T in 2 = h + h + (1 h )/ dV 2 res res res dt @t T @t n o (16) R R @T 2 1 @T = h dV + c h + (1 h )/ dV res v res 2 res T @t @t V V dP 2 @T in ( ) In the case that h + 1 h / < 0 for both time and space coordinates then < 0; res res @t dt one could anticipate that a minimum is reached for the entropy production of the whole system. In the 2 @T ( ) opposite case dP a maximum in entropy production could be attained if h + 1 h / > 0 in res res @t > 0 dt R R n o @T 2 1 @T and h dV + c h + (1 h )/ dV > 0 res v res 2 res T @t @t V V 3. Conclusions Comprehensive criteria for the extremum value of entropy production in the linear region of Extended Thermodynamics Framework were developed in this work. The above criteria are based on the calculation of the entropy production rate derivative with respect to time by assuming that heat flux could be approximated by the local heat transfer coecients. These criteria were applied for the heat transfer by conduction in a uniform body under constant volume in the Extended Thermodynamics Framework. For this purpose, the residual heat transfer coecients are introduced in the field and the temperature gradient is replaced by the dierence in temperature. It is believed that the results of this work might be used to further elucidate irreversible processes. Criteria for the extremum value of entropy production for heat transfer in the linear region of Extended Thermodynamics Framework were developed Introduction of local heat transfer coecients in the field The advantages of using local heat transfer coecients are (a) calculation of heat flux without resorting to analytical or numerical solutions and (b) temperature gradient is replaced by temperature dierence. Funding: This research received no external funding. Acknowledgments: The author is thankful to K. Somerscales for her help in the preparation of the manuscript. Conflicts of Interest: The author declares no conflict of interest. Nomenclature a thermodynamic parameter c specific heat capacity for constant volume h close to equilibrium local heat transfer coecient loc h residual local heat transfer coecient res J flux k thermal conductivity L phenomenological coecients relating fluxes with thermodynamic driving forces P entropy production rate inside the whole system in s entropy Axioms 2020, 9, 113 6 of 7 T absolute temperature T reference absolute temperature t time V volume x space coordinate X thermodynamic driving force Greek Letters constant chemical potential of i-th substance mass density entropy production rate per unit volume relaxation time References 1. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. Transport Phenomena, 2nd ed.; Wiley: New York, NY, USA, 2002; pp. 764–798. 2. de Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics, 2nd ed.; Dover Publications: New York, NY, USA, 1984; pp. 11–82. 3. Kuiken, G.D.C. Thermodynamics of Irreversible Processes—Applications to Diusion and Rheology., 1st ed.; Wiley: New York, NY, USA, 1994; pp. 1–135. 4. Demirel, Y. Nonequilibrium Thermodynamics, Transport and Rate Processes in Physical, Chemical and Biological Systems, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2007; pp. 1–144. 5. Kondepudi, D.; Prigogine, I. Modern Thermodynamics: From Heat Engines to Dissipative Structures, 2nd ed.; Wiley: New York, NY, USA, 2015; pp. 333–455. 6. Kondepudi, D.P.; Prigogine, I. Thermodynamics Nonequilibrium, in Encyclopedia of Applied Physics, 2nd ed.; Trigg, G.L., Ed.; Wiley-VCH: Weinheim, Germany, 2003. 7. Demirel, Y.A.; Sandler, S.I. Nonequilibrium Thermodynamics in Engineering and Science. J. Phys. Chem. B 2004, 108, 31–43. [CrossRef] 8. Michaelides, E.E. Transport Properties of Nanofluids. A Critical Review. J. Non Equil. Thermodynamics 2013, 38, 1–79. [CrossRef] 9. Onsager, L. Reciprocal Relations in Irreversible Processes I. Phys. Rev. 1931, 37, 405–426. [CrossRef] 10. Onsager, L. Reciprocal Relations in Irreversible Processes II. Phys. Rev. 1931, 37, 2265–2279. [CrossRef] 11. Onsager, L.; Fuoss, R.M. Irreversible Processes in Electrolytes: Diusion, Conductance, and Viscous Flow in Arbitrary Mixtures of Strong Electrolytes. J. Phys. Chem. 1932, 36, 2659–2778. [CrossRef] 12. Onsager, L. Theories and Problems of Liquid Diusion. Ann. N. Y. Acad. Sci. 1945, 46, 241–265. [CrossRef] [PubMed] 13. Ziegler, H. An Introduction to Thermomechanics; North-Holland: Amsterdam, The Netherlands, 1983. 14. Martyushev, L.M.; Seleznev, V.D. Maximum Entropy Production Principle in Physics, Chemistry and Biology. Phys. Rep. 2006, 426, 1–46. [CrossRef] 15. Truesdell, C.A. Rational Thermodynamics; Springer: New York, NY, USA, 1984. 16. Muschik, W.; Ehrentraut, H.; Papenfuss, C. Mesoscopic Continuum Mechanics. In Geometry, Continua and Microstructure; Collection Travaux en Cours 60; Maugin, G., Ed.; Herrman: Paris, France, 1999. 17. Öttinger, H.C. Beyond Equilibrium Thermodynamics; Wiley: Hoboken, NJ, USA, 2005. 18. Sieniutycz, S.; Farkas, H. Variational and Extremum Principles in Macroscopic Systems; Elsevier: Oxford, UK, 2005. 19. Bejan, A. Entropy Generation Minimization; CRC Press: Boca Raton, FL, USA, 1995. 20. Bejan, A. Shape and Structure, from Engineering to Nature; Cambridge University Press: Cambridge, UK, 2000. 21. Bejan, A.; Lorente, S. The Constructal Law of Design and Evolution in Nature. Phil. Trans. R. Soc. B 2010, 365, 1335–1347. [CrossRef] [PubMed] 22. Andresen, B. Finite-time thermodynamics. In Finite-Time Thermodynamics and Thermoeconomics, Advances in Thermodynamics; Sieniutycz, S., Salamon, P., Eds.; Taylor and Francis: New York, NY, USA, 1990. 23. Jou, D.; Lebon, G.; Casas-Vazquez, J. Extended Irreversible Thermodynamics, 3rd ed.; Springer: Berlin, Germany, 2001; pp. 39–70. Axioms 2020, 9, 113 7 of 7 24. Grinstein, G.; Linsker, R. Comments on a Derivation and Application of the “Maximum Entropy Production” Principle. J. Phys. A Math. Theory 2007, 40, 9717–9720. [CrossRef] 25. Jaynes, E.T. The Minimum Entropy Production Principle. Ann. Rev. Phys. Chem. 1980, 31, 579–601. [CrossRef] 26. Martyuchev, L.M. Entropy and Entropy Production: Old Misconceptions and New Breakthroughs. Entropy 2013, 15, 1152–1170. [CrossRef] 27. Di Vita, A. Maximum or Minimum Entropy Production? How to Select a Necessary Criterion of Stability for a Dissipative Fluid or Plasma. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2010, 81, 041137. [CrossRef] [PubMed] 28. Lucia, U.; Grazzini, G. The Second Law Today: Using Maximum-Minimum Entropy Generation. Entropy 2015, 17, 7786–7797. [CrossRef] 29. Slattery, J.C. Advanced Transport Phenomena, 1st ed.; Cambridge University Press: Cambridge, UK, 1999; pp. 283–366. 30. Li, S.-N.; Cao, B.-Y. Generalized variational principles for heat conduction models based on Laplace transform. Int. J. Heat Mass Transf. 2016, 103, 1176–1180. [CrossRef] 31. Li, S.-N.; Cao, B.-Y. On Entropic Framework Based on Standard and Fractional Phonon Boltzmann Transport Equations. Entropy 2019, 21, 204. [CrossRef] © 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png
Axioms
Multidisciplinary Digital Publishing Institute
http://www.deepdyve.com/lp/multidisciplinary-digital-publishing-institute/comprehensive-criteria-for-the-extrema-in-entropy-production-rate-for-Pg0k30hyc1