Next Article in Journal
On Epistemics in Expected Free Energy for Linear Gaussian State Space Models
Next Article in Special Issue
Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle
Previous Article in Journal
History Marginalization Improves Forecasting in Variational Recurrent Neural Networks
Previous Article in Special Issue
Performance Analysis and Four-Objective Optimization of an Irreversible Rectangular Cycle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities

by
David Pérez-Gallego
1,
Julian Gonzalez-Ayala
1,2,*,
Antonio Calvo Hernández
1,2 and
Alejandro Medina
1,2
1
Departamento de Física Aplicada, Universidad de Salamanca, 37008 Salamanca, Spain
2
Instituto Universitario de Física Fundamental y Matemáticas (IUFFyM), Universidad de Salamanca, 37008 Salamanca, Spain
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(12), 1564; https://doi.org/10.3390/e23121564
Submission received: 22 October 2021 / Revised: 17 November 2021 / Accepted: 19 November 2021 / Published: 24 November 2021

Abstract

:
A model for a pumped thermal energy storage system is presented. It is based on a Brayton cycle working successively as a heat pump and a heat engine. All the main irreversibility sources expected in real plants are considered: external losses arising from the heat transfer between the working fluid and the thermal reservoirs, internal losses coming from pressure decays, and losses in the turbomachinery. Temperatures considered for the numerical analysis are adequate for solid thermal reservoirs, such as a packed bed. Special emphasis is paid to the combination of parameters and variables that lead to physically acceptable configurations. Maximum values of efficiencies, including round-trip efficiency, are obtained and analyzed, and optimal design intervals are provided. Round-trip efficiencies of around 0.4, or even larger, are predicted. The analysis indicates that the physical region, where the coupled system can operate, strongly depends on the irreversibility parameters. In this way, maximum values of power output, efficiency, round-trip efficiency, and pumped heat might lay outside the physical region. In that case, the upper values are considered. The sensitivity analysis of these maxima shows that changes in the expander/turbine and the efficiencies of the compressors affect the most with respect to a selected design point. In the case of the expander, these drops are mostly due to a decrease in the area of the physical operation region.

1. Introduction

Increasing production of electric energy from renewable sources like wind or Sun makes essential energy storage technologies due to the inherent intermittency of such natural resources [1]. Substantial efforts are being devoted to the so-called Pumped Thermal Energy Storage (PTES) systems [2]. They are aimed at storing energy during the hours with an excess of, for instance, wind or photovoltaic production. When electric energy is subsequently required, heat is transformed again in electricity through some thermodynamic cycle. Required storage periods could be not very long, from four to eight hours. At grid-scale, the only proven competitors of this concept are Pumped Hydro Energy Storage and Compressed Air Energy Storage [3]. Nevertheless, these technologies require very particular geographical and geological scenarios that make difficult their implementation where required [4]. On the other hand, electric batteries, as with lithium ones, are still very expensive at grid-scale.
Recently, Dumont et al. [5] have proposed the term Carnot batteries to those systems that are used to store energy by establishing a temperature difference between two media (low and high temperature reservoirs). In the charge period, energy is stored (usually by means of a thermal heat pump). And after a storage period, electric energy is produced by a heat engine that works between those reservoirs. Several configurations and thermodynamic cycles can be considered [6]. To date, only a few prototype installations were built [7], and it is still necessary a research work at a fundamental level to provide optimal plant configurations and adequate intervals for the most significant design parameters. This work is focused on Brayton-like cycles [8] and constant temperature heat reservoirs as packed bed solid reservoirs at equilibrium.
The working fluids in Brayton-like thermodynamic cycles are gases stable at high temperatures, chemically inert, cheap, and environmentally friendly [9]. These requirements limit the options to Ar and N 2 (or air). Ar is actually the preferred option, because it can reach higher temperatures for the same storage pressure ratios [5], although there is more experimental background in the turbomachinery working with air. In any case, the ideal gas equation of state is sufficient for making efficiency and cost estimations. Moreover, in the case of Ar, specific heats are temperature independent. Packed bed thermal storage systems have been widely analyzed in the literature. A recent review on types of packed bed systems, the most relevant experimental works in the last years, and their numerical modeling are due to Esence et al. [10].
The use of packed bed sensible heat storage makes it possible to reach high temperatures, and there is a direct heat transfer between the working fluid in the thermodynamic cycle and the reservoir itself [5]. This avoids the use of intermediate heat exchangers, and overall efficiency is supposedly higher. Moreover, only one deposit at a high temperature is required at difference with the case of liquid storage. This kind of thermal storage is especially interesting for not overly long storage periods [5], for instance, of the order of several hours, and for Brayton cycles, whose efficiency when working as a heat engine is large for high turbine inlet temperatures.
When modeling this type of PTES configuration, several key points should be taken into account: the losses in the heat transfer between the gas working fluid and the solid particulates, the pressure decay in the heat exchanges, the efficiencies of expanders and compressors, and the heat leakage through the containers walls during storage. It has been demonstrated that in optimized designs, the round-trip efficiency is more sensitive to the performance of expanders, turbines, and compressors than to the losses coming from the packed bed thermal storage towards the surroundings [5]. As stated by McTigue et al. [11] heat leakage can be reduced to any desired level with an appropriate investment in storage insulation. Numerically a 1% heat leakage from each reservoir could be reflected in a reduction of about 2% in round-trip efficiency.
The objective of the present study is to propose a thorough thermodynamic model for the overall round-trip performance of a PTES system composed of a Brayton heat pump and a Brayton heat engine during charge and discharge periods, respectively, and considering solid storage reservoirs at an approximately constant temperature. All internal and external irreversibilities in heat exchanges and during compression and expansion processes will be considered. The system will be assumed to be at stationary conditions with fixed temperature reservoirs to obtain analytical equations allowing for detailed sensitivity analyses of the influence of all possible losses and different temperature levels in system performance. In previous studies by our group [12,13], PTES Brayton-like systems with liquid storage and their optimization were developed. This work is intended to present a complete and flexible thermodynamic model capable to include all the most relevant peculiarities and thermodynamic details emerging from the combination of Brayton-like heat devices and solid thermal reservoirs. Special emphasis is placed on the links between charge and discharge cycles, pressure drops during heat exchanges, and the losses in turbines, expanders, and compressors. The relative importance of external and internal irreversibilities is analyzed, as well as the influence of each of them in the optimal values of objective functions, such as round-trip efficiency or power output.

2. Thermodynamic Model

The PTES storage system is made up of two irreversible Brayton engines, one working as a heat pump and another as a heat engine. The first operates in a charge stage, from which heat is transferred to the hot reservoir and energy is stored. This is used later as an input heat for the discharge stage performed by the heat engine. The working fluid interacts with two reservoirs (hot or cold storage) at a constant temperature. In this process, energy is used or produced. The combination in a sequence of both cycles with a time delay between them constitutes the energy storage system. That time delay is equivalent to the intended storage period.
In the model developed it will be considered that heat reservoirs operate at a constant temperature, i.e., after a non-stationary transient period they reach an equilibrium temperature. This simplified stationary model will allow for detailed sensitivity analysis of the key thermodynamic parameters of the system and to optimize the elected objective functions. Thus, thermodynamic optimum pre-design intervals for the most influential parameters will be proposed. Another key hypothesis is that working fluid specific heats are approximately constant in the temperature interval considered. This is true for Ar or other noble gases and is not a bad supposition for diatomic gases such as air. In this way, the heat capacity of the working fluid per unit time, C w , that depends on the mass flow, m ˙ , and the specific heat at constant pressure, c p , w , can be written as: C w = m ˙ c p , w .

