IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE)

e-ISSN: 2278-1684,p-ISSN: 2320-334X, Volume 13, Issue 3 Ver. VII (May- Jun. 2016), PP 17-22
www.iosrjournals.org

CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat

Exchanger with Helical-Tap Inserts at Annulus of Inner Pipe
1
Bharat Bhushan Verma, 2Saurabh Kumar
1 2
ME Student, Associate Prof., Mechanical Department, Raipur Institutes of Technology, Raipur, India

Abstract: In heat exchanger for increasing the rate of heat transfer with helical-tape inserts have useful. Generally, helical-
tape inserts causes swirl flow introduces at outside of inner tube which continuously disrupts the thermal boundary layer of
fluid on the tube. This analysis has done on a single unit of pipe of heat exchanger those are used in a heat exchanger and
this investigation is useful to increase the thermal characteristics of a heat exchanger. To analyze the characteristics of
helical-tap inserts at annulus of inner pipe, a 3-D analytical model has been developed. From the analysis it is cleared that
there are good relation between Nusselt no. and friction factor for enhancing the heat transfer. SST k-ω turbulent model is to
be selected for simulation because it gives better turbulent model. From analysis, helical-tape inserts increase the heat
transfer rate with expectation of pressure drop. In this work an analysis has been done on heat transfer phenomenon of
helical-tap inserts at annulus of the inner pipe.
Keywords: Double pipe heat exchanger, Heat transfer augmentation techniques, Helical-tape insert, Pitch length,
Numerical investigation, Computational fluid dynamics, Turbulence modeling, Friction factor, Nusselt number.

I. Introduction
The procedure of designing the heat exchanger is difficult, because it required correct analysis of heat transfer rate,
flow rate, drop in pressure and other factor calculation that give the result for long term established and economically
preference of the equipment. The challenge is to designing a heat exchanger for better performance and makes it to compact
to gain a more heat transfer. Normally a heat exchanger used in power plants, chemical plants, AC equipments, Freezer and
other plants which provide or remove heat from different types of fluid. This involves huge investments annually for both
operation and capital costs. It is required to reduce the overall dimensions and characteristics of heat exchanger. The need to
analyze and optimize the heat exchanger to conserved and developed it.
Heat exchanger with helical-tap inserts at annulus of inner pipe, from the velocity vector, it is observed that the
flow of water in plain tube is straight line but in case of helical-tap flow is swirl, swirl flow in tube causes the surface area of
effected heat transfer in increases, thereby rate of heat transfer increases but pressure drop also increases because of flow
abstraction. This pressure drops varies along with Reynolds no. on helical-tap. To predict the performance heat exchanger is
relating to the governing parameters such as surface area for transferring heat overall heat transfer coefficient and
temperature difference. Assuming there is no transfer of heat to the surrounding and negligible KE and PE.
mh - mass of hot fluid entering (kg), mc - mass of cold fluid entering (kg), Ch - Sp. heat of hot fluid entering
(J/kg K),Cc – Sp. heat of cold fluid entering (J/kg K), th1 - temperature of hot fluid entering (K), tc1 - temperature of cold fluid
entering (K), th2 - temperature of hot fluid exit (K), tc2 - temperature of cold fluid exit (K), dh- hydraulic diameter.
Heat rejected by hot fluid–
𝑄𝑕 = 𝑚𝑕 𝐶𝑕 (𝑡𝑕1 − 𝑡𝑕2 ) (1.1)
Heat absorbed by cold fluid
𝑄𝑐 = 𝑚𝑐 𝐶𝑐 (𝑡𝑐1 − 𝑡𝑐2 ) (1.2)
Heat exchange by two fluids
𝑄 = 𝑈𝐴𝜃𝑚 (1.3)
Where, U- Overall heat transfer coefficient, A- Effective heat transfer area, θm - appropriate means temperature
difference across heat exchanger. As there is variation in temperature difference of hot and cold fluids point to point, so that
by the concept of mean temperature difference the term ϴm has introduced which is appropriate mean temperature difference
across heat exchanger or known as log mean temperature difference.
For parallel flow log mean temperature difference is given by
θ2− θ1
θm = θ2 (1.4)
𝐼𝑛 ( θ1 )
For counter flow log mean temperature difference is given by
θ1− θ2
θm = θ1 (1.5)
𝐼𝑛 ( θ2 )
Where, θ1 = th1 - tc1, θ2 = th2 - tc2

