Heat Exchanger Network Synthesis Integrated with Compression–Absorption Cascade Refrigeration System

: Compression–absorption cascade refrigeration system (CACRS) is the extension of absorption refrigeration system, which can be utilized to recover excess heat of heat exchanger networks (HENs) and compensate refrigeration demand. In this work, a stage-wise superstructure is presented to integrate the generation and evaporation processes of CACRS within HEN, where the generator is driven by hot process streams, and the evaporation processes provide cooling energy to HEN. Considering that the operating condition of CACRS has signiﬁcant e ﬀ ect on the coe ﬃ cient of performance (COP) of CACRS and so do the structure of HEN, CACRS and HEN are considered as a whole system in this study, where the operating condition and performance of CACRS and the structure of HEN are optimized simultaneously. The quantitative relationship between COP and operating variables of CACRS is determined by process simulation and data ﬁtting. To accomplish the optimal design purpose, a mixed integer non-linear programming (MINLP) model is formulated according to the proposed superstructure, with the objective of minimizing total annual cost (TAC). At last, two case studies are presented to demonstrate that desired HEN can be achieved by applying the proposed method, and the results show that the integrated HEN-CACRS system is capable to utilize energy reasonably and reduce the total annualized cost by 38.6% and 37.9% respectively since it could recover waste heat from hot process stream to produce the cooling energy required by the system.


Introduction
Exchanging heat between cold and hot process streams through heat exchangers is an effective form of waste heat recovery. A great many of researches have been undertaken, such as the optimization design of heat exchanger, and the synthesis of heat exchanger networks. For the optimization design of heat exchanger, the nanomaterials were mixed in the pure fluid since the nanoparticles have higher thermal conductivity and can give greater characteristics of carrier fluid [1]. In addition, the turbulent flow of nanofluid in a pipe has been modeled in order to reach better design in view of second law [2]. Also, the researches on heat exchanger network (HEN) synthesis have received considerable development, and still bolstered by new ideas. Mathematical programming approach based on superstructure is one of the most widely used methods for HEN synthesis. In recent years, the original stage-wise superstructure proposed by Yee and Grossmann [3] has been improved and perfected. Such as the assumption of isothermal mixing was shattered [4], and the superstructure was enhanced to contain more possibilities [5], as well as to consider the parallel connection of utilities in stream split branches [6] and series connection of utilities in substages [7].
The challenge of HEN synthesis with CACRS integrated is how to describe the coupling relationship between HEN and CACRS quantitatively, for which the COPs affected by operating condition are the crucial parameters to quantify the relationship. This paper handles the problem by means of discretizing the evaporation temperatures of two subsystems. At each evaporation temperature, the quantitative relationships between COPar 1 and heat source temperature, the quantitative relationships between COPec 2 and condensation temperature, are both obtained by simulating and data fitting. Therefore, the CACRS can be integrated into HEN in the form of a more detailed mathematic formulation. With this method, the optimal integrated system can be obtained with better performance and higher energy efficiency, large amount of excess heat can be recovered to minimize the economic cost.
This paper is organized as follows. The problem to be addressed in this paper is stated as Section 2. The simulation and model formulation are conducted in Section 3, wherein the two subsystems of CACRS are simulated, and the quantitative relationship between COPs and the operating variables of CACRS is determined with the simulation results, and the synthesis method for HEN integrated with CACRS is presented, which includes the superstructure composed of plentiful network alternatives and the detailed mathematical description for optimization purpose. The proposed method is demonstrated through two case studies in Section 4, and the discussion is presented in Section 5. Finally, conclusions are drawn in Section 6. Figure 1 shows the schematic representation of the CACRS synthesized with heat source. The system is composed of an absorption subsystem and a single stage electrical compression subsystem. These two subsystems share a heat exchanger which operates as the evaporator of absorption section meanwhile the condenser of compression section simultaneously. The generator of absorption subsystem is heated by heat sources either waste heat of hot process streams or external steam, and the evaporators of two subsystems provide process streams with cooling energy. The absorber and condenser are cooled down with cooling water. Here it should be noted that, when the required refrigerating temperature, lower than room temperature while no less than 0 • C, can be achieved by the evaporator of absorption subsystem, electrical compression subsystem is not needed any more. When the refrigerating temperature is too low to be achieved by ARC, below 0 • C, the demand of refrigeration is satisfied by compression subsystem's evaporation process of CACRS. subsystem and 2 COPec for the compression subsystem are used as indicators to evaluate the performance) of CACRS. Therefore, in this paper, CACRS and HEN are considered as a whole system to recover heat for refrigeration, where the operating parameters and COPs of CACRS and the structure of HEN are optimized simultaneously.