2.1. Thermodynamic Analysis: Heat Pump Cycle

A complete (round-trip) cycle begins with the charge cycle. It is considered as an inverse Brayton thermodynamic cycle. The working fluid absorbs heat from a cold storage at temperature T L , and using a network from the external energy source (wind, PV, or whichever) through a compressor, it is compressed by increasing its pressure and temperature. In this way, it can transfer heat to a hot storage at T H . In the process, energy is then released in an expander, and gas temperature and pressure decrease in an approximately isentropic path. A components scheme of this cycle is shown in Figure 1.
The charge cycle, in a T S diagram, is shown in the inset of Figure 1. It runs in an anti-clockwise direction. The compressor increases working fluid pressure (and temperature) from state 3 to state 2 ( 3 2 ). Ideally, it is a reversible isentropic compression, but the model considers an entropy increase through the consideration of isentropic compressor efficiency, ϵ c H P 1 , defined as:
ϵ c H P = W ˙ c , i d e a l H P W ˙ c H P = T 2 S T 3 T 2 T 3 ,
where T 2 S represents the temperature the working fluid eventually would reach after an ideal compression. Thus, ϵ c H P , is the ratio between the work required in ideal conditions and the one consumed in real conditions.
At state 2, the working fluid reaches the maximum pressure and temperature in the cycle. Then, the working fluid releases heat to the hot storage (process 2 1 ). This process is isobaric in the ideal case. To account for the real pressure losses, Δ p H , a pressure decay parameter is defined [12]:
ρ H H P = p 1 p 2 κ = p 2 Δ p H p 2 κ ,
where κ ( γ 1 ) / γ is defined by convenience and γ is the adiabatic coefficient of the working fluid that is considered approximately constant in the working temperature intervals.
At state 1, the working fluid has reduced its temperature, due to the transfer of heat to the hot storage. Subsequently, the working fluid enters the expander. Its temperature and pressure decrease up to state 4 ( 1 4 ). Like in compression, the isentropic expander efficiency is considered to account for losses:
ϵ t H P = W ˙ t H P W ˙ t , i d e a l H P = T 1 T 4 T 1 T 4 S .
Finally, and to close the thermodynamic cycle, the fluid absorbs heat from a cold storage. This storage is at a lower temperature than the hot storage, but at a higher temperature than the working fluid in state 4. The temperature of the cold reservoir is denoted as T L . The heat absorption in step 4 3 ideally is an isobaric process, but pressure drops in the real process are quantified through:
ρ L H P = p 3 p 4 κ = p 4 Δ p L p 4 κ ,
For the correct operation of the charge cycle, the temperatures of each state must comply in a sequential order. The temperature of each state is limited as follows.
0 T 4 H P T 3 H P T L T H T 1 H P T 2 H P .
The work per unit time consumed by the compressor and the work per unit of time released at the expander are:
W ˙ c H P = C w ( T 2 T 3 ) ; W ˙ t H P = C w ( T 1 T 4 ) ,
so, the power output required to run the pump can be written as:
P i n = W ˙ i n = W ˙ c H P W ˙ t H P = C w ( T 2 T 3 T 1 + T 4 ) ,
where W ˙ c H P > W ˙ t H P , in absolute value. The heat per unit time that the working fluid transfers to the hot storage in step 2 → 1, is Q ˙ H H P = C w ( T 2 T 1 ) . Similarly, the working fluid absorbs heat from a cold storage in step 4 → 3, Q ˙ L H P = C w ( T 3 T 4 ) . The pressure ratios of the compressor and expander can be expressed as: r c H P = p 2 / p 3 and r t H P = p 1 / p 4 . These parameters give us the relationship between the pressure of the working fluid in the inlet and the outlet of the compressor and expander. In terms of the temperatures it is easy to show that:
a c H P = T 2 S T 3 = r c κ = p 2 p 3 κ = p 2 p 4 Δ p L κ ,
a t H P = T 1 T 4 S = r t κ = p 1 p 4 κ = p 2 Δ p H p 4 κ ,
and so, there is a connection among these parameters:
a t H P = a c H P ρ H H P ρ L H P .
The new parameters were introduced to allow us to rewrite the work consumed by the compressor and the work done by the turbine. Using Equations (1), (3), (8) and (9):
W ˙ c H P = C w ( T 2 T 3 ) = C w T 3 a c 1 ϵ c ,
W ˙ t H P = C w ( T 1 T 4 ) = C w T 1 ϵ t a t 1 a t .
It is considered that there are losses in the heat transfers between the working fluid and the reservoirs. To quantify these external irreversibilities, the maximum heats which the working fluid would exchange under ideal conditions, reaching the temperatures of the hot storage ( T H ) and the cold storage ( T L ) would be:
Q ˙ H , m a x H P = C w ( T 2 T H ) ,
Q ˙ L , m a x H P = C w ( T L T 4 ) .
The parameters which characterize these external irreversibilities can be defined as the ratios between the actual heat transferred and the maximum heat transferred under ideal conditions.
ϵ H H P = Q ˙ H H P Q ˙ H , m a x H P = T 2 T 1 T 2 T H ,
ϵ L H P = Q ˙ L H P Q ˙ L , m a x H P = T 3 T 4 T L T 4 .
In a heat pump, the coefficient of performance (COP), ν , quantifies the quotient between the heat which the fluid releases to the hot storage and the work input used for that task:
ν Q ˙ H H P W ˙ i n = Q ˙ H H P Q ˙ H H P Q ˙ L H P = T 2 T 1 T 2 T 1 T 3 + T 4 .
In the last equality, it was used that the internal energy variation in a cyclic process is zero. All of the temperatures and, therefore, the COP can be written analytically in terms of geometric parameters (pressure and temperature ratios, and the adiabatic coefficient of the working fluid) and the parameters defined to account for both internal ( ρ H H P , ρ L H P , ϵ c H P , ϵ t H P ) and external irreversibilities ( ϵ H H P , ϵ L H P ). From Equations (1), (3), (10), (15) and (16) the four temperatures and a t H P can be expressed in terms of the irreversible parameters and a c H P , resulting in:
T 1 = a c H P ρ H H P ρ L H P T L ( ϵ H H P 1 ) ϵ L H P ( a c H P + ϵ c H P 1 ) T H ϵ c H P ϵ H H P D H P ,
where D H P is use to simplify the expression, being
D H P = ϵ t H P ( ϵ H H P 1 ) ( ϵ L H P 1 ) ( a c H P + ϵ c H P 1 ) ( a c H P ρ H H P ρ L H P 1 ) a c H P ρ H H P ρ L H P ϵ L H P ( ϵ H H P 1 ) ( a c H P + ϵ c H P 1 ) a c H P ϵ H H P + a c H P ϵ c H P ϵ H H P + ϵ H H P 1 ,
in the same way, the other temperatures are
T 2 = ( a c H P + ϵ c H P 1 ) T H ϵ H H P ( ϵ L H P 1 ) a c H P ρ H H P ρ L H P ( ϵ t H P 1 ) ϵ t H P + a c H P ρ H H P ρ L H P T L ϵ L H P D H P ,
T 3 = ϵ c H P T H ϵ H H P ( ϵ L H P 1 ) a c H P ρ H H P ρ L H P ( ϵ t H P 1 ) ϵ t H P + a c H P ρ H H P ρ L H P T L ϵ L H P D H P ,
T 4 = a c H P ρ H H P ρ L H P ( ϵ t H P 1 ) ϵ t H P T H ϵ c H P ϵ H H P T L ( ϵ H H P 1 ) ϵ L H P ( a c H P + ϵ c H P 1 ) a c H P ρ H H P ρ L H P ϵ c H P + ( ϵ H H P 1 ) ( ϵ L H P 1 ) ( a c H P + ϵ c H P 1 ) a c H P ρ H H P ρ L H P ( ϵ t H P 1 ) ϵ t H P ,