Of the many enhancement techniques that can be employed, swirl flow generation by means of full-length helical-
tape inserts is found to be extremely effective [1], [2]. Significant heat transfer improvement can be obtained, particularly in
laminar flows. Other examples of techniques that promote swirl flows include curved ducts, tangential fluid injection, and
twisted or convoluted ducts. Their thermal-hydraulic characteristics, heat transfer improvement potential, and typical
applications have been outlined.
Helical-tap – we know that the heat transfer is increase considerably for flow is well mixed and stirred. Because of this
principle of development for equipment to enhancement technique of heat exchanger helical-tap inserts are used at annulus
of inner pipe.

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 17 | Page

 CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap..

Heat transfer increase because of:

a. Tap reduce the hydraulic diameter, cause affect and enhance the coefficient of heat transfer.
b. The helical-tap causes a tangential velocity component hence speed of flow of fluid is increase near the wall surface.
c. There may be heat transfer from the tape.

Figure 1: CFD model of helical tape at annulus of inner pipe

Figure 1.1 shown of helical-tap, is described by the helical twisting nature of tap providing the fluid longer flow
region or greater time for transferring heat, the helical force for bulk flow are forcing for generation a secondary circulation
because of well mixed swirl flow increase the convective heat transfer[5].

The swirl flow which is in fully developed nature, performance of helical-tap and functional relation are given as below

𝑓 = ∅ 𝑅𝑒, 𝑦, 𝛿/𝑑 (1.6)

𝑁𝑢 = ∅ 𝑅𝑒, 𝑦, 𝛿/𝑑 , 𝑃𝑟 (1.7)
Based on a fundamental balance between inertia, viscous, and tape-geometry helical curvature induced forces, it is
proposed that tape-induced swirl flows can be scaled by a swirl parameter defined as [1]
𝑆𝑤 = 𝑅𝑒/ 𝑦 (1.8)
Where, the swirl Reynolds number is based on the swirl velocity, and
𝑅𝑒𝑠 = 𝜌𝑉𝑠 𝑑/𝜇 (1.9)

II. Numerical Investigation Of Heat Exchanger With Helical Tape Inserts At Annulus Of Inner Pipe
For analysis of our work the input data and boundary condition will be taken in [1] which is experimentally investigate the
heat transfer phenomenon.

Table I Input data for double pipe heat exchanger

Length of tube 2.2 m
Inner diameter of inner pipe, di 0.022 m
Outer diameter of inner pipe, do 0.026 m
Inner diameter of outer pipe, Di 0.054 m
Outer diameter of outer pipe, do 0.058 m
Material of pipe Copper
Inner pipe fluid Cold water(300K)
Annulus fluid Hot water(353k)
Table II Properties of water
Density, ρ 998.2 kg/m3
Specific Heat Capacity, Cp 4182 J/kg K
Thermal Conductivity, k 0.6 W/m K
Viscosity, µ 1.003×10-3 kg/m s

Table III Boundary condition for inner fluid

Inlet condition Velocity Inlet (0.376 m/s)
Outlet condition Pressure Outlet
Inlet Temperature 300 K

Table IV Boundary condition for annulus fluid

Inlet condition Velocity Inlet (Varies from 0.127 to 0.557 m/s)
Outlet condition Pressure Outlet
Inlet Temperature 353

III. Data Reduction Equations

The area weighted average temperature and static pressure were noted at the inlet and outlet surfaces of the pipe.
After setting all the above mentioned parameters simulation has been done and after simulation Friction Coefficients for all
the mentioned or considered Reynolds Number have been calculated as per the following relation and verified with the
findings of Patnala shankra Rao [1].
The correlation for the Calculation of „Friction Factor (f)‟ Colburn‟s Equation [1] is:
𝑓 = 0.046(𝑅𝐿 −0.2 )…………………………………… (4.1)

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 18 | Page

 CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap..

Reynolds No. RL is:

𝜌𝑉𝑑
𝑅𝐿 = ………………...…………………………… (4.2)
µ
Heat transfer rate is given by:
𝑄1 = 𝑚 × 𝐶𝑃 × (𝑇𝑐𝑜 − 𝑇𝑐𝑖 )………………………...… (4.3)
𝑄2 = 𝑚 × 𝐶𝑃 × 𝑇𝑕𝑖 − 𝑇𝑕𝑜 ………………………..… (4.4)
𝑄 +𝑄
𝑄𝐴𝑣𝑔 = 1 2 ……………...……………………..…… (4.5)
2
Area:
𝐴 = 𝜋 × 𝑑𝑕 × 𝐿………..……………………………... (4.6)
LMTD is given by:
𝑇𝑕 𝑖 −𝑇𝑐𝑜 −(𝑇𝑕 𝑜 −𝑇𝑐𝑖 )
𝐿𝑀𝑇𝐷 = 𝑇 𝑕 𝑖 −𝑇 𝑐𝑜
……..…………………… (4.7)
ln ( )
𝑇 𝑕 𝑜 −𝑇 𝑐𝑖
Convective heat transfer coefficient is:
𝑄
𝑕= ………………………………………….... (4.8)
𝐴×𝐿𝑀𝑇𝐷
Nusselt No. is given by:
𝑕×𝑑 𝑕
𝑁𝑢 = ……………………………………….......... (4.9)
𝑘

