Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 163839, 7 pages

http://dx.doi.org/10.1155/2013/163839

Research Article

CFD Simulation of Heat Transfer and Friction Factor Augmentation in a Circular Tube Fitted with Elliptic-Cut Twisted Tape Inserts

Sami D. Salman,1,2 Abdul Amir H. Kadhum,¹ Mohd S. Takriff,¹ and Abu Bakar Mohamad¹

¹Department of Chemical and Process Engineering, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia

²Biochemical Engineering Department, Al-Khwarizmi College of Engineering, University of Baghdad, Baghdad 47024, Iraq

Received 30 April 2013; Accepted 10 June 2013

Academic Editor: Tirivanhu Chinyoka

Copyright © 2013 Sami D. Salman et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents the application of a mathematical model for simulation of the swirling flow in a tube induced by elliptic-cut and classical twist tape inserts. Effects of the twist ratio (.93, 3.91, and 4.89) and cut depth (, 0.8, and 1.4 cm) on heat transfer enhancement (Nu) and friction factor (f) in laminar flow are numerically investigated. The simulation is carried out using commercial CFD package (FLUENT-6.3.26) to grasp the physical behaviour of the thermal and fluid flows of a constant heat-fluxed tube fitted with elliptic-cut twist tape in the laminar flow regime for the Reynolds number ranging from 200 to 2100. The simulated results matched the literature correlations of plain tube for validation with 8% variation for Nusselt number and 10% for friction factor. The results show that the heat transfer rate and friction factor in the tube equipped with elliptic-cut twist tape (ECT) are significantly higher than those fitted with classical twist tape (CTT). Moreover the results also reveal that the Nusselt number and the friction factor in the tube with elliptic-cut twisted tape (ECT) increase with decreasing twist ratios (y) and cut depths (w).

1. Introduction

The heat transfer augmentation techniques are widely utilized in many applications in the heating process to enable reduction in weight and size or increase the performance of heat exchangers. These techniques are classified as active and passive techniques. The active technique required external power such as surface vibration and electric or acoustic fields, whereas the passive techniques required fluid additives, special surface geometries, or swirl/vortex flow devices, that is, twisted tape inserts. The passive techniques are advantageous compared with the active techniques because the swirl inserts manufacturing process is simple and can be easily employed in an existing heat exchanger. Moreover the passive techniques can play an important role in the heat transfer augmentation if a proper configuration of the insert is being selected depending on working conditions that have been reported in the literature [1–9]. Due to advances in computer software, the Computational Fluid Dynamics (CFD) modelling technique was developed as a powerful and effective tool for more understanding the hydrodynamics of heat transfer when using twist tape inserts.

Zhang et al. [10] investigated the heat transfer characteristics of a helically baffled heat exchanger combined with a finned tube using CFD modeling. Sivashanmugam et al. [11] reported the modeling of heat transfer augmentation in circular tube fitted with helical twist insert in a laminar and turbulent flows using CFD.

Eiamsa-ard et al. [12] have conducted a numerical study on a tube equipped with loose-fit twisted tapes using four turbulence models. The results showed that the prediction obtained by SST k-ω turbulence model has a better agreement with measurement results compared with other models.

Nagrajam and Sivashanmugam [13] conducted CFD simulations of heat transfer characteristics of Al₂O₃ nanofluid in a circular tube fitted with helical twist inserts under constant heat flux using Fluent version 6.3.26. Different concentrations of Al₂O₃ nanoparticles (0.5%, 1.0%, and 1.5%) and twist tape inserts with different twist ratios ( = 2.93, 3.91, and 4.89) have been used for the simulation. The results showed that the Nusselt number and the friction factor are increased with increasing the nanofluid concentration and decreasing the twist ratios ().

Pathipakka and Sivashanmugam [14] proposed CFD simulation of heat transfer and friction factor behavior for the circular tube fitted with right-left helical twist insert with 100 mm spacer. The simulated Nusselt number and friction factor werecompared with the experimental data and observed to have good agreement.

Shabanian et al. [15] have carried out an experimental and CFD modeling on heat transfer and friction factor characteristics in air cooled heat exchanger using butterfly twist tape insert. They found that the insert configuration has a main effect on the Nusselt number, friction factor, and thermal performance factor and the maximum thermal performance factor was obtained by the configuration of butterfly insert with an inclined angle of 90°. The results verified that there was a good agreement between the predicted and measured values of Nusselt number and friction factor values.

