Temperature Disturbance Management in a Heat Exchanger Network for Maximum Energy Recovery Considering Economic Analysis

: The design of heat exchanger networks (HEN) in the process industry has largely focused on minimisation of operating and capital costs using techniques such as pinch analysis or mathematical modelling. Aspects of operability and flexibility, including issues of disturbances affecting downstream processes during the operation of highly integrated HEN, still need development. This work presents a methodology to manage temperature disturbances in a HEN design to achieve maximum heat recovery, considering the impact of supply temperature fluctuations on utility consumption, heat exchanger sizing, bypass placement and economic performance. Key observations have been made and new heuristics are proposed to guide heat exchanger sizing to consider disturbances and bypass placement for cases above and below the HEN pinch point. Application of the methodology on two case studies shows that the impact of supply temperature fluctuations on downstream heat exchangers can be reduced through instant propagation of the disturbances to heaters or coolers. Where possible, the disturbances have been capitalised upon for additional heat recovery using the pinch analysis plus-minus principle as a guide. Results of the case study show that the HEN with maximum HE area yields economic savings of up to 15% per year relative to the HEN with a nominal HE area.


Introduction
Heat integration has been a well-established energy saving technique for the chemical process industry since the global energy crisis in the 1970s [1]. There has also been extensive development of pinch analysis methodologies for industrial heat exchanger network (HEN) design focusing on minimisation of operating and capital costs [2]. The developed methodologies typically assume that process parameters such as supply/target temperatures and stream flowrates are fixed [3]. In practice, these process parameters may fluctuate due to plant start-ups/shut-downs, changes in feed or product demand as well as quality, changes in environmental conditions and other operational disturbances. The impact of these parameter changes influences energy related-decision making as a step for efficient energy management in the industry [4,5]. The extent to which a HEN is able to cope with disturbances is known as flexibility [3]. Previous works have included flexibility, safety, controllability and operability in the pre-design stage.
Marselle et al. [6] pioneered the study of operability considerations for HENs. They proposed a manual combination of several resilient designs with maximum energy efficiency under different worst-case scenarios. Hafizan et al. [7] proposed the controllability of HEN under uncertainty by having a large heat exchanger area. This method is able to recover maximum heat with minimum utility. The paper was also partly presented at the 13th SDEWES conference and proceedings. Linnhoff and Kotjabasakis [8] introduced the concept of downstream paths to identify disturbance propagation paths through HEN. Supply temperature and feed flowrate disturbances occurring at the feed stream of heat exchangers can affect the target temperatures of processes located downstream of the heat exchanger path. HEN modification, therefore, may need to be performed to reject the process disturbances.
Escobar et al. [9] introduced a computational framework for the synthesis of controllable and flexible HENs over a specified range of inlet temperatures and flow rate variations using a decentralised control system. This framework comprised two stages-the HEN design and the operability analysis and adjustment of control variables during operation in the presence of uncertain parameters. Hafizan et al. [10] developed a pinch analysis-based methodology, which considered both safety and operability aspects in HEN synthesis. The concept of downstream paths suggested by Linnhoff and Kotjabasakis [8] was used for the assessment of flexibility and structural controllability. Čuček and Kravanja [11] proposed a novel three-step methodology for HEN retrofit with fixed and flexible designs for large-scale total sites (TS). This method modified the HEN by forming profitable heat exchanger matches with improved utility consumption as well as by proposing intermediate utility production.
Most recent works have developed a multi-period formulation for the synthesis of flexible HENs. Isafiade and Short [12] proposed a three-step approach for improving the degree of flexibility of a multi-period HEN with unequal periods. However, this method is unable to cater for all the possible uncertainties of process parameters. Miranda et al. [13] recently proposed three sequential steps represented by linear programming (LP), multi-linear programming (MILP) and non-linear programming (NLP) for the synthesis of multi-period HEN. The same heat transfer can be operated under different operating conditions in the multi-period [13]. Kang and Liu [14] developed a threestep method for designing a flexible multi-period HEN when the disturbances of operating parameters occurred in sub-periods. The flexibility of nominal multi-period HEN is first determined and analysed prior to identifying the bottlenecks. As a next step, the design of flexible multi-period HEN is finalised by solving the sub-period debottlenecking model. Several authors have considered the operability issues at the early stages of the process design. The need for this s widely accepted and has motivated the integration of process design and control (IPDC). Narraway and Perkins [15] and Walsh and Perkins [16] were the among the earliest to take into account the general mathematical programming techniques for the simultaneous design and control problem using dynamic process models. Recently, Abu Bakar et al. [17] introduced a new model-based IPDC of HEN which is decomposed into four hierarchical sequential stages. The proposed methodology recommends a solution that satisfies the design, control and economic criteria.
The optimal HEN with unclassified hot/cold process streams was discussed in the work of Kong et al. [18], Quirante et al. [19] and Onishi et al. [20]. All these works depend on the process operating conditions to finalise the classification of process streams. Kong et al. [18] presented mixed-integer nonlinear programming (MINLP) for the heat integration model which accounts for unclassified process streams and variable stream temperatures and flowrates. Quirante et al. [19] later extended the disjunctive model of the pinch location method proposed by Quirante et al. [21] and work by Kong et al. [18] for the simultaneous process optimisation and heat integration. This work is also extended for the isothermal process streams, multiple utilities and area estimation of HEN. The area estimation is done based on the vertical heat transfer between the hot and cold balanced composite curves. Onishi et al. [20] proposed an optimisation model to enhance the work and HENs energy efficiency and cost-effective synthesis considering the unclassified process streams. It combined the methods of mathematical programming and pinch location while adjusting the pressure and temperature of unclassified streams.
State-of-the-art studies on HEN flexibility, show that there are a few key limitations associated with the existing methods. In these previous works, the optimal design of HEN synthesis operated under uncertain operating parameters were considered with an appropriate strategy for control and operation. The control variables were assumed to be adjusted during the operation to improve the flexibility of HEN. However, several possible HEN configurations were needed for each of the scenarios in order to increase the flexibility and some exchanger area adjustment was required during the operation. Besides that, heuristics to guide heat exchanger sizing and bypass placement, and which can be applied in all cases have not been introduced. The understanding of how the temperature fluctuation in HEN affects the amount of heat recovery is still not clearly understood.
This work presents a methodology to manage temperature disturbances in HEN design for maximum heat recovery. The impact of the supply temperature fluctuations on utility consumption, heat exchanger sizing and bypass placement are studied to ensure the target temperatures of affected streams are achieved. At the same time, reducing the impact of the fluctuation on downstream heat exchangers and the immediately propagation of the disturbances to heaters or coolers is desired. Where possible, taking advantage of the disturbances for additional heat recovery is also desirable. The plus-minus principle for process changes are used and new heuristics are introduced to guide heat exchanger sizing and bypass placement. Linnhoff and Vredeveld [22] have introduced the plusminus principle for visualising the impact of process modifications on the minimum utility target using the composite curves (CC). Chew et al. [23] applied the plus-minus principle for process modification aimed at maximising energy savings for total site heat integration (TSHI). This methodology enabled designers to identify the potential process changes to maximise energy recovery and reduce utility consumption. Song et al. [24] further implemented the plus-minus principle for inter-plant heat integration (IPHI) for case studies involving threshold problems. The proposed methodology provides a simple technique of rapidly assessing the effect of supply temperature (TS) fluctuations in heat recovery and utility reduction without the need for detailed process simulation.

