Open AccessArticle

Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes

Marco A. Barranco-Jiménez

^1,†,

Norma Sánchez-Salas

^2,*,† and

Israel Reyes-Ramírez

^3,†

Departamento de Formación Básica, Escuela Superior de Cómputo, Instituto Politécnico Nacional, Av. Miguel Bernard Esq. Juan de Dios Bátiz, U.P. Zacatenco, Mexico D.F. 07738, Mexico

Escuela Superior de Física y Matemáticas, Instituto Politécnico Nacional, Edif. 9 UP Zacatenco, Mexico D.F 07738, Mexico

Unidad Profesional Interdisciplinaria en Ingeniería y Tecnologías Avanzadas, Instituto Politécnico Nacional, Av. IPN 2580, L. Ticomán, Mexico D.F. 07340, Mexico

Author to whom correspondence should be addressed.

^†

These authors contributed equally to this work.

Entropy 2015, 17(12), 8019-8030; https://doi.org/10.3390/e17127860

Submission received: 18 October 2015 / Revised: 16 November 2015 / Accepted: 24 November 2015 / Published: 4 December 2015

(This article belongs to the Special Issue Exploring the Second Law of Thermodynamics)

Download

Browse Figures

Figure 1
Schematic representation of an endoreversible heat engine. "> Figure 2
Profit functions for <math display="inline"> <mrow> <mi>β</mi> <mo>=</mo> <mn>0</mn> </mrow> </math> and <math display="inline"> <mrow> <mi>β</mi> <mo>=</mo> <mn>5</mn> </mrow> </math> versus η with <math display="inline"> <mrow> <mi>R</mi> <mo>=</mo> <mn>1</mn> </mrow> </math> and <math display="inline"> <mrow> <mi mathvariant="sans-serif">τ</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </math>. "> Figure 3
The steady-state efficiencies working under maximum power output (<math display="inline"> <msub> <mover> <mi>η</mi> <mo>¯</mo> </mover> <mrow> <mi>M</mi> <mi>P</mi> </mrow> </msub> </math>), maximum-efficient power (<math display="inline"> <msub> <mover> <mi>η</mi> <mo>¯</mo> </mover> <mrow> <mi>E</mi> <mi>P</mi> </mrow> </msub> </math>) and maximum ecological function (<math display="inline"> <msub> <mover> <mi>η</mi> <mo>¯</mo> </mover> <mi>E</mi> </msub> </math>) conditions. "> Figure 4
Plot of relaxation times under maximum power conditions versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1). "> Figure 5
Plot of relaxation times under maximum efficient power versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1). "> Figure 6
Plot of relaxation times under maximum ecological function conditions versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1). "> Figure 7
Relaxation times in the endoreversible case (<math display="inline"> <mrow> <mi>R</mi> <mo>=</mo> <mn>1</mn> </mrow> </math>) versus fractional fuel cost for several values of τ for (a) Maximum efficient power conditions and (b) Maximum ecological function. "> Figure 8
Ratio of the relaxation times versus τ for a value of the fractional fuel cost and several values of the parameter R (cases (a), (c) and (e)), and for the endoreversible case <math display="inline"> <mrow> <mi>R</mi> <mo>=</mo> <mn>1</mn> </mrow> </math>, for different values of the fractional fuel cost, (cases (b), (d) and (f)). ">

Versions Notes

Abstract

A recent work reported a local stability analysis of a thermo-economical model of an irreversible heat engine working under maximum power conditions. That work showed that after small perturbations to the working temperatures, the system decreases exponentially to the steady state characterized by two different relaxation times. This work extends the local stability analysis considering other performance regimes: the Maximum Efficient Power (MEP) and the Ecological Function (EF) regimes. The relaxation time was shown under different performance regimes as functions of the temperature ratio

τ = T_{2} / T_{1}

, with

T_{1} > T_{2}

, the fractional fuel cost f and a lumped parameter R related to the internal irreversibilities degree. Under Maximum Efficient Power conditions the relaxation times are less than the relaxation times under both Maximum Ecological function and Maximum Power. At Maximum Power Efficient conditions, the model gives better stability conditions than for the other two regimes.

Keywords:

thermo-economics; local stability; irreversible heat engine

1. Introduction

