Numerical evaluation of laminar heat transfer enhancement in nanofluid flow in coiled square tubes

Convective heat transfer can be enhanced by changing flow geometry and/or by enhancing thermal conductivity of the fluid. This study proposes simultaneous passive heat transfer enhancement by combining the geometry effect utilizing nanofluids inflow in coils. The two nanofluid suspensions examined in this study are: water-Al2O3 and water-CuO. The flow behavior and heat transfer performance of these nanofluid suspensions in various configurations of coiled square tubes, e.g., conical spiral, in-plane spiral, and helical spiral, are investigated and compared with those for water flowing in a straight tube. Laminar flow of a Newtonian nanofluid in coils made of square cross section tubes is simulated using computational fluid dynamics (CFD)approach, where the nanofluid properties are treated as functions of particle volumetric concentration and temperature. The results indicate that addition of small amounts of nanoparticles up to 1% improves significantly the heat transfer performance; however, further addition tends to deteriorate heat transfer performance.


Introduction
Convective heat transfer can be enhanced by active as well as passive methods. While the former usually provide better enhancement, it requires additional external forces and/or equipment which can increase the complexity, capital, and operating costs of the system. In contrast, passive heat transfer enhancement can be achieved by changing flow geometry or modifying thermo-physical properties of working fluid. Hence, it is generally a more desirable approach when compared to an active method. In our previous study [1][2][3] (Sasmito AP, Kurnia JC, Mujumdar AS: Numerical evaluation of transport phenomena in a T-junction micro-reactor with coils of square cross section tubes, submitted), we have shown that coiled tubes provide better heat transfer performance relative to straight tubes under certain conditions. In this study, the potential application of coiled tubes using nanofluids to improve heat transfer performance is investigated.
Coiled tubes have been known as one of the passive heat transfer enhancement techniques in heat and mass transfer applications due to the presence of secondary flows which improve heat and mass transfer rates. They have been widely used in process industries, e.g., heat exchangers and chemical reactors, due to their compact design, high heat transfer rate, and ease of manufacture. Aside from their industrial applications, studies of the transport phenomena in coiled duct have also attracted many attention from engineering researchers. The presence of secondary flows induced by coil curvature and the complex temperature profiles caused by curvatureinduced torsion are among significant phenomena which can be observed in coiled tubes. Numerous experimental [4][5][6][7][8] and numerical [1][2][3][9][10][11][12][13] investigations on heat transfer and flow characteristics inside coiled tubes have already been reported. Furthermore, reviews on the flow and heat transfer characteristics and potential application of coiled tubes in process industries and heat transfer application can be found in [14,15].
It is well known that conventional heat transfer fluids including water, oil, and ethylene glycol mixtures have poor heat transfer rate due to their low thermal conductivity. Therefore, over the past decade, extensive research have been conducted to improve thermal conductivity of these fluids by suspending nanoparticles of diverse materials in heat transfer fluids, called nanofluids [16]. Modern technology provides opportunities to process and produce particles below 50 nm. It is also expected that nanofluids should provide not only higher heat transfer rate, but also good stability of the suspension by eliminating possible agglomeration and sedimentation to permit long-term application [17]. To date, several experimental (see for example [18][19][20][21][22][23]) and numerical (see for example [24][25][26][27][28]) investigations to characterize heat transfer performance of nanofluids have been already reported. Choi et al. [18] showed that addition of small amounts of less than 1% nanoparticles can double the thermal conductivity of working fluids. Vajjha et al. [24] showed that heat transfer rate increases up to 94% by adding 10% Al 2 O 3 nanofluid and increase up to around 89% by adding 6% CuO nanofluid. In addition, the comprehensive reference on nanofluids can be found in the book of Das et al. [29], while several reviews of nanofluids are available in the literature [30][31][32][33][34][35][36][37][38][39][40][41][42].
It has been shown that coiled tubes geometry and nanofluids can passively enhanced heat transfer performance. Now, to maximize the advantages of the heat transfer enhancement, we propose to combine both techniques simultaneously; i.e., employing the combination of coiled tubes filled with nanofluids. Therefore, the aim of the study presented here is threefold: (i) to investigate the heat transfer performance of various configurations of coils of square tubes, e.g., conical spiral, inplane spiral, and helical spiral, relative to the straight pipe; (ii) to evaluate simultaneous passive heat transfer enhancement-channel geometry and fluid thermo-physical properties-in coiled tubes filled with nanofluids; (iii) to study the heat performance of two different nanofluids, water-Al 2 O 3 and water-CuO, in coiled tubes at various nanoparticle concentrations. The most significant aspect of this study is to determine the potential advantages and limitations of heat transfer enhancement of coiled of square tubes filled with nanofluids and provide design guidelines for their applications through mathematical modeling.
The layout of the article is as follows. First, the mathematical model is introduced; it comprises conservation equations for mass, momentum, and energy. The nanofluid thermo-physical properties are treated as functions of particle volumetric concentration and temperature. The mathematical model is then solved numerically utilizing finite-volume-based CFD software Fluent 6.3.26, the User-Defined Function written in C language is used extensively to capture the nanofluid properties. The model is further validated against experimental data by Anoop et al. [19] in terms of heat transfer performance for both base-fluid and nanofluid. Fluid flow and heat transfer performance of various coiled tube designs filled with nanofluids is evaluated in terms of a figure of Merit Defined later. Parametric studies for particle concentration and nanofluid type are then carried out. Finally, conclusions are drawn and possible extensions of the study are highlighted.

