Next Article in Journal
Studying the Bond Performance of Full-Grouting Rock Bolts Based on the Variable Controlling Method
Next Article in Special Issue
Study of a Gas Turbine Cycle to Boost the Autonomy of Electric Cars
Previous Article in Journal
Impact of Gas Saturation and Gas Column Height at the Base of the Gas Hydrate Stability Zone on Fracturing and Seepage at Vestnesa Ridge, West-Svalbard Margin
Previous Article in Special Issue
Condensation Flow and Heat Transfer Characteristics of R410A in Micro-Fin Tubes and Three-Dimensional Surface Enhanced Tubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthesis of Heat-Integrated Water Networks Using a Modified Heat Exchanger Network Superstructure

1
Faculty of Technology, University of Tuzla, Urfeta Vejzagića 8, 75000 Tuzla, Bosnia and Herzegovina
2
Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova ulica 17, 2000 Maribor, Slovenia
3
Department of Chemical Engineering, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
*
Author to whom correspondence should be addressed.
Energies 2022, 15(9), 3158; https://doi.org/10.3390/en15093158
Submission received: 10 March 2022 / Revised: 16 April 2022 / Accepted: 20 April 2022 / Published: 26 April 2022

Abstract

:
This work presents the synthesis of heat-integrated water networks (HIWNs) by using mathematical programming. A new superstructure is synthesised by combining a water network and a modified heat exchanger network. Based on the proposed superstructure, a mixed-integer nonlinear programming (MINLP) model is developed. The model is solved by using a one-step solution strategy enabling different initialisations and the generation of multiple solutions, from which the best one is chosen. The results show that the proposed model can be effectively used for solving HIWN problems of different complexities, including large-scale problems.

1. Introduction

The sustainable utilisation of natural resources, including water and energy, is very important in our daily life as well as within the industrial sector. This sector uses large amounts of these resources and generates waste streams and emissions that are discharged into the environment. The increased consumption of natural resources, their future scarcity, greenhouse gas emissions and environmental pollution will influence the environment and climate changes. It is necessary to be in line with the adopted targets and sustainable development goals to contribute to the climate and energy framework to achieve economic efficiency and environmental sustainability.
Water and energy are used for various purposes within industrial processes. The minimisation of the consumption of these resources can be achieved by systematically exploring their interconnections in combined water and energy networks or heat-integrated water networks (HIWNs). Various systematic methods (conceptual, mathematical programming and their combination) have been used to achieve this goal. The conceptual methods are based on pinch analysis (PA), and the mathematical programming (MP) methods are based on superstructure optimisation. These methods have been applied to minimise water and energy consumption in various processes. The MP methods, compared to the PA methods, can better address HIWNs, including large-scale problems, multiple contaminants, multiple freshwater sources and gains and losses of water and heat. MP methods have been successfully applied to find the best trade-offs between investment and operating costs in HIWNs.
A work by Budak Duhbacı et al. [1] presented a review of papers considering water and energy minimisation and the improvements achieved in different industrial processes by the application of MP methods. Ahmetović et al. [2] presented a review of systematic methods (PA and MP) along with the the results of case studies in the Kraft pulp mills and reported typical savings in freshwater and energy consumption. The reader is also referred to review papers by Zhang et al. [3], Kermani et al. [4] and Ahmetović et al. [5] for more information about the various methods used for solving HIWN problems.
The focus of several papers have been related to proposing superstructures, models and efficient solution strategies that can be used for solving various HIWN problems, including large-scale ones. Leewongtanawit and Kim [6] proposed a simultaneous design methodology using an MINLP model and solved a large-scale HIWN problem with multiple contaminants and multiple freshwater sources of different temperatures and quality levels. Yang and Grossmann [7] proposed a simultaneous water- and heat-targeting model, combining a water network with a heat-t-argeting model [8] for flowsheet optimisation. This targeting model can be applied to HIWN problems including large-scale problems. Liu et al. [9] proposed a generalised disjunctive programming model (GDP) for the simultaneous inegration of water and energy in HIWNs. This model enables the trade-offs between freshwater and utility consumption and non-isothermal mixing, reducing the complexity of the subsequent heat exchanger network (HEN) design. Ibrić et al. [10] proposed a compact superstructure and an MINLP model for the simultaneous synthesis of HIWNs. The proposed heuristic rules are used for superstructure simplification. The proposed model is solved using a two-step solution strategy that includes targeting and design steps. Hong et al. [11] proposed a systematic simultaneous optimisation approach including a three-step solution strategy that can be used for solving large-scale industrial HIWN problems. In this strategy, a combination of NLP and MINLP models is used to design the HIWN. Ibrić et al. [12] proposed a modified compact superstructure of their previous work [10], an MINLP model and a three-step solution strategy for solving large-scale problems of HIWNs. The modified superstructure includes a simplified heat integration block, enabling non-isothermal mixing and indirect heat exchange within an overall water network with reduced complexity. Recently, Dong et al. [13] proposed a superstructure considering splitters and mixers for every water-using operation and heat exchange unit to simultaneously design HIWNs. The results are presented for a different number of heat exchange units starting from the minimum number of heat exchange units and increasing this number with each iteration.
Based on the published papers, it can be concluded that the majority of them have considered small- and medium-scale HIWN problems, and only several works addressed large-scale HIWN problems. In addition, based on the results reported in the literature for HIWN problems (small-, medium- and large-scale), it can be concluded that the proposed optimal designs have a relatively small number of heat exchangers, compared to traditional HEN designs in which hot and cold process water streams are involved in heat integration. In most cases, there are less than ten heat exchangers in the optimal HIWN designs reported in the literature. One of the reasons for this is that hot and cold water streams in HIWNs can be involved in indirect heat transfers through heat exchangers as well as direct heat transfers by mixing various water streams. Accordingly, due to the mixing of various hot and cold water streams (direct heat transfer), the number of heat exchangers in a final HIWN design can be significantly smaller than a classical HEN design, in which hot and cold process streams are mainly involved in an indirect heat integration.
This work presents a new superstructure and a simultaneous optimisation MINLP model of a heat-integrated water network. In this model, a water network (WN) is combined with a modified HEN model instead of considering a stage-wise HEN superstructure model [14,15] to reduce the overall network complexity and the number of continuous and discrete variables. The stage-wise HEN superstructure [14,15] consists of heat integration stages in which the corresponding stream is split and directed to an exchanger for each potential match between each hot or cold stream. The outlet streams from the exchangers are mixed and the stream is directed to the next stage. The number of stages is usually equal to the maximum number of hot or cold streams. However, the number of possible heat exchange matches increases significantly with the increased number of hot and cold streams. Thus, the overall complexity of the combined HIWN model depends on the identification of hot and cold streams within the WN. In previous research [10], hot and cold streams within the WN were defined with the heating and cooling stages, with a limited number of hot and cold streams directed to the HEN superstructure [14]. The model is later modified to reduce the number of heating/cooling stages [12] using the same HEN superstructure. However, even with the reduced number of hot and cold streams (e.g., a maximum of five), there is still a large number of possible heat-exchange opportunities. Large-scale problems are still difficult to solve and special solution strategies are required, including separate targeting and design steps [12]. The new superstructure proposed in this paper includes a new HEN design consisting of a limited number of heat exchangers and a corresponding number of heaters and coolers, with their interconnections. Additional opportunities for the splitting and mixing of hot and cold water streams are enabled as well as various patterns of heat exchangers (serial, parallel and their combinations) and bypasses. The MINLP model is solved by a one-step solution strategy to minimise the total annualised cost (TAC) of the overall HIWN. This strategy also enables the random generation of initial points, an iterative solution of the MINLP model and the generation of multiple solutions, from which the best one can be chosen. The results obtained for various HIWN problems, including large-scale problems, are in good agreement with the reported literature results.