In 2001, Santillán et al. [1] studied the local stability analysis of an endoreversible Curzon-Ahlborn heat engine [2] operating under maximum power conditions. Later, Guzmán-Vargas et al. [3] investigated the effect of the heat transfer laws and the thermal conductances on the local stability of an endoreversible heat engine. Recently, Páez-Hernández et al. [4], analyzed the local stability of a non-endoreversible Curzon-Ahlborn (CA) cycle taking into account implicity the engine’s time delays operating under maximum power regime. However, the local stability analysis described in these works have not considered economical effects. Within the context of Finite-Time Thermodynamics (FTT), the economical aspects were introduced earlier by De Vos [5] to study the thermo-economic model performance of a Novikov type power plant [6,7]. Later, Sahin and Kodal [8] studied the thermo-economics of an endoreversible heat engine in terms of maximization profit function defined as quotient between the power output and the annual investment and fuel consumption costs. This thermo-economic performance analysis [9] consists of maximizing a benefit function in terms of the power output and the cost involved in the power plant performance. All these thermo-economic studies have shown that the inclusion of costs of performance have an important impact in the trade off cost-benefit of the corresponding power plant models [5,10,11]. Recently, Barranco-Jiménez et al. [12,13] reported a local stability analysis of a thermo-economic model of an irreversible heat engine working under maximum power conditions. In those studies, they used two different heat transfer laws, the Newtonian [12] and the Dulong-Petit [13] ones. In this work, the local stability analysis is extended considering other performance regimes: The Maximum Efficient Power [14,15] and the Ecological Function regime [16,17]. The relaxation time shown under maximum efficient power conditions is less than the relaxation times under both maximum power and maximum ecological function; that is, under maximum efficient power conditions, there is better stability conditions than for the other two regimes. The paper is organized as follows: Section 2 presents the thermo-economic analysis of the irreversible heat engine under different performance criteria. In Section 3, the local stability analysis of the irreversible heat engine is shown. Finally, in Section 4, the conclusions are given.

2. Thermoeconomic Optimization of a Curzon-Ahlborn Engine Model at Different Regimes of Performance

In Figure 1, a schematic diagram of the irreversible heat engine (Curzon-Ahlborn model) is shown. This engine consists of a Carnot-like thermal engine that works in irreversible cycles and exchanges heat with external thermal reservoirs at temperatures

T_{1}

and

T_{2}

(

T_{1} > T_{2}

). In the steady state, the temperatures of isothermals of the Carnot-like cycle are

\bar{x}

and

\bar{y}

; here overbars are used to indicate the corresponding steady-state value. The steady-state heat flows as are shown in Figure 1 are denoted as

\bar{J_{1}}

and

\bar{J_{2}}

, respectively.

Figure 1. Schematic representation of an endoreversible heat engine.

Applying the Clausius theorem and using the fact that the inner Carnot-like engine works in irreversible cycles, the following inequality is obtained,

\frac{\bar{J_{1}}}{\bar{x}} - \frac{\bar{J_{2}}}{\bar{y}} < 0,

(1)

this expression can be transformed into an equality by introducing a parameter R leading to,

\frac{\bar{J_{1}}}{\bar{x}} = R \frac{\bar{J_{2}}}{\bar{y}} .

(2)

The lumped parameter R, which, in principle, is within the interval

0 \leq R \leq 1

(

R = 1

for the endoreversible limit), can be seen as a measure of the departure from the endoreversible regime, it has been used for non endoreversible thermal heat engine models as a way to include the global internal irreversibilities [18,19]. Assuming that the heat flows from

T_{1}

\bar{x}

and from

\bar{y}

T_{2}

are of the Newton type, then

\bar{J_{1}} = α (T_{1} - \bar{x}),

(3)

\bar{J_{2}} = α (\bar{y} - T_{2}),

(4)

where α is the thermal conductance. For simplicity of the calculations, it is assumed that the heat exchanges take place in conductors with the same thermal conductance α; that is, the materials of both conductors are the same. Applying the first and second laws of thermodynamics, the system’s steady-state power output and the efficiency can be written as,

\bar{P} = \bar{J_{1}} - \bar{J_{2}},

(5)

and

\bar{η} = \frac{\bar{P}}{\bar{J_{1}}} = 1 - \frac{\bar{J_{2}}}{\bar{J_{1}}} = 1 - \frac{1}{R} \frac{\bar{y}}{\bar{x}} .

(6)

By combining Equations (2)–(4) and (6), gives the steady-state temperatures

\bar{x}

and

\bar{y}

and the power output in terms of

T_{1}

T_{2}

, R and

\bar{η}

as [1,3],

\bar{x} = \frac{T_{1}}{1 + R} (1 + \frac{τ}{1 - \bar{η}}),

(7)

\bar{y} = \frac{R}{1 + R} T_{1} (1 + \frac{τ}{1 - \bar{η}}) (1 - \bar{η}),

(8)

\bar{P} = \frac{α}{1 + R} T_{1} \bar{η} (R - \frac{τ}{1 - \bar{η}}),

(9)

where

τ = T_{2} / T_{1}

. The De Vos thermoecomical analysis considers a profit function F, which is maximized [5]. This profit function is given by the quotient of the power output (

\bar{P}

) and the total cost involved in the performance of the power plant (

{\bar{C}}_{t o t}

), that is,

