Modern Mechanical Engineering, 2012, 2, 130-136 http://dx.doi.org/10.4236/mme.2012.24017 Published Online November 2012 (http://www.SciRP.org/journal/mme)

Numerical Study on Application of CuO-Water Nanofluid in Automotive Diesel Engine Radiator

Navid Bozorgan1*, Komalangan Krishnakumar2, Nariman Bozorgan1 1Department of Mechanical Engineering, Abadan Branch, Islamic Azad University, Abadan, Iran

2Department of Mechanical Engineering, College of Engineering, Trivandrum, India Email: *N.Bozorgan@gmail.com

Received August 15, 2012; revised September 17, 2012; accepted September 29, 2012

ABSTRACT

Application of CuO-water nanofluid with size of the nanoparticles of 20 nm and volume concentrations up 2% is nu- merically investigated in a radiator of Chevrolet Suburban diesel engine under turbulent flow conditions. The heat transfer relations between airflow and nanofluid coolant have been obtained to evaluate local convective and overall heat transfer coefficients and also pumping power for nanofluid flowing in the radiator with a given heat exchange ca- pacity. In the present study, the effects of the automotive speed and Reynolds number of the nanofluid in the different volume concentrations on the radiator performance are also investigated. The results show that for CuO-water nanofluid at 2% volume concentration circulating through the flat tubes with Renf = 6000 while the automotive speed is 70 km/hr, the overall heat transfer coefficient and pumping power are approximately 10% and 23.8% more than that of base fluid for given conditions, respectively. Keywords: Automotive Diesel Engine Radiator; CuO-Water Nanofluid; Turbulent Flow; Heat Transfer; Pumping

Power

1. Introduction

Cooling is one of the top technical challenges to obtain the best automotive design in multiple aspects (perform- ance, fuel consumption, aesthetics, safety, etc.). Automo- tive radiator is an important part of the engine cooling system. Due to limited space at the front of the engine, the size of the radiator is restricted and cannot be essen- tially increased. Therefore, it is necessary to increase the heat transfer capabilities of working fluids such as water and ethylene glycol in radiators because of their low thermal conductivity. A recent advancement in nanote- chnology has been the introduction of nanofluids, i.e. colloidal suspensions of nanometer-sized solid particles instead of common working fluids. Nanofluids were first innovated by Choi [1] in 1995 at the Argonne National Laboratory, USA. Compared with traditional solid-liquid suspensions containing millimeter or micrometer sized particles, nanofluids as coolants in the heat exchangers have shown better heat transfer performance because of small size of suspend solid particles. It causes that nano- fluids have a behavior similar to base liquid molecules. Because of their excellent thermal performance, nano- fluids have attracted attention as a new generation of heat transfer fluids in building heating, in heat exchangers, in

chemical plants and in automotive cooling applications. Recently there have been considerable research findings highlighting superior heat transfer performances of nan- ofluids. Peyghambarzadeh et al. [2] investigated experi- mentally the convective heat transfer enhancement of water and EG based nanofluids consisting of Al2O3 nanoparticles up 1% volume concentration as the cool-ants inside flat aluminum tubes of the car radiator under laminar and turbulent flows. The results show that the heat transfer enhances about 40% compared with the base fluids in the best conditions. Ollivier et al. [3] in-vestigated the use of nanofluids as a jacket water coolant in a gas spark ignition engine. They numerically simu-lated the unsteady heat transfer through the cylinder and inside the coolant flow and reported that because of higher thermal diffusivity of nanofluids, the thermal sig-nal variations for knock detection increased by 15% over the predicted using pure water. Vajjha et al. [4] numeri-cally investigated the heat transfer augmentation by ap-plication of two different nanofluids consisting of Al2O3 and CuO nanoparticles in an ethylene glycol and water mixture circulating through the flat tubes of an automo-bile radiator. The numerical results showed that at a Reynolds number of 2000, the percentage increase in the average heat transfer coefficient over the base fluid for a 10% Al2O3 nanofluid was 94% and that for a 6% CuO *Corresponding author.

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL. 131