In the present study, a new configuration of twist tape inserts is presented using CFD simulation to predict the Nusselt number and friction factor in laminar flow regimes based on experimental data listed in [14]. This study can be used as guideline for experimental works.

2. Physical Model

The geometry of the elliptic-cut twisted tape (ECT) insert is illustrated in Figure 1. Twist tape with thickness () of 0.08 cm and width () of 2.45 cm and relative twisted ratios ( = 2.93, 3.91, and 4.89) fits in a tube with a diameter () of 2.54 cm and length () of 180 cm. Different cut depth ( = 0.4, 0.8, and 1.4 cm) is used for twisted tape with twist ratio . Steel and aluminium were selected as the material of construction of the tube and twisted tape. Water was selected as the working fluid, and the thermophysical properties were assumed to be temperature independent. The thermophysical properties of water and materials used for simulation are listed in Table 1.

Table 1: Thermophysical properties of materials at 298 K.

Materials
Density (Kg/m³)
Specific heat (J/kg K)
Thermal conductivity (W/m K)
Viscosity (Pa s)

Water
998.2
4182
0.6
0.001003

Steel
8030
502.48
16.27
—

Aluminium
2719
871
202.4
—

Figure 1: Elliptic-cut twisted insert.

The Reynolds number (), the Nusselt number (), and the friction factor () are defined by the following equations:

where is the density, is velocity, is dynamic viscosity, and is the thermal conductivity of the fluid. is the heat transfer coefficient, and is the inner diameter of the tube.

3. Numerical Simulations

Three-dimensional numerical simulation of the conjugate heat transfer was conducted using the CFD code FLUENT 6.3.26. The modeling was carried out in order to predict and explain the effect of elliptic-cut configuration on the Nusselt number and friction factor at steady-state laminar flow for constant heat-fluxed tube. The CFD modeling involves numerical solutions of the conservation equations for mass, momentum, and energy. These equations for incompressible flows can be written as follows.

3.1. Continuity Equation for an Incompressible Fluid

It is given as follows:

3.2. Conservation of Momentum

 It is given as follows:

3.3. Conservation of Energy

It is given as follows:

4. Geometry and Grid Arrangement

The geometry and the gird of plain tube and elliptic-cut twisted tape were generated in GAMBIT and the CFD code FLUENT was used for simulation. The geometry consists of a cylindrical tube of diameter 25.54 cm and length of 180 cm. Figure 2 shows the grid for the plain tube configuration using a quadrilateral face mesh over the volume of the cylinder. Different types of meshing twisted tapes are available to mesh the volume but tetrahedral/hybrid and T grid type elements were the best mesh for irregular shapes. The grid generated in the tube fitted with ECT inserts is shown in Figure 3. The boundary conditions for the geometry were defined for inlet, outlet, wall of twist tape and cylinder wall, the continuum volume of fluid was defined as water. The mesh file is exported to FLUENT successfully.

Figure 2: Grid for the plain tube.

Figure 3: Grid for the plain tube with elliptic-cut twisted tape (ECT) insert.

5. Modeling Parameters

Numerical values of experimental data mentioned in [14] were used in a number of simulations thet are given in Table 2.

Table 2: Numerical values of the parameters used for simulations [14].

Flow rate (kg/sec)
0.003
0.005
0.006
0.008
0.01
0.116
0.0133
0.015
0.0166
0.02
0.03

Heat flux (w/m²)
240.03
340.77
459.15
563.38
749.85
1001.57
1363.26
1512.68
1893.56
2130.05
2445.25

6. Numerical Method

The commercial CFD package of Fluent 6.3.26 software was chosen as the CFD tool to solve the Navier-Stokes equation in common with energy equation companied with boundary conditions. Solution sequential algorithm (segregated solver algorithm) has been used with the following settings; implicit formulation, steady (time-independent) calculation, viscous laminar model, and energy equation.

7. Results and Discussion

7.1. Validation of Plain Tube Simulation Results

The simulated results of Nusselt number and friction factor for plain tube were compared with the correlations developed by Sieder and Tate [16] for validation. The predicted values of Nusselt number and friction factor are demonstrated in Figures 4 and 5; apparently, the present results reasonably agree well with the available correlations with 8% variation for Nusselt number and 20% for friction factor.

Figure 4: Plain tube simulated Nusselt number versus literature data.

Figure 5: Plain tube simulated friction factor versus literature data.

