AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. AUGMENTATION OF HEAT TRANSFER IN CORRUGATED TUBE USING FOUR-START SPIRAL WALL Dr. Ghassan Fadhil SMAISIM Mechanical Engineering Department, Kufa University, Iraq Ghassan.Smaisim@Uokufa.Edu.Iq mail: -E Received on 7August 2017 Accepted on 30 October 2017 Published on 15 January 2018 DOI: 10.30772/qjes.v10i4.493 Abstract: This article dealt with an important heat transfer field, which is passive heat transfer technique represented by corrugated tube. The study conducted numerically by ANSYS Fluent 14. The motivation behind the current study was to clarify the characteristics and merits of such tube geometry in terms of heat transfer and pressure drop. The obtained results reported under constant heat flux, temperature-dependent thermo-physical properties and Reynolds number range of 300-1500. Results show good heat transfer enhancement of 6.15-33.24% in spite of an increase in friction factor of 1.80- 2.93 times the smooth values. The corrugated tube with φ =4.76×10 -2 has the top thermal 1.16-1.25 for all Reynolds number. The most important finding is that the gained heat transfer is much more than the increase in pressure loss until a certain threshold of Reynolds number about 900, later, the pressure loss increase and dominates over the heat gained. The effect of Prandtl number on the heat transfer for three types of fluids produced too. A criterion correlation equation of Nusselt number developed to describe the cases of four starts spirally corrugated tubes by deviation ±3% compared with the simulation results. Keywords: Corrugated Tubes; Four-Starts; Severity Index; Heat Transfer Enhancement; Friction Factor. INTRODUCTION Different types of the heat exchangers have been used in the electrical, chemical, automotive and industrial applications. One of the important types in a co-axial heat exchanger is the spirally corrugated tube, which can improve the overall heat transfer performance of the heat exchangers [1]. Corrugated tubes have been developed from single start to multiple starts. Corrugation augments the heat transfer due to the mix of fluid, which produced during the flow separation and then re-attaches stages [2]. It is imperative to find out a way for sharing resources and saving energy due to horrible energy exhausting. Axiomatically, energy resources are running short with, and the globe on the verge of imminent lack [1]. Heat transfer field holds the promises; it is the most important field that could be helpful in such issue. It can reduce the wasted energy significantly in the case of employing an effective and well functional technique in many applications, such as electronic device cooling, heat exchangers, and nuclear reactor cooling [3-7]. In general, the desired magnification of gained heat could be achieved by primarily two techniques. First one; produces a good outcome, but it needs external power and it is called "active technique"[8]. The second one; does not need any external power and it could be achieved by adding something to the fluid or modify the surface, this method is called "passive technique" [9]. There is another method which called ''compound technique '' which involves the two previous types (active and passive) like using an electrostatic field or vibrated flow which is active types in addition to use the nanofluid or vortex mailto:Ghassan.Smaisim@Uokufa.Edu.Iq mailto:Ghassan.Smaisim@Uokufa.Edu.Iq AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. generators which are passive types [8,9]. Therefore, the active technique is not desirable, while the passive is suitable for many applications. In the last few years, passive heat transfer enhancement became very interested due to economical (saving energy) and environmental (weather mitigation) reasons [10]. Hence, many researchers investigate the passive heat transfer enhancement. One of the gurus in this field is Saha [11,12], who paid great effort and conducted many studies especially laminar flow in rippled tubes. The studies were conducted at different boundary condition with various corrugation styles. The reported results showed a remarkable magnification in heat transfer with an increase in pressure drop for the flow in a tube with specific corrugation style. Kareem et. al. [2] introduces an extensive survey review article for the flow in corrugated tubes. Most of the whole studies regarding corrugation were taken into account. The master key of the corrugated tube is the corrugation style and profile as the mentioned review article stated. Moreover, the harmony of the flow and the enhancement of heat transfer directly proportion to the corrugation style and characteristics. In another study, the using innovative corrugation style can maximize the heat transfer and induce better flow pattern with pressure drop reduction as Kareem et. al. [13] stated. Chen et al. [14] carried out the fluid flow and heat characteristics of the four-start rippled pipe. They studied the influence of the corrugation angles, Reynolds number and pressure drop along the longitudinal direction. They concluded that as the angle ratio increases, the friction loss increases. Zimparov [15] investigated experimentally the fluid flow and heat transfer characteristics of a three-start rippled tube and studied the effect of variable relative pitches for twisted tape inserts in the tube. The results proved that the friction coefficients for whirl flow enlarged at the relative pitch H/Di reduced. Ahn [16] investigated experimentally the effect of four-start spirally corrugated tube on the heat transfer performance. The results illustrated that there is a good enhancement of heat transfer than the use of the smooth tube .Liu et al. [17] presented an experiment and numerical analysis of the heat transfer and hydrodynamic performance in a spirally corrugated cooler and studied the influence the depth of the corrugation with pitch ratio on the heat transfer characteristics, the results showed that corrugated tube augments the heat transfer efficiency. Then Liu et al. [18] investigated numerical model on the heat transfer of rippled tubes with many numbers of start spirally tubes (one start to four-starts). Results showed that the Nusselt number for four starts is higher than the one, two or three starts. Promthaisong et al. [19] carried out numerical simulation on heat transfer characteristics of forced convection in turbulent flow regime with the influence of Reynolds number and pitch ratio on the overall heat performance, they concluded that spirally tubes generate vortex to mix the flow during the tube that helps to enhance the heat transfer. Accordingly, the current article deals with an innovative corrugated tube to enhance heat transfer in the laminar flow region using four spirally corrugated tube by study the effect of the five values of corrugation height e, nominal diameter Dn , severity index φ , fluid properties and Reynolds number on the heat and pressure drop characteristics. In order to mimic the real case of the problem, the thermo-physical parameters of the working fluid were taking as a temperature-dependent under constant heat flux condition. 1. MATHEMATIC MODEL First of all, one smooth and five corrugated tubes were created by utilizing SolidWorks software package. The tubes material was selected to be steel of 2 mm thickness and 2 m length. The bore diameter Db of all tubes remain constant at 10mm, on the other hand, the envelope diameter De is varying from 12 to 14 mm with increment step of 0.5mm. The primary parameters that usually used to describe the rippled tube are the corrugation height e, corrugation pitch p, nominal diameter Dn and the severity index φ as shown in Figure 1 . AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 1: Corrugated Tube Geometry Instead of discussing the effect of both corrugation height e and corrugation pitch p individually, the expression of "severity index φ" proposed which include the effect of both corrugation height and pitch. The severity index simply indicates how corrugation is severe and calculates from the following equation [2]: φ ( )⁄ (1) Where, Dn represents the nominal diameter and it often used in the corrugation tubes in order to denote to the "hydraulic diameter" due to irregular tube geometry. The nominal diameter could be calculated from the following equation [10]: ( ) ⁄ (2) The tubes characteristic and dimensions that employed in the current simulations described in the Table 1 and the corresponding cross sections for the cases appear clearly in Figure 2. TABLE 1: TUBES CHARACTERISTICS No. Db De Dn e P φ 1 - - 10.00 - - Smooth Tube 2 10 12.0 11.00 2.0 28 1.29×10 -2 3 10 12.5 11.25 2.5 28 1.98×10 -2 4 10 13.0 11.50 3.0 28 2.79×10 -2 5 10 13.5 11.75 3.5 28 3.72×10 -2 6 10 14.0 12.00 4.0 28 4.76×10 -2 φ =1.29×10 -2 φ =1.98×10 -2 φ =2.79×10 -2 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 151Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 2: Case Studies of the Corrugated Tubes 1.1. GOVERNING EQUATIONS The discretization of the governing equations should be made prior to simulation. For the flow of fluid in the tube, the following three equations are the main governing equations [8]. -Continuity Equation: (3) -Momentum Equation in z-direction (axial direction): ( ) ( ) (4) -Energy Equation: ( ) * ( ) + (5) The other important heat transfer parameter can be compute by using analytical equation, such as local heat transfer coefficient (h), which can be computed as below:- ( ) ( ) [ ( ) ( )] (6) where, both Tw(z) and Tb(z) represent local wall and bulk temperatures respectively. While q s (z) is the locally applied heat flux, which can be determined by:- ⁄ (7) (8) On the other hand, both of local and average Nusselt number determined from the following equations:- ( ) ( ) ( ) (9) ( ) ( ) (10) Eventually, Reynolds number was determined from the following expression (11) 1.2. NUMERICAL PROCEDURE AND BOUNDARY CONDITIONS ANSYS MESHING ICEM CFD 14 is the tool that employed for this purpose. Mesh process acting as an essential character in accuracy and reliability of the results. Tetrahedral elements of constant stress element owning four nodes were applied. For each node, the elements were represented in 3-Dim space φ =3.72×10 -2 φ =4.76×10 -2 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 155Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. possessing three degrees of freedom in X, Y and Z, directions with take care about the meshing qualities like mesh skewness, aspect ratios and mesh smooth that collected to be 0.75%, 1% and 20% respectively. For consistency in the discretization scheme, the density of the grid plays an important parameter to get the right solution [20]. It is known that any increases in the density about the spatial grids will minimize the errors. Element number has an important parameter to get mesh dependence results as shown in Figure 3. In the current study, the grid density limits were 2763566–3094583 tetrahedral cells. Grid Convergence Index (GCI) was performed to verify the meshing quality [21]. Mesh spaces were tested for each case 0.1, 0.2, 0.3, 0.4 and 0.5 mm respectively, where the GCI's were 0.107, 0.131, 0.169 and 0.205 mm respectively. Therefore, the mesh spacing of 0.1 mm was selected for the whole cases as shown in Figure 4. Moreover, the reduction in mesh spacing gives more stability and accurate outcome values, where the bulk fluid temperature becomes reasonable and has precise value, as the mesh space gets finer as shown in both Figures 3 and 4. A commercial code of ANSYS FLUENT 14.0 [22] used to construct the fluid flow model of smooth and corrugated pipes. Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm produced to solve the Navier Stokes equations by coupling the pressure and velocity. In discrete arrangement, second order upwind chosen for computing momentum and energy while standard format chosen for pressure. Convergence conditions absolute value of energy residual is setting less than 1×10 -9 , beside the value of other quality residuals are less than 1×10 -6 . For all cases, the simulations conducted the stainless steel pipe under constant heat flux of 10 KW/m 2 , with steady state and temperature-dependent thermo-physical properties. The working fluid is Newtonian fluid as pure water enters at temperature 300K with Reynolds number range was 300-1500. Figure 3: Element Number vice Bulk Temperature 311 312 313 314 315 316 317 318 319 320 0.5 1 1.5 2 2.5 3 3.5 T e m p e r a tu r e T b Element Number ( x106 ) AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 4: Grid Refinery Process 2. RESULTS AND DISCUSSIONS 2.1. VALIDATIONS Firstly, validation was made to confirm the results of smooth tube; this verification will give a confidence for the final numerical results of the corrugated tubes. The validation was present numerically for fluid flow in the smooth tube of stainless steel tube has 2 mm thickness, nominal diameter Dn of 10 mm and 2 m length. Nusselt number for the smooth pipe compared with both Churchil and Ozoe equation [9] beside Shah and London equation [23]. Churchil and Ozoe [9] prescribed an equation that covers the entrance and fully develops region confirmed for Pr about 0.7 to 10, and has reliable intimate of small Gz and large Pr, as:    (12) /29.6G1)(Pr/0.02071 04.19/ 1*))6.29(1(364.4 3 1 2 3 3 1 2 Z 2 1 3 2 6 1 2                                        Z Zx G GNu Where: (13) )]/(4/Pr).Re[( 4 * nDZ Dx x G n    Moreover, Shah and London equation [23] was: { ( ) ( ) ( ) (14) AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 5 shows a good agreement with the current outcomes with acceptable deviations of ±4.7% and 3.8% with the data reported by [9] and [23] respectively. This deviation came from both assumption approximations and numerical error. On the other hand, the friction factor is habitually used for the validation of pressure drop, where the validation achieved by comparing the numerical results of smooth tube with the Darcy friction factor as in the following Figure 6 and equation (15) with maximum deviations of ±5.87% occurs at Re=1500. f=64/Re (15) Figure 5: Validation of Nusselt Number for Smooth Tube. Figure 6: Friction Factor Validation of Smooth Tube 2.2. HEAT TRANSFER RESULTS The reasonable matching of both Nusselt number and friction factor validation for both the current numerical results give an indication and corroboration about the fidelity of the current study. For further reliability of the ANSYS simulation, the calculated Nusselt number of Severity Index φ=1.98×10 -2 is compared with the resultant values of previous correlation equations for all Reynolds number values. Srinivasan and Christensen [24] developed a correlation for wide range of spirally tube factors with 0.1≤e/Dvi≤0.4, 0.4≤p/Dvi≤7.3, 30≤ ≤58 where Dvi is the volume-based tube inner diameter as:         5000D ReFor Pr) 90 ()()(Re 014.0 (16) 5000D ReFor Pr) 90 ()()(Re 064.0 v i 4.07 0 5.02 9 3.00 6 7.00 .8 4 2 Dv i v i 4.05 9 9.01 0 8.02 4 2.00 .7 7 3 Dv i   Dvi p Dvi e Dvi p Dvi eDvi Nu Also Barba et al [25] developed a correlation equation for laminar flow of 100< Re <800 for ethylene glycol fluid flow as follows: 0 5 10 15 20 25 0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 N u ss le t N u m b e r N u Lenght of the Tube (m) Num. Smooth Tube Shah and London Eq. Churchil and Ozoe Eq. 0 0.05 0.1 0.15 0.2 0.25 0 300 600 900 1200 1500 F r ic ti o n F a c to r f Reynolds Nummber Nu Darcy Eq. Num. Smooth Tube AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Nu=0.7345 Re 0.5689 Pr -0.0374 (17) Abdel-Kariem and Fletcher [26] suggested a laminar correlation equation for fluid flows through corrugation tubes; valid to the range of Reynolds number at 70≤Re≤1200 at different value of , where is the helix corrugation angle as in the expression below: Nu= 0.777 Re 0.444 ( /45) 0.67 (18) Figure 7, shows the numerical Nusselt number versus Reynolds number for difference values in the severity index as well as the validation of the Nusselt number with theoretical values. As Reynolds number increases, Nusselt number Nu for both corrugated pipe and smooth tubes increase steadily. Obviously, the Nusselt number supposed to increase as both corrugation height and pitch increases, in other words, the heat transfer enhancement gets an increase as the severity index increase. As it was expected, the severity index helps to mix the fluid layers next to the tube wall while keeping the core flow intact. This behavior has two merits. First, enhancing the heat transfer by the disturbance near the wall; second; it keeps the flow laminar in the core zone, which minimizes the pressure loss. The mentioned contours figure strongly indicate that the coefficient of the heat transfer gets enhances through the growth of severity index due to the decrease of the overall wall temperature with increase of severity index. This attributes to the turbulence in the secondary flow region that helps to mix the fluid layers and lead to increase the rate of heat transferred from the wall toward the fluid. The differences of the water properties in this study and the ethylene glycol properties in Barba et al [25] beside their fixed parameters like e=1.5mm and p=11.5mm are the reasons of the dissimilar values of Nusselt number. Moreover the variations in helical corrugation angles were the reasons of the difference between the values of this study compared with Abdel-Kariem and Fletcher [26]. Srinivasan and Christensen [24] obtained the maximum divergence was ±4.91% so the numerical study for the corrugated tubes in different values of severity index can be acceptable. Figure 7: Nusselt Number versus Reynolds Number for Different Severity Index. Figure 8 shows the temperature distribution through the corrugated tube for φ =1.29×10 -2 , 2.79×10 -2 , 4.76×10 -2 at Reynolds Number of 300 to 1500 for the corrugated tubes. For the same Reynolds number, and once the severity index increases, the mixing of the fluid develops according to eddies generated due to the increase of the twisted wall compared with low severity index. Although of the figures declared that, the highest temperature increases with rising of the severity index, this is because of a miniature amount of fluid compared with the main flow encloses near the corners which has a high local temperature compared with the mean fluid flow temperature through the pipe. As the inlet velocity increases, Reynolds number increases, and mass flow rate will enlarge, therefore the overall heat transfer will enhance too thus the surface temperature of the tube decreases. 6 8 10 12 14 16 18 20 22 24 300 600 900 1200 1500 N u ss e lt N u m b e r N u Reynolds Number Re ϕ=1.29x10-2 ϕ=1.98x10-2 ϕ=2.79x10-2 ϕ=3.72x10-2 ϕ=4.76x10-2 Srinivasan and Christensen Eq. [24] Barba et al Eq. [25] Abdel-Kariem and Fletcher Eq. [26] AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 154Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. φ =1.29×10 -2 φ =2.79×10 -2 φ =4.76×10 -2 R e = 3 0 0 R e = 6 0 0 R e = 9 0 0 R e = 1 2 0 0 R e = 1 5 0 0 Figure 8: Temperature Contour with Severity Index of 1.29×10 -2 , 2.79×10 -2 , 4.76×10 -2 at Reynolds Number of 300 to 1500 2.3.FRICTION FACTOR RESULTS As it commonly known, when the inlet velocity increases, means Reynolds number increase and the mass flow rate increase too so the velocity contour increases as it clear shown in the Figure 9, for the corrugate tubes at Re=300,900 and 1500. As the severity index rises, eddies will generated gradually, the local velocity will decrease, and the local pressure will increase which is satisfied the Bernoulli Equation. While Figure 10 shows the pressure contours of the integer corrugation height at Re= 300,900 and 1500. It clearly indicates the actual case of losing energy (pressure loss) as both of corrugation height and Re increase. Pressure loss reaches the maximum value at the highest value of φ =4.76×10 -2 and Re=1500. This figure also shows that the maximum pressure occurring at the corrugation gap region due to the induced turbulence caused by the corrugation and there was an increase in pressure loss of range 20.6-39.1%. AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. φ =1.29×10 -2 φ =2.79×10 -2 φ =4.76×10 -2 R e = 3 0 0 R e = 6 0 0 R e = 9 0 0 R e = 1 2 0 0 R e = 1 5 0 0 Figure 9: Velocity Contours with Severity Index of 1.29×10 -2 , 2.79×10 -2 , 4.76×10 -2 at Re=300 to1500 φ =1.29×10-2 φ =2.79×10-2 φ =4.76×10-2 R e = 3 0 0 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. R e = 6 0 0 R e = 9 0 0 R e = 1 2 0 0 R e = 1 5 0 0 Figure 10: Pressure Contours with Severity Index of 1.29×10 -2 , 2.79×10 -2 , 4.76×10 -2 at Re=300 to1500 As it commonly known, the energy dissipation increases as the flow rate increase. Analogously, the increase in Reynolds number causes an increase in the pressure loss as shown in Figure 11. The desired heat gained is accompanied with appreciable pressure loss and pumping power. The figure shows the of increasing Reynolds number on the friction factor of the numerical study values at the severity index φ=1.98×10 -2 and compared with another previous correlation equations. Srinivasan and Christensen [24] offered a correlation equation as:                     )19( 1500D ReFor 1500D ReFor v i 3 8 4.0 )(0.4 5 Dv i Re6 4 .0 5 5 4.0 4 2.2 ) 9 0 ( )/0 8 3.24 5 4.1( )( * v i )/0 5.02 6.1( )( 0 .2 6 1- Dv i Re 2 0 9.1 )/6 7.36 6 9.2( ) 9 0 ( )/0 3 3.26 6.1( )( * Dvi e Dvie Dvi p Dvip Dvi e DvieDvie Dvi p f   AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. While Barba [25] manifested another correlation equation for ethylene glycol fluid flow tested to 100< Re <800 as follows: f=61.639 Re -0.8602 (20) Meng et al [27] presented correlation equation for oil flow in circular and elliptical tubes by p/D=2–3 for 500≤Re≤50000 as: f= 1.54 Re -0.32 (21) The deviation of Barba [25] and Meng et al [27] equations were as the different of the thermo physics properties between the water for this study and the other types of the work fluids. Srinivasan and Christensen equation obtained a maximum divergence of ±3.1% so that it can be conclude that the other values of the Severity Index are acceptable. Figure 11: Friction Factor against Reynolds Number for Various Severity Index 2.4.THERMAL PERFORMANCE FACTOR Figure12 shows Nusselt number ratio Nuc/Nus between the corrugated tubes with the smooth tube against Reynolds number. The Nu ratio indicates that the corrugated pipe is better than the smooth pipe that means that the heat transfer performance of the corrugated pipe is more efficient and the largest value in Nuc/Nus was reach 1.560. In additional the resistance coefficient ratio; fc/fs through the corrugated tube increases directly with booth the Reynolds number and Severity Index along with maximum value 1.80-2.93 times the smooth tube for Reynolds number from 300-1500 respectively as shown in Figure 13. Figure 12: Nusselt Number Ratios with Reynolds Number 0.03 0.06 0.09 0.12 0.15 0.18 0.21 0.24 0.27 0.3 0.33 0.36 0.39 0.42 0.45 300 600 900 1200 1500 F r ic ti o n F a c to r f Reynolds Number Re ϕ=1.29x10-2 ϕ=1.98x10-2 ϕ=2.79x10-2 ϕ=3.72x10-2 ϕ=4.76x10-2 Srinivasan and Christensen Eq.[24] Barba Eq.[25] Meng et al Eq. [27] 1 1.1 1.2 1.3 1.4 1.5 1.6 300 600 900 1200 1500 N u c/ N u s Reynolds Number Re ϕ=1.29x10-2 ϕ=1.98x10-2 ϕ=2.79x10-2 ϕ=3.72x10-2 ϕ=4.76x10-2 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 13: Friction Factor Ratios with Reynolds Number This factor was proposed to describe the overall performance of an enhanced tube or tube with the promoter. The gained heat transfer ratio to the gained pressure loss simply represents the thermal performance factor as in the following equation [20]. η= (Nuc / Nus) / (fc / fs) 1/3 = f(Res) (22) The thermal performance of the current cases shows in Figure 14. The figure shows the tube of highest severity index is the best one in term of overall thermal performance, while the tube of lowest severity index is the worst one and the maximum enhancement was in the range of Reynolds numbers 900. Figure 14: Thermal Performance Factor against Reynolds Number for Various Severity Index 2.5 ANALYSIS OF DIFFERENT FLUID TYPES Air, water and engine oil investigated to explore the effects of the medium properties on the overall characteristics of the heat transfer for smooth and corrugated tube. By varying the medium type, the fluid properties change and the relation between fluid property like Pr and Nu could be produced by the numerical method. Physical characteristics of the fluids at temperature 350 K exposed in Table 2 [28]. Table 2: PHYSICAL PROPERTY OF THE WORKING FLUIDS [28] Pr K W/m. 0 C μ N.s/m 2 Cp KJ/Kg.K ρ kg/m 3 Fluid at(350K) 0.7 0.03003 0.00002082 1.005 0.995 Air 2.33 0.668 0.000372 4.191 973.7 Water 546 0.138 0.0356 2.118 583.9 Engine Oil 0.5 1 1.5 2 2.5 3 300 600 900 1200 1500 fc /f s Reynolds Number Re ϕ=1.29x10-2 ϕ=1.98x10-2 ϕ=2.79x10-2 ϕ=3.72x10-2 ϕ=4.76x10-2 1 1.05 1.1 1.15 1.2 1.25 1.3 300 600 900 1200 1500 T h e rm a l P e rf o rm a n ce F a ct o r η Reynolds Number Re ϕ=1.29x10-2 ϕ=1.98x10-2 ϕ=2.79x10-2 ϕ=3.72x10-2 ϕ=4.76x10-2 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 141Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. As it concluded in the water cases study, the best corrugation type to enhance the heat transfer depends on the Severity Indexes occurs at φ=4.76×10 -2 , so that this severity index will used to get the evaluations for the other fluid mediums. The distribution of augmentation ratio of the heat transfer Nuc/Nus that evaluated for the three fluids through the smooth and corrugated tube as shown in Figure 15. It can be conclude that the Nusselt number Nu of four start corrugate pipe is better than the smooth pipe and for the greater Prandtl number working fluid like engine oil, Nusselt number will be better than the lower Prandtl number working fluid like in air. so that It can be establishing that the progression of Nusselt number ratio Nuc/Nus from low to high as: air, water, engine oil. Largest ratio of Nusselt number Nuc/Nus could be reach 5.778, which points to that for the engine oil medium with a high Pr cause a higher value in the augmentation of the heat transfer. Figure 15: Nusselt Number Ratios with Reynolds Number for All Fluid Media Figure 16 shows the friction factor ratio for the different working mediums (air, water and engine oil) at Severity Index φ=4.76×10 -2 with all range of Reynolds numbers. It can be shown the curves act as the same behavior of the water and the effect of the fluid properties on the resistance of the fluid inside the pipe is not large so that the curves of the friction factor ratio fc/fs is similar. It's clear that the friction factor for all fluid mediums for the corrugated tube is higher than the smooth tube. The curve of the engine oil medium is higher the others as the high viscosity compared with the other working fluids so the pumping power will be more than the other fluids. Figure 16: Friction Factor Ratio with Reynolds Number for All Working Fluid Figure 17 shows the thermal performance factor against the Reynolds number for the corrugated tube with Severity Index φ=4.76×10 -2 for the air, water and engine oil. It can be noted that the engine oil has a higher thermal performance than the other fluids as the high values of the Nusselt number ratio Nuc/Nus that acts in significant outcome more than the effect of its friction factor ratio fc/fs. The maximum values of the Thermal performance ratio was not in the same Reynolds number but its change in the range of Re=800- 1200 from the low Prandtl number to the high Prandtl number as the fluid properties are changed. 1 2 3 4 5 6 300 600 900 1200 1500 N u c/ N u s Reynolds Number Re Air Water Engine Oil 1 1.5 2 2.5 3 3.5 300 600 900 1200 1500 fc /f s Reynolds Number Re Air Water Engine Oil AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 145Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Figure 17: Thermal Performance Factor against the Reynolds Number for All Fluids As shown before, there are many parameters effects on the heat transfer phenomena as the geometry factors (corrugation depth e and pitch p), the properties of the fluid flow and Reynolds number Re, so that the non dimensional factors like e/dn, p/dn, Re and Pr are confirmed by Zhi-Jiang Jin [1] as parameters of the correlation criteria. The general form can expressed as: Nu=F(e/dn,p/dn, Re,Pr) (23) Zhi-Jiang Jin [1] mentions that the correlation equation could be adjust in power function as: Nu=e m1 (e/dn) m2 (p/dn) m3 Re m4 Pr 1/3 (24) The statistic program SPSS was used to interpret the output results and get the correlation equation using Multiple Regression Analysis. The criteria of the analysis of multiple regression results were; the Multiple R is 0.999949, R Square is 0.999899, Adjusted R Square is 0.952265 and the Standard Error is 0.021933. The criterion correlation of the heat transfer through five cases of four-start spirally corrugated pipes is: Nu= e 0.38887 (e/dn) 0.03644 (p/dn) -0.13939 Re 0.26243 Pr 1/3 (25) This criterion correlation is satisfied for water flow inside corrugated tube with 0.1818≤e/dn≤0.333, 2.333≤p/dn≤2.545, 300≤Re≤1500, 500≤L/e≤1000 and 166.667≤L/dn≤181.82. Then a comparison of Nu between simulation and calculation results is shows in Figure 18. The deviation between these results is inside ±3%, so that the correlation for the calculations of heat transfer through the corrugated tube is reliable. Figure 18: Comparison Nusselt Number Results CONCLUSION Four-start corrugated pipes simulated under constant heat flux condition to clarify the heat transfer enhancement in such passive technique. The outcomes strongly indicated that obtained heat transfer is in 0 1 2 3 4 5 300 450 600 750 900 1050 1200 1350 1500 T h e rm a l P e rf o rm a n ce F a ct o r η Renolds Number Re Air Water Engine Oil 5 7.5 10 12.5 15 5 7.5 10 12.5 15 C a lc u la ti o n N u Simulation Nu ±3% Deviation AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. directly proportion with the severity index, also with increase the Reynolds number and the maximum Nusselt number ratio was 1.56 at the highest Reynolds number through φ=4.76×10 -2 . Another conclusion is the temperature gradient at the center of pipe is miniature whereas a notice change of the temperature was near the wall that means the main heat transfer procedures near the pipe wall. At constant Reynolds number Re, as the ratio e/dn increase, the Nusselt number Nu increases. It was found that although there are an increase in friction factor of range 1.80-2.93 times the smooth tube values, the enhancement of heat transfer range was about 6.15-33.24%. As Reynolds number Re increases, the ratios Nuc/Nus and fc/fs increases. The most important finding is that the increases in the heat transfer is much more than the increase in pressure loss until a certain threshold (Re=900) and later, the pressure loss increase and dominate over the heat gained. It can be concluded a similar behavior for all working fluids, where Nusselt number increase as Reynolds number increases for both the smooth tube and corrugated tube and its values in the corrugated tube are always more the Nusselt number for the smooth tube. Moreover the heat transfer performance of the fluids that has high Prandtl number is best than the lower Prandtl number and the thermal performances was 1.18 times the smooth tube values for air at Re=800 and 4.05 times the smooth tube values for engine oil at Re=1200. Finally, the heat transfer calculation presented as a criterion correlation for corrugated tubes and it validated and proposed to be reliable with ±3% with the simulation results. REFERENCES 1. Jin, Zhi-Jiang; Chen, Fu-Qiang; Gao, Zhi-Xin; Gao ,Xiao-Fei, Effects of Pitch and Corrugation Depth on Heat Transfer Characteristics in Six-Start Spirally Corrugated Tube, Jin-yuan Qian, International Journal of Heat and Mass Transfer, 108, 2017, pg. 1011–1025. 2. Kareem ,Zaid S.; Jaafar, M.N. Mohd; Lazim ,Tholudin M.; Shahrir Abdullah; Abdulwahid ,Ammar F., Passive Heat Transfer Enhancement Review in Corrugation, Experimental Thermal and Fluid Science, 68, 2015, pg. 22–38. 3. Smaisim, G. F.; Fatla, O. M.H.; Medina, A. Valera; Rageb, A.M.; Syred, N., Investigation of Heat Transfer and Fluid Mechanics across a Heated Rotating Circular Cylinder in Crossflow, AIAA Science and Technology Forum and Exposition, San Diego, U.S.A. 4-8 Jan, 2016. 4. Balla ,Hyder H., Enhancement of Heat Transfer in Six-Start Spirally Corrugated Tubes, Case Studies in Thermal Engineering, 9, 2017, pg. 79–89. 5. Liang, Caihang; Tong , Xiaoman; Lei , Tengyue; Li ,Zhenxing; Wu ,Guoshan, Optimal Design of an Air-to-Air Heat Exchangerwith Cross-Corrugated Triangular Ducts by Usinga Particle Swarm Optimization Algorithm, Applied Science 7,2017,pg.554-562. 6. Ardila Marin, J.; Hincapie Zuluaga, D.,Casas Monroy, J., Comparison and Validation of Turbulence Models in The Numericalstudy of Heat Exchangers, TECCIENCIA, Vol. 10 No. 19., 2015 ,pg.49-60. 7. Lazim ,Tholudin Mat; Kareem, Zaid S.; Jaafar,M. N. Mohd; Shahrir Abdullah; Abdulwahid, Ammar F., Heat Transfer Enhancement in Spirally Corrugated Tube, International Review on Modelling and Simulations (I.RE.MO.S.), 7, (6), 2014, pg. 123-132. 8. Webb,Ralph L. and Kim, Nae-Hyun,'Principles of Enhanced Heat Transfer', 2Nd Edition, Taylor & Francis Routledge, 2006. 9. Churchill, S. W., and Ozoe, H., 'Correlations for Forced Convection with Uniform Heating in Flow Over a Plate and in Developing and Fully Developed Flow in a Tube', J. Heat Transfer, 95, 78–84, 1973. 10. Kareem, Zaid S.; Shahrir Abdullah; Lazim, Tholudin M.; Jaafar, M. N. Mohd; AbdulWahid, Ammar F., Heat Transfer Enhancement in Three-Start Spirally Corrugated Tube: Experimental and Numerical Study, Chemical Engineering Science,134,2015, pg.746–757. 11. Saha, S. K., Thermohydraulics of Laminar Flow of Viscous Oil Through a Circular Tube Having Axial Corrugations and Fitted With Centre-Cleared Twisted-Tape, Exp.Thermal Fluid Sci.,38,2012,pg.201–209. 12. Saha, S. K.,Thermohydraulics of Laminar Flow Through a Circular Tube Having Integral Helical Corrugations and Fitted with Helical Screw-Tape Insert, Chem. Eng. Commun. 200, 2013, pg.418–436. 13. Kareem, Zaid S.; Jaafar, M.N. Mohd; Lazim, Tholudin M.; Shahrir Abdullah; AbdulWahid, Ammar F., Heat Transfer Enhancement in Two-Start Spirally Corrugated Tube, Alexandria Engineering Journal 54,2015, pg. 415-422. AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 10, No. 4 ISSN: 1998-4456 144Page Copyright @ 2017 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. 14. Chen X.D.; Xu,X.Y.; Nguang,S.K., Characterization of The Effect of Corrugation Angles on Hydrodynamic and Heat Transfer Performance of Four-Start Spiral Tubes, J. Heat Transfer,123,(6),2001, pg.1149–1158. 15. Zimparov, V., Enhancement of Heat Transfer by a Combination of Three-Start Spirally Corrugated Tubes With a Twisted Tape, Int. J. Heat Mass Transf.,44, (3),2001, pg. 551–574. 16. Ahn,S.W.,Experimental Studies on Heat Transfer in The Annuli With Corrugated Inner Tubes, KSME Int. J. 17, (8),2003, pg.1226–1233. 17. Liu ,L.; Ling, X.; Peng,H. Analysis on Flow and Heat Transfer Characteristics of EGR Helical Baffled Cooler With Spiral Corrugated Tubes, Exp. Thermal Fluid Sci.,44, 2013, pg. 275–284. 18. Liu, J.J.; Liu, Z.C.; Liu, W., 3D Numerical Study on Shell Side Heat Transfer and Flow Characteristics of Rod-Baffle Heat Exchangers With Spirally Corrugated Tubes, Int. J. Therm. Sci., 89, 2015, pg.34–42. 19. Promthaisong ,P.; Jedsadaratanachai ,W.; Eiamsa-ard, S., 3D Numerical Study on The Flow Topology and Heat Transfer Characteristics of Turbulent Forced Convection in Spirally Corrugated Tube, Numer. Heat Trans. Part A: Appl. Int. J. Comput. Methodol. 69, (6), 2016, pg. 607–629. 20. Han, H., Li ,B., Wu, H.; Shao ,W., Multi–Objective Shape Optimization of Double Pipe Heat Exchanger With Inner Corrugated Tube Using RSM Method, Int. J. Therm.Sci. 90, 2015, pg. 173–186. 21. Versteeg,H.K. ; Malalasekera, W., An Introduction to Computational Fluid Dynamics the Finite Volume Method, Pearson Education Lmited, 2nd Edition, 2007. 22. ANSYS FLUENT User's Guide. Release 14.0. ANSYS, Inc. November 2011. 23. Shah, R.K.; London, A.L., Laminar Flow Forced Convection in Ducts, 1st Ed., Academic Press, Massachusetts, New York, 1978. 24. Srinivasan, V. and Christensen, R.N., Experimental Investigation of Heat Transfer and Pressure Drop Characteristics of Flow Through Spirally Fluted Tubes, Exp. Thermal Fluid Sci., 5, 1992, pg. 820–827. 25. Barba,A.; Rainieri, S.; Spiga,M., Heat Transfer Enhancement in Corrugated Tube, Int. Commun. Heat Mass Trans., 29, 2002, pg. 313–322. 26. Abdel-Kariem A.H.; Fletcher L.S., A Comparative Analysis of Heat Transfer and Pressure Drop in Plate Heat Exchangers, Proceedings of the Fifth ASME/JSME Thermal Engineering Conference, USA, 1999. 27. Meng,J.; Liang X.; Chen ,Z.; Li, Z., Experimental Study on Convective Heat Transfer in Alternating Elliptical Axis Tubes, Exp. Therm. Fluid Sci.,29, 2005, pg. 457–465. 28. Theodore L. Bergman;Adrienne S. Lavine; Frank P. Incropera; David P. Dewitt, Fundamentals of Heat and MassTransfer. 7th Ed, JOHN WILEY & SONS,2011. Symbols Subscripts e Corrugation Height r Radial Direction b Bore f Friction Factor Re Reynolds Number B Bulk GCI Grid Independence Index T Temperature c Corrugated h Coefficient of Heat Transfer u Fluid Velocity e Envelope k Thermal Conductivity z Axial Direction n Nominal L Length of Tube Greek Symbols s Smooth Nu Nusselt Number φ Severity Index z Local P Pressure Density of Fluid p Corrugation Pitch ν Kinematic Viscosity of Fluid Pr Prandtl Number Dynamic Viscosity q Heat Flux Angular Direction