Problem Statement
The challenge of HEN synthesis with CACRS integrated is how to describe the coupling relationship between HEN and CACRS quantitatively, for which the COPs affected by operating condition are the crucial parameters to quantify the relationship. This paper handles the problem by means of discretizing the evaporation temperatures of two subsystems. At each evaporation temperature, the quantitative relationships between 1 COPar and heat source temperature, the quantitative relationships between 2 COPec and condensation temperature, are both obtained by simulating and data fitting. Therefore, the CACRS can be integrated into HEN in the form of a more detailed mathematic formulation. With this method, the optimal integrated system can be obtained with better performance and higher energy efficiency, large amount of excess heat can be recovered to minimize the economic cost.
This paper is organized as follows. The problem to be addressed in this paper is stated as Section 2. The simulation and model formulation are conducted in Section 3, wherein the two subsystems of CACRS are simulated, and the quantitative relationship between COPs and the operating variables of CACRS is determined with the simulation results, and the synthesis method for HEN integrated with CACRS is presented, which includes the superstructure composed of plentiful network alternatives and the detailed mathematical description for optimization purpose. The proposed method is demonstrated through two case studies in Section 4, and the discussion is presented in Section 5. Finally, conclusions are drawn in Section 6. Figure 1 shows the schematic representation of the CACRS synthesized with heat source. The system is composed of an absorption subsystem and a single stage electrical compression subsystem. These two subsystems share a heat exchanger which operates as the evaporator of absorption section meanwhile the condenser of compression section simultaneously. The generator of absorption subsystem is heated by heat sources either waste heat of hot process streams or external steam, and the evaporators of two subsystems provide process streams with cooling energy. The absorber and condenser are cooled down with cooling water. Here it should be noted that, when the required refrigerating temperature, lower than room temperature while no less than 0 °C, can be achieved by the evaporator of absorption subsystem, electrical compression subsystem is not needed any more. When the refrigerating temperature is too low to be achieved by ARC, below 0 °C, the demand of refrigeration is satisfied by compression subsystem's evaporation process of CACRS.  The COPs are crucial parameters to establish the relationship between the cold energy requirements in the evaporator and heat required by the generator (or the electricity consumed by the compressor). The mathematical expressions of COPar 1 for ARC and COPec 2 for ECR are in Equations (1) and (2), respectively. Whereas, COPar 1 and COPec 2 are reflected by the operating condition of ARC and ECR

Problem Statement
The synthesis problem to be addressed in this paper can be stated as follows.
A set of hot process streams (some are below room temperature or below 0 • C), and a set of cold process streams are given with their heat capacity flow rates, supply and target temperatures. Parameters of utilities, such as their inlet temperatures, outlet temperatures and unit costs are also given. Other parameters given include the cost coefficients of heat transfer units, the minimum approach temperatures and so on. Then, the problem consists in determining the optimal design of HEN coupled with CACRS including the optimal design of HEN structure and the operating condition of CACRS, as well as the trade-off between the efficiency and amount of heat recovery, with the aim of minimizing the total annual cost (TAC). Thus, to solve the problem, the relationship between COPs and operating condition is simulated in this paper and based on which, the mathematical model of HEN integrated with CACRS is established for optimal design.

