SciELO - Scientific Electronic Library Online

 
vol.63 issue6Glucose in aqueous solution thermal characterization by photopyroelectric techniquesEnergy on a timelike particle in dynamical and electrodynamical force fields in De-Sitter space author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.63 n.6 México Nov./Dec. 2017

 

Research

The universal optimal relations of the allocation and effectiveness of the heat exchangers for power plants with n Carnot-like cycles

G. Aragón-González1 

A. León-Galicia1 

1PDPA. UAM- Azcapotzalco. Av. San Pablo # 180. Col. Reynosa. Azcapotzalco, Ciudad de México, México 02200. Tel. and FAX: (55) 5318-9057. *e-mail: gag@correo.azc.uam.mx


Abstract

A model of irreversible Carnot-like power plant with n Carnot-like cycles is optimized. The irreversibilities of each cycle are: finite rate heat transfer between the working fluid and the external heat sources, internal dissipation of the working fluid, and heat leak between reservoirs; is extended to two or more of this combined model. Applying the Bellman’ Principle, we find the optimal recurrence relations for the allocation of the heat exchangers for thes power plants. The optimal allocation or effectiveness is determined by two design rules, applied alternatively: internal thermal conductance fixed or areas fixed. The optimal recurrence relations obtained for this combined model are invariant to the power and efficiency and to the heat transfer law.

Keywords: Irreversibilities; Carnot; allocation; conductance; effectiveness; area; optimal

PACS: 01.40G; 05.70.-a; 64.70.F

Resumen

Se optimiza un modelo de planta de potencia irreversible con n ciclos tipo Carnot. Las irreversibilidades de cada uno de los ciclos son: transferencia finita de calor entre el fluido de trabajo y los depósitos de temperatura, disipación del fluido de trabajo y fuga de calor entre ambos depósitos; es extendido a dos o más de este ciclo combinado. Aplicando el Principio de Bellman, encontramos las relaciones de recurrencia óptimas para la dimensión y efectividad de los intercambiadores de calor para modelo combinado. La óptima dimensión y efectividad se determina mediante dos reglas de diseño, aplicadas alternativamente: conductancias térmicas restringidas ó área total constante. Las relaciones óptimas para la dimensión ó efectividad obtenidas para este modelo combinado son invariantes a la potencia y eficiencia y a la ley de transferencia de calor.

Descriptores: Irreversibilidades; Carnot; dimensi´on; efectividad; conductancia; área; óptima

PACS: 01.40G; 05.70.-a; 64.70.F

1. Introduction

Recently in1, a methodology of optimization was applied to an irreversible Carnot-like power plant, where the characteristic parameters were: the allocation or effectiveness of the heat exchangers of this plant. Although this methodology is applicable to any Carnot-like power plant, a standard irreversible Carnot-like cycle was chosen because of its simplicity to account for the main irreversibilities that usually arise in real heat engines2: “finite rate heat transfer between the working fluid and the external heat sources, internal dissipation of the working fluid, and heat leak rate between reservoirs”. The above standard irreversible Carnot-like cycle has been studied at length for many objective functions, different transfer heat laws and several characteristic parameters (see2-24 for more details).

The maximum power and efficiency have been obtained in3-6. In general, these optimizations were performed with respect to only one characteristic parameter: internal isentropic temperature ratio (x) including the time too18. However, Bejan and Lewins15,20 have considered the optimization with respect to other parameters (see also the reviews of17-18 for more details): “the allocation, cost and effectiveness of the heat exchangers of the hot and cold sides”. On the other hand, effects of heat transfer laws or when a property is independent of the heat transfer law for this Carnot-like cyclic model, have been discussed in several works,18-23, and so on.

On the other hand, n-stage combined Carnot cycle has been presented, by the law of heat conduction, in24-26 optimizing the specific power and efficiency for the characteristic parameters: isentropic temperature ratios , and effectiveness (other combined cycles can be found in27-36, coupled heat devices in37-38; see also the references therein included).

