1. Introduction
During more than four decades, a branch of irreversible thermodynamics called finite-time thermodynamics (FTT) has been developed [
1,
2,
3,
4,
5]. Within the context of FTT, it has been possible to elaborate detailed models of heat engines working at finite-time and with entropy production [
6,
7,
8,
9,
10]. One of the applications of FTT has been in the field of thermoeconomics [
11,
12,
13,
14]. In 1995, Alexis De Vos introduced for the first time a thermoeconomic study for the Novikov plant model [
15]. In his study, De Vos takes into account two costs: a cost related to the investment and another cost associated to the fuel consumption. From his thermoeconomic study, De Vos found that the optimal thermoeconomic efficiency under maximum power conditions is between the Curzon–Ahlborn (CA) efficiency and Carnot efficiency [
16]. Here, De Vos refers to the power production rate
of the Novikov plant model. Later, Barranco-Jiménez and Angulo-Brown [
17,
18] extended the De Vos thermoeconomic approach, but using the so-called ecological function [
19,
20], and they also showed that the economical efficiency under this regime of performance lies between the CA efficiency and the Carnot efficiency. In 2009, Barranco-Jiménez [
21] studied the thermoeconomics of a nonendoreversible Novikov engine. He used different heat transfer laws: The Newton’s law of cooling, the Stefan–Boltzmann radiation law [
22], the Dulong–Petit’s law of cooling [
23] and a phenomenological heat transfer law [
24,
25] (the usage of each transfer law depends on the kind of heat engine model to describe [
22,
23,
24,
25,
26,
27]). He calculated the thermoeconomic optimal efficiencies under two regimes of performance, namely, the maximum power regime and the ecological regime, and he found that when the Novikov model maximizes the ecological function, it reduces the rejected heat to the environment down to about 50% with respect to the heat rejected under maximum power conditions. Recently, Pacheco-Paez et al. [
28,
29] made a thermoeconomic analysis of a Novikov plant following the De Vos approach (with a linear heat transfer law) but including a new cost associated to the maintenance of the power plant. The aim of this paper is to extend the thermoeconomic analysis of the Novikov plant model (see
Figure 1) by considering different heat transfer laws. Besides, in the optimization of the benefit functions, we take into account a new cost associated to the maintenance of the power plant [
28,
29]. In our study, we also consider three regimes of performance: maximum power (MP) [
3], maximum efficient power (EP) [
30,
31] and maximum ecological function (ME) [
19,
20]. The methodology used in this paper is within the context of finite-time thermodynamics. Apart from the De Vos method, there exist other approaches for the thermoeconomic analysis of power plant models. Interestingly, the thermoeconomic analysis based on exergy has played an important role in this discipline [
32,
33,
34,
35,
36]. The work is organized as follows: In
Section 2, we present the thermoeconomic analysis of the Novikov power plant under different criteria of performance. In
Section 3, we analyze the environmental impact of the thermoeconomic model, and finally in
Section 4, we present our conclusions.
2. Thermoeconomic Analysis
Applying the first law of thermodynamics to
Figure 1, the power output is given by,
where
and
are the heat transfers between the thermal baths and the working fluid (
,
and
are quantities per cycle’s period). On the other hand, the internal efficiency of the engine is given by [
21],
where
,
and the parameter
is the non-endoreversibility parameter [
37,
38] (which characterizes the degree of internal irreversibility that comes from the Clausius inequality),
being the change of the internal entropy along the hot isothermal branch and
the entropy change corresponding to the cold isothermal compression. This parameter in principle is within the interval
(
for the endoreversible limit). If we consider that the heat transferred between the hot source and the working fluid obeys a generalized law of the type,
where the exponent
can be any real number different from zero and it takes the values
for a Newtonian heat transfer law (N) and
for a Dulong–Petit heat transfer law (DP) [
23,
27] and
is the thermal conductance (see
Figure 1). From Equations (2) and (3), the power output
results in
The De Vos thermoeconomical analysis considers a profit function
, which is maximized [
11]. This profit function is given by the quotient of the power output
(that is, the power production rate
of the Novikov plant model) and the total costs involved in the performance of the power plant, that is
Apart from this De Vos method, there exist other approaches for the thermoeconomic analysis of power plant models. Significantly, the thermoeconomic analysis based on exergy has played an important role in this matter [
32,
33,
34,
35,
36]. In his study, De Vos assumed that the running costs of the plant consist of two parts: a capital cost which is proportional to the investment and, therefore, to the size of the plant, and a fuel cost that is proportional to the fuel consumption and, therefore, to the heat input rate. In this work, we also take into account a cost associated to the maintenance of the power plant which is taken as proportional to the power output of the plant. This choice is based on an intuitive idea that has to do with the wheatering of the engine depending on its excessive usage. Therefore, the running costs of the plant exploitation are now defined as,
where the proportionality constants
a,
b and
c have units of $/Watt, and
is the maximum heat that can be extracted from the heat reservoir without supplying work in the same manner previously considered by De Vos [
11] (see
Figure 1). Equation (6) is taken as a linear relation inspired in the linear character of the first law of thermodynamics [
39]. Using Equations (3) and (4), the running costs of the plant are given by,
where
and
(
and
). After substituting Equations (4) and (7) into Equation (5), the dimensionless objective profit function is given by,
which is valid for the transfer laws of Newton (
), and Dulong–Petit (
). Now, we define two additional objective functions, one in terms of the so-called modified ecological function (
with
the rate of total entropy production of the Novikov model and
a parameter which depends on the heat transfer law used in the Novikov model [
19,
20]). This function is a slight modification of the former ecological function [
19], which preserves the original physical meaning of that function; that is, it represents a good trade-off between high power output and low dissipation. On the other hand, we also use a profit function based on the efficient power [
30,
31]
, which is an objective function of the class of compromise functions looking for a good trade-off between power output and efficiency. Both objective functions are divided by the total costs involved in the performance of the plant, that is,
and
, respectively. Analogously to Equation (8), by using Equations (4), (5) and (7),
and
are given by,
and,
The maximization of the objective functions given by Equations (8)–(10) by means of
gives the corresponding optimal efficiencies. Therefore, taking the derivatives of
,
and
with respect to
and setting them equal to zero, we obtain the optimal efficiency
that maximizes Equations (8)–(10), respectively. In
Figure 2, for the case of a Newton heat transfer law (
in Equation (3)), we show the behavior of the three objective functions versus the internal efficiency
for the following arbitrary values
(endoreversible case),
and
. We observe in
Figure 2 that there exists an efficiency value
which depends on the parameter
. The benefits function diminishes as the parameter
increases. Besides, for the case of maximum power conditions, the maximum value of
is obtained for practically the same value of
for the two values of
considered.
In
Figure 3, we can observe for the three benefit functions how the optimal efficiency tends to the reversible value (Carnot efficiency) when
, such as it occurs in the De Vos approach [
11,
21].
The optimal efficiencies are obtained in terms of the dimensionless parameter
, but this parameter is not easy to evaluate [
11] and it is more convenient to work in terms of the so-called fuel fractional cost (
) [
11,
17], defined as the ratio between the fuel cost and the total costs of the plant, that is,
which is in the interval
. From the previous equation and by using Equations (2)–(4), we get
From Equations (8)–(10), we can numerically find the optimal efficiencies by solving polynomial equations of the form:
. However, some analytical expressions can be obtained for
in the cases of both Newton (
) and Dulong–Petit (
) laws previously reported in [
21]. For the case of maximum power output, the optimal efficiencies are given by [
20]
As it was previously mentioned, under the maximum power conditions, the optimum value of
does not appreciably change the maximum of the benefit function, that is, the optimal efficiency is practically independent of the parameter
, which is associated to the maintenance costs of the plant (see
Figure 2). Interestingly, the inclusion of this new parameter has a remarkable effect on the optimal efficiencies for the cases of ME and EP maximizations. In fact,
and
increase with respect to the results without maintenance costs (see
Figure 2); that is, better maintenance is reflected as better efficiency. For the case of maximum efficient power and maximum ecological function (with
), the optimal efficiencies are given by,
where
. For the cases
and
, from Equations (13a), (14a) and (15a), we obtain the optimal efficiencies
,
and
previously reported by Curzon-Ahlborn [
16], Yilmaz [
30] and Arias-Hernández and Angulo-Brown [
20], for the Curzon-Ahlborn heat engine under maximum power output, maximum efficient power and maximum ecological function conditions, respectively.
Now, we consider two other heat transfer laws between the hot source and the working fluid as follows,
Here, the exponent
takes the values of
for a phenomenological heat transfer law (
) [
24,
25] and
for a Stefan–Boltzmann (SB) radiation law ([
22,
26]) and the Sign function leads to
if
, and
if
. In this case, the objective functions under the three regimes considered result in,
and,
In
Figure 4 and
Figure 5, we show the behavior of
,
and
for the case of a phenomenological (
) heat transfer law (
in Equations (18)–(20)) versus the internal efficiency for
,
and
(a similar behavior is found for the benefit functions for the case of the Stefan–Boltzmann law). We can see in
Figure 5 and
Figure 6 that these profit functions have a similar behavior to the Newton and the Dulong–Petit cases; that is, the benefits function also diminishes as the parameter
increases. In [
11], it was shown that the maxima points of the MP profit function tend to the reversible limit (Carnot point) when the parameter
tends to infinite. This same behavior occurs for the EP and ME profit functions. In
Figure 6 we only show, as an example, the case of the EP profit function for several increasing values of the parameter
, and it is clear how the curve that contains the maxima points tend to the Carnot efficiency (
, in this case).
In a similar way to the Newton and the Dulong–Petit cases, we can obtain the optimal efficiencies in terms of the parameter β (or ).
Analogously to Newton and Dulong–Petit laws, we can also obtain some analytical expressions for
for both Ph (
) and SB radiation (
) laws. Therefore, under maximum power output, maximum efficient power and maximum ecological function conditions, we obtain, respectively,
and,
Similarly as for the Newton heat transfer law, for a phenomenological heat transfer law when
and
, from Equation (21a–c), we obtain the optimal efficiencies
previously reported by De Vos [
24],
and
(these last two expressions were previously published by Barranco-Jiménez [
21]), for the Curzon-Ahlborn heat engine under maximum power output, maximum efficient power and maximum ecological function conditions, respectively.
3. Numerical Results: Environmental Impact
When we calculate the derivative of Equations (8)–(10), that is,
, for
, we obtain equations of the form
. Therefore, we have to solve numerically polynomial equations of the form,
These previous expressions are not shown in explicit form because they are very large. Finally, using the expression for the fractional fuel cost
[
11,
17] (see Equation (12)), we obtain the numerical optimal efficiencies in terms of the parameter
, the parameter
and the economic parameters
f and
. In
Figure 7, we show the optimal numerical efficiencies under maximum power, maximum efficient power and maximum ecological function conditions, respectively. In this Figure (see solid lines), we also show the optimal efficiencies given by Equations (13a), (14a) and (15a).
On the other hand, in
Figure 8 we depict the corresponding optimal efficiencies against
f for the Newtonian case, but with
R = 0.9 < 1 (this case is a representative one of any nonendoreversible situation). As we can see in this case, the
optimal values are all smaller than those corresponding to
Figure 7. This behavior for
R < 1 is present for any heat transfer law and any operating regime [
21,
40].
We can also see that the optimal efficiencies under the maximum ecological function regime (for all values of the parameters
,
and
for the case of Newton heat transfer law, see
Figure 7 and
Figure 8) are greater than the optimal efficiencies under both maximum power output and maximum efficient power conditions, respectively. For the case of both phenomenological and Stefan–Boltzmann laws, the optimal efficiencies have a similar behavior with respect to the parameters
,
and
. Besides, all of these optimal efficiencies satisfy the following inequality:
The previous inequality was recently obtained by Barranco-Jiménez et al. for the case of both endoreversible and nonendoreversible models of heat thermal engines [
21,
40]. In addition, we can see in
Figure 2 and
Figure 3 how the optimal efficiency under the maximum ecological function regime leads to a high reduction on the profits in comparison with those at maximum power output. However, this reduction is concomitant with a better efficiency for each value of the parameter
(see Table 1 in [
11] for fractional fuel costs for several kinds of fuels). This value of the efficiency provides a diminution of the energy rejected to the environment by the power plant. Therefore, the ecological criterion becomes an appropriate procedure for the search for insight into an engine’s performance for providing a less aggressive interaction with the environment. Therefore, we propose the following simplified analysis:
From the first law of thermodynamics we get,
where
is the heat rejected to the environment by the power plant. For the cases of heat transfer of the Newton and Dulong–Petit types, when the power plant works under maximum power conditions by using Equations (3), (4) and (25), we get
In a similar way, the expressions when the power plant works under maximum ecological function and maximum efficient power conditions are, respectively,
The solutions of Equations (26)–(28) were obtained in a numerical manner similarly to the equations in
Section 2. In our optimization, we consider different values for the maintenance of the power plant (parameter
) and we can observe for each thermodynamic regime that while the
parameter increases, the heat rejected to the environment slightly diminishes, as it is shown in
Figure 9.
In
Figure 9, we show the amount of heat rejected to the environment for each one of the three regimes of performance. As we can see, the heat rejected to the environment under the maximum ecological function presents a less aggressive impact to the environment in comparison with the other two regimes (see
Figure 10a,b). We can also observe in
Figure 10 that when the cost of the maintenance is taken into account, the rejected heat slightly diminishes in the three regimes of performance considered.
In
Figure 11, we show the ratio between the amounts of heat rejected for a heat transfer law of the Dulong–Petit type.
Furthermore, by applying the second law of thermodynamic (see
Figure 1) we get,
Using Equations (1), (3) and (4) and after some algebraic manipulations, for the cases of heat transfer of the Newton and Dulong–Petit types, the entropy production under maximum power conditions can be calculated by the following expression:
Similarly to Equation (30), the expressions when the power plant works under maximum ecological function and maximum efficient power conditions are, respectively,
It is important to mention that all of the calculations were made in a numerical form. The behavior of the above equations can be observed in
Figure 12. As can be seen in
Figure 12b, the entropy production of the power plant when it is working under maximum ecological conditions is less than when it is working under maximum power conditions.
In
Figure 13a, we show the ratio between the amounts of entropy production for two regimes: maximum ecological function and maximum power. In a similar manner, in
Figure 13b, we show the ratio between the amounts of entropy production for maximum efficient power conditions and maximum power conditions. In addition, in
Figure 14, we show the ratios of entropy production for a heat transfer of the Dulong–Petit type. In all of our numerical calculations we have attempted to use for the distinct involved parameters only typical values reported in the FTT literature.