nanofluid was 89%. Also the average heat transfer coef- ficient increases with the Reynolds number and also with the particle volumetric concentration. Leong et al. [5] have studied the application of nanofluids as working fluids in shell and tube heat recovery exchangers in a biomass heating plant and showed that about 7.8% of the heat transfer enhancement could be achieved with the addition of 1% copper nanoparticles in ethylene glycol based fluid at 26.3 kg/s and 111.6 kg/s mass flow rate for flue gas and coolant, respectively. Ijam et al. [6] theo- retically analyzed a minichannel heat sink with a 20 × 20 cm bottom for SiC-water and TiO2-water nanofluids as the coolants through hydraulic diameters under turbulent flow. Their results showed a maximum enhancement of 12.44% in thermal conductivity for SiC-water and 9.99% for TiO2-water at 4% of volume fraction. Also the maxi- mum improvement in heat flux by using SiC-water and TiO2-water nanofluids at 4% of volume fraction for inlet velocities of 2 and 6 m/s is calculated by ~7.63%, 12.43% and 7.25%, 12.43%, respectively. In 2012 Saeed- inia et al. [7] applied CuO-base oil particles varying in the range of 0.2% - 2% inside a circular tube and showed that the CuO nanoparticles suspended in base-oil in-creases the heat transfer coefficient even for a very low particle concentration of 0.2% volume concentration. They found a maximum heat transfer coefficient en- hancement of 12.7% for 2% CuO nanofluid. Shafahi et al. [8] used a two-dimensional analysis to study the thermal performance of a cylindrical heat pipe utilizing Al2O3, CuO and TiO2 nanofluids. Their results confirmed that the thermal performance of a heat pipe is improved and temperature gradient along the heat pipe and thermal resistance across the heat pipe are reduced and maximum capillary heat transfer of the heat pipe is observed when nanofluids are utilized as the working fluid.

In the present paper, 20 nm-CuO/water nanofluid with concentration up 2 vol% has been numerically studied as a coolant in a radiator of Chevrolet Suburban diesel en- gine with a given heat exchange capacity. It shall be noted that metal oxides such as CuO nanoparticles are chemically more stable than their metallic counterparts.

2. Analysis

In order to investigate the heat transfer performance of nanofluids and use them in practical applications, it is necessary first to study their thermophysical properties such as density, specific heat, viscosity and thermal con- ductivity. In this study, to validate the numerical results, the methodology employed to determine the efficiency of CuO-water nanofluid as coolant in a radiator of Chevro- let Suburban diesel engine with a given heat exchange capacity is based on new rheological and thermophysical properties correlations developed from experiments.

2.1. Prediction of Thermophysical Properties

At the first step, the heat characteristics of the nanofluid have been evaluated and at the next step the application of nanofluid coolant has been investigated to improve the heat transfer performance of the radiator. This radiator consists of 644 brass flat tubes and 346 continuous cop- per fins [9,10]. The main characteristics of radiator are listed in Table 1 that useful for assessing the radiator performance in this work. However, following assump- tions are made: The flow is an incompressible, steady state and tur-

bulent. The effect of body force is neglected. The thermophysical properties of nanofluids are con-

stant. The characteristics of nanoparticles and base fluid

used in this study are summarized in Table 2. The ne- cessary thermophysical properties in this paper are den-sity, viscosity, specific heat and thermal conductivity.

In this paper, density (ρnf) and special heat capacity (cp,nf) of CuO/water nanofluid have been calculated based on empirical correlations proposed by Pak [11] and Xuan [12] as follows:

1nf bf p (1)

,,

1 bf p bf p p pp nf

nf

c cc

, (2)

where

is nanoparticle volume concentration and ρp, ρbf and cp,p, cp,bf are the densities and the specific heats of

Table 1. Radiator specifications.

Description Air Coolant

Radiator frontal area, (h × w) 0.6 m × 0.5 m

Fin thickness, t 0.01 cm

Fin length for symmetric heating from primary to midpoint between plates, L

0.4 m

Hydraulic diameter, Dh 0.351 cm 0.373 cm

Heat transfer area/total volume, α 886 138

Fin area/total area, (Af/A) 0.845

Table 2. Thermophysical properties of base fluid and nanoparticles.

Property CuO Water

cp [J·kg−1·K−1] 540 4197

ρ [kg·m−3] 6510 971

k [Wm-1·K−1] 18.0 0.669

dp (nm) 20 0.384

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL. 132