Mathematical model
The physical model (see Figure 1) comprises four tube designs, e.g., straight pipe, conical spiral, in-plane spiral, and helical spiral, filled with two different nanofluids (water-Al 2 O 3 and water-CuO). We assume that the low particle volumetric concentration of nanoparticles (less than 3%) in the base-fluid makes it behave like a singlephase fluid and there is no agglomeration or sedimentation which occurs inside the tubes. A constant wall temperature is prescribed along all sides of the channel wall; the nanofluid is assumed incompressible and Newtonian. Furthermore, to ensure fidelity of the comparison of heat transfer performance for each tube design, the total length of each tube design is kept constant. Since this study relates only to laminar flow, a precise numerical solution is adequate to simulate reality very closely.
In the above equations, r nf is the nanofluid fluid density, u is the fluid velocity, p is the pressure, μ nf is the dynamic viscosity of the nanofluid, c p,nf is the specific heat of the nanofluid and k nf is thermal conductivity of the nanofluid.

Thermo-physical properties of nanofluids
The thermo-physical properties of nanofluid are functions of particle volumetric concentration and temperature. The nanofluid density is given by [24,29] where r np and r w is the nanoparticle density and water density, respectively, while j is the particle volumetric concentration. The nanofluid viscosity is estimated by [24] μ nf = C 1 exp (C 2 φ)μ w , where C 1 and C 2 are constants (summarized in Table  1), and μ w is the viscosity of base-fluid.
The specific heat of nanofluid is assumed to be a weighted average of the base-fluid and the nanoparticles, e.g., where c p,np and c p,w are the specific heats of nanoparticle and water, respectively. In this model, the thermal conductivity considers a combination of the static part of Maxwell's theory and the dynamic part taking the contribution of the Brownian motion of nanoparticles, defined as [24] where d np is the nanoparticle diameter, k 1 is the Brownian motion constant, k np and k w are thermal conductivity of nanoparticle and water, respectively. Here, the effect of temperature and particle volumetric concentration is taken into account in the Brownian motion from empirical data given by [24] where b 1 , b 2 , c 1, c 2, c 3 and c 4, are constants (see Table 1).

Thermo-physical properties of base-fluids
The base-fluid considered in this article is water. Thermo-physical properties of water were obtained as polynomial functions of temperature [43]; the water density is defined by while the water viscosity is given by and the thermal conductivity of water is calculated from The specific heat of water is considered constant at Properties of nanoparticles are given in Table 1.

Heat transfer performance
The heat transfer performance of the cooling channel is discussed in terms of the figure of merit, FoM, which is defined as where W pump is the pumping power required to drive the fluid flow through the channel. It is given by Here, h pump is the pump efficiency (assumed to be 70%), W is the total heat transfer rate, and Δp is the pressure drop in the cooling channel. The total heat transfer rate is given as whereṁ is the mass flow rate and T m,in and T m,out are mixed mean temperature at the inlet and outlet, respectively. The mixed mean temperatures is calculated as where A c is the cross section area of the channel and V is the mean velocity given by

