Nonlinear Thermodynamic Analysis and Optimization of a Carnot Engine Cycle

Michel Feidt 1, Monica Costea 2,*, Stoian Petrescu 2 and Camelia Stanciu 2 1 Laboratoire d'Énergétique et de Mécanique Théorique et Appliquée, UMR 7563, Université de Lorraine, 2, avenue de la Forêt de Haye, Vandoeuvre-lès-Nancy 54518, France; michel.feidt@univ-lorraine.fr 2 Department of Engineering Thermodynamics, University Politehnica of Bucharest, Splaiul Independentei 313, Bucharest 060042, Romania; stoian.petrescu@yahoo.com (S.P.); camelia.stanciu10@yahoo.com (C.S.) * Correspondence: monica.costea@upb.ro; Tel.: +40-21-402-9339


Introduction
The oldest and simplest model of a thermo-mechanical engine is due to Carnot [1].This model supposes two thermal reservoirs of infinite capacity (thermostats) at the source and sink, that implies constant temperatures T SH (source) and T SL (sink).
The thermal balance and the entropy balance equations applied to the reversible Carnot direct cycle allow finding the classical result, generally called "Carnot cycle efficiency".This efficiency corresponds to the upper limit of a thermo-mechanical engine efficiency, named First Law efficiency, that is defined as: mechanical work output energy consumption Thus, one recalls the result known as the Carnot cycle efficiency: It is often believed that the Carnot model that relies on Equilibrium Thermodynamics was firstly reconsidered by Curzon and Ahlborn [2], which is not quite precise.Nevertheless, the Finite Time Thermodynamics (FTT) proposes to examine the operation of an endoreversible engine in contact with two thermostats.In this case, the transfer irreversibility appears by the temperature differences at the hot-end and cold-end of the engine (see Figure 1a), respectively, ∆T H = T SH ´TH , ∆T L = T SL ´TL .The configuration of the reversible Carnot engine is optimized with respect to the mechanical work or power output.Thus, the expression of the First Law efficiency results, corresponding to maximum of work or power output as: This result has been previously obtained by Novikov [4].Hence it corresponds to what one may call the Carnot-Novikov-Curzon-Ahlborn engine.Equation ( 2) is often named "nice radical".
A complementary scheme was established in the same time by Chambadal [5,6].The assumption of a steady-state regime of a finite sensible heat source was applied, the source being characterized by a heat rate capacity where THSi is the inlet temperature of the source fluid in the hot-end heat exchanger.The maximum power output of the thermo-mechanical converter is easily obtained [5,6] for the case when considering the sink as thermostat, namely, when TL = TSL (perfect thermal contact at sink).The efficiency associated to maximum power output arises as: One recalls here that the hot source is characterized by a finite heat capacity rate SH C  .
By using the method of irreversibility ratio [13] instead of the entropy analysis approach, the efficiency corresponding to maximum mechanical work output of the irreversible Carnot engine becomes: with I = ΔSL/ΔSH (see Figure 1b).Another alternative to the Finite Time Thermodynamics is the Finite Speed Thermodynamics (FST) proposed by Petrescu and Costea.Reference [14] reports a complete synthesis in this domain.Also, results relative to irreversible Carnot cycle analysis in the frame of FST provide a corrected The configuration of the reversible Carnot engine is optimized with respect to the mechanical work or power output.Thus, the expression of the First Law efficiency results, corresponding to maximum of work or power output as: This result has been previously obtained by Novikov [4].Hence it corresponds to what one may call the Carnot-Novikov-Curzon-Ahlborn engine.Equation ( 2) is often named "nice radical".
A complementary scheme was established in the same time by Chambadal [5,6].The assumption of a steady-state regime of a finite sensible heat source was applied, the source being characterized by a heat rate capacity .

