arxiv: v1 [cond-mat.stat-mech] 21 Mar 2016

Transcription

1 Efficiencies of power plants, quasi-static odels and the geoetric-ean teperature Raandeep S. Johal arxiv: v1 [cond-at.stat-ech] 21 Mar 2016 Departent of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, S.A.S. Nagar, Manauli PO , Punjab, India Observed efficiencies of industrial power plants are often approxiated by the square-root forula: 1 T /T +, where T + (T ) are the highest (lowest) teperature achieved in the plant. This expression can be derived within finite-tie therodynaics, or, by entropy generation iniization, based on finite rates of processes. A closely related quantity is the optial value of the interediate teperature for the hot strea, which is given by the geoetric-ean value: T + T. It is proposed to odel the operation of plants by quasi-static work extraction odels, with one reservoir (source/sink) as finite, while the other as practically infinite. No siplifying assuption is ade on the nature of the finite syste. This description is consistent with two odel hypotheses, each yielding a specific value of the interediate teperature. We show that the expected value of the interediate teperature, defined as the arithetic ean, is very closely given by the geoetric-ean value. The definition is otivated as the use of inductive inference in the presence of liited inforation. Keywords: Power plants; Heat Engines; Theral efficiency; Inference PACS nubers: a, Ce, Ln rsjohal@iiserohali.ac.in

2 2 Introduction: In recent years, there has been a great interest in extending therodynaic odels to justify the observed perforance of industrial power plants [1 4]. observed efficiencies are usually uch less than the Carnot liit The η C = 1 T T +, (1) evaluated on the basis of the highest (T + ) and the lowest (T ) teperatures, for the particular plant. Naturally, real achines operate under irreversibilities caused by various factors, like finite rates of heat transfer and fluid flow, internal friction, heat leakage and so on, unlike the idealized quasi-static processes of a reversible cycle. Thus the analysis of irreversible odels with finite-rate processes sees a reasonable goal to pursue. One often-studied easure is the efficiency at axiu power of an irreversible odel which is then copared with the observed efficiency of these plants. The earliest known such odel is ascribed to Reitlinger [5], which involved a heat exchanger receiving heat fro a finite hot strea fed by a cobustion process. An analogous odel was applied to a stea turbine by Chabadal [6]. The considered heat exchanger in these odels is effectively infinite. Novikov [7] considered the heat transfer process between a hot strea and a finite heat exchanger with a given heat conductance. Two, siple but significant, assuptions enter these odels: i) constant specific heat of the inlet hot strea and ii) validity of Newton s law for heat transfer. Further, there appears a floating, teperature variable in between the highest and the lowest values, such as the teperature of the exhaust war strea, over which the output power can be optiized. This yields an optial value of the interediate teperature which is usually found to be T + T, the geoetric ean of T + and T. Related to this fact, is the conclusion that the efficiency at axiu power is given by an elegant expression: T η CA = 1. (2) T + Due to historical iperative [8], the above expression ay be called Reitlinger-Chabadal- Novikov efficiency. However, ore recently it was rediscovered by Curzon and Ahlborn (CA) [1]. Thus in the physics literature, it is ore popularly addressed as CA-value. This latter odel considered finite rates of heat transfer at both the hot and the cold contacts, but also explicitly considered the ties allocated to these contacts. The average power per cycle ay be optiized over these ties [1], or alternately, over the interediate teperature

3 3 variables [9]. This odel spawned uch activity and the new area borne thereforth was tered Finite-tie Therodynaics [10]. In the engineering literature, the corresponding approach is called Entropy Generation Miniization [2, 11]. A positive indication for the siple therodynaic approach is that the actually observed efficiencies of industrial plants happen to be quite close to CA-value. Fig.1 shows this coparison as tabulated in Table I. Although the agreeent is close, the observed values can be higher, or lower, than CA-value. This apparent discrepency has stiulated further extensions of odels, for instance using the low-dissipation assuption [3], which predict the efficiency at axiu power, to be bounded as: η C 2 η η C 2 η C. (3) It is then realized that ost of the observed efficiencies fall within these bounds [3, 4]. Η CA o FIG. 1. Data on the observed efficiencies (η o ) of power plants, plotted against the CA-value, η CA = 1 T /T + for respective plants as given in Table I. A point lying on the straight line has the observed efficiency equal to CA-value. So for points above the line, η o is below the CA-value, while the converse is true for the points below this line. Naturally, the question of the actual working constraints and the real optiization targets for each plant, is also relevant. Still, the effectiveness of these siple odels, in reproducing the gross features of diverse plants, cannot be denied. Apart fro finite-tie odels, the irreversibilities reducing the efficiency to a lower-than- Carnot value, ay also be treated within a quasi-static work extraction odels with finite source/sink [13 16]. Soe of these studies also derive the optial value of an interediate