2. Materials and Methods

2.1. Problem Formulation

Process industries usually require process water that is used in different process operations (process water-using units). Process operations require water of different temperatures and quality, and the water within process operations becomes contaminated in direct contact with the process streams. The wastewater leaving process is discharged into the environment, satisfying a specified temperature constraint. Figure 1 shows a generic representation of the HIWN problem.
The synthesis problem of HIWNs to be solved is stated as follows. The process consists of process water-using units ( p P ) requiring water of different quality and temperature. Within unit p, water becomes contaminated in direct contact with a process stream and, thus, the contaminant ( c C ) is transferred to the water stream. The water is supplied from external freshwater sources ( s S ). Freshwater sources can have different temperatures and contaminant concentrations. Water can be reused from the internal water sources (from process unit P 1 to process unit P p and otherwise) within the water network (WN). To supply water at the required temperatures for the process units p, water heating and cooling is required within the HEN. The HEN consists of heat exchangers ( e E ) enabling the maximum heat recovery between hot and cold water streams; thus, the minimum utility consumption at the specified exchanger minimum approach temperature (EMAT) is enabled also. The hot and cold utilities are available for additional heating and cooling to satisfy the target temperatures of hot and cold streams. The wastewater is discharged into the environment at the specified temperature. It is required to determine the optimum design of the HIWN to minimise the TAC. The network design should include the optimal connections between freshwater sources, process water-using units and heat exchangers and the optimal number of heat exchangers.

2.2. Motivation Example

As a motivation example, a simple problem is presented consisting of two process water-using units and a single contaminant. Freshwater is at temperature of 20 °C and wastewater is discharged into the environment at 30 °C. The data for this example are presented in Table 1, including the mass load of contaminant c in the process unit p ( L p , c ), the maximum inlet x p , c ( i n , m a x ) and outlet x p , c ( o u t , m a x ) concentrations of contaminant c and the operating temperature t p of the process water-using unit p. The cost data required for this problem and all other problems solved in this paper are the same.
If water integration is not considered, freshwater is used in P 1 and P 2 in a total amount of 87.5 kg/s. The water required by water-using units is at 100 °C and 75 °C, and the heating of freshwater is necessary by using external hot utility, if no heat integration is considered. Moreover, to satisfy the temperature constraint of the wastewater, the streams leaving the process water-using units need to be cooled. Thus, the consumption of hot and cold utilities is 25,462.5 kW and 21,787.5 kW, respectively. The basic design, without the water and heat integration, is given in Figure 2.
To reduce freshwater use, water reuse opportunities between process water-using units need to be considered. The water leaving P 1 with contaminant concentration 100 ppm can be partially reused in P 2 when mixed with freshwater to reduce contaminant concentration to the 50 ppm required at the P 2 inlet. It is clear from Figure 2 that the energy of the hot streams leaving the process units can be used to a certain degree, depending on the temperature approach, to preheat cold freshwater streams. Thus, the heat integration opportunities can also be considered to maximise heat recovery. Performing an optimisation of the HIWN using mathematical programming by minimising the TAC, an optimal design shown in Figure 3 can be obtained with reduced freshwater consumption (70 kg/s vs. 87.5 kg/s). The freshwater consumption is reduced due to water reuse P 1 P 2 in the amount of 20 kg/s. The utilities’ consumption is significantly lower for the hot (2940 kW vs. 25,462.5 kW) and the cold utility (0 kW vs. 21,787.5 kW), due to heat integration. The energy of the effluent stream (14,112 kW) is utilised to preheat freshwater from 20 °C to 68 °C. Additionally, the energy of the stream leaving P 1 (3780 kW) is used to additionally preheat freshwater from 68 °C to 86 °C. The hot utility (2940 kW) is required to achieve the operating temperature of P 1 (100 °C). An investment cost of heat exchangers for the integrated design is somewhat higher (USD 248,188/year vs. USD 229,751/year) than the base case due to higher heat exchange areas resulting from a lower temperature difference at the hot and the cold sides of the heat exchangers. The TAC of the integrated design is USD 2,112,569/year, obtained by the later-proposed model and solution approach. The solution obtained in the literature [16] by using the MINLP model (538 equations, 427 continuous variables, 92 discrete variables) proposed by Ahmetović and Kravanja [17] and the global optimisation solver BARON exhibits the same network design. The proposed model for this example consists of 118 single equations, 228 continuous variables and 9 discrete variables.

2.3. Superstructure

Figure 4 shows a generic design of the proposed HIWN superstructure. The superstructure consists of:
  • freshwater sources s S = S 1 , S 2 , S s ;
  • proces water-using units p P = P 1 , P 2 , P p ;
  • heat exchanger units e E = E 1 , E 2 , E e , with the corresponding heaters and coolers;
  • wastewater, discharged into the environment.
The elements of the superstructure within each of the blocks consist of mixers and/or splitters that enable their connections within and between the blocks. These connections enable multiple options for water reuse. Options for heat exchange include direct heat exchange through the non-isothermal mixing, as well as indirect heat exchange within heat exchangers, heaters and coolers. To each heat exchanger e E , one heater and one cooler are assigned for additional heating and cooling of the water streams.

2.4. Model

This section describes the proposed model for the HIWN superstructure shown in Figure 4. The model consists of the mass and heat balance of the WN and HEN superstructure elements as well as additional heat transfer constraints within the HEN. The model is developed under the following assumptions:
  • Temperature and contaminant concentrations of the freshwater source s are constant;
  • The temperature at the inlet/outlet of the process unit p is constant;
  • The specific heat capacity of the water streams is constant;
  • The flowrate of water through the process unit p is constant, under the assumption that the mass of the contaminants in the water streams is much lower than the overall mass of the water streams;
  • The process unit p operates at constant mass load of contaminants transferred to the water stream;
  • The individual heat transfer coefficients of hot and cold water streams and utilities are constant;
  • Heat losses to the environment are neglected;
  • Single hot and cold utilities are available at constant inlet/outlet temperatures;
  • Counter-current heat exchange is assumed;
  • Water is at liquid state, and no phase change occurs for the water streams;
  • The process operates continuously in steady state conditions.