C SH "
. m SH ¨cp SH , with .m SH -mass flow rate of the source fluid, c p SH -mass specific heat at constant pressure assumed independent of temperature.For a perfect heat transfer process to which an effectiveness of the hot-side heat exchanger ε H = 1 corresponds, the heat rate capacity .
C SH provides the maximum heat transfer rate to the endoreversible Carnot engine: .
where T HSi is the inlet temperature of the source fluid in the hot-end heat exchanger.The maximum power output of the thermo-mechanical converter is easily obtained [5,6] for the case when considering the sink as thermostat, namely, when T L = T SL (perfect thermal contact at sink).The efficiency associated to maximum power output arises as: One recalls here that the hot source is characterized by a finite heat capacity rate .
By using the method of irreversibility ratio [13] instead of the entropy analysis approach, the efficiency corresponding to maximum mechanical work output of the irreversible Carnot engine becomes: with I = ∆S L /∆S H (see Figure 1b).Another alternative to the Finite Time Thermodynamics is the Finite Speed Thermodynamics (FST) proposed by Petrescu and Costea.Reference [14] reports a complete synthesis in this domain.Also, results relative to irreversible Carnot cycle analysis in the frame of FST provide a corrected expression of the "nice radical" as a function of the piston speed, correlations between the maximum cycle efficiency and the entropy generation rate, or maximum power output/efficiency, and the piston speed [15][16][17].
This paper proposes reconsidering the Carnot engine approach by taking into account the finite physical dimensions.Heat transfer conductances and areas in regard to the finite speed of the piston, will be particularly considered.
An overview of the possible objective functions aims to give an idea of the complexity of the optimization problem involving as constraints finite temperature difference at source and sink and finite dimensions (heat conductances or areas) of the heat exchangers.However, only two of the listed objective functions will be discussed hereinafter, namely, maximum power output and minimum entropy generation rate.
By using the sequential optimization method, a three-step optimization is achieved corresponding to progressive consideration of different constraints.Hence, the modeling allowed to provide optimum temperatures of the working gas at the hot-end and cold-end of the engine, optimal allocation of heat transfer areas and optimum characteristic speed of the system associated to these objective functions.
Also, it will be shown (which is for the first time, to our knowledge) how the optimization of the heat exchanger finite surfaces is connected to the characteristic speed of the engine and to the energy dissipation, which also depends on the speed.
A fundamental equation is proposed showing the optimal connection between temperatures at the source and sink, irreversibilities of the converter, and characteristic speed, as well as the link to the system size.It is illustrated in a study case, leading to an analytical solution.These original results are discussed and the consequences are shown, as well as the extensions of this work.

Model of Carnot Cycle Engine
Although the case of the reverse cycle machines could be simultaneously approached, it seems preferable to limit the reported study to the case of direct cycles for reasons of results accessibility and demonstrations clarity.
An irreversible engine (converter) is considered characterized by irreversible heat transfer at the contact with the source SH or sink SL.
The source and sink are modeled as thermostats with constant temperatures T SH at the hot-end respectively, T SL at the cold-end.
The study is reported in the stationary dynamic operating regime so that the converter receives the heat flow .Q H > 0 from the source and rejects to the sink the heat flow .Q L < 0, then it produces mechanical power .W > 0 (see Figure 1a).The group of source-converter-sink is assumed without any heat loss (adiabatic system).The energy balance equation of the engine is given by: .
The First Law efficiency of the engine is expressed as: 2.2.Optimization of the CNCA Machine by Using the Entropic Ratio Method

