Document Type: Full Lenght Research Article
Authors
^{1} Faculty of Mechanical Engineering, Semnan University, P.O. Box 3513119111, Semnan, Iran
^{2} Faculty of Mechanical Engineering, Semnan University, P.O. Box 3513119111, Semnan, Iran.
Abstract
Keywords
Main Subjects
A B S T R A C T Many studies have been performed by researchers about Shell and Tube Heat Exchanger but the MultiObjective Big BangBig Crunch Algorithm (MOBBA) technique has never been used in such studies. This study presents application of thermaleconomic multiobjective optimization of shell and tube heat exchanger using MOBBA. For an optimal design of a shell and tube heat exchanger, it was first modelled thermally using eNTU method while BellDelaware procedure was applied to estimate its shell side heat transfer coefficient and pressure drop. The MOBBA method was applied to obtain the maximum effectiveness (heat recovery) and the minimum total cost as two objective functions. The results of the optimal designs are a set of multiple optimum solutions, called ‘Pareto optimal solutions'. In order to show the accuracy of the algorithm, a comparison is made with the Nondominated Sorting Genetic Algorithm (NSGAII) and MOBBA which are developed for the same problem. © 2016 Published by Semnan University Press. All rights reserved. 
History: Submitted 14 January 2014 Revised 24 July 2014 Accepted 19 June 2016
Keywords: Shell and Tube Heat Exchanger MultiObjective Big BangBig Crunch Algorithm (MOBBA) NonDominated Sorting Genetic Algorithm (NSGAII) Effectiveness Total cost.

Thermaleconomic optimization of shell and tube heat exchanger using a new multiobjective optimization method
Mohammad Sadegh Valipour^{1}, Mojtaba Biglari^{1}, Ehsanolah Assareh^{2}^{, *}
^{1}Faculty of Mechanical Engineering, Semnan University, P.O. Box 3513119111, Semnan, Iran ^{2,*}Department of Mechanical Engineering, Dezful Branch, Islamic Azad University, Dezful,Iran

Journal of Heat and Mass Transfer Research
Journal homepage: http://jhmtr.journals.semnan.ac.ir 


1. Introduction
Shell and tube heat exchangers are probably the most common type of heat exchangers applicable for a wide range of operating temperatures and pressures. Shell and tube heat exchangers are widely used in refrigerating, power generating, heating and air conditioning, chemical processes, manufacturing and medical applications.
Address of correspondence author: Ehsanolah Assareh, Department of Mechanical Engineering, Dezful Branch, Islamic Azad University, Dezful,Iran. Email: E.assareh@semnan.ac.ir

