SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Nanofluid impingement jet heat transfer

Obida Zeitoun* and Mohamed Ali

Author affiliations

King Saud University, Mechanical Engineering Department, King Abdullah Institute for Nanotechnology, 800 Riyadh, 11421, Saudi Arabia

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:139  doi:10.1186/1556-276X-7-139


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/139


Received:24 November 2011
Accepted:17 February 2012
Published:17 February 2012

© 2012 Zeitoun and Ali; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Experimental investigation to study the heat transfer between a vertical round alumina-water nanofluid jet and a horizontal circular round surface is carried out. Different jet flow rates, jet nozzle diameters, various circular disk diameters and three nanoparticles concentrations (0, 6.6 and 10%, respectively) are used. The experimental results indicate that using nanofluid as a heat transfer carrier can enhance the heat transfer process. For the same Reynolds number, the experimental data show an increase in the Nusselt numbers as the nanoparticle concentration increases. Size of heating disk diameters shows reverse effect on heat transfer. It is also found that presenting the data in terms of Reynolds number at impingement jet diameter can take into account on both effects of jet heights and nozzle diameter. Presenting the data in terms of Peclet numbers, at fixed impingement nozzle diameter, makes the data less sensitive to the percentage change of the nanoparticle concentrations. Finally, general heat transfer correlation is obtained verses Peclet numbers using nanoparticle concentrations and the nozzle diameter ratio as parameters.

Keywords:
nanofluid; heat transfer; impingement jet

Introduction

Fluid heating and cooling play very important roles in many industries including power generation, production processes, transportation and electronics. Heat transfer can be enhanced using different methods such as extended surfaces (fins), vibration of the heated surfaces, injection or suction of the fluid and applying electrical or magnetic fields. Nanofluid heat transfer is an innovative technology which can be used to enhance the heat transfer. The term nanofluid refers to a new kind of fluid produced by suspending nanoparticles in the base fluid.

Impinging liquid jet is an established technique to provide high local heat transfer coefficients between the impinged liquid and a surface. This cooling technique is considered as an attractive cost effective method of cooling [1,2]. Combining the liquid jet impingement and the nanofluid technologies is thought to capture the advantages of both and consequently enhances the heat transfer significantly. Enhancing the heat transfer means compact size and low weight which reduces the cooling system capital cost.

Liquid jets can be classified as submerged or free surface. A submerged jet is formed when a liquid jet is discharged into the same liquid medium. A free surface jet is formed when a liquid jet is discharged into a gas medium. For free surface jet, the liquid jet impingement has demonstrated high cooling capacity as reported by Liu and Lienhard [1], Lienhard and Hadeler [2], and Bergles [3]. The Jet impingement or free surface flow can be classified according to jet orientation, surface type or flow type as vertical or horizontal, flat or curved, and single or two phase flow, respectively.

Upon liquid jet impinging on a horizontal plate, the liquid spread due to inertia and gravity forces. A hydraulic jump occurs when the flow suddenly changes from shooting flow, before the hydraulic jump, to streaming flow after the hydraulic jump. This jump is accompanied by a sudden increase in the liquid film thickness. This change is accompanied by a significant decrease in the heat transfer. Jambunathan et al. [4] reviewed previous works of heat transfer between circular liquid jet and horizontal surface. Review of previous investigations reveals very high heat transfer rates at the stagnation point of the jet due to the very thin boundary layer at the stagnation point. However, at a distance of two to three nozzle diameters from the stagnation point, the cooling rate is less than half that of the stagnation value.

The pioneer work of Watson [5] has divided the flow field of an impinging jet over a horizontal circular disk before the hydraulic jump into two regions. The first region is near the center of the jet, where he assumed a boundary layer type flow. The second part was regarded as a free surface flow up until the hydraulic jump position. The region before the hydraulic jump is characterized by a high heat transfer capability. Chaudhury [6] followed Watson [5] and solved the energy equation, based on integral analysis, taking into consideration the viscous dissipation. The thermal boundary condition at the circular plate was assumed to be of the constant wall temperature type.

Experimental investigations of impingement jet heat transfer include works of Ishigai, Zhao and Masuoka, Baonga et al. and Teamah and Farahat [7-10]. Ishigai. [7] studied experimentally the flow and heat transfer of an impinging round jet over a horizontal plate. They compared their experimental measurements of the film thickness with the theoretical predictions of Watson [5]. There was a good agreement between the two solutions near the center of the jet, but as the radial location increases, the difference becomes wider. Zhao and Masuoka [8] have investigated flow and heat transfer due to liquid jet impingement on a circular surface. They have studied the heat transfer between small jets of 0.9 and 2 mm and a disk of 10 mm diameter. Baonga et al. [9] investigated liquid film, hydraulic jump and local heat transfer distributions along the radial direction of a circular disk. For jet Reynolds number in the range of 1,050 to 9,000, and for each nozzle diameter, the difference between the stagnation and average Nusselt numbers decreases significantly for higher Reynolds number. When the jet Reynolds number increases, the average heat transfer coefficient increases because of the increase in the liquid flow rate. Furthermore, Teamah and Farahat [10] have investigated both the heat transfer and fluid flow due to the impingement of vertical circular water jet on a horizontal heated surface numerically and experimentally. However, the hot surface used in their experiment was square of 0.95 m side.

