Graphical Abstract Figure
Graphical Abstract Figure
Close modal

Abstract

Cascade utilization is the most promising solution for efficient recovery of industrial waste heat. In this paper, a combined cooling and power (CCP) system including an organic Rankine cycle (ORC) and an absorption refrigeration cycle (ARC) is proposed, aiming to realize the cascade utilization of efficient utilization of low-grade flue gas. With the objective of maximizing the equivalent specific net power output, a particle swarm optimization (PSO) algorithm is employed to optimize the operation parameters of the proposed system and the composition and mass fractions of zeotropic working fluid. The results indicate that the proposed system can realize efficient cascade utilization for flue gas waste heat energy. The optimized proposed system achieves a specific net power output of 6.09 kW/kgFG of a specific refrigeration capacity of 27.65 kW/kgFG. Additionally, heat exchanger exergy losses account for 73.7% and 53.6% of the total exergy losses in the ORC and ABR subsystems, respectively. Therefore, optimizing the heat exchanger equipment is a feasible approach to further enhance the thermodynamic performance of the proposed system.

1 Introduction

Low-grade thermal energy (below 230 °C) accounts for about 60% of industrial waste heat and has significant recycling potential [1]. However, single utilization methods such as power generation, refrigeration cannot achieve deep utilization of this low-grade heat energy [2]. Therefore, studying the cascade utilization of industrial waste heat, particularly low-grade waste heat, is highly significant.

Efficient approaches to harness low-grade industrial waste heat include enhancing heat upgrading, generating electricity, and refrigeration [3]. However, these techniques are often limited by the temperature range of the waste heat, making complete utilization a challenge. Dividing the heat source into distinct temperature regions and utilizing them in a cascading manner is a common approach to address suboptimal utilization of waste heat and incomplete energy conversion [4]. For example, Mohammadi et al. [5] established a combined cooling, heating, and power system comprising a gas turbine, an organic Rankine cycle (ORC), and an absorption refrigeration cycle (ARC), achieving an exergy efficiency of 68%. Zhou et al. [6] introduced a cascade system combined with an ARC and ORC, investigating how critical parameters impact the system's thermodynamic performance. However, most of these studies focused only on the effect of key parameters on the system, neglecting the optimization of the system as a whole.

Additionally, since electricity is the highest-grade energy source, ORC power generation is highly valued in the field of industrial waste heat recovery. For ORC systems, most of the studies have focused on selecting the optimal working fluid. Emadi et al. [7] utilized genetic algorithms and neural networks to optimize 20 working fluids for a solid oxide fuel cell waste heat recovery system, revealing that R601 emerged as the optimal choice. Wang et al. [8] compared ORC and proton exchange membrane fuel cells using different working fluids, finding that a zeotropic mixture of R245fa/R123 (0.6/0.4) demonstrated the best thermodynamic performance. Wang et al. [9] employed a multi-objective genetic algorithm to optimize the zeotropic working fluids for a dual-loop ORC system using flue gas waste heat. The results indicated that cyclohexane/cyclopentane (0.2/0.8) and butane/pentane (0.65/0.35) performed best in the high-temperature and low-temperature loops, respectively. In these studies, zeotropic working fluids are considered an effective option for improving the thermodynamic performance of ORC systems. However, most existing studies focus on optimizing the mass fraction of a zeotropic working fluid with a predetermined composition. Few studies have simultaneously optimized both the composition and the mass fraction of the zeotropic working fluid.

In this study, a combined cooling and power (CCP) system is proposed for the deep utilization of low-grade industrial waste heat, including an ORC with zeotropic working fluids and an ABR. A particle swarm optimization (PSO) algorithm is employed to optimize its operating parameters as well as the composition and mass fraction of the ORC working fluid. Based on this, the thermodynamic performance of the optimized system is analyzed.

2 System Description

The proposed system is illustrated in Fig. 1. It comprises an ORC subsystem for power generation and an ABR subsystem for refrigeration. In this system, the heat source (flue gas) is divided into two sections: Sec. 1, the high-temperature section, serves as the heat source for the ORC subsystem, while Sec. 2, the low-temperature section, is used as the heat source for the ABR subsystem.

Fig. 1
Conceptual diagram of the CCP system
Fig. 1
Conceptual diagram of the CCP system
Close modal

