## Abstract

A heat exchanger is a unit operation used to transfer heat between two or more fluids at different temperatures. There are many different types of heat exchangers that are categorized based on different criteria, such as construction, flow arrangement, heat transfer mechanism, etc. Heat exchangers are optimized based on their applications. The most common criteria for optimization of heat exchangers are the minimum initial cost, minimum operation cost, maximum effectiveness, minimum pressure drop, minimum heat transfer area, minimum weight, or material. Using the data modeling, the optimization of a heat exchanger can be transformed into a constrained optimization problem and then solved by modern optimization algorithms. In this chapter, the thermal design and optimization of shell and tube heat exchangers are presented.

### Keywords

- log-mean temperature difference (LMTD)
- effectiveness-number of transfer units (ε-NTU)
- genetic algorithm (GA)
- particle swarm optimization (PSO)

## 1. Introduction

Heat exchangers are systems used to transfer heat between fluids with different temperatures. These devices have vast applications in many areas, such as refrigeration, heating, and air conditioning systems, power plants, chemical processes, food industry, automobile radiators, and waste heat recovery units. Heat exchangers can be classified according to different criteria such as construction, flow arrangement, heat transfer mechanism, etc [1]. The heat exchanger design can be divided into two main categories, thermal and hydraulic design and mechanical design. In thermal and hydraulic design, the focus is on calculating an adequate surface area transfer a certain amount of heat, pressure dope, pumping power work, etc. The goal of the mechanical design is to design the mechanical integrity of the exchanger, as well as designing various pressure and non-pressure components. In this chapter, the thermal and hydraulic design of heat exchangers is presented. To achieve better performance of heat exchangers, they optimize based on their application. Heat exchanger optimization can be performed using different optimization algorithms. Since most heat exchanger optimization problems are nonlinear, using traditional methods such as linear and dynamic programming and steepest descent may not lead to the desired solution and may even fail. Also, most traditional methods need gradient information to solve an optimization problem. On the other hand, advanced optimization algorithms are developed, which are gradient-free. Several advanced optimization methods, such as genetic algorithm, non-nominated sorting genetic algorithm, bio-geography-based optimization, particle swarm optimization, Jaya algorithm, and teaching-learning-based optimization, can be more efficient in solving an optimization problem. However, each of these methods has its advantages and disadvantages, which are discussed in the optimization section. In this chapter, genetic algorithm and particle swarm optimization are discussed in detail due to the vast applications that arise from their acceptable accuracy, as well as short computational time [2].

## 2. Basic equation of heat transfer

In most heat transfer problems, hot and cold fluids are divided by a solid wall. In this case, the mechanism of heat transfer from hot fluid to the cold fluid can be categorized into three steps:

Heat transfer from the hot fluid to the wall by convection.

Heat transfer through the wall by conduction.

Heat transfer from the wall to the cold fluid by convection.

Figure 1 shows a schematic of heat transfer between two fluids. As it can be seen, thermal resistance (* R*) is present at each stage of the transfer. Thermal resistance is a thermal (physical) property that indicates the resistance of each material to heat transfer due to temperature differences that can be calculated from [3]:

Where * L*is the thickness of the wall,

*is the cross-sectional area in which heat transfer occurs, and*A

*and*K

*are conduction and convection heat transfer coefficient, respectively.*h

Heat transfer in each stage can be calculated as follows [4]:

Where:

* R* = thermal resistance for convection in the hot side.

_{c,H}

* R* = fouling resistance of hot side.

_{f,H}

* R* = wall resistance.

_{w}

* R* = fouling resistance of cold side.

_{f,C}

* R* = thermal resistance for convection in cold side.

_{c,C}

Table 1 shows the fouling resistance of the most common fluids used in heat exchangers. The overall heat transfer coefficient can be obtained from Eq. (2) as follows [4]:

Gas and vapors | Fouling Factor | Liquids | Fouling Factor | |||
---|---|---|---|---|---|---|

Manufactured Gas | 0.01 | Industrial organic heat transfer media | 0.001 | |||

Engine Exhaust Gas | 0.01 | Refrigerating liquids | 0.001 | |||

Steam (non-oil bearing) | 0.0005 | Molten heat transfer salts | 0.0005 | |||

Exhaust Steam (oil bearing) | 0.001 | Hydraulic fluid | 0.001 | |||

Refrigerant Vapors (oil bearing) | 0.002 | |||||

Compressed Air | 0.002 | Fuel oil | 0.005 | |||