\bar{F} = \frac{\bar{P}}{{\bar{C}}_{t o t}} .

(10)

In his early study, De Vos assumed that the running cost of the plant consists 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

\bar{J_{1}}

. Assuming that

{\bar{J}}_{m a x}

is an appropriate measure for the size of the plant, the running costs of the plant exploitation are defined as [5],

{\bar{C}}_{t o t} = a {\bar{J}}_{m a x} + b {\bar{J}}_{1} = a α T_{1} [(1 - τ) + β (1 - \frac{\bar{x}}{T_{1}})]

(11)

where the proportionality constants a and b have units of

$ / J o u l e

β = b / a

and

{\bar{J}}_{m a x} = α (T_{1} - T_{2})

is the maximum heat that can be extracted from the heat reservoir without supplying work (see Figure 1). By using Equations (3), (6), (7), (10) and (11), the profit function can be written as,

a {\bar{F}}_{M P} = \frac{\bar{η} [\frac{1}{1 + R} (R - \frac{τ}{1 - \bar{η}})]}{(1 - τ) + \frac{β}{1 + R} (R - \frac{τ}{1 - \bar{η}})} .

(12)

If we calculate the derivative of

a \bar{F}

with respect to

\bar{η}

and we solve for the efficiency

{\frac{d (a {\bar{F}}_{M P})}{d \bar{η}}|}_{\bar{η} = {\bar{η}}^{*}} = 0

, we obtain

\bar{η} = {\bar{η}}^{*} (β, τ, R)

[20]. However, instead of expressing the optimal efficiency in terms of the parameter β, a number that is difficult to obtain in the literature [5], we can also express it in terms of the fractional fuel cost, which is defined as [5],

f = \frac{β {\bar{J}}_{1}}{{\bar{J}}_{m a x} + β {\bar{J}}_{1}} .

(13)

The fractional fuel costs (f) for various technologies were reported by De Vos for different energy sources; that is, for example; renewable energy

f = 0

, for Uranium

f = 0.25

, for Coal

f = 0.35

, and for natural gas

f = 0.5

[5]. By using Equations (3), (7) and (13), the parameter β in terms of the fractional fuel cost can be written as,

β = \frac{f}{1 - f} \frac{(1 + R) (1 - τ)}{R - \frac{τ}{1 - \bar{η}}} .

(14)

Therefore, the efficiency that maximizes the profit function is given by [20],

{\bar{η}}_{M P} (f, τ, R) = 1 - \frac{f}{2 R} τ - \frac{\sqrt{4 (1 - f) τ R + f^{2} τ^{2}}}{2 R} .

(15)

Equation (15) represents the optimal steady-state efficiency (

\bar{η}

) as a function of τ, f and R for a non endoreversible Novikov-Curzon-Alhborn heat engine working in the maximum-power regime.

Analogously to Equation (10), for our thermo-economic optimization approach, two objective functions, the so-called Efficient Power [14,15] and the so-called Ecological Function [16,17] are defined, both divided by the total cost. The Maximum Efficient Power performance [14,15] for heat engines was studied for Yilmaz [14] and previously defined by Stucki [21] in 1980 as the product of power output (P) by the efficiency (η) in the context of the first order irreversible thermodynamics. The ecological optimization criterion for the FTT-thermal cycles was proposed by Angulo-Brown [16]. This criterion considers the maximization of a function E which represents a compromise between high power output (P) and low entropy production Σ. The E function is given by,

E = P - T_{2} Σ

where P is the power output of the cycle, Σ the total entropy production per cycle and

T_{2}

is the temperature of the cold reservoir. One of the most important characteristics of a CA engine operating under maximum-E conditions is that it produces around 80% of the maximum power and only 30% of the entropy produced in the maximum power regime [16]. Another interesting property of the maximum-E regime is that the CA-engine’s efficiency in this regime, is given by

η_{E} \approx (η_{C} + η_{C A}) / 2

, where

η_{C}

is the Carnot efficiency and

η_{C A}

the Curzon-Ahlborn efficiency. For our thermo-economical approach, these objective functions are given by

{\bar{F}}_{E P} = \frac{η \bar{P}}{{\bar{C}}_{t o t}}

and

{\bar{F}}_{E} = \frac{\bar{P} - T_{2} Σ}{{\bar{C}}_{t o t}}

, respectively. In the same way that Equation (12) the profit functions can be written as,

a {\bar{F}}_{E P} = \frac{{\bar{η}}^{2} [\frac{1}{1 + R} (R - \frac{τ}{1 - \bar{η}})]}{(1 - τ) + \frac{β}{1 + R} (R - \frac{τ}{1 - \bar{η}})},