There are few models available for the average liquid jet impingement Nusselt number along circular disks; such models are those given by Zhao et al. [8] and Teamah and Farahat [10]. Integral analysis method was used by [8] to solve the flow along the radial direction of the circular disk. The film along the radial direction was divided into stagnation, near impingement point, boundary layer and similarity regions. Based on integral analysis results, Zhao et al. [8] have developed the following model:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M1">View MathML</a>

(1)

where Nuj is the Nusselt number based on the jet nozzle diameter Dj, D is the disk diameter; Pr is the Prandtl number and Rej is the Reynolds number based on jet nozzle diameter. For the complete list of the symbols used in all equations, please refer to the 'Appendix' section. Teamah and Farahat [10] developed a mathematical two-dimensional model in radial and vertical directions. They solved the governing equations numerically up to r /Dj = 50 in the radial direction. However, they integrated heat transfer results in steps along the radial direction, from r = 0 into r to obtain average Nusselt number as a function of the radial positions. Based on their numerical data, they have obtained two correlations based on jet Reynolds number. For Rej range from 5,000 to 20,000:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M2">View MathML</a>

(2)

Different jet orientations were also investigated by Silverman and Nagler, Tong and Rahman et al. [11-13]. Silverman and Nagler [11] investigated heat transfer between a horizontal water jet and a vertical surface. Tong [12] investigated the effect of liquid inclination angle on the hydrodynamics and heat transfer of the impingement of a liquid jet on a horizontal surface. The locations of the maximum Nusselt number as well as maximum pressure on the surface were found to be identical with the geometric jet impingement point. Rahman et al. [13] solved free surface flow in the presence or absence of gravity. The distribution of film height was found to be strongly affected by the magnitude and orientation of gravity.

Liquid properties effect on free surface flow was investigated by Sun et al. and Liu et al. [14,15]. Sun et al. [14] carried out local measurements to investigate the characteristics of heat transfer from small heater to liquid jet of Prandtl number between 7 and 262. The Nusselt number dependence of Pr1/3 was testified by their experimental data, as well as the data of water and heavy electrochemical liquid from other resources. Bula et al. [16] investigated effect of high Prandtl number fluid for jet impinging perpendicularly on a solid substrate of finite thickness containing small discrete heat sources. It was found that the local heat transfer coefficient had maximum value at the center of the disk and decreases gradually with radius as the flow moves downstream. Liu et al. [15] focused on heat transfer at the stagnation point considering in their analysis the effect of surface tension.

Carper et al. [17] investigated experimentally the heat transfer for impinging liquid jet on a rotating disk. They used petroleum oil as the working medium. The heat transfer coefficients increase approximately with the square root of the rotational velocity of the disk. For two phase flow applications, Robidou et al. [18] investigated the boiling on a hot plate cooled by a water jet. In the forced convection regime, they found that heat fluxes increased with the increase of the subcooling, jet velocity and decrease of the distance from the stagnation line. Boiling first starts in the parallel flow region and propagates in the direction of the jet. In the fully developed nucleate boiling regime, no influence of jet velocity, subcooling and the nozzle-to-plate spacing on the heat flux in the parallel flow region were observed.

Liquid jet and spray impingement cooling were studied experimentally by Oliphant et al. [19]. The comparison of the two cooling techniques revealed that spray cooling can provide the same heat transfer coefficient as jets at a substantially lower mass flux. It was concluded that the effective cooling of non-boiling sprays was primarily due to the unsteady boundary layer resulting from droplet impact and secondarily from evaporative cooling. For nanofluids, considerable works are concentrated on the thermal property measurements and modeling including the study of Williams et al. [20]. Experimental investigations have revealed that nanofluids have remarkably higher thermal conductivities.

Comprehensive reviews of nanofluid heat transfer have been presented by Trisaksri and Somchai [21] and Wang and Mujumdar [22]. Trisaksri and Somchai [21] have concluded that, among a lot of models for thermal conductivity it is not clear yet which model can be considered as the best. They also added that the suspended nanoparticles remarkably increased the forced convective heat transfer of the base fluid. In addition to that, they also noted that the heat transfer of the nanofluid increased as the volume fraction increased for fixed Reynolds number. Wang and Mujumdar [22] reported that the use of nanofluids in a wide range of applications appears promising, but the development of this field faces several challenges: (1) the lack of agreement between experimental results from different groups, (2) the often poor performance of suspensions (3) and lack of the theoretical understanding of the mechanisms. Further theoretical and experimental research investigations were needed to understand the heat transfer characteristics of nanofluids.

For nanofluid flow inside circular tube, experimental and numerical results have indicated that increasing the nanoparticle concentration enhanced the heat transfer coefficient on the cost of the pressure drop [20,23] and [24]. It has been found by Torii [23] that significant enhancement of heat transfer performance was observed in comparison with pure water. This enhancement was intensified with the increase in the Reynolds number and the nanoparticle concentration.