Industrial Organic Heat Transfer Media | 0.001 | Engine lube oil | 0.001 | |||

Transformer oil | 0.001 | |||||

Acid Gas | 0.001 | Quench oil | 0.004 | |||

Solvent Vapors | 0.001 | Vegetable oils | 0.003 | |||

Stable Overhead Product | 0.001 | |||||

Temperature of heating medium | 240°F≤ | 240–400°F | ||||

Atmospheric Tower Overhead Vapors | 0.001 | Temperature of water | 125°F≤ | 125°F> | ||

Light Naphthas | 0.001 | velocity | 3 ft.≤ | 3 ft.> | 3 ft. ≤ | 3 ft. > |

Vacuum Overhead Vapors | 0.002 | Seawater | 0.0005 | 0.0005 | 0.001 | 0.001 |

Natural Gas | 0.001 | Brackish water | 0.002 | 0.001 | 0.003 | 0.002 |

Overhead Products | 0.001 | Distilled water | 0.0005 | 0.0005 | 0.0005 | 0.0005 |

Coke Unit Overhead Vapors | 0.002 | Boiler blowdown | 0.002 | 0.002 | 0.002 | 0.002 |

or:

Where:

* h*and

_{H}

*= convection heat transfer coefficient of the hot and cold sides, respectively.*h

_{C}

* A*and

_{H}

*= and surface area of wall in the hot and cold side, respectively.*A

_{C}

The * r*can be calculated for flat wall and cylindrical walls using Eqs. (5) and (6), respectively.

_{w}

Where * d*is the thickness of the wall, and

_{w}

*and*r

_{o}

*are the outside and inside diameter of the wall, respectively.*r

_{i}

Total thermal resistance can be expressed as [7]:

The rate of heat transfer (* Q*) can be determined from

Where * U*is the overall heat transfer coefficient.

Where * A*is a reference area. If the heat transfer is carried out over a pipe, the inside and outside surface areas of the pipe are not equal. Hence, the

_{ref}

*must be determined (The outer surface of the pipes is usually selected).*A

_{ref}

## 3. Thermal design of heat exchangers

The thermal design of heat exchangers can be performed by several methods. The most commonly used methods are log-mean temperature difference (LMTD) and effectiveness-number of transfer units (** -**NTU) [8]. The LMTD was used to calculate heat transfer when the inlet and outlet temperatures of fluids are specified. When more than one inlet and/or outlet temperature of the heat exchanger is unknown, LMTD may be calculated by trial and errors solution. In this case, the ε

**NTU method is commonly used [3].**-

### 3.1 The log-mean temperature difference (LMTD) method

As mentioned earlier, by determining the temperature difference between hot and cold fluids, the amount of heat transfer can be calculated from Eq. (3). Figure 2 shows temperature changes of hot and cold fluids along with a heat exchanger with different types of flow configuration. As it can be seen the temperature difference between hot and cold fluids can vary along with the heat exchanger. Terminal temperatures of hot and cold fluids (* T*anf

_{H,out}

*) are very effective factors in a heat exchanger design. If*T

_{C,out}

*is lower than*T

_{C,out}

*for countercurrent flow, temperature approach occurs. In contrast, if*T

_{H,out}

*is higher than the*T

_{C,out}

*for countercurrent flow, temperature cross happens [9]. But if*T

_{H,out}

*is equal to*T

_{C,out}

*, temperature meet takes place. Based on the second law of thermodynamics, temperature cross can never take place for heat exchangers with co-current flow configuration [10]. In 1981, Wales [11] proposed a parameter G, which can be used to determine the temperature conditions in the heat exchangers. Eq. (10) defines the*T

_{H,out}

*parameters that can be between −1 and 1.*G

Since the temperature difference between hot and cold streams varies along with the heat exchanger, the basic question is which temperature difference should be considered to calculate the heat transfer rate. To answer this question, consider Figure 3 which shows heat transfer between two parallel and co-current fluids. Based on Figure 3, the heat transfer for the specified heat transfer area can be written in the form:

It can be said that the amount of heat transferred is reduced from the hot fluid and added to the cold fluid. Therefore [3]:

Where * C*and

_{p,H}

*are specific heat capacities of hot and cold fluids, respectively. The temperature difference can be written as below:*C

_{p,C}

By the combination of Eqs. (12) and (13) with Eq. (14):

By defining

Assuming the * M*is constant along with the heat exchanger:

On the other hand, by placing Eq. (16) in Eq. (11):

