Unification perspective of finite physical dimensions thermodynamics and finite speed thermodynamics

The paper studies the possibility of unifying the two branches of the irreversible engineering thermodynamics, namely finite physical dimensions thermodynamics (FPDT) and finite speed thermodynamics (FST), aiming to take into account their benefits and successes and to eliminate as much as possible their disadvantages. Actually, the two branches have the same goal, that of optimizing the performance of thermal machines and they were developed almost in parallel. Analysis of thermal machines cycles using the FPDT is based on the first and second law of thermodynamics, in the presence of the external irreversibility generated by the heat transfer at finite temperature difference at the thermal reservoirs and internal irreversibility, using the internal source of entropy considered as parameter or function to be specified. The FST is based on the mathematical expression of the first law for process with finite speed that involves three causes of internal irreversibility, namely the finite speed of the piston, internal friction and throttling. The direct method is used in the analysis of thermal machines cycles to provide analytical expression of the machine performance (efficiency and power) as a function of the speed of the process. The significant progress of these two branches of irreversible engineering thermodynamics makes their unification a desirable outcome. We hope that the new model yielded from this study will provide an even more important tool for engineers that will help their attempt to a better design and optimization of thermal machines.

Abstract The paper studies the possibility of unifying the two branches of the irreversible engineering thermodynamics, namely finite physical dimensions thermodynamics (FPDT) and finite speed thermodynamics (FST), aiming to take into account their benefits and successes and to eliminate as much as possible their disadvantages. Actually, the two branches have the same goal, that of optimizing the performance of thermal machines and they were developed almost in parallel. Analysis of thermal machines cycles using the FPDT is based on the first and second law of thermodynamics, in the presence of the external irreversibility generated by the heat transfer at finite temperature difference at the thermal reservoirs and internal irreversibility, using the internal source of entropy considered as parameter or function to be specified. The FST is based on the mathematical expression of the first law for process with finite speed that involves three causes of internal irreversibility, namely the finite speed of the piston, internal friction and throttling. The direct method is used in the analysis of thermal machines cycles to provide analytical expression of the machine performance (efficiency and power) as a function of the speed of the process. The significant progress of these two branches of irreversible engineering thermodynamics makes their unification a desirable outcome. We hope that the new model yielded from this study will provide an even more important tool for engineers that will help their attempt to a better design and optimization of thermal machines.
Keywords Finite physical dimensions thermodynamics Á Finite speed thermodynamics Á Internal and external irreversibility Á Source of internal entropy Á Thermal machines

Introduction
Two branches of the engineering irreversible thermodynamics, namely finite physical dimensions thermodynamics (FPDT) and finite speed thermodynamics (FST), have been developed almost in parallel in recent years. They have the same goal, to model and optimize the performance of thermal machines. The analysis and optimization of thermal machine cycles by FPDT are based on the first and second law of thermodynamics, where the external irreversibility generated by the heat transfer at finite temperature difference is considered, together with internal irreversibility introduced by the internal source of entropy. The finite dimension considered in the models could be thermal conductance, heat transfer area in the machine's heat exchangers, or cost of heat transfer area. Several models have been developed for direct and reverse cycle thermal machines as Stirling, Carnot, Otto, Diesel, Ericsson, Brayton-Joule, etc. [1][2][3][4][5][6][7].
FST is based on the expression of combined first and second law of thermodynamics for processes with finite speed that involves internal irreversibilities due to finite piston speed, internal friction and throttling of the gas in valves, and also external irreversibilities due to heat transfer at finite temperature difference at the source and the sink. Analysis of thermal machine cycles is made by the direct method that is able to derive analytical relationships of machine performance (efficiency/COP and power) depending on the speed of the process and other variables [8][9][10][11][12][13][14]. Also, the entropy generation in the cycle and system can be expressed as a function of speed and many constructive and operational parameters, such as temperature of the thermal reservoir, the compression ratio and the dimensions.
The significant progress of these two branches of engineering irreversible thermodynamics justifies an investigation into the feasibility of their unification. The present paper is aimed at establishing the framework of this unification tentative, highlighting the contribution of each branch.
Analyzing the number of variables and the relations between them, we find that the unification between FPDT and FST is not possible through a single ''connection point'' (the internal entropy source), but it has to be done through two ''connection points'' (the two heat fluxes, which have to be equal between the two theories).