the nanoparticles and base fluid, respectively. Also, thermal conductivity (knf) and viscosity (μnf) for

nanofluid have been estimated based on two semi-em- pirical equations presented by Corcione [13] in 2011 on the basis of a wide variety of experimental data available in the literature as following equations:

10 0.03

0.4 0.66 0.661 4.4Re Preff pbf

bf fr bf

k kT

k T k

(3)

0.3 1.031 34.87

bfnf

p bfd d

(4)

where kbf is the thermal conductivity of the base fluid, Re is the nanoparticle Reynolds number, Prbf is the Prandtl number of the base fluid, T is the nanofluid temperature, Tfr is the freezing point of the base fluid, kp is the thermal conductivity of the nanoparticles, μbf is the dynamic vis- cosity of the base fluid, dp is the diameter of the nanopar- ticles and dbf is the equivalent diameter of a base fluid molecule. The standard deviations of error of Equations (3) and (4) are 1.86% and 1.84%, respectively.

In more detail, the Reynolds number of the suspended nanoparticles can be calculated as follows [13]:

2

2Re

πbf b

bf p

k T

d

(5)

where in kb = 1.38066 × 10−23 J/K is the Boltzmann con- stant.

2.2. Heat Transfer and Pressure Drop Modeling

2.2.1. Heat Transfer Modeling The rate of heat transferred between nanofluid coolant and airflow in the radiator can be written as follows:

, 1 2 air ,air 2 1nf p nf pQ m c T T m c t t (6)

where nf and air denote the relevant parameters of nan- ofluid coolant and airflow. Here, T1 = 110˚C and T2 = 49˚C are the inlet and outlet temperatures of nanofluid coolant, t1 = 32˚C is the inlet air temperature and outlet air temperature (t2) can be calculated as:

2 1air ,airp

Qt t

m c

(7)

The mass flowrates are calculated based on the pump flowrate for water and the speed and frontal area for the air as follows:

.nf nf vol p pm Q (8)

air air Automotivem V h W (9)

Here, the Bosch Cobra pump is used with flowrate of Qp = 317 gallons/hour and volumetric efficiency (vol.p)

by 80%. In Equation (9), VAutomotive is the automotive speed.

The total heat transfer area of the radiator, A, is com- puted from the following equation:

QA

U LMTD

(10)

1 1 2 2

1 1

2 2

ln

T t T tLMTD

T t

T t

(11)

In Equation (10), U is the overall heat transfer coefficient which is expressed in the following form [9]:

1

air

air

1 1

nfonf

U Rh

h

f (12)

where hair is the heat transfer coefficient of air, o is the total surface temperature effectiveness, hnf is the heat transfer coefficient of nanofluid coolant and Rf is the fouling resistance that is assumed to be 5 × 10−4 m2·KW−1. In the present study, the heat transfer coefficient of air is calculated as:

0.14

,air 2 3

,

Re

Pra a a p a nf

wnfa h a

J Ch

D

(13)

In the above equation,

0.14

nf

wnf

is the viscosity cor-

rection factor that is defined as the ratio of viscosity of the nanofluid at the mean temperature of inlet and outlet conditions to that one at the mean temperature of wall tube and also cp,air = 1007 J·kg−1·K−1, Pra = 0.7057, μa = 0.00001889 (Nsm−2)

In Equation (13), Reynolds number and Colburn factor for the air are calculated with considering the air proper- ties as follows:

,Re a h aa

a

G D

(14)

0.383

0.174

Reaa

J (15)

where Ga is the mass flow rate of air. In Equation (12), the total surface temperature effectiveness is expressed as:

1.0 1.0 fo f

A

A (16)

where f is the fin efficiency of plate fin which is defined as follows:

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL. 133

,πt h nfA D H (27) air

air

2tan

2f

hh L

kt

hL

kt

(17)

Here, k = 377 Wm−1·K−1 is thermal conductivity of copper fin.

The heat transfer coefficient of the nanofluid coolant (hnf) flowing in the flat tubes can be calculated with con- sidering the turbulent Nusselt number presented by Li and Xuan [14] as follows:

, 0.6886 0.001

0.14

0.9238 0.4

0.0059 1.0 7.6286 Pe

Re Pr

nf h nfnf d

nf

nfnf nf

wnf

h DNu

k

(18)