Where * U*is the mean overall heat transfer coefficient which can be defined as below:

_{m}

The heat transfer rate can be written as flow:

Where

The simplicity of the LMTD method has led to its use in the design of many heat exchangers by introducing a correction factor, * F*, according to Eq. (23) [12]. The F is generally expressed in terms of two non-dimensional parameters, thermal effectiveness (

*), and heat capacity ratio (*P

*). The*R

*and*P

*are defined as Eqs. (24) and (25), respectively [13]. Figure 4 shows the correction factor for common shell and tube heat exchangers.*R

### 3.2 Effectiveness-number of transfer units (ε-NTU)

When more than one of the inlet and outlet temperatures of the heat exchanger is unknown, LMTD may be calculated by trial and errors solution. Another approach to calculating the rate of heat transfer is the effectiveness number of transfer units (* ε*-NTU) method. The

*-NTU can be expressed according to Eq. (26) where*ε

*is the minimum value between the heat capacity of cold fluid (*Cp

_{min}

*) and hot fluid (*Cp,

_{C}

*). The effectiveness (*Cp,

_{H}

*) can be defined as the ratio of the actual heat transfer rate (*ε

*) and the maximum possible heat transfer rate (*q

*) according to Eq. (27).*q

_{max}

Where:

The heat transfer rate using the * ε*-NTU method can express as [14]:

Tables 2 and 3 show the effectiveness and NTU relations for heat exchangers, respectively. It should be noted that * C*is the capacity and it can be defined as follows:

_{r}

Heat exchanger type | Effectiveness relation |
---|---|

Double pipe:Parallel-flow Counter flow with C _{r}< 1Counter flow with C _{r}= 1
| |

Sell and tube:On-pass and 2, 4, … tube
| |

Cross flow (single-pass):C _{max}and C_{min}unmixedC _{max}mixed and C_{min}unmixedC _{max}unmixed and C_{min}mixed
| |

All heat exchangers with C_{r} = 0 |

Heat exchanger type | NTU relation |
---|---|

Double pipe:Parallel-flow Counter flow with C _{r}< 1Counter flow with C _{r}= 1
| |

Sell and tube:On-pass and 2, 4, … tube
| |

Cross flow (single-pass):C _{max}mixed and C_{min}unmixedC _{max}unmixed and C_{min}mixed
| |

All heat exchangers with C_{r} = 0 |

## 4. Thermal and hydraulic design of shell and tube heat exchanger

Heat exchangers can be classified according to different criteria such as construction, flow arrangement, heat transfer mechanism, etc [1]. Shell and tube heat exchangers are some of the most convenient heat exchangers due to their versatility, wide operating range, and simplicity [16]. Hence, this chapter focuses on the design of this type of heat exchanger. In the design of shell and tube heat exchangers, a lot of consideration including the number of shells and tubes, tube pitch and layout, tube passes, baffles, etc., should be taken into account. In this case, there are some methods such as Kern and Bell-Delaware to design a heat shell and tube exchanger design. Since Kern’s method offers the simplest route, this chapter is focused on this method.

### 4.1 Kern’s method

Kern’s method is based on experimental data for typical heat exchangers. In this method, it is assumed the shell flow is ideal, and leakage and bypass are negligible. Based on this flow model, only a single stream flows in the shell that is driven by baffles. This can lead to a very simple and rapid calculation of shell-side coefficients as well as pressure drop [17]. Figure 5 shows a schematic of a shell and tube heat exchanger.

### 4.2 Number of tubes

The number of tubes (* N*) can be calculated as follows:

_{t}

Where * ṁ*is the flow rate of fluid inside the tube,

_{t}

*is the density of the fluid inside the tube,*ρ

_{t}

*is the velocity of the fluid inside the tube,*ν

_{t}

*is the cross-sectional area of the tube, and*A

_{t}

*is the tube inside diameter.*d

_{i}

### 4.3 Tube-side heat transfer coefficient

The heat transfer coefficient for the tube side (* h*) is calculated as follows:

_{t}

Where * Nu*is the Nusselt number for the tube-side fluid and

_{t}

*is the thermal conductivity of the tube-side fluid. The*k

_{t}

*is a function of Reynolds number (*Nu

_{t}

*) and Prandtl number (*Re

*).*Pr

*and*Re

*can be obtained by the following:*Pr

Where * μ*is the dynamic viscosity of the tube-side fluid,

_{t}

*is the heat conductivity coefficient, and*K

*is the heat capacity of the tube-side fluid. The*C

_{p}