The model
As in Ref. [4], either other processes inside the engine are adiabatic (Carnot engine) or the heat exchanges are entirely recuperated inside the engine (Stirling engine with perfect thermal regeneration).
The following analysis depends crucially on the number of parameters, the number of variables and the number of relations, so it is important to enumerate them clearly. The quantities of interest are given as follows: 1. The heat flux _ Q H transferred from the source, 2. The heat flux _ Q L transferred to the sink, 3. The mechanical power _ W, 4. The efficiency g ¼ _ W= _ Q H , 5. The internal entropy source _ S M , 6. The total entropy source _ S T , and 7. The total cost C T (the cumulated cost of the heat exchangers H and L).
Besides these, we also consider four absolute temperatures: T SH and T SL for the sources (constant parameters) and T H and T L for the internal system (variables). FST introduces another variable, the piston's average speed w.
The quantities should satisfy these inequalities: The first law of thermodynamics The systems H and L do not exchange mechanical work and do not change state, so for them the first law is simply _ Q in ¼ _ Q out . For the system M (the engine), the first law is given as follows: We use the general sign convention: the work _ W done by the system is positive, the received heat _ Q H is positive and the released heat _ Q L is negative.

The second law of thermodynamics
We use the same convention for entropy: the entropy going into the system is positive, the entropy going out of the system is negative. We define the following entropy sources: for the system H: for the system M: for the system L: In the above equations _ S M is the internal entropy source-it represents the entropy generated due to all the irreversibilities inside the engine because of finite speed (friction, throttling, spatial non-uniformity of pressure, etc.). _ S H and _ S L are the external entropy sources, associated with the heat transfer in finite time and through finite area surfaces. _ S T is the total entropy source of the whole system, which accounts for all the irreversibilities.
For the entropy sources we obtain immediately the formulae: The cost We consider that the cost of a heat exchanger is a function of: the heat flux _ Q that the heat exchanger must provide, the temperature T that the heat exchanger must provide, and the temperature T S of the source to which the heat exchanger is connected. In fact, we consider that the cost is directly proportional to the heat flux: This hypothesis is compatible with any of the heat transfer laws considered in [4]: where K is a constant.
Remark The cost might be a financial cost and could also be any other parameter that is additive and should be kept constant, e.g., it could be the area of the heat exchanger. This means that the present approach unifies all the conditions that before were treated separately: constant thermal conductance, constant heat transfer area, or constant cost of heat transfer area. Any of these can be treated in the same manner, just by changing the function f.
Considering that the total cost C T regards only the heat exchangers H and L, we have: Remark When the temperature is higher than that of the source (as it is the case for T L ), the heat flux changes its sign; for the cost to remain positive, the function f also needs to change its sign. Therefore, f L (T L ) is always negative.

Degrees of freedom
The system is completely characterized by 11 variables: seven quantities of interest plus four temperatures. The sources' temperatures are constant, so nine variable quantities remain ( Fig. 1).
There are five relations linking them: the power formula (the first law), the definition of efficiency, the definition of the internal entropy source (the second law for the internal system), the definition of the total entropy source (the second law for the total system) and the definition of the cost. This means that each quantity can be expressed as a function of four others (usually the two internal temperatures plus two other quantities).
To these relations one can add relations concerning the heat transfer (which relate heat fluxes to temperatures). The total system In Fig. 2, we show the dependencies between quantities of interest for the total system. In the total system we have only five variable quantities (the internal temperatures, the internal entropy source and the total cost vanish) and three relations (the definitions of the internal entropy source and of the total cost vanish). It follows that any quantity can be expressed as a function of other two. We choose one of these two to be the total entropy source _ S T . Independently of other conditions, from the relations valid for the total system follow the dependencies on the total entropy source shown in Table 1.
These formulae allow one to study the influence of the total entropy source when one of the quantities is held constant. For instance, in the first row one can see that at constant g the power _ W is directly proportional with _ S T , both reaching synchronized maximum/minimum values. In the second row, one can see that at constant _ W the efficiency g decreases with _ S T , reaching a maximum when _ S T is minimum and vice versa. How each quantity varies with _ S T is synthesized in Table 2.
Remark In this analysis, the internal temperatures were not taken into account: all the above results regarding the maximization/minimization of some quantities in relation to the total entropy source are valid generally, regardless of any other restrictions and of their influence on the internal temperatures. Particularly, it is true that at constant g, the maxima/minima of _ W is synchronized with those of _ S T ; while at constant _ W, g is maximum when _ S T is minimum (and vice versa)-regardless of other restrictions (such as a cost restriction, for example) and regardless of the heat exchange laws. Figure 3 illustrates the dependencies between the quantities of interest and the variables (in black circles) for the internal system. The novelties are the internal temperatures, together with the quantities depending on them: the internal entropy source _ S M and the total cost C T . On the diagram there is no direct connection between the entropy sources; in fact, this connection exists, but just as a consequence of the fact that they both depend on the heat fluxes. We cannot eliminate the heat fluxes between them, but we can obtain a relation that contains only their ratio (expressed through the efficiency g):

