Next Article in Journal
Life’s a Gas: A Thermodynamic Theory of Biological Evolution
Next Article in Special Issue
Thermodynamic Analysis of Closed Steady or Cyclic Systems
Previous Article in Journal
The Intrinsic Cause-Effect Power of Discrete Dynamical Systems—From Elementary Cellular Automata to Adapting Animats
Previous Article in Special Issue
Non-Fourier Heat Transfer with Phonons and Electrons in a Circular Thin Layer Surrounding a Hot Nanodevice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of the Changing Phase Fluid in a Carnot Type Engine for the Recovery of a Given Waste Heat Source

1
Laboratoire d'Énergétique et de Mécanique Théorique et Appliquée, UMR 7563, Université de Lorraine, 54518 Vandoeuvre-lès-Nancy, France
2
Agence de l'Environnement et de la Maîtrise de l'Energie (ADEME), 20, avenue du Grésillé-BP 90406, 49004 Angers Cedex 01, France
*
Author to whom correspondence should be addressed.
Entropy 2015, 17(8), 5503-5521; https://0-doi-org.brum.beds.ac.uk/10.3390/e17085503
Submission received: 10 July 2015 / Revised: 24 July 2015 / Accepted: 27 July 2015 / Published: 31 July 2015
(This article belongs to the Special Issue Selected Papers from 13th Joint European Thermodynamics Conference)

Abstract

:
A Carnot type engine with a changing phase during the heating and the cooling is modeled with its thermal contact with the heat source. In a first optimization, the optimal high temperature of the cycle is determined to maximize the power output. The temperature and the mass flow rate of the heat source are given. This does not take into account the converter internal fluid and its mass flow rate. It is an exogenous optimization of the converter. In a second optimization, the endogenous optimization, the isothermal heating corresponds only to the vaporization of the selected fluid. The maximization of the power output gives the optimal vaporization temperature of the cycled fluid. Using these two optima allows connecting the temperature of the heat source to the working fluid used. For a given temperature level, mass flow rate and composition of the waste heat to recover, an optimal fluid and its temperature of vaporization are deduced. The optimal conditions size also the internal mass flow rate and the compression ratio (pump size). The optimum corresponds to the maximum of the power output and must be combined with the environmental fluid impact and the technological constraints.

1. Introduction