Magïa et al. [25,26], Nguyen et al. [27-29] and Gherasim et al. [30] carried out experimental investigations to study the heat transfer performance of nanofluid (Al2O3) for confined and submerged impinging jets. Experimental data, obtained for turbulent flow regime, have clearly shown that the inclusion of nanoparticles into distilled water has produced a considerable enhancement of the convective heat transfer coefficient. For a particular nanofluid with 6.8% particle volume concentration, the heat transfer coefficient has been found to increase as much as 40% compared to that of the base fluid. Data of Nguyen et al. [28] indicate that the use of a nanofluid can provide a heat transfer enhancement. It has been observed that the highest surface heat transfer coefficients can be achieved using an intermediate nozzle-to-surface distance of 5 mm and a 2.8% nanoparticle volume fraction. Nanofluids with 6% or higher particle volume fraction have been found not appropriate for the heat transfer enhancement purpose under the confined impinging jet configuration.

Manca et al. [31] and Feng and Kleinstreuer [32] numerically investigated the coffined nanofluid jet heat transfer. The numerical data of Manca et al. [31] indicate that Nusselt number increases for increasing particle concentrations and Reynolds numbers. A maximum increase of 18% is detected at a concentration of 6%. The required pumping power as well as Reynolds number increase as the particle concentration grows, which is almost 4.8 times greater than the values calculated for the case of base fluid.

For nanofluid jet boiling, Chakraborty et al. [33] have investigated the cooling of hot steel plate using TiO2 nanofluid jet. Their results of water based-TiO2 nanofluid showed significantly higher cooling rate as compared to the water as a coolant. Chakraborty et al. [33] concluded that convective heat transfer by jet boiling of the nanofluid together with the lower surface tension and higher viscosity of the nanofluid could be important factors leading to this faster cooling.

Nanofluid spray boiling is investigated by Chang et al. [34] and Bellerová et al. [35]. Chang et al. [34] found that the optimal heat transfer performance is obtained using Al2O3 particle volume fraction of 0.001%. In spray cooling with high-volume-fraction nanofluids (0.025 to 0.05%), they found that the nanoparticles were easily deposited on the heated surface, thereby, reducing the number of active nucleation sites and hindering the convection heat transfer mechanism between the surface and the nanofluid. Low-volume-fraction nanofluids (0.001%) yield a significant improvement in the spray cooling efficiency since most of the nanoparticles rebound from the heated surface directly or are washed away by subsequently arriving droplets. The results of Bellerová et al. [35] have found, by comparing the nanofluid results with that of pure water, that an approximately 45% decrease of heat transfer coefficient of spray cooling with the volume fraction of the nanoparticle suspension increasing from 0 to 16.45%. The reduction of heat transfer coefficient caused by the change of the spraying impact duration due to the presence of nanoparticles.

This paper shows the effect of using a vertical alumina-water nanofluid jet with different concentrations on the cooling of horizontal heated circular disks with different diameters in terms of dimensionless parameters. This paper is focused on the effect of disk circular size and its relation to the nanofluid concentration on the cooling process.

Experimental details

The experimental test rig of the current experiment consists of a liquid (water or nanofluid) circuit and a hot circular disk as shown in Figure 1a. Water (or nanofluid) flows from the collecting tank (8) through a heat exchanger (9) using a pump (10). The flow rate is controlled and measured by a valve (11) and a turbine flow meter (12), respectively. The water jet emerges from the nozzle and impinges on a horizontal hot circular disk (1). The falling liquid is collected in a collecting tank (8). A heat exchanger (9) with fan is used to cool the fluid. The circular disk is heated using an electric heater as shown in Figure 1b. The power of the heaters is controlled using a variac transformer (6) (Slideup type made by Voltac Co. Ltd, Tokyo, Japan). The voltage across the heater, the current through the heater and the electric power are measured by a power meter (7) (Metrix PX110; Metrix Instruments, Paris, France). The uncertainities in volt, current and power are ± 0.5, ± 0.7 and ± 1.5 of readings, respectively. The temperatures of the hot disk (1), the jet and the bakalite insulation are measured using k-type thermocouples. The thermocouples (3) are connected to a data acquisition system (4) where the data is collected and stored on a labtop computer (5) for further analyses. The data acquisition type is OMB-DAQ-56 made by Omega Engineering, Inc. (Stamford, CT, USA).

thumbnailFigure 1. Experimental setup. (a) Schematics of test rig and (b) heating plate and (c) thermocouple locations.

Three nozzles of different diameters are used in the current investigation 0.0039, 0.0055 and 0.0082 m. The nozzle length to diameter ratio is kept above 20 to avoid the effect of entrance on jet exit. The height of the nozzle above the hot disk is maintained constant at 50 mm during this investigation. In addition to that, four circular heating plates are used with diameters 0.080, 0.100, 0.115 and 0.133 m. The first two are made of copper and the others of aluminum, respectively. The disks were nickel-electroplated to avoid rusting. A Bakelite insulating housing (3) is built to enclose the circular disk as shown in Figure 1b. The thickness of the hot disks is 25 mm. The electric heaters used are spiral stainless steel heating elements of 0.008-m diameter. The heater is placed in a circular gap of 0.015-m height as shown in Figure 1b. High thermal conductivity cement is used to hold the heater in its cylindrical gap housing. K-type thermocouples made by Thermocoax (Suresnes Cedex, France) are used to measure the disk surface temperatures. The thermocouples are of insulated junction and sheathed by stainless steel of 0.0005-m outside diameter. Those thermocouples are located along the radial direction of the disk each 0.005 m apart as shown in Figure 1c. The thermocouples are impeded in 0.002-m diameter holes in the disk using Omega high thermal conductivity cement (Omega Engineering, Inc., Stamford, CT, USA). The holes are 0.001 m below the disk surface. Two more thermocouples are used to measure the temperature difference across the Bakelite insulation. Water jet temperature is measured using a thermocouple before jet exit.

