razmi2020
razmi2020
Energy
journal homepage: www.elsevier.com/locate/energy
a r t i c l e i n f o a b s t r a c t
Article history: In the present article, thermodynamic, exergoeconomic, economic, and sustainability investigations of a
Received 25 April 2020 recently developed environmentally friendly hybrid absorption/recompression refrigeration cycle is
Received in revised form proposed to evaluate its feasibility for decision making and marketing. The proposed system is a novel
6 July 2020
hybridization of the conventional vapor compression and absorption cycles, wherein a booster
Accepted 28 July 2020
Available online 15 August 2020
compressor has been used between the generator and condenser of the single-effect absorption cycle to
improve its performance. Also, two separate multi-objective optimization models are developed using a
combination of the nondominated-storing-genetic algorithm (NSGA-II) and artificial neural network
Keywords:
Thermodynamic and thermoeconomic
(ANN) to address the optimum performance values concerning the objective functions. The obtained
analysis results approve that the proposed cycle is a promising concept from both thermodynamic and economic
Absorption/recompression refrigeration viewpoints. The results indicate that the system presents a coefficient of performance and exergy effi-
cycle ciency of 4.88 and 37.43% under the optimum working conditions. The overall rate of exergy destruction
Sustainability index of the system is 20.23 kW, and a sustainability index of around 1.53 can be achieved at a cooling capacity
Artificial neural network of 150 kW. The economic results indicate that the reference system has a payback period of 4.17 years,
Multi-objective optimization which is reduced to less than 4 years after doing the optimizations.
ANN
© 2020 Elsevier Ltd. All rights reserved.
NSGA-II
https://doi.org/10.1016/j.energy.2020.118559
0360-5442/© 2020 Elsevier Ltd. All rights reserved.
2 A.R. Razmi et al. / Energy 210 (2020) 118559
traditional absorption systems. (c) An improvement in COP, exergy advantage of this configuration, efficient operation at low heat
efficiency, and primary energy ratio. (d) A reduction in the initial source temperatures, indicating that the required heat source
investment cost in comparison to absorption cycles. temperature is reduced from 125 to 130 C in the conventional NH3/
The first configuration is when the compressor comes between H2O absorption cycles to 110e115 C in the proposed AR cycle [21].
the absorber and the evaporator. This design is the best option for In the second configuration, i.e. when a compressor is located
performance improvement and application of traditional absorp- between the absorber and generator, the cycle is driven by elec-
tion cycles in the lower heat source and ambient temperatures. Wu tricity rather than a heat source. This arrangement was initially
et al. [18] theoretically compared various working fluids for this presented by Osenbruck [22]; then, the cycle was developed by
arrangement, concluding that the maximum performance is ach- Altenkirch [23]. The generator of this configuration is capable of
ieved by using NH3/LiNO3 working fluid. Then, they substantiated extracting heat from low temperature fluids, while the absorber has
their findings with an experimental setup, representing that a COP the potential of releasing the heat of the absorption process to a
of 1.37 can be achieved at the evaporator and heat source tem- high-temperature fluid. Some of the major profits of this configu-
peratures of 25 C and 130 C by utilizing NH3/H2O as working ration are [24]: (a) Easy adjustment of heating and cooling capac-
fluid in the reference system [19]. As a result, this configuration can ities by solution composition; (b) A reduction in the unavoidable
be considered as an applicable heat pump for cold climates. heat losses in various processes; (c) A reduction in discharge tem-
Eventually, they investigated the annual performance of their perature and pressure of the compressor because of the lower
proposed system from technical and economic points of view un- saturation pressure of solution in comparison to the pure
der various heat source temperatures and compression ratios to refrigerant.
achieve an optimum operating condition [20]. It was concluded The third configuration uses a compressor between the low- and
that the obtained heat is markedly increased by employing a high-pressure stages. To develop this configuration, Ayala et al.
booster compressor between the absorber and the evaporator. [25,26] proposed theoretical and experimental analyses using NH3/
Accordingly, in spite the additional budget for purchasing the vapor LiNO3 working fluid, representing that the cooling and heating
compressor, the total capital cost of the system was remarkably capacities of this configuration can be simply and effectively
decreased due to a reduction in the size of the system. The second modulated. Their results disclosed that the provided heat at the
experimental analysis was done by We et al. to reveal another compressor exit is high enough to run an absorption subcycle.
A.R. Razmi et al. / Energy 210 (2020) 118559 3
Table 1
Comparison of different configurations of AR cycles.
X X
m_ i Xi ¼ m_ e Xe (2)
X X
m_ i hi þ Q_ ¼ m_ e he þ W_ (3)
X X
_
m_ i exi þ Ex _ W þ Ex
m_ e exe þ Ex _ D
Q ¼ (4)
here, subscripts i and e show the inlet and outlet streams of the
control volume and Ex_ D represents the exergy destruction of each
Fig. 2. Qualitative temperature-entropy diagram of the hybrid absorption/recom-
_ , and, Ex
control volume. Also, ex, Ex _ W are respectively specific,
pression system. Q
thermal, and power exergy rates, defining as follows [40]:
the refrigerant is saturated vapor and superheated vapor at the inlet ex ¼ ðh h0 Þ T0 ðs s0 Þ (5)
and outlet of the compressor (streams 1 and 2). Besides, the
refrigerant at the condenser and evaporator outlets (streams 3 and
_ ¼ 1 T0 Q_
Ex (6)
5) is respectively saturated liquid and vapor. Q
Tk k
The following assumptions are supposed to simulate the pro-
posed AR cycle [32e36]:
_ W ¼W
Ex _ (7)
k
1. The understudy cycle operates at a steady state condition.
where subscript 0 relates to the reference condition. It is worthy to
2. Variations of kinetic and potential energies are negligible.
3. The isentropic efficiencies of the compressor and pump are mention that the kinematic and potential exergies have been
ignored in the formulation [41]. Since chemical reactions do not
constant.
4. The condenser coil outlet (point 3) is a saturated liquid. happen within the cycle components, absorber fluid’s chemical
state is constant. Furthermore, considering the insignificant rate of
5. The outlet stream of the evaporator (point 5) is supposed to be
saturated vapor. chemical exergy in comparison to physical exergy, it is discarded in
the present research like most of the publications in the literature
6. The pressure losses associated with frictional effects in pipes
[42e45]. To provide a more comprehensive perspective from the
and between components are ignored.
exergetic performance of the components, second law efficiency
are expressed as follows [46]:
3. System simulation _
Ex _
Ex
hex;k ¼ _ p ¼ 1 _ D (8)
Exf Exf
The governing equations for thermodynamic and exer-
goeconomic analyses of the reference cycle are proposed in this _ p andEx
_ denote the exergy rates of the product and fuel in
here Ex f
section. The primary assumptions and design conditions used for
the components. The product exergy is the desired exergy flow
the simulation of the current study are outlined in Table 2.
coming out from each component, whereas the fuel is the
consumed exergy within the component for providing that prod-
Table 2 uct. The governing energy balance and fuel-product equations for
The assumed design condition for the system simulation [30,37,38]. different equipment of the proposed cycle are tabulated in Tables 3
Foundation parameter Value Unit and 4.
To assess the proposed system from the sustainable perspective,
Absorber temperature 35 C
Ambient pressure 1.01 bar a sustainability index (SI) is considered in this study. In fact, SI is a
Ambient temperature 25
C determinative parameter that connects the exergy analysis to the
Annual operating hours 4320 Hr environmental aspect. It states that how decreasing the total exergy
Cost of electricity [37] 0.033 $/kWh destruction of the system in the energy conversion processes re-
Evaporator capacity 150 kW
Evaporator temperature 10
C
sults in a drop in the environmental impact of system [47,48]:
Generator temperature 65 C
Heat exchanger effectiveness 0.9 eee 1
SI ¼ (9)
Inlet temperature of the absorber coil 30 C DP
Inlet temperature of the evaporator coil 18 C
Interest rate 15 %
here DP is an indicator of depletion number, defining as the pro-
Isentropic efficiency of compressor 88 %
Isentropic efficiency of pump 85 % portion of total exergy destruction to the inlet exergy. Hence, the
Maintenance factor 1.06 eee more the sustainability index, the less the environmental effect of
Outlet temperature of the absorber coil 40 C the reference system.
Outlet temperature of the evaporator coil 13 C The COP and exergy efficiency of the reference system are ob-
Service year of the system 20 eee
tained as follows [30,49,50].
A.R. Razmi et al. / Energy 210 (2020) 118559 5
Table 3
Energy balance equations of the cycle components.
Component Equation
purchase cost of the kth component, and tis the annual working
Table 4
Fuel-product definition of the cycle components.
time of the cycle per hour. CRF is defined as [54]:
Table 5
Cost functions and purchase cost of the cycle components [32,38].
Table 6
Cost balance and auxiliary equations of the cycle components.
Absorber C_ d C_ c ¼ C_ 5 þ C_ 11 C_ 6 þ Z_ Abs C_ 5 þ C_ 11 C_
cc ¼ 0 ; ¼ 6
E_ 5 þ E_ 11 E_ 6
Compressor C_ 2 C_ 1 ¼ C_ W;Com þ Z_ Com cw ¼ 0:033 $=kWh[37]
Evaporator C_ b C_ a ¼ C_ 4 C_ 5 þ Z_ Evp c4 ¼ c5 ; ca ¼ 0
Expansion valve 1 C_ ¼ C_ þ Z_ Ev
4 3 1
are considered for minimizing the total cost rate (C_ tot ). Total cost 4.1. Model validation
rate is achieved as the sum of the capital investment, operating, and
maintenance costs, as well as the cost related to the exergy To verify the simulation results for the hybrid AR cycle, a com-
destruction [58]: parison of the obtained data is considered with those of reported by
X X Riffat et al. [29]. This comparison has been outlined in Table 7. As it
C_ tot ¼ Z_ k þ C_ D;k (22) can be seen, there is a reasonable deviation between the current
study and the similar publication in the literature, meaning that the
model is reliable and authentic.
The net present value (NPV) method is applied to analyze the The thermodynamic properties of the streams in the reference
economic feasibility of the present cycle, which can address the cycle are listed in Table 8. A computer program was developed in
payback period of the invested capital budgeting and the cycle the EES software and reported data in Table 8 are extracted from
profit during its economic book life [59]. According to NPV method, the system thermodynamic simulation via simultaneous solving of
the net cash flow at the end of year is converted to its present value energy, exergy, and mass conversion equations.
based on the interest rate as follows [60]: The results associated with the exergy analysis, as well as the
specific exergy costs of product and fuel, cost of exergy destruction,
X
20 X
NPV ¼ Yð1 þ iÞN CCapital (23)
N¼0 Table 7
Validation of the simulation results with the literature.
here, N and i express the year and interest rate. Y indicates the net
Parameters Present work Ref. [29] Deviation (%)
cash flow at the end of the year. According to Table 1, an interest
rate of 15% and 20 years of service time have been considered for COP 4.17 4.21 0.95
the proposed system. QGen-Con (kW) 129.7 128.5 0.93
QEvp (kW) 100 100 0.00
QAbs (kW) 124.6 123.7 0.73
4. Results and discussion TGen ( C) 60 60 0.00
Wcom (kW) 23.98 23.74 1.01
The simulation results of the hybrid absorption/recompression Xlean (%) 53.28 52.21 2.05
Xrich (%) 58.85 57.21 2.86
cycle are reported in this section.
A.R. Razmi et al. / Energy 210 (2020) 118559 7
Table 8
The result of thermodynamic simulation for the system streams.
purchased equipment costs, exergoeconomic factor, and relative the exergy destruction streams of equipment.
cost difference for the system components are outlined in Table 9. The performance indicators and technical parameters of the
The most quantity of capital investment cost rate associates with system obtained from the thermodynamic and economic analyses
the absorber, followed by the evaporator, and the generator- are presented in Table 10. A vapor compressor with a compression
condenser components. Among the all, the compacted generator- ratio of 7.52 is needed to address the required heat power of the
condenser, owing to internal irreversibilities arisen at direct heat compacted generator-condenser (188 kW) for reaching the
exchange between the water vapor and LiBr-water solution, is designed generator temperature (65 C). As a result, 31.05 kW po-
responsible for about 52% of total exergy destruction in the system. wer is consumed in the compressor to deliver 150 kW refrigerating
As mentioned in Table 6, the specific exergy cost of the inlet water capacity in the evaporator with a COP of around 4.83. Moreover,
to the evaporator and the absorber (streams a and c) is supposed to 47.2 and 181 kW heat are respectively transferred in the heat
be zero. The specific exergy cost of cooling as the system product is exchanger and the absorber. The exergy analysis reveals that about
103.1 $/GJ. Also, the exergoeconomic factor of the second expansion 20.23 kW exergy is destructed in different processes and compo-
valve and the pump have the maximum values, meaning that they nents and the proposed cycle operates with an exergy efficiency of
have negligible cost of destruction. On the other hand, the lowest 34.84%. Also, the cycle has a sustainability index of around 1.53
values of exergoeconomic factor belong to the first expansion valve, which is a reasonable rate for refrigeration cycles. The economic
generator-condenser, compressor, and absorber. The calculated analysis denotes that an investment cost of around 24,586 $ is
values for the mentioned equipment indicate that the rate of cost needed to install a combined absorption/recompression refrigera-
destruction is higher than the capital cost rate. The total cost rates tion system with a refrigerating capacity of 150 kW. The reference
of the investment and the exergy destruction of the reference cycle cycle has a payback period of 4.17 years and 27,702 $ can be saved
are 0.96 and 3.35 $/hr, respectively, resulting in a total exer- during its service time (20 years). Finally, the exergoeconomic
goeconomic factor of 22.34%. It means that the cost of exergy analysis shows that the total cost rates of destruction and invest-
destruction is responsible for 77.66% of system cost rate, revealing ment are respectively 3.35 and 0.96 $/hr, leading to a total cost rate
the importance of applying the exergoeconomic analysis for the of 4.31 $/hr.
reference cycle. Thus, the use of expensive equipment that reduces
the cost destruction rate and enhances the investment cost rate,
4.3. Parametric analysis
improves the cycle performance from exergoeconomic perspective.
In addition, the second expansion valve has a low relative cost
In order to make a thorough assessment, the effect of the sig-
difference that indicates the negligible difference between the cost
nificant parameters on the system performance is examined in this
per exergy unit of its product and fuel, as well as the superior
section. Among all, the generator and evaporator temperatures are
benefit.
the most crucial factors that can substantially affect system oper-
A Sankey exergy flow diagram is also prepared and illustrated in
ation. Thus, the adjustment of the generator and evaporator tem-
Fig. 3 to give a better understanding of cycle operation for readers.
peratures are considered as the base of parametric analysis.
It is worthy to mention that the width of the vectors is comparative
The impact of the generator and evaporator temperatures on the
to the exergy values of the cycle streams and the dark arrows show
heat transfer rate of the compacted generator-condenser and the
Table 9
The result of exergy and exergoeconomic analyses for system components.
Component _ ½kW
Exf
_ p ½kW
Ex _ D ½kW
Ex hex ½% cf ½$=GJ cp ½$=GJ cD ½$=hr Z_ ½$=hr f ½% r ½%
Absorber 8.485 5.962 2.523 70.26 112.5 176.6 1.0220 0.3543 25.74 36.30
Compressor 31.05 29.08 1.974 93.64 9.167 10.03 0.0651 0.0256 28.21 8.639
Evaporator 7.867 4.859 3.008 61.76 74.69 102.9 0.8116 0.3207 65.33 37.45
Expansion valve 1 0.867 0.783 1.651 90.42 67.68 74.96 0.4022 0.0003 0.075 10.75
Expansion valve 2 28.11 28.11 0.000 100 54.78 54.81 0.0000 0.0039 100 0.071
Generator-condenser 31.04 20.53 10.51 66.15 67.68 15.99 2.5600 0.1317 4.894 76.37
Heat exchanger 3.789 3.220 0.569 84.97 54.78 75.05 0.1123 0.1227 52.21 27.01
Pump 0.002 0.000 0.002 2.86 9.167 eeee 0.0001 0.0046 98.44 99.95
8 A.R. Razmi et al. / Energy 210 (2020) 118559
Fig. 3. The Sankey exergy flow diagram for the hybrid absorption/recompression refrigeration system.
Table 10
The important result of the reference system.
Cr 7.52 e
C_D; tot
3.35 $/hr
C_ tot 4.31 $/hr
cp 76.49 $/GJ
COP 4.83 eee
_ D; tot
Ex 20.23 kW
hex 34.84 %
ftot 22.34 %
Q_Gencon
188 kW
Q_ Evp 150 kW
Q_Abs
181 kW
Q_ HXCH 47.2 kW
SI 1.53 eee
NPV20 27,702 $
tp 4.17 year
W _ 31.05 kW
Com
W _ 0.002 kW
P
XLean 52.2 %
XRich 56.3 %
Ztot 24,586 $ Fig. 4. Variation in the required heat and power of the generator and compressor with
Z_ tot 0.96 $/hr the generator and evaporator temperatures.
Fig. 5. Variation in the COP and compression ratio of the compressor with the
generator and evaporator temperatures.
Fig. 7. Variation in the exergy destruction ratio of the components with the generator
temperature.
Fig. 9. Variation in the sustainability index and total exergoeconomic factor with the
generator and evaporator temperatures.
5. Multi-objective optimization
subject to:
50 TGen 85
5 TEvp 10 (25)
32 TAbs 38
Fig. 11. Variation in the payback period and cost of product with the generator and
evaporator temperatures.
5.1. Optimization methodology and accuracy check
Fig. 15. Distribution of the Pareto frontier diagram for COP and payback period (model
(1)).
The Pareto frontier solution of the model (1) for optimizing the
COP and payback period is shown in Fig. 15.
Although all of the points shown in the Pareto frontier diagram
reveal the optimal solutions, point B is the ideal solution of the
multi-criteria optimization for model (a), considering both the
optimized solutions of the payback period and COP, simultaneously.
At this ideal point, the COP and payback period are respectively
4.88 and 3.97 years. Moreover, point A with a COP of 4.44 and a
payback period of 3.84 years indicates the best optimal solution
Fig. 14. Accuracy verification of the ANN model for a) COP, b) exergy efficiency, and c) from the economic standpoint. Alternatively, point C shows the
payback period. best optimal solution from the perspective of the COP, where the
A.R. Razmi et al. / Energy 210 (2020) 118559 13
Fig. 17. Distribution of the Pareto frontier diagram for the exergy efficiency and
payback period (model (2)).
However, Fig. 16 (b) represents that the optimal points for the
generator temperature are located between 68TGen72; most of
which are at the generator temperature of 71 C. Finally, as dis-
cussed in section 4.3, Fig. 16 (c) confirms the better performance of
the cycle in higher evaporator temperatures from the viewpoint of
the COP and economic. It is worthy to mention that designing a
refrigeration cycle with the evaporator temperature of more than
10 C is not so reasonable and not considered in the optimization.
Also, the Pareto frontier solution and scatter distribution of the
model (b) for optimizing the exergy efficiency and payback period
are presented in Figs. 17 and 18. Like the previous model, three
different points D, E, and F are considered in the Pareto frontier of
the model (2) to analyze the alterations of the exergy efficiency and
the payback period. Between all optimum points, the point E cor-
responds the best results for both exergy efficiency and the payback
period. Also, the point F can be preferred when the exergy effi-
ciency is the only objective function. In contrary, the point D rep-
resents the ideal solution for a single objective optimization based
on the payback period. Contrary to the scatter distribution of the
model (1), Fig. 18 (a) exhibits that there is a brilliant agreement for
the absorber temperature in the model (2), which should be
retained at the highest value (38 C). The more the generator
temperature, the less the payback period. Also, the exergy effi-
ciency is greater at the generator temperatures of higher than the
critical point (65 C); so, there is a good agreement between the
objective functions of the model (2) in higher generator tempera-
tures, as shown in Fig. 18 (b). In contrast to COP, the exergy effi-
ciency is greater at lower evaporator temperatures, whereas the
payback period increases. As a result, the scatter distribution of the
evaporator for the model (2) is distributed in a wide range and with
a more focus on the temperature of around 9.5 C.
Finally, the values of the optimum parameters and objectives for
Fig. 16. Scatter distribution for the important parameters of model (1); a) absorber
both models and points A, B, C, D, E, and F are indicated in Table 11.
temperature, b) generator temperature, and c) evaporator temperature.
6. Conclusions
payback period and COP are respectively 4.08 and 4.99. For a better
understanding of the key factors, the scatter distribution of model
Refrigeration systems play a crucial role in buildings’ energy
(a) is shown in Fig. 16.
saving and greenhouse gas emission reduction. The vapor
Fig. 16 (a) shows that the optimum points for the absorber
compression and absorption systems are the conventional kinds of
temperature are distributed at all ranges between 32 C and 38 C.
chillers, operating with opposite properties in terms of
14 A.R. Razmi et al. / Energy 210 (2020) 118559
Table 11
The optimum values for six Pareto frontier optimal solutions.
waste-driven CHPeORC plant with exhaust heat recovery. Energy Convers geostatistical modeling using genetic algorithm and direct search methods.
Manag 2019;187:512e22. J Petrol Sci Eng 2008;63:34e42.
[60] Najafi A, Jafarian A, Darand J. Thermo-economic evaluation of a hybrid solar- [64] Sohani A, Sayyaadi H, Hoseinpoori S. Modeling and multi-objective optimi-
conventional energy supply in a zero liquid discharge wastewater treatment zation of an M-cycle cross-flow indirect evaporative cooler using the GMDH
plant. Energy Convers Manag 2019;188:276e95. type neural network. Int J Refrig 2016;69:186e204.
[61] Esfahani IJ, Yoo CK. Exergy analysis and parametric optimization of three [65] Jiang H, Xi Z, Rahman AA, Zhang X. Prediction of output power with artificial
power and fresh water cogeneration systems using refrigeration chillers. neural network using extended datasets for Stirling engines. Appl Energy
Energy 2013;59:340e55. 2020;271:115123.
[62] Sohani A, Sayyaadi H, Mohammadhosseini N. Comparative study of the con- [66] Behzadi A, Arabkoohsar A, Gholamian E. Multi-criteria optimization of a
ventional types of heat and mass exchangers to achieve the best design of biomass-fired proton exchange membrane fuel cell integrated with organic
dew point evaporative coolers at diverse climatic conditions. Energy Convers rankine cycle/thermoelectric generator using different gasification agents.
Manag 2018;158:327e45. Energy 2020;201:117640.
[63] Maschio C, Campane A, Jose D. A framework to integrate history matching and