Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Optimization of Heat Transfer Using CFD Simulation For Concentric Helical Coil Heat Exchanger For Constant Temperature Outer Wall

Download as pdf or txt
Download as pdf or txt
You are on page 1of 45

Optimization of heat transfer using CFD simulation

for concentric helical coil heat exchanger


for constant temperature outer wall

A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR


THE DEGREE OF

BACHELOR OF TECHNOLOGY
IN
MECHANICAL ENGINEERING
BY
SAGAR DAS
110ME0287

DEPARTMENT OF MECHANICAL ENGINEERING

NATIONAL INSTITUTE OF TECHNOLOGY

ROURKELA – 769008

(2013-2014)

1|Page
Optimization of heat transfer using CFD simulation
for concentric helical coil heat exchanger
for constant temperature outer wall

A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR


THE DEGREE OF

BACHELOR OF TECHNOLOGY
IN
MECHANICAL ENGINEERING
BY
SAGAR DAS
110ME0287

Under the guidance of


Prof. ASHOK K. SATAPATHY

DEPARTMENT OF MECHANICAL ENGINEERING

NATIONAL INSTITUTE OF TECHNOLOGY

ROURKELA – 769008 (2013-2014)

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.

Prof. Ashok K. Satapathy


Dept. of Mechanical Engineering
National Institute of Technology
Rourkela – 769008

3|Page
ACKNOWLEDGEMENT

I wish to express my profound gratitude and indebtedness to Prof. Ashok K. Satapathy,


Department of Mechanical Engineering, NIT, Rourkela for introducing the present topic and for
his inspiring guidance, endless support and valuable suggestion throughout the project work. I
am extremely fortunate to be involved in an exciting and challenging research project like
“Optimization of heat transfer using CFD simulation for concentric helical coil heat
exchanger for constant temperature outer wall”. It has enriched my life, giving me an
opportunity to work in a new environment of Fluent. This project increased my thinking and
understanding capability as I started the project from scratch. He has not only been a wonderful
supervisor but also a genuine person. I consider myself extremely lucky to be able to work under
the guidance of such a dynamic personality. Actually he is one of such genuine person for whom
my words will not be enough to express.
I am also thankful to Dr. K. P. Maity, H.O.D of Department of Mechanical Engineering,
National Institute of Technology, Rourkela for his constant support and encouragement.
Last, but not the least I extend my sincere thanks to all faculty members of Mechanical
Engineering department for making my project a successful one, for their valuable advice in
every stage and also giving me absolute working environment where I unlashed my potential. I
would like to thank all whose direct and indirect support helped me completing my thesis in
time. I want to convey my heartiest gratitude to my parents for their unfathomable
encouragement.

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:

Table No. Caption Page No.


1 Names of parts of the body 21
2 Input and Output Variables 24
3 Boundary Conditions 25
4 Residuals Variable 26

6|Page
ABSTRACT

Thermodynamic Optimization in heat transfer of a concentric coiled tube-in-tube heat exchanger


under constant wall temperature condition, based on Fluid–Fluid heat transfer is focused in this
paper. The parameters which influence the nature of flow in a helical coil are pitch coil diameter,
pitch and tube diameter in helical coils. 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. Thus, Nusselt number and friction factor were
also calculated. Graphs were plotted between Nusselt number, friction factor, pressure drop and
power loss with Reynolds number. The point where the friction fraction intersects with the
Nusselt number is the point where the heat transfer is optimum, corresponding to that Reynolds
number. Beyond that Reynolds number, the friction factor decreases rapidly, hence the pressure
drop increases and so the power loss also increases. Various velocity and temperature contours
were also obtained. Hence, we found the optimum value of Reynolds number for the
corresponding tube-to-coil diameter ratios. It thus minimizes the degradation of thermal energy
and viscous dissipation of mechanical energy.