7.2. Effect of Twist Ratio on Heat Transfer and Friction Factor

The simulated data of the Nusselt number and friction factor and their variation with a Reynolds number of cut twisted tape inserts with twist ratios (, 3.91, and 4.89) are shown in Figures 6 and 7. Figure 6 indicates that the Nusselt number increases with Reynolds number increasing and the heat transfer rate is higher for the twist tape set than for the plain tube because of strong swirl flow in the presence of the twist tape. It is found that the heat transfer rate with the twist ratio () is higher than those with other ratios ( and 4.89); this means that the turbulent intensity obtained from the lower twist ratio is higher than those from higher ratios (). Figure 7 shows the variation of friction factor with a Reynolds number for different twist ratios ( = 2.93, 3.91, and 4.89). The friction factor obtained from the tube with twisted tape insert is significantly higher than the plain tube. Moreover, the use of smaller twist ratio leads to higher tangential contact between the swirling flow and the tube surface. Therefore, the twisted tape with twist ratio () has a maximum friction factor.

Figure 6: Simulated Nusselt number for plain tube and twisted tape with different twist ratios ().

Figure 7: Simulated friction factor for plain tube twisted tape with different twist ratios ().

7.3. Effect of Cut Depth on Heat Transfer and Friction Factor

The Nusselt number and friction factor and their variation with a Reynolds number of ECT inserts with twist ratio = 2.93 and cut depth ( = 0.4, 0.8, and 1.4 cm) are shown in Figures 8, 9, 10, and 11. For a given Reynolds number, Nusselt number and friction factor are increased with decreasing the cut depth; this is mainly due to the combined effect of common swirling flow by the twisted tape and turbulence generated by the alternative cuts along the edge of the twisted tape. Subsequently this leads to the destruction of the thermal boundary layer and creating better flow mixing between the fluids at the core and heating wall surface.

Figure 8: Simulated Nusselt number for elliptic-cut twisted tape with and different cut depth ().

Figure 9: Simulated friction factor for elliptic-cut twisted tape with and different cut depth ().

Figure 10: Simulated Nusselt number for elliptic-cut twisted tape (ECT) and classical twist tape (CTT).

Figure 11: Simulated friction factor for elliptic-cut twisted tape (ECT) and classical twist tape (CTT).

7.4. Velocity Field

Figures 12 and 13 show the velocity field of classical and Elliptic-cut twisted tape inserts with a twisted ratio () and cut depth () for . The figures illustrate that the predicted velocity in tube fitted with elliptic-cut twisted tape insert is higher than the classical twist tape insert and the vortices generated by the alternative cuts along the edge of the twisted tape are higher than those of classical twist tape. This can be the reason for more heat transfer rate obtained by the Elliptic-cut twisted tape as compared to classical twist tape insert.

Figure 12: The velocity field of classical twist tape in an axial slice inside the tube.

Figure 13: The velocity field of elliptic-cut twist tape in an axial slice inside the tube.

8. Conclusion

In the present study, a circular tube inserted by Elliptic-cut twisted tape and classical twist tape insert with twisted ratios ( = 2.93, 3.91, and 4.89) and cut depths ( = 0.4, 0.8, and 1.4 cm) in laminar flow conditions has been simulated using fluent version 6.23.26. The data obtained by simulation match the literature value for plain tube with a maximum discrepancy of 8% for Nusselt number and 10% for friction factor. The simulated results show that the elliptic-cut twisted tape with twist ratio () and cut depth ( cm) offered higher heat transfer rate and friction factor compared to the plain tube and other twisted tape. Furthermore, the influence of the cut depth ( cm) was more dominant than that of the cut depths ( = 0.8 and 1.4 cm) for all the Reynolds number.

Nomenclature

: | Energy component in energy equation
: | Force component in momentum equation, N
: | Fanning friction factor
: | Acceleration due to gravity, m/s²
: | Thermal conductivity in energy equation, W/m K
: | Mass flow rate of fluid, kg/s
: | Reynolds number based on internal diameter of the tube, dimensionless
: | Nusselt number, dimensionless
: | Pressure component in momentum equation, N/m²
: | Accumulation of mass, Kg
: | Accumulation of energy, J
: | Temperature,  °C
: | Velocity component in momentum equation, m/s
: | Twist ratio (length of one twist (360°)/diameter of the twist), dimensionless.

Greek Symbols

: | Density component in governing equations
: | Stress component in momentum equation, N/m².

Acknowledgments