The thermocupoles and data acquisition system are calibrated at ice and boiling points. At the ice point, the error in thermocople readings are in the range from -0.2 to 0.1°C. At the boiling point, the error in the thermocople readings are in the range from 0.15 to 0.4°C. The flow meter used is a turbine flow meter (FTB602) made by Omega Engineering, Inc. (Stamford, CT, USA). The flow meter is calibrated using a balance scale and a stop watch. The calibration data for water and nanofluid are shown in Figure 2. The difference between the reading of the flow meter and that of the measured flow rates is -2% at most.

thumbnailFigure 2. Flow meter calibration.

Nanofluid preparation

Alumina nanofluid (10-nm particles of Al2O3 dispersed in water) of 20% mass concentration supplied by the Nanostructured & Amorphous Materials, Inc. (Houston, TX, USA) is used in the present investigation. The manufacturer information about this nanofluid is that the average nanoparticle size = 10 nm gamma and it looks as a transparent liquid. The manufacturer 20%-mass concentration nanofluid was diluted to different concentrations using distilled water. These new diluted solutions are ultrasonic vibrated for about 6 h to insure complete dispersions of the nanoparticles. Still camera is used to check any precipitation. Photographs are taken just after ultrasonic vibration and 48 h later. No precipitation is observed during this period. Two mass concentrations, 6.6% and 10%, are used in the current investigation.

Nanofluid density is measured using DMA 35N density meter. Figure 3 shows the effect of mass concentration on the measured nanofluid density at ambient temperature of 25.5°C. The nanofluid viscosity at ambient temperature 24.5°C is measured using Viscolite 700. The Effect of nanoparticle concentration on viscosity is shown in Figure 4. As seen in this figure, the viscosity value is doubled as the concentration increases from 0% to 10%. The nanofluid thermal conductivity is measured using KD2Pro. The Effect of of nanoparticle concentration on thermal conductivity is shown in Figure 5. As shown in this figure, the thermal conductivity increased about 15% as the concentration increases from 0% to 10%.

thumbnailFigure 3. Nanofluid density at ambient temperature of 25.5°C.

thumbnailFigure 4. Nanofluid viscosity at 24.5°C.

thumbnailFigure 5. Nanofluid thermal conductivity at 26.5°C.

Comparison with the viscosity and the thermal conductivity models available in the literature revealed that the models proposed by Williams et al. [20], represented by Equations 3 to 6, can reasonably present the current measurements:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M3">View MathML</a>

(3)

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M4">View MathML</a>

(4)

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M5">View MathML</a>

(5)

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M6">View MathML</a>

(6)

where φ is the nanoparticle volumetric concentration, ρ, C, k and μ are the nanofluid density, specific heat, thermal conductivity and viscosity, respectively. For the based fluid, ρb, Cb, kb and μb are the density, specific heat, thermal conductivity and viscosity, respectively. For the nanoparticles, ρp, Cp and kp are the density, specific heat and thermal conductivity, respectively. The relation between mass concentration × and the volumetric concentration φ can be estimated from the following equation and Equation 3:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M7">View MathML</a>

(7)

where X is the nanofluid particle mass concentration; ρp, fluid density of a nanoparticle; ρb, fluid density base. The measurred viscosity falls within ± 20% of Equation 6, and the measurred thermal conductivity falls within ± 5% of Equation 5.

Experimental procedure

Experiments are carried out to examine the effect of nanofluid flow rates and concentrations on steady state cooling of horizontal circular disks. Results are obtained for mass concentrations of 0.0%, 6.6% and 10.0%, mass flow rate in the range 0.006 to 0.075 kg/s and for heat flux in the range 60 to 100 kW/m2. Water is used as liquid reference (0 concentrations) in the current investigation since it is the base fluid of the used nanofluid. Three nozzles are used to examine jet size effect on heat transfer from hot disks. Four disks are used to examine effect of disk size on heat transfer.

At the beginning of the experiment, the control valve is used to establish the required flow rate through the nozzle. Then, the heater is turned on where the electric power is adjusted using the variac transformer and recorded. The data acquisition starts to collect thermocouple readings and saves them in the computer. It should be noted that the system reaches steady state in 30 min where the temperatures are recorded through the data acquisition system. The experiments are done first for pure water only, then for nanofluid with different concentrations.

Data analysis

The average temperature Twam of the hot surface is calculated as:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M8">View MathML</a>

(8)

where Twk is the local wall temperature for the element area Ak. Twam is the average temperature of the local wall. The heat loss Qloss through the Bakelite housing is estimated as:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M9">View MathML</a>

(9)

where kins is the insulation thermal conductivity, Ains is the insulation disk area, ΔTins is the temperature difference across the insulation and δins is the insulation thickness. The heat transfer rate Q between the jet and the hot surface is calculated from the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M10">View MathML</a>

(10)