Simulation of Compression-Absorption Cascade Refrigeration System
The CACRS is modelled and simulated with the process modelling software Aspen Plus. Two steady state simulation models for the two subsystems have been structured and implemented based on sequential modular approach in the computer program. The simulation used LiBr-H 2 O fluid pairs within absorption subsystem and R134a within compression subsystem. The ELECNRTL and STEAMNBS property methods in Aspen Plus are chosen for LiBr-H 2 O solution and pure water. The RK-SOAVE property method is used for R134a. It is assumed that the pressure and heat losses are neglectable, working fluid at the exits of generator, absorber, evaporator, cascaded evaporator-condenser and condenser are all in saturation state. The model is validated by comparing with the data presented by Cimsit et al. [21].
In this study, the condensation temperature is set as 45 • C and the concentration of solution is set as 0.58 for absorption section, referring to Wang's work [29]. A series of evaporation temperatures alternatives (5 • C, 7 • C, 9 • C, 11 • C, 13 • C, 15 • C) are given for choosing, and the heat source temperature will be varied in certain range (from 105 • C to 175 • C). For the compression section, the optional evaporation temperatures are optimized in a series of values (−25 • C, −21 • C, −17 • C, −13 • C, −9 • C, −5 • C) according to refrigeration demands, and the corresponding condensation temperature are 10 • C, 12 • C, 14 • C, 16 • C, 18 • C, and 20 • C when the temperature difference of cascaded evaporator-condenser is supposed to be 5 • C. Figure 2a,b present the variation of COP with heat source temperature and condensation temperature at certain evaporation temperature for absorption and compression section separately. Obviously, COPar 1 is improved with the increase of heat source temperature (t) and evaporation temperature. Higher evaporation temperature (te) and lower condensation temperature (tc) contribute to higher COPec 2 . The corresponding functional relationship between COPar 1 and heat source temperature can be obtained in form of fitting curve, which is performed within MATLAB. The values calculated according to fitting function are compared with the simulated value to examine the accuracy, and the maximum error is less than 2%. In the same way, the matrix COPec 2 determined by evaporation and condensation temperatures can be also established. the simulated value to examine the accuracy, and the maximum error is less than 2%. In the same way, the matrix 2 COPec determined by evaporation and condensation temperatures can be also established.

Superstructure Presentation
A stage-wise superstructure is proposed to present the configuration of HEN incorporated with the generation and evaporation processes of CACRS in Figure 3. The cooperation of HEN and CACRS can be performed in two parts: (i) Hot process streams provide heat to drive the generators of absorption subsystem at the inner stage, (ii) low-temperature cooling energy produced in evaporators of absorption and compression subsystems are used to cool hot process streams to subambient temperature or below 0 °C at the stream ends after the use of cooling water.
As indicated, the generation (CACRSGE) and evaporation processes of absorption (ARCEV) and compression (ECREV) subsystems are implemented cooperatively with the heat exchange of mentioned process streams. Also, LPS can be used as heat source to motivate CACRS for producing cooling energy when hot process streams can't provide enough heat. The corresponding generation and evaporation process are LPSCACRSGE, LPSARCEV, and LPSECREV. Each hot process stream can be cooled down by cold process stream, cooling water and/or the generation process (CACRSGE) and evaporation processes (ARCEV, ECREV, LPSARCEV, and LPSECREV) of CACRS, and each cold process stream can be heated by hot process stream, low-pressure steam (LPS) and/or high-pressure steam (HPS). In addition, a two-stage structure of heat exchange between hot process stream and the evaporation processes of CACRS are performed at the stream end in order to include more matching possibilities.

Superstructure Presentation
A stage-wise superstructure is proposed to present the configuration of HEN incorporated with the generation and evaporation processes of CACRS in Figure 3. The cooperation of HEN and CACRS can be performed in two parts: (i) Hot process streams provide heat to drive the generators of absorption subsystem at the inner stage, (ii) low-temperature cooling energy produced in evaporators of absorption and compression subsystems are used to cool hot process streams to sub-ambient temperature or below 0 • C at the stream ends after the use of cooling water.