(16)

and

a {\bar{F}}_{E} = \frac{(2 \bar{η} + τ - 1) [\frac{1}{1 + R} (R - \frac{τ}{1 - \bar{η}})]}{(1 - τ) + \frac{β}{1 + R} (R - \frac{τ}{1 - \bar{η}})} .

(17)

In Equation (17), the second law of thermodynamics was applied to calculate the total entropy production given by

Σ = \frac{{\bar{J}}_{2}}{T_{2}} - \frac{{\bar{J}}_{1}}{T_{1}}

(see Figure 1). Figure 2 shows the behavior of the three objective functions versus internal efficiency. As can be seen from Figure 2 there exists an optimal efficiency value which depends on the parameter R, the parameter τ and the economic parameter β. In addition, Figure 2 shows for the three performance regimes how the optimal efficiency tends to the Carnot value when

β ⟶ \infty

. Analogously to Equation (12), calculating the derivatives of

a {\bar{F}}_{E P}

and

a {\bar{F}}_{E}

with respect to

\bar{η}

, and solving for the efficiency the following two equations

{\frac{d (a {\bar{F}}_{E P})}{d \bar{η}}|}_{\bar{η} = {\bar{η}}^{*}} = 0

and

{\frac{d (a {\bar{F}}_{E})}{d \bar{η}}|}_{\bar{η} = {\bar{η}}^{*}} = 0

and using Equation (14) we get,

{\bar{η}}_{E P} (f, τ, R) = 1 - \frac{1 + f}{4 R} τ - \sqrt{\frac{(1 - f) τ}{2 R} + \frac{{(1 + f)}^{2} τ^{2}}{16 R^{2}}} .

(18)

{\bar{η}}_{E} (f, τ, R) = 1 - \frac{f}{2 R} τ - \sqrt{\frac{(1 - f) (1 + τ) τ}{2 R} + \frac{f^{2} τ^{2}}{4 R^{2}}} .

(19)

Figure 2. Profit functions for

β = 0

and

β = 5

versus η with

R = 1

and

τ = 0.5

Figure 2. Profit functions for

β = 0

and

β = 5

versus η with

R = 1

and

τ = 0.5

Equations (18) and (19), represent the steady-state efficiencies working both under maximum-efficient power (

{\bar{η}}_{E P}

) and maximum ecological function conditions (

{\bar{η}}_{E}

), respectively. In analogous way to Equation (15), for the endoreversible case (

R = 1

) from Equations (18) and (19), when

f = 0

, (

β = 0

) we obtain

\bar{η} = 1 - τ / 4 - \sqrt{τ (8 + τ)} / 4

and

\bar{η} = 1 - \sqrt{τ (τ + 1) / 2}

, respectively, which were previously obtained by Yilmaz [14], and by Angulo-Brown [16,22,23], for the case of a Curzon-Ahlborn heat engine, working at maximum efficient power and maximum ecological function conditions, respectively.

We can see, in Figure 3, how the optimal efficiencies smoothly increase from the maximum efficiency point,

f = 0

(in each regime of operation), corresponding to energy sources where the investment is the preponderant cost up to the Carnot value for

f = 1

, that is, for energy sources where the fuel is the predominant cost [5], and we can also observe that the following inequality,

{\bar{η}}_{E} > {\bar{η}}_{E P} > {\bar{η}}_{M P}

is hold [10,11].

Figure 3. The steady-state efficiencies working under maximum power output (

{\bar{η}}_{M P}

), maximum-efficient power (

{\bar{η}}_{E P}

) and maximum ecological function (

{\bar{η}}_{E}

) conditions.

Figure 3. The steady-state efficiencies working under maximum power output (

{\bar{η}}_{M P}

), maximum-efficient power (

{\bar{η}}_{E P}

) and maximum ecological function (

{\bar{η}}_{E}

) conditions.

3. Local Stability Analysis

In this section, the local stability theory (see Appendix) is applied to the thermo-economical heat engine model mentioned previously. Following Santillán et al. [1], due to x and y are macroscopic objects (the working substance at the isothermal branches of the cycle) with heat capacity C [24,25], their temperatures change according to the following differential equations:

\frac{d x}{d t} = \frac{1}{C} [α (T_{1} - x) - J_{1}],

(20)

\frac{d y}{d t} = \frac{1}{C} [J_{2} - α (y - T_{2})],

(21)

where

J_{1}

and

J_{2}

are the heat flows from x to the working substance and from the Carnot engine to y, respectively. According to the non-endoreversibility hypothesis [19],

J_{1}

and

J_{2}

are given by

J_{1} = \frac{R x}{R x - y} P,

(22)

and

J_{2} = \frac{y}{R x - y} P .