*can be calculated according to the type of flow as follows:*Nu

_{t}

Where * L*is the length of the tube and

*is the friction factor of the tube side, which can be calculated from*f

_{t}

The convection heat transfer coefficient in the tube is obtained based on the value of the Re_{t} from [19]:

Where * μ*is the dynamic viscosity of the tube-side fluid at the wall temperature and λ is the Darcy friction coefficient which can be defined as [19]:

_{w,t}

The tube-side pressure drop is calculated by the following:

Where * N*is the tube passes.

_{p}

### 4.4 Shell diameter

Inside sell diameter (* D*) is calculated as follows:

_{s}

Where * A*is the projected area of the tube layout expressed as an area corresponding to one tube and can be obtained from Eq. (45). Also,

*is tube pitch and*P

_{t}

*is the tube layout constant. Figure 6 depicts two common tube layouts, square pitch and triangular pitch. The*CL

*is the tube count calculation constant that accounts for the incomplete coverage of the shell diameter by the tubes, due to necessary clearances between the shell and the outer tube circle and tube omissions due to tube pass lanes for multitude pass design [15]. Eq. (46) shows the*CTP

*for different tube passes.*CTP

Combining Eq. (44) with Eq. (45) as well as defining tube pitch ratio as P_{r}, one gets:

Where * d*is tube outside the diameter. Eq. (47) can be written as follows:

_{o}

Where * A*is the outside heat transfer surface area based on the outside diameter of the tube and can be calculated from:

_{o}

The shell side flow direction is partially along the tube length and partially across to tube length or heat exchanger axis. The inside shell diameter can be obtained based on the cross-flow direction and the equivalent diameter (* D*) is calculated along the long axes of the shell. The equivalent diameter is given as follows:

_{e}

From Figure 6 the equivalent diameter for the square pitch and triangular pitch layouts are as below:

Reynolds number for the shell-side (* Re*) can be obtained as follows:

_{s}

Where * ṁ*is the flow rate of shell-side fluid,

_{s}

*is the viscosity of the shell-side fluid, and*μ

_{s}

*is the cross-flow area at the shell diameter which can be obtained as below:*A

_{s}

Where * B*is the baffle spacing and

*is the clearance between adjacent tubes. According to Figure 6*C

_{t}

*is expressed as follows:*C

_{t}

The shell-side mass flow rate (* G*) is found with:

_{s}

In Kern’s method, the heat transfer coefficient for the shell-side (* h*) is estimated from the following:

_{s}

Where * k*is the thermal conductivity of the shell-side fluid. The tube-side pressure drop is calculated by the following:

_{s}

Where * N*is the number of baffles,

_{b}

*is the density of the shell-side fluid,*ρ

_{s}

*is the viscosity of the shell-side fluid at bulk temperature, and*μ

_{b}

*is the viscosity of the tube-side fluid at wall temperature. The*μ

_{w,s}

*is the friction factor for the shell and can be obtained as follows:*f

_{s}

The wall temperature can be calculated as follows:

According to Eq. (21), the heat transfer surface area (* A*) of the shell and tube heat exchanger is obtained by the following:

The required length of the heat exchanger can be calculated based on the heat transfer surface area as follows:

## 5. Optimization of heat exchangers

The applications of heat exchangers are very different. Therefore, they are optimized based on their application. The most common criteria for optimizing heat exchangers are the minimum initial cost, minimum operating cost, maximum effectiveness, minimum pressure drop, minimum heat transfer area, minimum weight or material, etc. These criteria can be optimized individually or in combination. It is clear from the above that the optimal design of heat exchangers is based on many geometrical and operational parameters with high complexity. So it is difficult to design a cheap and effective heat exchanger. The optimization techniques are usually applied to ensure the best performance as well as lower the cost of the heat exchanger. The optimization is carried out using different techniques. Traditional techniques such as linear and dynamic programming as well as steepest descent usually fail to solve nonlinear large-scale problems. The need for gradient information is another drawback of traditional techniques. Therefore, it is not possible to solve non-differentiable functions using these methods. To overcome these difficulties, advanced optimization algorithms are developed which are gradient-free. Several advanced optimization methods, such as genetic algorithm (GA) [21], non-nominated sorting GA (NSGA-II) [22], bio-geography-based optimization (BBO) [23], particle swarm optimization (PSO) [24], Jaya algorithm, and teaching–learning-based optimization (TLBO) [25], had been used for the optimization of heat exchangers by many researchers each of which has its advantages and disadvantages. Using GA, it is possible to solve all optimization problems, which can be described with the chromosome encoding and solves problems with multiple solutions. But in order to use GA, it is necessary to set a number of specific algorithmic parameters such as jump probability, selection operator, cross probability. NSGA-II has an explicit diversity preservation mechanism and elitism prevents an already found Pareto optimal solution from being removed, but crowded comparison can limit the convergence and it needs the tuning of algorithmic-specific parameters including mutation probability, crossover probability, etc. The optimization using BBO is also effective and it inhibits the degradation of the solutions, but poor exploiting the solutions is the main drawback of this method. The PSO is a heuristic and derivative-free technique that has the character of memory but it needs the tuning of algorithmic specific parameters and plurality of the population is not enough to achieve the global optimal solution. Similarly, TLBO and Jaya need the tuning of their own algorithmic-specific parameters [26].