Designing the heat exchangers is a complex procedure and it requires a good knowledge of thermodynamics, fluid dynamics, cost estimation and optimization. The objectives involved in the design optimization of heat exchangers are thermodynamic (i.e. maximum efficiency) and economic (i.e. minimum cost). The conventional design approach for heat exchangers involves rating a large number of different exchanger geometries to identify those designs that satisfy a given heat duty and a set of geometric and operational constraints. This approach is timeconsuming, and does not guarantee an optimal solution. [3]
2. Literature Review
Several studies are presented to propose optimization of shell and tube heat exchanger.Min Zhao and Yanzhong Li used an effective layer pattern optimization model for multistream platefin heat exchanger using genetic algorithm. In this study, an effective layer pattern optimization model using Genetic Algorithm (GA) was developed in detail [4]. Suxin Qian et al. presented applicability of entransy dissipation based thermal resistance for the design optimization of twophase heat exchangers. In this study, the evaluation of twophase entransy was achieved by optimizing one tubefin heat exchanger and one micro channel heat exchanger based on a validated heat exchanger modeling tool [5]. Khaled Saleh et al. applied approximation assisted optimization of headers for a new generation of aircooled heat exchangers. In this study an online multiobjective approximation assisted optimization approach was used to design optimum headers for compact aircooled heat exchangers [6]. Amin Hadidi and Ali Nazari used the design and economic optimization of shellandtube heat exchangers using BiogeographyBased Optimization (BBO) algorithm. In this research, a new shell and tube heat exchanger optimization design approach was developed based on biogeographybased optimization (BBO) algorithm. The BBO algorithm has some good features in reaching the global minimum in comparison with other evolutionary algorithms [7]. SalimFettaka et al. developed the design of shellandtube heat exchangers using multiobjective optimization. In this study, a multiobjective optimization of the heat transfer area and pumping power of a shellandtube heat exchanger was presented to provide the designer with multiple Paretooptimal solutions which captured the tradeoff between the two objectives [8]. Sreepathi and G.P. Rangaiah applied an improved heat exchanger network retrofitting using exchanger reassignment strategies and multiobjective optimization. In this study, several ERS (exchanger reassignment strategies) for HEN retrofitting were proposed and tested by performing single objective optimization on retrofit problems using IDE (integrated differential evolution) [9]. Viviani C. Onishiet al. developed a mathematical programming model for heat exchanger design through optimization of partial objectives. Mathematical programming can be used for the optimal design of ShellandTube Heat Exchangers (STHEs). This study proposed a Mixed Integer NonLinear Programming (MINLP) model for the design of STHEs, following rigorously the standards of the Tubular Exchanger Manufacturers Association (TEMA). Bell–Delaware Method was used for the shellside calculations [10]. Jiang feng Guo et al. used multiobjective optimization of heat exchanger based on entransy dissipation theory in an irreversible Brayton cycle system. In this study a multiobjective optimization of the main heat exchanger in a regenerative Brayton cycle system was carried out based on entransy dissipation. The best tradeoff between the entransy dissipation numbers caused by heat transfer and fluid friction is achieved in the Pareto optimal solutions and a decrease in the entransy dissipation related to heat transfer inevitably leads to an increase in the entransy dissipation due to fluid friction, and vice versa. The entransy dissipation plays a decisive role in the network output due to heat transfer rather than due to fluid friction [11]. R. Venkata Rao and Vivek Patel et al. presented multiobjective optimization of heat exchangers using a modified teachinglearningbased optimization algorithm. In this study, a modified version of the TLBO algorithm was introduced and applied for the multiobjective optimization of heat exchangers. The platefin heat exchanger and the shell and tube heat exchanger were considered for the optimization [12]. Ming Pan et al. used optimization for the retrofit of large scale heat exchanger networks with different intensified heat transfer techniques. In this study, an optimization method was developed for dealing with the retrofit of large scale HENs in which the location of intensified heat transfer within the network and its degree of intensification were identified systematically, given the objective function and design constraints, including the topological limitation in the existing heat recovery systems[13]. Juan I. Manassaldi et al. presented optimization mathematical model for the detailed design of air cooled heat exchangers. This study presents a disjunctive mathematical model for the optimal design of air cooled heat exchangers. The model involves seven discrete decisions which are related to the selection of the type of the finned tube, the number of tube rows, the number of tube per row,the number of passes, the fins per unit length, the mean fin thickness and the type of the flow regime. Each discrete decision is modeled using disjunctions, boolean variables and logical propositions. The main continuous decisions are fan diameter, bundle width, tube length, pressure drops and velocities in both sides of the ACHE, heat transfer area and fan power consumption[14].Vivek Patel and VimalSavsani presented optimization of a platefin heat exchanger design through an improved MultiObjective TeachingLearning Based Optimization (MOITLBO) algorithm. In this study, the MultiObjective Improved TeachingLearningBased Optimization (MOITLBO) algorithm was introduced and applied for the multiobjective optimization of platefin heat exchangers[15]. Costa and Queiroz developed design optimization of shellandtube heat exchangers [16]. Caputo et al. presented heat exchanger design based on the economic optimization [17]. Fesanghary et al. applied a harmony search algorithm to design optimization of shell and tube heat exchangers [18]. Hilbert et al. developed parallel genetic algorithms to shape optimization of a heat exchanger [19]. Sanaye and Hajabdollahi used multiobjective optimization of shell and tube heat exchangers [20]. PonceOrtega et al. used the genetic algorithms for the optimal design of shellandtube heat exchangers [21]. Jie Yang et al. developed optimization of shellandtube heat exchangers using a general design approach motivated by constructal theory [22]. Daniël Walraven et al. used optimum configuration of shellandtube heat exchangers in lowtemperature organic Rankine cycles [23]. Jie Yang et al. developed optimization of shellandtube heat exchangers conforming to TEMA standards with the designs motivated by constructal theory [24]. Mohsen Amini and Majid Bazargan used twoobjective optimization in shellandtube heat exchangers using genetic algorithm [25]. Literature review also indicates that BBA algorithm has never been used for such a study.
2. Mathematical Model for Optimization of Shell and Tube Heat Exchanger Mathematical model
Based on the work of Sanaye and Hajabdollahi [20], effectiveness of the standards of the Tubular Exchanger Manufactures Association (TEMA) Etype Shell and Tube Heat Exchanger (STHE) is given as the following[20]:
(1)