Optimization General Approach
The ratio method is especially used in Finite Speed Thermodynamics and Direct Method [14], which commonly express the mechanical work of any thermal piston machine as: with sign (+) for compression and sign (´) for expansion processes.Hence, there are three ratios in Equation ( 8) that represent dissipative mechanisms accounting for pressure losses due to the finite speed of the piston, friction and throttling of the working fluid [14].Actually, these losses take into account internal irreversibilities of the engine.They can be expressed as functions of the piston speed, design and operational parameters, respectively, and will be considered in further development of the model.
The step proposed here is a global and generic one since it involves only one ratio, I.This ratio appears in the entropy balance equation in the form (9a) or (9b): The case with I H = I L = 1 corresponds to an endoreversible system.Hereinafter, it is chosen to work with the ratio I H , that can be either a parameter (constant on a certain domain), or a function to be identified [18], or sought using the Direct Method, for example, by assuming degradation mechanisms [15].
Most conventional optimizations for an engine are related to the following objective functions: In the present study, η lim is the Carnot cycle efficiency, η C , given by: Note that generally, η II is not identical to the exergetic efficiency, η ex .Actually, η ex is expressed as: Thus, it is obvious that η ex is identical to η II , if T SL = T 0 .This is often close to reality, except certain cases, as for example CHP (Combined Heat and Power) systems.
To the four preceding objective functions given by Equations ( 6), ( 7), ( 10) and (12), it is possible to add those resulting from the energy triangle that is representative of any system or process (Figure 2).This triangle has on the top the Useful Effect provided by the system as a natural optimization choice, while in the bottom the unavoidable Consumption (energy in our case) and Rejects (energy or pollutants) appear.Obviously, the efficiency is inside the triangle, as it is determined by two of the previous fluxes.Then, the environment has to be mentioned, since the ambient properties (mainly temperature) are accounted for (see Equation ( 12)). rejects, R, that may be recoverable (RR, recyclable waste), but most often pollutants for the Environment (RP, pollutant rejects); accordingly, OF 6 namely MIN (R).
The efficiency of the system is represented by the various efficiencies introduced by dimensionless Equations ( 7), ( 10)- (12).
Many other optimization criteria could be considered, including environmental and economic ones [19].These will not be discussed here.
The presentation is mainly concerned with irreversibility in light of its fundamental importance in Thermodynamics.For the case proposed here, the irreversibility of the converter appears as a ratio, IH, which is a parameter, depending on the problem (Equation (9a)).For example, IH = fI(TL, TH).
Hence, Equation (9a) acts as a constraint (C1) between natural intensive variables of the problem TL, TH.The choice of an objective function (OF) of these variables TL, TH (control variables), with C1 as constraint, provides one degree of freedom (DF) leading to the sequential optimization [6].In this case, KH, KL, TSH, TSL are parameters (imposed) of the study, i.e., they appear in OF 6 and constraint, C1.
Furthermore, the objective function, OF 7, can be added.It is related to ISH, for which the minimum, min ISH, will be sought.Note that ISH represents the entropy ratio of the studied system and, therefore, it satisfies the system entropy balance equation: This general optimization approach has reviewed several possible objective functions of the optimization procedure among which only two will be further discussed, namely, maximum power output, OF 1, and minimum entropy production, OF 7.

Maximization of the Power Output of CNCA Engine
Supposing that ISH = f(TH, TL), the Lagrangian function of the optimization can be expressed from Equations ( 6) and (9a) as follows: The necessary condition for the optimum is derived as:

W)
‚ rejects, R, that may be recoverable (RR, recyclable waste), but most often pollutants for the Environment (RP, pollutant rejects); accordingly, OF 6 namely MIN (R).
The efficiency of the system is represented by the various efficiencies introduced by dimensionless Equations ( 7), ( 10)- (12).
Many other optimization criteria could be considered, including environmental and economic ones [19].These will not be discussed here.
The presentation is mainly concerned with irreversibility in light of its fundamental importance in Thermodynamics.For the case proposed here, the irreversibility of the converter appears as a ratio, I H , which is a parameter, depending on the problem (Equation (9a)).For example, Hence, Equation (9a) acts as a constraint (C1) between natural intensive variables of the problem T L , T H .The choice of an objective function (OF) of these variables T L , T H (control variables), with C1 as constraint, provides one degree of freedom (DF) leading to the sequential optimization [6].In this case, K H , K L , T SH , T SL are parameters (imposed) of the study, i.e., they appear in OF 6 and constraint, C1.
Furthermore, the objective function, OF 7, can be added.It is related to I SH , for which the minimum, min I SH , will be sought.Note that I SH represents the entropy ratio of the studied system and, therefore, it satisfies the system entropy balance equation: This general optimization approach has reviewed several possible objective functions of the optimization procedure among which only two will be further discussed, namely, maximum power output, OF 1, and minimum entropy production, OF 7.