The first thermal model applied to an engine was developed by Carnot [1]. He describes a series of ideal thermodynamic transformations to convert heat into mechanical energy. The main features of his concept are the isothermal heating and the isothermal cooling. He chooses a vapor engine with a phase change and an air piston engine as examples. Carnot introduces a fundamental efficiency limit for converting thermal energy into mechanical energy. The maximal efficiency of a Carnot engine is given by the expression η C = 1 T L S / T H S , where T H S and T L S are, respectively, the temperature of the heat source and the heat sink. He clearly states that the nature of the fluid and its state do not affect this theoretical limit. However, in this approach, the working fluid is supposed to be at the heat source temperature during the isothermal expansion and at the cold sink temperature during the isothermal compression. These assumptions require thermal equilibrium between the working fluid and the heat source or sink during the thermal transfers (Figure 1a). These assumptions imply infinitely slow transformations and, under these conditions, zero power output.
Figure 1. Entropy diagram of a Carnot cycle; (a) for an ideal Carnot engine; (b) for an ideal Carnot engine with thermal transfer consideration.
Figure 1. Entropy diagram of a Carnot cycle; (a) for an ideal Carnot engine; (b) for an ideal Carnot engine with thermal transfer consideration.
Entropy 17 05503 g001
Curzon and Ahlborn consider another fundamental limitation on efficiency that is caused by the rate at which heat can be exchanged between the working fluid and the heat reservoirs [2] (Figure 1b). They introduce an efficiency of the Carnot engine at maximum power output, which is given by the expression η C N C A = 1 T L S / T H S . Historically, this theoretical limitation was introduced by Novikov [3] and, simultaneously, by Chambadal in France [4]. Some recent articles [5,6] refer to an older paper of Yvon [7]. In this paper, Yvon introduced the optimal high temperature, which maximizes the mechanical power output as if it was a widely known and accepted result. These considerations have led to Finite Time Thermodynamics [8]. The Carnot ideal engine was extended to an irreversible one with a finite heat source and heat sink reservoirs (Figure 2) [9,10,11,12]. This system approach allows sizing an optimal Carnot engine independently of the working fluid and its state. The external fluids (heat source and cold sink) are considered limiting. This optimization is described as exogenous, that is, external to the converter. Finite Time Thermodynamics evolved to Optimal Thermodynamics in Finite Physical Dimensions [13].
Figure 2. Entropy diagram of an irreversible Carnot engine with finite heat source and heat sink reservoirs.
Figure 2. Entropy diagram of an irreversible Carnot engine with finite heat source and heat sink reservoirs.
Entropy 17 05503 g002
The working fluid choice is important and has been the subject of many articles, particularly for Rankine cycles [14,15,16,17,18,19,20]. The ideal working fluid must have good thermo-physical properties, but also adequate chemical stability in the desired temperature range. The fluid selection considerably affects the operating conditions, the environment impact, and the cost. In practice, a fluid that satisfies all the ideal fluid criteria does not exist. In this article, an endoreversible changing phase Carnot engine is studied, as shown in Figure 3. The vaporization and the condensation at a constant temperature are means to implement a Carnot engine. There is no technical description of a specific Carnot type engine in the present work; only a proposition of a general design process. It can be adapted to the Rankine engine and also to other technologies operating according to a cycle similar to the Carnot cycle [21].
Figure 3. Changing phase Carnot engine.
Figure 3. Changing phase Carnot engine.
Entropy 17 05503 g003
The model used is an endoreversible engine (Figure 4). The hot fluid enters at temperature T H S i and exits at temperature T H S o . The cold fluid enters at temperature T L S i and exits at temperature T L S o . Heat transfer is taken into account here at the hot side only and not at the cold side. This simplifying assumption allows avoiding the temperature gap at the cold sink. It is equivalent to assuming known the condensation temperature or the condensation pressure in the condenser. Both physical quantities are strongly related to the environment. That is why two cases are covered in this work. Firstly, the low temperature is imposed equal to the ambient temperature T 0 . This corresponds to assuming that T L = T L S i = T L S o = T 0 , which corresponds to Chambadal model [4]. Secondly, the low pressure is imposed equal to the atmospheric pressure. The cycle is optimized to maximize the mechanical power output. The parameters are the fluid characteristics, which are the temperature at the critical point, the temperature at the triple point, and the latent heat of vaporization at the triple point. An identified expression of the latent heat of vaporization is given depending on the parameters and on the variable T H , the temperature of vaporization. This optimization is described as endogenous, that is, internal to the converter and independent of the heat source.
Figure 4. Endoreversible heat engine.
Figure 4. Endoreversible heat engine.
Entropy 17 05503 g004
The last specific feature of this work is a way of linking exogenous and endogenous optima. This allows determining, for a given heat source temperature level and mass flow rate, the adapted working fluid, its optimal temperature of vaporization, and its optimal mass flow rate. The chosen criterion is thermodynamic. This optimum constitutes a tool for choosing the working fluid, even if it cannot be considered independently of environmental, technological, and economic criterions.

2. General Endoreversible Carnot Engine Model

2.1. Assumptions and Balances

The engine diagram is given in Figure 4. The heat source is considered as a finite heat reservoir, the hot fluid enters at T H S i and exits at T H S o . Symmetrically, the cold source is considered as a finite heat reservoir, the cold fluid enters at T L S i and exits at T L S o . The cycled fluid is at T H at the hot side and at T L at the cold side. To simplify, the assumption is added that T L is known. Thus, T L is a parameter as mentioned in Section 1. The steady state hypothesis is used so the entering heat rate is Q ˙ H , the outgoing heat rate is Q ˙ L and the outgoing mechanical power is W ˙ . The endoreversible case is considered, which corresponds to the case where the internal entropy production of the converter is neglected. The thermodynamic convention is used: each quantity going into the system is positive and the quantity going out is negative.
The energy balance is:
W ˙ = Q ˙ H + Q ˙ L
The entropy balance of the converter in the endoreversible case is:
Q ˙ H T H + Q ˙ L T L = 0

2.2. Exogenous Optimization

