Abstract

The additives of solid nanoparticles to liquids are significant enhancement of heat transfer and hydrodynamic flow. In this study, the thermal properties of three types of nanoparticles (Al2O3, TiO2, and SiO2) dispersed in water as a base fluid were measured experimentally. Forced convection heat transfer turbulent flow inside heated flat tube was numerically simulated. The heat flux around flat tube is 5000 W/m2 and Reynolds number is in the range of to . CFD model by finite volume method used commercial software to find hydrodynamic and heat transfer coefficient. Simulation study concluded that the thermal properties measured and Reynolds number as input and friction factor and Nusselt number as output parameters. Data measured showed that thermal conductivity and viscosity increase with increasing the volume concentration of nanofluids with maximum deviation 19% and 6%, respectively. Simulation results concluded that the friction factor and Nusselt number increase with increasing the volume concentration. On the other hand, the flat tube enhances heat transfer and decreases pressure drop by 6% and −4%, respectively, as compared with circular tube. Comparison of numerical analysis with experimental data available showed good agreement with deviation not more than 2%.

1. Introduction

The increasing demand for more efficient heat transfer fluids in many applications led to enhance heat transfer to meet the cooling challenge necessary such as the electronics, photonics, transportation, and energy supply industries [1]. The thermal conductivity of metallic liquids is much greater than that of nonmetallic liquids. Therefore, the thermal conductivities of fluids that contain suspended solid metallic particles could be expected to be significantly higher than those of conventional heat transfer fluids as shown in Table 1. Sen Gupta et al. [2] carried out experimental study to measure thermal conductivity in graphene nanofluids using the transient hot-wire method. They also studied carbon nanotube (CNT) and graphene oxide nanofluids. The magnitude of enhancement was between CNT and metallic/metal oxide nanofluids. Multicurrent hot-wire technique was used by Peñas et al. [3], to measure thermal conductivity of SiO2 and CuO in water and ethylene glycol nanofluids with concentrations up to 5% in mass fraction. They found a good agreement within 2%, with published thermal conductivities of the pure fluids. Optimizing thermal conductivity of nanofluids proposed by [4], there were various nanoparticles involving Al2O3 with different sizes, SiC with different shapes, MgO, ZnO, SiO2, Fe3O4, TiO2, diamond, and carbon nanotubes. The base fluids used deionized water, ethylene glycol (EG), glycerol, silicone oil, and the binary mixture of DW and EG. Results showed that the thermal conductivity enhancements of nanofluids could be influenced by multifaceted factors including the volume fraction of nanoparticles, the tested temperature, thermal conductivity of the base fluid, nanoparticles size, pretreatment process, and the additives of the fluids. The viscosity and specific heat of silicon dioxide (SiO2) nanoparticles with various diameters (20, 50, and 100 nm) in a 60 : 40 (by weight) ethylene glycol and water mixture were investigated experimentally by Namburu et al. [5]. Results showed new correlation from experimental data, which related viscosity with particle volume percent and nanofluid temperature. Also, specific heat of SiO2 nanofluid for various particle volume concentrations was presented. Bahiraei et al. [6] examined the effect of temperature and volume fraction on the viscosity for TiO2-water nanofluid. Results were recorded and analyzed within a temperature range of 25 to 70°C and volume fraction 0.1, 0.4, 0.7, and 1%. Viscosity was measured experimentally by [7] using a rheometer. It was obtained as a function of the nanoparticles mass fraction and shear rate. Water was already used as a base fluid and two different materials single wall carbon nanohorn (SWCNH) and titanium dioxide (TiO2). Results proposed empirical correlation equations of viscosity.