Particular Results for the Objective Functions MAX
.

Maximization of the Power Output of CNCA Engine
Supposing that I SH = f(T H , T L ), the Lagrangian function of the optimization can be expressed from Equations ( 6) and (9a) as follows: The necessary condition for the optimum is derived as: This system of two equations with two unknowns needs numerical resolution to get T H , T L corresponding to the general case.
Remark 1.The function, I SH = f(T H , T L ), could be expressed by using FST that is able to provide the entropy generation as a function of the piston speed and temperatures.This improvement makes part of further development of the modeling.
In the present approach, the particular case involves f I = I = ct, f I, H " Then simple analytical calculations provide: leading to: .
Equation (19) shows that the entropy ratio is an increasing function of the converter irreversibility, the minimum ratio of the system being imposed by the temperatures of the thermostats.
Contrarily, the efficiency corresponding to maximum power (Equation ( 20)) is a decreasing function of the converter irreversibility.Equation (18) appears as a function of the conductances, K H and K L .
Taking account of the finite dimension of the conductances K H and K L as a constraint C2 that imposes: one can find by a second sequential optimization that the second order maximum of .
W is associated to Lagrangian L 2 (K H , K L ): Some calculations will lead to the next analytical results: T H ˚" a I T SH ?T SH `?T SL 1 `?I Consequently, Equation ( 18) becomes: Note that Equations ( 19) and (20) are not affected by this second optimization.The second order maximum of .
W is an increasing function of the size of the system by K T , of the temperature difference (T SH -T SL ), but a decreasing function of the irreversibility, I.The result associated with the endoreversible engine is found.