where Q is heat transfer and Qele, electrical heat transfer. The average temperature Twam of the hot surface is obtained by correcting Twam due to the effect of metal thickness, δth = 0.001 m between the hot surface and the thermocouple locations as seen in Figure 1c:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M11">View MathML</a>

(11)

where km is the local temperature of metal and q is heat flux. The average heat transfer coefficient based on temperature difference between the wall average temperature Twa and the water jet temperature Tj (before striking the hot surface) can be calculated as:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M12">View MathML</a>

(12)

where h is the average heat transfer coefficient. The Nusselt number, NuD based on disk diameter is given by:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M13">View MathML</a>

(13)

where hD is the average heat transfer coefficient of the disk diameter. The jet velocity at the nozzle exit, Vj, is obtained using the continuity equation at the nozzle exit. The jet velocity at impingement point can be obtained by applying the energy equation between the nozzle exit and just before the impingement point as:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M14">View MathML</a>

(14)

where Zo is the nozzle height above the horizontal disk plate, Vi is the jet velocity at impingement and Vj is the velocity at nozzle exit. It was maintained at 50 mm in this study. Furthermore, the jet diameter Di just before the impingement is calculated from the continuity equation just before impingement:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M15">View MathML</a>

(15)

where Dj is the nozzle diameter. The jet Reynolds number Rei based on the jet velocity and the diameter at the impingement point is given by:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M16">View MathML</a>

(16)

It should be noted that the physical properties of water used in the above equations are estimated at the film temperature Tf:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M17">View MathML</a>

(17)

where Tj is the temperature of the jet.

Results and discussion

Uncertainty propagation technique is used to calculate how the uncertainties in each of the measured variables propagate into the value of the calculated quantities. The method for determining this uncertainty propagation is described in Taylor and Kuyatt [36]. The uncertainties in the measured quantities are considered as follows: in temperature, ± 0.4°C; in length and diameter, ± 0.05 mm; in nanofluid mass, ± 0.1 g; in electric heating power, ± 1.5% and in flow rate, ± 2%. The estimated uncertainity in the calculated Peclet and Nusselt numbers falls within ± 2% and ± 8.7%, respectively.

Samples of the measured wall temperature along the radial direction are shown in Figure 6. As seen in this figure, the lowest wall temperature is observed at the impingement point (radial distance = 0). For the same Reynolds number, the wall temperature increases along the radial direction. This indicates that the heat transfer coefficient decreases along the radial direction. The reduction in the heat transfer is caused by the increase in the liquid film thickness and the decrease in the film velocity along the radial direction. As the jet Reynolds number increases, the wall temperature decreases as seen in Figure 6. These results reveal that the jet Reynolds number plays a major role in the heat transfer process.

thumbnailFigure 6. Measured local wall temperature along the radial direction. (a) Data for heating disk of diameter D = 0.08 m, (b) 0.10 m, (c) 0.115m and (d) 0.133 m.

Comparison between the current data for water jet and the model of Zhao et al. [8] (Equation 1) and Teamah and Farahat [10] (Equation 2) are shown in Figures 7 and 8, respectively. The current data fall within ± 20% of Zhao et al.'s model [8]; however, Teamah and Farahat's model [10] under-predicts the current data by about 25%. Nevertheless, these comparisons indicate that the current data fall within the predicted models.

thumbnailFigure 7. Comparison with the model of Zhao and Masuoka [8]. Black triangles represent current data of water jets.

thumbnailFigure 8. Comparison with the model of Teamah and Farahat [10].

Preliminary investigation has shown that the best presentation of the data is to be in terms of the Nusselt number NuD based on disk diameter versus the jet Reynolds number Rei based on jet velocity and diameter at the impingement point. It should be noted that using the Vi at the impingement point in estimating the Reynolds number takes into account the effect of the nozzle height. Figures 9,10,11 show the data of the mean Nusselt numbers versus the Reynolds number Rei for the nanoparticle concentrations of 0%, 6.6% and 10%, respectively, for different jets Dj = 0.0039, 0.0055 and 0.0082 m and heated disks of diameters D = 0.080, 0.100, 0.115 and 0.133 m. The general trend drawn from these figures is that as the Reynolds number increases the Nusselt number increases too. It should be noted that using Rei in the data presentation eliminates the effect of jet diameter and height Zo on NuD profiles as shown in this figures. Figure 9a,b,c,d show that all the data for pure water are heating disk-independent and, therefore, collapse on one curve. As the nanoparticle concentration percentage increases, the absolute viscosity increases which leads to reduced Reynolds number as seen in Figure 10a,b,c,d and Figure 11a,b,c,d for 6.6% and 10% concentrations, respectively. In spite of using different jet diameters, the experimental data collapse on one curve for each heater disk. It was also found that presenting the data in terms of Reynolds number at the impingement disk can take into account the effect of nozzle height and jet diameter. In other words, using this choice of Reynolds number makes NuD less sensitive to nozzle diameter change.

thumbnailFigure 9. Nusselt numbers versus Reynolds numbers for water. (a) Data for heating disk of diameter D = 0.08 m, (b) 0.10 m, (c) 0.115m and (d) 0.133 m.

