 Open access peer-reviewed chapter

# Heat Exchanger Design and Optimization

Written By

Shahin Kharaji

Submitted: July 30th, 2021Reviewed: September 14th, 2021Published: November 20th, 2021

DOI: 10.5772/intechopen.100450

From the Edited Volume

## Heat Exchangers

Edited by Laura Castro Gómez, Víctor Manuel Velázquez Flores and Miriam Navarrete Procopio

Chapter metrics overview

View Full Metrics

## 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)

## 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 :

R=LKA,forconductionR=1hA,forconvectionE1

Where Lis the thickness of the wall, Ais the cross-sectional area in which heat transfer occurs, and Kand hare conduction and convection heat transfer coefficient, respectively.

Heat transfer in each stage can be calculated as follows :

Q=T1T2Rc,H=T2T3Rf,H=T3T4Rw=T4T5Rf,C=T5T6Rc,CE2

Where:

Rc,H = thermal resistance for convection in the hot side.

Rf,H = fouling resistance of hot side.

Rw = wall resistance.

Rf,C = fouling resistance of cold side.

Rc,C = thermal resistance for convection in cold side.

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 :

Gas and vaporsFouling Factor
hrft2FBtu
LiquidsFouling Factor
hrft2FBtu
IndustrialIndustrial liquids
Manufactured Gas0.01Industrial organic heat transfer media0.001
Engine Exhaust Gas0.01Refrigerating liquids0.001
Steam (non-oil bearing)0.0005Molten heat transfer salts0.0005
Exhaust Steam (oil bearing)0.001Hydraulic fluid0.001
Refrigerant Vapors (oil bearing)0.002Industrial oils
Compressed Air0.002Fuel oil0.005
Industrial Organic Heat Transfer Media0.001Engine lube oil0.001
Chemical ProcessingTransformer oil0.001
Acid Gas0.001Quench oil0.004
Solvent Vapors0.001Vegetable oils0.003
Petroleum ProcessingTemperature of heating medium240°F≤240–400°F
Atmospheric Tower Overhead Vapors0.001Temperature of water125°F≤125°F>
Light Naphthas0.001velocity3 ft.≤3 ft.>3 ft. ≤3 ft. >
Natural Gas0.001Brackish water0.0020.0010.0030.002

### Table 1.

Fouling factors for different types of fluid [5, 6].

Q=T1T6Rc,H+Rf,H+Rw+Rf,C+Rc,CE3

or:

Q=T1T61hHAH+rf,HAH+rwAH+rfCAC+1hCACE4

Where:

hHand hC = convection heat transfer coefficient of the hot and cold sides, respectively.

AHand AC = and surface area of wall in the hot and cold side, respectively.

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

rw=dwKA,forflatwallE5
rw=lnrori2πLK,forcylindricalwallE6

Where dwis the thickness of the wall, and roand riare the outside and inside diameter of the wall, respectively.

Total thermal resistance can be expressed as :

Rt=1hHAH+rf,HAH+rwAH+rf,CAC+1hCACE7

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

Q=UAΔTE8

Where Uis the overall heat transfer coefficient.

U=11hHAHAref+rf,HAHAref+rwAHAref+rf,CACAref+1hCACArefE9

Where Arefis 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 Arefmust be determined (The outer surface of the pipes is usually selected).

## 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) . 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.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 (TH,outanf TC,out) are very effective factors in a heat exchanger design. If TC,outis lower than TH,outfor countercurrent flow, temperature approach occurs. In contrast, if TC,outis higher than the TH,outfor countercurrent flow, temperature cross happens . But if TC,outis equal to TH,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 . In 1981, Wales  proposed a parameter G, which can be used to determine the temperature conditions in the heat exchangers. Eq. (10) defines the Gparameters that can be between −1 and 1. Figure 2.Temperature changes of hot and cold fluids along with a heat exchanger with different types of flow configuration.
G=TH,outTC,outTH,inTC,in,G>0temperatureapproachG=0temperaturemeetG<0temperaturecrossE10

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:

dQ=UΔTdAE11

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

dQ=Cp,HdTHE12
dQ=Cp,CdTCE13

Where Cp,Hand Cp,Care specific heat capacities of hot and cold fluids, respectively. The temperature difference can be written as below:

dΔT=dTHdTCE14

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

dΔT=dQCp,HdQCp,CE15

By defining 1M=1Cp,C+1Cp,H, Eq. (15) can be written as follow:

dQ=MdΔTE16

Assuming the Mis constant along with the heat exchanger:

0QdQ=MΔTinΔToutdΔTQ=MΔToutΔTinE17

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