Where the number of transfer units (NTU) and the heat capacity ratio (C) are defined as follows [20]:
(2) 

(3) 
Where A_{t} is the total tube outside heat transfer surface area and U_{o} is the overall heat transfer coefficient which is computed by the following equation [20]:
(4) 
Where L, N _{t}, d_{i}, d_{o}, R_{i,f}, R _{o,f}, and k_{w }are the tube length, tube number, tube inside and outside diameter, tube and shell side fouling resistances and thermal conductivity of tube wall respectively.
(5) 
The tube side heat transfer coefficient is estimated by the following equation [20]:
(6) 
Where k_{t} and Pr_{t} are the tube side fluid thermal conductivity and Prandtl number respectively. Also Re_{t} is the tube flow Reynolds number which is defined as follows [20]:
(7) 
Where m_{t} is the mass flow rate and A_{o,t} is the tube side flow cross section area per pass estimated as [20]:
(8)

And n_{p }is the number of tube passes. Furthermore, the tube side pressure drop is also estimated by the following formula:
(9)

Where ΔP_{t }includes the pressure drop due to the flow contraction, acceleration, friction, and expansion, four terms in Eq. (9). K_{c} and K_{e }are the tube entrance and exit pressure loss coefficients. Furthermore f_{t} is the tube side friction factor estimated as the following [20]:
(10)

The shell diameter is estimated by the following equation:
(11)

The shell side Reynolds number is as follows:
(12)

Where p_{t} is the tube pitch and CL is the tube layout constant that has a unit value for 45 and 90 tube arrangements and 0.87 for 30 and 60 tube arrangements. Also, CTP is the tube count constant which is 0.93, 0.9, and 0.85 for the single pass, two passes and three passes of tubes, respectively.
BellDelaware method is used in this study to compute the shell side heat transfer coefficient and pressure drop. For more information, readers are referred to ref. [20].
3. Entropy Generation Number
The irreversibility losses in the heat exchanger are evaluated in terms of entropy generation. The entropy generation rate caused by the finite temperature difference can be written as follows [26]:
(13)

For an incompressible fluid under nonadiabatic condition the entropy generation rate caused by fluid friction is expressed as follows [26]:
(14)

The total entropy generation rate in heat exchanger can be written as follows [26] :
(15)

When the heat capacity rate of the hot fluid is larger than that of the cold fluid, the outlet temperature of both fluids can be calculated as follows [20]:
(16)


(17)