thumbnailFigure 10. Nusselt numbers versus Reynolds numbers for 6.6% nanofluid concentration. (a) Data for heating disk of diameter D = 0.08 m, (b) 0.10 m, (c) 0.115m and (d) 0.133 m.

thumbnailFigure 11. Nusselt numbers versus Reynolds numbers for 10% nanofluid concentration. (a) Data for heating disk of diameter D = 0.08 m, (b) 0.10 m, (c) 0.115m and (d) 0.133 m.

Cross plots of Figures 9 to 11 for fixed disk diameter D = 0.100 m and for different concentrations are shown in Figure 12a,b,c. Figure 12a shows the cross plot for Dj = 0.0039 m, (b) for Dj = 0.0055 m and (c) for 0.0082 m. Figure 12a,b,c shows increase in NuD as the concentration increases. For example, at Rei = 8,000 the Nusselt number NuD increases from 500 at zero concentration to 700 at 6.6% concentration up to 1,000 at 10% concentration. For the same Reynolds number, the experimental data show an increase in Nusselt number that can reach up to 100% at higher concentration. This result indicates that using nanofluid as a heat transfer carrier can enhance the heat transfer process. It should be noted that, as the concentration increases, Reynolds number decreases as mentioned earlier, and the data for 10% concentration is shifted to the left and that for pure water is on the right. It means that the increase in Nusselt number is on the account of force required to push the flow. Figure 13 is constructed to show the effect of nanofluid concentration on the Nusselt number. The data in this figure is extracted from Figure 12. As seen in Figure 13, the increase in Nusselt number can reach up to 40% and 75% at 6.6% and 10% concentrations, respectively.

thumbnailFigure 12. Comparison between Nusselt numbers for water and nanofluid for 100-mm diameter heater. (a) Cross plots for Dj = 0.0039 m, (b) for Dj = 0.0055 m and (c) for Dj = 0.0082 m.

thumbnailFigure 13. Enhancement in Nusselt number due to an increase in the nanofluid concentration. (a) Data for heating disk of diameter D = 0.08 m and (b) 0.10 m.

The effect of disk diameter and noanofluid concentration is shown in Figure 14. In this figure, the heat transfer coefficient is presented in terms of NuD /D or (h /k) versus disk diameter D. This ratio was chosen such that Nusselt number be independent of disk diameter D. These data indicate that heat transfer coefficient decreases as heating disk diameter increases. This decrease could be attributed to the fact that, as the liquid film thickness increases along the radial direction, the film velocity decreases. Therefore, as the heating disk diameter increases, the local heat transfer coefficient along the radial direction of the disk decreases. Consequently, the average heat transfer coefficient decreases. Furthermore, the opportunity of the flow to turn from shooting before the hydraulic jump to a streaming flow after the hydraulic jump increases as the disk diameter increases. The sudden increase in the film thickness after the hydraulic jump leads to a significant decrease in the heat transfer coefficient. The effect of nanofluid concentration on the heat transfer is also shown in Figure 14. This figure shows that increasing the concentration enhances the heat transfer coefficient. Using 10% concentration can increase the heat transfer (NuD /D or h /k) about 60%. Moreover, the actual increase in the heat transfer coefficient should be more due to the increase in thermal conductivity k of the nanofluid.

thumbnailFigure 14. The effect of disk diameter and nanofluid concentration on heat transfer.

Figure 15a,b,c presents the profiles of Nusselt number versus Peclet number, Pei, for disk diameter D = 0.100 m and for different concentrations. Figure 15a shows the data for jet diameter Dj = 0.0039 m, (b) for 0.0055 m and (c) for 0.0082 m, respectively. As shown in this figure, the data of various concentrations are collapsed on one curve when Pei is used in presenting the data. The nanoparticle concentration has a strong impact on viscosity, and consequently, on Prandtl number. Therefore, as the nanoparticle concentration increases from 05 to 10%; the viscosity increases about 150%. However, using Peclet number in Figure 15a,b,c ensures the effect of Prandtl number which makes the data to collapse on one curve.

thumbnailFigure 15. Nusselt numbers versus Peclet numbers for 100 mm diameter heater. (a) Data for jet diameter Dj = 0.0039 m, (b) 0.0055 m and (c) 0.0082 m.

Different correlations were examined to fit the data in Figure 15a,b,c. It is found that the best correlation trend is obtained by correlating the experimental data in terms of Pei, nanofluid concentration represented by (1-X) where X is the percentage of nanoparticle concentration, and heating disk to jet diameter ratio D /Di. The effect of diameter ratio is chosen to be similar to that of Zhao et al.'s model [8] in Equation 1. The correlation obtained is:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/139/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/139/mathml/M18">View MathML</a>

(18)

where 0 ≤ × ≤ 0.1 with a correlation coefficient R = 96.5%. Comparison between the predictions of this correlation and the experimental data are shown in Figure 16. This figure shows that 96% of the data fall within error bands of ± 20%.

thumbnailFigure 16. Comparison between the predicted and the experimental Nusselt numbers. Black squares represent predicted Nusselt number using Equation 18 at experimental conditions of Nusselt number shown on x-axis.

Conclusions