2.2. Thermodynamic Analysis: Heat Engine Cycle

The discharge cycle after the storage period is also considered as a Brayton-like thermodynamic cycle, in this case corresponding to a heat engine. This cycle is analogous to the charge one, but the working fluid runs in the opposite direction in the T S diagram, in a clockwise direction (see Figure 2).
In the discharge cycle, the working fluid enters the compressor at state 4. Working fluid is compressed up to state 1 ( 4 1 ). Losses in the compressor are accounted for through its isentropic efficiency, ϵ c H E . Afterward, in the heat input process ( 1 2 ) the working fluid is heated up to the turbine inlet temperature. Pressure decay in the heat input is considered through the parameter ρ H H E . Irreversibilities in the heat transfer from the reservoir to the working fluid are quantified with the parameter ϵ H H E . Subsequently, the working fluid enters the turbine at its maximum temperature, T 2 . The expansion to state 3 ( 2 3 ) is also considered as non-isentropic. This internal irreversibility is considered through ϵ t H E . Finally, a heat release to the cold storage allows for a cyclic process for the working fluid. The corresponding external irreversibilities are represented by ϵ H H E . For the correct operation of the discharge cycle, the temperatures of each state comply with an order analogous to the charge cycle. The temperature of each state is limited as follows:
T L T 4 H E T 3 H E T 2 H E T H ,
T L T 4 H E T 1 H E T 2 H E T H .
Internal and external irreversibilities parameters are defined next:
ϵ c H E = W ˙ c , i d e a l H E W ˙ c H E = T 1 S T 4 T 1 T 4 ; ϵ t H E = W ˙ t H E W ˙ t , i d e a l H E = T 2 T 3 T 2 T 3 S
ρ H H E = p 2 p 1 κ = p 1 Δ p H p 1 κ ; ρ L H E = p 4 p 3 κ = p 3 Δ p L p 3 κ
ϵ H H E = Q ˙ H H E Q ˙ H , m a x H E = T 2 T 1 T H T 1 ; ϵ L H E = Q ˙ L H E Q ˙ L , m a x H E = T 3 T 4 T 3 T L .
The work per unit time required by the compressor and the one produced by the turbine are:
W ˙ c H E = C w ( T 1 T 4 ) ; W ˙ t H E = C w ( T 2 T 3 ) ,
so, the power output is:
P o u t = W ˙ o u t = W ˙ c H E W ˙ t H E = C w ( T 1 T 4 T 2 + T 3 ) ,
The heat per unit time the working fluid absorbs from the hot storage in step 1 → 2, is defined as follows:
Q ˙ H H E = C w ( T 2 T 1 ) .
The temperature of the working fluid at state 2 must be lower than the temperature of the hot storage because the heat is transferred from the storage to the working fluid. In the same way, the working fluid releases heat to a cold storage in step 3 → 4. The temperature of the working fluid in state 4 must be higher than the temperature of the cold reservoir.
Q ˙ L H E = C w ( T 3 T 4 ) .
The pressure ratios of the compressor and the turbine are:
r c H E = p 1 p 4 ; r t H E = p 2 p 3 ,
and the temperature ratios,
a c H E = T 1 S T 4 = r c κ ; a t H E = T 2 T 3 S = r t κ ,
so,
a t H E = a c H E ρ H H E ρ L H E .
The new parameters introduced allow us to rewrite the work consumed by the compressor and the work done by the turbine as:
W ˙ c H E = C w T 4 a c 1 ϵ c ; W ˙ t H E = C w T 2 ϵ t a t 1 a t .
Heat engine efficiency ( η ) quantifies the amount of useful energy in relation to the amount of heat input from the hot storage
η W ˙ o u t Q ˙ H H E = P o u t H E Q ˙ H H E = T 1 T 4 T 2 + T 3 T 2 T 1 .
All the cycle temperatures and therefore its efficiency can be expressed in terms of pressure and temperature ratios, and the parameters defined to account for internal (pressure losses and non-ideality of compressor and turbine) and external irreversibilities (heat transfer from the reservoirs to the working fluid). The following analytical expressions for cycle temperature are found from Equations (25), (27) and (34), in this case:
T 1 = ( a c H E + ϵ c H E 1 ) T H ϵ H H E ( ϵ L H E 1 ) a c H E ρ H H E ρ L H E ( ϵ t H E 1 ) ϵ t H E + a c H E ρ H H E ρ L H E T L ϵ L H E D H E
where, once again, the term D H E is introduced to simplify the expression, being
D H E = ϵ t H E ( ϵ H H E 1 ) ( ϵ L H E 1 ) ( a c H E + ϵ c H E 1 ) ( a c H E ρ H H E ρ L H E 1 ) a c H E ρ H H E ρ L H E ϵ L H E ( ϵ H H E 1 ) ( a c H E + ϵ c H E 1 ) a c H E ϵ H H E + a c H E ϵ c H E ϵ H H E + ϵ H H E 1 ,
in the same way, the other temperatures are
T 2 = a c H E ρ H H E ρ L H E T H ϵ c H E ϵ H H E T L ( ϵ H H E 1 ) ϵ L H E ( a c H E + ϵ c H E 1 ) D H E ,
T 3 = a c H E ρ H H E ρ L H E ( ϵ t H E 1 ) ϵ t H E T H ϵ c H E ϵ H H E T L ( ϵ H H E 1 ) ϵ L H E ( a c H E + ϵ c H E 1 ) D H E ,
T 4 = ϵ c H E T H ϵ H H E ( ϵ L H E 1 ) a c H E ρ H H E ρ L H E ( ϵ t H E 1 ) ϵ t H E a c H E ρ H H E ρ L H E T L ϵ L H E a c H E ρ H H E ρ L H E ϵ c H E + ( ϵ H H E 1 ) ( ϵ L H E 1 ) ( a c H E + ϵ c H E 1 ) a c H E ρ H H E ρ L H E ( ϵ t H E 1 ) ϵ t H E ,

2.3. Round-Trip Efficiency