Boundary conditions
The boundary conditions for the flow inside the channel are prescribed as follows • Inlet At the inlet, we prescribe inlet mass flow rate and inlet temperature.
• Outlet At the outlet, we specify the pressure and streamwise gradient of the temperature is set to zero; the outlet velocity is not known a priori but needs to be iterated from the neighboring computational cells.
• Walls At walls, we set no slip condition for velocities and constant wall temperature.
In this article, a constant mass flow rate at a Reynolds number (Re = rUD h /μ) of approximately 1000 is prescribed at the inlet for comparison purposes.

Numerics
The computational domains (see Figure 2) were created in AutoCAD 2010; the commercial pre-processor software GAMBIT 2.3.16 was used for meshing, labeling boundary conditions and determines the computational domain. Three different meshes, 1 × 10 5 , 2 × 10 5 , and 4 ×10 5 , were tested and compared in terms of the local pressure, velocities, and temperature to ensure a mesh independent solution. It is found that mesh number of around 2 × 10 5 gives about 1% deviation compared to mesh size of 4 × 10 5 ; whereas the results from mesh number of 1 × 10 5 deviate by up to 8% compared to those from the finest one. Therefore, a mesh of around 2 × 10 5 (20 × 20 × 500) elements was considered sufficient for the numerical investigation purposes; a fine structured mesh near the wall to resolve the boundary layer and an increasingly coarser mesh in the middle of the channel to reduce the computational cost.
Equations 1-3 together with appropriate boundary conditions and constitutive relations comprising of five dependent variables, u, v, w, p, and T, were solved using the finite volume solver Fluent 6.3.26. User-Defined functions (UDF) were written in C language to account for particle volumetric concentration and temperature-dependence of the thermo-physical properties of the nanofluids. The equations were solved with the well-known Semi-Implicit Pressure-Linked Equation (SIMPLE) algorithm, first-order upwind discretization and Algebraic Multi-grid (AMG) method. As an indication of the computational cost, it is noted that on average, around 200-500 iterations and 500 MB of Random Access Memory (RAM) are needed for convergence criteria for all relative residuals of 10 -6 , this takes 5-30 min on a workstation with a quadcore processor (1.83 GHz) and 8 GB of RAM.

Results and discussion
The numerical simulations were carried out for four different tube geometries, four different nanofluid concentrations, and two different nanofluid suspensions. The base-case conditions together with the physical parameters are listed in Table 1, while the geometric details can be found in Table 2.

Validation
When developing and implementing mathematical model to predict the behavior of nanofluid heat transfer, one needs to pay special attention to validation of the model due to inherent complexity of coupled physical phenomena and interaction between base-fluid and nanoparticle. In this study, we aim to validate our model with an experimental nanofluid heat transfer by Anoop et al. [19], which has error of approximately 4%. The heat transfer performance of nanofluid flows in circular tube with diameter 4.75 × 10 -3 m and length of 1.2 m is approximated with 2D axisymmetric model, see Anoop et al. [19] for details of the experimental setup.
The validation is initiated with heat transfer performance of water flowing at a constant Reynolds approximately 1580; after which, the heat transfer performance of 4 wt% of water-Al 2 O 3 nanofluid with nanoparticle size 45 nm flows at Reynolds approximately 1588 is compared, as depicted in Figure 3. It is found that the model predictions agree well with the heat transfer performance from  This implies that the model correctly accounts for the fundamental physics associated with nanofluid heat transfer.

Effect of geometry Base-fluid
One of the key factors that determine the heat transfer performance is the cross-sectional tube geometry. This study examines four different square cross section tubes geometries: straight, conical spiral, in-plane spiral, and helical spiral with water as the base working fluid. Since the convective heat transfer inside the tube is directly linked to flow behavior, it is of interest to investigate the flow patterns inside the tubes. In our previous studies [1][2][3], albeit using air as working fluid, showed that the presence of centrifugal force due to curvature leads to significant radial pressure gradients in flow core region. In the proximity the inner and outer walls of the coils, however, the axial velocity and the centrifugal force will approach zero. Hence, to balance the momentum transport, secondary flow should develop along the outer wall. This is indeed the case, as can be seen in Figure 4, where the secondary flow with higher velocities is generated in the outer wall region of coiled tubes (see Figure 4b,c,d). However, this is not the case for the straight tube ( Figure 4a) as a fully developed flow exists inside the tube. It is noted that at this particular Reynolds number (approximately 1000), the secondary flows appear as one-pair for conical spiral and helical spiral tubes; whereas in the in-plane spiral tube, the secondary flows appeared as two-pairs. The presence of secondary flow with high velocities is expected to have direct impact on the heat transfer rate. This can be inferred from Figure 5 which presents temperature distribution over the cross sections of various tube designs. As can be seen from Figure 5, temperatures in coiled tubes are higher than in straight tube at the same axial distance which indicates that coiled tubes have higher heat transfer rate when compared to that of the straight tube due to the presence of secondary flows. It is also worth noting that the higher intensity of secondary flow will tend to lead to higher heat transfer rate. Now looking at the mixed mean temperature and total heat transfer variation along the tube length (see dotted line in Figure 6), it is noted that coiled tubes have superior heat transfer performance when compared to that of the straight tube; the total heat transfer rate can be up to  Figure 3 Comparison of heat transfer coefficient between simulation (lines) and experimental data [19] (symbols) for water and nanofluid.
almost three times higher than that for the straight tube. In the near-inlet region, the heat transfer performance of in-plane spiral yields the best result among others, followed by conical spiral and helical spiral; whereas, in the near-outlet region, the helical coil performs the best followed by in-plane spiral and conical spiral. This indicates that, for water as working fluid, in-plane spiral is more effective to be used in short tube applications, while the helical spiral is more effective for long tube applications in terms of amount of heat transferred.