Generally, an optimization design starts by selecting criteria (quantitatively) to minimize or maximize, which is called an objective function. In an optimization design, the requirements of a particular design such as required heat transfer, allowable pressure drop, limitations on height, width and/or length of the exchanger are called constraints. Several design variables such as operating mass flow rates and operating temperatures can also participate in an optimization design [15]. The single target optimization can be expressed as [27]:

Where * f(x)*is the objective function,

*≥ 0 is the inequality constraint, and*g

_{i}(x)

*≥ 0 is the equality constraint. Multi-objective combination optimization can be indicated as [27]:*h

_{i}(x)

Using the data modeling, the optimization of a heat exchanger can be transformed into a constrained optimization problem and then solved by modern optimization algorithms. In this chapter, the focus is on GA and PSO because many researchers mentioned that these algorithms lead to remarkable savings in computational time and have an advantage over other methods in obtaining multiple solutions of the same quality. So it gives more flexibility to the designer [2].

### 5.1 Genetic algorithm

GA is a search heuristic that is inspired by Charles Darwin’s theory of survival of the fittest, which explains inferior creatures pass away and superior creatures remain [28]. In GA, sets of design variables are codified by sequences with fixed or variable lengths, similar to chromosomes or individuals in biological systems. Each chromosome is formed of several design variables, which are known as genes. In repetitive processes such as GA, each repetitive stage is a generation and a collection of solutions associated with each generation is a population. Generally, the initial population is generated randomly [29]. In GA statistical methods are used to achieve optimum points. In the process of natural selection, populations are selected based on their fitness. A new population is formed using genetic operations containing selection, crossover, mutation, etc. This cycle continues until a certain result is achieved or the stop criterion is satisfied [30]. Figure 7 shows the flowchart of GA and the steps of binary GA are discussed below [32].

The initial population of GA includes binary numbers generated randomly which are chromosomes or GA strings, consisting of bits called genes. Actually, the initial population is the probable solution to the optimization problem. The number of gens (* n*) assigned to represent a variable in the chromosome depends on the precision ϵ and the range of the variable [

_{g}

*], and is given by*x

_{min}, x

_{max}

The fitness value of each GA string is examined by first determining the decoded values of the variables * D*, and next the corresponding real values are obtained as follows:

The fitness function values are then computed knowing the real values of design variables.

In this step, chromosomes with better fitness values to participate in the crossover are selected. Several selection modes such as roulette wheel selection or proportionate selection, rank-based selection, and tournament selection can be used in this step [33]. In proportionate selection, the probability of a chromosome to be selected is directly proportional to its fitness value. Hence, the chromosome having a better fitness value has a higher chance of selection for reproduction. This may result in premature convergence of the solution because there is a chance of losing diversity. The tournament selection is faster compared with the other two selection methods. In this method, * n*chromosomes are randomly picked from the population of solutions, where n represents the tournament size. The chromosome having the best fitness value is copied to the mating pool and all the n GA strings are returned to the population. This process is repeated for obtaining all the individuals of the mating pool.

The genes are exchanged between two-parent chromosomes in the crossover step, which leads to a new set of solutions, called children. The crossover operation represents the selection pressure or exploitation of fit chromosomes for even better solutions. The crossover probability (* P*) specifies the number of individuals taking part in the crossover operation, and this control parameter value is optimally chosen as nearly equal to 1.0. Several schemes of crossover such as single-point crossover, two-point crossover, multipoint crossover, and uniform crossover can be used in this step. A comparison of these methods is given in the literature [34].

_{c}