4 4 Industrial Plant T + T η o η CA T ( ) T (+) T (av) G(T +, T ) Alaraz II (Nuclear, Spain) [2] Calder Hall (Nuclear, UK) [2] CANDU (Nuclear, Canada) [1] Cofrentes (Nuclear, Spain) [2] Doel 4 (Nuclear, Belgiu) [2] Heysha (Nuclear, UK) [2] Larderello (Geotheral, Italy) [1] Sizewell B (Nuclear, UK) [2] West Thurrock (Coal, UK) [1] Pressurized water nuclear reactor [12] Boiling water nuclear reactor [12] Fast neutron nuclear reactor [12] (Stea/Mercury, US) [2] (Stea, UK) [2] (Gas Turbine, Switzerland) [2] (Gas Turbine, France) [2] TABLE I. Observed efficiencies (η o ) of power plants working between teperatures T + and T, copared with η CA = 1 T /T +. The effective teperatures are defined as: T (+) = T + (1 η 0 ), T ( ) = T /(1 η 0 ), T (av) = (T (+) + T ( ) )/2, and G(T +, T ) = T + T. teperature which is the geoetric-ean value, and consequently the efficiency at axiu work equals the CA-value. However, here again, the siplifying assuption of a constant heat capacity, say, of the source or the working ediu, heavily deterines the conclusion. More recently, it was observed [17] that within linear response theory, the bounds such as Eq. (3), also follow within a quasi-static odel of work extraction fro finite-sized heat source and sink. Ref. [17] akes no specific assuption about the nature of the heat source/sink. The present study is otivated by this analysis, and ais to study the perforance of real plants within a quasi-static approach. More precisely, we address an inverse question. Instead of finding the optial interediate teperature and then fro it, the efficiency of the process, we use the inforation on the highest and the lowest tepera-

5 5 tures, along with the value of observed efficiency, to infer the interediate teperature at axiu work. The correspondence between this teperature and efficiency is as follows. If the forer is exactly equal to the geoetric-ean value, then the efficiency is equal to CA-value. The odel of work extraction is based on one syste as a finite source, or sink, and the other as an infinite reservoir, which allows for two alternate scenarios: i) when the source is finite and ii) when the sink is finite. The liited inforation on the actual situation being realized out of these two possibilities, otivates to do an inference analysis. Thus we estiate an expected value of the interediate teperature scale. Interestingly, this value is found to be very close to the geoetric ean of T + and T. We also present a quantitative arguent for this proxiity to the geoetric-ean value. Thus our analysis indicates the role of geoetric ean fro a novel angle which has two distinctive features: first, it is based on axiu work approach, and second, we do not per se assue a specfic nature of the finite source or sink. Our starting point is the reversible cycle operating between two infinite heat reservoirs, for which the efficiency is the Carnot value η C. As a first step arking a deviation fro reversibility, we consider one reservoir as finite, while the other reservoir reains very large copared to the forer, or practically infinite. Now, we first assue that syste A acts as a finite heat sink at teperature T, relative to a very large heat source at teperature T +. The two are coupled by an ideal engine which delivers work to a reversible work source, by running in infinitesial heat cycles, that successively increase the teperature of A, till A coes in equilibriu with the source, see Fig.1 (i). At this point, the extracted work is axiu. Suppose that in this process, the syste A oves fro an equilibriu state characterized by energy U and entropy S to an equilibriu state with the corresponding values of U + and S +. The heat reoved fro the source is Q + = T + (S + S ). The heat rejected to the finite sink is q + = U + U. The work extracted is given by W + = Q + q +, or W + = T + (S + S ) (U + U ). (4) Then the efficiency at axiu work, η + = W + /Q +, is evaluated to be: where we have defined η + = 1 T (+), (5) T + T (+) = U + U S + S, (6)