7|Page
NOMENCLATURE

 f = Darcy’s friction factor


 m. = mass flow rate

 A = area of heat transfer (m2), A = (π*d2)/4, m2

 D = helical Coil diameter


 h = heat transfer coefficient (W m-2 K-1)
 H = helical pitch (m)
 K’ = thermal conductivity (W m-1 K-1)
 L= length of the pipe (m)
 Nu = Nusselt number
 Pr = Prandtl number
 q = heat transferred (W)
 R i = radius of the inner tube (m)
 R out = radius of the outer tube (m)
 R = resistance the flow of thermal energy (W-2 m2 K)
 Rc = pitch circle radius of the pipe (m)
 Re= Reynolds number
 v = velocity (ms-1)
 U = overall heat transfer coefficient (Wm-2 K-1)
 V = volume (m3)
 O = helix angle (rad)
 δ = curvature ratio = D/d
 Δ = temperature difference (K)
 μ = viscosity (kg m-1 s-1)
 ρ = density (kg m-3)
 LM = logarithmic mean
 w = power loss

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.

1.2. Types of Heat Exchangers [8]

1. Recuperators or Transfer type heat exchanger


2. Regenerators or Storage type heat exchanger
3. Mixers or Direct contact type heat exchanger

1.3. Flow arrangements in Recuperative heat exchangers [8]

1. Parallel flow heat exchanger


2. Counter flow heat exchanger
3. Cross flow heat exchanger

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:

1. Shell and tube heat exchanger


2. Double pipe heat exchanger
3. Spiral tube heat exchanger

1.5. Helical tube-in-tube Heat exchanger:

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.

Fig. Geometry of a helical coil

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:

Advantages of helical coil:

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:

1. Design of helical tube-in-tube heat exchanger is complex.


2. Cleaning of the tubes are more difficult than jackets and shells.
3. For high reactive or corrosive fluids, coils can’t be used, instead jackets are preferred.
4. Coils generally plays a major role in the selection of any agitation system, as densely
packed coils may create some unmixed regions by any interference with the fluid flow.

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:

Q = (T1 – T0) / ( R ins – R conv )

= (T1 – T0) / [(ln (R2/R1) / 2* 3.14*L k) + 1/(h*3.14*R2L)]

We know, that the critical radius (Rcr) of insulation depends on:


1. The thermal conductivity of insulation ‘k’ and
2. Convection heat transfer coefficient ‘h’.
Rate of heat transfer increases from the cylindrical surface with the increase of
insulation for R2 < Rcr, which reaches a maximum value when r2 = rcr, and then starts to
decrease from r2 >rcr. It is thus clear that the value of critical radius Rcr , will be large when k is
high and h is less.

So, for a cylindrical surface the following relation holds:

Rcr, max = kmax, insulation / hmin

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.

To optimize the heat transfer coefficient by plotting graphs between:

a) Friction factor and Reynolds number


b) Nusselt number and Reynolds number
c) Decrease in pressure and Reynolds number
d) Power loss and Reynolds number
e) Temperature profile and velocity profile

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]

N. Ghorbani, H. Teherain, M. Gorhi and H. Mirgolbabaei worked on the


Experimental study of mixed convection of heat exchanger in vertically concentric helical coil
heat exchanger. Mixed convection heat transfer in a coil in cell heat exchanger with varying
Reynolds no. with varying tube-to-coil diameter ratios and coil pitch were investigated
experimentally, for both laminar and turbulent flow. The effects of the coil pitch and diameter of
tube on shell-side heat transfer coefficient were studied of the helically coil heat exchanger.
Nusselt number correlation with variable coil parameters were analyzed to best fit the data. It is
compared with similar studies with specific boundary conditions. It, concluded that, the tube
diameter has negligible influence on heat transfer on shell- side, coil surface has –ve effect on h0.
The overall heat transfer increases with h0. [4]
Timothy J. Rennie, Vijaya G.S. Raghavan worked on the Numerical studies of
a double pipe helical heat exchanger. A double-pipe helical heat exchanger was modeled
numerically for numerical flow and the heat transfer characteristics were studied for different
fluid flow rates and tube sizes. The overall heat transfer coefficient for both counter and parallel
flow were calculated. Simulations were validated by comparing the Nusselt numbers. Greatest
thermal resistance were found around the annular region. A correlation was found between the
annulus Nusselt no. with a modified dean noumber, which gave a strong linear relationship. [5]
Akiyama, M. and Cheng, K.C. worked on the Boundary vorticity method for
laminar forced convection heat transfer in curved pipes. Based on their theoretical work he
proposed that the Vortex formation and conversion from laminar to turbulent flow occurs due to
the secondary flow by the curvature of the helically coil, which produces a centrifugal force and
was analyzed by the CFD using Fluent and various correlations were developed. [6]