The authors would like to thank Universiti Kebangsaan Malaysia and MOSTI of Grant Research (FRGS/1/2013/TK07/UKM/01/1) for financial support.

References

L. Wang and B. Sundén, “Performance comparison of some tube inserts,” International Communications in Heat and Mass Transfer, vol. 29, no. 1, pp. 45–56, 2002.

A. Ujhidy, J. Németh, and J. Szépvölgyi, “Fluid flow in tubes with helical elements,” Chemical Engineering and Processing, vol. 42, no. 1, pp. 1–7, 2003.

S. Jaisankar, T. K. Radhakrishnan, and K. N. Sheeba, “Experimental studies on heat transfer and friction factor characteristics of thermosyphon solar water heater system fitted with spacer at the trailing edge of twisted tapes,” Applied Thermal Engineering, vol. 29, no. 5-6, pp. 1224–1231, 2009.

J. Guo, A. Fan, X. Zhang, and W. Liu, “A numerical study on heat transfer and friction factor characteristics of laminar flow in a circular tube fitted with center-cleared twisted tape,” International Journal of Thermal Sciences, vol. 50, no. 7, pp. 1263–1270, 2011.

K. Wongcharee and S. Eiamsa-ard, “Friction and heat transfer characteristics of laminar swirl flow through the round tubes inserted with alternate clockwise and counter-clockwise twisted-tapes,” International Communications in Heat and Mass Transfer, vol. 38, no. 3, pp. 348–352, 2011.

E. Z. Ibrahim, “Augmentation of laminar flow and heat transfer in flat tubes by means of helical screw-tape inserts,” Energy Conversion and Management, vol. 52, no. 1, pp. 250–257, 2011.

S. K. Saha, U. N. Gaitonde, and A. W. Date, “Heat transfer and pressure drop characteristics of turbulent flow in a circular tube fitted with regularly spaced twisted-tape elements,” Experimental Thermal and Fluid Science, vol. 3, no. 6, pp. 632–640, 1990.

R. M. Manglik and A. E. Bergles, “Heat transfer and pressure drop correlations for twisted-tape inserts in isothermal tubes—part II: transition and turbulent flows,” Journal of Heat Transfer, vol. 115, no. 4, pp. 890–896, 1993.

V. K. Naumov and N. N. Semashko, “Analytical model of estimating thermo-physical and strength parameters of cooled pipe with the twisted tape under asymmetric heating by a pulse of external heat flux,” Plasma Devices and Operations, vol. 3, no. 4, pp. 267–280, 1994.

Z. Zhang, D. Ma, X. Fang, and X. Gao, “Experimental and numerical heat transfer in a helically baffled heat exchanger combined with one three-dimensional finned tube,” Chemical Engineering and Processing, vol. 47, no. 9-10, pp. 1738–1743, 2008.

P. Sivashanmugam, P. K. Nagarajan, and S. Suresh, “Experimental studies on heat transfer and friction factor characteristics of turbulent flow through a circular tube fitted with right and left helical screw-tape insert,” Chemical Engineering Communications, vol. 195, no. 8, pp. 977–987, 2008.

S. Eiamsa-ard, K. Wongcharee, and S. Sripattanapipat, “3-D numerical simulation of swirling flow and convective heat transfer in a circular tube induced by means of loose-fit twisted tapes,” International Communications in Heat and Mass Transfer, vol. 36, no. 9, pp. 947–955, 2009.

P. K. Nagrajam and P. Sivashanmugam, “CFD Simulation of heat transfer augmentation in circular tube filled with right-left helical inserts with spacer,” International Journal of Chemical Engineering Research, vol. 6, no. 1, pp. 1–11, 2009.

G. Pathipakka and P. Sivashanmugam, “Heat transfer behaviour of nanofluids in a uniformly heated circular tube fitted with helical inserts in laminar flow,” Superlattices and Microstructures, vol. 47, no. 2, pp. 349–360, 2010.

S. R. Shabanian, M. Rahimi, M. Shahhosseini, and A. A. Alsairafi, “CFD and experimental studies on heat transfer enhancement in an air cooler equipped with different tube inserts,” International Communications in Heat and Mass Transfer, vol. 38, no. 3, pp. 383–390, 2011.

E. N. Sieder and C. E. Tate, “Heat transfer and pressure drop of liquids in tubes,” Industrial & Engineering Chemistry, vol. 28, no. 12, pp. 1429–1435, 1936.