6 6 (i) A W (ii) B W FIG. 2. Scheatic of the engine between a finite syste and a heat reservoir, for a given pair of initial teperatures (T +, T ): (i) Syste A as a finite sink at T and an infinite source at T +, coupled with a reversible work source. Work extraction W +, Eq. (4), is copleted when the teperature of A becoes T +. (ii) Syste B as a finite source at T + and an infinite sink at T. Maxiu extracted work is W, Eq. (8), when the teperature of B becoes T. a quantity characteristic of syste A. Further T (+) ay be regarded as an effective teperature of an infinite reservoir [18], such that the present situation of an infinite source and a finite sink, at teperatures T + and T respectively, is equivalent to axiu work extraction in a reversible cycle fro two infinite reservoirs at teperatures T + and T (+) (< T + ). In the latter case, the extracted work per cycle is: W + = (T + T (+) )(S + S ), with (Carnot) efficiency 1 T (+) /T +, which is the sae as Eqs. (4) and (5). Now, we assue that the inforation on the nature of syste A is not available, or, in other words, T (+) is not given. But if the quasi-static odel with an infinite source and a finite sink, is a good odel for the observed perforance of an industrial plant, then we ay infer the relevant value of the effective teperature T (+), when the theoretical efficiency (η + ) for the odel is set equal to the observed value η o. This iplies that we write T (+) = T + (1 η o ). (7) Such values for soe industrial plants, based on the observed values of their efficiencies, are tabulated in Table I, and also depicted in Fig. 3, in coparison to geoetric ean G(T +, T ) = T + T. The latter value is chosen siply because when T (+) = G, the observed

7 G, T FIG. 3. The effective teperature T (+) plotted against the geoetric ean G(T +, T ) = T + T, the two quantities being equal along the straight line. efficiency is equal to the CA-value. Thus the spread of observed efficiency values around the CA-value in Fig. 1, is translated here, to a spread in the effective teperatures around the geoetric-ean values. More precisely, T (+) G indicates η o η CA, as is evident fro Table I. It is to be noted that inferring the effective teperature fro the observed efficiency does not deterine the nature of the syste A or the for of fundaental relation U(S). The latter inforation is not at our disposal within the assuptions of the odel. 550 G, T FIG. 4. The effective teperature T ( ) plotted against the geoetric ean G(T +, T ) = T + T, the two quantities being equal along the straight line. However, the scenario of work extraction fro a finite syste coupled to an infinite reservoir, via a reversible process, is consistent with an alternate picture too. This involves

8 8 that, for the sae initial teperatures T + and T, we consider a finite source (B) at T +, coupled to an infinite sink at T, see Fig. 2 (ii). Here again, we can extract axial work by utilizing the teperature gradient between B and the reservoir, till B coes to be at teperature T [19]. Assuing that the syste goes fro soe equilibriu state (U +, S +) to another one (U, S ), the heat reoved fro the source is Q = U + U and the aount rejected to sink is q = T (S + S ). So the extracted work is W = Q q, or W = (U + U ) T (S + S ). (8) This is called exergy in the engineering literature [20]. The efficiency η = W /Q is given by where η = 1 T T ( ), (9) T ( ) = U + U. (10) S + S It is clear fro the expressions for W and η, that an equivalence exists between the above odel and that of work extraction in a reversible cycle fro two infinite reservoirs at T and T ( ) (> T ). Now again, to apply the above case to odel an industrial plant, we can equate the observed efficiency to the theoretical efficiency, as η o = η and infer the corresponding effective teperature T ( ) It is clear that T + fro Eq. (9), as T ( ) = T 1 η o. (11) > T ( ) > T. The calculated values of T ( ), based on the observed efficiencies of the plants, are also tabulated in Table I, and shown graphically in Fig. 4 in coparison to G(T +, T ). The following reark sees appropriate here. An analogy ay be ade between the above odel and the earlier irreversible odels such as proposed by Chabadal [2, 6], in which an interediate teperature T w (of the war exit strea) enters the analysis and the theoretical efficiency of the odel is 1 T /T w (copare with Eq. (9)). For an optial value T w = T + T, the power output becoes optial with the corresponding efficiency equal to CA-value. The crucial differences, fro our odel, are twofold: a) we study quasistatic axiu work extraction process b) the nature of syste B, the finite source, is not specified.