Model Formulation
A MINLP model consisting of HEN and CACRS is formulated according to the superstructure proposed in this paper. The following sets are defined before presenting the model formulation: I, J, ST, ST' represent the sets for hot process streams, cold process streams, stage numbers of heat exchange inside the superstructure, and stage numbers of heat exchange at end of the superstructure, S and L are the sets for temperature grades of evaporation processes for CACRS.

•
Total energy balances for process streams As indicated, the generation (CACRSGE) and evaporation processes of absorption (ARCEV) and compression (ECREV) subsystems are implemented cooperatively with the heat exchange of mentioned process streams. Also, low-pressure steam (LPS) can be used as heat source to motivate CACRS for producing cooling energy when hot process streams can't provide enough heat. The corresponding generation and evaporation process are LPSCACRSGE, LPSARCEV, and LPSECREV. Each hot process stream can be cooled down by cold process stream, cooling water and/or the generation process (CACRSGE) and evaporation processes (ARCEV, ECREV, LPSARCEV, and LPSECREV) of CACRS, and each cold process stream can be heated by hot process stream, LPS and/or high-pressure steam (HPS).
In addition, a two-stage structure of heat exchange between hot process stream and the evaporation processes of CACRS are performed at the stream end in order to include more matching possibilities.

Model Formulation
A MINLP model consisting of HEN and CACRS is formulated according to the superstructure proposed in this paper. The following sets are defined before presenting the model formulation: I, J, ST, ST' represent the sets for hot process streams, cold process streams, stage numbers of heat exchange inside the superstructure, and stage numbers of heat exchange at end of the superstructure, S and L are the sets for temperature grades of evaporation processes for CACRS.

Model for Heat Exchanger Network Synthesis
• Total energy balances for process streams Equations (3) and (4) It should be note that all equal to t i,k , t in j,i,k and t LPSin j,k equal to t j,k+1 , on account that the inlet temperature of stage k is the outlet temperature after mixing in previous stage.
• Non-isothermal mixing for each inner stage For process streams, the mass balances and energy balances for internal stages of the superstructure are shown as Equations (10)- (13). With these equations, the outlet temperature of stage k after mixing can be calculated.
• Energy balances for matches of evaporation processes at the stream ends The heat removed by the evaporation processes of CACRS in any stage k' which is at the ends of hot streams are given by Equations (14)- (17). , but the evaporation process belongs to the compression section of CACRS.
• Energy balance of non-isothermal mixing for evaporation processes at the stream ends • Energy balances for utilities at the stream ends T out j − t HPSin At the end of the superstructure, the heat removed by cooling water is determined by Equation (20), and the heat removed by the evaporation processes of CACRS are given by Equation (21). With Equations (22) and (23), the heating requirement of each cold process stream at the stream end are satisfied by LPS and HPS are determined.
• Feasibility constraints of temperature Equation (24) denotes the temperature of hot process stream i will decrease after exchanging heat with M, M is the set including cold process stream, generator of CACRS, cooling water, and evaporators of CACRS. Similarly, Equation (25) denotes the temperature of cold process stream j will increase after exchanging heat with N which contains hot process stream, LPS and HPS. Equations (26) and (27) are used to ensure the monotonic decrease of temperature from the left side to the right side of the superstructure. The constraints on binary variable of CACRS are given by Equations (28)-(32). Equation (28) means that for the generation process of CACRS, there is at most one corresponding evaporation temperature of certain grade. Equations (29)-(32) indicate that for evaporation processes of CACRS, only one grade of evaporation temperature is considered so as to simplify the model and shorten computing time.

Model for CACRS
In the CACRS, condensation temperature and absorption temperature can be regarded as constants since the absorber and condenser are cooled down with cooling water which is used as cold utility in HEN with given inlet and outlet temperatures.
• The determination of matrix COP In the case that condensation and absorption temperatures of absorption section are defined, COPar 1 is only related to the inlet temperature of generator (heat source temperature) t CACRSGEin i,k and the evaporation temperature T ARCEV for ARC, COPec 2 can be determined by the condenser temperature (equal to T ARCEV + ∆T min2 ) and evaporation temperature T ECREV for ECR. Evaporation temperature of two subsystems are graded by discretizing into several values. Once the grades of these two evaporation temperatures are determined, copar 1 (s) can be expressed as the function of t CACRSGEin  (35) give the energy balances of the CACRS, which means that the evaporator of absorption subsystem get heat from hot process stream and the condenser of compression subsystem.