In this study, the working fluid of the ORC subsystem is considered as a binary zeotropic mixture. The specific components and mass fractions are obtained from the optimization (specific optimization method is given in Sec. 3), and the thermophysical properties of each alternative component (from Refprop 9.1) are listed in Table 1. Ammonia–water solution and lithium bromide–water solution are widely employed as working fluids in ABR. Compared to lithium bromide–water solution, ABR with an ammonia–water solution as the working fluid achieves lower refrigeration temperatures. Therefore, ammonia–water solution is selected as the working fluid for the ABR subsystem.

Table 1

Thermophysical properties of alternative components

First componentsCritical temperature (K)Critical pressure (MPa)Second componentsCritical temperature (K)Critical pressure (MPa)
R601469.73.37R12385.124.14
R601a460.353.38R124395.433.62
R113487.213.39R1234ze382.513.63
R123456.833.66R152a386.414.52
R141b477.54.21R236fa398.073.20
R245ca447.573.94R227ea374.902.93
First componentsCritical temperature (K)Critical pressure (MPa)Second componentsCritical temperature (K)Critical pressure (MPa)
R601469.73.37R12385.124.14
R601a460.353.38R124395.433.62
R113487.213.39R1234ze382.513.63
R123456.833.66R152a386.414.52
R141b477.54.21R236fa398.073.20
R245ca447.573.94R227ea374.902.93

3 Methodology

3.1 Evaluation Indicators.

The thermodynamic model of the proposed system is based on the energy conservation, as detailed in Refs. [10,11], and some evaluation indicators are listed in this section. The thermal efficiency of the ORC subsystem is calculated by
(1)
where Wexp and Wp denote the power output of the turbine and the power consumption of the pump, respectively, and Qeva denotes the heat exchange rate of evaporator.
The specific net power output of the ORC is defined as
(2)
where m˙FG denotes the mass flowrate of flue gas.
For ARC, the refrigeration efficiency is the ratio between the unit cooling capacity of the evaporator and the unit heat load of the generator:
(3)
where Qrf and Qgen are refrigeration capacity and the heat exchange rate of generator, respectively.
In addition, in order to evaluate the proposed system in a combined manner, the refrigeration capacity of the ABR subsystem is converted to electricity as follows:
(4)
where COPVC is the refrigeration efficiency of a standard electrically driven vapor compression system.
Therefore, the equivalent specific net power output of the proposed system is
(5)
The exergy represents the portion of energy within a system that is available for power generation, while also indicating the level of irreversibility present in the system. The exergy of a substance at any state can be mathematically expressed as follows:
(6)
where m˙ denotes the mass flowrate, h and s are the specific enthalpy and specific entropy, respectively, and h0, T0, and s0 are the specific enthalpy, temperature, and specific entropy at ambient conditions, respectively.
The balance equation for the exergy of a system or equipment can be expressed as
(7)
where Exin, Exout, and Exloss are the exergy input, exergy output, and exergy loss of the equipment, respectively.
Exergy efficiency represents the ratio of the output exergy of the system to the input exergy of the system:
(8)
where Exin,system and Exout,system are the exergy input and exergy output of the system, respectively.

3.2 Optimization Method.

In this study, the PSO algorithm is employed to achieve global optimization of the proposed CCP system. The objective of the optimization is to maximize the equivalent net power output of the system, and objective function is described as
(9)
where X denotes the decision-making vector composed of decision-making variables, and it can be described as
(10)
where yAi and yBi denote the first and second components of the ORC working fluid, respectively, and z represents the mass fraction of the first component. As an example, the vector of the first alternative component is [R601, R601a, R113, R123, R141b, R245ca]. When yAi = 1, the first component is determined to be R601. The same applies to the second component, allowing for the simultaneous optimization of the zeotropic working fluid components and mass fractions. With this clarification, the range of values for decision-making vector X is listed in Table 2.
Table 2

The lower and upper bounds of the decision variables