(23)

On the other hand, Equations (6), (15), (18) and (19) are used to construct the expressions for τs which relate the internal variables x and y, to the external temperatures

T_{1}

and

T_{2}

, under maximum power, maximum efficient power and maximum ecological function conditions, respectively,

τ_{M P} = \frac{{\bar{y}}^{2}}{{\bar{x}}^{2} (1 - f) R + \bar{x y} f} .

(24)

In addition, for both efficient power and ecological function conditions performance, we obtain

τ_{P_{E P}} = \frac{2 {\bar{y}}^{2}}{{\bar{x}}^{2} (1 - f) R + \bar{x y} (1 + f)},

(25)

τ_{E} = \frac{(1 - f) R {\bar{x}}^{2} + 2 f \bar{x y} + \sqrt{8 (1 - f) R {\bar{x}}^{2} {\bar{y}}^{2} + {\bar{x}}^{2} {((f - 1) R \bar{x} - 2 f \bar{y})}^{2}}}{2 (f - 1) R {\bar{x}}^{2}}

(26)

The functional dependence of the internal and external temperatures in each regime of operation is necessary to analyze the local stability. Using the assumption [1] that out of the steady state but not too far away, the power output of a Curzon-Ahlborn heat engine depends on x and y in the same way that it depends on

\bar{x}

and

\bar{y}

at the steady-state

(P (\bar{x}, \bar{y}, f, R) ⟶ P (x, y, f, R))

; that is, this assumption is applicable only in the vicinity of the steady state, the dynamical equations for x and y are written as follows:

\frac{d x}{d t} = \frac{1}{C} [α (T_{1} - x) - \frac{R x}{R x - y} P (x, y, f, R)],

(27)

\frac{d y}{d t} = \frac{1}{C} [\frac{y}{R x - y} P (x, y, f, R) - α (y - T_{2})] .

(28)

To analyze the system stability near the steady state, we proceed following the steps described in the Appendix. For the case of maximum efficient power we get,

h (x, y, f, R) = \frac{1}{C} [α (T_{1} - x) - \frac{R x}{R x - y} \frac{(f - 1) {(y - R x)}^{2}}{(f - 1) R x - (2 R + f + 1) y}],

(29)

g (x, y, f, R) = \frac{1}{C} [\frac{y}{R x - y} \frac{(f - 1) {(y - R x)}^{2}}{(f - 1) R x - (2 R + f + 1) y} - α (y - T_{2})],

(30)

where

h (x, y, f, R)

and

g (x, y, f, R)

are defined in the Appendix. The case of Maximum Power conditions was reported in [12]. Similar expressions are obtained for the case of maximum ecological function, but they are quite lengthy. After solving the corresponding eigenvalue equation, we find that both eigenvalues (

λ_{1}

and

λ_{2}

) are function of α, C, τ, f and R. The final expression and the algebraic details are not shown but this can be easily reproduced with the help of a symbolic algebra package. Moreover, our calculations show that both eigenvalues are real and negative. Thus, the steady state is stable because any perturbation would decay exponentially. For the case

f = 0

, expressions for the eigenvalues previously obtained by Santillán et al. [1] and Guzmán-Vargas et al. [3] are recovered.

In Figure 4, Figure 5 and Figure 6, the relaxation times are plotted against τ for different values of fractional fuel cost f, for a fixed value of R (

R = 1

), that is, the endoreversible case. It is observed that

t_{1}

(Equation (A.8)) is a decreasing function of τ. This relaxation time decreases as the fuel cost increases, indicating a faster decay as

f ⟶ 1

. For

t_{2}

(see Equation (A.9)), this relaxation time remains almost constant for

f = 0

. As the fractional fuel cost f increases,

t_{2}

slowly increases too. In the limit

f ⟶ 1

, both relaxation times tend to be closer each other, but there is a stronger inequality

t_{2} < t_{1}

in the interval

0 < τ < 1

. These figures also show the relaxation times as a function of τ, for several values of the parameter R, and for a fixed value of the fractional fuel cost f. These figures show that

t_{1}

is a decreasing function of τ and decreases as the parameter R decreases and

t_{2}

remains almost constant when the irreversibility parameter changes. From the findings shown in Figure 4, Figure 5 and Figure 6, it is concluded that the system is stable for

τ > 0

. As the fractional fuel cost f increases,

t_{1}

decreases whereas

t_{2}

increases, for a given value of R. In contrast, for a given value of f, as the irreversibility parameter R decreases,

t_{1}

decreases, whereas

t_{2}

increases (see phase portrait [3] where the phase diagram are showed in details). The power output and the efficiency depend on τ for the cases analyzed here, and both energetic quantities are decreasing functions of this parameter, that is, the system’s stability moves in the opposite direction to that of the steady-state as P, η and τ varies.