The internal system
In this formula one can verify that if the internal entropy source is zero, then the engine efficiency has to be the Carnot efficiency between the internal temperatures. And if the total entropy source is zero, then the engine efficiency has to be the Carnot efficiency between the sources' temperatures.
Besides the efficiency, the formula (10) contains both the internal temperatures. As a consequence, all the relations regarding the internal entropy source also contain them. With the help of formula (10) we get the formulae relating _ S M to other quantities when keeping one parameter constant, presented in Table 3.
We can also write formulae for the total cost (they would make an additional column, which do not fit in the table): at _ W ¼ const:: at _ Q L ¼ const:: From these we obtain the formulae for constant total cost (they would have made an additional row in the table, but it does not fit there): The presence of the internal temperatures in all the above formulae prevents us from studying the quantities' variation with the internal entropy source _ S M , as it was possible for the total system. It appears that it would be necessary to know how the internal entropy source _ S M is related to the internal temperatures-not through the heat fluxes (the definition of _ S M , which we already have), but independently of them. This is the motivation of the present study: researching the possibility that FST and FPDT can provide a formula for the internal entropy source _ S M -a formula dependent on the internal temperatures and on other variables or parameters, but independent of the heat fluxes.  Table 3 Internal system's quantities dependence on the internal entropy source

Finite speed thermodynamics
Approaching the problem from the inside of the engine, FST focused on the internal irreversibilities: spatial non-uniformity of pressure, various frictions and throttling in valves. All these phenomena are related to the average speed of the engine's piston. Evaluating these losses for each transformation and adding them, FST succeeds in providing an analytical formula for the engine's power and efficiency, both as functions of the piston's average speed. In [15], an analytical formula for the internal entropy source is provided, which seems to be exactly what we are looking for. The calculation scheme is in each step, the internal temperatures T H and T L are considered given, and based on them, on the speed w and on other constructive parameters, FST computes the heat; then, it computes the cycle duration and thus finds the heat fluxes _ Q H and _ Q L ; from them, the power _ W and the efficiency g follows immediately (Fig. 4). Knowing the heats and the temperatures, the entropy source _ S M follow from its definition. Coupling the two diagrams (Figs. 3, 4), we get Fig. 5. Now we can compute the entropy sources, based on the heat fluxes provided by FST. For each triple of values given to the variables (w, T H , T L ) we get the heat fluxes, the power, the efficiency, the two entropy sources and the total cost. The determining factor is not one (the internal entropy source), but two: the heat fluxes-which, once computed, determine uniquely all the other quantities of interest.
FST is not only able to provide the internal entropy source, but in the process also determines uniquely the heat fluxes, which cannot be varied any more by FPDT. We introduced a new variable (w) and two relations (the heat fluxes provided by FST). If before the unification we could express each quantity as a function of four quantities (two temperatures and two other quantities- Table 3), now we can express each quantity as a function of three variables: the speed w and the two internal temperatures. But this was anyway done by FST alone, so we gained nothing-we did nothing more than to ''unblock'' the temperatures, which were constant in FST and now are variable.
Unfortunately, this is not enough for optimization. Because there will be no optimum point-in any regime, if the internal temperatures get closer to the external ones we will always get more power. FST does find an optimum speed w at which the power is maximum, but this result is relative to the internal temperatures-increasing their ratio even by a tiny bit will increase the power, so that we will always discover that ''we can do better''. In other words, the absolute maximum of power will be that calculated by FST when the internal temperatures equal the sources' temperatures.
To be able to optimize something, we have to introduce either some more relations or a restriction. They have to depend on the internal temperatures, but in such a way that they introduce a new relation, not already present. There are two options, both based on heat transfer: (a) We add specific relations for the heat transfer (thus linking heat fluxes to temperatures through known, fixed conductances), or (b) We impose a fixed total cost of the heat exchangers (the cost being based on heat transfer laws, but with unknown conductances).
Both these options lead to a problem that can be solved through the Lagrange multipliers method: find the extrema of some objective function O with the restrictions stated in

Heat transfer laws
Unfortunately, in our case it is not simple to add the heat transfer laws. Because we already have connections between the temperatures and the heat fluxes: those given by FST. If we add formulae for heat transfer, we get the situation in Fig. 6. The arrows marked with ''!'' show a double determination: starting from the same temperatures, FST determines some heat fluxes, and the heat transfer laws determine other heat fluxes.
We can overcome the difficulty by considering the two additional relations as restrictions and using the Lagrange multipliers method.
Let us consider that the heat transfer is described by one of the following formulae: where the constant K Q encapsulates the properties and the area of the heat exchange surface. In general we can write Denoting by O the chosen objective function, we write the Lagrange function: We nullify the partial derivatives: From the first two equations we find the Lagrange multipliers, which we introduce in the third equation. We eventually get a system of three equations with three unknowns: The last two equations are the initial heat equalities. One cannot solve them analytically, because the temperatures are too deeply buried into the formulae-they show up simultaneously in exponentials, in logarithms and under the radical (see Sect. 6).
The system can be solved numerically by choosing a suitable algorithm. Now let us explore solution (b), choosing the total cost as a constant parameter.