Nanofluids
Four square cross section tube geometries were examined for flow of nanofluid suspensions of water-Al 2 O 3 with nanoparticle concentration of 1%. The results are depicted in Figure 6 where the mixed mean temperature and total heat transfer of base-fluid and nanofluids are shown. It is noted that adding 1% concentration of Al 2 O 3 in water improves the heat transfer performance. The total heat transfer for straight tube increases up to 50% as compared to that for water, whereas for coiled tubes, the heat transfer improves by about 50% in the near-inlet region and then decreases toward the outlet. Furthermore, among the coiled tube geometries, in-plane spiral gives the highest heat transfer improvement, followed by helical spiral and conical spiral tubes. This implies that in-plane spiral tube may have potential application to be used along with nanofluid due to its higher heat transfer performance. Therefore, the most of the following results refer to in-plane spiral coils.

Effect of nanoparticle concentration
The amount of nanoparticles suspended in the basefluid plays a significant role in deter-mining heat transfer performance. Intuitively, adding larger amount of nanoparticles in the base-fluid increases thermal conductivity of the nanofluid; however, care has to be taken as it also increases the friction factor and may reduce the stability of nanofluids due to agglomeration and sedimentation. To study the impact of these factors, we investigated four different nanoparticle concentrations: 0, 1, 2, and 3% of Al 2 O 3 in the base-fluid (water). Figure 7 displays the velocity profiles for the in-plane spiral tube for various nanoparticle concentrations. Interestingly, the velocity profiles are not strongly affected by the additional nanoparticle suspension, especially at low concentrations. We note that at 1 and 2% of Al 2 O 3 concentration, there is no significant difference on the secondary flow development inside the tube; whereas, at 3% Al 2 O 3 concentration, the effect of nanofluid suspension becomes stronger: the secondary flow appears in two-pairs as compared to that in one-pair at lower nanoparticle concentrations. A plausible explanation is the fact that nanofluid suspension does not significantly change viscosity of the fluid. Conversely, this is not the case for thermal conductivity of the nanofluid, as mirrored in Figure 8, where the addition of small amount of nanoparticle (1%) drastically changes the temperature profiles inside the tube. Furthermore, the temperature profiles for higher amount of nanoparticle concentration (2 and 3%) also slightly change, but they are mainly affected by the hydrodynamics (secondary flows).
Proceeding to the local mixed mean temperature and total heat transfer along the tube, as illustrated in Figure 9, it is clearly seen that additional small amounts     of nanoparticles improves the heat transfer performance significantly, especially in the near-inlet area. How-ever, increase in nanoparticle concentration leads to a reduction of total heat transfer along the tube by approximately 5%. It is noteworthy that adding large amounts of nanoparticles in the suspension is not effective in enhancing heat transfer. Moreover, low nanoparticle concentration also has advantages of better stability of the suspension as it minimizes agglomeration and sedimentation.