2.4.1. Freshwater Splitter

The freshwater splitter with the corresponding streams is shown in Figure 5. The mass balance of the freshwater splitter s S is given by Equation (1).
m s = p m s , p + e m s , e ( h s ) + e m s , e ( c s ) + e m s , e ( H ) + e m s , e ( C ) , s S

2.4.2. Water-Using Unit

The process water-using unit p P , with the corresponding water streams, is shown in Figure 6. The overall mass balance and the mass balance of contaminant c C for the process unit p P mixer is given by Equations (2) and (3).
m p ( i n ) = s m s , p + e m e , p ( h s ) + e m e , p ( c s ) + p p p m p , p + e m e , p ( C ) + e m e , p ( H ) , p P
m p ( i n ) x p , c ( i n ) = s m s , p x s , c + e m e , p ( h s ) x e , c ( h s ) + e m e , p ( c s ) x e , c ( c s ) + p p p m p , p x p , c ( o u t ) + e m e , p ( C ) x e , c ( C ) + e m e , p ( H ) x e , c ( H ) , p P , c C
The heat balance of the process unit p mixer is given by Equation (4).
m p ( i n ) t p ( i n ) = s m s , p t s + e m e , p ( h s ) t e ( h s , o u t ) + e m e , p ( c s ) t e ( c s , o u t ) + p p p m p , p t p ( o u t ) + e m e , p ( C ) t e ( C , o u t ) + e m e , p ( H ) t e ( H , o u t ) , p P
Under the assumption that the mass of contaminants in the water streams is much lower than the overall mass, the flow-rate of water at the inlet and outlet of the process water-using unit p is equal, as given by Equation (5). However, the concentration of contaminant c is increased as the mass load of contaminant L p , c is transferred to the water stream within the process unit p, as given by Equation (6).
m p ( i n ) = m p ( o u t ) , p P
m p ( i n ) x p , c ( i n ) + L p , c = m p ( o u t ) x p , c ( o u t ) , p P , c C
The mass balance of the process unit p splitter is given by Equation (7).
m p ( o u t ) = p p p m p , p + e m p , e ( h s ) + e m p , e ( c s ) + e m p , e ( C ) + e m p , e ( H ) + m p ( e n v ) , p P

2.4.3. Heat Exchangers

Figure 7 shows a counter-current heat exchanger with the corresponding hot and cold streams. The exchanger e E consists of an exchanger hot stream mixer and splitter as well as a cold stream mixer and splitter. The selection of mixers in front of heat exchangers enables the possibility that any water stream within the network can be directed towards any heat exchanger as a hot and/or cold stream. The number of heat exchangers is set arbitrarily and the existence of multiple heat exchangers enables many options for heat exchange (serial, parallel, bypass and their combination), as there are interconnections between individual heat exchangers. The overall mass balance and the mass balance for each contaminant c C for the exchanger e E hot stream mixer is given by Equations (8) and (9).
m e ( h s ) = p m p , e ( h s ) + s m s , e ( h s ) + e e e m e , e ( h s ) + e e e m e , e ( c s , h s ) , e E
m e ( h s ) x e , c ( h s ) = p m p , e ( h s ) x p , c ( o u t ) + s m s , e ( h s ) x s , c + e e e m e , e ( h s ) x e , c ( h s ) + e e e m e , e ( c s , h s ) x e , c ( c s ) , e E , c C
The heat balance of the exchanger e E hot stream mixer is given by Equation (10).
m e ( h s ) t e ( h s , i n ) = p m p , e ( h s ) t p ( o u t ) + s m s , e ( h s ) t s + e e e m e , e ( h s ) t e ( h s , o u t ) + e e e m e , e ( c s , h s ) t e ( c s , o u t ) , e E
The mass balance of exchanger e E hot stream splitter is given by Equation (11).
m e ( h s ) = m e ( h s , e n v ) + p m e , p ( h s ) + e e e m e , e ( h s ) + e e e m e , e ( h s , c s ) + e m e , e ( h s , C ) , e E
The overall mass balance and the mass balance for each contaminant c C for the exchanger e E cold stream mixer is given by Equations (12) and (13). The heat balance of the exchanger e E cold stream mixer is given by Equation (14).
m e ( c s ) = p m p , e ( c s ) + s m s , e ( c s ) + e e e m e , e ( c s ) + e e e m e , e ( h s , c s ) , e E
m e ( c s ) x e , c ( c s ) = p m p , e ( c s ) x p , c ( o u t ) + s m s , e ( c s ) x s , c + e e e m e , e ( c s ) x e , c ( c s ) + e e e m e , e ( h s , c s ) x e , c ( h s ) , e E , c C
m e ( c s ) t e ( c s , i n ) = p m p , e ( c s ) t p ( o u t ) + s m s , e ( c s ) t s + e e e m e , e ( c s ) t e ( c s , o u t ) + e e e m e , e ( h s , c s ) t e ( h s , o u t ) , e E
The mass balance of exchanger e E cold stream splitter is given by Equation (15).
m e ( c s ) = m e ( c s , e n v ) + p m e , p ( c s ) + e e e m e , e ( c s ) + e e e m e , e ( c s , h s ) + e m e , e ( c s , H ) , e E

2.4.4. Coolers and Heaters

According to the proposed superstructure shown in Figure 4, to each heat exchanger e E , one heater (H) is assigned for additional heating of the cold streams and one cooler (C) for additional cooling of the hot streams. Figure 8 shows the cooler and heater with their corresponding streams. The overall mass balance and the mass balance for each contaminant c C for the cooler e E mixer is given by Equations (16) and (17). The heat balance of the cooler e E mixer is given by Equation (18).
m e ( C ) = p m p , e ( C ) + s m s , e ( C ) + e e e m e , e ( C ) + e m e , e ( h s , C ) , e E
m e ( C ) x e , c ( C ) = p m p , e ( C ) x p , c ( o u t ) + s m s , e ( C ) x s , c + e e e m e , e ( C ) x e , c ( C ) + e m e , e ( h s , C ) x e , c ( h s ) , e E , c C
m e ( C ) t e ( C , i n ) = p m p , e ( C ) t p ( o u t ) + s m s , e ( C ) t s + e e e m e , e ( C ) t e ( C , o u t ) + e m e , e ( h s , C ) t e ( h s , o u t ) , e E
The mass balance of the cooler e E splitter is given by Equation (19).
m e ( C ) = p m e , p ( C ) + m e ( C , e n v ) + e e e m e , e ( C ) , e E
The overall mass balance and the mass balance for each c C for the heater e E mixer is given by Equations (20) and (21). The heat balance of the heater e E mixer is given by Equation (22).
m e ( H ) = p m p , e ( H ) + s m s , e ( H ) + e e e m e , e ( H ) + e m e , e ( c s , H ) , e E
m e ( H ) x e , c ( H ) = p m p , e ( H ) x p , c ( o u t ) + s m s , e ( H ) x s , c + e e e m e , e ( H ) x e , c ( H ) + e m e , e ( c s , H ) x e , c ( c s ) , e E , c C
m e ( H ) t e ( H , i n ) = p m p , e ( H ) t p ( o u t ) + s m s , e ( H ) t s + e e e m e , e ( H ) t e ( H , o u t ) + e m e , e ( c s , H ) t e ( c s , o u t ) , e E
The mass balance of the heater e E splitter is given by Equation (23).
m e ( H ) = p m e , p ( H ) + m e ( H , e n v ) + e e e m e , e ( H ) , e E