ItemLower boundUpper bound
The first component of working fluid, yAi16
The second component of working fluid, yBi16
Mass fraction of the first component, z0.020.98
Evaporation temperature (ORC, °C), Teva80150
Condensation temperature (ORC, °C), Tcon2040
Superheat (ORC, °C), ΔTsh0.55
Refrigeration temperature (ABR, °C), Trf−205
Generator inlet temperature (ABR, °C), Tfeed2560
Flue gas outlet temperature (°C), Tout5080
Generation temperature (ABR, °C), Tgen90107
ItemLower boundUpper bound
The first component of working fluid, yAi16
The second component of working fluid, yBi16
Mass fraction of the first component, z0.020.98
Evaporation temperature (ORC, °C), Teva80150
Condensation temperature (ORC, °C), Tcon2040
Superheat (ORC, °C), ΔTsh0.55
Refrigeration temperature (ABR, °C), Trf−205
Generator inlet temperature (ABR, °C), Tfeed2560
Flue gas outlet temperature (°C), Tout5080
Generation temperature (ABR, °C), Tgen90107
Additionally, for practical purposes, all heat exchangers in the proposed system must have a pinch point temperature difference greater than 3 °C. This requirement is incorporated into the optimization process through a penalty function denoted as
(11)

The whole flowchart of the PSO is shown as Fig. 2.

Fig. 2
Flowchart of PSO

4 Results and Discussion

4.1 Optimization Performance.

The base design parameters of the proposed system are listed in Table 3. The variation curves of the PSO fitness values under these parameter settings are shown in Fig. 3, and the optimal design parameters are listed in Table 4.

Fig. 3
System single objective global optimization process convergence diagram: (a) ORC and (b) ARC
Fig. 3
System single objective global optimization process convergence diagram: (a) ORC and (b) ARC
Close modal
Table 3

Base design parameters of the proposed system

ItemsUnitValues
Flue gas inlet temperature°C150
Flue gas middle temperature (ABR inlet)°C110
Cold water temperature°C20
Cold water temperature rise°C15
Flue gas mass flowratekg/s1
Isentropic efficiency of the pump0.75 [12]
Isentropic efficiency of the turbine0.8 [12]
ItemsUnitValues
Flue gas inlet temperature°C150
Flue gas middle temperature (ABR inlet)°C110
Cold water temperature°C20
Cold water temperature rise°C15
Flue gas mass flowratekg/s1
Isentropic efficiency of the pump0.75 [12]
Isentropic efficiency of the turbine0.8 [12]
Table 4

The optimal design parameters of the proposed system

ORCARC
ParametersValuesParametersValues
The first componentR141bRefrigeration temperature (°C)−5.9
The second componentR124Feed temperature (°C)27.3
The first component's mass fraction62.5%Flue gas outlet temperature (°C)51.4
Evaporation temperature (°C)136.3Generation temperature (°C)93.4
Evaporation pressurea (kPa)2283.1Ammonia's concentration in strong solutiona48.3%
Condensation temperature (°C)26.8The mass fraction of ammonia in dilute solutiona32.8%
Condensation pressure (kPa)188.8
Degree of superheat (°C)3.8
ORCARC
ParametersValuesParametersValues
The first componentR141bRefrigeration temperature (°C)−5.9
The second componentR124Feed temperature (°C)27.3
The first component's mass fraction62.5%Flue gas outlet temperature (°C)51.4
Evaporation temperature (°C)136.3Generation temperature (°C)93.4
Evaporation pressurea (kPa)2283.1Ammonia's concentration in strong solutiona48.3%
Condensation temperature (°C)26.8The mass fraction of ammonia in dilute solutiona32.8%
Condensation pressure (kPa)188.8
Degree of superheat (°C)3.8
a

Not the decision-making variables.

Based on the optimization results, the most suitable working fluid for the ORC subsystem is the R141b/R124 with a mass fraction of 0.625:0.375. Under this configuration, the ORC subsystem exhibits a specific net power output of 6.09 kW/kgFG (net power output per unit mass flowrate of flue gas), accompanied by a thermal efficiency of 14.07%. Additionally, when operating at a refrigeration temperature of −5.9 ℃, the ARC subsystem demonstrates a specific refrigeration capacity of 27.65 kW/kgFG and a COP of 0.4453. Consequently, the overall system showcases an exceptional equivalent specific output power amounting to 11.45 kW/kgFG.

In addition, to demonstrate the advantages of zeotropic working fluids for the ORC subsystem of the proposed system, the ORC subsystem with zeotropic working fluid was compared with the ORC subsystem using a single-component working fluid, as shown in Fig. 4. It is worth noting that the ORC subsystem with a single-component working fluid employs the same optimization methodology (PSO) to ensure a fair comparison. Using zeotropic working fluids in the ORC results in an increase in net power output by 0.43 kW/kgFG and thermal efficiency by 0.99% compared to using R141b alone. Furthermore, when compared to an ORC employing R124 as its working fluid, the use of zeotropic working fluid leads to a significant improvement in net power output by 1.53 kW/kgFG and thermal efficiency by 3.53%. These findings demonstrate that zeotropic mixtures substantially enhance the thermodynamic performance of the ORC subsystem.