Effect of nanofluid type
So far, the simulated nanofluid type chosen was water-Al 2 O 3 ; it is, therefore, of interest to see the heat transfer performance for a different nanofluid. In this study, we compare the performance of water-Al 2 O 3 and water-CuO nanofluids. Note that other types of nanofluid suspensions can be easily simulated within the framework of this model once their properties are known. Figure 10 shows temperature profiles for an in-plane spiral tube flowing through with water (Figure 10a), 1% of Al 2 O 3 nanofluid (Figure 10b) and 1% of CuO nanofluid (Figure 10c). We note that the temperature profiles for both nanofluids (Figure 10b,c) are much higher than that of water (Figure 10a). Closer inspection reveals that a slightly larger area of higher temperature exists for the Al 2 O 3 suspension (Figure 10b) as compared to that for CuO suspension (Figure 10c). This is attributed to the stronger secondary flow observed in Al 2 O 3 nanofluid when compared to that of the CuO nanofluid (not shown here due to page limitation).
The heat transfer performance of two different nanofluid types is further evaluated in terms of the local mixed mean temperature and total heat transfer. As seen in Figure 11, the mixed mean temperature for the nanofluid is around 15% higher than that of water. There is no discernible difference between Al 2 O 3 and CuO suspensions in terms of the mixed mean temperature. For total heat transfer, Al 2 O 3 gives somewhat higher heat transfer (approximately 5%) when compared to the CuO nanofluid. Therefore, it can be deduced that Al 2 O 3 nanofluid performs better heat transfer performance than that of CuO nanofluid, but not significantly. The stability and cost would decide the selection between these two nanofluids.

Overall heat transfer performance
A summary of heat transfer performance for all cases considered in this article is presented in Figure 12. Here several features are apparent; foremost among them is that the coiled tubes provide significantly higher heat transfer than that of straight tube, and addition of a small amount of nanoparticles in the base-fluid enhances heat transfer further (see Figure 12a). It is   noted that the maximum heat transfer performance is achieved at 1% nanoparticle concentration, decreasing with higher amounts of nanoparticles. Aside from higher heat transfer performance, keeping pressure drop at a minimum is of interest for reducing the operating cost and saving energy. Figure 12b shows a summary of the pressure drop required for all cases studied. Note that the mass flow rate is kept constant in all cases; hence, it can be used directly to represent the pumping power required. The straight channel requires the lowest pressure drop among all cases; whereas the coiled tube designs require more than double the pressure drop of the straight channel. Among the coiled tubes, helical spiral tube needs the highest pressure drop, followed by in-plane spiral and conical spiral tubes. An interesting phenomenon is observed at a nanofluid concentration of 1% when the pressure drop for coiled tubes is slightly lower than that for water. This is due to the fact that at low particle concentrations, the particle volumetric concentration affects the nanofluid viscosity negligibly while the effect of temperature increases in the nanofluid thermo-physical properties.
With respect to the heat transfer performance and pressure drop required in the system, the "Figure of Merit" concept is introduced as a measure of the heat transferred per unit pumping power (see Equation 14 for details). Figure 12c presents the computed figures of merit for various tube geometries, nanofluid concentrations and nanofluid types. It is found that apart from the higher heat transfer rate, the coiled tubes have lower figures of merit than those of the straight tube. This can be explained by the higher pressure drops required in the coiled systems (see Figure 12b). Among all coiled tubes tested, the conical spiral tube gives the highest figure of merit, followed by in-plane spiral and conical spiral tubes. Furthermore, for the straight tube, the addition of nanoparticles improves the figure of merit significantly, albeit it decreases with increasing concentration. For coiled tubes filled with nanofluids, on the other hand, the improvement of figure of merit is only shown at low particle concentration of 1% and then it drops lower than that of water when more nanoparticles are added. Clearly, these results suggest that one can add nanoparticle up to 1% volumetric concentration to water to enhance heat transfer performance in coiled tubes; higher nanoparticle concentrations are not recommended.

Concluding remarks
A computational study was conducted to investigate the laminar flow heat transfer performance of square cross section tubes, i.e., straight, conical spiral, in-plane spiral, and helical spiral, with water and two nanofluids. It is found that adding 1% nanoparticle volumetric concentration improves heat transfer performance and the figure of merit for all tubes. However, higher amounts of nanoparticles is not recommended. In-plane spiral tubes give better performance than other coiled tubes for nanofluids. Furthermore, Al 2 O 3 nanofluid gives slightly better heat transfer performance than CuO nanofluid in coiled tubes. Future study will evaluate various modeling approaches for nanofluid heat transfer, e.g., single-phase, two-phase mixture, Euler-Euler, and Euler-Lagrange models, in coils with respect to the effect of secondary flow to the nanoparticle concentration.