Wide variety of practical and industrial applications of forced convection in a flat tube as heat exchangers, heating processes, power generation, chemical processes, microelectronics, and cooling processes led to interest in this type of studies [816]. Duangthongsuk and Wongwises [17] found the heat transfer coefficient and friction factor of the nanofluid TiO2-water flowing in a horizontal double-tube counterflow heat exchanger under turbulent flow conditions experimentally and showed that the heat transfer coefficient of nanofluid is higher than that of the base liquid. Forced convection turbulent flow of Al2O3-water nanofluid inside an annular tube with variable wall temperature was investigated experimentally by Prajapati [18]; the results showed the enhancement of heat transfer due to the nanoparticle presence in the fluid. The forced convection flow between two corrugated cylinders was studied by Kittur [19]; results found friction factor and heat transfer on the boundaries. Horizontal double-tube heat exchanger counter turbulent flow was studied numerically by Bozorgan et al. [20]; Al2O3-water nanofluid of 7 nm with volume concentrations up to 2% was selected as a coolant use; the results showed that the pressure drop of nanofluid is slightly higher than water and increases with the increase in volume concentrations. Forced convection of nanofluid of a double-tube counterflow heat exchanger using CFD simulation FLUENT software was investigated by Demir et al. [21]. A double-tube coaxial heat exchanger heated by solar energy using Aluminum oxide nanofluid was presented experimentally and numerically by Luciu et al. [22]; results showed that nanofluids have a higher performance of heat transfer than base fluid.

In this study thermal conductivity and viscosity of three types of nanofluids (Al2O3, SiO2, and TiO2) in water are measured experimentally. Simulation study by using commercial software to solve governing equation of forced convection heat transfer through flat tube was conducted. The data measured was used as inputs parameters, and friction factor and heat transfer coefficient represented output parameters. Results were compared with experimental data available in the literature.

2. Thermal Properties Measurements

2.1. Density and Specific Heat

Many of researchers used regression equations of density () and specific heat () [23] as

2.2. Thermal Conductivity

Nanofluids were prepared in thermodynamics laboratory of Mechanical Engineering Faculty University Malaysia Pahang. Nanopowders were purchased from US Research Nanomaterials, Inc. (NovaScientific Resources (M) Sdn. Bhd.). They represented three types of commercial nanoparticles (Al2O3, TiO2, and SiO2) as shown in Figure 1 dispersed in water as a base fluid. Water was prepared in a laboratory by double distillation before using it for the experiments. The nanoparticles have size diameters of 13, 30, and 30 nm, respectively. Measured quantities of nanoparticles were dispersed in distilled water to obtain mass concentration nanofluids. The aid of mechanical stirrers was used to achieve a homogenously dispersed solution. This method depended on [2427] and was then subjected to ultrasonics for at least 3 hrs to break up any residual agglomerations. Mass of nanoparticles () and water () was measured with the accuracy of 0.001 g to estimate weight percentage () using [25]

Equation (3) was used to estimate the volume concentration of nanofluid depending on nanoparticles density () and base fluid density () at 25°C:

Sedimentation of nanoparticles at the bottom of the samples led to change of physical properties of the bulk nanofluids with time [26, 27]. The measurement of the thermal properties of the nanofluids requires many individual measurements for at least one month, so should take to check sample stability. Samples were checked after finishing each test but no visible sedimentation was found.

The transient hot-wires method as shown in Figure 1 was used to measure thermal conductivity of nanofluids experimentally. The wire placed along the axis of the container which was surrounded by the fluid whose thermal conductivity is to be measured. Platinum has high electrical resistivity that is,  Ωm (at 20°C) and order of magnitude higher than that of other metals. Also, it has a temperature coefficient of resistance of 0.0003925°C−1 (for pure platinum), which is much higher than that of other metals chosen as material of wire. The wire is to be used as a line heat source, so the wire diameter was kept within 100 μm. The length of the wire was kept to just a few centimeters, which compared to the wire diameter represents an infinitely long line heat source, assuring one directional (radial) heat transfer. Calibration method was used with standard fluid (glycerin) which brought with devices already; the error among reading data and the standard are 0.0023. After that verification was performed by using the pure liquid (water) and comparing with standard and the error among them is 0.0014. Thermal conductivity values were estimated by [2830] where represents the shape factor which was equal to 3 for spherical nanoparticles. The standard deviation corresponding to the series of individual data measured for each nanofluid was in all cases less than 0.15%.