Methodology
This section describes the methodology that was developed to manage supply temperature disturbances through modification of a conventional heat exchange network. By planning the right size for the heat exchangers, and by utilising the bypass streams, a HEN can be designed with the flexibility to cope with supply temperature disturbances. An illustrative case study of a HEN experiencing supply temperature disturbances on each process stream is used to demonstrate the applicability of the methodology. For the case study, disturbances are assumed to occur at all supply temperatures (TSH1, TSH2, TSH3, TSC1, TSC2 and TSC3) with a deviation of ± 5 ℃ . The manipulated variables for process control include heat exchangers, bypasses and utility flowrates of coolers and heaters. The controlled variables are all the target temperatures (TtH1, TtH2, TtH3, TtC1, TtC2 and TtC3). Table 1 shows the stream data which are used to illustrate the effect of disturbances on maximum energy recovery HEN. There are three hot streams (H1, H2 and H3) and three cold streams (C1, C2 and C3) involved in the process. The required data for the pinch study includes the supply temperature, Ts; target temperature, Tt; heat capacity flowrate, FCp; enthalpy, ∆H and the supply temperature fluctuation temperature range. The minimum temperature difference, ∆Tmin is set as 10 ℃. The maximum energy recovery (MER) targets are determined for the nominal case (without disturbances) by using pinch analysis targeting methods such as the problem table algorithm or composite curves by Linnhoff and Flower [25] or streams temperature vs enthalpy plot (STEP) by Wan Alwi and Manan [26]. The targeted minimum hot utility requirement QHmin is 450 MW and the minimum cold utility requirement QCmin is 180 MW. The hot pinch temperature, TPinch, hot is at 250 ℃ and the cold pinch temperature, TPinch, cold is at 240 ℃. Figure 1 shows the composite curves of the nominal case.