Figure 4. Plot of relaxation times under maximum power conditions versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1).

Figure 5. Plot of relaxation times under maximum efficient power versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1).

Figure 6. Plot of relaxation times under maximum ecological function conditions versus τ for (a) several values of the endorreversibility parameter and a value of the fractional fuel cost and (b) for several values of the fractional fuel cost f in the endoreversible case (R = 1).

Additionally, in Figure 7, for the cases of maximum efficient power and maximum ecological function conditions, show the relaxation times versus fuel fractional cost for several values of τ. In these cases, it can be seen, how the fast (slow) relaxation time slightly increases (decreases) as f changes from 0 to 1.

Figure 7. Relaxation times in the endoreversible case (

R = 1

) versus fractional fuel cost for several values of τ for (a) Maximum efficient power conditions and (b) Maximum ecological function.

Figure 7. Relaxation times in the endoreversible case (

R = 1

) versus fractional fuel cost for several values of τ for (a) Maximum efficient power conditions and (b) Maximum ecological function.

Finally, Figure 8 shows the ratio of the relaxation times previously calculated for the three regimes of operation considered. In this Figure, the inequality

τ_{E P} > τ_{E} > τ_{M P}

holds.

Figure 8. Ratio of the relaxation times versus τ for a value of the fractional fuel cost and several values of the parameter R (cases (a), (c) and (e)), and for the endoreversible case

R = 1

, for different values of the fractional fuel cost, (cases (b), (d) and (f)).

Figure 8. Ratio of the relaxation times versus τ for a value of the fractional fuel cost and several values of the parameter R (cases (a), (c) and (e)), and for the endoreversible case

R = 1

, for different values of the fractional fuel cost, (cases (b), (d) and (f)).

4. Conclusions

