Cost-Optimal Heat Exchanger Network Synthesis Based on a Flexible Cost Functions Framework

In this article an approach to incorporate a flexible cost functions framework into the cost-optimal design of heat exchanger networks (HENs) is presented. This framework allows the definition of different cost functions for each connection of heat source and sink independent of process stream or utility stream. Therefore, it is possible to use match-based individual factors to account for different fluid properties and resulting engineering costs. Layout-based factors for piping and pumping costs play an important role here as cost driver. The optimization of the resulting complex mixed integer nonlinear programming (MINLP) problem is solved with a genetic algorithm coupled with deterministic local optimization techniques. In order to show the functionality of the chosen approach one well studied HEN synthesis example from literature for direct heat integration is studied with standard cost functions and also considering additional piping costs. Another example is presented which incorporates indirect heat integration and related pumping and piping costs. The versatile applicability of the chosen approach is shown. The results represent designs with lower total annual costs (TAC) compared to literature.


Introduction
The application of heat integration strategies can have a significant impact on reducing the amount of utility used by a process and thus improve its economic performance.Against the background of increasing global competition, environmental specifications, climate change and assumedly increasing energy costs heat integration using heat exchanger networks (HENs) have a significant importance [1].
Heat integration strategies have been developed to reduce both, capital and operating costs since decades by now.Pinch technology [2] and mathematical programming [3] have been the two main approaches and have been improved numerous times by many researchers [4].Not only single processes but also total site heat integration has been considered.Initial works have been carried out by Dhole and Linnhoff [5].Due to new challenges, in recent years publications have covered relevant practical issues in a higher degree of detail.As a consequence the problem complexity increased.The main issues influencing practical implementation of total site integration have been formulated by Chew et al. [6].The consideration of further impact factors like safety related issues has been a major topic in heat integration during the last years [1].The identification of critical risk equipment and respective streams for total site heat integration was developed by Liu et al. [7].Nemet et al. [8,9] developed approaches for including risk assessment already during the HEN synthesis.Multiperiod HEN synthesis as well as controllability and disturbance propagation have been studied [10][11][12].Due to operational issues and safety concerns direct heat integration is not always practical to realize [13,14].Therefore, Wang et al. [14] developed a graphical methodology Energies 2019, 12, 784 2 of 18 to investigate different connection patterns for total site heat integration.This methodology was developed further by applying mathematical models to determine the optimal solution for multi-plant heat integration [15].Multi-plant heat integration has been further considered by Chang et al. [13,16].The consideration of plant layout issues is an important factor during optimization.Liew et al. [17] introduced an improved heat cascade algorithm considering pressure drop and heat loss for utility targeting in total site heat integration.Pouransari and Maréchal [18] took into account individual priority levels for different possible connections and Souza et al. [19] included pressure drops in piping as well as in heat exchangers.
The review given above shows that the various demands on HEN optimization in the literature are manifold and a huge variety of optimization models are used.In this work our aim is to incorporate a flexible consideration of cost functions into the cost-optimal HEN synthesis to account for various fields of application.The most important part is the definition of individual factors for each possible connection of heat source and sink.Therefore it is possible to represent a multitude of practical implementation requirements with the same mixed integer nonlinear programming (MINLP) model.For example, it is possible to consider individual cost functions for different types of heat exchangers required for different operation conditions and properties of the involved process streams.Depending on the properties of the process streams the material and thus the costs of the installed heat exchanger can vary significantly.Peripheral equipment, layout constraints as well as cost of premises can be taken into account.These cost functions can be directly incorporated without changing the model or the solution algorithm itself.Concerning the algorithm performance it is the clear aim to be able to generate valid network structures that are competitive to the best solutions published in literature by now.Therefore, different approaches are combined.This model was primarily developed for direct heat integration [1], but is also shown to be applicable for the cost optimization of indirect heat integration problems.

Methodology
The utilized simultaneous cost optimization model is mainly based on a superstructure MINLP formulation and a hybrid genetic algorithm developed by Luo et al. [20] to solve the problem resulting from this formulation.A common HEN optimization problem statement is described by a set of N h hot process streams, N c cold process streams, their respective heat capacity flow rates .W and heat transfer coefficients h.The supply and target temperatures (T and T out ) of the process streams are given as well as the hot utility HU and cold utility streams CU with their respective temperature levels.Furthermore, the cost parameters should be known to perform a cost optimization.