where Ped is the nanofluid Peclet number and is defined in the following form:

Pe nf pd

nf

u d

(19)

where αnf is the nanofluid thermal diffusivity which is defined as follows:

,

nfnf

nf p nf

k

c

(20)

The Reynolds and Prandtl numbers in Equation (18) are calculated with considering the nanofluid properties as follows:

,Re nf nf h nfnf

nf

u D

(21)

,Pr p nf nfnf

nf

c

k

(22)

2.2.2. Pressure Drop Modeling In the present paper, the fanning friction factor (fn), the pressure drop (Δpnf) and pumping power (PP) for CuO- water nanofluid coolant flowing in the flat tubes are cal- culated as follows [10]:

0.250.079 Renf nff

(23)

0.252

,

2 nf nf nfnf

nf h nf wnf

G f Hp

D

(24)

p nfPP Q p (25)

In Equation (24), H is the total nanofluid flow length which is obtained as:

where At is the total heat transfer area of flat tubes.

3. Simulation Results and Discussion

3.1. Local Heat Transfer

In Figure 1 the ratio of thermal conductivity and heat transfer coefficient of nanofluid in comparison with the base fluid (kr and hr) emphasize the local heat transfer augmentation. This analysis is carried out at Renf = 6000. This figure reveals that as the concentration increases, the effect of increasing nanoparticles concentration on changing the thermal conductivity is lower than changing the heat transfer coefficient. As an example, for a volume concentration of 2 vol%, the local heat transfer coeffi- cient ratio is about hr = 2.05 while the thermal conductiv- ity ratio is kr = 1.55. Figure 2 shows the effect of nanoparticles concentration on the heat transfer coeffi- cient and Nusselt number. Results show that the heat transfer coefficient and Nusselt number can be enhanced by adding nanoparticles to the base fluid. Enhancement of heat transfer by the nanofluid may be resulted from the following two aspects: first the incremented particles which increase the thermal conductivity of the mixture; second, chaotic movement of ultrafine particles acceler- ates the energy exchange between the fluid and the wall of the heat exchanger. However, it should be noted that increasing the particles concentration raises the fluid viscosity and decreases the Reynolds number and con- sequently decreases the local heat transfer behavior. But the results shown in Figures 1 and 2 indicate that in- creasing in particles concentration raises the heat transfer coefficient. This is due to the fact that the change in the thermal conductivity is more effective than the change in

Figure 1. Comparison between kr and hr for CuO-water nanofluid coolant in the radiator at different volume frac- tions. t fA A A (26)

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL. 134

Figure 2. The Local heat transfer coefficient and Nusselt number for CuO-water nanofluid coolant at various con- centrations. the fluid viscosity on heat transfer enhancement of nano- fluid at higher concentration.

3.2. Overall Heat Transfer

Figures 3-5 show the results of simulating the effects of using CuO-water nanofluid coolant with different vol- ume fractions flowing in the flat tubes in the various ranges of the automotive speed and Reynolds number of the nanofluid on the overall heat transfer coefficient. The results of this study show that nanofluid enhances the overall heat transfer coefficient in the radiator and there- fore the total heat transfer area of the radiator can be re- duced.

Figure 3 shows the overall heat transfer coefficient for CuO-water nanofluid coolant in an automotive radiator that has been calculated by Equation (12). In this analysis, Reynolds number of the nanofluid is 6000 and the automotive speed is assumed to be 70 km/hr. As shown in this figure, the overall heat transfer coefficient is high when the probability of collision between nanoparticles and the wall of the flat tubes has increased under higher concentration conditions. It confirms that nanofluids have considerable potential to use in the automotive ra- diator. A further inspection of Figure 3 shows that the overall heat transfer coefficient of the CuO-water nano- fluid for volume concentrations in the range of 0.1% to 2% is about 122 and 134.27 W/m2k, respectively.

Figure 3. Variations of overall heat transfer coefficient for nanofluid at different volume fractions.

0%0.5%1.0%1.5%2.0%

Figure 4. Effect of coolant Reynolds number on overall heat transfer coefficient in the radiator.

0%0.5%1.0%1.5%2.0%

The coolant Reynolds number has an important role in