Once the charge and discharge phases of the system have been modeled separately, a global or round-trip cycle can be defined and analyzed. The global cycle starts with the charging of heat towards the hot storage. The energy is stored for the required period of time, and later on it is converted into electric energy again through the discharge heat engine cycle. Each thermodynamic charge or discharge phase can be carried out one or multiple times depending on the energy demand and the dimensions of the PTES system. In this work, the PTES system will be analyzed over a single charge and discharge sequence. Figure 3 contains two representative examples of the coupling of charge and discharge cycles.
Overall performance of the energy storage cycle is usually evaluated through the round-trip efficiency, Φ [14]. It is defined as the ratio between the net energy obtained from storage during discharge and the net energy introduced into the system during the charge.
Φ = W o u t W i n = P o u t P i n ,
where W o u t refers to the useful energy obtained from storage, W i n is the work input in the charge stage. P o u t and P i n are the electrical net powers from the discharge and charge phases. For the last equality identical charge and discharge times are assumed in this model. Recently, the analysis of the case in which the charge and discharge times are different has been studied in an endoreversible machine model, studying the time influence on the performances for different operating regimes [15].
The performance of heat devices is obviously limited by the maximum theoretical efficiency of Carnot. This efficiency depends on the temperature of the storages between which the cycle operates. The individual charge and discharge efficiencies are linked to this theoretical restriction. However, the round-trip efficiency quantifies the energy obtained from the stored energy. In this process, the heat used in the discharge phase comes from the energy stored in the charge phase, Q H H P . The round-trip efficiency of the overall cycle, Φ is not limited by the theoretical Carnot efficiency. Using the definitions of COP in the charging mode and efficiency in the discharge phase the round-trip efficiency can be written as,
Φ = η ν Q ˙ H H E Q ˙ H H P .
If the energy stored in the hot storage is the same that the heat used for the discharge phase (without heat leak storage losses) then Q ˙ H H P = Q ˙ H H E and in consequence the round-trip efficiency is:
Φ = η ν .
The condition Q ˙ H H P = Q ˙ H H E constitutes an idealization in this analysis, providing upper values for the coupled system performance. This constraint establishes a relationship between the compression ratios of the working fluid at the inlet and outlet of the compressor in the charge/discharge phases, from which a c H P = a c H P ( a c H E ) .

2.4. Endoreversible Limit

In the framework of the model considered, it is possible to recover the endoreversible limit as a particular case where the only losses come from the external irreversibilities due to the coupling of the working fluid with the heat storage, i.e., in this limit compressors, turbines, and expanders isentropic efficiencies are considered as 1, as well as the parameters ρ H and ρ L for charge and discharge (no pressure decays in heat inputs or releases). Thus, the overall cycle is internally reversible, and so the compression and expansion steps occur without heat exchanges between the working fluid and the surroundings. Thus, processes 3 2 and 1 4 of the charge phase or 2 3 and 4 1 in discharge phase correspond to a constant entropy vertical line in the T S diagram. In the same way, as there are no internal irreversibilities, the pressure drops are zero in processes 1 2 and 3 4 in charge or in the steps 2 1 and 4 3 in discharge. They are isobaric processes for the working fluid that is assumed to behave as an ideal gas. Figure 4 depicts the typical shape of a complete cycle in the endoreversible limit. In this example, the temperatures of the hot and cold storages are set at realistic values for solid thermal storage: T H = 1000 K and T L = 300 K.
Analytical expressions for P o u t and Φ were calculated under endoreversible conditions with the aim of comparing our results with previous ones in the literature as a validation. The following expressions are found:
Φ = a c H P ( a c H E 1 ) ( T H a c H E T L ) ( a c H P 1 ) a c H E ( a c H P T L T H ) ,
P ¯ o u t = ( a c H E 1 ) ϵ H ϵ L ( a c H E T L T H ) a c H E ( ϵ H 1 ) ϵ L a c H E ϵ H .
Particularly, the round-trip efficiency in conditions of maximum net output power (power at the discharge), Φ m a x P o u t , was obtained in two steps. First, the expression of the net output power in the discharge with respect to the parameter a c H E is maximized. And second, the value obtained for a c H E is substituted in Φ . The results are shown in the next equations:
a c , m a x P o u t H E = T H T L ,
m a x P ¯ o u t = ϵ H ϵ L ( T H + T L 2 T H T L ) ϵ H + ϵ L ϵ H ϵ L ,
Φ m a x P o u t = 2 T H T L 2 T H + T L .
These equations are the same as those provided in Refs. [15,16,17]. Although in those works the values were obtained in the context of an irreversible Carnot model, it should be noted that these types of models are equivalent to Brayton-type cycles for both engines and refrigerators, as has been demonstrated in [18,19,20] for linear transfer laws and constant heat capacities. Figure 5 depicts the evolution of Φ m a x P o u t with the temperature of the hot reservoir. A monotonic increasing behavior is observed with numerical values for the round-trip efficiency relatively large, for instance about 0.57 for T H = 1000 K.
It is also possible to show that at endoreversible conditions, the value of a c leading to maximum round-trip efficiency is identical to that giving maximum power output, i.e., Equation (47), so maximum round-trip efficiency and power output at maximum round-trip efficiency are the same as those obtained at maximum power output, Equations (48) and (49) respectively. In summary, the proposed thermodynamic model, at the endoreversible limit, recovers previous results in the literature ensuring its reliability at least under these conditions. Nevertheless, the model is capable to dig in specifically in all the main irreversibility sources. This is the objective of the next section.

3. Results

3.1. Physically Acceptable Regions