Experimental investigation to study the heat transfer between a vertical round water jet having nanoparticles of aluminum oxide, and a horizontal circular round surface is carried out for different jet flow rates, jet diameters, nanoparticle concentrations and heating surface diameters. Three different nanoparticle concentrations 0%, 6.6% and 10% are used in the current investigation. For the same Reynolds number, the experimental data show an increase in the Nusselt number that can reach up to 100% for some higher concentrations. This result indicates that using nanofluid as a heat transfer carrier can enhance the heat transfer process. It was also found that presenting the data in terms of Reynolds number at impingement jet diameter can take into account the effect of jet heights and nozzle diameters. The data have also indicated that increasing heating disk diameter decreases the heat transfer coefficient. Experimental data was correlated in terms of jet Peclet number Pei, nanofluid concentration presented by (1-X) and diameter ratio (D /Di) with a correlation coefficient of R = 96.5%. Predictions of the obtained correlation fall within error bands of ± 20% of the experimental data.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

OZ designed and installed experimental rig, carried out experimental runs, analyzed experimental data and wrote and prepared manuscript. MA participated in installing experimental rig, data analysis, preparing and writing the manuscript. All authors read and approved the final manuscript.

Appendix

Nomenclature

A Area of disk, πD2/4, square meter

Ak Area element assigned to each thermocouple, square meter

Aj Area of jet nozzle, πDj2/4, square meter

C Specific heat, joules per kilogram Kelvin

D Disk diameter, meter

Di Jet diameter at impingement, meter

Dj Nozzle diameter, meter

h Average heat transfer coefficient, watts per square meter Kelvin

k Thermal conductivity, watts per meter Kelvin

NuD Nusselt number, hD /k

Nuj Nusselt number, hDj /k

Pei Peclet number, Rei Pr

Rei Jet Reynolds number at impingement, ρVi Di /μ

Rej Nozzle Reynolds number, ρVj Dj

T Temperature, Kelvin

Vi Jet velocity at impingement, meter per second

Vj Velocity at nozzle exit, meter per second

Q Heat transfer, watt

q Heat flux, Q /A, watt per square meter

X Nanofluid particle mass concentration, percent

Zo Height of jet above surface disk, meter

Greek symbols

ν kinematic viscosity, square meter per second

μ viscosity, Pascal second

ρ fluid density, kilogram per cubic meter

φ volumetric concentration, percent

Subscripts

a average

b base fluid

ele electrical

i impingement

ins insulation

j jet

k local temperature

loss losses

m metal

p nanoparticle

w wall

Acknowledgements

This experimental investigation is supported by King Abdullah Institute for Nanotechnology at King Saud University under the project no. 42/1429. This support is highly appreciated and acknowledged.