The results of the general optimization are taken here for a Carnot engine. We can express the entering heat rate as:
Q ˙ H = C ˙ H ( T H S i T H S o )
C ˙ H is the heat source calorific rate:
C ˙ H = m ˙ H S c p H S
m ˙ H S is the mass flow rate of the heat source and c p H S is the calorific capacity, assumed not variable during the heat transfer.
The exchanger is supposed to be an ideal one, so the output temperature is equal to the high temperature of the internal fluid. The exchanger effectiveness is supposed to equal 1. We can write:
Q ˙ H = C ˙ H ( T H S i T H )
Injecting (2) and (3) in (1), the power output is given by:
W ˙ = C ˙ H ( T H S i T H ) ( 1 T L T H )
For C ˙ H , T H S i , and T L given, the optimal high temperature at maximum power output is the solution of d ( W ˙ ) / d T H = 0   [5]:
T H o p t e x o = T L T H S i
The corresponding maximum power output is:
W ˙ e x o   = C ˙ H ( T H S i T L ) 2
The entering heat rate becomes:
Q ˙ H e x o = C ˙ H T H S i ( T H S i T L )
The outgoing heat rate is:
Q ˙ L e x o = C ˙ H T L ( T L T H S i )
These three energy flows (Equations (8) to (10)) can be made dimensionless:
T H S i = t h s i T 0   ;   T L = t l T 0     ;   Q ˙ = q ˙ C ˙ H T 0   ;   W ˙ = w ˙   C ˙ H T 0
The expression of the optimal dimensionless temperature is:
t h o p t e x o = t l t h s i
The expression of the dimensionless power output is:
w e x o   = ( t h s i t l ) 2
The dimensionless entering heat rate is:
q h e x o = t h s i ( t h s i t l )
The dimensionless outgoing heat rate is:
q l e x o = t l ( t l t h s i )
The optimal dimensionless temperature, the dimensionless power output, the dimensionless entering heat rate, and the dimensionless outgoing heat rate are plotted as a function of t h s i in Figure 5a. The low temperature is fixed at an ambient temperature, assumed equal to 300   K . This choice of ambient temperature must be adapted to the considered environment.
Figure 5. (a) Optimal high temperature, power output, entering heat flow rate and outgoing heat flow rate; (b) first law and second law efficiency behavior; for inlet temperature heat source variable and low temperature imposed equal to T 0 = 300   K .
Figure 5. (a) Optimal high temperature, power output, entering heat flow rate and outgoing heat flow rate; (b) first law and second law efficiency behavior; for inlet temperature heat source variable and low temperature imposed equal to T 0 = 300   K .
Entropy 17 05503 g005
The first law efficiency corresponding to the power output divided by the entering heat rate is:
η I = W ˙ Q ˙ H = 1 T L T H S i
The second law quality factor is:
f I I = η I 1 T L T H S i = 1 1 + T L T H S i
The evolution of the first law efficiency and the second law quality factor is given in Figure 5b for a Chambadal model ( T L = T 0 ) [4].

3. Changing Phase Carnot Engine

3.1. Evolution of the Latent Heat of Vaporization with the Temperature