Fundamentals
The model for HEN optimization in this work is based on counterflow heat exchangers for the process-to-process heat exchangers as well as the utility heat exchangers.The heat load .Q in kW is calculated as the product of overall heat transfer coefficient U in kW/(m 2 K), the heat exchanger area A in m 2 and the logarithmic mean temperature difference (LMTD) ∆T m in K: .
The overall heat transfer coefficient is calculated via the individual heat transfer coefficients h h and h c in kW/(m 2 K) of the connected hot and cold process streams.The thermal resistance of the wall is neglected: Energies 2019, 12, 784 3 of 18 Due to the consideration of counterflow heat exchangers, the LMTD is used as temperature driving force for heat transfer: The correlations for hot and cold outlet temperatures for the heat exchangers are as follows: and: In order to represent different HEN solutions, a stage-wise superstructure proposed by Yee et al. [21] is used.As an example, a superstructure for two hot process streams j and three cold process streams k with stages i is shown in Figure 1 including the nomenclature used in the following equations.The correlations for hot and cold outlet temperatures for the heat exchangers are as follows: and: In order to represent different HEN solutions, a stage-wise superstructure proposed by Yee et al. [21] is used.As an example, a superstructure for two hot process streams  and three cold process streams  with stages  is shown in Figure 1 including the nomenclature used in the following equations.In each stage of the superstructure each hot stream is possibly connected with each cold stream.At the end of the streams a utility heat exchanger may be placed to ensure the desired target temperatures to get reached.Within the superstructure every heat exchanger has a specific index  depending on the positioning: Rearranging Equations ( 1)-( 5) according to the outlet temperatures leads to a formulation which allows for the numerical calculation of the HEN temperatures: with  , being: and: In each stage of the superstructure each hot stream is possibly connected with each cold stream.At the end of the streams a utility heat exchanger may be placed to ensure the desired target temperatures to get reached.Within the superstructure every heat exchanger has a specific index ijk depending on the positioning: Rearranging Equations (1)-( 5) according to the outlet temperatures leads to a formulation which allows for the numerical calculation of the HEN temperatures: Energies 2019, 12, 784 4 of 18 with R h,ijk being: and: The corresponding non-isothermal mixing temperature for each process stream after each stage is calculated based on an energy balance [22].The maximum number of stages to represent a HEN during optimization is set to N s = max{N h , N c } [22].

Objective Function
The objective function is structured in such way that a flexible cost functions framework can get incorporated in order to consider different implementation related factors like piping and pumping.A different cost function for every coupling of possible hot and cold process stream matches can be defined.In addition, different cost functions can get added for each utility heat exchanger.If necessary they can be further refined for every stage in the superstructure representation of the HEN.The objective function is based on the proposed structure by Rathjens and Fieg [1] and was further developed for the present work: Hot utility and cold utility are used to reach the desired target outlet temperatures T out,n and cause costs based on their individual supply costs C HU and C CU .During optimization hot and cold utility can potentially be applied to each process stream depending on the maximum positive or negative deviation from the specified bounds of target temperatures.The respective utility which remains unused will make no contribution to the total annual costs (TAC).The binary variables z state the existence of a heat exchanger within the superstructure.The match-based costs X are calculated based on different cost function formulations with a wide range of possible correlations that can get implemented.The common cases of using a power function with the coefficients a 0 , a 1 and a 2 as well as a general polynomial representation up to 4th degree (coefficients b 0 , b 1 , b 2 , b 3 and b 4 ) are included.More specialized dependencies are covered using exponential (coefficients c 0 and c 1 ) or logarithmic expressions (coefficients d 0 and d 1 ): The different cost function formulations given in Equation (11) are chosen because they were able to represent the cost correlations for practical implementation that were faced when working together with our industrial partners.In order to be used in the objective function, the dependencies given in Equation (11) have to get expressed through the individual heat exchanger areas A, the temperature levels T and heat capacity flow rates .W of the involved streams and the heat load .

Q:
Energies 2019, 12, 784 The expression of cost dependencies through distinct variables was chosen due to ensuring the universal applicability of the shown approach without extending the programming work for each different problem.Furthermore these variables can be used with minimal computational overhead.Because ξ can be chosen out of an amount of alternatives given in Equation ( 12) the combinatorial possibilities of representing different cost dependencies are immense.The costs of utility heat exchangers X HU,n and X CU,n are structured equivalently to X ijk .