4. MultiObjective Optimization Problems
Many optimization problems can be presented by the following general mathematical model [27]:
Max. /Min. F (x) xϵ X 
(18) 
F: X → R 
(19) 
In which f(x) is the objective function, x is the decision variables vector, X is the decision variables space, and R is the set of real numbers.
Figure (1a) shows a schematic of transferring data from a two dimensional decision variable space to a decision space in single objective problems. In multiobjective optimization problems, optimizing a vector of objectives is considered instead of satisfying a single objective.
Generally, in mathematical form, multiobjective optimization problems can be deﬁned as follows [27]:
Max. /Min f(x) = (f_{1}(x), f_{2}(x)… f_{m}(x)) xϵ X 
(20) 
F: X → R^{m}, mϵ M 
(21) 
Where m denotes the set of natural numbers and N denotes the number of objectives. There are techniques such as the weighting method and Ɛconstraint method which transfer multiobjective problems to a singleobjective one using different combinations of a weighting vector and constraints. Thus, each optimal solution can be assigned to a speciﬁc combination of weighting vector and constraint. Hence, in each run of the algorithm, a single point (solution) can be achieved. However, multiobjective evolutionary algorithms are capable of ﬁnding almost all candidate solutions (Pareto) in a single run.
Figure (1b) presents a schematic of transferring data from two dimensional decision variables space to the decision space in a twoobjective problem.
As it is shown, each set of decision variables has been related to a couple of objectives. Note that none of the solutions dominates the others. In other words, if all the objective values of a solution dominate the corresponding values of another solution, the former will be a dominated solution and the latter will be removed. Otherwise, both solutions will be located in the nondominated set.
The dominated and nondominated relations between objectives values in a biobjective problem are shown in Fig.2.Both objectives are minimized. In this ﬁgure, the solutions labeled by 1 or 2 have nondominated conditions individually. We should note that the set labeled 1 dominates the set labeled 2. In the optimization procedure, the best set of nondominated solutions is called Paretofront. Thus, there are two Pareto in Figure (2), and the one which is labeled 1 is the Paretofront. [27]
(a) 
(b) 
Fig. 1 The schematic of decision and decision variable spaces in: (a) singleobjective and (b) multiobjective problems. 
Fig.2 The schematic of dominated and nondominated conditions of solutions in a projective problem. 
The analogy of Big Bang  Big Crunch Algorithm (BBBCA) with evolutionary algorithms makes it evident that using a Pareto ranking scheme [28] could be the straightforward way to extend the approach to handle the multiobjective optimization problems. The historical record of the best solutions found by a particle (i.e., an individual) could be used to store nondominated solutions generated in the past (this would be similar to the notion of elitism used in evolutionary multiobjective optimization). The use of global attraction mechanisms combined with a historical archive of previously found nondominated vectors would motivate convergence toward globally nondominated solutions.
5. Big BangBig Crunch Algorithm (BBA)
In the BBBC algorithm proposed by Erol and Eksin[29], the initial bigbang is identical to the first step of the other evolutionary methods in that an initial population of candidate solutions is randomly generated over the entire search space. Erol and Eksin compared this random nature of the Big Bang to the energy dissipation or the transformation from an ordered state (a convergent solution) to a chaotic state (generation of new set of candidate solutions). In the Big Crunch phase following the Big Bang, a contraction operation is applied to randomly distributed candidate solutions. The contraction operator takes the current positions (represented by the values of the design variables) of each candidate solution in the population and its corresponding penalized fitness function value to compute a center of mass. The center of mass is the weighted average of the candidate solution positions where the weight associated with the position of each candidate solution is the inverse of the corresponding penalized fitness function. The averaging is done with respect to the inverse of the penalized fitness function values. The center of mass is computed as the following [29]:
(22)

Where is the position of candidate k in an ndimensional searchspace, Fk is the penalized fitness function value of candidate k, and NC is the candidate population size. New positions of the candidate solutions for the next iteration of the Big Bang are obtained using the following equation [29]:
(23) 
Where is the standard deviation of the normal distribution. In the BBBC algorithm, the standard deviation is related to asubset of search space and is obtained by the following equation [29]:
(24) 
Where r is a random number from a standard normal distribution, is a parameter limiting the size of the search space, and are the upper and lower limits of the values of the design variables respectively, and is the number of Big Bang iterations.
In order to improve the computational efficiency of BBBC algorithm, Kaveh and Talatahari use the social behaviour of bird flocking and fish schooling model in particle swarm optimization. The swarm’s movement is directed by both their own experience and the population’s experience. For every iteration, a particle moves towards a direction computed from the local best solution and the global best solution. This concept is used in this study where the BBBC algorithm not only utilizes the center of mass but also employs the global best solution to generate the new solution.
A modified version of Eq. (23) is given as follows [29]:
(25)