2.4.5. Wastewater Mixer

The streams entering the wastewater mixer are those leaving all the splitters in the water network except the freshwater splitter, as shown in Figure 9. The overall mass balance of the wastewater mixer and the mass balance for the contaminant c C are given by Equations (24) and (25). The heat balance of the wastewater mixer is given by Equation (26).
m ( e n v ) = p m p ( e n v ) + e m e ( c s , e n v ) + e m e ( h s , e n v ) + e m e ( C , e n v ) + e m e ( H , e n v ) , e E
m ( e n v ) x c ( e n v ) = p m p ( e n v ) x p , c ( o u t ) + e m e ( c s , e n v ) x e , c ( c s ) + e m e ( h s , e n v ) x e , c ( h s ) + e m e ( C , e n v ) x e , c ( C ) + e m e ( H , e n v ) x e , c ( H ) , e E , c C
m ( e n v ) t ( e n v ) = p m p ( e n v ) t p ( o u t ) + e m e ( c s , e n v ) t e ( c s , o u t ) + e m e ( h s , e n v ) t e ( h s , o u t ) + e m e ( C , e n v ) t e ( C , o u t ) + e m e ( H , e n v ) t e ( H , o u t ) , e E

2.4.6. Global Mass Balance

According to Karuppiah and Grossmann [18], the global mass balance equations can improve the lower bounds of the relaxed model. The overall mass balance and the contaminant mass balance are given by Equations (27) and (28).
s m s = m ( e n v )
s m s x s , c + p L p , c = m ( e n v ) x c ( e n v ) , c C

2.4.7. Heat Transfer Constraints

Within the heat exchanger e E , heat transfer occurs between the hot and cold water streams. The heat allocated from the hot water stream ( q e ) is given by Equation (29). The same amount of heat ( q e ) is allocated to the cold water stream as the heat losses are neglected, as given by Equation (30).
q e = m e ( h s ) c p t e ( h s , i n ) t e ( h s , o u t ) , e E
q e = m e ( c s ) c p t e ( c s , o u t ) t e ( c s , i n ) , e E
The heat balance of the cooler and heater are given by Equations (31) and (32).
q e ( C ) = m e ( C ) c p t e ( C , i n ) t e ( C , o u t ) , e E
q e ( H ) = m e ( H ) c p t e ( H , o u t ) t e ( H , i n ) , e E
By assuming counter-current heat exchangers with the upper bound on temperature differences Γ , the temperature differences at the hot and cold sides of the heat exchangers, the coolers and the heaters are given by Equations (33)–(38).
Δ t h e t e ( h s , i n ) t e ( c s , o u t ) + Γ 1 z e , e E
Δ t c e t e ( h s , o u t ) t e ( c s , i n ) + Γ 1 z e , e E
Δ t h e ( C ) t e ( C , i n ) t ( C U , o u t ) + Γ 1 z e ( C ) , e E
Δ t c e ( C ) t e ( C , o u t ) t ( C U , i n ) + Γ 1 z e ( C ) , e E
Δ t h e ( H ) t ( H U , i n ) t e ( H , o u t ) + Γ 1 z e ( H ) , e E
Δ t c e ( H ) t ( H U , o u t ) t e ( H , i n ) + Γ 1 z e ( H ) , e E
Equations (39)–(42) describe the feasibility of temperatures of water streams entering and leaving the heat exchangers, heaters and coolers.
t e ( h s , i n ) t e ( h s , o u t ) , e E
t e ( c s , i n ) t e ( c s , o u t ) , e E
t e ( H , i n ) t e ( H , o u t ) , e E
t e ( C , i n ) t e ( C , o u t ) , e E
With the binary variables z e , z e ( C ) , z e ( H ) denoting the existence of heat exchangers, coolers and heaters, the heat transfer is constrained by the upper bounds by the following Equations (43)–(45).
q e q e ( u p ) z e 0 , e E
q e ( C ) q e ( C , u p ) z e ( C ) 0 , e E
q e ( H ) q e ( H , u p ) z e ( H ) 0 , e E

2.4.8. Objective Function of the Model

The objective function of the model minimises the TAC of the network consisting of freshwater cost, hot and cold utility cost and the investment cost of HEN, as given by Equation (46)
min T A C = s m s C F W s f + e q e ( H ) C H U + e q e ( C ) C C U + e a z e + b A e n + e a z e ( C ) + b A e ( C ) n + e a z e ( H ) + b A e ( H ) n
where the heat exchanger and the cooler and heater heat exchange area parameters are given as follows:
A e = q e 1 h e ( h s ) + 1 h e ( c s ) Δ t h e Δ t c e Δ t h e + Δ t c e 2 1 / 3
A e ( C ) = q e ( C ) 1 h e ( h s ) + 1 h C U Δ t h e ( C ) Δ t c e ( C ) Δ t h e ( C ) + Δ t c e ( C ) 2 1 / 3
A e ( H ) = q e ( H ) 1 h e ( c s ) + 1 h H U Δ t h e ( H ) Δ t c e ( H ) Δ t h e ( H ) + Δ t c e ( H ) 2 1 / 3
Chen’s [19] approximation of the logarithmic mean temperature difference was used in Equations (47)–(49) to avoid numerical problems with the logarithmic term. Chen’s approximation slightly overestimates the heat exchange area requirements [20], on average by 0.025% in the studied examples.

2.5. Solution Strategy

The model is implemented in the General Algebraic Modelling System (GAMS) v. 24.6.1 [21] and solved on a laptop computer with a 2.80 GHz processor (4 cores) and 16 GB RAM. The proposed model is solved by using a one-step iterative strategy. For the solution of the MINLP model, the SBB solver is used with CONOPT as NLP solvers for the root node and sub-nodes. Multiple solutions are generated by solving the same model iteratively by generating different initialisation points for all variables on each solve run. Initialisation points are randomly generated for variables between zero and the upper bound for each variable using the GAMS function uniform and generating different initial values in each solve run (execseed = gmillisec(jnow)). Generic upper bounds are obtained based on the initial data from the problem formulation, as proposed by Ibrić et al. [10].