Genetic Algorithm
The genetic algorithm uses the common genetic operations like selection, crossover and mutation.As optimization variables the parameters (U A) ijk , .W h,ijk and .W c,ijk are used like proposed by Fieg et al. [22].The selection is based on the fitness value F. The fitness value represents the quality of an individual with respect to the objective function.It is based on the relative relation towards the average costs among all individuals and the minimum TAC of all individuals.The calculation of the fitness value is given in the following equation: The selection is carried out according to roulette wheel selection.Crossover is split in parameter crossover and structure crossover.Mutation is carried out on the variable parameters stated above.For details about the explicit formulations please refer to Fieg et al. [22].The probability for the application of a crossover operation is 89% in this work and the probability for parameter crossover was chosen as 23% according to Brandt [23].The probability for the general mutation is 1%, for parameter mutation it is 50% with a respective gene mutation probability of 1% [23].
Because of the optimization of the parameter (U A) ijk using the explicit temperature solution depicted in Equation ( 7) instead of the heat loads of the heat exchangers, every solution is thermodynamically feasible.Otherwise occurring temperature or heat load constraints can be omitted.The binary variable handling and corresponding constraints considering the continuous variables are adopted from Luo et al. [20].During optimization the strategy of excessive use of utilities is used [24].If the specified hot utility temperature is higher than any target temperature of all process streams and cold utility temperature is lower than any cold target temperature of all process streams, each network generated throughout optimization is feasible.As a result no outlet temperature constraints are considered here.Because .W h,ijk and .W c,ijk are also used as optimization variables, these values have to get constrained to be within a feasible region to hold the constraints given in Equations ( 14) and ( 15): During optimization .
W * h,ijk and .
W * c,ijk can occur that violate the above mentioned constraints.In order to keep these constraints, the normalization strategy from Fieg et al. is applied [22], which is shown by the following equations: . .
Furthermore, the structural control strategy proposed by Luo et al. [20] was adapted to ensure heterogeneity in the population and avoid local optimal solutions.
Structurally forbidden matches are removed automatically from a possible occurrence in the superstructure by setting the respective (U A) ijk to zero.A simple assignment of high costs for the specific match is not efficient against the background of fast convergence.This is particularly true for the second example shown below.