Where , are thecontrolling parameters, is the local best solution,and is the global best solution. Since normally distributednumbers can exceed ±1, it is necessary to limit candidate positions to the prescribed search space boundaries. As a result of this contraction, there is an accumulation of candidate solutions at the search space boundaries.
In Eq. (25) , and are the adjustable parameters that control the influence of the local and the global best solutions on the new positions of candidate solutions.
In addition, the BBBC algorithm employs a multiphase search. In a twophase search, the BBBC algorithm is initially applied to the entire search space. After the convergence of Phase I, the Phase II search starts in a reduced search around from Phase I. The search space is reduced to 20% of the original search space [30].
6. EconomicModelling
The total cost of the STHE heat exchanger is given by the following [20]:
(26)

Where C_{in} and C_{op} are the total investment cost and the operating cost of the heat exchanger respectively and are defined as follows [20]:
(27)


(28) 
Where n_{y} is the equipment life in year, i is the annual discount rate and k is the depreciation time in year. Also, ΔP and m are the pressure drop and the mass flow rate of the fluid. Cop is the annual operating cost and it is calculated as follows:
(29)

Where k_{el}, τ and ɳ_{p} are the electricity unit cost, the operation hours of the heat exchanger per year and the pump efficiency respectively. The detailed calculation of different parameters of the above mentioned equations is given in Refs. [20, 26].
7. Case Study
The analysis of the case study has been performed by Sanaye et al. [20] using NSGAII approach according to the literature [2]. The original design specifications, shown in Table (1), are fed as inputs into the MOBBA algorithm.
In Figure (3), the optimal heat exchanger architectures obtained by MOBBA are compared with those of Sanaye et al. [20] using the NSGAII approach and also with the original design solution given by Shah et al. [2].
The application example of STHE was also analysed to ensure the capability of MOBBA which was earlier analysed using NSGAII by Sanaye and Hajabdollahi [20]. With the help of fresh water, the oil cooler STHE used to lower the temperature of the oil was designed and optimized for providing the maximum effectiveness and the minimum total cost. The high temperature of oil with mass flow rate of 8.1 kg/s was poured into the shell side of STHE at a temperature of 78.3 ºC. Also, the fresh water with mass flow rate of 12.5 kg/s at a temperature of 30 ºC was poured into the tube side of STHE. The thermo physical properties of both the fluids were taken from Ref. [20].
The following inequality constraints which are bound by lower and upper limits of the design variables are considered in the present study and three different tube arrangements (30º, 45º and 90º) are also considered for the optimization.
0.0112 ≤ di ≤ 0.0153 
1.25 ≤ pt/do ≤ 2

3 ≤ L ≤ 8

100 ≤ Nt ≤ 600

0.19 ≤ bc/Ds ≤ 0.32

0.2 ≤ bs/Ds ≤ 1.4

3 ≤ L/Ds ≤ 12

In this study, the following assumptions were made:
Table1. The process input and physical properties for different case studies" 

Properties 
Pr 
K(W/m^{2}.K) 
Cp(j/kg.K) 
µ(Pa.s) 
R_{f} 
Shell side 
966 
0.14 
2115 
0.0643 
0.00015 
Tube side 
4.77 
0.634 
4120 
0.000695 
0.000074 
Fig.3 Theoptimized point's distribution in pareto by using MOBBA 
8. Distance Criterion
The distance criterion measures (d) cover and give distribution to the proposed space. In this study, the distance criterion is used to compare the results of two algorithms, MOBBA and NSGAII. The result is discussed in section 9.
(30)