This subsection presents an analysis of the links and constraints among the different variables of the overall system and overall charge/discharge cycle when all the losses that the model considers are taken into account together. It is a key point in a complex thermodynamic system like the one considered as to which are the physically acceptable regions and which are the constraints for the main parameters in order to manage an adequate design.
In each subsystem (heat pump or engine) there are six irreversibility parameters (pressure losses, irreversibilities in the heat transfer from the working fluid to the reservoirs, and losses in expanders/turbines or compressors) and one operation variable, a c H E . Moreover, due to the different options for possible storage media, it will be of interest to study the influence of the maximum temperature at which the storage device can operate. This adds up T H as a variable parameter. As it has been addressed in the context of energy storage in liquid media [12] there are physical constraints that define operation regions, that is, the possible values that the operation variables can acquire are restricted and the physical region lastly depends on the values of the irreversibility parameters. The combination of these values might prevent the system from reaching certain operation regimes. So far, the joint variation of such parameters has been analyzed in a multiobjective optimization scheme, obtaining a Pareto front for a given irreversibilities configuration [13]. Here, the analysis will be focused on the influence of each parameter in a certain operation regime, where the operation regime refers to the maximum available value of an objective function and not the global maximum, which might not be available.
Table 1 summarizes the constraints among the temperatures of the working fluid and pressure ratios that must be satisfied to obtain physically acceptable values of the heat transfers and thermodynamic figures of merit (COP values, heat engine efficiency, and round-trip efficiency).
Figure 6 shows the acceptable region considering the specific constraints and for a representative combination of parameters. The intersection is the physically adequate region to perform optimization and sensitivity analyses. The set of parameters characterizing irreversibilities range from typically (and optimistic) values in the literature and are considered for the so-called reference design point in the subsequent analyses [21,22,23]. They are given by: ϵ c H P , H E = ϵ t H P , H E = 0.95 , ϵ H H P , H E = ϵ L H P , H E = 0.95 , and ρ H H P , H E = ρ L H P , H E = 0.98 . Further on, a sensitivity analysis showing the importance of these parameters will be shown. In the figure, it was taken T L = 300 K and the working fluid was considered as monoatomic (for instance Ar) in which refers to its adiabatic coefficient. In other figures, a c H E will be taken as the independent variable. This avoids the election of a particular numerical value for the adiabatic coefficient. The figure shows the allowable correlations between T H and r c H E for the specified set of realistic values of irreversibility parameters. Numerical intervals for both T H and r c H E are realistic and acceptable. Other combinations of irreversibility parameters will lead to slightly different acceptable regions. This will be commented on in detail afterward. Once the possible physical regions are delimited, it is possible to analyze meaningful functions such as round-trip efficiency, net power output, or efficiency.
Figure 7 shows the possible attainable values that can be obtained from the physical region for the power output, HE efficiency, and round-trip efficiency, and also for the compressor temperature ratio of the heat pump. The latter linearly increases with T H as shown in Figure 7a. The net power output obtained in the discharge, P ¯ o u t , (see Figure 7b) increases with T H and has a parabolic behavior with a c H E , i.e., there is an optimum value of a c H E maximizing P ¯ o u t for each value of the temperature of the hot reservoir. The round-trip efficiency ( Φ ) and the discharge engine efficiency ( η ) have non-trivial behaviors with the temperature T H and with a c H E (see Figure 7c,d). In general, increasing the temperature T H (with constant a c H E ) implies increasing these functions, but after a particular temperature the physically acceptable region is narrower and both efficiencies start to decrease. On the other hand, if T H is set constant and a c H E changes, the functions will present a maximum value for a specific a c H E . So, optimizing the parameter a c H E and raising the temperature T H seems to be the most beneficial option for not particularly high temperatures as can also be seen in Figure 6, where an increase of T H reduces the acceptable interval for a c H E , which is an important consideration to account for. The optimization of these objective functions will be analyzed in the next section.
The thermodynamic functions of interest generally exhibit non-trivial behaviors in the acceptable physical region defined by the considered independent variables and, thus, at least two detailed analyses are in order: an optimization survey for the chosen independent variables for a particular set of irreversibility parameters and a sensitivity analysis of the optimum values of the objective functions with respect to irreversibilities. These studies will be shown in the next subsections.

3.2. Optimization of Overall Performance

The influence of irreversibilities in the coupled cycle will be performed in the so-called design point (see the caption of Figure 6), a representative configuration of irreversibilities that serves as a reference to compare the influence of each parameter. The value for the temperature T H is chosen such that a c H E has a large action region, allowing it to perform at several operation regimes and, of course, it should be a realistic temperature for solid heat reservoirs. A value T H = 1000 K fulfills these requirements.
Figure 8a represents the input/output powers as a function of a c H E . Notice that P i n is a monotonic decreasing function whereas P o u t is parabolic and displays a maximum slightly below a c H E = 1.8 . For a monoatomic working fluid as Ar, it corresponds to a pressure ratio r c H E 4.3 that is a realistic value for actual compressors. For air, that value of r c H E would be about 7.8, which is also an acceptable value. Figure 8b depicts the round-trip efficiency, the discharge cycle efficiency ( η ), and the charge cycle coefficient of performance (COP) (minus unity, to scale the graphic and compare the functions). The black point in all curves in Figure 8 corresponds to the maximum round-trip efficiency and the red point corresponds to maximum P o u t ( m a x P o u t ). As it can be seen, unlike in the endoreversible limit, both maxima are reached at different a c H E values, which nonetheless is in agreement with previously analyzed PTES systems [12,13]. It is interesting to note that (see Figure 8b) the maximum round trip leads to a high value of η , not exactly the corresponding maximum, but a numerical value close to it. Numerically, the maximum round-trip is obtained at a pressure ratio r c H E 6.4 for monoatomic gases and r c H E 13.4 for diatomic ones, values that are above those leading to maximum power output.
Figure 9 shows the parametric plots of round-trip versus several other thermodynamic functions ranging on all possible a c H E values. Maximum P ¯ o u t and max Φ are indicated by black and red dots, the region in between is usually considered as the optimal design region of the device. As it will be seen later, the physical constraints might forbid the system from reaching such operational states.

3.3. Irreversibility Analysis

For Brayton-like HP and HE cycles, optimistic values for the irreversibilities can be considered in the ranges depicted in Table 2. The sensitivity of the main operation regimes to different parameters and for charge/discharge stages has not yet been carefully addressed and, thus, it is the objective of this subsection.
The effects on the power output P o u t and in the round-trip ( Φ ) of the coupled cycle stemming from these parameters is presented in Figure 10. In Figure 10a–c the influence of the pressure drops is shown, as are the expander and compressor efficiencies, and the effectiveness of heat transfers, respectively, over the charge phase. Similarly, in Figure 10d–f it is shown the influence of the pressure drops; turbine and compressor efficiencies, and the effectiveness, respectively, over the discharge phase. In each figure, solid curves represent the power and round-trip values in the physically acceptable region. Each curve starts at P o u t = 0 with the minimum value of a c H E . As a c H E increases, the curve goes through the parametric loop in a clockwise direction. At a certain point, a maximum value of a c H E is reached (in the sense of Figure 6), which leads to the end of the physical region of our cycle (dashed lines begins) and the mentioned physical and practical constraints are not fulfilled. In other words, for each combination of irreversibility parameters, an intersection area as that shown in Figure 6 is built. Dashed lines in the figures correspond to points outside that region, i.e., are not acceptable.
As it can be seen in a general way in Figure 10, the overall cycle is more sensitive to irreversibilities in the discharge cycle than those of the charge (left versus right column). The irreversibilities in the components of the charge cycle (compressor, expander, and heat exchangers) do not directly affect the thermodynamic functions of the discharge cycle, such as the maximum power output. However, they produce a large reduction in the acceptable physical region, especially for large irreversibilities. This reduction in the acceptable physical region is shown as an increase in the length of the dashed lines and a reduction in the length of the solid lines. When the reduction of the acceptable physical region becomes very noticeable, the values of the parameter a c H E which maximize the power output and the HE efficiency might no longer be included in the physical region, decreasing the maximum available P o u t and Φ . Therefore, indirectly and through the decrease in the acceptable physical region, the irreversibilities of the HP affect the discharge thermodynamic functions. In Figure 10b it is observed how ϵ c , t H P reduce the acceptable physical region more abruptly than the external irreversibilities, ϵ H , L H P , shown in Figure 10c. This same behavior is reproduced for the HE in the 2nd column.

3.4. Sensitivity Analysis