3. Results and Discussion

To demonstrate the capabilities of the proposed model, five examples are solved in increasing order of complexity, concerning a number of process water-using units, contaminants and freshwater sources. The following data are used for the solutions of the problems.
Freshwater: The freshwater is at a temperature of either 20 °C or 30 °C, with a unit cost of USD 0.375/t. The specific heat capacity of all water streams is 4.2 kJ/(kgK). The freshwater is free of contaminants.
Utility: The hot utility is fresh steam at either 120 °C or 150 °C, with unit costs of USD 377/(kW-year) and USD 388/(kW-year). The inlet and outlet temperatures of the hot utility are the same. The cold utility is cooling water at an inlet temperature of 10 °C and an outlet temperature of 20 °C. The unit cost of cooling the water is USD 189/(kW-year).
Heat exchangers: The individual heat transfer coefficients for the hot and cold water streams and the hot and cold utilities are 1 kW/(m2K). The heat exchangers are counter-current with annualised investment cost (USD/year) given as a function of the heat exchange area (m2) as follows: 8000 + 1200 ( A ) 0.6 .
The plant operates continuously 8000 h/year.

3.1. Example 1

This example presents a single-contaminant problem with four process water-using units previously studied in the literature [17,22,23,24]. The freshwater is at 20 °C without contaminants and the wastewater is discharged into the environment at 30 °C. The data for the process water-using units are given in Table 2.
The optimal network design, shown in Figure 10, exhibited the minimum freshwater consumption (25 kg/s) and consumption of only the hot utility (1050 kW). The optimal HEN design consists of two heat exchangers and one heater, with the HEN investment cost of USD 131,525/year. This is a somewhat improved solution compared to the solution obtained by using a more complex stage-wise HEN superstructure combined with the water network superstructure [23]. The same network design and TAC were obtained by Yan et al. [24] by using the NLP model and the global optimisation solver BARON. A comparison of the results with those from the literature for Example 1 is given in Table 3.

3.2. Example 2

In this example, a problem originally solved as an isothermal water network problem [18] is modified as a non-isothermal water network problem, as reported in the literature [25]. The problem includes four process water-using units with two contaminants and a single freshwater source at 20 °C without contaminants. The wastewater is discharged into the environment at 30 °C with no restrictions on contaminant concentrations in the effluent stream. Data for the process water-using units are given in Table 4. The EMAT is 1 °C.
Ibrić et al. [25] reported a solution with the freshwater consumption of 19.444 kg/s and the hot-utlity consumption of 816.67 kW. The network exhibited a design comprising of two heat exchangers and one heater, with the HEN investment cost of USD 188,801/year. The TAC of the network was USD 706,687.7/year. The optimal design shown in Figure 11 exhibited slightly increased consumption of freshwater compared to the literature [25] (20.698 vs. 19.444 kg/s). Consequently, the consumption of hot utility is also increased (869.33 vs. 816.67 kW). However, as a trade-off exist, the HEN investment cost is reduced by ≈41% (USD 111,207 vs. USD 188,801/year). The TAC of the network (USD 662,489/year) is reduced by 6.25% compared to the literature [25].

3.3. Example 3

This example is slightly more complex than the previous one, considering a number of contaminants. The problem consists of four process water-using units and three contaminants originally solved in the literature [22] and later studied by other authors [10,24,26]. The data for the process water-using units are given in Table 5. The freshwater is available at 20 °C without contaminants and the wastewater is discharged into the environment at 30 °C. The EMAT is 1 °C.
The optimal network design obtained by using the model proposed in this paper is shown in Figure 12. Compared to the network design in the literature [10], the basic design, the freshwater and hot utility consumption and the heat integration opportunities are the same. The only difference is related to the water reuse opportunities related to P 4 , due to the increased flowrate in P 4 . Table 6 shows a comparison of the results with those from the literature for Example 3.

3.4. Example 4

The complexity of the studied problems of HIWNs significantly increases with larger number of process water-using units. This example considers a single-contaminant HIWN problem with eight process units. The data for the example are taken from the literature [27] and presented in Table 7. The problem is solved for an EMAT of 1 °C.
In a recent publication, Ibrić et al. [12] solved this example by using a three-step solution strategy. However, the HEN superstructure still exhibits a large number of possible heat-integration opportunities, even with a limited number of hot and cold streams. With the new and modified HEN superstructure proposed in this paper, the HEN combinatorial problem is reduced and a good solution is obtained. Dong et al. [13] proposed a solution to this problem by fixing a number of heat exchangers, heaters and coolers and solving problem for their different numbers. In this way, a set of solutions is obtained with different numbers of exchangers. An optimal network design obtained in this work is shown in Figure 13. The network design exhibited the same heat integration opportunities as in the literature [13] and also the same freshwater and utilities consumption. However, there is a difference in the internal distribution of water streams with very similar network designs. A comparison of the results with those from the literature for is presented in Table 8.

3.5. Example 5

In the previous example, only the hot utility was required; thus, the problem was considered a threshold problem concerning utility consumption. This problem occurs when the temperature of the freshwater source ( t s ) is lower than the temperature of the wastewater ( t ( e n v ) ) and the maximum temperature of the process water-using units ( t p ( i n , m a x ) ). If t s t ( e n v ) , a cooling utility is also required and, thus, the problem is considered to be a pinched HIWN problem. In this example, a single-contaminant pinched problem is considered with fifteen process water-using units. The temperature of the freshwater source is 30 °C and is the same as the temperature of the wastewater discharged into the environment. The data for the process water-using units are taken from the literature [9] and given in Table 9. The hot utility temperature is 150 °C. The EMAT is 10 °C.
Figure 14 shows an optimal network design for Example 5. The optimal design exhibited a minimum freshwater consumption of 100 kg/s as well and a minimum consumption of hot (4200 kW) and cold (4200 kW) utility. When compared to a recent work [12], the hot and cold utility consumption is reduced (4200 kW vs. 4246.2 kW) at the expense of increasing heat exchanger investment costs (USD 352,273/year vs. USD 339,573/year). Compared to the solution presented by Dong et al. [13], the freshwater and utilities consumption are the same. Moreover, the heat exchangers’ heat load and heat exchange areas are identical, as well as for an additional cooler. However, the distribution of heat loads and heat exchange areas for the heaters are different, causing a decrease in HEN investment (USD 352,273/year vs. USD 360,448/year). A comparison of the results with those from the literature for Example 5 is presented in Table 10.

3.6. Model Statistics