Consideration of System Irreversibility
Equation (27) shows that the maximum power of the system is obtained for the endoreversible converter, but as it has been seen (Equation ( 19)), the irreversibility ratio of the system, I SH , integrates the irreversibility of the converter by the term, ?I (internal irreversibility), and the irreversibility associated with heat transfer, through a T SH {T SL .Both ratios are greater than 1, but yet different.

Condition of transfer entropies equipartition
It is proposed here to study the flow of entropy produced at the hot-end, cold-end, and inside the system.Firstly, it is necessary to clearly differentiate the entropy of transit at the hot-end of the converter, and the entropy of transit at the cold-end compared to the converter entropy generation.Actually, the heat transfer rate at the hot-end can be written as: .
Then the entropy production rate due to heat transfer from the source to the converter, .
S i H , results as: The transposition of the reasoning from the converter to the cold-end leads to: . .

S H ď
. S L , hence there is no equipartition of transit entropies at the two ends of the converter, except for the endoreversible case (trivial).
Remark 2. The equipartition of entropy [20,21] emerged from the linear thermodynamics of the irreversible processes approach when the various fluxes are taken as linear functions of the driving forces.On these bases, structural properties of the optimized processes are established.
However, the entropy production distribution between source and sink provided by the optimization procedure should be analyzed in connection with the model assumptions, as this may differ with their change.
For the particular case of the linear heat transfer law, the expressions of .
S i H and .

S i L result as follows:
. .
The equipartition of the transfer entropies in the general case involves equality of Equations ( 33) and (34) under the condition (9a), leading to a system of two equations with two unknowns to be solved: There is a physical solution to the problem in the endoreversible case: This solution also makes T L eq = T SL (T H eq = T SH ), which is equivalent to the total reversibility condition (equilibrium thermodynamics).
Note that T L eq , T H eq , T H , T L are values of temperatures at transfer entropies equipartition.

Minimization of the entropy production rate of the system
The entropy production rate in the system is expressed as: . S i S " .
In the case of the linear heat transfer law, Equation (39) is written: .
The optimization considering the objective function OF 7, .
S i S , given by Equation (40), is always subject to the constraint, C1 (Equation (9a)).The Lagrangian results as: Problem solving requires a system of two equations with two unknowns, T H , T L .In the particular case where I H = I = constant, we find that the condition minimizing the system entropy production rate is the thermodynamic equilibrium, whatever K H , K L are.Hence, .S i S is an increasing function of T L and a decreasing one of T H .One finds a result obtained for the particular case of a thermoelectric system [22].

General Statement of the Optimization
This entropy production rate method has been preferred in our research for a long time because it is considered appropriate to the identification of an entropy generation function .S relative to the whole or only a part of the system or process (which is the case in Section "Consideration on System Irreversibility").Also, it can lead to particular expressions of functions .S by dissipation mechanisms models, thus joining the Direct Method previously cited [14].
Again, the most comprehensive method is adopted, so that the entropy balance equation of the converter will get a unique form that is equivalent to Equation (9a) or (9b): .
By preserving the same assumptions as above, namely, linear heat transfer laws, the energy balance equation given by Equation ( 6) remains the same.The constraint given by the entropy balance equation of the converter becomes: Actually, it has been shown that .
S C is a function of T H , T L .Moreover, it also depends on the characteristic speed of the system, w.This was experimentally confirmed for reverse cycle machines (heat pumps) [21], and it should also apply to engines.

Entropy production rate of the converter
The entropy production rate of the converter is generally denoted as: .
The simplest mathematical forms of Equation ( 44) that can be used are: -.

S C "
.

S C " C ¨ln T H
T L (logarithmic dependence on temperature ratio) Experiments done on heat pumps have shown that the coupling to the speed is made by the linear law through S', such that:

Heat transfer laws at the source and sink
The most commonly used heat transfer law is the linear one, expressed as: .
The conductances, K i , were supposed constant in Section 2.2, where they were introduced.A step further in the development of a more general model is based on the fact that the literature provides correlations expressed as: Therefore, the functions f i from Equation (46) could be identified as nonlinear heat transfer laws, namely:
with K i -generalized conductances.

W
The Lagrangian function leading to the maximization of the power output of the engine L 4 (T H , T L , w) is composed as follows: There are three variables (T H , T L , w) for one constraint, thus the model has two independent variables (2 degrees of freedom).
The system of three equations to be solved is studied in the particular case where the dependence of K H , K L on temperature is negligible, so that f H (T H , w) = f H (w); f L (T L , w) = f L (w).Furthermore, if one supposes identical modes of heat transfer to the hot-and cold-ends, it comes: This corresponds to the particular case of the dominant convective heat transfer laws (not radiative, nor conductive).
In the particular case of convective heat transfer laws of the Newton type (linear) satisfying the condition (50) and a linearized law of the converter energy dissipation rate settled by Equation (44) [27], it results: f C pT H , T L , wq " S C pwq pT H ´TL q (51) In the sequential optimization relative to the control variables, T H , T L , the following required condition is derived: Note that in Equation ( 52) the variables can be separated.Therefore, by considering the intermediate variable, α, the resolution is reduced to a second order equation as following: Equation ( 54) can be numerically solved, except for the endoreversible case.
As the purpose here is not to present a complete solution, but the existence of an extensive combined solution, the results are illustrated in the case close to endoreversibility.

Sequential Optimization of the Engine Power Output
Asymptotic Solution for the Case S C (w) Smaller than k H A H and k L A L For this case one gets: Furthermore, the maximum power output of the engine is derived: The power output of the engine expressed by Equation (56) depends on k H , k L , A H , A L , S C , T SH , T SL , which are parameters.Optimum temperatures, T H * , T L * ,are obtained by the asymptotic forms of the following expressions: Optimization of the Physical Geometric Dimensions A H , A L , of the System The heat transfer areas, A H , A L , are considered as variables, subject to a constraint of finite dimensions given by the A T parameter.The optimization aims to find the optimal distribution of the heat transfer area between the hot-and cold-ends of the system: Remark 3. Apparently, Equations ( 21) and (59) provide similar finite dimension constraints for the optimization problem.However, they are quite different from the optimal conditions yielded from the optimization procedure.Besides, the finite constraint on heat conductances is simpler and less precise but gives easier analytical results.Thus, both of them are important in the optimization modeling.
The variational calculus allows to find a solution close to the limit and improvable by iteration, such as: with Then, by assuming the formal symmetry of the heat transfer laws (k H (w)= k L (w) = k(w)) and by limiting S C (w) to the first order form, the expression of the maximum power MAX 2 ˇˇ.W ˇˇin the limiting case results: Knowing that Bk Bw " k , w and BS C Bw " S C, w are positive, Equation (64) proves by derivation with respect to w that there exists an optimum power with respect to a characteristic speed of the system, w.It results w*, the optimum characteristic speed that satisfies the equation: This fundamental equation is obtained for the first time to our knowledge.It shows the optimal connection between heat transfers to the source and sink, with the irreversibilities of the converter, as well as the link to the system size (A T ) and temperatures at the source and sink (T SH , T SL ).

Discussion-Partial Conclusion
The particular forms of Equation (64) provide some interesting trends that are discussed hereafter.Note that in the case of purely convective heat transfers, k(w) takes the general form: (a) For a dissipation law, similar to that reported in [18], Equation (51) becomes: Then, from Equation (64), the optimum speed can be expressed as: Remark 4. Generally, the value of n is close to 3, according to the correlations indicated in the literature and Dulong-Petit, a law mainly used by Angulo Brown's group [23,24], for which Equation (68) yields: it is known that w* increases with both T SH and T SL , as well as with the ratio S C '/k', but decreases with A T .Small machines have larger optimum speeds.
(b) For a dissipation law as an mth power function of the speed: Equation (64) leads to the following expression of the optimum speed: Therefore, there is a confirmation of the trends identified on Equation (67).Note that the irreversibility function (Equation ( 68)) satisfies the Second Law of Thermodynamics since the entropy flux vanishes at zero speed (corresponding to a quasi-static process that is characteristic to the Equilibrium Thermodynamics).
Equation (66) that was obtained by identification to experiments can only be a particular relationship, valid in a speed-limited physical domain.However, it accounts well for the irreversibility associated with temperature differences between input and output of the converter (T H -T L ).
It is assumed that a more comprehensive relationship must exist to account for all of irreversibility coupled mechanisms (including ∆P with mass transfers).

Optimization Statement on an Entropic Base
This approach differs a little from the previous one (Section 2.3).It is based on the following fundamental relationship: .
Relation (66) preserves the heat transfer assumption used in Sections 2.2 and 2.3.Note that here, .
S i is the entropy rate associated with transit (or transfer) processes, not the entropy production rate, as it will be seen below.The shaping is almost identical to previous except that in the calculations, a factor, a(w), appears in Equations ( 57) and ( 58) such that: This arises from the variable conductance reasoning, K(w), and not transfer coefficient, k(w), but with the assumption of equipartition of transfer conductances (K H (w) = K L (w) = K(w)).
Then the resolution leads to an asymptotic solution: and it results for MAX ˇˇ.W ˇˇ: This asymptotic solution confirms the results obtained in Section "Optimization of the Physical Geometric Dimensions A H , A L , of the System" and also the observed trends, in addition to the existence of an optimal characteristic speed providing maximum engine power in the presence of irreversibility.

Particularization to Calculate the Various Entropy Production Rates
To illustrate the results, the case of the endoreversible engine ( .S C = 0) is chosen, with K H and K L parameters.One finds analytically: The temperatures corresponding to the optimum power of engine result as: and the maximum power of the endoreversible engine is: Then, it is easy to show that when variable conductances fulfilling the condition K H + K L = K T are considered, equipartition of conductances leads to maximum power, expressed by: This result is in full agreement with that one deduced from Equation (63).It also allows to examine the distribution of the entropy production rates in the most favorable case (endoreversible).Hence, the entropy production rate at the hot-end is: .
The entropy production rate at the cold-end of the endoreversible engine is: .
Concluding, equipartition of the entropy production of heat transfer processes can be encountered only if: On the other side, the total entropy production rate of the endoreversible engine, .

S i T corresponds to:
.
This expression shows that when K H + K L = K T (finite physical dimensions constraint), the total entropy production rate of the endoreversible engine passes through a maximum for equipartition of heat transfer conductances (K H = K L = K T /2).
Thereby, one can conclude that for the endoreversible engine the equipartition of heat transfer conductances corresponds to the maximum mechanical power of the engine, and also to the maximum entropy production rate.This confirms other works' results [7,22,[28][29][30][31] and shows that the entropy generation minimization (EGM) recommended in quite an important number of works [32][33][34] does not have the character of generality we would like to give it.It was widely used in the 1990's as a powerful optimization tool.For some particular constraints, the optimization results have shown equipartition of heat transfer conductances associated with maximum power output.Actually, the results can be completely different when the constraints associated with an objective function are changed.Thus it is necessary to properly raise the issue.Generally, the minimum entropy generation rate is associated with high values of or maximum efficiency (if it exists) or to an economical operation regime of the system.
The incidence of converter irreversibilities affects the results more.Other variables, like the characteristic speed, w, can further modify the results (Section 2.3.2).

Conclusions
This work provides models of the irreversible Carnot engine with finite dimensions (heat transfer conductances or surfaces, characteristic speed of the engine), and local description of the engine irreversibilities through (1) the entropy ratio method of irreversibility and (2) the method of entropy production rate.The aim of the study was to develop the models by considering the engine irreversibility introduced by the two methods, to prove that they can provide optimal regimes and to derive analytical results, at least for particular cases.
Even if a general approach is proposed by the formal optimization based on the variational calculus, the objective functions are focused on the maximum power output of the engine and the minimum entropy production rate in steady state operation conditions.
The method of sequential optimization used in the modeling allowed to obtain the temperatures associated to optimum (T H *, T L *), afterwards the optimal allocation of heat transfer areas (A H *, A L *), and finally w*, the engine speed at the optimum operation point.The main findings related to maximizing the engine power extend previous results reported in the literature.Thus, a new generalized form of the nice radical has been proposed in the paper.
The models proposed herein also allow to demonstrate the relationship between the physical dimensions of the system (T i , K i , A i ) and a characteristic speed of the engine for the first time, to our knowledge.This connection appears highly nonlinear owing to the coupling dimension-speed.The asymptotic behavior near the endoreversibility can provide meaningful analysis results corresponding to realistic tendencies.Primarily, there is an optimum speed associated with the best compromise between transfer and conversion of energy.
Another important conclusion is obtained through the analysis of the associated entropy production rate.Regardless of the method used to characterize the irreversibility (entropy ratio method, or entropy production rate method), the result of the study shows that:

‚
the maximum power supplied by the engine does not correspond systematically to the minimum entropy production rate; ‚ equipartition of entropy production rate is not associated to the maximum power delivered by the engine (endoreversible or not).
Therefore, these results open up a new horizon for the future of Engineering Thermodynamics, also from a fundamental thermodynamic point of view.
The further development of modeling regarding the other objective functions considered in the general optimization procedure, the replacement of the entropic ratio used here as a parameter with an expression derived by using FST for the entropy generation rate as a function of the piston speed, as well as the comparison of the two methods (ratio method, or entropy production rate method), is in progress.
m  -mass flow rate of the source fluid, SH p c -mass specific heat at constant pressure assumed independent of temperature.For a perfect heat transfer process to which an effectiveness of the hot-side heat exchanger εH = 1 corresponds, the heat rate capacity SH C  provides the maximum heat transfer rate to the endoreversible Carnot engine:

Figure 2 .
Figure 2. Energy Triangle connecting the Useful Effect (UE), the Energy Consumption (EC) and the Reject (R).

Figure 2 .
Figure 2. Energy Triangle connecting the Useful Effect (UE), the Energy Consumption (EC) and the Reject (R).
2.3.2.General and Required Condition for the Existence of Optimum Corresponding to MAX