Optimization of Heat Transfer Using CFD Simulation For Concentric Helical Coil Heat Exchanger For Constant Temperature Outer Wall
Optimization of Heat Transfer Using CFD Simulation For Concentric Helical Coil Heat Exchanger For Constant Temperature Outer Wall
Optimization of Heat Transfer Using CFD Simulation For Concentric Helical Coil Heat Exchanger For Constant Temperature Outer Wall
BACHELOR OF TECHNOLOGY
IN
MECHANICAL ENGINEERING
BY
SAGAR DAS
110ME0287
ROURKELA – 769008
(2013-2014)
1|Page
Optimization of heat transfer using CFD simulation
for concentric helical coil heat exchanger
for constant temperature outer wall
BACHELOR OF TECHNOLOGY
IN
MECHANICAL ENGINEERING
BY
SAGAR DAS
110ME0287
2|Page
CERTIFICATE
This is to certify that the thesis entitled “Optimization of heat transfer using CFD simulation
for concentric helical coil heat exchanger for constant temperature outer wall” submitted by
Sagar Das (Roll no. 110ME0287) in partial fulfillment of the requirements for the award of
Bachelor of Technology degree in Mechanical Engineering at the National Institute of
Technology, Rourkela is an authentic work carried out by him under my supervision and
guidance.
To the best of my knowledge, the matter embodied in the thesis has not been submitted to any
other University/Institute for the award of any Degree or Diploma.
3|Page
ACKNOWLEDGEMENT
Sagar Das
110ME0287
Bachelor of Technology,
Mechanical Engineering Dept.
4|Page
CONTENTS
Certificate 3
Acknowledgement 4
List of Tables 6
Abstract 7
Nomenclature 8
1. Introduction 9
1.1. Heat Exchanger 10
1.2. Types of Heat Exchanger 10
1.3. Flow arrangement in recuperative Heat Exchanger 10
1.4. Tubular Heat Exchanger 11
1.5. Helical Tube-in-Tube Heat Exchanger 11
1.6. Characteristics of Helical Coil 12
1.7. Advantages and Disadvantages 13
1.8. Critical Radius of Insulation of Cylindrical Surface 14
1.9. Aim of the present work 15
1.10. Methodology 15
2. Literature Survey 16
3. CFD Modeling 19
3.1.Geometry 20
3.1.1. Sketching 20
3.1.2. Sweeping 21
3.1.3. Boolean Operation 21
3.1.4. Merging 21
3.2. Mesh 23
3.2.1. Mapped Face Meshing 23
5|Page
3.2.2. Edge Sizing 23
3.2.3. Named Selection 23
3.3. Solution 23
3.3.1. Problem Set-up 23
3.3.2. Models 24
3.3.3. Material 24
3.3.4. Cell Zone Conditions 25
3.3.5. Boundary Conditions 25
3.3.6. Reference Values 26
3.3.7. Solution Methods 26
3.3.8. Solution Control and Initialization 26
3.3.9. Measure of Convergence 27
3.3.10. Run Calculation 27
4. Results and Discussion 28
4.1. – 4.5. Tabulation and Plots 29
5. Conclusions 44
6. References 45
List of Tables:
6|Page
ABSTRACT
7|Page
NOMENCLATURE
8|Page
CHAPTER 1
INTRODUCTION
9|Page
1.1. Heat Exchangers
Heat transfer between the flowing fluids is an important physical process to be focused upon,
thus heat exchangers plays a vital role in our day to day life. They have a wide range of
application like they can be used in a varied type of installations, food processing, transportation,
domestic applications, process industries, nuclear power plant, air conditioning, compact heat
exchangers, recovery process, HVACs, refrigeration, etc.
The main purpose of the heat exchanger is to frame an efficient method of heat transfer/exchange
from one fluid to another, by simple direct or indirect contact. The heat transfer mostly occurs by
three principles i.e. convection, conduction and radiation. The heat transfer through radiation in a
heat exchanger is generally not taken into account, as it is comparatively negligible to the heat
exchange by conduction and convection. The process conduction occurs when there is a
temperature gradient between the solid wall. It can be maximized by selecting a critical radius of
insulation of the wall and a high conductive material. Convection plays an important role in the
heat exchanger performance. Forced convection accelerates the heat transfer in a heat exchanger
from a moving stream of fluid to the wall of pipe or vise-versa.
In the applications of heat exchanger, improvement is focused on the efficiency, substantial cost,
material saving and space.
10 | P a g e
1.4. Tubular Heat exchanger
They are mainly of circular cross section. The circular cross section provides flexibility in the
design parameters like length, diameter, thickness of tube and their arrangement can be modified
easily and can be arranged to multiple complex shapes. It is mostly used for single phase i.e.
liquid-to-liquid heat transfer. It is further classified into:
The tube-in-tube helical heat exchanger consists of one tube oriented concentrically inside
another with a greater radius. The flow configuration may be of the following two types i.e.
parallel or counter flow. It can also be arranged in a lot of parallel and series configurations to
meet different requirements of heat transfer. The helical arrangement stands out to be used in
various industrial applications. Though this configuration has been widely used, through
knowledge should be focused upon the heat transfer coefficient, temperature gradient,
pressure drop with various flow patterns are of much importance. The curvature forms a
secondary flow in the tubes, which is just normal to the direction of flow to the primary axis.
Heat transfer occurring between the wall and the fluid in increased substantially by this
secondary flow, which also offers a greater area for the heat transfer within a small compact
space, with higher heat transfer coefficient. Types of flow in curved pipes has been focused and
effect of Prandtl number and Reynolds number has been related on the following flow patterns
and also on Nusselt number.
11 | P a g e
1.6. Characteristics of Helical Coil [1]
The helical coil has a pitch of height ‘H’, diameter of tube as ‘2r’, the coil diameter ‘2Rc’,
curvature ratio as ‘i’ i.e. ratio of tube and coil diameter as ‘r/Rc’. The helix angle ‘2a’ is the
angle between its projection on a surface and measuring angle between the coil, shown by ‘v’
shape.
The centrifugal force is governed by the curvature of the tubular coil, while torsion occurs due to
the helix angle or pitch. The centrifugal force leads to the development of secondary flow in the
helical heat exchanger.
12 | P a g e
1.7. Advantages & Disadvantages:
1. Helical tubes have high heat transfer co-efficient compared to straight tubes and
are more compact.
2. It increases the contact area and time for heat exchange between the two fluids,
which leads to higher rate of heat exchange.
3. The tube curvature causing secondary flow pattern, which is perpendicular to the
main stream flow acts as an additional convective heat transfer mechanism.
4. Coils generally give better heat transfer performance, as they do have higher process side
coefficient and lower wall resistance.
5. The entire surface area of the curved helical tube is exposed to moving fluid, which thus
eliminates any dead-zones, which is a common drawback in the shell-tube type heat
exchanger.
6. The spring-like coil in helical heat exchanger eliminates any thermal expansion and
thermal shock problems, which has its application in high pressure operations.
7. Fouling is comparatively less in helical coil type than shell and tube type because of
greater turbulence created inside the curved pipes.
Disadvantages of coils:
13 | P a g e
1.8. CRITICAL RADIUS OF INSULATION OF CYLINDRICAL
SURFACE [9]
Heat transfer decreases when we add insulation to the tubular pipe. Heat transfer
rate decreases to a greater extent, as we increase the thickness of insulation. If the surface or area
of heat transfer is kept constant and insulation is added, then it gradually increases the thermal
resistance, thus increases thermal resistance.
The resistance of conduction increases as the addition of the insulation layer, but
at the same time the convection resistance decreases because of the notifying increase in the
outer surface area for heat transfer due to convection. The heat transfer from the pipe may
decrease or increase, depending on which effect, i.e. convection or conduction dominates. The
rate of heat transfer to the surrounding atmosphere from the insulated pipe is expressed as:
14 | P a g e
1.9. Aim of the Present Work:
Numerical analysis were carried out to determine heat transfer characteristics for the
tube-in-tube helical heat exchanger. CFD analysis was carried out and their variation on thermal
and hydraulic characteristics were analyzed, with varying Reynolds number (hot fluid) and
varying tube-to-coil diameter ratios for a given flow velocity of cold fluid. The analysis was
carried with Ansys 13.0 Fluent, for turbulent counter-flow with fluid water. The correlations for
heat transfer and drop in pressure were analyzed. The objective of the project is also to obtain a
better and quantitative insight to the heat transfer process that generally occurs when a fluid
flows inside a helical coiled tube. The materials were decided and fluid taken was water and the
material for the pipe was taken as copper because it has better conducting properties.
1.10. Methodology:
Literature Review
Project definition and description
Selection of variation of parameters of tube.
Selection of material along with specifications such as thickness and size.
Boundary conditions: Inlet (velocity, mass flow rate, hydraulic diameter for backflow
parameters), Outlet ( Pressure), wall ( constant temperature and no slip condition)
Observing and analyzing the characteristics of flow.
15 | P a g e
.
CHAPTER 2
LITERATURE REVIEW
16 | P a g e
S.D. Sancheti & DR.P.R .Suresh have worked on the Experimental and CFD
estimation of the heat transfer in helically coiled heat exchanger. His work focused on the fluid –
to – fluid heat transfer. He validated the basic methodology of CFD analysis in a heat exchanger,
without considering actual properties of fluid, a constant value was established instead. For
various boundary conditions, the heat transfer characteristics was compared for a helical coil. He
found that specification for constant heat flux and constant temperature boundary condition
doesn’t yield desired modeling for an actual possible heat exchanger. So, heat exchanger was
analyzed considering conjugate temperature dependent and heat transfer properties. The
fabrication of an experimental set up for the heat transfer characteristics was developed.
Experimental results were compared with the results of CFD calculation using CFD package i.e.
FLUENT 6.2. Use of constant values for the basic thermal and transport properties in the heat
exchanger resulted in prediction of inaccurate heat transfer coefficients. From the results
obtained from experiment a correlation was developed for the calculation of inner heat transfer
coefficient in a helical coil heat exchanger. [1]
Rahul Kharat, Nitin Bhardwaj and R. S Jha worked on the Development of heat
transfer coefficient correlation for concentric helical coil heat exchanger. The existing correlation
was found to result in a large discrepancies with increase in the gap between the two concentric
coils when they were compared with the experimental results. In their work, CFD simulations
and with Fluent 6.3.26 was compared with the experimental data, which was used to develop
improved correlation for the heat transfer coefficient. Mathematical model was also developed
for analyzing the data, which were obtained from experimental results and CFD to accounts for
the effects, by different parameters and functional variables like coil gap, coil diameter and tube
diameter. Using numerical technique Optimization was done for the heat transfer coefficient by
the new correlation, which fits the experimental data within an error band of 3-4%. [2]
J.S. Jayakumar , S.M Mahajani, J.C Mandal, Kannan N. Iyer and P.K. Vijayan
worked on the Thermal hydraulic characteristics of air-water two-phase flows in helical pipes.
He worked on the Two fluid Eulerian-Eulerian calculation in Fluent 6.3 for the analysis. The
parameters which influence the nature of flow are pitch coil diameter, pitch and diameter in
helical coils. CFD analysis was carried out and their variation on thermal and hydraulic
characteristics by changing the inlet void fraction for a given flow velocity. The correlations for
17 | P a g e
heat transfer and drop in pressure were analyzed. Estimation of inner heat transfer coefficient by
changing the void fraction and flow velocity, results in reduction in ‘h’ is less below 5% and
significant above 15% void fraction [3]
18 | P a g e
CHAPTER – 3
CFD MODELING
19 | P a g e
Computational fluid dynamics (CFD) is the study of system, which starts with
the basic construction of the desired geometry and then meshing for modeling the desired
dominion. Generally, the geometry is thus simplified for CFD studies. Mesh is the
discretization/division of the domain into further small volumes where then the equations are
simplified and solved by the help of simple iterative methods. Modeling starts with describing
the boundary and the initial conditions for the domain and leads to the modeling of the total
system. It is thus followed by analysis of the obtained results, discussions, calculations and
conclusions.
3.1. Geometry:
Heat exchanger with a concentric tube-in-tube is created in the ANSYS 13.0 workbench design
module. It is a heat exchanger, with counter-flow. Initially, the fluid flow (fluent) module is
selected from the workbench. The design model opens up as a new window on clicking the
geometry option. Scale is taken in mm.
3.1.1. Sketching
Initially, there are 3 reference planes, i.e. XY-plane, YZ-plane and ZX-plane. Any of the plane is
selected say the XY-plane is selected. Then the sketch button which appears on the toolbar menu
is selected 5 times as 5-sketches are to be drawn and then integrated to form the geometry, which
comprises of 4 circles and a line with required boundary conditions to constraint the sketch.
So, the first sketch is selected from the tree outline and then selected the
sketching button and a circle is selected with the center and radii needed to constraint it. The
circle is drawn on the +ve, X-axis and then selected the modeling option to constraint it by
giving the radii as 5 mm and distance from the center of circle to the origin as 50 mm (taken D/d
= 10) i.e. ratio of coil diameter and tube diameter as 10 in the first case.
This process is followed by the next 3 concentric circles with radii as 6 mm,
11 mm and 12 mm. The 5-6 mm is the thickness of the inner tube and 11-12 mm is the thickness
of outer tube. The 5th sketch drawn is a line on the +ve Y-axis starting from the origin with the
height = (length of pitch of helix * no. of turns). So, here 1.5 turns and pitch = 30 mm is taken,
so height = 45 mm is taken. Now, generate button is clicked. So, the basic sketch is created.
20 | P a g e
3.1.2. Sweeping
The Sketch 1, 2, 3 & 4 are swept along the path, i.e. line drawn in sketch 5, by editing the details
view on the tray by changing the operation to “add frozen” to construct the 3D model. The
helical sweep is of 1.5 turns as the twist specification which is defined as number of turns. Then
the helical tubes are created by generating the sweep.
The Boolean operation is selected from the create toolbar 3 times as four volumes are to be
created without any volume being intersected. So, the 1st Boolean operation is selected and
subtract operation is selected with 4th volume – 3rd volume, with 3rd volume retained. Likewise
the other Boolean operations are performed.
3.1.4. Merging
After sweep and Boolean operation, it shows that the model has 4 parts and 4 bodies. So, in the
merge operation, all the 4 parts are then selected using control and selecting new parts and thus
merged as 1 part. Then, it will show1 part and 4 bodies. The 4 bodies obtained within 1 part are
thus named and edited the type fluid/solid as follows:
Naming of various
parts of the body Part Of The Model State Type
with state type
Part number
1 Inner_Fluid Fluid
2 Inner_Pipe Solid
3 Outer_Fluid Fluid
4 Outer_Pipe Solid
Thus, the geometry is created which is then generated, refreshed and then updated.
Save the project.
21 | P a g e
Geometry of Helical coil with 1.5 turns
22 | P a g e
Meshing of helical coil
3.2. Meshing :
This mesh generally contains mixed cells (i.e. Tetra and Hexahedral cells) with both triangular
faces and quadrilateral faces at boundaries. Care should be taken to use the structured hexahedral
cells to an extent as high as possible. A fine mesh is thus generated. For this fine meshing, the
edges and the regions of preferably high temperatures and pressure gradient are meshed finely.
The mesh can be applied by creating or clicking on mapped face meshing 4 times, to create 4
mapped face mesh. The opposite faces of the inlet and outlet are selected and applied after one
another.
3.2.2. Sizing
Then the edge sizing is defined for the no. of divisions we want to divide the entire volume into
small parts to perform the finite element analysis. Each edge is defined with no. of even divisions
increasing with the edge diameter.
Different surfaces i.e. inlet, outlet, boundary of the solid are then named as per the required inlet and
outlet for cold and hot fluids. The outermost wall is named as constant temperature surface. Generate
the mesh and then close the window. Refresh and then update the project on the workbench. Then,
open the setup. ANSYS Fluent Launcher will then open in a small window. Set the dimension as 3-D,
option: Double Precision and click OK. The Fluent window will gradually open.
3.3.Solution:
The mesh is automatically checked and the quality is obtained, if it is compatible. The type
analysis is then changed - Pressure Based type. The velocity is changed to absolute and the time
to steady state.
23 | P a g e
3.3.2. Models
Energy is set - ON position. As we prefer to analyze turbulent flow, so the viscous model is
selected: “k-ε” model (2 equations).
3.3.3. Materials
The edit option is then clicked to add the water-liquid (for the hot and cold fluid) and copper
(tubes) to the list - fluid and solid resp. from the specified fluent database.
24 | P a g e
No. of Coils n -
Pitch of helix P mm
Coil dia D mm
Boundary conditions are used according to the need of the model. The inlet and outlet conditions
are defined as velocity inlet and pressure outlet. As this is a counter-flow with two tubes so there
are two inlets and two outlets. The walls are separately specified with respective boundary
conditions. No slip condition is considered for each wall. Except the tube walls each wall is set
to zero heat flux. The outer wall constant temperature is taken as 313 K.
25 | P a g e
3.3.6 Reference Values
The inner fluid inlet is selected from the list of “compute from”. The values are:
Scheme = Simplec
Gradient = Least Square Cell Based
Pressure = Linear
Momentum = Power Law
Turbulent Kinetic Energy = Power Law
Turbulent Dissipation Rate = Power Law
Energy = Second Order Upwind
26 | P a g e
The solution initialization method, is then set to Standard Initialization where the reference frame
is set to the Relative cell zone. The inner fluid inlet is selected from the list and solution is
initialized.
For a stable convergence through the simulation, the criteria is made strict to get an accurate
result. For this, the residuals are edited as per the given table:
Residual
Table 4 Residuals Variable
x-velocity 10-4
y-velocity 10-4
z-velocity 10-4
Continuity 10-4
Specific dissipation energy/ dissipation energy 10-5
Turbulent kinetic energy 10-5
Energy 10-4
The number of iterations is set to 750 and the solution is thus calculated and different contours,
vectors and plots are thus obtained.
27 | P a g e
CHAPTER- 4
28 | P a g e
4.1. D/d = 10 ( d = 10 mm )
So, for Turbulent flow, v > 0.936 m/s. So, the velocity of cold fluid is taken 2.5 m/s and velocity
of hot fluid is taken from 1 m/s, increased in steps of 0.2 m/s to 2 m/s.
Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )
Tf vs Nodes
338.8
338.6
338.4
338.2
0 10000 20000 30000 40000 50000 60000 70000 80000 90000
29 | P a g e
For Divisions (30, 32, 60, 62), it is Grid Independent, The graphs corresponding to it:
Nu vs Re
180
160
140
120
100
80 Nu vs Re
60
40
20
0
0 5000 10000 15000 20000 25000
f vs Re
0.054
0.053
0.052
0.051
0.05
0.049
f vs Re
0.048
0.047
0.046
0.045
0.044
0 5000 10000 15000 20000 25000
180 0.054
160 0.053
140 0.052
0.051
120
0.05
100
0.049
80
0.048
60
0.047
40 0.046
20 0.045
0 0.044
0 5000 10000 15000 20000 25000
Nu variation friction factor variation
30 | P a g e
" P in-P out vs Re"
4500
4000
3500
3000
2500
2000 " P in-P out vs Re"
1500
1000
500
0
0 5000 10000 15000 20000 25000
w vs Re
0.0008
0.0007
0.0006
0.0005
0.0004
w vs Re
0.0003
0.0002
0.0001
0
0 5000 10000 15000 20000 25000
The Optimum heat transfer rate and minimum power loss or dissipation of mechanical energy at
D/d =10 occurs corresponding to Reynolds no. 12250
31 | P a g e
4.2. D/d = 15 ( d = 10 mm )
So, for Turbulent flow, v > 0.817 m/s. So, the velocity of cold fluid is taken 2.5 m/s and velocity
of hot fluid is taken from 1 m/s, increased in steps of 0.2 m/s to 2 m/s.
Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )
Tf vs Nodes
337.2
337.1
337
336.9
336.8
336.7
0 10000 20000 30000 40000 50000 60000 70000 80000 90000
32 | P a g e
For Divisions (30, 32, 60, 62), it is Grid Independent, The graphs corresponding to it:
Nu vs Re
180
160
140
120
100
80 Nu vs Re
60
40
20
0
0 5000 10000 15000 20000 25000
f vs Re
0.046
0.045
0.044
0.043
0.042
f vs Re
0.041
0.04
0.039
0.038
0 5000 10000 15000 20000 25000
Chart Title
180 0.046
160 0.045
140 0.044
120
0.043
100
0.042
80
0.041
60
40 0.04
20 0.039
0 0.038
0 5000 10000 15000 20000 25000
33 | P a g e
P in - P out vs Re
6000
5000
4000
3000
P in - P out vs Re
2000
1000
0
0 5000 10000 15000 20000 25000
w vs Re
0.001
0.0009
0.0008
0.0007
0.0006
0.0005
0.0004
0.0003
0.0002
0.0001
0
0 5000 10000 15000 20000 25000
w vs Re
The Optimum heat transfer rate and minimum power loss or dissipation of mechanical energy at
D/d =15 occurs corresponding to Reynolds no. 12750
34 | P a g e
4.3. D/d = 20 ( d = 10 mm )
So, for Turbulent flow, v > 0.74734 m/s. So, the velocity of cold fluid is taken 2.5 m/s and
velocity of hot fluid is taken from 1 m/s, increased in steps of 0.2 m/s to 2 m/s.
Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )
Tf vs Nodes
336
335.8
335.6
335.4
335.2
335
334.8
334.6
0 20000 40000 60000 80000 100000 120000
35 | P a g e
For Divisions (30, 32, 60, 62), it is Grid Independent, The graphs corresponding to it:
Nu vs Re
180
160
140
120
100
80 Nu vs Re
60
40
20
0
0 5000 10000 15000 20000 25000
f vs Re
0.045
0.04
0.035
0.03
0.025
0.02 f vs Re
0.015
0.01
0.005
0
0 5000 10000 15000 20000 25000
180 0.045
160 0.04
140 0.035
120 0.03
100 0.025
80 0.02
60 0.015
40 0.01
20 0.005
0 0
0 5000 10000 15000 20000 25000
36 | P a g e
P in - P out vs Re
7000
6000
5000
4000
3000 P in - P out vs Re
2000
1000
0
0 5000 10000 15000 20000 25000
w vs Re
0.0012
0.001
0.0008
0.0006
w vs Re
0.0004
0.0002
0
0 5000 10000 15000 20000 25000
The Optimum heat transfer rate and minimum power loss or dissipation of mechanical energy at
D/d =20 occurs corresponding to Reynolds no. 10000
37 | P a g e
4.4. D/d = 25 ( d = 10 mm )
So, for Turbulent flow, v > 0.698 m/s. So, the velocity of cold fluid is taken 2.5 m/s and velocity
of hot fluid is taken from 1 m/s, increased in steps of 0.2 m/s to 2 m/s.
Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )
Tf vs Nodes
335
334.5
334
333.5
333
332.5
0 20000 40000 60000 80000 100000 120000 140000
38 | P a g e
For Divisions (32, 34, 68, 70), it is Grid Independent, The graphs corresponding to it:
f vs Re
0.04
0.039
0.038
0.037
0.036
f vs Re
0.035
0.034
0.033
0.032
0 5000 10000 15000 20000 25000
Nu vs Re
180
160
140
120
100
80 Nu vs Re
60
40
20
0
0 5000 10000 15000 20000 25000
Chart Title
180 0.04
160 0.039
140 0.038
120
0.037
100
0.036
80
0.035
60
40 0.034
20 0.033
0 0.032
0 5000 10000 15000 20000 25000
39 | P a g e
P in - P out vs Re
9000
8000
7000
6000
5000
4000 P in - P out vs Re
3000
2000
1000
0
0 5000 10000 15000 20000 25000
w vs Re
0.0014
0.0012
0.001
0.0008
0.0006 w vs Re
0.0004
0.0002
0
0 5000 10000 15000 20000 25000
The Optimum heat transfer rate and minimum power loss or dissipation of mechanical energy at
D/d =25 occurs corresponding to Reynolds no. 12500
40 | P a g e
4.5. D/d = 30 ( d = 10 mm )
So, for Turbulent flow, v > 0.6612 m/s. So, the velocity of cold fluid is taken 2.5 m/s and velocity
of hot fluid is taken from 1 m/s, increased in steps of 0.2 m/s to 2 m/s.
Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )
Tf vs Nodes
333.6
333.4
333.2
333
332.8
332.6
332.4
0 20000 40000 60000 80000 100000 120000 140000 160000 180000
41 | P a g e
For Divisions (32, 34, 68, 70), it is Grid Independent, The graphs corresponding to it:
f vs Re
0.038
0.037
0.036
0.035
0.034 f vs Re
0.033
0.032
0.031
0 5000 10000 15000 20000 25000
Nu vs Re
180
160
140
120
100
80 Nu vs Re
60
40
20
0
0 5000 10000 15000 20000 25000
180 0.038
160
0.037
140
0.036
120
100 0.035
80 0.034
60
0.033
40
0.032
20
0 0.031
0 5000 10000 15000 20000 25000
42 | P a g e
P in - P out vs Re
10000
8000
6000
4000 P in - P out vs Re
2000
0
0 5000 10000 15000 20000 25000
w vs Re
0.0016
0.0014
0.0012
0.001
0.0008
w vs Re
0.0006
0.0004
0.0002
0
0 5000 10000 15000 20000 25000
The Optimum heat transfer rate and minimum power loss or dissipation of mechanical energy at
D/d =30 occurs corresponding to Reynolds no. 13000
43 | P a g e
CHAPTER- 5
5. CONCLUSION
44 | P a g e
6. REFERENCS
1. Experimental and CFD estimation of heat transfer in helically coiled heat enchangers – S.D.
Sancheti, Dr. P.R. Suresh, 2012.
2. Development of heat transfer coefficient for concentric helical coil heat exchanger – Rahul
Kharat, Ntin Bhardwaj, R.S. Jha. International Journal of Thermal Sciences 48, 2009 (2300-
2308).
3. Thermal hydraulic characteristics of air-water two-phase flows in helical pipes. - J.S. Jayakumar ,
S.M Mahajani, J.C Mandal, Kannan N. Iyer and P.K. Vijayan , 2010 (501-512).
4. Experimental study of mixed convection of heat exchanger in vertically concentric helical coil
heat exchanger - N. Ghorbani, H. Teherain, M. Gorhi and H. Mirgolbabaei, 2010 , Experimental
thermal and Fluid science 34 (900-905).
5. Numerical studies of a double pipe helical heat exchanger - Timothy J. Rennie, Vijaya G.S.
Raghavan, Applied thermal engineering 26, 2006 (1266-1273).
6. Boundary vorticity method for laminar forced convection heat transfer in curved pipes - Akiyama,
M. and Cheng, K.C., 14(10): (1659–1675).
7. Thermodynamic optimization of a coiled tube heat exchanger under constant wall heat flux
condition – Ashok k. Satapathy - Science direct (1122-1126 )
8. Tata McGraw-Hill Education, 1978 – Heat and mass transfer – P.K. Nag, 3rd edition.
9. Tata McGraw-Hill Education, 1989 – Heat transfer – Yunus A. Cengel, 3rd edition.
45 | P a g e