In this section, a changing phase Carnot engine is considered (Figure 3). The transformation heat is supposed to be the latent heat of vaporization at T H . It imposes knowing the way that the latent heat of vaporization varies with T H . For a given fluid, the latent heat of vaporization is equal to zero at the critical temperature T C and is supposed to be known at the temperature of the triple point T T :
L ( T T ) = L T T
The latent heat of vaporization decrease with an increased of temperature can be modeled in a reasonable way using the following analytical equation:
L ( T ) = L T T ( T C T T C T T ) 1 n
For a linear approximation of the latent heat of vaporization, n is equal to 1. In this approximation, if the saturation dome of the considered fluid is symmetric, the latent heat of vaporization is shown in Figure 6. The increase of the n value is also shown in the figure, in this case there is an adapted n value of between 2 and 3. A large selection of fluids with their main properties is given in Table 1. The properties are taken from the software Refprop. The considered fluids are chosen because their temperatures at the triple point are below 300  K . They are also chosen for their variable critical temperatures, ranked in increasing order, in the table. The knowledge of the latent heat of vaporization variation with temperature for each fluid allows an optimal estimation of n   by a least square regression, of which output is plotted together with the reference data. The comparisons for inorganic pure fluids are given in Figure 7. Figure 7a is the latent heat of vaporization of water and Figure 7b that of ammonia. A large selection of hydrocarbons is represented in Figure 8 ((a) butane, (b) isobutene, (c) cyclopentane, and (d) cyclohexane) and in Figure 9 ((a) propane, (b) benzene, (c) decane, (d) toluene, and (e) acetone). Halogenated hydrocarbons are also considered. Some HCFC (HydroChloroFluoroCarbons) are given in Figure 10 ((a) R22 and (b) R123) and some HFC (HydroFluoroCarbons) are given in Figure 11 ((a) R32, (b) R125, (c) R134a, (d) R143a, and (e) R152a) and Figure 12 ((a) R227ea and (b) R245fa). Approximation (21) with different n values is acceptable for all of the selected fluids. The relative bias is given by:
e L ( T H ) = | L a p p r o x i m ( T H ) L e x p e r i m e n t a l ( T H ) | L a p p r o x i m ( T H )
Figure 6. Approximation of the latent heat of vaporization for a symmetric dome of saturation.
Figure 6. Approximation of the latent heat of vaporization for a symmetric dome of saturation.
Entropy 17 05503 g006
Table 1. Properties of the fluid selection.
Table 1. Properties of the fluid selection.
SubstanceMolecular formulaPhysical dataApproximation data
TC (K)TC (K) L T T (J/kg) n value e L m a x (%) e L a v e (%)
R125 C H F 3 C F 3 3391731.90 × 1052.62.20.5
R143a C 2 H 2 F 3 3461612.67 × 1052.61.60.4
R32 C H 2 F 2 3511364.63 × 1052.62.70.5
R22 C H C l F 2 3691163.03 × 1052.52.21.1
propane C 3 H 8 370865.63 × 1052.61.71
R134a C 2 H 2 F 4 3741702.63 × 1052.53.00.6
R227ea C 3 H F 7 3751461.75 × 1052.52.31.6
R152a C 2 H 4 F 2 3861554.06 × 1052.51.70.6
ammonia N H 3 4051951.48 × 1062.64.41.4
isobutane C 4 H 10 4081144.81 × 1052.61.50.7
butane C 4 H 10 4251354.96 × 1052.61.30.7
R245fa C 3 H 3 F 5 4271712.59 × 1052.3163
R123 C 2 H C l 2 F 3 4571662.24 × 1052.48.91.8
acetone C H 3 C O C H 3 5081796.49 × 1052.451
cyclopentane C 5 H 10 5121804.87 × 1052.5252.2
cyclohexane C 6 H 12 5542794.04 × 1052.51.50.7
benzene C 6 H 6 5622794.48 × 1052.51.20.6
toluene C 7 H 8 5921784.92 × 1052.35.82.2
decane C 10 H 22 6182443.93 × 1052.32.71.4
water H 2 O 6472732.50 × 1062.9121.7
The maximum average bias is 3% and relative to R245fa, which is acceptable (Table 1). The maximum relative error is for the cyclopentane (25%), which is high but is very close to the critical point value. At this point, the latent heat of vaporization is close to zero, so a small absolute error becomes a very large relative error. The other maximum relative bias has the same cause.
Figure 7. Approximated latent heat of vaporization for (a) water and (b) ammonia.
Figure 7. Approximated latent heat of vaporization for (a) water and (b) ammonia.
Entropy 17 05503 g007
Figure 8. Approximated latent heat of vaporization for (a) butane; (b) isobutane; (c) cyclopentane; (d) cyclohexane.
Figure 8. Approximated latent heat of vaporization for (a) butane; (b) isobutane; (c) cyclopentane; (d) cyclohexane.
Entropy 17 05503 g008aEntropy 17 05503 g008b
Figure 9. Approximated latent heat of vaporization for (a) propane; (b) benzene; (c) decane; (d) toluene; (e) acetone.
Figure 9. Approximated latent heat of vaporization for (a) propane; (b) benzene; (c) decane; (d) toluene; (e) acetone.
Entropy 17 05503 g009
Figure 10. Approximated latent heat of vaporization for (a) R22 and (b) R123.
Figure 10. Approximated latent heat of vaporization for (a) R22 and (b) R123.
Entropy 17 05503 g010
Figure 11. Approximated latent heat of vaporization for (a) R32; (b) R125; (c) R134a; (d) R143a; and (e) R152a.
Figure 11. Approximated latent heat of vaporization for (a) R32; (b) R125; (c) R134a; (d) R143a; and (e) R152a.
Entropy 17 05503 g011
Figure 12. Approximated latent heat of vaporization for (a) R227ea and (b) 245fa.
Figure 12. Approximated latent heat of vaporization for (a) R227ea and (b) 245fa.
Entropy 17 05503 g012

3.2. Endogenous Optimization