MdΔTΔT=UdA=AUdAAE18
A0AUdAA=MΔTinΔToutdΔTΔTUmA=MlnΔToutΔTinE19

Where Umis the mean overall heat transfer coefficient which can be defined as below:

Um=0AUdAAE20

The heat transfer rate can be written as flow:

Q=UmAΔTLMTDE21

Where ΔTLMTDis the logarithmic mean temperature difference (LMTD) which can be defined as below:

ΔTLMTD=ΔToutΔTinlnΔToutΔTinE22

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) . The F is generally expressed in terms of two non-dimensional parameters, thermal effectiveness (P), and heat capacity ratio (R). The Pand Rare defined as Eqs. (24) and (25), respectively . Figure 4 shows the correction factor for common shell and tube heat exchangers. Figure 4.Correction factor for common shell and tube heat exchangers . (a) One-shell pass and 2, 4, 6, etc. (any multiple of 2), tube pass. (b) Two-shell pass and 4, 8, 12, etc. (any multiple of 4), tube pass. (c) Single-pass cross-flow with both fluids unmixed. (d) Single-pass cross-flow with one fluid unmixed and other unmixed.
Q=UAFΔTLMTD,0<F<1E23
P=TC,outTC,inTH,inTC,inE24
Rcr=TH,inTH,outTC,outTC,inE25

### 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 Cpminis the minimum value between the heat capacity of cold fluid (Cp,C) and hot fluid (Cp,H). The effectiveness (ε) can be defined as the ratio of the actual heat transfer rate (q) and the maximum possible heat transfer rate (qmax) according to Eq. (27).

NTU=UAmCpminE26
ε=qqmaxE27

Where:

q=Cp,HTH,inTH,out=Cp,CTC,outTC,inE28
qmax=CpminTH,inTC,inE29

The heat transfer rate using the ε-NTU method can express as :

Q=εCPminTH,inTC,inE30

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

Heat exchanger typeEffectiveness relation
Double pipe:
• Parallel-flow

• Counter flow with Cr < 1

• Counter flow with Cr = 1

ε=1expNTU1+Cr1+Cr
ε=1expNTU1Cr1CrexpNTU1Cr
ε=NTU1+NTU
Sell and tube:
• On-pass and 2, 4, … tube

ε=21+Cr+1+C21+expNTU1+C21expNTU1+C21
Cross flow (single-pass):
• Cmax and Cmin unmixed

• Cmax mixed and Cmin unmixed

• Cmax unmixed and Cmin mixed

ε=1exp1CrNTU0.22expCrNTU0.781ε=1Cr1expCr1expNTU
ε=1exp1Cr1expCrNTU
All heat exchangers with Cr = 0ε=1expNTU

### Table 2.

Effectiveness relation for heat exchangers .

Heat exchanger typeNTU relation
Double pipe:
• Parallel-flow

• Counter flow with Cr < 1

• Counter flow with Cr = 1

NTU=ln1ε1+Cr1+Cr
NTU=1Cr1lnε1εCr1
NTU=ε1ε
Sell and tube:
• On-pass and 2, 4, … tube

NTU=11+C2ln2ε1Cr1+C22ε1Cr+1+C2
Cross flow (single-pass):
• Cmax mixed and Cmin unmixed

• Cmax unmixed and Cmin mixed

NTU=ln1+ln1εCrCr
NTU=lnCrln1ε+1Cr
All heat exchangers with Cr = 0NTU=ln1ε

### Table 3.

NTU relation for heat exchangers .

Cr=CPminCPmaxE31

## 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 . Shell and tube heat exchangers are some of the most convenient heat exchangers due to their versatility, wide operating range, and simplicity . 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 . Figure 5 shows a schematic of a shell and tube heat exchanger. Figure 5.Schematic of a shell and tube heat exchanger a) fixed-tube b) floating-head c) removable U-tube .

### 4.2 Number of tubes

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

Nt=4mt.ρtνtπdi2E32

Where tis the flow rate of fluid inside the tube, ρtis the density of the fluid inside the tube, νtis the velocity of the fluid inside the tube, Atis the cross-sectional area of the tube, and diis the tube inside diameter.

### 4.3 Tube-side heat transfer coefficient

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

ht=NutktdiE33

Where Nutis the Nusselt number for the tube-side fluid and ktis the thermal conductivity of the tube-side fluid. The Nutis a function of Reynolds number (Re) and Prandtl number (Pr). Reand Prcan be obtained by the following:

Ret=ρtνtdtμtE34
Prt=CpμtKE35