9 9 Now as observed in Figs. 3 and 4, the values of effective teperatures see to be distributed, apparently in a rando fashion, about the geoetric-ean value. However, it is rearkable to note that for a given plant, the calculated values of T (+) and T ( ) are such that, to a high accuracy, they are equidistant fro the geoetric ean G. More precisely, if we define an average value of teperature T (av), as T (+) T (av) = T (av) value T (av) T ( ), then this is very close to the G value for that situation. In other words, we define an average scale of teperature as the arithetic ean of the two inferred teperatures T (±), and so given by T (av) = 1 [ T + (1 η o ) + T ]. (12) 2 1 η o Then the above average value is found to be closely approxiated by G(T +, T ) for ost observed cases of industrial plants, see Fig. 5, as well as Table I. Now, we turn to a ore quantitative characterization of the above observation. It is G, T av FIG. 5. The average effective teperature T (av) plotted against the geoetric ean G(T +, T ) = T+ T, the two quantities being equal along the straight line. T (av) values calculated fro the observed data on ost of the power plants, is rearkably close to G(T +, T ). easy to see that T (av) takes a iniu value of T + T, w.r.t to η o (dt (av) /dη o = 0, and d 2 T (av) /dηo 2 > 0), at η o = 1 T /T +. This iplies that any possible value of T (av) is equal to or greater than T + T, and so will lie on or below the straight line plotted in Fig. 5. The second aspect is related to the observation ade earlier that the deviations in values of effective teperatures T (±) fro the corresponding G-values, reflect the fact how the observed efficiency deviates fro the CA value. However, the deviations fro G-values, are

10 10 suppressed considerably in case of T (av). This ay be argued by considering sall deviations (ɛ) in efficiency fro the CA value, and expanding T (av) in powers of ɛ. Let η o = η CA + ɛ. Then we see that up to second order: T (av) = T + T + 1 T+ 3 ɛ 2 + O(ɛ 3 ). (13) 2 T Thus, the first non-zero correction fro the geoetric-ean value is of second order in ɛ, while it is of the first order for η o, or T (±). Clearly, the agnitude of fluctuations is suppressed in the case of T (av). Finally, we address the eaning of the average effective teperature. If we again consider the two extree situations envisaged in Fig. 2, then they are utually exclusive, or one ay say, they are counterfactual. The average teperature is not necessarily seen in an actual realization, except for the special case T (+) = T ( ) = T + T. In this sense, the eaning one can attach to the definition of T (av), can be given best in the language of inductive inference [21, 22]. In latter ters, the average teperature represents an expected scale of the effective teperature, in view of our inability to choose between the two alternatives (i) and (ii), where each represent our hypothesis for the odel of work extraction applicable to the plant. In inference, when one of, say, two utually exclusive hypotheses cannot be given a preference over the other, then we ust assign equal weights to the inferences derived fro each [22, 23]. The effective teperatures T (±) are our inferred values fro the given data on the observed efficiency, and, so T (av), defined with equal weights assigned to both the inferences, represents the expected value of the effective teperature, coensurate with the inforation at our disposal. Any deviation fro equal weights would iply that we have soe extra piece of inforation about the process, and thus would be inconsistent [24]. Concluding, we have proposed to odel the observed efficiencies, using the quasi-static odels of work extraction where one reservoir (hot or cold) is finite while the other is practically infinite. This brings an extra scale of interediate teperature into the analysis. We note two consistent odels depending on which reservoir is taken as finite, and consequently, we get two possible values of the interediate teperature. The liited inforation on the working conditions does not allow to prefer one odel hypothesis over the other and so an equally-weighted average represents the expected scale of interediate teperature. For ost of the power plants, this inferred value is found to be quite close to the geoetric ean of the highest and the lowest teperatures. Thus we rediscover the significance of

11 11 the geoetric-ean teperature, which was ephasized in earlier irreversible odels of plants [5 7, 25], but there it was often based on siplifying assuptions such as constant heat capacity for the hot strea, and Newton s law for heat transfer. Such assuptions are not iportant in our analysis and we do not ake use of a particular nature of the finite reservoir. In the present approach, if the deviations in the observed efficiency fro the CA-value are sall, then the geoetric-ean teperature appears as a rational estiate for the interediate teperature. An inference based approach has been used earlier to study odels of therodynaic achines with liited inforation. The eergence of CA efficiency fro inference has been noted in quantu heat engine [26] and esoscopic odels like Feynan s ratchet [27]. For classical odels of work extraction fro two finite reservoirs, the results for efficiency at axiu work are reproduced through inference, beyond linear response [28]. Further, reversible odels with liited inforation have been related to irreversible odels through inference based reasoning [29]. In the present context, it is rearkable how siple, but general considerations can lead to estiates close to the geoetric-ean value for the interediate teperature. It is indeed curious to ask, how the inference approach leads to the sae conclusions as the siple irreversible odels with constant heat capacity. Does it indicate a connection between the siplicity of objective odelling and assuing inial inforation in inference? More research is needed for a deeper understanding of the connection between the use of liited inforation and therodynaic odelling. The purpose of inductive inference is not just to provide an acadeic and philosophical point of view which ay soeties rival the ore objective odelling we usually resort to in science, but also seek how ethods of inference can give insight into the actual state of affairs, while incorporating the prior inforation, norally not considered in therodynaic odels. In this context, we note that although ost of the data on plants yield an expected value of the interediate teperature close to geoetric-ean value, still there are a couple of not-so-insignificant deviations in Fig. 5, near the top right corner. A valid question is, why do these exaples differ fro the rest of the cases? Does it indicate other easures of optiization being used in the actual operation? As far as the available data is concerned, we note these plants operate under higher teperature gradients than ost other plants. Thus, one ay argue that our odel and the assuptions ay not be good approxiations for large teperature differences. In any case, once we have noticed deviations fro the