2.3. Viscosity

To evaluate the thermal properties of nanofluids the viscosity is a significant indication. A commercial Brookfield DV-I prime viscometer was used to measure viscosity at different temperatures and rotor RPMs was shown in Figure 1. Base fluid (water) was used to measure viscosity for calibration; after that nanofluids were used to measure viscosity. Viscosity of nanofluid () was determined from the well-known Einstein equation for estimating viscosity, which is validated to spherical particles and volume concentration less than 5.0 vol. % and defined as [3133]

3. Computational Method

3.1. Physical Model

Cylindrical geometry coordinates of problem undertaken in Figure 2 are shown. Dimensions of the flat tube are major and minor diameter ( mm, mm); the length () and hydraulic diameter () of the flat tube are 500 mm and 4.68 mm. Reynolds number was calculated regarding hydraulic diameter () as [34]

The problem undertaken was assumed to be two-dimensional, steady, incompressible and Newtonian turbulent fluid flow, constant thermophysical properties of nanofluid, no effect of gravity and heat conduction in the axial direction and wall thickness of tubes neglected.

3.2. Governing Equations

Infinitesimal (less than 100 nm) solid particles were assumed to be able to use single phase approach, so single phase approach was adopted for nanofluid modeling. For all these assumptions, the dimensional conservation equations for steady state mean conditions are as follows: continuity, momentum, and energy equations Bejan [35]

High Reynolds number as input parameter was estimated; pressure treatment adopted SIMPLE scheme and turbulent viscous model was employed; converged solutions were considered for residuals lower than 10−6 for all the governing equations. The results of simulation for nanofluid were compared with the equations of Blasius (8) for friction factor and Dittos-Boelter equation (9) for Nu as [36]:

3.3. Boundary Conditions

Volume concentration nanofluids (1, 1.5, 2, and 2.5%) at 25°C base temperature were used for nanofluid as input. For comparison purposes, water was employed as the working fluid. CFD studies were carried out with uniform velocity profile at the inlet and pressure outlet condition used at the outlet of the flat tube. Turbulent intensity () was specified for an initial guess of turbulent quantities ( and ). The turbulent intensity calculated for each case based on the formula

The walls of tube were assumed to be perfectly smooth and the constant heat flux condition were specified on the inside tube wall with a value of 5000 W/m2. Reynolds number varied from to at each step of iterations as input data. The friction factor and Nu were introduced as output data.

3.4. Grid Independence Test

Grids independence in GAMBIT software for flat tube as cells and , subdivisions in the axial length, and surface face, respectively, were tested. To find the most suitable size of mesh faces, grid independent test was performed for the physical model. In this study, rectangular cells were used to mesh the surfaces of the tube wall but triangular cells were used to mesh the surfaces of gap as shown in Figure 3. The grid independence was checked by using different grid systems and four mesh faces were considered (, and ) for pure water. Nusselt number was estimated for all four mesh faces and results were proper. However, any number of mesh faces for these four cases can be used, but in this study, mesh faces with were adopted as the best in terms of accuracy as shown in Figure 4.

3.5. CFD Simulation

CFD simulations used FLUENT software with solver strategy. Gambit software was used to analyze problems. To make possible numerical solution of governing equations, single phase conservation equations were solved by control volume approach then by converting them to a set of algebraic equations. Simulation results were tested by comparing the predicted results [8, 17, 49] that used circular heated tube in experimental work. FLUENT software was used in CFD analysis in the literature and a detailed description of the mathematical model can be found in the FLUENT User’s Guide [50]. CFD modeling region could be classified into few major steps: preprocess stage, the geometry of problem undertaken is constructed as flat narrow and computational mesh generated in GAMBIT. It followed by the physical model that boundary conditions and other parameters were appropriately defined in models setup and solving stage. All scalar values and velocity components of the problem were calculated at the center of control volume interfaces where the grid schemes were used intensively. Throughout the iterative process accurate monitoring of the residuals was done. When the residuals for all governing equations were lower than 10−6, all solutions were assumed to be converged. Finally, the results could be obtained when FLUENT iterations lead to converged results defined by a set of converged criteria. The friction factor and Nusselt number inside elliptical tube could be obtained throughout the computational domain in the postprocess stage.