In this section an analysis of the sensitivity of the maximum round-trip, discharge performance, power output and pumped heat with respect to losses will be performed. This will provide a landscape for irreversibilities relevance/influence for each figure of merit. For this, a reference configuration will be considered (design point shown in caption of Figure 6 and T H = 1000 K). In Figure 11 such effects on the maxima are depicted. The parameter α I r r * I r r × 100 is introduced as a measure of the irreversibilities variation with respect to the design point. I r r refers to each parameter and I r r * is the value at the reference/design point. In this way, at the design point: α = 0 ; and according Table 2 all ϵ coefficients vary from 1 ( α = 5 ) to 0.9 ( α = + 5 ). For the pressure drops ( ρ ) will be varied from 1 ( α = 2 ) to 0.95 ( α = + 3 ). Figure 11a shows the maximum round-trip efficiency when pressure drops vary from ideal conditions, down to 0.95 . Figure 11b shows the effects on the maximum round-trip efficiency due to the increment or decrement of all the efficiencies. In Figure 11c the maximum efficiency is depicted for variation of pressure drops and in Figure 11d the same is done for variations of each ϵ . This sensitivity analysis is done also for the maximum power density, P ¯ o u t , in Figure 11e for variation in the ρ ’s and in Figure 11f for variations in all ϵ ’s. Finally, the sensitivity of the maximum pumped heat is depicted in Figure 11g for variations of each ρ and in Figure 11h for variations in each ϵ .
It was found that the irreversibilities in the expander/turbine, ϵ t , are the most determinant in the drop of the maximum values for Φ , η and P o u t (see Figure 11b,d,f). Especially the losses in the expander of the charge cycle, ϵ t H P , are more determinant than those of the turbine in the discharge, ϵ t H E whose influence stem from the operation region reduction. It is also noticeable that the reduction in the maximum values of those functions with ϵ t H P over design point values is clearly not linear (see purple curves). From the Figure 11a,b, it can be seen that by reducing the irreversibilities with respect to the considered design point, maximum round-trip efficiencies over 0.4 could be obtained, which is a good result compared to previously reported values.
Notice that the maximum heat pumped strongly depends on the external irreversibilities of the discharge cycle, ϵ H H E and ϵ L H E (see black and blue curves in Figure 11h). Meanwhile, the influence of ϵ H H P is small, and as expected, there is no influence from ϵ L H P . The relevance of the HE in the pumped heat is due to the linking between both subsystems, stemming from the constraint Q H H P = Q H H E . This behavior is reproduced in the case of pressure drops; those of the HP do not influence the maximum Q H , but the influence of the HE is significant.

4. Discussion and Conclusions

In the present work, a PTES system based on a coupled Brayton-like HE and HP has been analyzed. This includes the analysis of 12 parameters associated with irreversibilities, eight for internal irreversibilities (pressure losses and expanders/turbines and compressors losses), and four for external ones (heat transfers with the reservoirs). One variable, a c H E , the pressure ratio of the compressors in the discharge cycle, is used as an independent variable to establish the operation regime. Meanwhile, a c H P is given by the constraint Q H H P = Q H H E .
Unlike in the case of energy storage in liquid media with variable temperature in the TES [13], here it is shown that in the present model it is not possible to simultaneously optimize power output and round-trip efficiency (see Figure 8 and Figure 9), mostly because the physically acceptable region is wider. In the case reported in [13], the physical region was strongly constrained by the thermal properties of the molten salts (lower and upper limits for the temperatures of the reservoirs). Here the absence of this constraint allows for a wider range of operation regimes in the Φ - P o u t curve. Additionally, increasing the temperature of the hot storage ( T H ) allows for better results in P, η , and Φ altogether. However, above a certain value, the increase of T H reduces the acceptable physical region, restraining the range of acceptable values for a c H E .
The internal irreversibilities, when compared with the external ones, show the following remarkable double effect:
  • The values of the internal irreversibilities considerably reduce the maximum achievable power output and round-trip efficiency with respect to the same values for the external irreversibilities.
  • The internal irreversibilities greatly reduce the range of acceptable values for a c H E ; affecting the available operation regimes. In this way, the irreversibilities of the discharge can influence significantly the thermodynamic functions of the charge by severely reducing its acceptable physical region.
Some points can also be remarked from the performed sensitivity analysis:
  • The expander/turbine is the most determinant element, especially the expander of the charge cycle is more decisive than the turbine of the discharge cycle.
  • Not all irreversibilities exhibit an effect on all thermodynamic functions. The irreversibilities of the charge cycle do not influence the functions of the discharge other than in the case when they considerably reduce the physical region.
  • By reducing irreversibilities it is possible to reach maximum values of the round-trip efficiency around 40% or even above. This is a significant result that should be considered for further studies on this type of model.
By establishing the most relevant components affecting the main energetic functions, this analysis might be used to focus on specific improvements in the turbomachinery. The cost of improving specific components might be weighted by its energetic relevance.
Future perspectives include analyzing the constraints imposed by the external sources (wind, photovoltaic, or others) in the inlet power in the charge phase, studying in detail the operation strategies considering realistic storage times, and including a more detailed analysis of the peculiarities of storage reservoirs and transient regimes.

Author Contributions

The authors contributed equally to the paper. All authors have read and approved the final manuscript.

Funding

This research was funded by University of Salamanca: 0218 463AB01; University of Salamanca: 18.KB.YF/463A.C.01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

J.G.A. acknowledges financial support from University of Salamanca, with contract number 0218 463AB01; and A.M. with contract number 18.KB.YF/463A.C.01.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Acronyms
COPCoefficient of performance
HEHeat engine
HPHeat pump
PTESPumped Thermal Energy Storage
RTERound-trip efficiency
TESThermal Energy Storage
Symbols
aTemperature ratio
CHeat capacity (J/Ks)
Q ˙ Heat flow (J/s)
WWork (J)
PPower (J/s)
P ¯ o u t Specific Power P / C w (K)
κ Dimensionless factor γ 1 γ
m ˙ Mass flow (kg/s)
pPressure (atm)
rPressure ratio
TTemperature (K)
Greek Letters
γ Adiabatic coefficient of ideal gas
ϵ Efficiency parameter
η Discharge efficiency
ν Coefficient of performance
ρ Pressure drops coefficient
Φ Round-trip efficiency
Subscript
wWorking fluid
cCompressor
HHot storage
LCold storage
tExpander/turbine