In this paper we extend the results of allocation or effectiveness for one cycle to two or more combined cycles of the same irreversible Carnot-like power plant using iteratively the Bellman’ Principle39 which has been successfully applied in Refs. 13 and 14. We found optimal recurrence relations remarkable for the irreversible power plant with n Carnot-like cycles for two constraints (design rules): constrained internal thermal conductance or fixed total area of the heat exchangers from hot and cold side.

This paper is organized as follows. In Sec. 2 we present the irreversible power plant model and the functional relation between power and efficiency which extends to the corollary presented in1 to the irreversible power plant. Section 3 presents the optimization of the allocation and effectiveness, by two design rules, of heat exchangers for the Carnot-like 𝑛-cycles of the power plant. Section 4 is devoted Conclusions.

2. Power plant with n Carnot-like cycles

The power plant with n Carnot-like cycles is shown in Fig. 1. Each cycle satisfies the conditions expressed in2,4: leak heat rate Q and finite heat transfer rates Qi and internal dissipations of the working fluid expressed by constants Ii(i=1n) such that Ii=(ΔSi+1/ΔSi)1;i=1n 4 which make the Claussius inequality to become equality (for more details see1,24). Each cycle of the power plant consists of two isothermal and two adiabatic processes. Denoting, for each cycle, the temperatures of the working fluid in the hot and cold isothermal processes as T2i-1 and T2i(i=1n), respectively, and the end temperatures as TH and TL.

Figure 1 A power plant with n Carnot-like cycles, with heat leak rate, finite heat transfer rates, and internal dissipations of the working fluid in each cycle. 

Following to1,25, the thermal efficiency of the power plant is given by (this functional form for only a cycle has appeared in Refs. 1, 2, and 17):

η=Pf(x1,.....,xn)P+Q (1)

where i=T2i/T2i-1 corresponds to the isentropic temperature ratio for each cycle i(i=1,....,n), P is the power of the plant, Q is the leak heat rate and the function

f(x1,.....,xn)=11-(i=1nIixi) (2)

The Eq. (1) is obtained of the following way: according to the second Law of the Thermodynamics, for each i(=1,2,...,n), we have (cfr. With25),

Qi+1=IixiQi (3)

Thus,

Qn+1=InxnQn=....=IixiQi=......=I1x1Q1 (4)

Now, using the Eqs. (3) and (4), the efficiency is given by,

(5)

From the Eq. (5), Q1 is equal to:

Q1=P1-(i=1nIixi) (6)

And from the Eq. (6), (1) is obtained.

The Eq. (1) extends to the Eq. (6) of1. Therefore, we can seek, as in1, invariant optimal relations, by following the spirit of Carnot’s work, for other characteristic parameters, which can be independent from the heat transfer law, of the irreversible power plant.

Also, the Eq. (1) say us that it is enough analyzes only for power and efficiency (for other operation regimes, e.g. algebraic combination of power and/or efficiency that have thermodynamic meaning and satisfy imposed power conditions: see Eq. (7) below and for more details see1). For instance, the functional expression of the efficiency, power and the heat transfer, without require of algebraic explicit expressions, have been used in2, for two constraints (design rules): constrained internal thermal conductance or fixed total area of the heat exchangers from hot and cold side. The isentropic temperature ratio do not has this condition1,2,3,13,24-26.

Now, if we take zxi(i=1,2....n) as other independent variable the Corollary given at1 can be paraphrased as:

“The power 𝑃 achieves a maximum value in 𝑧 𝑚𝑝 if and only if the efficiency 𝜂 achieves a maximum value in the value: zmp=zme

Thus, in the optimization of power plant, with respect to 𝑧 it is enough to find the maximum power by,

Pz|zmp=zme=0and2Pz2|zmp=zme<0 (7)

The optimization performed, with respect to will be: a property independent from the heat transfer law. This a remarkable conclusion of the Eq. (1) is that it can find the maximum for one and only one the power (see the criterion of2). For example the power, for power plant model, by

Pz|zmp=zme=0and2Pz2|zmp=zme<0 (8)

We show the above as follows: Let zxi(i=1,2....n) and the vector x=(x1,....,xn) and zmp the point in which the power achieves a maximum value, then,

Pz|zmp=zmp=0and2Pz2|zmp<0.

The power and the efficiency satisfy the functional relationship given by the Eq. (1). Deriving (1) with respect to:

ηz=Q(Pz)f(x)P+Q2 (9)

since, we can suppose that Q does not depend of the variable z.

Therefore,

ηz|zme=0=Pz|zmp;

where zme is the point in which the efficiency η achieves a maximum value which implies that their roots are the same zmp=zme. Also, it is easily seen that for zmp=zme, the efficiency reach a maximum since

2ηz2|zme=zme=Q2Pz2|zme=zmef(x)P+Q2<0 (10)

In all the above calculations none transfer heat law has been used.

Therefore, it is enough to have one, and only one, algebraic expression of the power and any heat transfer law. We choose for our convenience the power and the conduction heat transfer Law since they are algebraically the simplest. In the optimization of the power of the plant for two constraints (design rules): constrained internal thermal conductance (z=ϕ) or fixed total area (z=ψ) of the heat exchangers from hot and cold side of the heat exchangers from hot and cold side.

3. Universal optimal relations for the allocation and effectiveness of the heat exchangers of the power plant

The optimization of power and efficiency with respect to is well known for combined cycles13,24-26. Henceforth, x will be fixed and we will assume that the law of heat transfer rate can be any law, included also the heat leak rate. Essentially, we will treat the following two design rules: internal thermal conductance fixed; or areas fixed for heat exchangers; which will be applied alternatively.

The first design rule is that the internal conductance of the Carnot-like cycle is constrained to: i=1nαi=Γ; where Γ is a constant, which is applied to the allocation of the heat exchangers from hot and cold side with the same overall heat transfer coefficient U by unit of area A in both ends of the cycle i and αi,αi+1(I=1,2....n-1) are the thermal conductance correspondent to the finite heat transfers of the hot/cold sides for this cycle. Thus,

i=1nUAi=Γ (11)

where Ai,Ai+1 are heat transfer areas on hot/cold sides of the cycle i(i = 1,2…n - 1).

Alternatively, the second design rule is that the total area is constrained by: i=1nAi=A; where Ai, Ai+1 are heat transfer areas on cold and hot side for the cycle i. Now, the total area (A) is fixed, but when distributed it has distinct overall heat transfer coefficients and hence different effectiveness on each one of the hot and cold sides. As αi=UiAi (see 1 and 12), then,

A=i=1nAi=i=1nαiUi (12)

where Ui, Ui+1 are the overall heat transfer coefficients on cold/hot sides of the cycle i(i=1,2.....n). From the result of the Sec. 2, it is enough to found “The maximum power of each Carnot-like cycle i(i=1,2.....n) for the two design rules (11) and (12)”

Indeed, the dimensionless power output (Eq. (2.5) for one cycle in3) for each cycle i is given by:

pi=piαiT2i-1=h(xi)1α1+Iiα1+1 (13)

here h (xi) is a function which depends on the heat transfer law since xi=T2i/T2i-1 is the isentropic temperature ratio of the cycle i ([1,2,3,13,24-26]); and i, αi+1 i(i=1,2.....n-1) are the thermal conductance correspondent to the finite heat transfers of the cold/hot sides for this same cycle.

Next, applying the Bellman’ Principle39, in form iterative, we can obtain the optimal relations for the allocation of the heat exchangers of the power plant (Fig. 1). This principle will be applied as it is indicated in20: “to state that every part of an optimum path is optimal”.

3.1 Constrained internal thermal conductance

The 𝑛 thermal conductance can be written as: αi=UA1; i = 1…n + 1 where U is overall heat transfer coefficient and 𝐴 𝑖 ; i=1n+1 are the available areas for heat transfer. Thus, for the first optimization we can take for this first design rule:

α1+α2=Γ1

where Γ1=Γ-i=3nαi is supposed to be a constant. Equivalently,

α1+α2Γ1=1