(31) 
Also, using di , we can find the minimum point distance of the diagram based on the ideal point using such criterion [20].
9. Optimization Results
The numbers of iterations for finding maxima or minima, in the search space are almost 104 and the number of populations and repositories are 100. The results are shown in figure (4) via the Pareto diagram. The figure demonstrates clearly the presented differences between both objective functions.
Increasing heat transfer in the heat exchanger, the total cost would increase. So, we use multi objective optimization method for shelltube heat exchanger to decrease the cost. Figure (6) shows the minimum amounts of total cost with the least effectiveness. Here, the MOBBA algorithm results are compared with the results of the NSGA II algorithm previously proposed by Mr. Sanaye in 2010 [20].
We use distance criterion for comparing the results of two algorithms, MOBBA and NSGAII, which this scale is 0.3919 for the NSGAII algorithm and 0.3112 for the MOBBA algorithm. These amounts indicate that the point achieved by the MOBBA algorithm is closer to the ideal point rather than the point achieved by the NSGAII algorithm. Figure (4) shows a better performance of the MOBBA algorithm compared to the NSAG II algorithm.
For both algorithms in the figure (3), A and C points are happened when optimization is done as a single objective and so our objective functions are maximization of the effectiveness. For B and D points, our objective functions are minimization of the total cost. Figures (5) and (6) are related to the maximum effectiveness and minimum total cost, respectively. It is obvious that for the total cost, the results in the MOBBA algorithm are optimized to 7.36 % and so for Q maximum, the MOBBA algorithm has better results than the NSGAII.
10. External Temperatures
According to the equations (16) and (17) the relations for external temperatures of heat exchanger and the amounts of these temperatures are obtained based on the limitation of the total cost, effectiveness and decision variables in figure (3) for both algorithms in the figures (8) to (10).
These figures emphasize the problem that there is an external temperature in the tube for an external temperature in the shell with the optimized amounts of decision variables. The maximum amounts of effectiveness in (16) and (17) equations are accompanied with the minimum amounts for external temperature of the shell, and the maximum external temperature of the tube.
Also, the optimized amount of the external temperature in the tube is 63.85 for the MOBBA algorithm and for the external temperature of the shell is 44.440C
11. Entropy Balance
Figure (11) and figure (12) show the entropy amounts based on the effectiveness and also show the total cost for highlighted points in figure (3).
Figure (11) shows the entropy amount based on the output heat. More specifically, it shows that the amount of system entropy increases with increasing the output.
Also, figure (12) shows incurred cost increased with increasing the output. For the amount of the best optimized point, which is obtained by pareto diagram in figure (3) for the MOBBA algorithm, the amount of effectiveness is 0.7051 and the cost is 28184 $ and the amount of produce entropy cost is 29287.
12. Conclusion
Designing the heat exchanger can be a complicated task. Thus, the advanced optimization methods are useful for identifying the best and cheapest heat exchanger for a specific duty. In this study, the shell and tube heat exchanger design and the multiobjective optimization are performed with the MOBBA algorithm. The results are compared with those of Sanaye paper [20]. The MOBBA algorithm is a useful method for the shell and tube heat exchanger design and optimization. These results are demonstrated in similar conditions for twoobjective functions (effectiveness and total cost) and indicate that the MOBBA is a better and more effective method than the NSGAII algorithm. The algorithm proposed here can help the manufacturer and engineers to optimize heat exchangers in engineering applications.
Fig. 4 The related results to D criterion for the MOBBA & NSGAII.

Fig. 5 The comparison of the maximum effectiveness for the MOBBA& NSGAII[8].

Fig. 6 The comparison of the least total cost for MOBBA& NSGAII [8].

Fig. 7 The comparison of the maximum heat transfer for the MOBBA& NSGAII [8]. 
Fig. 8[S1] The external temperature based on effectiveness for the MOBBA& NSGAII [8].

Fig. 9 The comparison of the maximum tube external temperature for the MOBBA& NSGAII [8].

Fig. 10 The comparison of the least shell external temperature for the MOBBA & NSGAII [8].

Fig. 11 The entropy balance based on effectiveness efficiency for the MOBBA.

Fig. 12 The entropy balance based on total cost for the MOBBA.

Nomenclature 