The entering heat rate depends on the latent heat of vaporization L ( T H ) :
Q ˙ H = m ˙ L ( T H )
m ˙ is the cycled fluid mass flow rate.
The outgoing mechanical power can be expressed with (2) and (21) as:
W ˙ = m ˙ L ( T H ) ( 1 T L T H )
To maximize the outgoing mechanical power we have to solve:
( W ˙ ) T H | m ˙ = 0
The optimal high temperature, taking Equation (19) into account, is:
T H o p t e n d o = 1 2 T L ( ( n 1 ) 2 + 4 n T C T L ( n 1 ) )
n depends on the fluid and is between 2 and 3. T C is the critical temperature of the considered fluid and T L is assumed to be known. The result means that, for a given fluid, we know at which temperature the fluid must be evaporated to deliver the maximum power for an endoreversible cycle and a fixed condensation temperature. To view this result, the variation of the area of the cycle versus the evaporation temperature is shown in Figure 13 for a wet fluid. In Figure 13a, this area is represented in the entropy diagram. In Figure 13b, the area is directly plotted versus temperature T H . For T H = T L or for T H = T C , the area is null. The optimal temperature is the temperature at maximal power output, which is equivalent to the maximum area. It is also valid for dry and isentropic fluids because it is only the vaporization which fixes the length of the rectangle. The height of the rectangle depends on the temperature difference. The optimization corresponds to the optimal allocation of both sides to maximize the surface area of the rectangle. Recalled that the low temperature is fixed.
Figure 13. Specific work variation with the high temperature (a) representation in the entropy diagram; (b) corresponding behavior.
Figure 13. Specific work variation with the high temperature (a) representation in the entropy diagram; (b) corresponding behavior.
Entropy 17 05503 g013
The maximum power output is:
W ˙ e n d o = m ˙   L T T ( T C T H o p t e n d o T C T T ) 1 n ( 1 T L T H o p t e n d o )
The entering heat rate becomes:
Q ˙ H e n d o = m ˙   L T T ( T C T H o p t e n d o T C T T ) 1 n
The outgoing heat rate is:
Q ˙ L e n d o = m ˙   L T T T L T H o p t e n d o ( T C T H o p t e n d o T C T T ) 1 n

3.2.1. Temperature T L Equal to the Ambient Temperature T 0

For the first case studied, the low temperature T L is fixed and equal to the ambient one T 0 . The variation of the optimal dimensionless temperature t h o p t e n d o = T H o p t e n d o / T 0 with the increase of the critical dimensionless temperature t c = T C / T 0 is plotted in Figure 14 for different fluids. The vaporization temperature and the corresponding compression ratio for the fluid selection are listed in Table 2. Note that the condensation pressure at atmospheric temperature is very low for some fluids ( P L   < 5 × 10 3   Pa for water, decane and toluene), and quite low for acetone, cyclopentane, cyclohexane and benzene ( P L < 5 × 10 4   Pa ). This makes the implementation difficult because compression ratios are not physically acceptable. Moreover, for some fluids (R125, R134a, R32, R22, propane, R134a, R152, and ammonia) the low pressure is higher than 5 × 10 5   Pa . This situation is not an ideal one because it induces high evaporation pressure, which imposes some technological constraints (principally the choice of exchanger and the maintenance cost).
Figure 14. Optimal high temperature for variable critical temperature with T L = 300   K .
Figure 14. Optimal high temperature for variable critical temperature with T L = 300   K .
Entropy 17 05503 g014
Table 2. Optimal vaporization temperature and corresponding compression ratio for fluid selection.
Table 2. Optimal vaporization temperature and corresponding compression ratio for fluid selection.
SubstanceTL imposed at 300 KPL imposed at 1 × 105 Pa
TL (=T0)PL T H o p t e n d o PH τ TLPL (=P0) T H o p t e n d o PH τ
R1253001.45 × 106327 ± 22.80 × 1061.92251 × 105300 ± 11.46 × 10615
R143a3001.32 × 106332 ± 12.81 × 1062.12261 × 105305 ± 11.49 × 10615
R323001.77 × 106335 ± 14.17 × 1062.42211 × 105306 ± 12.08 × 10621
R223001.10 × 106347 ± 13.27 × 10632321 × 105320 ± 11.83 × 10618
propane3009.98 × 105348 ± 12.86 × 1062.92311 × 105321 ± 11.65 × 10616
R134a3007.03 × 105351 ± 12.49 × 1063.62471 × 105330 ± 11.55 × 10616
R227ea3004.81 × 105351 ± 11.77 × 1063.72571 × 105334 ± 11.21 × 10612
R152a3006.30 × 105358 ± 12.65 × 1064.22491 × 105339 ± 11.72 × 10617
ammonia3001.06 × 106371 ± 16.04 × 1065.72401 × 105346 ± 13.55 × 10635
isobutane3003.70 × 105372 ± 11.98 × 1065.42611 × 105357 ± 11.46 × 10615
butane3002.58 × 105383 ± 11.87 × 1067.32721 × 105372 ± 11.51 × 10615
R245fa3001.59 × 105382 ± 11.55 × 106102881 × 105377 ± 11.39 × 10614
R1233009.78 × 104401 ± 11.41 × 106143011 × 105401 ± 11.42 × 10614
acetone3003.33 × 104432 ± 11.31 × 106393291 × 105444 ± 11.68 × 10617
cyclopentane3004.55 × 104434 ± 11.44 × 106323221 × 105444 ± 11.71 × 10617
cyclohexane3001.41 × 104458 ± 11.04 × 106743531 × 105483 ± 11.56 × 10616
benzene3001.38 × 104462 ± 11.20 × 106873531 × 105488 ± 11.80 × 10618
toluene3004.18 × 103473 ± 17.67 × 1061843831 × 105515 ± 11.50 × 10615
decane3002.07 × 102492 ± 12.49 × 10612024471 × 105558 ± 18.90 × 1068.9
water3003.54 × 103507 ± 13.62 × 10610233731 × 105554 ± 16.51 × 10665