This work presented a local stability analysis of a thermo-economic model of an irreversible heat engine working under different performance regimes: The Maximum Power Output, Maximum Efficient Power and the Ecological Function regime considering a linear heat transfer law (the Newtonian law case). The relaxation times are shown as function of α, C, τ, f and R; that is, they depend on the materials that separate the working fluid from the reservoirs (through α); and on the working fluid (through C); on the reservoir temperatures (through τ); on fractional fuel costs (through f, which is associated with several energy sources from renewable energy to natural gas, (see Table 1. reported by De Vos [5]) and on the internal irreversibilities (through R). After a small perturbation, the system exponentially declines to the steady state characterized by two different relaxation times. In addition, the relaxation time under maximum efficient power conditions is less than the relaxation times under both maximum ecological function and maximum efficient power; that is, under maximum efficient power conditions, there is better stability conditions than for the other two regimes, although these differences are not large. Our cycle model is an FTT simple version of a Carnot-type engine, but, additionally, it considers the fractional fuel costs. Finally, the results confirm the local stability properties of the CA-engine for the three performance regimes.

Acknowledgments

The authors acknowledge CONACYT, EDI-IPN, and COFAA-IPN for supporting this work. The authors also thank Fernando Angulo-Brown for the fruitful comments and suggestions relating to this work. Israel Reyes-Ramírez expresses their thanks to the hospitality of Universidad de Salamanca, Spain.

Author Contributions

All three authors contributed equally to the present work. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix

A. Linearization and Stability Analysis

This section presents a brief description of the linear stability analysis of a two-dimensional system [26]. Consider the dynamical system,

\frac{d x}{d t} = h (x, y),

(A.1)

\frac{d y}{d t} = g (x, y),

(A.2)

where h an g are functions of x and y. Let (

\bar{x}, \bar{y}

) be a fixed point such that

h (\bar{x}, \bar{y}) = 0

and

g (\bar{x}, \bar{y}) = 0

. Consider a small perturbation around this fixed point and write

x = \bar{x} + δ x

and

y = \bar{y} + δ y

, where

δ x

and

δ y

are small disturbances from the corresponding fixed point values. By substituting into Equations (A.1) and (A.2), expanding

h (\bar{x} + δ x, \bar{y} + δ y)

and

g (\bar{x} + δ x, \bar{y} + δ y)

in a Taylor series, and using the fact that

δ x

and

δ y

are too small to neglect quadratic terms, the following equations are obtained for the perturbations:

\begin{matrix} (\begin{matrix} \frac{d δ x}{d t} \\ \frac{d δ y}{d t} \end{matrix}) = (\begin{matrix} h_{x} & h_{y} \\ g_{x} & g_{y}, \end{matrix}) (\begin{matrix} δ x \\ δ y \end{matrix}) \end{matrix}

(A.3)

where

{h_{x} = \frac{\partial h}{\partial x}|}_{\bar{x}, \bar{y}}

{h_{y} = \frac{\partial h}{\partial y}|}_{\bar{x}, \bar{y}}

{g_{x} = \frac{\partial g}{\partial x}|}_{\bar{x}, \bar{y}}

and

{g_{y} = \frac{\partial h}{\partial y}|}_{\bar{x}, \bar{y}}

. Equation (A.3) is a linear system of differential equations. Thus, the general solution of the system is of the form,

δ \vec{r} = e^{λ t} \vec{u},

(A.4)

with

δ \vec{r} = (δ x, δ y)

and

\vec{u} = (u_{x}, u_{y})

. Substitution of the solution

δ \vec{r}

into Equation (A.3) yields the following eigenvalue equation:

A δ \vec{r} = λ δ \vec{u},

(A.5)

where A is the matrix given by the first term on the right-hand side of Equation (A.3). The eigenvalues of this equation are the roots of the characteristic equation,

|A - λ I| = (h_{x} - λ) (g_{y} - λ) - g_{x} h_{y} = 0 .

(A.6)

λ_{1}

and

λ_{2}

are solutions of Equation (A.6), the general solution of the system is

δ \vec{r} = c_{1} e^{λ_{1} t} {\vec{u}}_{1} + c_{2} e^{λ_{2} t} {\vec{u}}_{2},

(A.7)

where

c_{1}

and

c_{2}

are arbitrary constants and

{\vec{u}}_{1}

and

{\vec{u}}_{2}

are the eigenvectors corresponding to

λ_{1}

and

λ_{2}

, respectively. To determine

{\vec{u}}_{1}

and

{\vec{u}}_{2}

we use Equation (A.5) again for each eigenvalue. Information about the stability of the system can be obtained from the eigenvalues

λ_{1}

and

λ_{2}

. In general,

λ_{1}

and

λ_{2}

are complex numbers. If both

λ_{1}

and

λ_{2}

have negative real parts, the fixed point is stable. Moreover, if both eigenvalues are real and negative, the perturbations decrease exponentially. In this last case, it is possible to identify relaxation times for each eigendirection as,

t_{1} = \frac{1}{|λ_{1}|},

(A.8)

t_{2} = \frac{1}{|λ_{2}|} .

(A.9)

References

Santillán, M.; Maya, G.; Angulo-Brown, F. Local stability analysis of and endoreversible Curzon-Ahlborn-Novikov engine working in a maximum-power like regime. J. Phys. D Appl. Phys. 2001, 34, 2068–2072. [Google Scholar] [CrossRef]
Curzon, F.L.; Ahlborn, B. Efficiency of a Carnot engine at maximum power output. Am. J. Phys. 1975, 43, 22–24. [Google Scholar] [CrossRef]
Guzman-Vargas, L.; Reyes-Ramirez, I.; Sánchez, N. The effect of heat transfer laws and thermal conductances on the local stability of an endoreversible heat engine. J. Phys. D Appl. Phys. 2005, 38, 1282–1291. [Google Scholar] [CrossRef]
Paéz-Hernández, R.; Angulo-Brown, F.; Santillán, M. Dynamic Robustness and Thermodynamic Optimization in a Non-Endoreversible Curzon-Ahlborn Engine. J. Non-Equilib. Thermodyn. 2006, 31, 173–188. [Google Scholar] [CrossRef]
De Vos, A. Endoreversible thermoeconomics. Energy Convers. Manag. 1995, 36, 1–5. [Google Scholar] [CrossRef]
Novikov, I.I. The efficiency of atomic power stations (a review). J. Nucl. Energy (1954) 1958, 7, 125–128. [Google Scholar] [CrossRef]
Chambadal, P. Les Centrales Nucléaires; Armand Colin: Paris, France, 1967. (In French) [Google Scholar]
Sahin, B.; Kodal, A. Performance analysis of an endoreversible heat engine based on a new thermoeconomic optimization criterion. Energy Convers. Manag. 2001, 42, 1085–1093. [Google Scholar] [CrossRef]
Wu, C.; Chen, L.; Chen, J. Recent Advances in Finite-Time Thermodynamics; Nova Science Publishers: New York, NY, USA, 1999. [Google Scholar]
Barranco-Jiménez, M.A.; Angulo-Brown, F. Thermoeconomic optimisation of Novikov power plant model under maximum ecological conditions. J. Energy Inst. 2007, 80, 96–104. [Google Scholar] [CrossRef]
Barranco-Jiménez, M.A.; Angulo-Brown, F. Thermoeconomic optimisation of endoreversible heat engine under maximum modified ecological criterion. J. Energy Inst. 2007, 80, 232–238. [Google Scholar] [CrossRef]
Barranco-Jiménez, M.A.; Páez-Hernández, R.T.; Reyes-Ramírez, I.; Guzmán-Vargas, L. Local Stability Analysis of a Thermo-Economic Model of a Chambadal-Novikov-Curzon-Ahlborn Heat Engine. Entropy 2011, 13, 1584–1594. [Google Scholar] [CrossRef]
Barranco-Jiménez, M.A.; Cervantes-Espinoza, L.; Hurtado-Aguilar, D.; Reyes-Ramirez, I.; Guzmán-Vargas, L. Local stability Analysis of thermoeconomic model of a Curzon-Ahlborn heat engine with a Dulong-Petit heat transfer law. In Proceedings of the 24th International Conference on ECOS 2011, Novi Sad, Serbia, 4–7 July 2011; pp. 410–418.
Yilmaz, T. A new performance criterion for heat engines: Efficient power. J. Energy Inst. 2006, 79, 38–41. [Google Scholar] [CrossRef]
Arias-Hernández, L.A.; Barranco-Jiménez, M.A.; Angulo-Brown, F. Comparative analysis of two ecological type modes of performance for a simple energy converter. J. Energy Inst. 2009, 82, 223–227. [Google Scholar] [CrossRef]
Angulo-Brown, F. An ecological optimization criterion for finite-time heat engines. J. Appl. Phys. 1991, 69, 7465–7469. [Google Scholar] [CrossRef]
Arias-Hernández, L.A.; Angulo-Brown, F. A general property of endoreversible thermal engines. J. Appl. Phys. 1997, 81, 2973–2979. [Google Scholar] [CrossRef]
Chen, J. The maximum power output and maximum efficiency of an irreversible Carnot heat engine. J. Phys. D Appl. Phys. 1994, 27. [Google Scholar] [CrossRef]
Arias-Hernández, L.A.; AresdeParga, G.; Angulo-Brown, F. On Some Nonendoreversible Engine Models with Nonlinear Heat Transfer Laws. Open Syst. Inf. Dyn. 2003, 10, 351–375. [Google Scholar] [CrossRef]
Barranco-Jiménez, M.A. Finite-time thermodynamics optimization of a non endoreversible heat engine. Rev. Mex. Phys. 2009, 55, 211–220. [Google Scholar]
Stucki, J.W. The Optimal Efficiency and the Economic Degrees of Coupling of Oxidative Phosphorylation. Eur. J. Biochem. 1980, 109, 269–283. [Google Scholar] [CrossRef] [PubMed]
Arias-Hérnandez, L.A.; Angulo-Brown, F. Thermodynamic optimization of endoreversible engines. Rev. Mex. Fis. 1994, 40, 866–877. [Google Scholar]
Calvo-Hernández, A.; Medina, A.; Roco, J.M.M.; White, J.A.; Velasco, S. Unified optimization criterion for energy converters. Phys. Rev. E 2001, 63, 037102. [Google Scholar] [CrossRef] [PubMed]
Wang, Y.; Tu, Z.C. Efficiency at maximum power output of linear irreversible Carnot-like heat engines. Phys. Rev. E 2012, 85, 011127. [Google Scholar] [CrossRef] [PubMed]
Gonzalez-Ayala, J.; Arias-Hernandez, L.A.; Angulo-Brown, F. Connection between maximum-work and maximum-power thermal cycles. Phys. Rev. E 2013, 88, 052142. [Google Scholar] [CrossRef] [PubMed]
Strogatz, S.H. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering; Westview Press: Boulder, CO, USA, 2000. [Google Scholar]

© 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons by Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Barranco-Jiménez, M.A.; Sánchez-Salas, N.; Reyes-Ramírez, I. Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes. Entropy 2015, 17, 8019-8030. https://doi.org/10.3390/e17127860

AMA Style

Barranco-Jiménez MA, Sánchez-Salas N, Reyes-Ramírez I. Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes. Entropy. 2015; 17(12):8019-8030. https://doi.org/10.3390/e17127860

Chicago/Turabian Style

Barranco-Jiménez, Marco A., Norma Sánchez-Salas, and Israel Reyes-Ramírez. 2015. "Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes" Entropy 17, no. 12: 8019-8030. https://doi.org/10.3390/e17127860

APA Style

Barranco-Jiménez, M. A., Sánchez-Salas, N., & Reyes-Ramírez, I. (2015). Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes. Entropy, 17(12), 8019-8030. https://doi.org/10.3390/e17127860

Article Menu

Local Stability Analysis for a Thermo-Economic Irreversible Heat Engine Model under Different Performance Regimes

Abstract

1. Introduction

2. Thermoeconomic Optimization of a Curzon-Ahlborn Engine Model at Different Regimes of Performance

3. Local Stability Analysis

4. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

Appendix

A. Linearization and Stability Analysis

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI