Abstract
This paper reports numerical investigations of heat transfer and friction factor characteristics in swirling flow conditions using CFD simulation. A commercial CFD package, FLUENT 6.3.26, was used in this study. 3D models for circular tube fitted with classical and quadrant-cut twisted tape (QCT) inserts with three twist ratios ( = 2.93, 3.91, and 4.89) and different cut depths ( = 0.5, 1.0, and 1.5 cm) were generated for the simulation. The data obtained from the CFD simulation were verified with the literature correlations of plain tube with the discrepancy of less than ±8% for Nusselt number and ±10% for friction factor. The results show that there was a significant increase in heat transfer coefficient and friction factor in the tube fitted with quadrant-cut twisted tape (QCT) with decreasing of twist ratio () and cut depth (). Furthermore, the configuration of QCT insert with a twist ratio of = 2.93 and a cut depth of = 0.5 cm offered higher heat transfer rate and friction factor than other twist ratios.
1. Introduction
The performance of heat exchangers especially for single-phase flows can be enhanced by many augmentation techniques. In general, heat transfer augmentation techniques can be divided into two categories: active techniques which require external power source and passive techniques which do not require external power source. Several experimental studies on heat transfer augmentation techniques using twisted tape as passive technique have been reported in the literature [1–15]. Thereafter and due to the advances in computer hardware and software and consequently increase in computation speed, the CFD modeling technique has been recognized as a powerful and effective tool to gain better understanding of the complex hydrodynamics in many industrial processes. Kumar et al. [16] studied the hydrodynamics and heat transfer characteristics of tube in a pilot scale helical heat exchanger fitted with semicircular plate in the annular area using a commercial CFD package to predict the flow and thermal properties in the heat exchanger tube. They found that the Nusselt number and friction factor in the inner and outer tubes from the CFD results coincided with the experimental data. Sivashanmugam et al. [17] reported the modeling of heat transfer augmentation in a circular tube fitted with a helical twist insert in a laminar and turbulent flow using CFD. Jayakumar et al. [18] offered a comparison study of CFD simulations and experiment on convective heat transfer in a double pipe helical heat exchanger and developed an empirical correlation for estimation of the inner heat transfer coefficient of a helical coil. Kharat et al. [19] developed heat transfer coefficient correlation for the flue gas side of concentric helical heat exchanger based on both the CFD simulation and experimental work. The correlation was used to analyze the data obtained from CFD and experimental results and to account for the effects of the gap between the concentric coil, tube diameter, and coil diameter on heat transfer coefficient. The predicted results of the correlations were in good agreement with the experimental results within an error band of 3-4%. Nagarajan and Sivashanmugam [20] reported the simulation of heat transfer augmentation and friction factor characteristics of a circular tube fitted with a right-left helical twist insert with a 100 mm spacer using CFD. The results of the simulated Nusselt number and friction factor for a given flow rate and twist ratio were compared with the experimental data and found a good agreement between the two. Shabanian et al. [21] conducted an experiment and (CFD) modeling on heat transfer, friction factor, and thermal performance of an air cooled heat exchanger equipped with three types of tube insert including butterfly, classic, and jagged twisted tape. They found that the predicted results in terms of turbulence intensity are in good agreement with measured values of Nusselt number and friction factor. Wang et al. [22] explicated the optimum configuration of regularly spaced short-length twisted tape in a circular tube by using computational fluid dynamics (CFD) modeling. Air was used as the test fluid with a turbulent flow with the configuration parameters including the free space ratio (), twist ratio (), and rotated angle (). The predicted results were in good agreement with the experimental data. The results showed that the larger rotation angle yields a higher heat transfer value and greater flow resistance, whereas the smaller twist ratio resulted in better heat transfer performance except for a larger rotation angle at high Reynolds number. Pathipakka and Sivashanmugam [23] performed a CFD simulation on the heat transfer characteristics of Al2O3 nanofluid in a circular tube fitted with helical twist inserts under constant heat flux using FLUENT version 6.3.26 in a laminar flow. The Al2O3 nanoparticles in water at different concentrations (0.5%, 1.0%, and 1.5%) and helical twist inserts with different twist ratios (, 3.91, and 4.89) were used in the simulation. The data obtained from the CFD simulation were compared with the reported values in the literature for water in a plain tube helical tape inserts. Salman et al. [24] report the application of a mathematical model of the heat transfer enhancement and friction factor characteristics of water in constant heat-fluxed tube fitted with elliptical cut twisted tape inserts using FLUENT version 6.3.26. Two types of swirl flow generator, classical and elliptical cut twisted tape with twist ratios (, 3.91, and 4.89) and different cut depths (, , and cm), were used in the CFD simulation. The results show that the enhancement of heat transfer rate and the friction factor induced by elliptical cut twisted tape inserts increases with Reynolds number and decreases with twist ratio. In addition, the results show that the elliptical cut twisted tape with twist ratio () and cut depth ( cm) offered higher heat transfer rate with significant increases in friction factor.
Salman et al. [25] numerically investigated the heat transfer of water in a uniformly heated circular tube fitted with V-cut twisted tape inserts in laminar flow using Fluent version 6.3.26. Classical and elliptical cut twisted tape with twist ratios (, 3.91, and 4.89) and different cut depths (, , and cm) was employed for the simulation. The results show that the V-cut twisted tape with twist ratio () and cut depth () generated the maximum heat transfer rate with significant increases in friction factor. In the present paper, the effects of quadrant-cut twist tape inserts with three twist ratios (, 3.91, and 4.89) and different cut depths (, 1.0, and cm) on heat transfer and friction factor characteristics in a constant heat-fluxed tube are numerically investigated. The numerical simulations were carried out by solving the mathematical model of swirling flow using a commercial CFD package (FLUENT 6.3.26). Results of interests such as Nusselt number, friction factor, and the twist tape configuration are reported.
2. Technical Details
The configuration of the quadrant-cut twisted tape (QCT) inserts is illustrated in Figure 1. An aluminium tape of 0.08 cm thickness and 2.45 cm width is uniformly winding over a length of 7.5, 10, and 12.5 cm to produce twist ratios of 2.93, 3.91, and 4.89, respectively. The twist ratio “” is defined as the ratio of the length of one full twist (360°) to the tape width. Three cut depths (, 1, and cm) were used for each twisted tape to produce a Parabolic-cut twisted tape. Steel tube with a diameter () of 2.54 cm and length () of 180 cm was used as test section and water was selected as the working fluid. The dimensions and the thermophysical properties of fluid and materials are summarized in Table 1.
3. Geometry and Grid Generation
GAMBIT, a preprocessor program in FLUENT-6.3.26 package, was used for generating the geometry and grids, whereas FLUENT version 6.3.26 was used in the module preprocessing. The geometry and the grid of the plain tube and quadrant-cut twisted tape inserts are shown in Figures 2 and 3. These geometries were created in GAMBIT and imported into FLUENT for the CFD simulation. The geometries of classical and Quadrant-cut twisted tape inserts with different twist ratios (, 3.91, and 4.89) and different cut depths (, 1, and cm) were made by winding a uniform strip of 25.54 mm width using the twist option in the sweeping of faces. The twist angle is 360° with lengths of 75, 100, and 125 mm for various twist ratios that are generated by using a perpendicular type of sweeping for the entire length of 1800 mm. The volume required for simulation is created by subtracting the twisted tape inserts from the plain tube geometry. The edged meshes were applied to each edge by using the particular interval count, whereas the front circular face is meshed by using tetrahedral and pave type meshing. The mesh face was swept over the entire volume using tetrahedral/hybrid elements and a T-grid type. The boundary conditions for the mesh volume are defined as inlet, outlet, wall, and type of fluid. Water was used as the fluid in the CFD simulation. Subsequently, the mesh file was exported to FLUENT for simulation.
4. Numerical Procedure
The commercial CFD package (FLUENT 6.3.26) was used to perform three-dimensional numerical calculations of the plain tube and twisted tape inserts in a constant heat-fluxed tube using the following governing equations.(1)Continuity equation for an incompressible fluid: (2)Conservation of momentum: (3)Conservation of energy:
5. Modeling Parameters
Numerical values of the mass flow rate and constant heat flux used for simulation are elaborated in Table 2. These data form the basis for the present numerical study. Steel and aluminium have been selected as the material of construction of the plain tube and twisted tapes, respectively. Water was used as the working fluids for all simulation cases.
6. Numerical Procedure
The commercial CFD package, Fluent 6.3.26, was used to solve the abovementioned governing equations based on the defined boundary conditions. The solution sequential algorithm (segregated solver algorithm) has been chosen for this study and solver setting includes implicit formulation, steady (time independent) calculation, viscous laminar model, and energy equation. The SIMPLE algorithm has been selected as the pressure-velocity coupling method and the first-order upwind scheme was used for the energy and momentum equations solution. The following equations are used to calculate the Nusselt number (Nu) and the friction factor () as a skin friction factor: where is the tube diameter, is the heat transfer coefficient, is the conductivity of water, is the heat flux on the tube, is the tube wall temperature, and is the bulk temperature of water , water inlet temperature, water outlet temperature, is the density, is dynamic viscosity, and is the water velocity.
7. Grid Dependence
A test on the grid dependence was conducted to evaluate the effects of grid sizes on the simulated results. In this study, the control volume was divided into a number of nodes (cell) and was tested for . Five mesh volumes were considered for plain tube and classical and quadrant-cut twist with twist ratio and cut depth 0.5 cm for the grid dependence test. It is observed that all the five mesh volumes have similar results of the Nusselt number with a small percentage of difference up to 0.5%. Hence, a domain with mesh volumes of 312624, 964970, and 985116 was chosen for the plain tube and classical and quadrant-cut twist tape, respectively, to reduce the computational time.
8. Results and Discussion
8.1. Model Validation
The numerical computation was first performed in a particular case of plain tube in order to validate the CFD simulation results with predicted data obtained by Sieder and Tate correlation [26]. The simulated Nusselt number and friction factor for the plain tube are presented in Figures 4 and 5, respectively. These figures show that the results of the present numerical study agree well with the available correlations within ±8% and ±10% for Nusselt number and friction factor, respectively.
8.2. Effect of Twist Ratio on Heat Transfer and Friction Factor
Figures 6 and 7 show the variation of the Nusselt number and friction factor with Reynolds number for the classical twisted tape (CTT) inserts with twist ratios (, 3.91, and 4.89), respectively. Figure 6 shows that the Nusselt number increases with increasing Reynolds numbers, and the heat transfer rate with the application of the twist tape is much higher than that of the plain tube. The higher heat transfer rate is most likely due to enhancement of mixing as a result of swirl flow generation. In addition, the twist tape with smaller twist ratio () provides higher heat transfer rate than others twist tapes ( and 4.89). This means that the smaller twist ratio produced higher turbulent intensity with longest flow path than those of larger twist ratios (). Figure 7 shows the variation of friction factor with Reynolds numbers for different twist ratios. It seems that the friction factor generated by twisted tape insert with smaller twist ratio is significantly higher than the plain tube and other twist ratios. Subsequently, the smallest twist ratio leads to higher tangential contact between the swirling flow and the tube surface.
8.3. Effect of Cut Depth on Heat Transfer and Friction Factor
Figures 8, 9, 10, and 11 present the variation of the Nusselt number and friction factor with Reynolds numbers for classical and quadrant-cut twisted tape inserts with twist ratio and different cut depths (, 1, and cm), respectively. These figures show that the Nusselt number and friction factor increased with decreasing cut depth at the same value of the Reynolds number. This phenomenon is due to the combined effects of the swirling flow and turbulence that occurred at smallest alternative cuts along the edge of the twisted tape resulting in higher destruction of the thermal boundary layer which leads to enhanced mixing between the fluid at the core and the heated surface.
8.4. Velocity Field
The velocity vectors predicted for the plain tube and the quadrant-cut twisted tape inserts with twist ratio () and different cut depths (, 1, and cm) are presented in Figure 12. This figure shows the longitudinal and radial vortices generated in the quadrant-cut twisted tape inserts. These vortices play a critical role for disturbing the boundary layer and uniform the temperature in the core flow. Consequently, enhanced mixing by extension heat transfer rate between the fluid at the core and the heated surface can be achieved.
8.5. Temperature Field
Figure 13 shows the contour of temperature field for plain tube and the quadrant-cut twisted tape inserts with twisted ratio () and cut depths (, 1, and cm). The figure shows that the QCT with a cut depth of cm produced better temperature distribution on tube wall than plain tube and other cut depths.
9. Conclusion
In the present study, the heat transfer and friction factor of a constant heat-fluxed tube fitted with quadrant-cut twisted tape inserts in laminar flow conditions were simulated using Fluent version 6.3.26. The following conclusions can be drawn from the results of the present study.(i)The data obtained by simulation are well matched with the literature correlations for plain tube with a discrepancy less than ±8% for Nusselt number and ±10% for friction factor, respectively.(ii)The presence of the tube with quadrant-cut twisted tape inserts yields a higher heat transfer rate (Nu) and friction factor of the plain tube due to the formation of vortices that disturbed the boundary layer, and, thus, enhanced mixing by extension heat transfer rate between the fluid at the core and the heated surface can be achieved.(iii)Over the range of Reynolds number considered, the tube equipped with the quadrant-cut twisted with and cut depth of 0.5 cm offered the best heat transfer performance as compared to the plain tube and other twisted tapes used. (iv)The proposed quadrant-cut twisted tape inserts can be used for heat transfer augmentation.
Nomenclature
| : | Energy component in energy equation |
| : | Force component in momentum equation, N |
| : | Fanning friction factor |
| : | Acceleration due to gravity, m/s2 |
| : | 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/m2 |
| : | 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. |
| : | Density component in governing equations |
| : | Stress component in momentum equation, N/m2. |
Acknowledgments
The authors would like to thank Universiti Kebangsaan Malaysia and Ministry of Higher Education for the financial support and Grant research (FRGS/1/2013/TK07/UKM/01/1).