Local Optimization
Due to the strong nonlinearity of the stated objective function, an algorithm relying only on genetic operations would take too long to find an acceptable solution.Therefore, deterministic methods are used for local optimization.In order to generate promising HEN structures in the initialization step the approach of enhanced vertical heat transfer proposed by Stegner et al. [25] is used.The method proposed by Stegner et al. [25] adopts ideas known from the conventional vertical heat transfer concepts from traditional pinch approach but a new form of graphic depiction was implemented.Loop breaking as one of the heuristic approaches known from pinch technology is considered as well.This approach was successfully combined with stochastic optimization by Brandt et al. [26].
Furthermore, Newton's method is utilized for local parameter optimization [22].Equation (18) shows the nomenclature in Newton's method for the variables (U The partial derivatives needed in Equation ( 18) are calculated numerically using a step size of 10 −4 .

Examples and Results
In order to show the application of the presented approach two examples have been chosen.Example 1 is a mid-sized optimization problem for direct heat integration.Example 2 is an optimization problem for indirect heat integration using a heat recovery loop (HRL).

Example 1
Example 1 was first stated by Pho and Lapidus [27] and was named 10SP1.It consists of five hot and five cold streams.In the first part of the analysis in our work we focus on the cost functions already used in other publications to achieve comparable results (see Section 3.1.1).In the second part we focus on additional costs caused by the consideration of piping (see Section 3.1.2).

Example 1a
The 10SP1 problem is well studied and has been used as an optimization case several times so far with decreasing costs throughout the years by many researchers (e.g., Nishida et al. [28] or Flower and Linnhoff [29]. During the last two decades the 10SP1 problem has still attracted a considerable amount of attention.The associated problem data is given in Table 1.Two different formulations for the cost calculation of the annual costs of heat exchangers for the 10SP1 problem are given in the literature.The first cost formulation is given in Equation ( 19): The second cost formulation is given in Equation ( 20): A summary of the achieved TAC is given in Table 2. Several different approaches were used to tackle the 10SP1 problem.Lewin et al. [24] used a two-level approach consisting of a genetic algorithm and a lower level transformation into a linear parametric optimization problem.Lewin [30] also used an alternative nonlinear programming approach for the lower level optimization.Lin and Miller [31] used a tabu search algorithm to tackle the 10SP1 problem.Pariyani et al. [32] used a randomized algorithm for problems with stream splitting and a modified version of a previous work from Chakraborty and Ghosh [33] for designing networks without stream splitting.Yerramsetty and Murty [34] used a differential evolution algorithm and Peng and Cui [35] utilized a two-level simulated annealing algorithm.The latest work of Aguitoni et al. [36] found the best solution for the cost formulation 1 with a combination of a genetic algorithm and differential evolution. 1Solution without stream splits. 2 Solution with stream splits. 3Revised by Pariyani et al. [32].
The optimization procedure with the chosen approach incorporates 375 optimization variables considering five hot streams, five cold streams and five stages for the superstructure.Each heat exchanger has the same area-related cost function.
In our work cost formulation 1 was used for the optimization.The minimum TAC found was 42,963 $/yr.The solution is depicted in Figure 2 and shows a configuration of two stream splits as well as the use of two utility heat exchangers.simulated annealing algorithm.The latest work of Aguitoni et al. [36] found the best solution for the cost formulation 1 with a combination of a genetic algorithm and differential evolution.
The optimization procedure with the chosen approach incorporates 375 optimization variables considering five hot streams, five cold streams and five stages for the superstructure.Each heat exchanger has the same area-related cost function.
In our work cost formulation 1 was used for the optimization.The minimum TAC found was 42,963 $/yr.The solution is depicted in Figure 2 and shows a configuration of two stream splits as well as the use of two utility heat exchangers.The cost formulation for the optimization is favoring maximum energy recovery (MER) networks.Therefore, the obtained network recovers the maximum amount of energy and only cold utility is used.The TAC obtained for example 1a is significantly lower than in the publications cited in relation to the progress in TAC made over the years.The HEN shown in Figure 2 was also evaluated using cost formulation 2 and the resulting costs are 43,321 $/a, which is the lowest TAC found for formulation 2 so far.The results show that the chosen approach is capable of generating competitive results compared to other approaches presented in literature.Therefore, the The cost formulation for the optimization is favoring maximum energy recovery (MER) networks.Therefore, the obtained network recovers the maximum amount of energy and only cold utility is used.The TAC obtained for example 1a is significantly lower than in the publications cited in relation to the progress in TAC made over the years.The HEN shown in Figure 2 was also evaluated using cost formulation 2 and the resulting costs are 43,321 $/a, which is the lowest TAC found for formulation 2 so far.The results show that the chosen approach is capable of generating competitive results compared to other approaches presented in literature.Therefore, the consideration of even more complex problem formulations should be manageable, which is shown in the following parts.

Example 1b
Example 1b is built on the problem definition from example 1a.In addition, piping costs are considered.The required coordinates for the corresponding streams are taken from Pouransari and Maréchal [18].They added arbitrary Cartesian coordinates to the 10SP1 problem stated above.For the Energies 2019, 12, 784 9 of 18 hot utility coordinates the values from Rathjens and Fieg [1] are used.The additional data is given in Table 3.The required inner pipe diameter D i in m for each possible connecting pipe between a heat source and heat sink is calculated with the correlation for optimal diameter given by Peters et al. [37] for turbulent flow and D i ≥ 0.0254 m: (21) with the volumetric flow rate .
V in m 3 /s, the fluid density of 983 kg/m 3 and a specific heat capacity of 4.18 kJ/(kgK).
Peters et al. [37] give a diagram for the correlation between pipe diameter and the related costs per meter of pipe.Due to manually selecting and reading out values from this diagram there might be a slight deviation from the book values which cannot be further specified.Considering stainless-steel welded pipe of type 304 [37] results in the following cost correlations with respect to the stream heat capacity flow rates and heat exchanger areas: and: with the distance L being: The number of optimization variables is 375 as in example 1a.However, the number of potentially different cost functions is 25 (each match between hot and cold stream).The number of potentially different cost functions for the hot and cold utility usage is 10 for each utility.The real number of different cost functions is lower due to the fact that some distances between multiple hot and cold streams as well as utilities are equal.
The relatively short distances between the heat sources and heat sinks are compensated by the annualization of the piping costs over a one-year period.As a result the portion of piping costs becomes around 35% of TAC.
The solution with the lowest TAC found in this work is shown in Figure 3a.It has a TAC of 68,476 $/yr.A structurally different solution which was found during optimization with a less complex structure but the same utility configuration is shown in Figure 3b.
The relatively short distances between the heat sources and heat sinks are compensated by the annualization of the piping costs over a one-year period.As a result the portion of piping costs becomes around 35% of TAC.
The solution with the lowest TAC found in this work is shown in Figure 3a.It has a TAC of 68476 $/yr.A structurally different solution which was found during optimization with a less complex structure but the same utility configuration is shown in Figure 3b.As opposed to the solution of example 1a only two stages within the superstructure are taken by heat exchangers.Furthermore, the number of utility heat exchangers increased to three compared to the case not considering piping.The results still represent MER networks.Figure 4 and Figure 5 show the relative orientation of the heat sources and heat sinks towards each other.As opposed to the solution of example 1a only two stages within the superstructure are taken by heat exchangers.Furthermore, the number of utility heat exchangers increased to three compared to the case not considering piping.The results still represent MER networks.Figures 4 and 5 show the relative orientation of the heat sources and heat sinks towards each other.
annualization of the piping costs over a one-year period.As a result the portion of piping costs becomes around 35% of TAC.
The solution with the lowest TAC found in this work is shown in Figure 3a.It has a TAC of 68476 $/yr.A structurally different solution which was found during optimization with a less complex structure but the same utility configuration is shown in Figure 3b.As opposed to the solution of example 1a only two stages within the superstructure are taken by heat exchangers.Furthermore, the number of utility heat exchangers increased to three compared to the case not considering piping.The results still represent MER networks.Figure 4 and Figure 5 show the relative orientation of the heat sources and heat sinks towards each other.The configuration containing three utilities was dominant throughout the optimizations.The consideration of piping costs and thus the changed objective function had a strong influence on the heat load distribution within the HENs.The comparison between the optimal results of example 1b and example 1a shows about 14.2% (4016 $/yr) less piping costs, which is counterbalanced by an increase in installed heat exchanger area of about 16.2% (39.4 m 2 ) and thus an increase in capital costs for the heat exchangers of 14.7% (1307 $/yr).The configuration containing three utilities was dominant throughout the optimizations.The consideration of piping costs and thus the changed objective function had a strong influence on the heat load distribution within the HENs.The comparison between the optimal results of example 1b and example 1a shows about 14.2% (4016 $/yr) less piping costs, which is counterbalanced by an increase in installed heat exchanger area of about 16.2% (39.4 m 2 ) and thus an increase in capital costs for the heat exchangers of 14.7% (1307 $/yr).
The favorable design of using cold utility for the streams H3 and H5 is still preferred when considering piping costs during optimization due to the coordinates chosen by Pouransari and Maréchal [18].Furthermore the optimization is dominated by the close arrangement of the streams H3 and C4 as well as the streams H4, H5 and C3.Despite that, the solutions found considering piping costs already during the optimization are structurally different.Rathjens and Fieg [1] used a different cost functions formulation and were able to show significantly more local clustering of heat exchange similar to Figure 5b.

Example 2
The second example is a heat integration case with indirect heat integration utilizing a HRL with an intermediate fluid.It was taken from Chang et al. [13].It is a case study of a heat integration project in the southern part of China.The example comprises two different plants with a distance  of 1000 m.One plant is the heat source plant and the other one the heat sink plant with seven process streams each.The problem data is given in Table 4.The favorable design of using cold utility for the streams H3 and H5 is still preferred when considering piping costs during optimization due to the coordinates chosen by Pouransari and Maréchal [18].Furthermore the optimization is dominated by the close arrangement of the streams H3 and C4 as well as the streams H4, H5 and C3.Despite that, the solutions found considering piping costs already during the optimization are structurally different.Rathjens and Fieg [1] used a different cost functions formulation and were able to show significantly more local clustering of heat exchange similar to Figure 5b.

Example 2
The second example is a heat integration case with indirect heat integration utilizing a HRL with an intermediate fluid.It was taken from Chang et al. [13].It is a case study of a heat integration project in the southern part of China.The example comprises two different plants with a distance L of 1000 m.One plant is the heat source plant and the other one the heat sink plant with seven process streams each.The problem data is given in Table 4. C HU = 20 $/(kWyr), C CU = 8 $/(kWyr) [38]; C HX = 11, 000 + 150•A HX [39].
Energies 2019, 12, 784 12 of 18 In addition to the utility costs and the capital costs for the heat exchangers, the costs for pumping and piping are considered.The detailed equations are given in the following part.
For the calculation of the inner pipe diameter D i for the HRL, Equation ( 21) is used.The fluid density of 960 kg/m 3 and a specific heat capacity flow rate of 4.2 kJ/(kgK) is assumed.
The outer diameter D out in m, the specific pipe weight wt pipe in kg/m and the resulting pipe capital costs Pcul in $/m for schedule 80 steel pipes are calculated according to Stijepovic and Linke [39]: The resulting capital costs for piping in $ are: For the calculation of the pumping capital costs the fanning friction factor f , the Reynolds number Re and the fluid velocity u in m/s are used to estimate the pressure drop ∆p in Pa and the resulting costs [13,39].The fluid velocity is calculated according to Equation (29) using the fluids mass flow rate .m in kg/s of the intermediate fluid: The Reynolds number is calculated assuming a dynamic viscosity µ of 0.0002834 Pa s: The resulting pressure drop ∆p is given by: with the fanning friction factor of: The cumulated capital costs for the two pumps C pump in $ for the HRL are calculated according to Equation (33), utilizing the corresponding parameters for centrifugal pumps given in Jabbari et al. [40]: (33) The operating costs for pumping C pump,op in $/yr are composed of the pump efficiency η pump of 0.7, the price for electricity C el of 0.1 $/(kWh) [13], an operation duration t of 8000 h/yr and the parameters calculated beforehand [13]: Energies 2019, 12, 784 13 of 18 The annualization of the capital costs for the heat exchangers, piping and the required pumps is done using an annualization factor AF in 1/yr.The annualization factor is calculated according to Chang et al. [13] with an operation time n of five years and a fractional interest rate I of 10% per year: In order to use the superstructure of the utilized optimization model in this work, two pseudo-streams are added to the problem definition to emulate the HRL.The respective heat transfer coefficients for the two streams are 1 kW/(m 2 K) each [13,14].The resulting superstructure has eight stages as well as eight hot streams and eight cold streams.In order to represent the cost correlations described beforehand, they are transferred to a dependency in heat capacity flow rate to get integrated in the cost functions framework proposed in this work.Therefore, the costs for piping and pumping are added up and calculated for different mass flow rates.The resulting costs are plotted against heat capacity flow rates calculated based on the respective mass flow rates.A curve fitting is carried out afterwards.The deviation from calculating the exact costs based on the equations given above is below 0.0012% considering rounded values for the piping and pumping costs.The resulting correlations are shown in Equation ( 36): The utility heat exchangers for the two plants are considered as already installed like it was done by Chang et al. [13].Therefore, X HU,n and X CU,n are zero for this example.The first and last stages are used to include the costs for pumping and piping.The areas of the heat exchangers to emulate the HRL (ijk = {57, 505}) are chosen to be near infinite by the algorithm to avoid utility usage, thus create a self-adjusting temperature level and balance the heat transferred from the heat source plant to the heat sink plant.In order to prohibit direct heat integration between the two plants, all direct matches are forbidden.Altogether the applied superstructure yields 1536 optimization variables (compared to 1029 for direct heat integration without the pseudo streams).Due to the large number of forbidden matches (426) in relation to the number of possible matches (512), the number of optimization variables considered is reduced equivalently.The whole problem is described with only three different cost functions (see Equation (36)).
Iteration over the mass flow rate yields the optimal configuration for the HRL for example 2 which is shown in Figure 6.The superstructure representation of the solution in Figure 6 is shown in Figure A1 in Appendix A.
The TAC achieved for the solution shown in Figure 6 is 1.54 M$/yr.The associated annualized capital costs for the heat exchangers are 0.574 M$/yr (37.2% of TAC) and the utility consumption causes annual costs of 0.603 M$/yr (39.1% of TAC).The overall piping costs are 0.326 M$/yr (21.1% of TAC) and pumping costs are 0.038 M$/yr (2.5% of TAC).The solution obtained in this work is therefore 4.4% cheaper than the solution reported by Chang et al. [13].These differences are explained by the non-isothermal mixing model used and the relaxed temperature constraints in this work.

Conclusion
The presented approach to incorporate a flexible cost functions framework to synthesize cost-optimal heat exchanger networks (HENs) was carried out successfully and showed promising results.The flexible structured objective function allows for the integration of individual, match-dependent cost functions.The introduction of pseudo streams in combination with the flexible cost functions framework allow for the application for various problems.Corresponding optimizations utilizing different parametrizations have been carried out successfully.The universal applicability was shown by the execution of optimizations with different areas of application.The presented approach is applicable for direct as well as indirect heat integration utilizing the same superstructure and the same genetic algorithm for solving the problems.The use for combined direct and indirect heat integration is possible if initially forbidden matches get extended with cost information.This work presents results with lower TAC than other results published in literature beforehand.
For the practical implementation it is advisable to incorporate factors like piping already during optimization.The consideration of these factors can have a huge impact of the overall HEN structure 6. Plant representation of the optimal HEN configuration found of example 2 including a HRL (TAC: 1.54 M$/yr).

Conclusions
The presented approach to incorporate a flexible cost functions framework to synthesize cost-optimal heat exchanger networks (HENs) was carried out successfully and showed promising results.The flexible structured objective function allows for the integration of individual, match-dependent cost functions.The introduction of pseudo streams in combination with the flexible cost functions framework allow for the application for various problems.Corresponding optimizations utilizing different parametrizations have been carried out successfully.The universal applicability was shown by the execution of optimizations with different areas of application.The presented approach is applicable for direct as well as indirect heat integration utilizing the same superstructure and the same genetic algorithm for solving the problems.The use for combined direct and indirect heat integration is possible if initially forbidden matches get extended with cost information.This work presents results with lower TAC than other results published in literature beforehand.
For the practical implementation it is advisable to incorporate factors like piping already during optimization.The consideration of these factors can have a huge impact of the overall HEN structure and can lead towards significantly different and more efficient solutions.Local clustering was observed for some solutions like already reported by Rathjens and Fieg [1].

Acknowledgments:
The authors gratefully acknowledge the support of our industrial partner, the "weyer group".

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 1 .
Figure 1.Stage-wise superstructure of a heat exchanger network (HEN) with two hot and three cold process streams.

Figure 1 .
Figure 1.Stage-wise superstructure of a heat exchanger network (HEN) with two hot and three cold process streams.

Figure 4 .
Figure 4. (a) 3D layout representation of the optimal HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load; and (b) 3D layout representation of the alternative HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load.

Figure 4 .
Figure 4. (a) 3D layout representation of the optimal HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load; and (b) 3D layout representation of the alternative HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load.

Figure 4 .
Figure 4. (a) 3D layout representation of the optimal HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load; and (b) 3D layout representation of the alternative HEN configuration of example 1b, the chosen line thickness is proportional to the corresponding heat exchanger heat load.

Figure 5 .
Figure 5. (a) 2D layout representation of the optimal HEN configuration of example 1b; (b) 2D layout representation of the alternative HEN configuration of example 1b.

Figure 5 .
Figure 5. (a) 2D layout representation of the optimal HEN configuration of example 1b; (b) 2D layout representation of the alternative HEN configuration of example 1b.

Figure 6 .
Figure 6.Plant representation of the optimal HEN configuration found of example 2 including a HRL (TAC: 1.54 M$/yr).

Funding:
Authors would like to acknowledge the financial support from German Federal Ministry for Economic Affairs and Energy through ZIM program (Zentrales Innovationsprogramm Mittelstand, project number: ZF4025905CL7).The publication was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)-project number: 392323616 and the Hamburg University of Technology (TUHH) in the funding program "Open Access Publishing".

Table 1 .
Problem data for example 1a.

Table 2 .
Results comparison for example 1a.

Table 2 .
Results comparison for example 1a.

Table 3 .
Additional Cartesian coordinates for the 10SP1 problem.

Table 4 .
Problem data for example 2.

Table 4 .
Problem data for example 2.