A 
heat transfer area 
(m^{2}) 
BC 
baffle cut 
(m) 
C_{p} 
specific heat in constant pressure 
(J/kg K) 
C 
heat capacity rate ratio 
(C_{min}/C_{max}) 
C_{in} 
total investment cost 
($) 
C_{op} 
total operating cost 
($) 
C_{o} 
annual operating cost 
($/yr) 
C_{total} 
total cost 
($) 
C_{L} 
tube layout constant 
() 
C_{TP} 
tube count calculation constant 
() 
d_{i} 
tube side inside diameter 
(m) 
d_{o} 
tube side outside diameter 
(m) 
D_{s} 
shell diameter 
(m) 
f 
friction factor 
() 
G_{s} 
fluid mass velocity based on the minimum free area 
(kg/s.m^{2}) 
h_{i} 
tube side heat transfer coefficient 
(W/m^{2} K) 
h_{o} 
shell side heat transfer coefficient 
(W/m^{2} K) 
i 
annual discount rate 
(%) 
j_{c} 
correction factor for baffle configuration 

j_{s} 
correction factor for bigger baffle spacing at the shell inlet and outlet sections 

j_{r} 
correction factor for the adverse temperature gradient in laminar flows 

j_{b} 
correction factor for bundle and pass partition bypass streams 

j_{l} 
correction factor for baffle leakage effect 

K_{c} 
entrance pressure loss coefficient 
() 
K_{e} 
exit pressure loss coefficient 
() 
k_{el} 
price of electrical energy 
($/kWh) 
k 
thermal conductivity 
(W/m k) 
L 
tube length 
(m) 
L_{bc} 
baffle spacing 
(m) 
m 
mass flow rate 
(kg/s) 
n_{p} 
number of tube pass 
()

N_{t} 
number of tube 
() 
NTU 
number of transfer units 
() 
P_{t} 
tube pitch 
(m) 
Q 
heat transfer rate 
(kW) 
Re 
Reynolds number 
() 
T 
temperature 
(^{o}C) 
U 
overall heat transfer coefficient 
(W/m^{2} K) 
σ 
ratio of minimum free flow area to frontal area 
() 
P 
pumping water 
(W) 
ratio of minimum free flow area to frontal area 
() 

ρ 
density

(kg/m^{3}) 
Subscripts

i= input O= output 