Objective Function
The heat exchanging area can be calculated with Equation (36), herein, h M and h N represent the film heat transfer coefficients, dt 1 and dt 2 are the temperature differences before and after heat exchanging, the mean temperature difference is approximated according to Chen [30]. Objective of the optimization is to minimize the TAC, considering both the expenditures of capital and operation for HEN and CACRS. The corresponding cost representations are given in Equations (38)-(40). HEN cost is the summation of capital cost for heat exchangers and operating cost for utility consumption. CACRS cost is deduced from the investment of absorption subsystem (containing generator, evaporator, condenser and absorber) and compression subsystem (cascaded evaporation-condenser, compressor and evaporator are included), as well as the cooling water consumption in the absorber and condenser which belong to absorption subsystem and the electricity consumption of compressor in compression subsystem.

Case Study
In this section, the proposed model is demonstrated through two case studies. To make comparison, a base case is set to cool down process streams by employing electrical compression refrigeration (R134a), which produces cooling energy by consuming electricity rather than process heat. Analogously, for the HEN integrated with ECR, each hot process stream is allowed to be cooled down by multiple parallel evaporation processes belonging to different grades after the using of cooling water. In the absorption subsystem of CACRS, LiBr/H 2 O is used as the working fluid. For single LiBr/H 2 O absorption refrigeration, the temperature of heat source is not allowed to be lower than 80 • C or higher than 180 • C. In this paper, the COPar 1 lower than 0.5 is not considered, so the designated range of heat source temperature is 105 • C to 175 • C. The minimum approach temperature for process stream matching and CACRS matching are stipulated to be ∆T min1 = 10 • C, ∆T min2 = 5 • C respectively. The inlet and outlet temperatures of cooling water are 30 • C and 40 • C, temperatures of LPS and HPS are 158 • C and 255 • C respectively. Similar to CACRS, the condenser of ECR is also cooled down by cooling water, and accordingly, the COP of ECR (COPec) can be determined by the selectable evaporation temperatures T ECREV (denoted by s and l) which are identical to CACRS.
The data of the involved process streams for the two cases are shown in Tables 1 and 2. In case 2, one more hot process stream is added into the HEN, which means more heat available, meanwhile more demand for cooling. In order to optimize the operating condition of CACRS, six evaporation temperature alternatives for absorption subsystem and compression subsystem are given and graded respectively according to the requirement of target temperature. The grading about evaporation temperatures for absorption subsystem is shown in Table 3. Similarly, the evaporation temperature for compression subsystem is also graded into six levels in the range of −25 • C to −5 • C, with l 1 to l 6 denoting −25 • C, −21 • C, −17 • C, −13 • C, −9 • C, and −5 • C respectively. Based on the simulation and data fitting procedures stated in Section 3, the functional relation of COPar 1 with t CACRSGEin i,k is obtained and presented as Equation (40), and the correlation coefficients are listed in Table 3. Corresponding cop ec 2 s,l in matrix COPec 2 are given in Equation (41) according to the general presentation of Equation (33). When only ECR is used to satisfy the refrigeration demand, the cop ec s,l in matrix COPec is given in Equation (42), wherein, the cop ec s used to produce cooling energy below room temperature is given as Equation (43), and the cop ec l used to produce cooling energy lower than 0 • C is given as Equation (44).   The formulated MINLP mathematical models for mentioned HEN-CACRS and HEN-ECR designs are coded and solved in GAMS with PC specification: CPU 3.30 GHz, 8 GB RAM. For each model, BARON is used as solver. To be clear, note that the global optimum cannot be guaranteed for MINLP problems. Yet, good initial values and boundaries as well as circular iteration are employed to facilitate the search process and improve the quality of locally-optimal solutions.
In HEN-CACRS result, the evaporation temperature of grade 's' produced in ARC subsystem can cool hot process stream to below room temperature, and the evaporation temperature of grade 'l', produced in ECR subsystem, can be used for cooling hot process stream to below 0 • C. For the integrated CACRS, (s, l) denotes that there is cooling energy produced at evaporation temperature of grade 'l' by the compression section, whose temperature of cascaded evaporator-condenser belongs to grade 's', while for HEN-ECR result, just the evaporation temperature grade of 's' or 'l' needs to be determined. Optimal structures of the integrated HENs with CACRS and ECR for two cases are presented in Figure 4a,b and Figure 5a,b, respectively.   Figure 4a shows the optimal structure of HEN integrated with CACRS for case 1, as indicated in which one set of CACRS and one set of ARC are implemented to absorb heat from process stream for refrigeration purpose. The CACRSGE (5 °C) is motivated by H1 and H2, and the cooling energy, 5 °C (s1) produced by the evaporation process of absorption section, and −13 °C (l4) produced in evaporator of compression section are used to cool down H1 (ARCEV(5 °C), and ECREV (5 °C, −13 °C)) successively. In addition to exchanging heat with cold process stream and the generator of CACRS, H2 heats the generator of the ARC (CACRSGE (9 °C)), then to be refrigerated by the evaporation  Figure 4a shows the optimal structure of HEN integrated with CACRS for case 1, as indicated in which one set of CACRS and one set of ARC are implemented to absorb heat from process stream for refrigeration purpose. The CACRSGE (5 • C) is motivated by H1 and H2, and the cooling energy, 5 • C (s 1 ) produced by the evaporation process of absorption section, and −13 • C (l 4 ) produced in evaporator of compression section are used to cool down H1 (ARCEV(5 • C), and ECREV (5 • C, −13 • C)) successively. In addition to exchanging heat with cold process stream and the generator of CACRS, H2 heats the generator of the ARC (CACRSGE (9 • C)), then to be refrigerated by the evaporation process of ARC (ARCEV(9 • C), (s 3 )) after the use of cooling water. In HEN-ECR result, as shown in Figure 4b, the mentioned refrigeration demands are all satisfied by ECRs. The evaporating temperature for refrigerating H1 and H2 are −13 • C (l 4 ).
For case 2, the optimal integrated structures of the two concerned design scenarios are shown in Figure 5a,b. For the scenario HEN integrated with CACRS, two ARC sets and two CACRS sets are involved to provide cooling energy to HEN. Since the heat of process streams is not sufficient, LPS is used to compensate the heat demand of CACRS. One ARC set is run by LPS (LPSCACRSGE (5 °C)), and the evaporator absorbs heat from H1 (LPSARCEV (5 °C), (s1)). The other ARC set gets heat from H3 in generator (CACRSGE (5 °C)) and the cooling energy produced in the evaporator is used to cool H2 and H3 (ARCEV (5 °C), (s1)). One CACRS is motivated by LPS (LPSCACRSGE (9 °C)), and the cooling energy produced in the evaporator of compression section is used to cool down H1 (LPSECREV (9 For the scenario HEN integrated with CACRS, two ARC sets and two CACRS sets are involved to provide cooling energy to HEN. Since the heat of process streams is not sufficient, LPS is used to compensate the heat demand of CACRS. One ARC set is run by LPS (LPSCACRSGE (5 • C)), and the evaporator absorbs heat from H1 (LPSARCEV (5 • C), (s 1 )). The other ARC set gets heat from H3 in generator (CACRSGE (5 • C)) and the cooling energy produced in the evaporator is used to cool H2 and H3 (ARCEV (5 • C), (s 1 )). One CACRS is motivated by LPS (LPSCACRSGE (9 • C)), and the cooling energy produced in the evaporator of compression section is used to cool down H1 (LPSECREV (9 • C, −25 • C), (s 3 , l 1 )) at the stream end. The other CACRS gets heat from H2 (CACRSGE (7 • C)) and the cooling energy produced in the evaporator of absorption and compression subsystems is used for cooling H2 (ARCEV (7 • C), (s 2 )) and H3 (ECREV (7 • C, −25 • C), (s 2 , l 1 )) to their target temperatures separately.

Discussion
From obtained results of the two cases, it can be observed that, in HEN-CACRS and HEN-ECR scenarios, the process streams with cooling requirements are always cooled to the lowest temperature which can be achieved by cooling water first, and then are refrigerated by the evaporation processes of refrigeration system subsequently. The cause of this state is that cooling energy is produced at the cost of consuming much heat and electricity which will increase the operating cost. What is more, by comparing the optimal structure of HEN-CACRS with HEN-ECR, it can be observed that HEN-ECR prefers equal or higher evaporating temperature grades than that of CACRS at the last stage of evaporation processes. The reason for this tendency is that, higher evaporating temperature is beneficial for higher COP, which will decrease the consumption of expensive electricity.
Higher COP is pursued by both ECR and CACRS, which means consuming less electricity and heat. But when CACRS is integrated into HEN, high COP cannot be reached all the time due to the limit of cost target and characteristics of process streams (inlet temperature, outlet temperature, heat load and so on), there are trade-offs between the quantity and efficiency of heat recovery and between the cost of hot utility and electricity consumption.
For the absorption section, the lower the evaporation temperatures, the higher the hot source temperature (inlet temperature of generator) will be needed for higher COPar 1 . Such as the evaporation temperature used to cool H2 in case1 is 9 • C (s 3 ), the corresponding heat source temperature is 139.33 • C, the evaporating temperature which is used to cool H1 is 5 • C (s 1 ), and the matching heat temperature is 145 • C. But there is also 5 • C (s 1 ) cooling energy produced, which is motivated the heat source of 135 • C, for recovering much heat. Likewise, for the compression subsystem of CACRS, in case 2, there are two matches of evaporation temperatures for ARC section and ECR section, (7 • C, −25 • C) and (9 • C, −25 • C), obviously, the match of (7 • C, −25 • C) has higher COPec 2 , but the match of (9 • C, −25 • C) is motivated by LPS, higher evaporation temperature of ARC section will lead to higher COPar 1 and lower COPec 2 at the same time, so there is a trade-off between LPS and electricity consumption costs.
For the two cases, the optimal integrated HEN-ECR system are achieved at the TAC of 2,595,319 $·y −1 and 4,780,274 $·y −1 , and the TAC of optimal integrated HEN-CACRS are 1,592,527 $·y −1 and 2,970,188 $·y −1 , which are 38.6% and 37.9% lower than HEN-ECR respectively. The concrete analysis is as follows. The TAC consists of the capital cost of HEN, capital cost of CACRS and operating cost of HEN-CACRS. The operating cost of HEN-CACRS is composed of electricity cost, hot utility cost (LPS and HPS) and cold utility cost (cooling water). The cost of each part can be seen in Figure 6. $·y −1 and 4,780,274 $·y −1 , and the TAC of optimal integrated HEN-CACRS are 1,592,527 $·y −1 and 2,970,188 $·y −1 , which are 38.6% and 37.9% lower than HEN-ECR respectively. The concrete analysis is as follows. The TAC consists of the capital cost of HEN, capital cost of CACRS and operating cost of HEN-CACRS. The operating cost of HEN-CACRS is composed of electricity cost, hot utility cost (LPS and HPS) and cold utility cost (cooling water). The cost of each part can be seen in Figure 6. Apparently, the capital cost of refrigeration forms a considerable part of TAC, and the cost of CACRS is close to that of ECR for the two cases studied. By comparing the cost difference of HEN-CACRS with HEN-ECR, it can be found that there is little difference in the capital cost of HEN, the hot and cold utilities costs of HEN-CACRS are higher than that of HEN-ECR, but the electricity cost of HEN-CACRS is much less than that of HEN-ECR. There are the detail explanations. In HEN-CACRS scenario, the integration of CACRS helps to promote energy efficiency of hot process streams, some waste heat is used to motivate the generator of CACRS, and the energy produced by the Apparently, the capital cost of refrigeration forms a considerable part of TAC, and the cost of CACRS is close to that of ECR for the two cases studied. By comparing the cost difference of HEN-CACRS with HEN-ECR, it can be found that there is little difference in the capital cost of HEN, the hot and cold utilities costs of HEN-CACRS are higher than that of HEN-ECR, but the electricity cost of HEN-CACRS is much less than that of HEN-ECR. There are the detail explanations. In HEN-CACRS scenario, the integration of CACRS helps to promote energy efficiency of hot process streams, some waste heat is used to motivate the generator of CACRS, and the energy produced by the evaporators of CACRS can be used to partly meet the refrigeration demand, which lead to sharp decrease of electricity cost. Even though the utilization of waste heat can decrease the consumption of cooling water for cooling hot process streams, the cooling water cost of HEN-CACRS is still a little higher than that of HEN-ECR, because the cooling water consumed by the refrigeration system is also considered. In addition, the integration of CACRS into HEN may lessen the heat load matching with cold process stream or increase the consumption of LPS as heat source to motivate the refrigeration system, and result in the increase of hot utility cost. But with comprehensive consideration of all costs, it can be concluded that the TAC of HEN-CACRS is still lower than that of HEN-ECR.

Conclusions
This paper has presented a mathematical programming model to consider the integration of HEN with CACRS. The coupling relationship of CACRS and HEN are quantitatively described through process simulating and data fitting. Then a detailed mathematical description for CACRS and HEN is formulated. Considering the mutual effect of HEN and CACRS, the structure of the HEN and the performance of CACRS are optimized simultaneously within the design. Compared with the results of base case where only electrical compression refrigeration is used for refrigeration purpose, it can be found that the integrated HEN-CACRS system is capable to utilize energy reasonably and reduce the total annualized cost by 38.6% and 37.9% respectively since it could recover waste heat from hot process stream to produce the cooling energy required by the system. On base of this study, future work will be concentrated on the optimal design of the HEN integrated with CACRS and ORC, and more internal parameters will be considered. H Y annual working hour, h C cw unit cost of cooling water, $/(kW·y) C lps unit cost of low-pressure steam, $/(kW·y) C hps unit cost of high-pressure steam, $/(kW·y) C elec unit cost of electricity, $/(kW·y) C comp capital cost of compressor, $·y −1 Superscripts CACRS compression-absorption cascade refrigeration system LPS/LPS' low-pressure steam LPSCACRSGE the evaporation process belongs to absorption subsystem of CACRS motivated by low-pressure steam CACRSGE the generation process of CACRS motivated by hot process stream ARCEV the evaporation process belongs to absorption subsystem of CACRS motivated by hot process stream ECREV the evaporation process belongs to compression subsystem of CACRS motivated by hot process stream LPSARCEV the evaporation process belongs to absorption subsystem of CACRS motivated by low-pressure steam LPSECREV the evaporation process belongs to compression subsystem of CACRS motivated by low pressure steam CW/CW' cooling water HPS high-pressure steam Subscripts i hot process stream j cold process stream k index for inner stage of HEN superstructure (1, . . . , NOK) NOK total number of inner stages k' index for evaporation stage at hot stream end s evaporation temperature grades of the CACRS's absorption subsystem l evaporation temperature grades of the CACRS's compression subsystem Variables t in inlet temperature, • C t out outlet temperature, • C fcp heat capacity flowrate for a stream branch, kW· • C −1 q heat exchange load, kW A heat exchange area, m 2 Wele the consumption of electricity, $ COPar 1 coefficient of performance for the absorption section of CACRS COPec 2 coefficient of performance for the compression section of CACRS h M /h N film heat transfer coefficient of the stream needs to be heated or cooled, kW·( • C·m 2 ) −1 dt 1 /dt 2 temperature difference before or after heat exchanging, • C