Table 11 shows the basic model statistics, including model sizes and computational times, for the given number of solve iterations. A number of iterations can be increased in search of a better solution if the best one cannot be found in a given smaller number of iterations. Ibrić et al. [12], solved Example 5 including fifteen process water-using units and their proposed model consisted of 645 equations, 1028 continuous variables and 96 discrete variables, with a total computational time of 3924 s in only 3 iterations. The proposed model in this paper, for the same example, requires a reduced number of equations with a slightly increased number of continuous variables (1254 vs. 1028). However, the number of discrete variables is reduced (15 vs. 96). The total computational time is 1596 s for a total of 1000 solve iterations, which is, on average, less than 2 s per iteration.

4. Conclusions

This paper presented a mathematical programming model for the synthesis of heat-integrated water networks with a modified heat exchanger network superstructure. The heat exchanger network consists of a given maximum number of heat exchangers with corresponding heaters and coolers including different heat exchange opportunities. A mixed-integer nonlinear programming model is solved by using a one-step solution strategy and a local solver SBB within GAMS software. Examples of different complexities are solved to demonstrate the capabilities of the model, and the results obtained are in good agreement with those in the literature. The model can be extended, considering practical constraints and different types of heat exchangers, to take into account the maximum and minimum heat exchange areas. Moreover, additional model extensions are necessary to include the piping design, so as to ensure practical water flowrates and to manage and control the network, considering the piping investment costs as well as the pumping costs. The model can be extended to include a set of treatment units to enable wastewater treatment for streams leaving the process water-using units. This introduces additional options for water regeneration reuse and recycling and reduces freshwater consumption. This has been rarely addressed in the literature for large-scale HIWN problems.

Author Contributions

Conceptualisation, N.I. and E.A.; methodology, N.I.; software, N.I.; validation, N.I., E.A., A.N., Z.K. and I.E.G.; formal analysis, N.I.; data curation, N.I.; writing—original draft preparation, N.I., E.A., A.N., Z.K. and I.E.G.; writing—review and editing, N.I., E.A., A.N., Z.K. and I.E.G.; visualisation, N.I.; supervision, E.A., Z.K. and I.E.G.; project administration, E.A. and Z.K.; funding acquisition, E.A. and Z.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency (P2-0414 and J7-1816).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the program for Scientific and Technological Cooperation between the Republic of Slovenia and Bosnia and Herzegovina in the period 2019–2022 (Project: “SISUMP”, No.: 01-1282-1/20). In addition, support from the Federation of Bosnia and Herzegovina (Decision No.: 05-35-3210-1/21) for the project “Optimisation of water and energy consumption in the industry” is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CPUCentral Processing Unit
EMATExchanger Minimum Approach Temperature
GAMSGeneral Algebraic Modelling System
GDPGeneralised Disjunctive Programming
HENHeat Exchanger Network
HIWNHeat-Integrated Water Network
MINLPMixed-Integer Nonlinear Programming
MPMathematical Programming
NLPNonlinear Programming
PAPinch Analysis
TACTotal Annualised Cost
WNWater Network
Indices
cContaminant
eHeat exchanger
pProces water-using unit
sFreshwater source
Sets
CContaminants
EHeat exchangers
PProces water-using units
SFreshwater sources
Superscripts
CCooler
C U Cold utility
c s Cold stream
e n v Environment
HHeater
h s Hot stream
H U Hot utility
i n Inlet
m a x Maximum
o u t Outlet
u p Upper bound
Parameters
aFixed cost for heat exchanger, USD/year
bArea cost coefficient for heat exchanger, USD/(m2-year)
C F W s Freshwater cost, USD/kg
C C U Cold utility cost, USD/(kW-year)
C H U Hot utility cost, USD/(kW-year)
c p Specific heat capacity of water, kJ/(kgK)
fPlant annual operating fraction, s/year
hIndividual heat transfer coefficient, kW/(m2K)
L p , c Mass load of contaminant c in process unit p, kg/s
nCost exponent for heat exchanger
t p ( i n ) Inlet temperature of of the process unit p, °C
t p ( i n ) Outlet temperature of of the process unit p, °C
t s Freshwater temperature of the source s, °C
t ( e n v ) Temperature of wastewater discarged into the environment, °C
x s , c Concentration of contaminant c in freshwater source s, ppm
x p , c ( i n , m a x ) Maximum concentration of contaminant c at the inlet to process
unit p, ppm
x p , c ( o u t , m a x ) Maximum concentration of contaminant c at the outlet of process
unit p, ppm
Γ Upper bound for temperature driving force, °C
Continuous variables
AHeat exchange area, m2
mWater mass flowrate, kg/s
xContaminant concentration, ppm
tStream temperature, °C
qHeat flowrate (load), kW
Δ t h Temperature difference at the hot side of heat exchanger, °C
Δ t c Temperature difference at the cold side of heat exchanger, °C
Δ t c Temperature difference at the cold side of heat exchanger, °C
T A C Total annualised cost, USD/year
Binary variables
zExistence of heat exhanger, -