IV. Result And Discussion

The analysis of heat exchanger is done by the using different pitch length of 50, 100, 150, 200 and 250 mm for 40
design points. This analysis is for hot water at velocity range 0.127 to 0.577 m/s and for cold water at constant velocity 0.367
m/s. The inlet temperature of cold water and hot water is respectively 300 K and 353 K [1].

Figure 2: CFD Model for heat exchanger with helical tap of different pitch length inserts.

Figure 3: Inlet velocity of hot domin

Figure 4: Inlet velocity of cold domin

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 19 | Page

 CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap..

Figure 5: Average outlet temperature of hot water

Figure 6: Average outlet temperature of cold water

Figure 7: Meshed view of helical tape

Figure 8: Friction Factor Vs Reynolds No. with helical tape

Figure 9: Graph represents Nusselt No. Vs Reynolds No. with different pitch length

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 20 | Page

 CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap..

35
30

Heat transfer,KW
25
20
15
10
5
0
0 50 100 150 200 250 300

Pitch length, mm

Figure 10: Graph represents optimum value of pitch length

It is cleared from the above result heat transfer and heat transfer rate are increases with decrease the pitch length.
The above result is performed for the different pitch length and plot the comparison of Reynolds No. and Nusselt Number
and Reynolds No. and Friction factor. It is cleared that the Friction factor decrease with increase in Reynolds Number and
Nusselt No. increase with Increase in Reynolds Number. From the above result is also optimizing the helical tape pitch
length for different velocity. The maximum heat transfer achieved at minimum pitch length which is 50 mm and maximum
velocity of 0.557 m/s. The present research also predicts that by increasing the mass flow rate of fluids, there is miner
variation in heat transfer rate.

V. Conclusion
It is clear that insertion of a helical tape in a plain tube increase the thermal performance of the tube and
furthermore if the pitch length of helical tape is reduce increase in surface it increases the tube‟s thermal performance more.
Area Weighted Average of the fluid‟s temperature at the outlet of the tube has been increased due to the insertion of a helical
tape in the tube. The reason for the increment of these parameters is that, due to the insertion of a helical tape a swirl flow is
created in the pipe which helps the fluid to take more and more heat from the tube wall.
So it may be concluded that, modifications should be done in such a way so that average temperature as well as
heat flux both increases. To do this, optimization procedure may be adopted to optimize different parameters to achieve the
desired goal.
Scope for future work: Further detailed studies can be carried out in this area either through experiments or with
the aid of software. Nusselt number and friction factor values can be obtained for helical tap with the same pitch at different
velocity and similarly for helical tap with the same velocity and different pitch in order to study the effect of helical tap pitch
length on Nusselt number and friction factor. Some other inserts may be used and similar investigations can be done and the
values compared to those of helical tap inserts.

References
[1]. Patnala Sankara Rao, K Kiran Kumar(2014), Numerical and Experimental Investigation of Heat Transfer Augmentation in Double Pipe Heat
Exchanger with Helical and Twisted Tape Inserts, Volume-4.
[2]. K.Sivakumar, K.Rajan(2014-2015) Performance Analysis of Heat Transfer and Effectiveness on Laminar Flow with Effect of Various Flow Rates,
Volume-7.
[3]. Amarjit Singh and Satbir S. Sehgal (2013), Thermohydraulic Analysis of Shell-and-Tube Heat Exchanger with Segmental Baffles.
[4]. Kamlesh R. Raut, Prof. H.S. Farkade (2014), Convective Heat Transfer Enhancement in Tube Using Insert – A Review, Volume-2.
[5]. Prof. Naresh B. Dhamane, Prof. Mathew V. Karvinkoppa, Prof. Murtuza S. Dholkhwala (2012), Heat Transfer Analysis of Helical Strip Insert with
Regularly Spaced Cut Sections Placed Inside a Circular Pipe, Volume-2.
[6]. Neeraj kumar Nagayach, Dr. Alka Bani Agrawal (2012), Review of Heat Transfer Augmentation in Circular and Non-Circular Tube, Volume-2.
[7]. Prof. P. B. Dehankar, Prof. N. S. Patil (2014), Heat Transfer Augmentation - A Review for Helical Tape Insert, Volume-3.
[8]. K.Sivakumar, K. Rajan, S. Murali, S. Prakash, V. Thanigaive, T. Suryakumar (2015), Experimental Investigation Twisted Tape Insert on Laminar
Flow With Uniform Heat Flux For Enhancement Of Heat Transfer.
[9]. P. B. Malwadkar, Lalit Pawar Pratik Satav (2014), Experimental Investigation of Heat Transfer Performance of Matrix Coil Wire Inserts using CFD ,
Volume-1.
[10]. Dr. D. S. Kumar “Heat and Mass Transfer”. S. K. Kataria and Sons. 7th Edition. Chapter no. 14, pp 681-725