Mutation means the change of a bit from 0 to 1 and from 1 to 0 in the solution chromosome, which is used for the exploration of new solutions. It helps to come out of the local basin and search for a global solution. The mutation probability * P*specifies the number of mutations and is commonly kept very low. Because if its value is high, the qualified solutions may be lost. The range of

_{m}

*is given as*P

_{m}

Where * l*represents the length of the GA string. Steps 2, 3, 4, and 5 are repeated until the termination criterion (maximum number of generations or desired precision of solution) is met.

### 5.2 Particle swarm optimization

The PSO is inspired by the way fish and birds swarm search for food [35]. In this method, each particle represents one solution to a problem and they aim to find optimum points in a search space. This method is also based on the behavior of birds that they use to find their orientation. Based on this direction, the collective location of the swarm and the best individual location of particles per time are calculated and a new search orientation is composed of these two orientations and the previous orientation. In a search space of the D dimension, the best individual location of a particle and the best location of the overall particle are defined as Eqs. (69) and (70), respectively.

The best location in the vicinity of each particle is given as below:

Displacement of particles after determination their velocity is as follows:

The best individual location of the particle and the best collective location of particles as two springs connected to the particle are used to model the force entered in the particle. The first spring is directed to the best individual experience and the second spring is directed to the best swarm experience. Eq. (74) shows the force entered in the particle.

Where * C*and

_{1}

*are acceleration coefficients. The particle velocity at dimension*C

_{2}

*and the next repetition can be obtained as follows [35]:*d (ν

_{id})

This shows the velocity of particle i at the star topology or global best. The rand_{1} and rand_{2} are random numbers that have a constant distribution in the range 0–1. Figure 8 shows the flowchart of PSO and the steps of PSO are discussed below [32].

The swarm of potential solutions is generated with random positions and velocities. The i* th*particle in D-dimensional space may be denoted as

*and*X

_{i}= (x

_{i1}, x

_{i2},...., x

_{id})

*, where*i = 1, 2,..., N

*denotes the size of the swarm.*N

The corresponding fitness values of the particles are evaluated.