Fixing the temperature Tj(j=3n+1) and applying only for the first cycle:

a1α1+α2Γ1=1;a1=1

In parameterizing,

ϕ=α1Γ1;1-a1ϕ=α2Γ1

According to the Eq. (13), in optimizing

1a1ϕ+I11-a1ϕ

with respect to ϕ; we obtain:

ϕ1=11+I1

since a1=1. Solving,

1-a1ϕϕ1=I1=α2α1;α2=I1α1.

Continuing of this way we arrive to the following optimal recurrence relation:

αnIn-1α1,n=2,3...... (14)

And

an=j=0n-1Ij;I0=1.

From the Eqs. (11) and (14) we have: the areas de transference of heat decreases from TH to TL:

A1A2........An-1An (15)

and the equality is fulfilled if there is not internal dissipation in the n cycles.

3.2 Constrained areas of heat exchangers

For simplicity, we suppose Ii=1(i=1....n), without internal irreversibilities. Applying the design second design rule (Eq. (12)),

A1+A2=A1

where A1=A-i=3nAi is supposed a constant and fixing the temperature Tj(j=3......n+1). How,

A1=α1U1;A2=α2U2

Then

A1=a1α1+u1α2

where u1=U1/U2; a1=1. In parameterizing,

ψ=α1A1;1-a1ψ=α2A1

From according to the Eq. (13):

1a1ψ+1(1-a1ψ)u1,

the first optimization with respect to ψ, gives:

ψ1=u11+a1u1

Then,

1-a1ψ1ψ1=1+u1(a1-1)u1.

Solving,

A2=A1a1u1

and, the results of1 are recovered.

Continuing of this way arrive us to the following optimal recurrence relation:

An=1+un-1(an-1-1)1+an-1un-1A1 (16)

With

un-1=U1/Un;an-1=a2+an-2,n>2;a2=1+1a1u1;a1=1.

3.3 Some special cases

To illustrate the optimal relations (14 and 16) we can choose a power plant with 2 or 3 like-like-Carnot cycles:

a). Case constrained internal thermal conductance

In this case, from the Eq. (14), we obtain the conductance: α2=I1α1; α3=I2α1 and α4=I3α1 for only three cycles. If Ij=1; j=1,2,3 (without internal irreversibilities). Then, the optimal conductance without internal irreversibilities are obtained: α2=α3=α4=α1=(1/4)Γ4. For two cycles: α3=α2=α1=(1/3)Γ3 only and so on. This last result was reported in13 using only the heat conduction law. Thus, the expressions obtained from (14) are extended to any heat transfer law.

b). Constrained areas:

In this case, from the Eq. (16), we obtain the areas, for two cycles,

A1=u11+u1A;A2=A1+u1;A3=u1+u1u2(u1+u2u1+u2)(1+u1).

And for three cycles,

A2=1u1A1;A3=u1+u2u1+u2u1+u2A1A4=u1+u2+u2u1+u3u1+u3u2u1+u2+u1+2u3u1+2u3u2+u3u2u1A1

2 where uj-1=U1/Uj; Ij=1; j=2,3,4; and so on. The expressions obtained by (16) are extended for any heat transfer law.

4. Conclusions

We have found and determined the optimal allocation of heat exchangers of a power plant with n Carnot-like cycles. We must note that the optimal recurrence relations (Eqs. (14 and 16)) found for this model are valid for any heat transfer law and have not been reported in the literature previously; specially the Eq. (16). Thus, following the spirit of Carnot’s work, we have found invariant optimal relations for power and efficiency maximum which are valid for any heat transfer law. Moreover, these relations can be satisfied for other operation regimes, e.g. algebraic combination of power and/or efficiency that have thermodynamic meaning and satisfy the conditions imposed to the power (Eq. (7)) which can be carried out as in1. Nevertheless, the optimal isentropic temperature ratios depend on the heat transfer law and the operation regime of the engine as is discussed in1,2,3,13,24-26.