Dr. Ashok K. Satapathy worked on the Thermodynamic Optimization of a


coiled tube heat exchanger under constant wall heat flux condition, based on Fluid – Fluid heat
transfer. Optimum diameter ratio of coil, minimizes the degradation of thermal energy and
viscous dissipation of mechanical energy. [7]

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.

3.1.3. Boolean Operation

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:

Table 1. Names of parts of the body

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.

3.2.1. Mapped Face Meshing

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.

3.2.3. Named Selection

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:

3.3.1. Problem Setup

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.

Table 2. Input and Output Variables:

(i) Thermodynamic Properties

Constant Notation Units


Outlet temperature Tout K
Outlet pressure Pout Pa
Outlet density ρ kg/m3
Mass flow rate M kg/s

(ii) Fluid Properties

Property Notation Units


Specific heat Cp -
viscosity μ -

(iii) Geometric Inputs

Component Dimension Notation Unit


Helical Tube Inner radius RI mm
Thickness radius Rth mm
Outer radius Rout mm

24 | P a g e
No. of Coils n -
Pitch of helix P mm
Coil dia D mm

(iv) Input fluid Conditions

Component Constant Notation Unit


Heat Hot fluid Inlet velocity vin, hot m/s
exchanger Cold fluid Inlet velocity vin, cold m/s
as double Cold fluid Inlet temperature Tin, cold Kelvin
helical tube Hot fluid Inlet temperature Tin, hot Kelvin

3.3.4. Cell zone conditions


The parts are then assigned as water-liquid and copper as per the fluid/solid parts.
3.3.5. Boundary Conditions

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.

Table 3. Boundary Conditions

Boundary Velocity Magnitude Turbulent Turbulent Temperature


Condition Type Kinetic Dissipation
Energy Rate
Inner fluid Inlet Velocity Inlet 1 m/s 0.01 m2/s2 0.1 m2/s3 343 K
Inner fluid Outlet Pressure Outlet - - - -
Outer fluid Inlet Velocity Inlet 2.5 m/s 0.01 m2/s2 0.1 m2/s3 293 K
Outer fluid Outlet Pressure Outlet - - - -

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:

 Density = 998.2 kg/m3


 Length 473.145 mm
 Temperature = 343 K
 Velocity = 1 m/s
 Viscosity = 0.001003 kg/m-s
 Ratio of specific heats = 1.4

3.3.7. Solution Methods

The solution methods are specified as follows:

 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

3.3.8. Solution Control and Initialization

Under the relaxation factors the parameters defined are

 Pressure = 0.3 Pascal


 Density = 1 kg/m3
 Body forces = 1 kg/m2 s2
 Momentum = 0.7 kg-m/s
 Turbulent kinetic energy = 0.8 m2/s2

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.

3.3.9. Measure of Convergence

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:

Table 4 Residuals Variable

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

3.3.10. Run Calculation

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

RESULTS AND ANALYSIS

28 | P a g e
4.1. D/d = 10 ( d = 10 mm )

For Turbulent Flow,


Re = 2300 [ 1+ 8.6 * ( d/D )0.45 ]

= 9318.2 ( for d/D = 1/10 )

Now, Putting, Re = density * v * d / viscosity

 9318.2 = 998.2 * v * 10-2 / 0.001003


 => v = 0.936 m/s

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.

Leff = n [ (3.14*D)2 + p2 ]0.5 = 473.145 m, ( where n = 1.5, p = 30 mm )

Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )

No. of divisions Nodes Elements Outlet Fluid Temp.


1 20, 24, 44, 48 32554 29904 338.26
2 24, 26, 50, 54 43598 106098 338.696
3 30, 32, 60, 62 64555 338407 338.894
4 40, 44, 80, 84 77314 357041 338.899

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 )

For Turbulent Flow,