Constant total cost
We introduce a cost function as shown in Fig. 7.
We consider that the cost of a heat exchanger is a function of the constant K Q from the heat transfer law: ''Cost'' might not only mean a financial cost, but could also mean size, or any other parameter which we wish to maintain constant.
We consider that the total cost consists solely of the costs of the heat exchangers H and L: We want to optimize an objective function O with the constant cost restriction, so we employ again the Lagrange multipliers method. The Lagrange function is B is the budget available for the heat exchangers. We nullify the partial derivatives: From the first equation, we express k and then plug it in into the other two. We obtain a system of three equations with three unknowns: . . . partial derivatives equation. . . . . . partial derivatives equation. . .
The last equation is the constant cost condition. It cannot be solved analytically, again because of its complexity. The system can be solved numerically, if we choose a suitable algorithm.
The irreversible heat fluxes and the entropy source provided by FST From Ref. [15] we have the formula of the internal entropy source for an irreversible Carnot engine, which was determined based on the heat fluxes. Adapting the formulae in [15] to the present notations and adding some new notations, the heat fluxes _ Q H;FST and _ Q L;FST depend on the average piston speed w, on the internal temperatures T H and T L , and on the following seven parameters:The formulae for the heat fluxes are where c 23 ¼ 1 þ ðc À 1ÞB ad;2À3ir ; ð30Þ a ¼ 1 and the numerical coefficients B are The maximum pressure p 1 is and the denominator p 3 (1 ? e L ) in the last formula is Remark In [15], this pressure is considered reversible (as if c 23 = c 41 = c); but this approximation is not necessary, since at this point we already know B ad,2-3ir and B ad,4ir-1 , so we know c 23 and c 41 .
In the expressions of B coefficients, y is the stroke of the respective transformation (isothermal or adiabatic). Taking all these strokes equal to the piston's semi-stroke, we approximate it reversibly: The expressions (28) and (29) of the irreversible heat fluxes are needed for optimizing the speed w together with the temperatures T H and T L . In those formulae, the temperatures appear simultaneously in exponentials, in logarithms and under the radical (and combined). It seems clear that they are too complicated to give an analytical solution and the resulting systems of equations need to be treated numerically.
Applying the definition of the internal entropy source to the formulae for the irreversible heat fluxes _ Q H;FST and _ Q L;FST , in [15] the following synthetic formula is obtained: This formula, regarded by itself, hides the fact that the internal entropy source was calculated based on the heat fluxes and creates the illusion that it could be used with arbitrary heat fluxes. But doing so would be a mistake: the formula (40) is valid only for the heat fluxes _ Q H;FST and _ Q L;FST used to calculate it and is not valid for some arbitrary heat fluxes _ Q H and _ Q L that may happen to have the same weighted sum, _ Q H =T H þ _ Q L =T L . When we write the formula (40) the heat fluxes are already determined, as well as all the other quantities in the system (power, efficiency, etc.).

Conclusions
FST can provide to FPDT a formula for the internal entropy source of the irreversible Carnot engine, but only through the heat fluxes. This means that the heat fluxes cannot be varied independently in FPDT: they are determined by FST and they determine all the other quantities. Analyzing the number of quantities, variables and relations between them, we conclude that the unification between FPDT and FST is not possible through a single connection (the relation for the internal entropy source provided by FST), but it needs to be done through two connections (the two relations for the heat fluxes provided by FST).
In this way, the average speed of the piston and the internal temperatures can be optimized simultaneously either by imposing heat transfer laws (conductive, convective, radiative, or a combination) or by maintaining the total cost constant. Even if obtaining analytical formulae is not possible, the resulting systems of equations can be solved using the Lagrange multipliers method and a suitable numerical algorithm.
Some of the studies at one constant parameter can be performed in a simpler way, independently of irreversibilities and of heat transfer laws. For instance, using only The First Law and The Second Law of Thermodynamics (and no other constraints, not even the heat transfer laws), we showed that: -at constant efficiency, the regime that maximizes the power of the irreversible Carnot engine is the one that also maximizes the total production of entropy, -at constant power, the regime that maximizes the efficiency of the irreversible Carnot engine is the one that minimizes the total production of entropy.