Finally, the above equations can be extended to n Carnot-like cycles presented here including other heat transfer rate laws by the substitution of either of (14) or (16) in the objective function (algebraic combination of power and/or efficiency) and optimize only for the isentropic temperature ratio. However, the latter requires a comprehensive study of the implications for the power plant considered. We will study such implications as future work.

REFERENCES

1. G. Aragón-González, A. Canales-Palma, A. Galicia-León, and J.M. Rivera-Camacho, J. Energy Inst. 85 (2012) 201-208. [ Links ]

2. G. Aragón-González, A. Canales-Palma, A. Galicia-León, and J.M. Rivera-Camacho, J. Phys. A: Math. Theor. 42 (2009) 425205. [ Links ]

3. G. Aragón-González, A. Canales-Palma, A. Galicia-León and I.R. Morales-Gómez, Braz. J. of Phys. 38 (2008) 543-550. [ Links ]

4. J. Chen, J. Phys. D: Appl. Phys. 27 (1994) 1144-1149. [ Links ]

5. Z. Yan and L. Chen, J. Phys. A: Math. Gen. 28 (1995) 6167-6175. [ Links ]

6. G. Aragón-González, A. Canales-Palma, A. Galicia-León and M. Musharrafie-Martínez. J. Phys. D: Appl. Phys. 36 (2003) 280-87. [ Links ]

7. N. Cegonho and T. Domingos, J. Non-Equilib. Thermodyn. 30 (2005) 113-127. [ Links ]

8. L.G. Chen, J. Zhu, F. Sun, C.A. Wu, Appl. Energ. (2004) 78 305-13. [ Links ]

9. L.G. Chen, J. Li, F.R. Sun, Appl. Energ. 85 (2008) 52-60. [ Links ]

10. J.M. Gordon, J. Appl. Phys. 69 (1991) 1-7. [ Links ]

11. A. Bejan, Entropy Generation through Heat and Fluid Flow (Wiley, New York 1982). [ Links ]

12. A. Bejan, Advanced Engineering Thermodynamic (Wiley, New York 1988). [ Links ]

13. A. Bejan, Int. J. Heat Mass Transfer. (1995) 38 433-44. [ Links ]

14. A. Bejan, Entropy Generation Minimization (CRC Press, Boca Raton, FL 1996). [ Links ]

15. F.L. Curzon and B. Ahlborn, Am. J. Phys. 43 (1975) 22-24. [ Links ]

16. M.H. Rubin, Phys. Rev. A 15 (1979) 2094-2102. [ Links ]

17. K.H. Hoffman, J.M. Burzler and S. Shuberth, J. Non-Equilib. Thermodyn. 22 (1997) 311-355. [ Links ]

18. A. Durmayaz, O.S. Sogut, B. Sahin and H. Yavuz, Progr. Energ. and Combus. Sci. (2004) 30 175-217. [ Links ]

19. G. Aragón-González, A. Canales-Palma, A. Galicia-León and J.R. Morales-Gómez. Rev. Mex. Fis. (2006) 52 309-314. [ Links ]

20. J.D. Lewins, Int. J. Mech. Eng. Edu. 28 (2000) 41-46. [ Links ]

21. A. De Vos A. Am. J. Phys. 53 (1985) 570-573. [ Links ]

22. A. Calvo, A. Medina, J.M. Roco, J.A. White and S. Velasco, Phys. Rev. E 63 (2001) 0371021-5. [ Links ]

23. G Aragón-González, A Canales-Palma and A Galicia-León. J. Phys. D: Appl. Phys. 33 (2000) 1403-10. [ Links ]

24. J.D. Lewins, Int. J. Mech. Eng. Edu. (1999) 27 91-101. [ Links ]

25. J. Chen, J. Phys. A: Math. Theor. 3 (1998) 3383-3394. [ Links ]

26. S. Bandyopadhyay, N.C. Bera and S. Bhattacharyya, Energ. Convers. Manag. 42 (2001) 359-371. [ Links ]

27. L.G. Chen, F.R. Sun and W.Z. Chen, Energy The International Journal 20 (1995) 1049-1053. [ Links ]

