MaximumPower Point Characteristics of Generalized Heat Engines with Finite Time and Finite Heat Capacities

We revisit the problem of optimal power extraction in four-step cycles (two adiabatic and two heat-transfer branches) when the finite-rate heat transfer obeys a linear law and the heat reservoirs have finite heat capacities. The heat-transfer branch follows a polytropic process in which the heat capacity of the working fluid stays constant. For the case of ideal gas as working fluid and a given switching time, it is shown that maximum work is obtained at Curzon-Ahlborn efficiency. Our expressions clearly show the dependence on the relative magnitudes of heat capacities of the fluid and the reservoirs. Many previous formulae, including infinite reservoirs, infinite-time cycles, and Carnot-like and non-Carnot-like cycles, are recovered as special cases of our model.


Introduction
Curzon-Ahlborn efficiency, η CA = 1 − T 2 /T 1 , where T 1 and T 2 are the reservoir temperatures [1], is regarded as a landmark result of finite-time thermodynamics.It models the effect of irreversibilities due to finite rate of heat transfer on the performance of heat engines.Such models are termed as endoreversible.The flow of heat between the working fluid and the reservoirs is assumed to be Newtonian.All the components maintain an internal equilibrium and losses due to friction, and heat leak are assumed to be negligible.Many authors [2][3][4][5][6][7][8][9][10] extended and clarified the scope of this model.In fact, various other models of heat engines also yield maximum power output at efficiencies very close to CA efficiency [11][12][13][14][15].This universal-like behaviour of efficiency has been recently analysed using Bayesian probabilities in the case of quantum models of engines [16].
In this paper, we revisit the problem of optimal performance with "linear" irreversibilities of finite time and finite heat reservoirs in classical models of engines.This question was addressed in [7] using a Lagrangian formalism, for a one-component working fluid without assuming an equation of state.However, the generic problem with finite reservoirs could not be solved in a closed form.Then Gordon [8] used an ancillary device of intermediate reservoirs to arrive at a closed form solution.In particular, it was shown [8] that for Carnot-like engines, finiteness of reservoirs has no effect on the efficiency at maximum power, and it is still at Curzon-Ahlborn efficiency.The case of other non-Carnot heat engines (such as Otto cycle, Joule-Brayton cycle, and Atkinson cycle) was discussed by using an infinite chain of Carnot cycles.
The present analysis of the problem is based on the following features: (a) the working fluid is a classical ideal gas; (b) the total cycle time as well as switching time is given.We consider a generic four-step cycle with two adiabatic and two heat-transfer branches, which follow a polytropic process with a constant heat capacity C > 0. The process is described by the relation TV x = k.Here, T is the temperature of the working fluid, V is the volume of the working fluid, and k is an arbitrary constant.Consequently, many popular heat engines like Carnot cycle (x = 0), Otto cycle (x → ∞), and Joule-Brayton (x = −1) can be incorporated in this study.
Our key results may be listed as follows.For ideal gas as the working fluid and the generic four-step cycle (including

Journal of Thermodynamics
Carnot and non-Carnot heat cycles), the power of the engine is maximum at CA efficiency, even for finite reservoirs.Many special cases like infinite time and infinite reservoirs [11], finite time and infinite reservoirs with isothermal branches [1] or polytropic branches can be described.We note that closed form expressions for work can be derived without invoking a device like an intermediate reservoir [8].Further, the heat capacities of the working fluid and the (finite) reservoirs enter explicitly in the expressions, and their relative sizes appear as the respective ratios of heat capacities.Also, the explicit expressions of work for these special cases clearly show that the irreversibilities of finite time and/or finite reservoirs reduce the maximum amount of work extracted as compared with infinite time and infinite reservoirs cases.