determining the radiator performance. The coolant Rey- nolds number must be properly controlled by thermo- stating pump because engine might be overcooled or overheated. Therefore, the effect of the coolant Reynolds number on the radiator performance at an automotive speed of 70 km/hr has been studied in this paper. Ac- cording to Figure 4, the increase in the coolant Reynolds number increases the overall heat transfer coefficient that

Figure 5. Effects of the automotive speed and CuO volume fraction on overall heat transfer coefficient in the radiator. the magnitude of this property for CuO-water nanofluid at higher concentration is more than that of a base fluid. For instance, the results show that the increase in overall heat transfer coefficient for 2% CuO-water nanofluid at

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL. 135

Renf = 8000 compared with Renf = 6000 is about 2.3%. As well as, the automotive speed plays an important

role in the radiator performance. As expected, the mass flowrate of the air increases with increasing the automo- tive speed and therefore Reynolds number and mass ve- locity for the air increase. Increasing the air Reynolds number raises the air mass velocity (refer to Equation (14)) and decreases the air Colburn factor (refer to Equa- tion (15)). But the results shown in Figure 5 indicate that increasing in automotive speed and air Reynolds number raises the overall heat transfer coefficient because of the increase in air heat transfer coefficient. Therefore, it can be concluded that the change in the air mass velocity is more than the change in Colburn factor with increasing the air Reynolds number. In the analysis shown in Figure 5, the overall heat transfer coefficient of the radiator is calculated for the fixed coolant Reynolds number of 6000 at different volume fractions. For instance, the re- sults show that the increase in overall heat transfer coef- ficient for 2% CuO-water nanofluid at automotive speed of 80 km/hr compared with 60 km/hr is about 16%.

3.3. Pumping Power

In order to apply the nanofluids for practical application, in addition to the heat transfer performance it is neces- sary to study their flow features. With increasing nano- particles loading in the base fluid, viscosity and density of the nanofluids increase and therefore the friction fac- tor and the pressure drop must be increased. Hence, nan- ofluids generally require the greater pumping power than their base fluid. In the present paper, the pumping power for CuO-water nanofluid coolant flowing in the flat tubes in the various ranges of the coolant Reynolds number (at VAutomotive = 70 km/hr) and the automotive speed (at Renf = 6000) are calculated and shown in Figures 6 and 7. For instance, Figure 6 shows that the increase in pumping power for 2% CuO-water nanofluid at Renf = 8000 com- pared with Renf = 6000 is about 61.75%. Also as an ex- ample, Figure 7 shows that the decrease in pumping power for 2% CuO-water nanofluid at the automotive speed of 80 km/hr compared with 60 km/hr is about 16.7%.

4. Conclusion

This paper presented a numerical investigation of the use of CuO-water nanofluid as a coolant in a radiator of Chevrolet Suburban diesel engine with a given heat ex- change capacity. The local convective and overall heat transfer coefficients and pumping power for CuO-water nanofluid at different volume fractions (0.1% - 2%) was studied under turbulent flow conditions. Also, the effects of the coolant Reynolds number and the automotive speed on the radiator performance are considered in this

0%0.5%1.0%1.5%2.0%

Figure 6. Effect of coolant Reynolds number on pumping power in the radiator.

0%0.5%1.0%1.5%2.0%

Figure 7. Effect of the automotive speed on pumping power for nanofluid with different volume fractions in the radia- tor. work. The simulation results indicate that the overall heat transfer coefficient of nanofluid is greater than that of water alone and therefore the total heat transfer area of the radiator can be reduced. However, the considerable increase in associated pumping power may impose some limitations on the efficient use of this type of nanofluid in automotive diesel engine radiators.

5. Acknowledgements

The authors would like to express their appreciation to the Islamic Azad University of Abadan Branch for pro- viding financial support.

REFERENCES [1] S. U. S. Choi and J. A. Eastman, “Enhancing Thermal

Conductivity of Fluids with Nanoparticles,” ASME In-ternational Mechanical Congress and Exposition, San

Copyright © 2012 SciRes. MME

N. BOZORGAN ET AL.

Copyright © 2012 SciRes. MME

136

Francisco, 12-17 November 1995.

[2] S. M. Peyghambarzadeh, S. H. Hashemabadi., S. M. Ho- seini and M. S. Jamnani, “Experimental Study of Heat Transfer Enhancement Using Water/Ethylene Glycol Based Nanofluids as a New Coolant for Car Radiators,” International Communications in Heat and Mass Trans- fer, Vol. 38, No. 9, 2011, pp. 1283-1290. doi:10.1016/j.icheatmasstransfer.2011.07.001