Step 3: Construct the Grid Diagram (GD)
The nominal maximum energy recovery (MER) network is designed using the pinch design method by Linnhoff and Flower [25] and drawn on the grid diagram (GD) shown in Figure 2. Based on the feasibility criteria proposed by Linnhoff and Hindmarsh [27], streams matching above the pinch have to meet the criterion of CPCOLD ≥ CPHOT and streams matching below the pinch have to meet the criterion of CPHOT ≥ CPCOLD, at the pinch location. Equations (1) and (2) are used to determine the enthalpy balances for hot and cold streams.

Step 4: Manage Fluctuating Ts (Ts Disturbances) in HEN to Achieve MER
The steps to manage fluctuating TS in HEN to achieve Tt and MER can be defined in three stages. The first stage describes the effect of increasing or decreasing TS on the hot and cold streams energy requirement. The second stage provides the fundamental theory of the plus-minus principle. The final stage explains the effect of TS fluctuation on utilities based on the plus-minus principle. Heuristics are proposed for each disturbance scenario that necessitates bypass placement and valve opening, as well as the correct heat exchanger sizing to maximise utility savings. Figure 3 shows the effect of increasing or decreasing TS for a hot stream. The target temperature, Tt is assumed to be maintained. An example of this situation is a fluctuating reactor exit temperature. Before the stream enters the next unit operation, the Tt should be maintained at the setpoint value. From Figure 3, it can be seen that increasing TS for hot stream results in an increase in enthalpy, ∆H, while decreasing TS results in a decrease in ∆H. On the other hand, Figure 4 shows the effect of increasing or decreasing TS on a cold stream. From the figure, it can be seen that an increase in TS results in an increase in ∆H while a decrease in TS results in a decrease in ∆H. The target temperature, Tt is also assumed to be maintained. An example of this situation is a feed stream coming from a storage tank experiencing temperature fluctuations due to changes in ambient conditions as a result of weather changes in a four-season country.  Figure 5 shows the effect of plus-minus principle on composite curves (CC), as explained by Linnhoff and Vredeveld [22]. Above the pinch, increasing the enthalpy of a hot stream and decreasing the enthalpy of a cold stream decreases the minimum hot utility target, QH,min. Doing the reverse above the pinch increases QH,min. Below the pinch, increasing the enthalpy of a cold stream and decreasing the enthalpy of a hot stream decreases the minimum cold utility target, QC,min. Doing the reverse below the pinch increases QC,min.