Where μtis the dynamic viscosity of the tube-side fluid, Kis the heat conductivity coefficient, and Cpis the heat capacity of the tube-side fluid. The Nutcan be calculated according to the type of flow as follows:

Nut=ft/2RetPrt1.07+12.7ft/21/2Prt2/31;for:104<Re<5×106&0.5<Pr<200E36
Nut=1.86RetPrtdiL1/3;for:RetPrtdiL1/3>2&0.48<Pr<16700E37

Where Lis the length of the tube and ftis the friction factor of the tube side, which can be calculated from

ft=1.58lnRet3.282E38

The convection heat transfer coefficient in the tube is obtained based on the value of the Ret from :

ht=ktdi3.657+0.0677RetPrtdiL1.31+0.1PrtRet+diL0.3;forRet<2300E39
ht=ktdiλ8Ret1000Prt1+12.7λ8Prt0.6711+diL0.67;for2300<Ret<10000E40
ht=ktdi0.027Ret0.8Prt13μtμw,t0.14;forRet>10000E41

Where μw,tis the dynamic viscosity of the tube-side fluid at the wall temperature and λ is the Darcy friction coefficient which can be defined as :

λ=1.82log10log10Ret1.642E42

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

ΔPt=4ftLNpdi+4Npρtμt22E43

Where Npis the tube passes.

### 4.4 Shell diameter

Inside sell diameter (Ds) is calculated as follows:

Ds=4ANtCTPπE44

Where Ais the projected area of the tube layout expressed as an area corresponding to one tube and can be obtained from Eq. (45). Also, Ptis tube pitch and CLis the tube layout constant. Figure 6 depicts two common tube layouts, square pitch and triangular pitch. The CTPis 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 . Eq. (46) shows the CTPfor different tube passes.

A=Pt2CLCL=1forsquarepitchlayoutCL=0.866fortriangularpitchlayoutE45
CTP=0.93foronetubepassCTP=0.9fortwotubepassCTP=0.85forthreetubepassE46

Combining Eq. (44) with Eq. (45) as well as defining tube pitch ratio as Pr, one gets:

Ds=4Prdo2CLNtCTPπE47
Pr=PtdoE48

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

Ds=4Pr2CLAodoCTPπ2LE49

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

Ao=πdoNtLE50

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 (De) is calculated along the long axes of the shell. The equivalent diameter is given as follows:

De=4×freeflowareawettedperimeterE51

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

De=4Pt2πdo24πdo;forsquarepitchtubeE52
De=4Pt234πdo28πdo2;fortriangularpitchtubeE53

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

Res=ms.AsDeμsE54

Where sis the flow rate of shell-side fluid, μsis the viscosity of the shell-side fluid, and Asis the cross-flow area at the shell diameter which can be obtained as below:

As=DsPtB×CtE55

Where Bis the baffle spacing and Ctis the clearance between adjacent tubes. According to Figure 6Ctis expressed as follows:

Ct=PtdoE56

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

Gs=ms.AsE57

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

hs=0.36ksDeRes0.55Prs1/3for2×103<Res=GsDeμs<1×106E58

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

ΔP=fsGs2Nb+1Ds2ρsDeμbμw,s0.14E59

Where Nbis the number of baffles, ρsis the density of the shell-side fluid, μbis the viscosity of the shell-side fluid at bulk temperature, and μw,sis the viscosity of the tube-side fluid at wall temperature. The fsis the friction factor for the shell and can be obtained as follows:

fs=exp0.5760.19lnRes;for400,Res<1×106E60

The wall temperature can be calculated as follows:

Tw=12TH,in+TH,out2+TC,in+TC,out2E61

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

A=QUmFΔTLMTDE62

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

L=AπdoNtE63

## 5. Optimization of heat exchangers

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 . The single target optimization can be expressed as :

minfxgix0,j=1,2,,Jhkx0,k=1,2,,KE64

Where f(x)is the objective function, gi(x) ≥ 0 is the inequality constraint, and hi(x) ≥ 0 is the equality constraint. Multi-objective combination optimization can be indicated as :

minfx=f1xf2xfmxgix0,j=1,2,,Jhkx0,k=1,2,,KE65

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 .

### 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 . 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 . 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 . Figure 7 shows the flowchart of GA and the steps of binary GA are discussed below .

Step 1: Initialization of population

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 (ng) assigned to represent a variable in the chromosome depends on the precision ϵ and the range of the variable [xmin, xmax], and is given by

ng=log2xminxmaxE66

Step 2: Fitness evaluation

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:

x=xmin+xmaxxmin2ng1E67

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

Step 3: Reproduction/selection

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 . 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, nchromosomes 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.