[3] E. Ollivier, J. Bellettre, M. Tazerout and G. C. Roy, “De- tection of Knock Occurrence in a Gas SI Engine from a Heat Transfer Analysis,” Energy Conversion and Man- agement, Vol. 47, No. 7-8, 2006, pp. 879-893. doi:10.1016/j.enconman.2005.06.019

[4] R. S. Vajjha, D. K. Das and P. K. Namburu, “Numerical Study of Fluid Dynamic and Heat Transfer Performance of Al2O3 and CuO Nanofluids in the Flat Tubes of a Ra- diator,” International Journal of Heat and Fluid Flow, Vol. 31, No. 4, 2010, pp. 613-621. doi:10.1016/j.ijheatfluidflow.2010.02.016

[5] K. Y. Leong, R. Saidur, T. M. I. Mahlia and Y. H. Yau, “Modeling of Shell and Tube Heat Recovery Exchanger Operated with Nanofluid Based Coolants,” International Journal of Heat and Mass Transfer, Vol. 55, No. 4, 2012, pp. 808-816. doi:10.1016/j.ijheatmasstransfer.2011.10.027

[6] A. Ijam and R. Saidur, “Nanofluid as a Coolant for Elec- tronic Devices (Cooling of Electronic Devices),” Applied Thermal Engineering, Vol. 32, 2012, pp. 76-82. doi:10.1016/j.applthermaleng.2011.08.032

[7] M. Saeedinia, M. A. Akhavan-Behabadi and M. Nasr, “Experimental Study on Heat Transfer and Pressure Drop of Nanofluid Flow in a Horizontal Coiled Wire Inserted Tube under Constant Heat Flux,” Experimental Thermal and Fluid Science, Vol. 36, 2012, pp. 158-168. doi:10.1016/j.expthermflusci.2011.09.009

[8] M. Shafahi, V. Bianco, K. Vafai and O. Manca, “An In- vestigation of the Thermal Performance of Cylindrical Heat Pipes Using Nanofluids,” International Journal of Heat and Mass Transfer, Vol. 53, No. 1-3, 2010, pp. 376- 383. doi:10.1016/j.ijheatmasstransfer.2009.09.019

[9] D. G. Charyulu, G. Singh and J. K. Sharma, “Perfor- mance Evaluation of a Radiator in a Diesel Engine—A Case Study,” Applied Thermal Engineering, Vol. 19, No. 6, 1999, pp. 625-639. doi:10.1016/S1359-4311(98)00064-7

[10] K. Y. Leong, R. Saidur, S. N. Kazi and A. H. Mamun, “Performance Investigation of an Automotive Car Radia- tor Operated with Nanofluid-Based Coolants (Nanofluid as a Coolant in a Radiator),” Applied Thermal Engineer- ing, Vol. 30, No. 17-18, 2010, pp. 2685-2692. doi:10.1016/j.applthermaleng.2010.07.019

[11] B. C. Pak and Y. I. Cho, “Hydrodynamic and Heat Transfer Study of Dispersed Fluids with Submicron Me- tallic Oxide Particles,” Experimental Heat Transfer, Vol. 11, No. 2, 1998, pp. 151-170. doi:10.1080/08916159808946559

[12] Y. Xuan and W. Roetzel, “Conceptions of Heat Transfer Correlation of Nanofluids,” International Journal of Heat and Mass Transfer, Vol. 43, No. 19, 2000, pp. 3701-3707. doi:10.1016/S0017-9310(99)00369-5

[13] M. Corcione, “Empirical Correlating Equations for Pre- dicting the Effective Thermal Conductivity and Dynamic Viscosity of Nanofluids,” Energy Conversion and Man- agement, Vol. 52, No. 1, 2011, pp. 789-793. doi:10.1016/j.enconman.2010.06.072

[14] Q. Li and Y. Xuan, “Convective Heat Transfer and Flow Characteristics of Cu-Water Nanofluid,” Science in China Series E: Technological Sciences, Vol. 45, No. 4, 2002, pp. 408-416. doi:10.1360/02ye9047