4. Results and Discussion

The thermal conductivity always increased when nanopowders addition on base fluid to prepare nanofluid. Thermal conductivity of (Al2O3, TiO2, and SiO2) nanofluids at 25°C were measured from 1% to 2.5% volume concentration. Figure 5 indicated thermal conductivity of nanofluid with volume concentration at temperature 25°C. Al2O3 nanofluid has the highest values of thermal conductivity followed by TiO2 and SiO2 suspended in water. The nanoparticles have significant thermal conductivity enhancement in water [51]. The theoretical data of thermal conductivity as (4) was indicated for pure water as solid black line to compare with data measured. The experimental data of other researchers of was indicated to validate data measured, and there is good agreement with deviation no more than 2% that may be related to different nanoparticle size [51]. The enhancement () of thermal conductivity represented percentage deviation of thermal conductivity of nanofluids against base fluid and not more than 20% as compared with base fluid. Experimental data of viscosity of nanofluids and base fluid was measured from 1% to 2.5% volume concentration at 25°C. The expected values of the shear rate range from 200 to 1300 s and all measured data is close to this range and the percentage average deviation is about 4%. The experimental data of viscosity was presented in Figure 6 with volume concentrations, theoretical data of viscosity as (5) showed for pure water as solid black line and the experimental data of other investigators to validate data measured. The viscosity decreased () represented percentage deviations among experimental data of nanofluids and base fluid in the the range 1.2 to 6%. The increase in viscosity with the increase of volume concentration was due to increase in percentage deviation which might be related to not using any surfactant or chemical additives when preparing nanofluids [52]. It seems Al2O3 that nanofluid has the highest values of viscosity followed by SiO2 and TiO2; finally pure water has the lowest viscosity. Nanofluids showed an increase in viscosity with the increase in volume concentration, similar behavior of base fluid. As compared to other investigators for both the thermal conductivity and viscosity data, the results indicated that there was difference among measured data and correlations; the reason may be due to various parameters such as particle preparation (different condition and accuracy) [53], particle size (the various size due to different thermal properties), the measurement technique, or even different particle sources [31, 32]. The effect of different volume concentration on friction factor is shown in Figures 7, 8, and 9 at the range of 1–2.5% nanoparticle volume concentration. The results revealed that the nanofluid with highest concentration of volume 2.5% has the highest friction factor at all Reynolds numbers. The reason of increasing of friction factor with the increasing of concentration of volume of nanofluid is the increasing of fluid viscosity which reduced the moving of fluid. The friction factor values of SiO2-water nanofluid appear higher than those of other types of nanofluid followed by Al2O3 and TiO2 in water. On the other side, the CFD analysis of pure water inside flat tube was compared with Blasius equation (8) for water inside circular tube which was indicated as black solid line. It appears that the values of friction factor at flat tube are less than that at circular tube, and this was achieved to use flat tube to decrease pressure drop. Figures 10, 11, and 12 show Nusselt number with Reynolds number at (1–2.5%) volume concentrations of nanofluid at 25°C base temperature. Results showed that the nanofluid has the highest Nusselt number with the highest volume concentration 2.5% at all Reynolds numbers. The reason of increasing of Nusselt number with the increasing of concentration of volume of nanofluid is the increasing of fluid thermal conductivity. The Nusselt number of SiO2 is the highest value followed by Al2O3 and TiO2 and the cause may be related to the fact that SiO2 nanofluid has the lowest thermal conductivity than other nanofluids and the highest average velocity among the fluids due to lowest density compared with the others. The fluid velocity plays an important role on the heat transfer in case of forced convection and it represents the main reason to give high heat transfer coefficient [49]. The irregular and random movements of particles increase the energy exchange rates in the fluid with penalty on the wall shear stress and consequently enhance the thermal dispersion of the flow. The CFD analysis of Nusselt number for pure water inside flat tube was compared with Dittus-Boelter equation (9) for water inside circular tube was shown as black solid line. It seems that the values of Nusselt number at flat tube are higher than circular tube referring to enhance heat transfer. The CFD analysis of friction factor and Nusselt number for Al2O3, TiO2, and SiO2 was compared with experimental data of Sunder and Sharma [8], Duangthongsuk and Wongwises [17], and Mohammed et al. [49], respectively, as shown in Figure 13. It shows good agreement among CFD analysis and experimental data. Figure 14 illustrated velocity and temperature profiles at different Reynolds number. It seems that the maximum and minimum velocity is 0.135 and 0.958 m/s at Reynolds numbers 5000 and 50000, respectively. Furthermore the minimum temperature is 43 and 54°C at Reynolds number is 5000 and 50000 respectively. The reason may be related to the low Reynolds number that is given more time to make deference temperature between input and output of tube. The regression equations among input and output parameters have been analyzed by using Minitab 16 software. The input parameters represented that Reynolds number as a velocity inlet and Prandtl number as a concentration of nanofluid; on the other hand, the output parameters represented friction factor and Nusselt number. Figure 15 showed friction factor and Nusselt number from CFD analysis against (8) and (9) data to correlate among parameters.