References

  1. Budak Duhbacı, T.; Özel, S.; Bulkan, S. Water and energy minimization in industrial processes through mathematical programming: A literature review. J. Clean. Product. 2021, 284, 124752. [Google Scholar] [CrossRef]
  2. Ahmetović, E.; Kravanja, Z.; Ibrić, N.; Grossmann, I.E.; Luciana, S. State of the art methods for combined water and energy systems optimisation in Kraft pulp mills. Optim. Eng. 2021, 22, 1831–1852. [Google Scholar] [CrossRef]
  3. Zhang, C.; Ren, C.; Liao, Z.; Sun, J.; Wang, J.; Yang, Y. Recent Progresses on Optimal Design of Heat Integrated Water Allocation Network. China Petrol. Process. Petrochem. Technol. 2021, 23, 69–75. [Google Scholar]
  4. Kermani, M.; Kantor, I.D.; Maréchal, F. Synthesis of Heat-Integrated Water Allocation Networks: A Meta-Analysis of Solution Strategies and Network Features. Energies 2018, 11, 1158. [Google Scholar] [CrossRef] [Green Version]
  5. Ahmetović, E.; Ibrić, N.; Kravanja, Z.; Grossmann, I.E. Water and energy integration: A comprehensive literature review of non-isothermal water network synthesis. Comput. Chem. Eng. 2015, 82, 144–171. [Google Scholar] [CrossRef]
  6. Leewongtanawit, B.; Kim, J.K. Synthesis and optimisation of heat-integrated multiple-contaminant water systems. Chem. Eng. Process. Process Intensif. 2008, 47, 670–694. [Google Scholar] [CrossRef]
  7. Yang, L.; Grossmann, I.E. Water Targeting Models for Simultaneous Flowsheet Optimization. Ind. Eng. Chem. Res. 2013, 52, 3209–3224. [Google Scholar] [CrossRef]
  8. Duran, M.A.; Grossmann, I.E. Simultaneous optimization and heat integration of chemical processes. AIChE J. 1986, 32, 123–138. [Google Scholar] [CrossRef]
  9. Liu, Z.; Luo, Y.; Yuan, X. Simultaneous integration of water and energy in heat-integrated water allocation networks. AIChE J. 2015, 61, 2202–2214. [Google Scholar] [CrossRef]
  10. Ibrić, N.; Ahmetović, E.; Kravanja, Z. Mathematical programming synthesis of non-isothermal water networks by using a compact/reduced superstructure and an MINLP model. Clean Technol. Environ. Policy 2016, 18, 1779–1813. [Google Scholar] [CrossRef]
  11. Hong, X.; Liao, Z.; Sun, J.; Jiang, B.; Wang, J.; Yang, Y. Energy and Water Management for Industrial Large-Scale Water Networks: A Systematic Simultaneous Optimization Approach. ACS Sustain. Chem. Eng. 2018, 6, 2269–2282. [Google Scholar] [CrossRef]
  12. Ibrić, N.; Ahmetović, E.; Kravanja, Z.; Grossmann, I.E. Simultaneous optimisation of large-scale problems of heat-integrated water networks. Energy 2021, 235, 121354. [Google Scholar] [CrossRef]
  13. Dong, X.; Zhang, C.; Peng, X.; Chang, C.; Liao, Z.; Yang, Y.; Sun, J.; Wang, J.; Yang, Y. Simultaneous design of heat integrated water allocation networks considering all possible splitters and mixers. Energy 2022, 238, 121916. [Google Scholar] [CrossRef]
  14. Yee, T.F.; Grossmann, I.E. Simultaneous optimization models for heat integration—II. Heat exchanger network synthesis. Comput. Chem. Eng. 1990, 14, 1165–1184. [Google Scholar] [CrossRef]
  15. Yee, T.F.; Grossmann, I.E.; Kravanja, Z. Simultaneous optimization models for heat integration—III. Process and heat exchanger network optimization. Comput. Chem. Eng. 1990, 14, 1185–1200. [Google Scholar] [CrossRef]
  16. Ahmetović, E.; Kravanja, Z. Effects of the Different Stages of Superstructure Development On the Efficiencies and Designs of Heat-Integrated Process-Water Networks. In Proceedings of the 2012 AIChE Annual Meeting, Pittsburgh, PA, USA, 2 November 2012. [Google Scholar]
  17. Ahmetović, E.; Kravanja, Z. Simultaneous synthesis of process water and heat exchanger networks. Energy 2013, 57, 236–250. [Google Scholar] [CrossRef]
  18. Karuppiah, R.; Grossmann, I.E. Global optimization for the synthesis of integrated water systems in chemical processes. Comput. Chem. Eng. 2006, 30, 650–673. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, J. Comments on improvements on a replacement for the logarithmic mean. Chem. Eng. Sci. 1987, 42, 2488–2489. [Google Scholar] [CrossRef]
  20. Chen, J. Logarithmic mean: Chen’s approximation or explicit solution? Comput. Chem. Eng. 2019, 120, 1–3. [Google Scholar] [CrossRef]
  21. GAMS Development Corporation. General Algebraic Modeling System (GAMS) Release 24.6.1; GAMS Development Corporation: Fairfax, VA, USA, 2016. [Google Scholar]
  22. Bogataj, M.; Bagajewicz, M.J. Synthesis of non-isothermal heat integrated water networks in chemical processes. Comput. Chem. Eng. 2008, 32, 3130–3142. [Google Scholar] [CrossRef]
  23. Ibrić, N.; Ahmetović, E.; Kravanja, Z. Two-step mathematical programming synthesis of pinched and threshold heat-integrated water networks. J. Clean. Product. 2014, 77, 116–139. [Google Scholar] [CrossRef]
  24. Yan, F.; Wu, H.; Li, W.; Zhang, J. Simultaneous optimization of heat-integrated water networks by a nonlinear program. Chem. Eng. Sci. 2016, 140, 76–89. [Google Scholar] [CrossRef]
  25. Ibrić, N.; Ahmetović, E.; Kravanja, Z. Simultaneous optimization of water and energy within integrated water networks. Appl. Thermal Eng. 2014, 70, 1097–1122. [Google Scholar] [CrossRef]
  26. Ahmetović, E.; Kravanja, Z. Simultaneous optimization of heat-integrated water networks involving process-to-process streams for heat integration. Appl. Thermal Eng. 2014, 62, 302–317. [Google Scholar] [CrossRef]
  27. Bagajewicz, M.; Rodera, H.; Savelski, M. Energy efficient water utilization systems in process plants. Comput. Chem. Eng. 2002, 26, 59–79. [Google Scholar] [CrossRef]
  28. Hong, X.; Liao, Z.; Jiang, B.; Wang, J.; Yang, Y. Simultaneous optimization of heat-integrated water allocation networks. Appl. Energy 2016, 169, 395–407. [Google Scholar] [CrossRef]
  29. Hong, X.; Liao, Z.; Jiang, B.; Wang, J.; Yang, Y. Targeting of heat integrated water allocation networks by one-step MILP formulation. Appl. Energy 2017, 197, 254–269. [Google Scholar] [CrossRef]