Appendix
Table V Numerical result at 50mm pitch
S Velocity Average Convective Reynolds Friction Nusselt
No. of heat heat No. factor No.
annulus transfer transfer RL F Nu
fluid KW coefficient
(m/s) h
(KW/m2K)
1 0.127 7.04181 1.15676 3538.98 0.00897 53.98
2 0.208 11.10660 1.78283 5796.12 0.00813 83.19
3 0.26 13.98420 2.32348 7245.16 0.00777 108.42
4 0.346 18.89485 3.11277 9641.63 0.00734 145.26
5 0.4 22.24153 3.75048 11146.40 0.00713 175.02
6 0.433 24.34318 4.14922 12065.97 0.00702 193.63
7 0.52 29.22962 5.02292 14110.47 0.00680 234.40
8 0.577 32.82157 5.72580 15657.20 0.00666 267.20

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 21 | Page

 CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap..

Table VI Numerical result at 100mm pitch

S Velocity Average Convective Reynolds Friction Nusselt
No. of heat heat No. factor No.
annulus transfer transfer RL f Nu
fluid KW coefficient
(m/s) h
(KW/m2K)
1 0.127 6.33410 0.97627 3538.98 0.00897 45.559
2 0.208 10.3044 1.58114 5796.12 0.00813 73.786
3 0.26 13.5127 2.13015 7245.16 0.00777 99.407
4 0.346 18.9805 3.04332 9641.63 0.00734 142.021
5 0.4 21.8896 3.57807 9952.14 0.00729 166.976
6 0.433 24.1505 4.00747 12065.97 0.00702 187.015
7 0.52 29.1781 4.87827 14490.32 0.00676 227.65
8 0.577 32.8194 5.55982 16078.68 0.00663 259.46

Table VII Numerical result at 150mm pitch

S Velocity Average Convective Reynolds Friction Nusselt
No. of heat heat No. factor No.
annulus transfer transfer RL F Nu
fluid KW coefficient
(m/s) h
(KW/m2K)
1 0.127 5.94475 0.8868 3538.98 0.00902 41.384
2 0.208 11.2788 1.8226 5796.12 0.00817 85.054
3 0.26 14.47137 2.3996 7245.16 0.00777 111.98
4 0.346 19.472 3.2434 9641.63 0.00734 151.35
5 0.4 22.6103 3.7846 11146.4 0.00713 176.61
6 0.433 24.9284 4.2096 12065.97 0.00702 196.44
7 0.52 29.4910 4.9762 14490.32 0.00676 232.22
8 0.577 32.6799 5.5098 16078.68 0.00662 257.12

Table VIII Numerical result at 200mm pitch

S Velocity Average Convective Reynolds Friction Nusselt
No. of heat heat transfer No. factor No.
annulus transfer coefficient RL f Nu
fluid KW h
(m/s) (KW/m2K)
1 0.127 6.5382 1.0267 3538.98 0.00897 40.91
2 0.208 10.7652 1.6922 5796.12 0.00813 78.96
3 0.26 13.7418 2.1890 7245.16 0.00777 102.15
4 0.346 18.4353 2.9554 9294.83 0.00739 137.91
5 0.4 21.5496 3.4813 11146.40 0.00713 162.46
6 0.433 24.0553 3.9819 11749.68 0.00705 185.82
7 0.52 29.9605 5.1237 14110.47 0.00680 239.10
8 0.577 31.6990 5.2577 16078.68 0.00662 245.35

Table IX Numerical result at 250mm pitch

S Velocity Average Convective Reynolds Friction Nusselt
No. of heat heat transfer No. factor No.
annulus transfer coefficient RL f Nu
fluid KW h
(m/s) (KW/m2K)
1 0.127 4.7826 0.7153 3538.98 0.00897 33.38
2 0.208 7.5287 1.1951 5796.12 0.00813 55.77
3 0.26 9.6671 1.4319 7245.16 0.00777 66.82
4 0.346 13.1439 1.9699 9294.83 0.00734 91.92
5 0.4 15.6422 2.3993 11146.40 0.00713 111.96
6 0.433 17.2848 2.6673 11749.68 0.00702 124.47
7 0.52 21.2171 3.3257 14110.47 0.00680 155.19
8 0.577 23.5976 3.5927 16078.68 0.00666 167.65

DOI: 10.9790/1684-1303071722 www.iosrjournals.org 22 | Page