28. L.G. Chen, C. Wu and F.R. Sun, Applied Thermal Engineering 17 (1997) 639-645. [ Links ]

29. L.G. Chen, F.R. Sun, C. Wu and J.Z. Gong, Applied Thermal Engineering 17 (1997) 629-637. [ Links ]

30. L.G Chen, H. Duan, F.R. Sun and C. Wu, Journal of Non-Equilibrium Thermodynamics 24 (1999) 280-290. [ Links ]

31. L.G. Chen, Y.H. Bi, F.R. Sun and C. Wu, International Journal of Thermal Science 38 (1999) 712-718. [ Links ]

32. L.G. Chen, S. Cao, C. Wu and F.R. Sun, International Journal of Power and Energy Systems (1999) 19 103-106. [ Links ]

33. L.G. Chen, F.R. Sun, C. Wu, N. Ni, International Journal of Energy, Environment and Economics, 9 (1999) 35-49. [ Links ]

34. L.G. Chen, C. Wu and F.R. Sun, Energy Conversion and Management, 39 (1998) 113-116. [ Links ]

35. K. Ma, L.G. Chen, and F.R. Sun , Academy Proceedings in Engineering Science, 34 (2009) 851-864. [ Links ]

36. L.G. Chen, D.L. Ni, Z.L. Zhang and F.R. Sun, Applied Thermal Engineering, 102 (2016) 447-453. [ Links ]

37. B. Jiménez de Cisneros and A Calvo Hernández. Phys. Rev. E 77 (2007) 041127 [ Links ]

38. B. Jiménez de Cisneros and A. Calvo Hernández, Phys. Rev. Letters 98 (2008) 130602. [ Links ]

39. R. Bellman, Applied dynamic Programming (Princeton 1962). [ Links ]

40. D. Gunderson, Handbook of Mathematical Induction: Theory and Applications CRC Press series on discrete mathematics and its applications (2011). [ Links ]

Appendix

For the second optimization of the first design rule, the constraint is now:

a2α1+α3Γ2=1

with a2=1+I1; Γ2=Γ-i=4nαi a constant and fixing the temperature Tj(j=4....n+1) and applying only for the second cycle. In parameterizing,

ϕ=α1Γ2;1-a2ϕ=α3Γ2

Similarly, in optimizing

1a2ϕ+I21-a2ϕ

with respect to ϕ; we obtain:

ϕ2=1a2+I2

Solving

1-a2ϕ2ϕ2=I2=α3α1

so,

α3=I2α1

and a3=1+I1+I2. Continuing of this form arrive us to the Eq. (14). We can show the relation (14) by Mathematical Induction

38. Indeed for n=2 is true. We suppose that for n = k is true and we must show for n=k+1. From the hypothesis of induction the following is fulfilled. For the k-th optimization, the constraint is now,

akα1+αk+1Γk=1

where ak=j=0k-1Ij. In parameterizing

ϕ=α1Γk;1-akϕ=αk+1Γk

In optimizing

1akϕ+Ik1-akϕ,

we obtain:

ϕk=1ak+Ik.

Solving,

1-akϕkϕk=Ik=αk+1α1.

Therefore,

αk+1=Ikα1

And

ak+1=j=0kIj

And the proof is complete. For the second optimization of the second design rule, the constraint is now,

a2A1+A3u2=A2

Where

a2=1+a1u1a1u1;u2=U1U3;A2=A-i=4nAi

is a constant and fixing the temperature Tj(j=4....n+1). In parameterizing,

ψ=A1A;1-a2ψ=A3A

Where

a2=1+1a1U1.

From (13)

1a2ψ+1(1-a2ψ)u2,

the optimization with respect to ψ, gives,

ψ2=u21+a2u2

And

1-a2ψ2ψ2=1+u2(a2-1)1+a2u2

Thus,

A3=1+u2(a2-1)u2A1

Continuing of this form arrive us to the Eq. (15). In analogous form, the proof of the relations (15) can be carry out by Mathematical Induction.

Received: April 05, 2017; Accepted: August 07, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License