Fig. 4
Performance comparison of the ORC subsystem using single-component working fluid and zeotropic working fluid: (a) net power output and (b) thermal efficiency
Fig. 4
Performance comparison of the ORC subsystem using single-component working fluid and zeotropic working fluid: (a) net power output and (b) thermal efficiency
Close modal

4.2 Exergy Analysis.

By focusing on the exergy flow within the system, the exergy analysis provides valuable insights into the distribution of “energy quality” losses in the thermodynamic system. The input exergy of the system originates from industrial flue gas and certain components, resulting in a total input exergy of 25.553 kW. The output exergy primarily comprises power generated by ORC and cooling capacity provided by ARC, totaling 11.230 kW. Consequently, the exergy efficiency stands at 43.9%. A comprehensive examination of exergy flow and loss within this system is presented below to provide a deeper understanding.

Figure 5 illustrates the exergy flow of the proposed CCP system. As depicted in the figure, the system's total input exergy amounts to 25.553 kW, with flue gas waste heat contributing 22.192 kW (86.8% of the total). The remaining exergy is supplied by pumps, distillers, and the non-freezing fluid in the ARC evaporator. After passing through the heat exchanger, the industrial flue gas delivers 9.535 kW of exergy to the power generation system and 8.816 kW to the refrigeration system. However, as this exergy traverses various heat exchangers, pumps, throttling valves, and other equipment, losses occur amounting to 14.313 kW (56.1% of the total input). The ORC and ARC exhibit exergy efficiencies of 53.2% and 35.7%, respectively.

Fig. 5
Exergy flow diagram of the proposed system
Fig. 5
Exergy flow diagram of the proposed system
Close modal

Figure 6 illustrates the proportion of exergy loss in the primary equipment of the system. The ORC contributes 39.15% of the overall exergy loss, while the ARC accounts for 60.85%. Notably, heat exchanger equipment plays a dominant role in causing exergy losses. Specifically, within the ORC, the condenser and evaporator contribute to 73.7% of the total exergy loss. Similarly, heat exchanger equipment in the CCP system is responsible for 56.0% of the total exergy loss. Therefore, further enhancements for the proposed system should primarily focus on optimizing and designing more efficient heat exchangers.

Fig. 6
Proportion of equipment exergy loss of the proposed system
Fig. 6
Proportion of equipment exergy loss of the proposed system
Close modal

5 Conclusions

This study introduces an innovative CCP system that effectively harnesses waste heat from industrial processes, specifically utilizing flue gas waste heat through a cascade utilization approach. A PSO algorithm is employed to optimize the proposed system, specifically focusing on the composition and mass fraction of the zeotropic working fluid in the ORC subsystem. The main findings are as follows:

  1. The optimized ORC subsystem of the proposed system achieves a specific net power output of 6.09 kW/kgFG and a thermal efficiency of 14.07%. Meanwhile, the ABR subsystem attains a specific refrigeration capacity of 27.65 kW/kgFG and a COP of 0.4453.

  2. The zeotropic working fluid significantly enhances the thermodynamic performance of the ORC subsystem. Compared to ORC systems using pure R141b and pure R124, the ORC subsystem of the proposed system experiences a net power output increase of 0.43 kW/kgFG and 1.53 kW/kgFG, respectively.

  3. In the ORC subsystem, heat exchanger exergy losses account for 73.7% of the total exergy losses. In the ABR subsystem, heat exchanger exergy losses account for 53.6% of the total exergy losses. Therefore, further enhancements for the proposed system should primarily focus on optimizing and designing more efficient heat exchangers.

Acknowledgment