References
[1] G.F.Hewitt, editor. Heat exchanger design handbook. New York: Begell House; (1998).
[2] R.K. Shah, Bell KJ. The CRC handbook of thermal engineering. CRC Press; (2000).
[3] V. Rao, V. Patel, Multiobjective optimization of heat exchangers using a modified teachinglearningbased optimization algorithm, Applied Mathematical Modeling, 37(3), 11471162, (2013).
[4] Min Zhao , Yanzhong Li,An effective layer pattern optimization model for multistream platefin heat exchanger using genetic algorithm, International Journal of Heat and Mass Transfer, 60, 480–489, (2013).
[5] Suxin Qian, Long Huang, Vikrant Aute, Yunho Hwang, ReinhardRadermacher, Applicability of entransy dissipation based thermal resistance for design optimization of twophase heat exchangers. Applied Thermal Engineering, 55, (1), 140–148, (2013).
[6] KhaledSaleh, Omar Abdelaziz, Vikrant Aute, ReinhardRadermacher, ShapourAzarm, Approximation assisted optimization of headers for new generation of aircooled heat exchangers. Applied Thermal Engineering, 61(2), 817–824, (2013).
[7] Amin Hadidi, Ali Nazari, Design and economic optimization of shellandtube heat exchangers using biogeographybased (BBO) algorithm. Applied Thermal Engineering, 51(1), 1263–1272, (2013).
[8] Salim Fettaka, Jules Thibault, Yash Gupta. Design of shellandtube heat exchangers using multiobjective optimization. International Journal of Heat and Mass Transfer, 60, 343–354, (2013).
[9] Bhargava Krishna Sreepathi, G.P. Rangaiah, Improved heat exchanger network retrofitting using exchanger reassignment strategies and multiobjective optimization. Energy Available online 24 February 2014 In Press, Corrected Proof.
[10] Viviani C. Onishi, Mauro A.S.S. Ravagnani, José A. Caballero, Mathematical programming model for heat exchanger design through optimization of partial objectives, Energy Conversion and Management, 74, 60–69, (2013)
[11] JiangfengGuo, XiulanHuai , Xunfeng Li, Jun Cai, Yongwei Wang. Multiobjective optimization of heat exchanger based on entransy dissipation theory in an irreversible Brayton cycle system. Energy, 63, 95–102, (2013).
[12] R. VenkataRao ,Vivek Patel, Multiobjective optimization of heat exchangers using a modified teachinglearningbased optimization algorithm. Applied Mathematical Modelling, 37(3), 1147–1162, (2013).
[13] Ming Pan, Igor Bulatov, Robin Smith, JinKuk Kim, Optimisation for the retrofit of large scale heat exchanger networks with different intensified heat transfer techniques. Applied Thermal Engineering, 53( 2) , 373–386, (2013).
[14] Juan I. Manassaldi , Nicolás J. Scenna, Sergio F. Mussati, Optimization mathematical model for the detailed design of air cooled heat exchangers. Energy, 64, 734–746, (2014)
[15] Vivek Patel, Vimal Savsani, Optimization of a platefin heat exchanger design through an improved multiobjective teachinglearning based optimization (MOITLBO) algorithm. Chemical Engineering Research and Design, Available online 11 February 2014 In Press, Accepted Manuscript.
[16] L.H. Costa, M. Queiroz, Design optimization of shellandtube heat exchangers, Applied Thermal Engineering, 28, 17981805,(2008).
[17] A. Caputo, P. Pelagagge, P. Salini, Heat exchanger design based on economic optimization, Applied Thermal Engineering, 28(10), 1151–1159, (2008).
[18] M. Fesanghary, E. Damangir, I. Soleimani, Design optimization of shell and tube heat exchangers using global sensitivity analysis and harmony search algorithm, Applied Thermal Engineering, 29, 1026–1031, (2009).
[19] R. Hilbert, G. Janiga, R. Baron, D.Venin, Multiobjective shape optimization of a heat exchanger using parallel genetic algorithms, International Journal of Heat and Mass Transfer,49(15–16), 2567–2577, (2006).
[20] S. Sanaye, H. Hajabdollahi, Multiobjective optimization of shell and tube heat exchangers, Applied Thermal Engineering, 30(14–15), 1937–1945,(2010).
[21] J. PonceOrtega, M. SernaGonzález, A. JiménezGutiérrez, Use of genetic algorithms for the optimal design of shellandtube heat exchangers, Applied Thermal Engineering, 29(2–3), 203–209,(2009).
[22] Jie Yang, SunRyung Oh, Wei Liu,Optimization of shellandtube heat exchangers using a general design approach motivated by constructal theory.International Journal of Heat and Mass Transfer,Volume 77, October 2014, Pages 1144–1154.
[23] DaniëlWalraven, Ben Laenen, William D’haeseleer,Optimum configuration of shellandtube heat exchangers for the use in lowtemperature organic Rankine cycles.Energy Conversion and Management. Volume 83, July 2014, Pages 177–187.
[24] Jie Yang, Aiwu Fan, Wei Liu, Anthony M. Jacobi.Optimization of shellandtube heat exchangers conforming to TEMA standards with designs motivated by constructal theory.Energy Conversion and Management.Volume 78, February 2014, Pages 468–476.
[25] Mohsen Amini, MajidBazargan.Two objective optimization in shellandtube heat exchangers using genetic algorithm.Applied Thermal Engineering.Volume 69, Issues 1–2, August 2014, Pages 278–285.
[26] R.K. Shah, P.Sekulic, Fundamental of Heat Exchanger Design, John Wiley & Sons, Inc, (2003).
[27] E. Mehdipour, O. Haddad, M. Tabari, M. Mari. Extraction of decision alternatives in construction management projects: Application and adaptation of NSGAII and MOPSO, Expert Systems with Applications, 39(3), (2012).
[28] D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning, Reading, MA: AddisonWesley, (1989).
[29] Osman K. Erol, Ibrahim Eksin,A new optimization method: Big Bang–Big Crunch, Advances in Engineering Software, 37( 2), February 2006, Pages 106–111.
[30] H. Tang , J. Zhou , S. Xue , L. Xie, Big BangBig Crunch optimization for parameter estimation in structural systems, Mechanical Systems and Signal Processing, 24, 2888–2897, (2010).