References

  1. Liu X, Lienhard J: Extremely high heat fluxes beneath impinging jets.

    J Heat Transfer 1993, 115:472-476. Publisher Full Text OpenURL

  2. Lienhard J, Hadeler J: High heat flux cooling by liquid jet-array modules.

    Chem Eng Technology 1999, 22(11):967-970. Publisher Full Text OpenURL

  3. Bergles A: High-flux processes through enhanced heat transfer. Rohsenow Symposium-MIT, Cambridge, Mass; 2003.

  4. Jambunathan K, Lai E, Moss M, Button B: A review of heat transfer data for single circular jet impingement.

    Int J Heat and Fluid Flow 1992, 13(2):106-115. Publisher Full Text OpenURL

  5. Watson E: The radial spread of a liquid jet over a horizontal plans.

    J Fluid Mech 1962, 20:part 3, 481-499. OpenURL

  6. Chaudhury Z: Heat Transfer in a Radial Liquid Jet.

    J Fluid Mech 1964, 20:part 3, 501-511. OpenURL

  7. Ishigai S: Heat transfer of the impinging round water jet in the interference zone of film flow along the wall.

    Bulletin of the JSME 1977, 20(139):85-92. Publisher Full Text OpenURL

  8. Zhao Y, Masuoka T: Conjugated heat transfer on a horizontal surface impinged by circular free-surface liquid jet.

    JSME International Journal Series B 2002, 45(2):307-314. Publisher Full Text OpenURL

  9. Baonga J, Louahlia-Gualous H, Imbert M: Experimental study of the hydrodynamic and heat transfer of free liquid jet impinging a flat circular heated disk.

    Applied Thermal Engineering 2006, 26:1125-1138. Publisher Full Text OpenURL

  10. Teamah M, Farahat S: Experimental and numerical heat transfer from impinging of single free liquid jet.

    Alexandria Engineering Journal 2003, 42(5):559-575. OpenURL

  11. Silverman I, Nagler A: High heat flux cooling with water jet impingement. Proceedings of 2004 ASME Heat Transfer/Fluids Engineering Summer Conference, Charlotte, North Carolina: 2004, July 11-15;

  12. Tong A: On the impingement heat transfer of an oblique free surface plane jet.

    International Journal of Heat and Mass Transfer 2003, 46:2077-2085. Publisher Full Text OpenURL

  13. Rahman M, Faghri A, Hankey W, Swanson T: Computation of the free surface flow of a thin liquid film at zero and normal gravity.

    Numerical Heat Transfer Part A 1990, 17:53-71. Publisher Full Text OpenURL

  14. Sun H, Ma C, Chen Y: Prandtl number dependence of impingement heat transfer with circular free-surface liquid jets.

    lnt J Heat Mass Transfer 1998, 41(10):1360-1363. Publisher Full Text OpenURL

  15. Liu X, Gabour L, Lienhard J: Stagnation-point heat transfer during impingement of laminar liquid jets: analysis including surface tension.

    J Heat Transfer 1993, 115:99-105. Publisher Full Text OpenURL

  16. Bula A, Rahman M, Leland J: Axial steady free surface jet impinging over a flat disk with discrete heat sources.

    International Journal of Heat and Fluid Flow 2000, 21:11-21. Publisher Full Text OpenURL

  17. Carper H, Savedra J, Suwanprateep T: Liquid jet impingement cooling of a rotating disc.

    Journal of Heat Transfer, Transactions of the ASME 1986, 108:540-544. Publisher Full Text OpenURL

  18. Robidou H, Auracher H, Gardin P, Lebouch M: Controlled cooling of a hot plate with a water jet.

    Experimental Thermal and Fluid Science 2002, 26:123-129. Publisher Full Text OpenURL

  19. Oliphant K, Webb B, McQuay M: An experimental comparison of liquid jet array and spray impingement cooling in the non-boiling regime.

    Experimental Thermal and Fluid Science 1998, 18:1-10. Publisher Full Text OpenURL

  20. Williams W, Buongiorno J, Hu L: Experimental investigation of turbulent convective heat transfer and pressure loss of alumina/water and zirconia/water nanoparticle colloids nanofluids in horizontal tubes.

    Journal of Heat Transfer 2008., 130(4)

    042412

    OpenURL

  21. Trisaksri V, Wongwises S: Critical review of heat transfer characteristics of nanofluids. Renewable and Sustainable Energy Reviews.

    2007, 11(3):512-523.

  22. Wang X, Mujumdar A: A review on nanofluids - part ii: experiments and applications.

    Brazilian Journal of Chemical Engineering 2008, 25(4):631-648. Publisher Full Text OpenURL

  23. Torii S: Experimental study on convective heat transfer of aqueous suspensions of nano-diamond particles.

    Proceedings of International Symposium on EcoTopia Science 2007 ISETS07 2007 OpenURL

  24. Ali M, Zeitoun O: Nanofluids forced convection heat transfer inside circular tubes.

    Int J Nanoparticles 2009, 2(1-6):164-171. OpenURL

  25. Maïga S, Palm S, Nguyen C, Roy G, Galanis N: Heat transfer enhancement by using nanofluids in forced convection flows.

    International Journal of Heat and Fluid Flow 2005, 26:530-546. Publisher Full Text OpenURL

  26. Maïga S, Nguyen C, Galanis N, Roy G: Heat transfer behaviours of nanofluids in a uniformly heated tube.

    Superlattices and Microstructures 2004, 35:543-557. Publisher Full Text OpenURL

  27. Nguyen C, Roy G, Gauthier C, Galanis N: Heat transfer enhancement using Al2O3-water nanofluid for an electronic liquid cooling system.

    Applied Thermal Engineering 2007, 27:1501-1506. Publisher Full Text OpenURL

  28. Nguyen C, Galanis N, Polidori G, Fohanno S, Popa C, Bechec A: An experimental study of a confined and submerged impinging jet heat transfer using al2o3-water nanofluid.

    International Journal of Thermal Sciences 2009, 48:401-411. Publisher Full Text OpenURL

  29. Nguyen C, Laplante1 G, Cury M, Simon G: Experimental investigation of impinging jet heat transfer and erosion effect using al2o3-water nanofluid. 6th IASME/WSEAS International Conference on Fluid Mechanics and Aerodynamics (FMA'08), Rhodes, Greece, August 20-22, 2008

  30. Gherasim I, Roy G, Nguyen CT, Vo-Ngoc D: Heat transfer enhancement and pumping power in confined radial flows using nanoparticle suspensions (nanofluids).

    International Journal of Thermal Sciences 2011, 50:369-377. Publisher Full Text OpenURL

  31. Manca O, Mesolella P, Nardini S, Ricci D: Numerical study of a confined slot impinging jet with nanofluids.

    Nanoscale Research Letters 2011, 6:188. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  32. Feng Y, Kleinstreuer C: Nanofluid convective heat transfer in a parallel-disk system.

    International Journal of Heat and Mass Transfer 2010, 53:4619-4628. Publisher Full Text OpenURL

  33. Chakraborty S, Chakraborty A, Das S, Mukherjee T, Bhattacharjee D, Rayisij R: Application of water based-tio2 nano-fluid for cooling of hot steel plate.

    ISIJ International 2010, 50(1):124-127. Publisher Full Text OpenURL

  34. Chang T, Syu S, Yang Y: Effects of particle volume fraction on spray heat transfer performance of Al2O3-water nanofluid.

    International Journal of Heat and Mass Transfer 2012, 55:10141021. OpenURL

  35. Bellerová H, Tseng AA, Pohanka M, Raudensky M: Spray cooling by solid jet nozzles using alumina/water nanofluids.

    International Journal of Thermal Sciences 2011, in press.

    available online

    OpenURL

  36. Taylor BN, Kuyatt CE: Guidelines for evaluating and expressing the uncertainty of nist measurement results. NIST Technical Note 1297, National Institute of Standards and Technology; 1994.