Step 4: Crossover

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 (Pc) 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 .

Step 5: Mutation

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 Pmspecifies 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 Pmis given as

0.1lPm1lE68

Where lrepresents 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 . 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.

P1=C1Pi1Pi2PiDE69
g1=C1g1g2gDE70

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

n=ni1ni2niDE71

Displacement of particles after determination their velocity is as follows:

xt=xt1+νtE72
νt=νt1+Ft1E73

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.

Fi1=C1Pi1xi1+C2gi1xi1E74

Where C1and C2are acceleration coefficients. The particle velocity at dimension d (νid)and the next repetition can be obtained as follows :

νidt=ωνidt1+C1rand1Pidt1xidt1+C2rand2nidt1xidt1E75

This shows the velocity of particle i at the star topology or global best. The rand1 and rand2 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 .

Step 1: Initialization

The swarm of potential solutions is generated with random positions and velocities. The ithparticle in D-dimensional space may be denoted as Xi = (xi1, xi2,...., xid)and i = 1, 2,..., N, where Ndenotes the size of the swarm.

Step 2: Fitness evaluation

The corresponding fitness values of the particles are evaluated.

Step 3: Determination of personal and global best

The best individual location of a particle (P1) ​is sorted, the particle having the best fitness value is determined for the current generation, and the best location (g1) is updated.

Step 4: Velocity and position update

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) . Their aim was to minimize the total annual cost (Ctot) 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 (Ci) of the heat exchanger and the cost of power consumption in the shell and tube heat exchanger (Cod). 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:

Shell side: methanol27.8954075028400.000340.190.000330.00038
Tube side: sea water68.9254099542000.00080.590.00020.00052

### Table 4.

Process input data and physical properties for three case studies .

ParametersLower valueUpper value
Tubes outside diameters(m)0.0150.051
Shell diameters(m)0.11.5
Central baffle spacing(m)0.050.5

### Table 5.

Bounds for design parameters .

Ctot=Ci+CodE76

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).

PSOGA
L (m)2.68713.9089
do(m)0.0150630.015
B (m)0.499670.49989
Ds(m)0.811430.74105
Nt(m)12381365.5
νt(m/s)0.833490.893
Ret27,90913,386
Prt5.695.69
ht(W/m2K)3740.84639.7
ft0.02120.0073
ΔPt(Pa)47305191.3
De(m)0.01410.0106
νs(m/s)0.4980.499
Res15,48911,716
Prs5.15.1
hs(W/m2K)9075.71648.9
fs0.3130.353
ΔPs(Pa)21,35518,033
U (W/m2K)900.98686.71
A (m2)198.78252.58
Ci($)44,11650,737 Co($)2561.51085.2
Cod($)23406685 Ctot($)46,45657422.51

### Table 6.

Optimal parameter of heat exchanger using GA and PSO algorithms .

## 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.

## Conflict of interest

The authors declare no conflict of interest.

## Nomenclature

A

Total heat transfer area

As

Cross-flow area at the shell diameter

B

Baffle spacing

Cp

Heat capacity

Ci

Capital investment cost

Co

Total operating cost

Cod

Total discounted operating cost

Cr

Capacity ratio

Ct

Ctot

Total annual cost

CL

Tube layout constant

CTP

Tube count calculation constant

dw

Wall thickness

Ds

Shell diameter

Fs

Friction factor for shell

Gs

Mass flow rate of shell-side fluid

h

Convection heat transfer coefficient

k

Thermal conductivity

K

Heat conductivity coefficient

L

Tube length

LMTD

Logarithmic mean temperature difference

ṁs

Flow rate of shell-side fluid

Nt

Number of tubes

NTU

Number of heat transfer unit

P

Thermal effectiveness

Pr

Tube pitch ratio

Pt

Tube pitch constant

Pr

Prandtl number

q

Actual heat transfer rate

qmax

Maximum possible heat transfer rate

Q

Heat transfer rate

ri

Inside diameter

ro

Outside diameter

R

Thermal resistance

Re

Reynolds number

T

Temperature

U

Overall heat transfer coefficient

Um

Mean overall heat transfer coefficient

## Greek symbols

Δ

Difference

ε

Effectiveness

ρ

Density

μ

Dynamic viscosity

ν

Velocity

## Subscripts

b

Bulk

C

Cold

c

Convection

H

Cold

i

Inside

max

Maximum

min

Minimum

o

Outside

s

Shell side

t

Tube side

Written By

Shahin Kharaji

Submitted: July 30th, 2021Reviewed: September 14th, 2021Published: November 20th, 2021