3.2.2. Pressure P L Equal to the Ambient Pressure P 0

A second optimization is possible. The low pressure is fixed to the atmospheric pressure P 0 and the low temperature is the temperature of condensation at the atmospheric pressure. T L is deduced from the nature of the fluid (especial temperature of vaporization at ambient pressure). The optimal results are given in Figure 15 and Table 2. The condensation temperature under these conditions is equal to T 0 for R123, both hypotheses are equivalent for this fluid. For a large range of selected fluids, the low temperature is lower than the ambient one (R125, R143a, R32, R22, propane, R134a, R227ea, R152a, ammonia, isobutene, butane, and R245fa). That is not physically acceptable. A compromise between low temperature and low pressure is necessary for getting physically and technologically realistic conditions. This compromise depends on the selected fluid.
Figure 15. Optimal high temperature for variable critical temperature with P L = 1.0 × 10 5   Pa .
Figure 15. Optimal high temperature for variable critical temperature with P L = 1.0 × 10 5   Pa .
Entropy 17 05503 g015

4. Endogenous and Exogenous Optimization of a Carnot Engine

The external fluid is considered as a limiting factor so the energy balances are the same as with the exogenous case (Equations (5) and (6)). The Objective Function ( O F ) to maximize is the mechanical power is:
O F = W ˙ = C ˙ H ( T H S i T H ) ( 1 T L T H )
There is a supplementary constraint C on the mass flow rate. The heat rate transferred by the heat source to the cycled fluid is the energy of vaporization. Equations (5) and (21) are equals:
C = C ˙ H ( T H S i T H ) m ˙ L ( T H ) = 0
L ( T H ) can be replaced by its value (Equation (19)). The system has an optimum, which maximizes the mechanical power output for the variables T H and m ˙ . The Lagrangien ( L g ) is:
L g = O F + λ C = C ˙ H ( T H S i T H ) ( 1 T L T H ) + λ [ C ˙ H ( T H S i T H ) m ˙ L ( T H ) ]
The system to solve is:
{ L λ = 0 L T H = 0 L m ˙ = 0
The optimal solutions are:
{ T H o p t = T H S i T L m ˙ o p t = C ˙ H T H S i ( T H S i T L ) L T T ( ( T C T H S i T L ) T C T T ) 1 n
The dimensionless mass flow rate m r o p t is:
m r o p t = m ˙ o p t L T T T 0 C ˙ H = T H S i ( T H S i T L ) T 0 ( ( T C T H S i T L ) T C T T ) 1 n
Recall that the assumption of an ideal exchanger is made, it means that the exchange area is infinite. This assumption affects the optimal mass flow rate.
The dimensionless high temperature t h o p t is:
t h o p t = T H o p t T 0 = T H S i T L T 0
The corresponding maximal power output, the entering heat rate, and the outgoing heat rate are the same as the exogenous optimization (Equations (8)–(10)). Coupling the endogenous and exogenous equation permits to add a constraint to optimize the mass flow rate of the cycled fluid.
To have the better fluid for a given heat source, we seek to determine the optimal critical temperature depending on the entering heat source temperature. For this, the optimization results are equalized to those of the endogenous optimization. Indeed, the high temperature must satisfy:
T H o p t = T H o p t e n d o
To satisfy (35), the heat source temperature T H S i and the critical temperature of the working fluid T C are linked by the relation:
T C = T H S i n ( T L ( n 1 ) + T H S i )
This relation depends on the fluid. The dimensionless critical temperature is:
t c = T C T 0
We can see in Figure 16a that the behavior of the dimensionless optimal critical temperature t c and the optimal high temperature of the cycled fluid t h o p t for variable entering heat source temperature t h s i = T H S i / T 0 , for different fluids and T L = 300   K . The corresponding dimensionless mass flow rate is given in Figure 16b. The corresponding values are listed in Table 3. For the case where P L is fixed to 1 × 10 5   Pa , the behavior of the results are given in Figure 17a,b.
Figure 16. (a) Optimal critical temperature and optimal vaporization temperature; (b) internal mass flow rate; for entering heat source temperature variable and T L = 300   K .
Figure 16. (a) Optimal critical temperature and optimal vaporization temperature; (b) internal mass flow rate; for entering heat source temperature variable and T L = 300   K .
Entropy 17 05503 g016
Figure 17. (a) Optimal critical temperature and optimal vaporization temperature; (b) internal mass flow rate; for entering heat source temperature variable and P L = 1.0 × 10 5   Pa .
Figure 17. (a) Optimal critical temperature and optimal vaporization temperature; (b) internal mass flow rate; for entering heat source temperature variable and P L = 1.0 × 10 5   Pa .
Entropy 17 05503 g017
Table 3. Optimal high temperature, heat source temperature and mass flow rate for different fluids.
Table 3. Optimal high temperature, heat source temperature and mass flow rate for different fluids.
SubstanceTL imposed at 300 KPL imposed at 1 × 105 Pa
T H o p t T H S i m r o p t T H o p t T H S i m r o p t
R125327 ± 23580.280300 ± 14010.589
R143a332 ± 13680.323305 ± 14120.635
R32335 ± 13760.368306 ± 14240.713
R22347 ± 14020.487320 ± 14430.786
propane348 ± 14040.504321 ± 14480.830
R134a351 ± 14100.467330 ± 14410.681
R227ea351 ± 14110.491334 ± 14360.675
R152a358 ± 14300.538339 ± 14620.756
ammonia371 ± 14600.593346 ± 15000.835
isobutane372 ± 14640.686357 ± 14890.863
butane383 ± 14910.757372 ± 15090.880
R245fa384 ± 14880.736377 ± 14940.783
R123402 ± 15370.890402 ± 15370.888
acetone432 ± 16171.14444 ± 15981.02
cyclopentane434 ± 16281.16445 ± 16151.06
cyclohexane458 ± 16981.22483 ± 16601.01
benzene462 ± 17121.26488 ± 16741.06
toluene478 ± 17501.57516 ± 16941.22
decane492 ± 17881.59558 ± 16961.02
water507 ± 18931.80554 ± 18241.45
For example, for a heat source of 1   kg / s of industrial flue gas at 410   K , the calorific rate is equal to 1000  W / K . For a low temperature equal to the ambient temperature at 300   K , the optimal high temperature is 351   K . The corresponding fluid is the R134a. The optimal mass flow rate of internal fluid is 0.537   kg / s .

5. Conclusion

A Carnot type engine is modeled and optimized. The heating corresponds to the vaporization of the internal working fluid and the cooling of the condensation. The temperature of condensation is supposed known and the internal dissipations (entropy production) are neglected. A known optimization of a Carnot engine with the Chambadal model is recovered and completed by an endogenous optimization. This second optimization takes into account the working fluid and the behavior of its latent heat of vaporization. Making these two optima equal allows linking the heat source to the internal cycle. For a given temperature, mass flow rate, and composition of the waste heat to recover, an optimal fluid and its temperature of vaporization is proposed. The optimal conditions size also the internal mass flow rate and the compression ratio (pump size). The optimum corresponds to the maximum of power output, and must be combined with the environmental fluid impact and the technological constraints. These results have been proposed for pure fluids. Extension to mixtures is envisioned.

Acknowledgments

I thank the company Hevatech (Malataverne, France) and the Agence de l'environnement et de la maîtrise de l'énergie (ADEME, Angers, France).

Author Contributions

Michel Feidt designed and supervised the research. Denis Maillet oversaw and improved the approach. Mathilde Blaise developed the calculations, plotted the results and wrote the article. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Carnot, S. Réflexions sur la puissance motrice du feu et sur les machines propres à développer cette puissance; Bachelier: Paris, France, 1824. (in French) [Google Scholar]
  2. Curzon, F.L.; Ahlborn, B. Efficiency of a Carnot engine at maximum power output. Am. J. Phys. 1975, 1, 22–24. [Google Scholar] [CrossRef]
  3. Novikov, I. The efficiency of atomic power stations (a review). J. Nucl. Energy 1958, 7, 125–128. [Google Scholar] [CrossRef]
  4. Chambadal, P. Le choix du cycle thermique dans une usine génératrice nucléaire. Revue Générale de l'Electricité 1958, 67, 332–345. (in French). [Google Scholar]
  5. Vaudrey, A.; Lanzetta, F.; Feidt, M.H.B. Reitlinger and the origins of the efficiency at maximum power formula for heat engines. J. Non-equilib. Thermodyn. 2014, 39, 199–203. [Google Scholar] [CrossRef]
  6. Moreau, M.; Gaveau, B.; Schuman, S. Stochastic dynamics, efficiency and sustainable power production. Eur. Phys. J. D 2011, 62, 67–71. [Google Scholar] [CrossRef]
  7. Yvon, J. The Saclay reactor: Two years experience on heat transfer by means of compressed gas. In Proceedings of the International Conference on Peaceful Uses of Atomic Energy, Geneva, Switzerland, 8–20 August 1955.
  8. Andresen, B. Finite-Time Thermodynamics; Physics Laboratory II, University of Copenhague: Copenhague, Denmark, 1983. [Google Scholar]
  9. Feidt, M.; Costea, M.; Petre, C.; Petrescu, S. Optimization of the direct Carnot cycle. Appl. Therm. Eng. 2007, 27, 829–839. [Google Scholar] [CrossRef]
  10. Chen, J.; Yan, Z.; Lin, G.; Andresen, B. On the Curzon–Ahlborn efficiency and its connection with the efficiencies of real heat engines. Energy Convers. Manag. 2001, 42, 173–181. [Google Scholar] [CrossRef]
  11. Kodal, A.; Sahin, B.; Yilmaz, T. A comparative performance analysis of irreversible Carnot heat engines under maximum power density and maximum power conditions. Energy Convers. Manag. 2000, 41, 235–248. [Google Scholar] [CrossRef]
  12. Agnew, B.; Anderson, A.; Forst, T.H. Optimization of a steady state flow Carnot cycle with external irreversibilities for maximum specific output. Appl. Therm. Eng. 1997, 17, 3–15. [Google Scholar] [CrossRef]
  13. Feidt, M. Thermodynamique Optimale en Dimensions Physiques Finies; Lavoisier: Paris, France, 2013. [Google Scholar]
  14. Tchanche, B.F.; Papadakis, G.; Lambrinos, G.; Frangoudakis, A. Fluid selection for a low-temperature solar organic Rankine cycle. Appl. Therm. Eng. 2009, 29, 2468–2476. [Google Scholar] [CrossRef]
  15. Chen, H.; Goswami, D.Y.; Stefanakos, E.K. A review of thermodynamic cycles and working fluids for the conversion of low-grade heat. Renew. Sustain. Energy Rev. 2010, 14, 3059–3067. [Google Scholar] [CrossRef]
  16. Hung, T.C.; Wang, S.K.; Kuo, C.H.; Pei, B.S.; Tsai, K.F. A study of organic working fluids on system efficiency of an ORC using low-grade energy sources. Energy 2010, 35, 1403–1411. [Google Scholar] [CrossRef]
  17. Badr, O.; Probert, S.D.; O’Callaghan, P.W. Selecting a working fluid for a Rankine-cycle engine. Appl. Energy 1985, 21, 1–42. [Google Scholar]
  18. Liu, B.-T.; Chien, K.-H.; Wang, C.C. Effect of working fluids on organic Rankine cycle for waste heat recovery. Energy 2004, 29, 1207–1217. [Google Scholar] [CrossRef]
  19. West, H.H.; Patton, J.M. Selection of working fluids for the organic Rankine cycle. In Proceedings of the First Industrial Energy Technology Conference, Houston, TX, USA, 22–25 April 1979; pp. 953–959.
  20. Le, V.L.; Feidt, M.; Kheiri, A.; Pelloux-Prayer, S. Performance optimization of low-temperature power generation by supercritical ORCs (Organic Rankine Cycles) using low GWP (Global Warming Potential) working fluids. Energy 2014, 67, 513–526. [Google Scholar]
  21. De Witt, B.; Hugo, R. Naturally-Forced Slug Flow Expander for Application in a Waste-Heat Recovery Cycle. Energies 2014, 7, 7223–7244. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Blaise, M.; Feidt, M.; Maillet, D. Optimization of the Changing Phase Fluid in a Carnot Type Engine for the Recovery of a Given Waste Heat Source. Entropy 2015, 17, 5503-5521. https://0-doi-org.brum.beds.ac.uk/10.3390/e17085503

AMA Style

Blaise M, Feidt M, Maillet D. Optimization of the Changing Phase Fluid in a Carnot Type Engine for the Recovery of a Given Waste Heat Source. Entropy. 2015; 17(8):5503-5521. https://0-doi-org.brum.beds.ac.uk/10.3390/e17085503

Chicago/Turabian Style

Blaise, Mathilde, Michel Feidt, and Denis Maillet. 2015. "Optimization of the Changing Phase Fluid in a Carnot Type Engine for the Recovery of a Given Waste Heat Source" Entropy 17, no. 8: 5503-5521. https://0-doi-org.brum.beds.ac.uk/10.3390/e17085503

Article Metrics

Back to TopTop