References

  1. Benato, A.; Stoppato, A. Pumped Thermal Electricity Storage: A technology overview. Therm. Sci. Eng. Prog. 2018, 6, 301–315. [Google Scholar] [CrossRef]
  2. Benato, A. Performance and cost evaluation of an innovative Pumped Thermal Electricity Storage power system. Energy 2017, 138, 419–436. [Google Scholar] [CrossRef]
  3. Gallo, A.B.; Simões Moreira, J.R.; Costa, H.K.M.; Santos, M.M.; Moutinho dos Santos, E. Energy Storage in the energy transition context: A technology review. Renew. Sust. Energy Rev. 2016, 65, 800–822. [Google Scholar] [CrossRef]
  4. Guney, M.S.; Tepe, Y. Classification and assessment of energy storage systems. Renew. Sust. Energy Rev. 2017, 75, 1187–1197. [Google Scholar] [CrossRef]
  5. Dumont, O.; Frate, G.F.; Pillai, A.; Lecompte, S.; De paepe, M.; Lemort, V. Carnot battery technology: A state-of-the-art review. J. Energy Storage 2020, 32, 101756. [Google Scholar] [CrossRef]
  6. Cabeza, L.F.; Solé, A.; Barreneche, C. Review on sorption materials and technologies for heat pums and thermal energy storage. Renew. Energy 2016, 110, 3–39. [Google Scholar] [CrossRef] [Green Version]
  7. Davenne, T.R.; Peters, B.M. An Analysis of Pumped Thermal Energy Storage With De-coupled Thermal Stores. Front. Energy Res. 2020, 8, 160. [Google Scholar] [CrossRef]
  8. Guo, J.; Cai, L.; Chen, J.; Zhou, Y. Performance evaluation and parametric choice criteria of a Brayton pumped thermal electricity storage system. Energy 2016, 113, 693–701. [Google Scholar] [CrossRef]
  9. Laughlin, R.B. Pumped thermal grid storage with heat exchange. J. Renew. Sust. Energy Rev. 2017, 9, 044103. [Google Scholar] [CrossRef] [Green Version]
  10. Esence, T.; Bruch, A.; Molina, S.; Stutz, B.; Fourmigué, J.F. A review on experience feedback and numerical modeling of packed-bed thermal energy storage systems. Sol. Energy 2017, 153, 628–654. [Google Scholar] [CrossRef]
  11. McTigue, J.D.; White, A.J.; Markides, C.N. Parametric studies and optimisation of pumped thermal electricity storage. Appl. Energy 2015, 137, 800–811. [Google Scholar] [CrossRef] [Green Version]
  12. Salomone-González, D.; González-Ayala, J.; Medina, A.; Roco, J.M.M.; Curto-Risso, P.L.; Hernández, A.C. Pumped heat energy storage with liquid media: Thermodynamic assessment by a Brayton-like model. Energy Conv. Manag. 2020, 226, 113540. [Google Scholar] [CrossRef]
  13. González-Ayala, J.; Salomone-González, D.; Medina, A.; Roco, J.M.M.; Curto-Risso, P.L.; Hernández, A.C. Multicriteria optimization of Brayton-like pumped thermal electricity storage with liquid media. J. Energy Storage 2021, 44, 103242. [Google Scholar] [CrossRef]
  14. Wang, L.; Lin, X.; Chai, L.; Peng, L.; Yu, D.; Chen, H. Cyclic transient behavior of the Joule-Brayton based pumped heat electricity storage: Modeling and analysis. Renew. Sust. Energy Rev. 2019, 111, 523–534. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Wang, Z. Comparative study on optimized round-trip efficiency of pumped thermaland pumped cryogenic electricity storages. Energy Conv. Manag. 2021, 238, 114182. [Google Scholar] [CrossRef]
  16. Thess, A. Thermodynamic Efficiency of Pumped Heat Electricity Storage. Phys. Rev. Lett. 2013, 111, 110602. [Google Scholar] [CrossRef] [Green Version]
  17. Chen, J.; Guo, J. Comment on “Thermodynamic Efficiency of PumpedHeat Electricity Storage”. Phys. Rev. Lett. 2016, 116, 158901. [Google Scholar] [CrossRef]
  18. 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] [Green Version]
  19. Gonzalez-Ayala, J.; Roco, J.M.M.; Medina, A.; Hernández, A.C. Carnot-Like Heat Engines Versus Low-Dissipation Models. Entropy 2017, 19, 182. [Google Scholar] [CrossRef] [Green Version]
  20. Gonzalez-Ayala, J.; Roco, J.M.M.; Hernández, A.C. Entropy generation and unified optimization of Carnot-like and low-dissipation refrigerators. Phys. Rev. E 2018, 97, 022139. [Google Scholar] [CrossRef] [Green Version]
  21. Roco, J.M.M.; Velasco, S.; Medina, A.; Hernández, A.C. Optimum performance of a regenerative Brayton thermal cycle. J. Appl. Phys. 1997, 82, 2735–2741. [Google Scholar] [CrossRef]
  22. Sánchez-Orgaz, S.; Medina, A.; Hernández, A.C. Thermodynamic model and optimization of a multi-step irreversible Brayton cycle. Energy Conv. Manag. 2010, 51, 2134–2143. [Google Scholar] [CrossRef]
  23. Ahmadi, M.H.; Ahmadi, M.A.; Pourfayaz, F.; Bidi, M. Thermodynamic analysis and optimization for an irreversible heat pump working on reversed Brayton cycle. Energy Conv. Manag. 2016, 110, 260–267. [Google Scholar] [CrossRef]