Control Mechanism Decision by Using the Plus-Minus Principle
The plus-minus principle can be used to describe the effect of disturbances on HEN operation, and how HEN can be controlled. As shown in Step 4(i), the increase or decrease in supply temperature, TS of the hot or cold stream results in an increase in either QH or QC; or a decrease in QH or QC. For effective control of Tt, utility units such as heaters and coolers have been widely used in process plants. Utility heaters/coolers are typically placed after a series of heat exchangers to supplement stream heating or cooling. A bypass, on the other hand, functions as a mechanism for controlling process stream parameters including for disturbance rejection.
A bypass placed on the hot stream side can be described using Equation (3) while a bypass placed on the cold stream side can be described using Equation (4) [28]. T represents the target temperature of the heat exchanger at hot the stream side, T represents the target temperature of the heat exchanger at a cold stream side; T represents the outlet temperature of the heat exchanger at hot stream side, T represents the outlet temperature of the heat exchanger at a cold stream side; T represents the supply temperature of the heat exchanger at hot stream side, T represents the supply temperature of heat exchanger at a cold stream side; u represents the bypass fraction; u represents the stream fraction at hot stream side, u represents the stream fraction at cold stream side. Application of Equations (1) and (2) is explained next, in the context of heuristics development.
New heuristics have been introduced in this work as guides for the appropriate placement of a bypass and sizing of the heat exchanger that can reduce hot and cold utilities in cases with recurring Ts disturbances. Each proposed heuristic shall refer to the plus-minus principle of process changes. Applicability of the heuristics is explained using a case study.
Below the pinch: Observation 1. Decreasing TS for a hot or cold stream below the pinch results in decreasing QC.
The first observation states that decreasing Ts for a hot or cold stream located below the pinch results in QC decreasing. For example, the TS of cold stream C1 decreases from 40 ℃ to 35 ℃. As shown in Figure 4, decreasing TS on cold stream C1 increases C1 enthalpy by 15 MW. Based on the plus-minus principle, since stream C1 is located below the pinch, it can be used to recover more energy. This results in a reduction of CU2 cold utility from 180 MW to 165 MW, which is desirable. In order to control the target temperature of stream C1 at 220 ℃, the C1 flowrate entering HE5 is selected as a manipulated variable. The duty of HE5 is increased from 20 MW to 35 MW to allow more heat to be exchanged. To be able to do this, HE5 should be designed with a bigger area to accommodate up to 35 MW heat duty.
Observation 2. Below the pinch, the point where QC equals zero is the limit for hot or cold stream TS to decrease. Further decrease in Ts leads to a penalty in QH.
It is observed that there exists a limit for TS to decrease for the hot or cold stream below the pinch, i.e., at the point of zero Qc. Any additional decrease in TS for hot or cold stream leads to a penalty of QH. Figures 6a and 6b illustrate Observation 2 involving the disturbance scenario for HE4. It can be seen for the case that when Ts of stream H2 decrease, the enthalpy decrease exceeds the MER of cold stream C1, leading to a penalty of hot utility (see Figure 6b). This situation is also illustrated using the plus-minus principle shown in Figure 7. Moreover, decreasing TS for the hot stream at the pinch point or where the heat exchanger inlet has ∆T = ∆Tmin also incurs a penalty at the other side stream. However, in this case, cold stream C1 does not end at the pinch and the difference of temperature between TS,H2 and Tt,C1 is 25 ℃, which is more than Tmin.   Based on this third observation, increasing TS for a hot or cold stream located below the pinch causes QC to increase. To illustrate this, consider an increase in TS for the cold stream C1 from 40 ℃ to 45 ℃ due to a disturbance. Figure 4 shows that the increase in TS for the C1 resulted in the enthalpy for C1 decreasing by 15 MW. Since the Ts of stream C1 is located below the pinch, based on the plus-minus principle, the increase in Ts resulted in an increase in CU2 from 180 MW to 195 MW as shown in Figure 8. For this scenario, the bypass stream is selected as the manipulated variable as this deviation means the heat duty of heat exchanger HE5 is not high enough to keep the Tt at 220 ℃. Equation (2) is used to calculate the bypass fraction for the bypass placed on the cold streamside. In this scenario, supply temperature, T of HE5 at the cold stream side is 45 ℃. The decreasing of the heat duty of HE 5 from 20 MW to 5MW caused the target temperature of HE5 at cold stream side, T increased to 46.667 ℃. As the outlet temperature, T of HE5 with maximum duty at cold stream C1 is 56.667 ℃, the bypass is placed on the cold stream side. Now the valve is opened at a bypass fraction, u of 0.857. The bypass is calculated by rearranging Equation (2). As the heat duty of heat exchanger HE5 is not high enough to achieve the final target temperature of H3, the cold utility of cooler C2 is increased to absorb the remaining heat. Observation 4. Below the pinch, size the heat exchanger to achieve the maximum energy recovery target when TS decreases, and the cooler to achieve the higher utility when Ts increases.
Observation 4 states that the heat exchanger size below the pinch should be designed to achieve the maximum energy recovery when TS decreases and the cooler size should be designed to achieve the higher utility when Ts increases. Previously, the impact of the increase or decrease of Ts on cold stream C1 was shown. In the case of decreasing Ts, more heat is allowed to be exchanged. It is preferable to design HE5 with a bigger area to accommodate up to at 35 MW heat duty instead of 20 MW heat duty for the nominal case. CU2 cold utility should also be designed with the duty of 195 MW instead of the nominal case with the duty of 180 MW in order to absorb the remaining heat when Ts increases. Figures 9a to 9c illustrate this situation by using the plus-minus principle. The bypass stream is used to control the duty of HE5 during the disturbances as shown in Figures 10a to 10c, while the CU2 cold utility is used to absorb the remaining heat. Above the pinch: Observation 5. Increasing TS for a hot or cold stream above the pinch results in decreasing QH. Observation 5 states that, increasing TS for hot or cold stream located below the pinch results in QH decreasing. For example, TS of hot stream H1 increases from 310 ℃ to 315 ℃. As shown in Figure  3, the increasing TS on hot stream H1 increases H1 enthalpy by 15 MW. Based on the plus-minus principle, since stream H1 is located above the pinch, it can be used to recover more energy. This would result in a reduction of HU1 hot utility from 130 MW to 115 MW, which is desirable (as shown in Figure 11. In order to control the target temperature of stream H1 at 270 ℃, the H1 flowrate entering HE2 is selected as a manipulated variable. The duty of HE2 is increased from 120 MW to 135 MW to allow more heat to be exchanged. In order to attain it, HE2 should be designed with a bigger area to accommodate up to 135 MW heat duty. Observation 6. Above the pinch, the point where QH equals zero is the limit for hot or cold stream TS to increase. Further increases in TS lead to a penalty of QC. As explained in Observation 2, the same situation can also occur above the pinch. There is a limit for TS to increase for hot or cold streams above the pinch, i.e., at the point of zero QH. Further increases in TS for the hot or cold stream lead to a penalty of QC. Besides that, increasing TS for the hot stream at the pinch or where the heat exchanger inlet has ∆T = ∆Tmin also leads to a penalty at the other side stream. Figure 12 shows that TS increases for cold stream C3 from 240 ℃ to 245 ℃. The nominal value of TS on cold stream C3 and the outlet temperature of HE1 on the hot stream side is at the pinch temperature. The ∆T of HE1 decrease to 5 ℃, which is less than the original ∆Tmin of 10 ℃. This cause the enthalpy of cold stream C3 reduced to 385 MW. HE1 with the duty of 100 MW led the hot utility HU2 minimised to 285 MW (see Figure 12). This situation is also illustrated by the plus-minus principle shown in Figures 13a to 13b. Figure 13a shows the nominal temperature of stream C3 matching with stream H3 with ∆Tmin of 10 ℃. Figure 13b shows that, increasing of TS on cold stream results in ∆T less than ∆Tmin.

Observation 7.
Decreasing TS for a hot or cold stream above the pinch results in increasing QH. As stated in observation 7, decreasing TS for hot or cold stream located above the pinch causes QH to increase. To illustrate this, consider a decrease in TS for the hot stream H1 from 310 ℃ to 305 ℃. Figure 3 shows that the decrease in TS for the H1 resulted in the enthalpy for H1 to decrease by 15 MW. Since the TS of stream H1 is located above the pinch, based on the plus-minus principle, the decrease in TS resulted in an increase in hot utility HU1 from 130 MW to 145 MW as shown in Figure  14. For this scenario, the bypass stream is selected as a manipulated variable as this deviation means the heat duty of heat exchanger HE2 is not high enough to keep the target temperature at 270 ℃. Equation (1) is used to calculate the bypass fraction for the bypass placed on the hot streamside. In this scenario, the supply temperature, T of HE2 at the hot stream side is 305 ℃. The decrease in heat duty of HE2 from 120 MW to 105 MW caused the target temperature of heat exchanger HE2 at hot stream side, T to decrease to 270 ℃. As the outlet temperature, T of HE2 with maximum duty at hot stream H1 is 260 ℃, the bypass is placed on the hot stream side. Then, the valve is opened at a bypass fraction, u of 0.286. The bypass is calculated by rearranging Equation (1). As the heat duty of heat exchanger HE2 is not high enough to achieve the final target temperature of C2, the hot utility of heater HU1 is increased in order to satisfy the remaining heat. Observation 8. Above the pinch, size the heat exchanger to achieve the maximum energy recovery when TS increase, and the heater to achieve the minimum utility when TS decrease.
Observation 8 states that the heat exchanger size above the pinch should be designed to achieve the maximum energy recovery when TS increases, and the heater size should be designed to achieve the higher utility when TS decreases. Previously, the impact of increase or decrease in Ts on hot stream H1 was shown. In the case of increasing Ts, more heat is allowed to be exchanged. It is preferable to design HE2 with a bigger area to accommodate up to 135 MW heat duty instead of 120 MW for the nominal case. HU1 hot utility should also be designed with the bigger capacity of 145 MW heat duty instead of the nominal case of 130 MW, in order to cater for the remaining heat when TS decreases. Figures 15a-c illustrate this situation by using the plus-minus principle. The bypass stream is used to control the duty of HE2 when a disturbance occurs as shown in Figures 16a-c, while HU1 hot utility is used is to satisfy the remaining heat.  (2) Heuristic 2: Heat exchanger sizing to cater for disturbances in MER HEN a. Size a heat exchanger to cater to the highest amount of energy to be exchanged, considering all disturbances scenario.
b. Size a utility heat exchanger to cater to the highest amount of utility needed, considering all disturbances scenario.
After all the heuristics have been applied to the HEN design to cater for all the possible scenarios, there are some issues that must be checked: (1) After adjusting heat exchanger duties, temperature feasibility test should be done for all the affected streams. If temperature infeasibility occurs, designers should consider redistributing the duty to the utilities.
(2) The bypass fraction should be calculated for all possible cases considering the biggest heat exchanger size. Table 2 summarises the effects of cold and hot streams' supply temperature disturbances on hot and cold utilities as explained in the heuristics. None Increase

Final Sizing and Bypass Design for MER HEN with Disturbance on TS
The heuristics proposed are applied for designing a HEN that is flexible to disturbances. Analysis of the effect of disturbances on each stream in the HEN is summarised in Table 3. It is shown that disturbances can cause either positive or negative impacts on the QH and QC of HEN. The bypass fraction is determined based on the heuristics proposed. Figure 17 shows the final HEN with the bypass placed. Each stream applied different heuristics according to the scenario described previously.

Case Studies
Illustrative case studies extracted from the literature are used to verify the applicability and accuracy of the proposed methodology for optimal HEN synthesis, considering the uncertainties at the supply temperature.

Case Study 1
The first example is based on a methanol synthesis process adapted from Kijevčanin et al. [29]. The Case Study 1 has ∆Tmin of 20 ℃. The pinch temperature is at 357.2 ℃. The minimum heating requirement (QH,min) is 1,953.88 kW while the minimum cooling requirement (QC,min) is 3,463.53 kW. Expected variations of ±10 ℃ in the inlet temperature of streams H1 and C3 are assumed to vary from their nominal values. The case study consists of eight hot streams and three cold streams as shown in Table 4. Four scenarios involving the variations of the inlet temperatures between streams H1 and C3 are observed (Table 5). All the scenarios are applied in designing a HEN with the nominal HE area (previous work) and the maximum HE area (this work). Figure 18 shows the HEN at the nominal condition with the bypass placement but with different HE area and bypass fraction. Analysis of the effect of disturbances on streams H1 and C3 for the nominal HE area (previous work) is shown in Table 6 while the maximum HE area (this work) is shown in Table 7. It is observed that the duty of HE2 increased to 319.70 kW when the supply temperature of hot stream H1 located across or above the pinch increased. This decreased the duty of hot utility HU2 to 181.80 kW. However, if the HE2 area is maintained at the nominal size, the duty of hot utility is also maintained at the nominal value. On the other hand, the duty of HE2 for both works decreased to 227.70 kW when the supply temperature of hot stream H1 is decreased. This led to an increase of hot utility HU2 to 273.80 kW in both works. For the cold stream C3 located below the pinch, it is observed that the duty of HE5 decreased to 186.60 kW when the supply temperature is increased for both works. Both designs required large HU2 utility loads to cope with the deficit enthalpy of HE5. On the contrary, the duty of HE5 can be increased to 256.60 kW, at the maximum HE area when the supply temperature of C3 decreased. Consequently, the minimum cold utility for CU2 can be obtained. The impact of the changes on the economics for the cases of nominal HE area and maximum HE area were analysed by comparing the annualised capital and utility costs for both HEN using Equation (5) [30]. The basic rule to target for cost-effective minimum utilities is to maximise the use of higher temperature cold utilities and lower temperature hot utilities. The type of utilities suggested based on temperature interval is shown in Table 8. The rates for the utilities refer to Sun et al. [31].
Annualised capital cost = Annualised factor × (1,300 + 1,000 A . ) (5) where the annualised factor is 0.298.  Table 9 compares the heat recovery and economic performance of HEN with nominal HE area and maximum HE area for all the scenarios. Results of this study show that the new heuristics can guide the user to manage temperature disturbances in HEN design for maximum heat recovery with minimum total costs. Although this work has a high annualised capital cost due to larger HE area compared to the previous work, the annualised total costs for scenarios A, B and C is still much cheaper than the previous work. For scenario D, the total annualised cost is 0.16% higher due to the large HE area, but the utility load remains the same.

Case Study 2
In this case study, the proposed methodology is applied to solve an illustrative example with three hot streams and three cold streams. The data used for this case study are adapted from the work of Escobar et al. [9]. The nominal data for the problem is listed in Table 10. The expected variations in the inlet temperatures are assumed ± 10 K with ∆Tmin of 10 K. In contrast to Escobar et al. [9], the nominal configuration of HEN is maintained. The maximum heat exchanger area approach is applied to increase the flexibility of HEN due to the uncertainty of operating conditions. The uncertain parameters considered in the design are given in Table 11.  Initially, the MER for the nominal case (without disturbances) is determined by using pinch analysis targeting methods. For this case study, the nominal condition corresponds to the first iteration. By using the same (nominal) HEN configuration, in order for the target temperature to be achieved, the enthalpy for the cold streams is required to be increased while the enthalpy for the hot streams is decreased. Thus, the utility consumption (cold utility) can be reduced. At the same time, the area of the heat exchangers is required to be at the maximum size for the HEN design to be feasible. Figure 19 shows the HEN at the nominal condition with the bypass placement. On the other hand, Escobar et al. [9] suggested two different HEN designs for each iteration as the nominal design is not feasible for the variations up to 10 K in the inlet temperatures. The steps where the critical point is added to the nominal conditions and the multi-period optimisation problem needs to be repeated until the flexibility is accomplished. The new HEN configuration is designed with high flexibility. The TAC is comprised of the annualised utility cost and annualised capital cost. The TAC of this work for the first iteration is much higher than the previous work as the same HEN design with maximum heat exchanger area is applied. However, the TAC for this work with the consideration of uncertainty is $25,986.93/y, comprising $3,028.80/y associated with operating expenses (utility consumption) and $22,958.13/y to capital investment. This work gives 4%/y lower TAC compared to the work of Escobar et al. [9]. Although it has a higher capital cost, the utility consumption is reduced. This method is able to give positive effects even though the HEN has uncertain operating parameters, as it provides better utility usage and would be a good approach for considering a reduction in the environmental impact associated with the use of fossil-based energy sources. The results of this work and those of Escobar et al. [9] are compared in Table 12.

Conclusions
Heat exchanger network configuration can influence disturbance propagation and process behaviour, as well as limit process controllability and operability. A systematic methodology has been developed in this work to manage temperature disturbances and make a heat exchanger network more flexible and operable toward achieving the maximum heat recovery. The method considers the impact of supply temperature fluctuations on utility consumption, heat exchanger sizing and bypass placement. The maximum energy recovery targets are determined for the nominal case (without disturbances) by using pinch analysis targeting methods. The steps to manage the fluctuating supply temperature in the heat exchanger network to achieve target temperature and maximum energy recovery are defined in three stages: (1) the effect of increasing or decreasing supply temperature on the hot and cold streams energy requirement; (2) the fundamental theory of the plus-minus principle; (3) the effect of TS fluctuation on utilities based on the plus-minus principle. From all these steps, key observations were made and new heuristics based on the plus-minus principle of pinch analysis have been introduced for heat exchanger network design in maximising utility savings. The heuristics were applied by simulating various scenarios of disturbances occurring in the process streams of the heat exchanger network. Guidelines on the sizing of heat exchangers and bypass locations and fractions have also been proposed. The grid diagram and temperature vs enthalpy plot were used to illustrate the effect of changes in design and operating parameters. This approach involves a trade-off between the utility cost and capital cost of the affected heat exchangers. Application of the method on two case studies showed that the configuration of the heat exchanger network is maintained for all the scenarios. In addition, the exchanger areas are designed at the maximum size with a bypass to handle the most critical uncertainty of operating conditions while minimising the utility usage. It showed that this method improved the annualised utility cost by up to 89%. Previous works in the literature required the exchanger area adjustment during the operation and have several heat exchanger network configurations for each scenario. Nevertheless, in industry, the configuration and area should be fixed during the operation even with unforeseen uncertainty. The modification of heat exchanger network design during the operation will be costly. This work is beneficial as it has a reduced utility usage and a higher capital cost.
Major novelties introduced by this work include:  The approach is able to handle the most critical uncertainty of operating conditions with the maximum heat exchanger areas while minimising the utility usage.
 The nominal configuration of HEN based on pinch analysis is maintained and controlled by the bypass.
 The effect of uncertain operating conditions on the heat exchange and utility consumption can be easily visualised through the plus-minus principle.
 The proposed new heuristics based on the plus-minus principle of pinch analysis can be applied for all problems.
However, overdesign factors could affect the annualised capital cost of heat exchanger networks. Future research should also include the probability of disturbances occurring in the heat exchanger network. This probability may influence the requirement of the maximum heat exchanger size in the design as well as the total annualised cost.