This work was supported by the State Key Laboratory of Low-Carbon Smart Coal-Fired Power Generation and Ultra-clean Emission and National Natural Science Foundation of China (Grant No. 52022020).

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Li
,
T. X.
,
Wang
,
R. Z.
,
Yan
,
T.
, and
Ishugah
,
T. F.
,
2014
, “
Integrated Energy Storage and Energy Upgrade, Combined Cooling and Heating Supply, and Waste Heat Recovery With Solid–Gas Thermochemical Sorption Heat Transformer
,”
Int. J. Heat Mass Transfer
,
76
, pp.
237
246
.
2.
Sun
,
W.
,
Yue
,
X.
, and
Wang
,
Y.
,
2017
, “
Exergy Efficiency Analysis of ORC (Organic Rankine Cycle) and ORC-Based Combined Cycles Driven by Low-Temperature Waste Heat
,”
Energy Convers. Manage.
,
135
, pp.
63
73
.
3.
Xia
,
L.
,
Liu
,
R.
,
Zeng
,
Y.
,
Zhou
,
P.
,
Liu
,
J.
,
Cao
,
X.
, and
Xiang
,
S.
,
2019
, “
A Review of Low-Temperature Heat Recovery Technologies for Industry Processes
,”
Chin. J. Chem. Eng.
,
27
(
10
), pp.
2227
2237
.
4.
Han
,
B. C.
,
Cheng
,
W. L.
,
Li
,
Y. Y.
, and
Nian
,
Y. L.
,
2017
, “
Thermodynamic Analysis of Heat Driven Combined Cooling Heating and Power System (CCHP) With Energy Storage for Long Distance Transmission
,”
Energy Convers. Manage.
,
154
, pp.
102
117
.
5.
Mohammadi
,
A.
,
Kasaeian
,
A.
,
Pourfayaz
,
F.
, and
Ahmadi
,
M. H.
,
2017
, “
Thermodynamic Analysis of a Combined Gas Turbine, ORC Cycle and Absorption Refrigeration for a CCHP System
,”
Appl. Therm. Eng.
,
111
, pp.
397
406
.
6.
Zhou
,
T.
,
Liu
,
J.
,
Ren
,
J.
, and
Yang
,
S.
,
2022
, “
Conceptual Design, Modelling and Optimization of an Integrated System by Combining Organic Rankine Cycle and Absorption Refrigeration Cycle for Efficient Energy Recovery
,”
J. Taiwan Inst. Chem. Eng.
,
133
, p.
104276
.
7.
Emadi
,
M. A.
,
Chitgar
,
N.
,
Oyewunmi
,
O. A.
, and
Markides
,
C. N.
,
2020
, “
Working-Fluid Selection and Thermoeconomic Optimisation of a Combined Cycle Cogeneration Dual-Loop Organic Rankine Cycle (ORC) System for Solid Oxide Fuel Cell (SOFC) Waste-Heat Recovery
,”
Appl. Energy
,
261
, p.
114384
.
8.
Wang
,
C.
,
Li
,
Q.
,
Wang
,
C.
,
Zhang
,
Y.
, and
Zhuge
,
W.
,
2021
, “
Thermodynamic Analysis of a Hydrogen Fuel Cell Waste Heat Recovery System Based on a Zeotropic Organic Rankine Cycle
,”
Energy
,
232
, p.
121038
.
9.
Wang
,
Z.
,
Xia
,
X.
,
Pan
,
H.
,
Zuo
,
Q.
,
Zhou
,
N.
, and
Xie
,
B.
,
2021
, “
Fluid Selection and Advanced Exergy Analysis of Dual-Loop ORC Using Zeotropic Mixture
,”
Appl. Therm. Eng.
,
185
, p.
116423
.
10.
Zhi
,
L. H.
,
Hu
,
P.
,
Chen
,
L. X.
, and
Zhao
,
G.
,
2020
, “
Performance Analysis and Optimization of Engine Waste Heat Recovery With an Improved Transcritical-Subcritical Parallel Organic Rankine Cycle Based on Zeotropic Mixtures
,”
Appl. Therm. Eng.
,
181
, p.
115991
.
11.
Zhou
,
T.
,
Liu
,
J.
,
Ren
,
J.
, and
Yang
,
S.
,
2022
, “
Comprehensive Assessment of a Coupled LiBr/H2O Absorption Refrigeration/ORC System for Low-Grade Residual Heat Recovery Based on Advanced Exergy and Exergoeconomic Analysis
,”
ACS Sustain. Chem. Eng.
,
10
(
18
), pp.
5825
5837
.
12.
Wang
,
S.
,
Liu
,
C.
,
Zhang
,
S.
,
Li
,
Q.
, and
Huo
,
E.
,
2022
, “
Multi-objective Optimization and Fluid Selection of Organic Rankine Cycle (ORC) System Based on Economic-Environmental-Sustainable Analysis
,”
Energy Convers. Manage.
,
254
, p.
115238
.