Figure 1. A schematic representation of the HIWN problem.
Figure 1. A schematic representation of the HIWN problem.
Energies 15 03158 g001
Figure 2. Base case for the motivation example without water and heat integration.
Figure 2. Base case for the motivation example without water and heat integration.
Energies 15 03158 g002
Figure 3. Optimal design of HIWN.
Figure 3. Optimal design of HIWN.
Energies 15 03158 g003
Figure 4. The proposed HIWN superstructure.
Figure 4. The proposed HIWN superstructure.
Energies 15 03158 g004
Figure 5. Freshwater splitter with corresponding streams.
Figure 5. Freshwater splitter with corresponding streams.
Energies 15 03158 g005
Figure 6. Process water-using unit with corresponding streams.
Figure 6. Process water-using unit with corresponding streams.
Energies 15 03158 g006
Figure 7. Heat exchanger unit with corresponding streams.
Figure 7. Heat exchanger unit with corresponding streams.
Energies 15 03158 g007
Figure 8. Cooler and heater with corresponding streams.
Figure 8. Cooler and heater with corresponding streams.
Energies 15 03158 g008
Figure 9. Wastewater mixer with corresponding streams.
Figure 9. Wastewater mixer with corresponding streams.
Energies 15 03158 g009
Figure 10. Optimal network design for Example 1.
Figure 10. Optimal network design for Example 1.
Energies 15 03158 g010
Figure 11. Optimal network design for Example 2.
Figure 11. Optimal network design for Example 2.
Energies 15 03158 g011
Figure 12. Optimal network design for Example 3.
Figure 12. Optimal network design for Example 3.
Energies 15 03158 g012
Figure 13. Optimal network design for Example 4.
Figure 13. Optimal network design for Example 4.
Energies 15 03158 g013
Figure 14. Optimal network design for Example 5.
Figure 14. Optimal network design for Example 5.
Energies 15 03158 g014
Table 1. Process water-using units data for motivation example.
Table 1. Process water-using units data for motivation example.
Unit L p , c (g/s) x p , c ( in , max ) (ppm) x p , c ( out , max ) (ppm) t p (°C)
P 1 550100100
P 2 305080075
Table 2. Process water-using units data for Example 1.
Table 2. Process water-using units data for Example 1.
Unit L p , c (g/h) x p , c ( i n , m a x ) (ppm) x p , c ( o u t , m a x ) (ppm) t p (°C)
P 1 2000010040
P 2 500050100100
P 3 30,0005080075
P 4 400040080050
Table 3. Comparison of the results with those from the literature for Example 1.
Table 3. Comparison of the results with those from the literature for Example 1.
ReferenceFreshwater (kg/s)Hot Utility (kW)No. of HEsHEN Investment (USD/year)TAC (USD/year)
Bogataj and Bagajewicz [22] *2510503146,748812,598
Ahmetović and Kravanja [17]2510502134,227800,077
Yan et al. [24]2510503131,525797,375
Ibrić et al. [23]2510502134,227800,077
This paper2510503131,527797,377
* Recalculated in Ahmetović and Kravanja [17].
Table 4. Process water-using units data for Example 2.
Table 4. Process water-using units data for Example 2.
Unit L p , c (kg/h) x p , c ( i n , m a x ) (ppm) x p , c ( o u t , m a x ) (ppm) t p (°C)
ABABAB
P 1 11.5002537.540
P 2 115050707050
P 3 11505066.66766.66775
P 4 21505078.57178.571100
Table 5. Process water-using units data for Example 3.
Table 5. Process water-using units data for Example 3.
Unit L p , c (kg/h) x p , c ( i n , m a x ) (ppm) x p , c ( o u t , m a x ) (ppm) t p (°C)
ABCABCABC
P 1 213015010010010040
P 2 50155010030100200250100
P 3 304010010010080075060075
P 4 4221740038025080080080050
Table 6. Comparison of the results with those from the literature for Example 3.
Table 6. Comparison of the results with those from the literature for Example 3.
ReferenceFreshwater (kg/s)Hot Utility (kW)No. of HEsHEN Investment (USD/year)TAC (USD/year)
Bogataj and Bagajewicz [22] *26.5351115.74199,861907,066
Ahmetović and Kravanja [26]26.5351114.54183,063889,772
Yan et al. [24]26.7161122.14183,107894,856
Ibrić et al. [10]26.5351114.53177,859884,595
This paper26.5351114.53177,859884,595
* Recalculated in Ahmetović and Kravanja [26].
Table 7. Process water-using units data for Example 4.
Table 7. Process water-using units data for Example 4.
Unit L p , c (g/s) x p , c ( i n , m a x ) (ppm) x p , c ( o u t , m a x ) (ppm) t p (°C)
P 1 2258040
P 2 2.882590100
P 3 42520080
P 4 35010060
P 5 305080050
P 6 540080090
P 7 240060070
P 8 1010050
Table 8. Comparison of the results with those from the literature for Example 4.
Table 8. Comparison of the results with those from the literature for Example 4.
ReferenceFreshwater (kg/s)Hot Utility (kW)No. of HEsHEN Investment (USD/year)TAC (USD/year)
Bagajewicz et al. [27]125.945289.612NANA
Hong et al. [28]127.7135363.944274,4713,651,247
Hong et al. [29]126.4265325.123274,4713,641,653
Ibrić et al. [23]125.9435289.64273,6603,628,022
Ibrić et al. [12]125.9435289.64271,9073,626,269
Dong et al. [13]126.335305.63257,2453,621,773
This paper126.3255305.63257,2473,621,775
Table 9. Process water-using units data for Example 5.
Table 9. Process water-using units data for Example 5.
Unit L p , c (g/s) x p , c ( i n , m a x ) (ppm) x p , c ( o u t , m a x ) (ppm) t p (°C)
P 1 1010050
P 2 2010080
P 3 40200110
P 4 4.5030090
P 5 0.50100100
P 6 225040070
P 7 2.25150300100
P 8 1.5200250120
P 9 130035080
P 10 1.835038095
P 11 1300350110
P 12 1.235040095
P 13 60200100
P 14 0.3380400100
P 15 0.835080060
Table 10. Comparison of the results with those from the literature for Example 5.
Table 10. Comparison of the results with those from the literature for Example 5.
ReferenceFreshwater (kg/s)Hot Utility (kW)Cold Utility (kW)No. of HEHEN Investment (USD/year)TAC (USD/year)
Liu et al. [9]100420042009492,9003,918,600
Hong et al. [28]100420042005369,3443,872,744
Hong et al. [29]100420042007447,9763,951,376
Ibrić et al. [12]1004246.24246.25339,5733,869,655
Dong et al. [13]100420042006360,4483,863,848
This paper100420042006352,2733,855,672
Table 11. Model statistics for the studied examples.
Table 11. Model statistics for the studied examples.
ExampleNo. of EquationsNo. of Continuous VariablesNo. of Discrete VariablesCPU Time (s)No. of Iterations
1165458124921000
21894411259100
321546612103200
422077815218100
5262125415922500
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ibrić, N.; Ahmetović, E.; Nemet, A.; Kravanja, Z.; Grossmann, I.E. Synthesis of Heat-Integrated Water Networks Using a Modified Heat Exchanger Network Superstructure. Energies 2022, 15, 3158. https://doi.org/10.3390/en15093158

AMA Style

Ibrić N, Ahmetović E, Nemet A, Kravanja Z, Grossmann IE. Synthesis of Heat-Integrated Water Networks Using a Modified Heat Exchanger Network Superstructure. Energies. 2022; 15(9):3158. https://doi.org/10.3390/en15093158

Chicago/Turabian Style

Ibrić, Nidret, Elvis Ahmetović, Andreja Nemet, Zdravko Kravanja, and Ignacio E. Grossmann. 2022. "Synthesis of Heat-Integrated Water Networks Using a Modified Heat Exchanger Network Superstructure" Energies 15, no. 9: 3158. https://doi.org/10.3390/en15093158

APA Style

Ibrić, N., Ahmetović, E., Nemet, A., Kravanja, Z., & Grossmann, I. E. (2022). Synthesis of Heat-Integrated Water Networks Using a Modified Heat Exchanger Network Superstructure. Energies, 15(9), 3158. https://doi.org/10.3390/en15093158

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

Article Metrics

Back to TopTop