The Model
2.1.Temperature Profiles.First of all, we calculate how the temperatures of the two bodies change, when kept in contact with one another for some time.The heat transfer between the bodies is Newtonian, and during heat transfer, their heat capacities remain constant.
Let the two bodies denoted as A and B have heat capacities, C A and C B , respectively.The bodies are kept in thermal contact with one another from time, t = 0 to t = t 1 .The temperatures of the two bodies at any time t are T A (t) and T B (t), respectively, and the coefficient of heat conductivity is K.
At time t, the heat flow between bodies A and B should be equal to heat gained (lost) by one or heat lost (gained) by the other.Therefore, Solving the above equations, temperature profiles of the two bodies at time t are calculated as where 2.2.The Cycle.We have two finite heat reservoirs at initial temperatures T 1 (0) and T 2 (0) and heat capacities are C 1 and C 2 , respectively.A working fluid is chosen to extract work from the reservoirs.The time period for one working cycle is t 2 .One working cycle consists of the following steps.
(i) Working fluid is brought in contact with the hot reservoir from time t = 0 to t 1 .The fluid expands along a path with heat capacity C. Consequently, temperature of the hot reservoir reduces from T 1 (0) to T 1 (t 1 ), and that of the working fluid rises from T(0) to T − (t 1 ).
(ii) The fluid is allowed to expand adiabatically.This step is assumed to be instantaneous and takes negligible time.As a result, the temperature of the working fluids jumps from T − (t 1 ) to T + (t 1 ).
(iii) The fluid is brought in contact with the cold reservoir from time t 1 to t 2 .The fluid contracts along a path with heat capacity C. Consequently, temperature of the cold reservoir rises from T 2 (t 1 ) to T 2 (t 2 ) and that of the working fluid reduces from T + (t 1 ) to T − (t 2 ).
(iv) Working fluid is allowed to contract adiabatically.This step is also assumed to take negligible time.As a result, the temperature of the working fluids jumps from T − (t 2 ) to T(0), completing the cycle.
Using ( 2)-( 6), the temperature profiles of the reservoirs and the working fluid in different steps can be described as follows: (i) 0 < t < t 1 : where Since T 2 (t 1 ) = T 2 (0), this implies  where )), The complete cycle is depicted in Figure 1.
We will now use the above temperatures to evaluate the work performed by the engine.Due to the cyclic process, the first law of thermodynamics implies that the net heat exchanged by the working fluid during the cycle equals the work performed.Thus CdT. ( Since specific heat is not a function of temperature, Using ( 9), (10), and ( 14), we obtain Substituting values for λ, λ 2 , β and β 1 from ( 8) and ( 13), we have Heat absorbed from the hot reservoir (Q 1 ) is given by There is no change in the entropy of the working fluid after one cycle.Since the change in the entropy of the working fluid occurs only in the nonadiabatic branches ((i) and (iii)), which yields Here, it should be noticed that no heat leakage or entropy production is considered in this model.

Case I: Finite Time Studies with Infinite Reservoirs
In this section, we will first optimize the work per cycle over the initial temperature T(0) of the working fluid and then with respect to the switching time t 1 .
Under the infinite reservoirs condition (C 1 → ∞, C 2 → ∞), the temperatures of hot and the cold reservoirs remain fixed at values T 1 (0) and T 2 (0), respectively.Using ( 9), (14), and ( 22) and substituting λ = K 1 /C and β = K 2 /C, we can calculate temperature of working fluid after the adiabatic step as Substituting the above expression for T + (t 1 ) in ( 19), the work obtained is given by where for convenience, we have defined On using (20) with the infinite reservoir condition, and (24), the efficiency of the heat engine (η = W/Q 1 ) is given by Now we will optimize the work per cycle with respect to initial temperature of the working medium, for a fixed switching time, and later, optimise with respect to the switching time also Solving the above equation, we observe that at maximum work, Substituting the optimum values of α and T(0) in (24), we obtain 1 − e −(K1/C)t1−(K2/C)(t2−t1) . (30) Using ( 26) and ( 28), the efficiency at optimum work is found to be Note that the efficiency is independent of the switching time t 1 .
We can further maximize the work with respect to the switching time: Solving the above equation, we get Though we have not obtained an explicit formula for t 1 , (33) can be solved numerically to get the optimum t 1 .

The Isothermal Limit:
When we take the nonadiabatic paths to be isothermal, that is, C → ∞, we obtain the following results for the maximum work: The temperature of the working fluid is equal to that given by (34) from time, t = 0 to t = t 1 and to that given by (35) from time t = t 1 to t = t 2 .Further, with the isothermal condition (C → ∞), we can simplify (33) to calculate the optimum switching time explicitly as Using this switching time and the isothermal condition, the maximum work that can be extracted from the heat engine is calculated to be The above results are exactly those obtained for an endoreversible cycle [1].We observe that in the isothermal limit, our heat engine becomes identical to the Curzon-Ahlborn cycle.

Case II: Finite Time Studies with Finite Reservoirs
In this section, we carry out the optimisation of work under the constraints of both finite reservoirs and finite cycle time.Solving similarly as in Section 3, the work becomes optimal when the initial temperature of the fluid is given by The optimal work is explicitly given as The impact of finiteness of time can be seen on the work obtained, in Figure 2. It is shown that when time is large enough, the maximum work obtained is approximately equal to that in infinite time case (43).The efficiency at optimal work is again found to be CA value.To the best of our knowledge, only for Carnot-like cycles with isothermal branches (C → ∞) [8], it has been shown that finiteness of the reservoirs does not affect the efficiency at maximum power.Here, we have seen that the generic heat cycle with two adiabatic and two polytropic branches obtains maximum power at the CA efficiency.This efficiency is still unaffected by the switching time, even with the finite reservoirs.We may further maximize the work with respect to the switching time: The optimal switching time is the solution to the following equation: where λ and β are given by ( 8) and (13).The impact of finiteness of reservoirs on the optimal switching time can be seen in Figure 3.The optimal switching time with infinite reservoirs (but finite C for working fluid) for the parameters in Figure 3 is approximately 0.792, whereas the corresponding Curzon-Ahlborn expression, (37) (infinite reservoirs and infinite C working fluid), will yield about 0.828.In short, from numerical results we observe that optimal switching time varies significantly from the above limiting values if the reservoirs and working fluid are finite.
4.1.Infinite Time Limit.Further, if in (33), we keep the constraints of finite reservoirs, but let the nonadiabatic branches proceed very slowly (t 1 → ∞ and (t 2 − t 1 ) → ∞), we get the following expression for the maximum possible work that can be extracted from the heat engine: The factor 1/(1 + (C/C 1 ) + (C/C 2 )) accounts for finiteness of the reservoirs but with an infinite cycle time.In the limit of both infinite reservoirs as well as infinitely slow cycles, the above expression for the maximum work goes back to the one derived by Leff [11] (particularly Otto cycle and Joule-Brayton cycle), where C is the heat capacity of the working fluid.

Summary
In this paper, the maximum power point characteristics of a generalized four-step heat engine are studied.The heat cycle has two adiabatic steps and in the remaining two steps the working fluid follows a polytropic process with a constant heat capacity.It is observed that some common heat cycles such as Otto cycle, Joule-Brayton cycle, and Carnot cycle can be incorporated as special cases of this model.Curzon and Ahlborn in their seminal paper derived the expression for efficiency at maximum power for a Carnotlike engine with infinite reservoirs.Later, Gordon proved that even with finite reservoirs, Carnot-like engines show this efficiency at maximum power.Here we have analysed this result for the generalised heat cycle and shown that CA efficiency is obtained by optimizing the initial temperature of the working fluid and is independent of the switching time.However, we can further maximize this work over the switching time.Finally, the effect of finiteness of cycle time alone, the finiteness of the reservoirs alone, and the finiteness of both the cycle time and the reservoirs together on the maximum work per cycle can be respectively attributed to the following factors: 1 − e −(K1/C)t1 1 − e −(K2/C)(t2−t1) 1 − e −(K1/C)t1−(K2/C)(t2−t1) , ( It can be easily shown (see the appendix) that all these factors are less than unity, and thus, the irreversibilities due to the finiteness of time and/or the reservoirs actually reduce the maximum work that can be performed by these model engines, over the infinite-time and infinite-reservoir models.
Heat lost to cold reservoir, with process with specific heat, C contraction Temperature = T + (t 1 ) Temperature = T − (t 2 ) Temperature = T(0)

Figure 1 :
Figure 1: Heat cycle of the working fluid.