12 12 expected behavior in these plants, further studies on the actual working conditions ay be undertaken, which ay then yield inforation to generalize the inference analysis, or ay help to iprove the perforance, at par with other plants. [1] F. L. Curzon and B. Ahlborn, A. J. Phys. 43, 22 (1975). [2] A. Bejan, Advanced Engineering Therodynaics (Wiley, New York, 1997) p [3] M. Esposito, R. Kawai, K. Lindenberg, and C. Van den Broeck, Phys. Rev. Lett. 105, (2010). [4] G. B. Moreau, M. and L. S. Schulan, Phys. Rev. B 85, (2012). [5] H. B. Reitlinger, Sur l utilisation de la chaleur dans les achines á feu ( On the use of heat in stea engines, in French). (Vaillant-Caranne, Liége, Belgiu, 2008). [6] P. Chabadal, Arand Colin, Paris, France 4, 1 (1957). [7] I. Novikov, J. Nucl. Energy II 7, 125 (1958). [8] A. Vaudrey, F. Lanzetta, and M. Feidt, J. Noneq. Ther. 39, 199 (2014). [9] G. Lebon, D. Jou, and J. Casas-Vázquez, Understanding Non-equilibriu Therodynaics: Foundations, Applications, Frontiers (Springer-Verlag, Berlin Heidelberg, 2008). [10] B. Andresen, Ange. Che. 50, 2690 (2011). [11] A. Bejan, J. Appl. Phys. 79, 1191 (1996). [12] M. Borlein, Kerntechnik (Vogel Buchverlag, Wurzburg, Gerany, 2009). [13] W. Thoson, Phil. Mag. 5, 102 (1853). [14] M. J. Ondrechen, B. Andresen, M. Mozurchewich, and R. S. Berry, A. J. Phys. 49, 681 (1981). [15] H. S. Leff, A. J. Phys. 55, 602 (1986). [16] B. H. Lavenda, A. J. Phys. 75, 169 (2007). [17] R. S. Johal and R. Rai, EPL (Europhys. Lett.) 113, (2016). [18] R. S. Johal, arxiv: [19] Y. Izuida and K. Okuda, Phys. Rev. Lett. 112, (2014). [20] M. J. Moran, H. N. Shapiro, D. D. Boettner, and M. B. Bailey, Fundaentals of Engineering Therodynaics, 7th ed. (Wiley, New York, 2010) chap. 7. [21] H. Jeffreys, Theory of Probability, 2nd ed. (Clarendon Press, Oxford, 1948).

13 13 [22] E. T. Jaynes, Probability Theory: The Logic of Science (Cabridge University Press, Cabridge, 2003). [23] P. S. Laplace, Meoir on the Probabilities of the Causes and Events, Stat. Sc. 1, (translated by S. M. Stigler, 1986). [24] E. T. Jaynes, The evolution of Carnot s principle, edited by G. J. Erickson and C. R. Sith, in Maxiu-Entropy and Bayesian Methods in Science and Engineering (Kluwer, Dordrecht, 1988). [25] J. Yvon, The Saclay Reactor: Two years experience in heat transfer by eans of copressed gas as heat transfer agent, in Proceedings of the International Conference on Peaceful Uses of Atoic Energy (Geneva, Switzerland, 1955). [26] R. S. Johal, Phys. Rev. E 82, (2010). [27] G. Thoas and R. S. Johal, J. Phys. A: Math. and Theor. 48, (2015). [28] P. Aneja and R. S. Johal, J. Phys. A: Math. and Theor. 46, (2013). [29] R. S. Johal, R. Rai, and G. Mahler, Found. Phys. 45, 158 (2015).