Re = 2300 [ 1+ 8.6 * ( d/D )0.45 ]

= 8147.7 ( for d/D = 1/15 )

Now, Putting, Re = density * v * d / viscosity

 8147.7 = 998.2 * v * 10-2 / 0.001003


 v = 0.817 m/s

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.

Leff = n [ (3.14*D)2 + p2 ]0.5 = 708.38 m, ( where n = 1.5, p = 30 mm )

Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )

No. of divisions Nodes Elements Outlet Fluid Temp.


1 20, 24, 44, 48 40392 30415 336.768
2 24, 26, 50, 54 52798 116098 336.98
3 30, 32, 60, 62 70052 390200 337.1246
4 40, 44, 80, 84 81128 395345 337.13

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

Nu variation friction factor variation

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 )

For Turbulent Flow,


Re = 2300 [ 1+ 8.6 * ( d/D )0.45 ]

= 7437.63 ( for d/D = 1/20 )

Now, Putting, Re = density * v * d / viscosity

 7437.63 = 998.2 * v * 10-2 / 0.001003


 v = 0.74734 m/s

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.

Leff = n [ (3.14*D)2 + p2 ]0.5 = 943.074 m, ( where n = 1.5, p = 30 mm )

Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )

No. of divisions Nodes Elements Outlet Fluid Temp.


1 20, 22, 44, 48 63289 88512 334.66
2 24, 26, 50, 54 85738 179275 335.251
3 30, 32, 60, 62 92423 516630 335.753
4 40, 44, 80, 84 108405 578459 335.79

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

Nu variation friction factor variation

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 )

For Turbulent Flow,


Re = 2300 [ 1+ 8.6 * ( d/D )0.45 ]

= 6946.79 ( for d/D = 1/25 )

Now, Putting, Re = density * v * d / viscosity

 6946.79 = 998.2 * v * 10-2 / 0.001003


 v = 0.698 m/s

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.

Leff = n [ (3.14*D)2 + p2 ]0.5 = 1178.36 m, ( where n = 1.5, p = 30 mm )

Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )

No. of divisions Nodes Elements Outlet Fluid Temp.


1 20, 22, 44, 48 63579 90599 332.897
2 24, 26, 50, 54 99512 192873 333.217
3 32, 34, 68, 70 116043 649525 334.5029
4 40, 44, 80, 84 129849 694785 334.61

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

Nu variation friction factor variation

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 )

For Turbulent Flow,


Re = 2300 [ 1+ 8.6 * ( d/D )0.45 ]

= 6580.768 ( for d/D = 1/30 )

Now, Putting, Re = density * v * d / viscosity

 6580.768 = 998.2 * v * 10-2 / 0.001003


 v = 0.6612 m/s

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.

Leff = n [ (3.14*D)2 + p2 ]0.5 = 1414.437 m, ( where n = 1.5, p = 30 mm )

Grid Independent Test: ( v = 1 m/s : hot fluid, cold fluid : v= 2.5 m/s )

No. of divisions Nodes Elements Outlet Fluid Temp.


1 20, 22, 44, 48 74563 102920 332.482
2 24, 26, 50, 54 124575 116870 332.781
3 32, 34, 68, 70 134678 717085 333.341
4 40, 44, 80, 84 157484 984747 333.348

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

Nu variation friction factor variation

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

Thermodynamic Optimization in heat transfer of a concentric coiled tube-in-tube heat exchanger


under constant wall temperature condition, based on Fluid–Fluid heat transfer is discussed. 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. Thus, Nusselt number and friction factor were also calculated. Graphs were
plotted between Nusselt number, friction factor, pressure drop and power loss with Reynolds
number. The point where the friction fraction intersects with the Nusselt number is the point
where the heat transfer is optimum, corresponding to that Reynolds number. Beyond that
Reynolds number, the friction factor decreases rapidly, hence the pressure drop increases and so
the power loss also increases. Various velocity and temperature contours were also obtained.
Hence, we found the optimum value of Reynolds number for the corresponding tube-to-coil
diameter ratios. It thus minimizes the degradation of thermal energy and viscous dissipation of
mechanical energy.

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

You might also like