These equations are concluded under and ; the regression equation of friction factor has been expected as Elsewhere, the regression equation of Nusselt number has been found as:

5. Conclusions

In the present study, thermal properties of three types of nanoparticles were suspended in water measured experimentally in thermal laboratory. Thermal conductivity and viscosity were measured experimentally between 1% and 2.5% volume concentrations at 25°C. Results of thermal conductivity showed that the increase in volume concentrations and temperatures is due to the increase in thermal conductivity of nanofluids with maximum deviation (19%) for nanoparticles dispersed in water. The study showed that Al2O3 nanofluid has the highest values of thermal conductivity followed by TiO2 and SiO2 suspended in water. The data measure also showed that the viscosity of nanofluids significantly increases with increasing particle volume concentration with maximum deviation (6%) as compared with water. The results showed that Al2O3 nanofluid has the highest values of viscosity followed by SiO2 and TiO2; finally pure water had the lowest viscosity. The proposed correlations for predicting the thermal conductivity and viscosity of nanofluids showed good agreement with the data measured and the experimental results of [32, 33]. Forced convection heat transfer under turbulent flow by numerical simulation with uniform heat flux boundary condition around flat tube was studied. The heat transfer enhancement resulting from various parameters such as nanoparticle concentration of volume and Reynolds number were reported. The finite volume methods was used to solve the governing equations with certain assumptions and appropriate boundary conditions. The Nusselt number and friction factor were obtained through the numerical simulation. The study concluded that the enhancement of friction factor and Nusselt number is −4% and 6% for enhanced tube more than that of the circular tube at all Reynolds numbers. The (2.5%) volume concentration of nanofluid has the highest friction factor values, followed by 2, 1.5, and 1%. The CFD analysis of pure water friction factor at flat tube has lower values than that of circular tube which estimated from Blasius equation. The Nusselt number of SiO2 is the highest value followed were by Al2O3 and TiO2. There is a good agreement among the CFD analysis of friction factor and Nusselt number of nanofluid with experimental data of Sunder and Sharma [8], Duangthongsuk and Wongwises [17], and Mohammed et al. [49] with deviation being not more than 2%. The regression equations were found under and , for friction factor and Nusselt number have been concluded as in (11) and (12), respectively, where the constants of friction equation (11) Nusselt number equation (12) were tabulated as shown in Tables 2 and 3.

Acknowledgment

The financial support to the authors by University Malaysia Pahang is gratefully acknowledged.