The best individual location of a particle (

The velocity and position of the ith particle are updated based on Eq. (75).

Here, an example of the design and optimization of a shell and tube heat exchanger is presented. This example was used by Karimi et al. (2021) [36]. Their aim was to minimize the total annual cost (* C*) for a shell and tube heat exchanger based on optimization algorithms. The total annual cost is the sum of the initial cost for the construction (

_{tot}

*) of the heat exchanger and the cost of power consumption in the shell and tube heat exchanger (*C

_{i}

*). Hence, the total annual cost was considered as an objective function that should be minimized using GA and PSO. Process input data and physical properties for this case study are presented in Table 4. Also, bounds for design parameters are listed in Table 5. The objective function can be written as follows:*C

_{od}

27.8 | 95 | 40 | 750 | 2840 | 0.00034 | 0.19 | 0.00033 | 0.00038 | |

68.9 | 25 | 40 | 995 | 4200 | 0.0008 | 0.59 | 0.0002 | 0.00052 |

0.015 | 0.051 | |

0.1 | 1.5 | |

0.05 | 0.5 |

The results show that the use of PSO has been led to lower Ctot, which means that the minimization of cost function was performed better using this algorithm. Also, the use of PSO resulted in lower Δp and A as well as higher U (Table 6).

PSO | GA | |
---|---|---|

2.6871 | 3.9089 | |

_{o} | 0.015063 | 0.015 |

0.49967 | 0.49989 | |

_{s} | 0.81143 | 0.74105 |

_{t} | 1238 | 1365.5 |

_{t} | 0.83349 | 0.893 |

_{t} | 27,909 | 13,386 |

_{t} | 5.69 | 5.69 |

_{t}^{2} | 3740.8 | 4639.7 |

_{t} | 0.0212 | 0.0073 |

_{t} | 4730 | 5191.3 |

_{e} | 0.0141 | 0.0106 |

_{s} | 0.498 | 0.499 |

_{s} | 15,489 | 11,716 |

_{s} | 5.1 | 5.1 |

_{s}^{2} | 9075.7 | 1648.9 |

_{s} | 0.313 | 0.353 |

_{s} | 21,355 | 18,033 |

^{2} | 900.98 | 686.71 |

^{2} | 198.78 | 252.58 |

_{i} | 44,116 | 50,737 |

_{o} | 2561.5 | 1085.2 |

_{od} | 2340 | 6685 |

_{tot} | 46,456 | 57422.51 |

## 6. Conclusion

This chapter has discussed the thermal design and optimization of shell and tube heat exchangers. The basic equations of heat transfer were investigated and log-mean temperature difference (LMTD) and effectiveness-number of transfer units (* ε*-NTU) were presented. The thermal design was focused on Kern’s method. In this method, it is assumed the shell flow is ideal, and leakage and bypass are negligible. Based on this flow model, only a single stream flows in the shell that is driven by baffles. This can lead to a very simple and rapid calculation of shell-side coefficients and pressure drop. The optimization of heat exchangers is presented based on the genetic algorithm (GA) and particle swarm optimization (PSO) due to the recommendation of these methods by man researchers because of quick convergence and obtaining multiple solutions.

## Nomenclature

Total heat transfer area

_{s}

Cross-flow area at the shell diameter

Baffle spacing

Heat capacity

_{i}

Capital investment cost

_{o}

Total operating cost

_{od}

Total discounted operating cost

_{r}

Capacity ratio

_{t}

Clearance between adjacent tubes

_{tot}

Total annual cost

Tube layout constant

Tube count calculation constant

_{w}

Wall thickness

_{s}

Shell diameter

_{s}

Friction factor for shell

_{s}

Mass flow rate of shell-side fluid

Convection heat transfer coefficient

Thermal conductivity

Heat conductivity coefficient

Tube length

Logarithmic mean temperature difference

Flow rate of shell-side fluid

_{t}

Number of tubes

Number of heat transfer unit

Thermal effectiveness

_{r}

Tube pitch ratio

_{t}

Tube pitch constant

Prandtl number

Actual heat transfer rate

_{max}

Maximum possible heat transfer rate

Heat transfer rate

_{i}

Inside diameter

_{o}

Outside diameter

Thermal resistance

Reynolds number

Temperature

Overall heat transfer coefficient

_{m}

Mean overall heat transfer coefficient

## Greek symbols

Difference

Effectiveness

Density

Dynamic viscosity

## Subscripts

Bulk

Cold

Convection

Cold

Inside

Maximum

Minimum

Outside

Shell side

Tube side

## References

- 1.
Ghalandari, M., et al., Applications of intelligent methods in various types of heat exchangers: a review. Journal of Thermal Analysis and Calorimetry, 2021. 145(4): p. 1837-1848 - 2.
Lahiri SK, Khalfe NM,Wadhwa SK. Particle swarm optimization technique for the optimal design of shell and tube heat exchangers. Chemical Product and Process Modeling. 2012; 7 (1), Article 14 - 3.
Ezgi CN. Basic Design Methods of Heat Exchanger. In: Heat Exchangers-Design, Experiment and Simulation. Rijeka: IntechOpen; 2017 - 4.
Chen Q, Ren J. Generalized thermal resistance for convective heat transfer and its relation to entransy dissipation. Chinese Science Bulletin. 2008; 53 (23):3753-3761 - 5.
Marner W, Suitor J. Survey of Gas-Side Fouling in Industrial Heat-Transfer Equipment. Final Report. Pasadena, CA (USA): Jet Propulsion Lab; 1983 - 6.
Whitaker S. Fundamental Principles of Heat Transfer. Oxford, United kingdom: Pergamon press, Elsevier; 2013 - 7.
Çengel YA, Ghajar AJ. Heat and Mass Transfer: Fundamentals [and] Applications. New York City, United States: McGraw-Hill Education; 2020 - 8.
Borjigin S et al. Coupling ε-NTU method for thermal design of heat exchanger in cabinet cooling system. Applied Thermal Engineering. 2019; 159 :113972 - 9.
Smith R. Chemical Process Design And Integration. New York City, United States: John Wiley & Sons; 2005 - 10.
Vengateson U. Design of multiple shell and tube heat exchangers in series: E shell and F shell. Chemical Engineering Research and Design. 2010; 88 (5–6):725-736 - 11.
Wales, R. and W. RE, Mean Temperature Difference in Heat Exchangers. New York, United state Chemical Engineering Journal; 1981 - 12.
Fakheri A. Log mean temperature correction factor: an alternative representation. In: Proceedings of the International Mechanical Engineering Congress and Exposition. New Orleans, Louisiana: Citeseer; 2002 - 13.
Gulyani BB, Jain A, Kumar S. Optimal synthesis of multipass heat exchanger without resorting to correction factor. International Journal of Mechanical, Aerospace, Industrial, Mechatronic and Manufacturing Engineering. 2011; 5 (5):898-904 - 14.
Iu I et al. Applying the effectiveness-NTU method to elemental heat exchanger models. ASHRAE Transactions. 2007; 113 (1):504-513 - 15.
Shah RK, Sekulic DP. Fundamentals of Heat Exchanger Design. New York City, United States: John Wiley & Sons; 2003 - 16.
Silaipillayarputhur K, Khurshid H. The design of shell and tube heat exchangers–A review. International Journal of Mechanical and Production Engineering Research and Development. 2019; 9 (1):87-102 - 17.
Alperen MA, Kayabaşi E, Kurt H. Detailed comparison of the methods used in the heat transfer coefficient and pressure loss calculation of shell side of shell and tube heat exchangers with the experimental results. Energy Sources, Part A: Recovery, Utilization, and Environmental Effects. 2019:1-20 - 18.
Standard I. Specification for Shell and Tube Type Heat Exchangers. New Delhi: BIS 2007; 2007 - 19.
Jalilirad S, Cheraghali MH, Ashtiani HAD. Optimal design of shell-and-tube heat exchanger based on particle swarm optimization technique. Journal of Computational Applied Mechanics. 2015; 46 (1):21-29 - 20.
Anescu, G., A DMSACO Approach to Economic Heat Exchanger Design. Politechnica University of Bucharest in Romania - 21.
Guo J, Cheng L, Xu M. Optimization design of shell-and-tube heat exchanger by entropy generation minimization and genetic algorithm. Applied Thermal Engineering. 2009; 29 (14–15):2954-2960 - 22.
Yin Q et al. Optimization design and economic analyses of heat recovery exchangers on rotary kilns. Applied Energy. 2016; 180 :743-756 - 23.
Hadidi A, Nazari A. Design and economic optimization of shell-and-tube heat exchangers using biogeography-based (BBO) algorithm. Applied Thermal Engineering. 2013; 51 (1–2):1263-1272 - 24.
Patel V, Rao R. Design optimization of shell-and-tube heat exchanger using particle swarm optimization technique. Applied Thermal Engineering. 2010; 30 (11–12):1417-1425 - 25.
Rao RV, Patel V. Multi-objective optimization of heat exchangers using a modified teaching-learning-based optimization algorithm. Applied Mathematical Modelling. 2013; 37 (3):1147-1162 - 26.
Rao RV et al. Design optimization of heat exchangers with advanced optimization techniques: A review. Archives of Computational Methods in Engineering. 2020; 27 (2):517-548 - 27.
Yao J. A review of industrial heat exchange optimization. In: IOP Conference Series: Earth and Environmental Science. JIOP Publishing Ltd in Bristol, England; 2018 - 28.
Haupt RL, Haupt SE. Practical Genetic Algorithms. New York City, United States: John Wiley & Sons; 2004 - 29.
Najafi H, Najafi B, Hoseinpoori P. Energy and cost optimization of a plate and fin heat exchanger using genetic algorithm. Applied Thermal Engineering. 2011; 31 (10):1839-1847 - 30.
Konak A, Coit DW, Smith AE. Multi-objective optimization using genetic algorithms: A tutorial. Reliability Engineering & System Safety. 2006; 91 (9):992-1007 - 31.
Sadeghzadeh H, Ehyaei M, Rosen M. Techno-economic optimization of a shell and tube heat exchanger by genetic and particle swarm algorithms. Energy Conversion and Management. 2015; 93 :84-91 - 32.
Swayamsiddha S. Bio-inspired algorithms: principles, implementation, and applications to wireless communication. In: Nature-Inspired Computation and Swarm Intelligence. Cambridge, Massachusetts, United States: Academic Press; 2020. pp. 49-63 - 33.
Blickle T, Thiele L. A Comparison of Selection Schemes used in Genetic Algorithms , TIK-Report No. 11. Computer Engineering and Communication Networks Lab (TIK) Zurich, Switzerland; 1995 - 34.
Spears WM, Anand V. A study of crossover operators in genetic programming. In: International Symposium on Methodologies for Intelligent Systems. Berlin, Heidelberg: Springer; 1991 - 35.
Shi Y. Particle swarm optimization: developments, applications and resources. In: Proceedings of the 2001 Congress on Evolutionary Computation (IEEE Cat. No. 01TH8546). Seoul, Korea (South); 2001 - 36.
Karimi, H., H. Ahmadi Danesh Ashtiani, and C. Aghanajafi, Optimization Techniques for Design a shell and tube heat exchanger from economic view. AUT Journal of Mechanical Engineering, 2021. 5 (2): 9-9