Figure 1. Components scheme of the Brayton-like HP charge cycle and the approximate T S diagram in the inset at the right. The processes 2 1 and 4 3 have been represented as isobaric (ideal), therefore they do not end at the point corresponding to the final state, because of the pressure drops. Hereinafter, realistic processes including losses will be shown.
Figure 1. Components scheme of the Brayton-like HP charge cycle and the approximate T S diagram in the inset at the right. The processes 2 1 and 4 3 have been represented as isobaric (ideal), therefore they do not end at the point corresponding to the final state, because of the pressure drops. Hereinafter, realistic processes including losses will be shown.
Entropy 23 01564 g001
Figure 2. Scheme of the discharge Brayton cycle and its corresponding T S diagram in the inset.
Figure 2. Scheme of the discharge Brayton cycle and its corresponding T S diagram in the inset.
Entropy 23 01564 g002
Figure 3. Two illustrative examples of the coupling of charge and discharge cycles in T S diagrams. Values considered for the temperatures of the reservoirs: T H = 1000 K, T L = 300 K. In (a) the case where a c H E = 1.8 and a c H P = 4.57, and in (b) the case a c H E = 1.2 and a c H P c = 5.17.
Figure 3. Two illustrative examples of the coupling of charge and discharge cycles in T S diagrams. Values considered for the temperatures of the reservoirs: T H = 1000 K, T L = 300 K. In (a) the case where a c H E = 1.8 and a c H P = 4.57, and in (b) the case a c H E = 1.2 and a c H P c = 5.17.
Entropy 23 01564 g003
Figure 4. Example of an endoreversible coupled cycle. The figures correspond to the following parameters: a c H E = 1.2 , a c H P = 5.47 , ϵ c H P , H E = ϵ t H P , H E = 1 , ρ H H P , H E = ρ L H P , H E = 1 , ϵ H H P , H E = ϵ L H P , H E = 0.95 , T H = 1000 K, and T L = 300 K.
Figure 4. Example of an endoreversible coupled cycle. The figures correspond to the following parameters: a c H E = 1.2 , a c H P = 5.47 , ϵ c H P , H E = ϵ t H P , H E = 1 , ρ H H P , H E = ρ L H P , H E = 1 , ϵ H H P , H E = ϵ L H P , H E = 0.95 , T H = 1000 K, and T L = 300 K.
Entropy 23 01564 g004
Figure 5. Round-trip efficiency in the endoreversible limit at maximum P o u t .
Figure 5. Round-trip efficiency in the endoreversible limit at maximum P o u t .
Entropy 23 01564 g005
Figure 6. Physically acceptable regions for the two parameters chosen as independent variables, r c H E and T H . Parameters quantifying irreversibilities at a representative design point are: [21,22,23]: ϵ c H P , H E = ϵ t H P , H E = 0.95 , ϵ H H P , H E = ϵ L H P , H E = 0.95 , and ρ H H P , H E = ρ L H P , H E = 0.98 . Green region: 0 Φ 1 . Blue: Combinations of T H and r c H E satisfying the temperature bounds in Table 1. The intersection between both represents the physically allowable region.
Figure 6. Physically acceptable regions for the two parameters chosen as independent variables, r c H E and T H . Parameters quantifying irreversibilities at a representative design point are: [21,22,23]: ϵ c H P , H E = ϵ t H P , H E = 0.95 , ϵ H H P , H E = ϵ L H P , H E = 0.95 , and ρ H H P , H E = ρ L H P , H E = 0.98 . Green region: 0 Φ 1 . Blue: Combinations of T H and r c H E satisfying the temperature bounds in Table 1. The intersection between both represents the physically allowable region.
Entropy 23 01564 g006
Figure 7. Evolution with the heat engine compressor temperature ratio, a c H E , and the equilibrium temperature of the hot storage, T H , of some key functions: (a) heat pump temperature ratio, (b) heat engine power output, P ¯ o u t P o u t / C w , (c) heat engine efficiency, η , and (d) round-trip efficiency, Φ . Irreversibilities parameters are like in Figure 6 and T L = 300 K.
Figure 7. Evolution with the heat engine compressor temperature ratio, a c H E , and the equilibrium temperature of the hot storage, T H , of some key functions: (a) heat pump temperature ratio, (b) heat engine power output, P ¯ o u t P o u t / C w , (c) heat engine efficiency, η , and (d) round-trip efficiency, Φ . Irreversibilities parameters are like in Figure 6 and T L = 300 K.
Entropy 23 01564 g007
Figure 8. (a) Input and output power versus discharge temperature ratio. (b) Round-trip efficiency, COP and discharge efficiency versus discharge temperature ratio. Black dots indicate the value of a c H E leading to maximum round-trip efficiency, the red dots those of maximum power output.
Figure 8. (a) Input and output power versus discharge temperature ratio. (b) Round-trip efficiency, COP and discharge efficiency versus discharge temperature ratio. Black dots indicate the value of a c H E leading to maximum round-trip efficiency, the red dots those of maximum power output.
Entropy 23 01564 g008
Figure 9. Parametric loops: in (a) round-trip efficiency versus  P ¯ i n , in (b) round-trip efficiency versus P ¯ o u t , in (c) round-trip efficiency versus the HP COP and in (d) round-trip efficiency versus the discharge efficiency. Black dots indicate the value of a c H E leading to maximum round-trip efficiency, and red dots those of maximum power output.
Figure 9. Parametric loops: in (a) round-trip efficiency versus  P ¯ i n , in (b) round-trip efficiency versus P ¯ o u t , in (c) round-trip efficiency versus the HP COP and in (d) round-trip efficiency versus the discharge efficiency. Black dots indicate the value of a c H E leading to maximum round-trip efficiency, and red dots those of maximum power output.
Entropy 23 01564 g009
Figure 10. Parametric plots of P ¯ o u t vs Φ over all available a c H E values ( a c H E increases clockwise). Solid lines indicate physically acceptable regimes meanwhile dashed curves denote the not available operation states. In each plot the corresponding irreversibility parameter varies taking the values shown in Table 2: (a) ρ H , L H P , (b) ϵ c , t H P , (c) ϵ H , L H P , (d) ρ H , L H E , (e) ϵ c , t H E , (f) ϵ H , L H E . Notice the limit of the physical region when the dashed and semi-transparent lines begin. In all figures T H = 1000 K and T L = 300 K.
Figure 10. Parametric plots of P ¯ o u t vs Φ over all available a c H E values ( a c H E increases clockwise). Solid lines indicate physically acceptable regimes meanwhile dashed curves denote the not available operation states. In each plot the corresponding irreversibility parameter varies taking the values shown in Table 2: (a) ρ H , L H P , (b) ϵ c , t H P , (c) ϵ H , L H P , (d) ρ H , L H E , (e) ϵ c , t H E , (f) ϵ H , L H E . Notice the limit of the physical region when the dashed and semi-transparent lines begin. In all figures T H = 1000 K and T L = 300 K.
Entropy 23 01564 g010
Figure 11. Sensitivity analysis of the maximum round-trip efficiency (a,b), HE efficiency (c,d), power output (e,f) and pumped heat (g,h) with respect to the reference/design irreversible configuration. The parameter α I r r * I r r × 100 measures the relative changes of each parameter with respect to the design point. See text for details.
Figure 11. Sensitivity analysis of the maximum round-trip efficiency (a,b), HE efficiency (c,d), power output (e,f) and pumped heat (g,h) with respect to the reference/design irreversible configuration. The parameter α I r r * I r r × 100 measures the relative changes of each parameter with respect to the design point. See text for details.
Entropy 23 01564 g011
Table 1. Bounds to be considered in the coupled charge-discharge system.
Table 1. Bounds to be considered in the coupled charge-discharge system.
Overall Cycle Constraints and Bounds
Heat conservation Q H H P = Q H H E T 2 H P T 1 H P = T 2 H E T 1 H E
Charge conditions 0 T 4 H P T 3 H P T L T H T 1 H P T 2 H P
Discharge conditions for T 3 H E T L T 4 H E T 3 H E T 2 H E T H
Discharge conditions for T 1 H E T L T 4 H E T 1 H E T 2 H E T H
Energy conservation 0 Φ 1 0 T 2 H E T 1 H E T 3 H E + T 4 H E T 2 H P T 1 H P T 3 H P + T 4 H P 1
Carnot maximum efficiency 0 η 1 T L / T H & 1 ν 1 T L / T H 1
Pressure and temperature ratios links a t H P = a c H P ρ H H P ρ L H P
Pressure and temperature ratios links a t H E = a c H E ρ H H E ρ L H E
Table 2. Typical values of irreversibilities.
Table 2. Typical values of irreversibilities.
IrreversibilitiesInterval
ϵ c , t H P , H E 0.85 , 1
ϵ H , L H P , H E 0.85 , 1
ρ H , L H P , H E 0.9 , 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pérez-Gallego, D.; Gonzalez-Ayala, J.; Calvo Hernández, A.; Medina, A. Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities. Entropy 2021, 23, 1564. https://doi.org/10.3390/e23121564

AMA Style

Pérez-Gallego D, Gonzalez-Ayala J, Calvo Hernández A, Medina A. Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities. Entropy. 2021; 23(12):1564. https://doi.org/10.3390/e23121564

Chicago/Turabian Style

Pérez-Gallego, David, Julian Gonzalez-Ayala, Antonio Calvo Hernández, and Alejandro Medina. 2021. "Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities" Entropy 23, no. 12: 1564. https://doi.org/10.3390/e23121564

APA Style

Pérez-Gallego, D., Gonzalez-Ayala, J., Calvo Hernández, A., & Medina, A. (2021). Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities. Entropy, 23(12), 1564. https://doi.org/10.3390/e23121564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop