Abstract
Experimental investigation to study the heat transfer between a vertical round aluminawater 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 jetIntroduction
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 [710]. 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:
where Nu_{j }is the Nusselt number based on the jet nozzle diameter D_{j}, D is the disk diameter; Pr is the Prandtl number and Re_{j }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 twodimensional model in radial and vertical directions. They solved the governing equations numerically up to r /D_{j }= 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 Re_{j }range from 5,000 to 20,000:
Different jet orientations were also investigated by Silverman and Nagler, Tong and Rahman et al. [1113]. 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 Pr^{1/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 nozzletoplate 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 nonboiling 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. [2729] and Gherasim et al. [30] carried out experimental investigations to study the heat transfer performance of nanofluid (Al_{2}O_{3}) 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 nozzletosurface 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 TiO_{2 }nanofluid jet. Their results of water basedTiO_{2 }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 Al_{2}O_{3 }particle volume fraction of 0.001%. In spray cooling with highvolumefraction 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. Lowvolumefraction 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 aluminawater 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 ktype 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 OMBDAQ56 made by Omega Engineering, Inc. (Stamford, CT, USA).
Figure 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 nickelelectroplated 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.008m diameter. The heater is placed in a circular gap of 0.015m height as shown in Figure 1b. High thermal conductivity cement is used to hold the heater in its cylindrical gap housing. Ktype 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.0005m 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.002m 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.
Figure 2. Flow meter calibration.
Nanofluid preparation
Alumina nanofluid (10nm particles of Al_{2}O_{3 }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%.
Figure 3. Nanofluid density at ambient temperature of 25.5°C.
Figure 4. Nanofluid viscosity at 24.5°C.
Figure 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:
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}, C_{b}, k_{b }and μ_{b }are the density, specific heat, thermal conductivity and viscosity, respectively. For the nanoparticles, ρ_{p}, C_{p }and k_{p }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:
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/m^{2}. 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 T_{wam }of the hot surface is calculated as:
where T_{wk }is the local wall temperature for the element area A_{k}. T_{wam }is the average temperature of the local wall. The heat loss Q_{loss }through the Bakelite housing is estimated as:
where k_{ins }is the insulation thermal conductivity, A_{ins }is the insulation disk area, ΔT_{ins }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:
where Q is heat transfer and Q_{ele}, electrical heat transfer. The average temperature T_{wam }of the hot surface is obtained by correcting T_{wam }due to the effect of metal thickness, δ_{th }= 0.001 m between the hot surface and the thermocouple locations as seen in Figure 1c:
where k_{m }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 T_{wa }and the water jet temperature T_{j }(before striking the hot surface) can be calculated as:
where h is the average heat transfer coefficient. The Nusselt number, Nu_{D }based on disk diameter is given by:
where hD is the average heat transfer coefficient of the disk diameter. The jet velocity at the nozzle exit, V_{j}, 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:
where Z_{o }is the nozzle height above the horizontal disk plate, V_{i }is the jet velocity at impingement and V_{j }is the velocity at nozzle exit. It was maintained at 50 mm in this study. Furthermore, the jet diameter D_{i }just before the impingement is calculated from the continuity equation just before impingement:
where D_{j }is the nozzle diameter. The jet Reynolds number Re_{i }based on the jet velocity and the diameter at the impingement point is given by:
It should be noted that the physical properties of water used in the above equations are estimated at the film temperature T_{f}:
where T_{j }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.
Figure 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] underpredicts the current data by about 25%. Nevertheless, these comparisons indicate that the current data fall within the predicted models.
Preliminary investigation has shown that the best presentation of the data is to be in terms of the Nusselt number Nu_{D }based on disk diameter versus the jet Reynolds number Re_{i }based on jet velocity and diameter at the impingement point. It should be noted that using the V_{i }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 Re_{i }for the nanoparticle concentrations of 0%, 6.6% and 10%, respectively, for different jets D_{j }= 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 Re_{i }in the data presentation eliminates the effect of jet diameter and height Z_{o }on Nu_{D }profiles as shown in this figures. Figure 9a,b,c,d show that all the data for pure water are heating diskindependent 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 Nu_{D }less sensitive to nozzle diameter change.
Figure 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.
Figure 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.
Figure 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 D_{j }= 0.0039 m, (b) for D_{j }= 0.0055 m and (c) for 0.0082 m. Figure 12a,b,c shows increase in Nu_{D }as the concentration increases. For example, at Re_{i }= 8,000 the Nusselt number Nu_{D }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.
Figure 12. Comparison between Nusselt numbers for water and nanofluid for 100mm diameter heater. (a) Cross plots for D_{j }= 0.0039 m, (b) for D_{j }= 0.0055 m and (c) for D_{j }= 0.0082 m.
Figure 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 Nu_{D }/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 (Nu_{D }/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.
Figure 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, Pe_{i}, for disk diameter D = 0.100 m and for different concentrations. Figure 15a shows the data for jet diameter D_{j }= 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 Pe_{i }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.
Figure 15. Nusselt numbers versus Peclet numbers for 100 mm diameter heater. (a) Data for jet diameter D_{j }= 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 Pe_{i}, nanofluid concentration represented by (1X) where X is the percentage of nanoparticle concentration, and heating disk to jet diameter ratio D /D_{i}. 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:
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%.
Figure 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 xaxis.
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 Pe_{i}, nanofluid concentration presented by (1X) and diameter ratio (D /D_{i}) 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, πD^{2}/4, square meter
A_{k }Area element assigned to each thermocouple, square meter
A_{j }Area of jet nozzle, πD_{j}^{2}/4, square meter
C Specific heat, joules per kilogram Kelvin
D Disk diameter, meter
D_{i }Jet diameter at impingement, meter
D_{j }Nozzle diameter, meter
h Average heat transfer coefficient, watts per square meter Kelvin
k Thermal conductivity, watts per meter Kelvin
Nu_{D }Nusselt number, hD /k
Nu_{j }Nusselt number, hD_{j }/k
Pe_{i }Peclet number, Re_{i }Pr
Re_{i }Jet Reynolds number at impingement, ρV_{i }D_{i }/μ
Re_{j }Nozzle Reynolds number, ρV_{j }D_{j }/μ
T Temperature, Kelvin
V_{i }Jet velocity at impingement, meter per second
V_{j }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
Z_{o }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

Liu X, Lienhard J: Extremely high heat fluxes beneath impinging jets.
J Heat Transfer 1993, 115:472476. Publisher Full Text

Lienhard J, Hadeler J: High heat flux cooling by liquid jetarray modules.
Chem Eng Technology 1999, 22(11):967970. Publisher Full Text

Bergles A: Highflux processes through enhanced heat transfer. Rohsenow SymposiumMIT, Cambridge, Mass; 2003.

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):106115. Publisher Full Text

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

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):8592. Publisher Full Text

Zhao Y, Masuoka T: Conjugated heat transfer on a horizontal surface impinged by circular freesurface liquid jet.
JSME International Journal Series B 2002, 45(2):307314. Publisher Full Text

Baonga J, LouahliaGualous 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:11251138. Publisher Full Text

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

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 1115;

Tong A: On the impingement heat transfer of an oblique free surface plane jet.
International Journal of Heat and Mass Transfer 2003, 46:20772085. Publisher Full Text

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:5371. Publisher Full Text

Sun H, Ma C, Chen Y: Prandtl number dependence of impingement heat transfer with circular freesurface liquid jets.
lnt J Heat Mass Transfer 1998, 41(10):13601363. Publisher Full Text

Liu X, Gabour L, Lienhard J: Stagnationpoint heat transfer during impingement of laminar liquid jets: analysis including surface tension.
J Heat Transfer 1993, 115:99105. Publisher Full Text

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:1121. Publisher Full Text

Carper H, Savedra J, Suwanprateep T: Liquid jet impingement cooling of a rotating disc.
Journal of Heat Transfer, Transactions of the ASME 1986, 108:540544. Publisher Full Text

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:123129. Publisher Full Text

Oliphant K, Webb B, McQuay M: An experimental comparison of liquid jet array and spray impingement cooling in the nonboiling regime.
Experimental Thermal and Fluid Science 1998, 18:110. Publisher Full Text

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

Trisaksri V, Wongwises S: Critical review of heat transfer characteristics of nanofluids. Renewable and Sustainable Energy Reviews.
2007, 11(3):512523.

Wang X, Mujumdar A: A review on nanofluids  part ii: experiments and applications.
Brazilian Journal of Chemical Engineering 2008, 25(4):631648. Publisher Full Text

Torii S: Experimental study on convective heat transfer of aqueous suspensions of nanodiamond particles.
Proceedings of International Symposium on EcoTopia Science 2007 ISETS07 2007

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

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:530546. Publisher Full Text

Maïga S, Nguyen C, Galanis N, Roy G: Heat transfer behaviours of nanofluids in a uniformly heated tube.
Superlattices and Microstructures 2004, 35:543557. Publisher Full Text

Nguyen C, Roy G, Gauthier C, Galanis N: Heat transfer enhancement using Al2O3water nanofluid for an electronic liquid cooling system.
Applied Thermal Engineering 2007, 27:15011506. Publisher Full Text

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 al2o3water nanofluid.
International Journal of Thermal Sciences 2009, 48:401411. Publisher Full Text

Nguyen C, Laplante1 G, Cury M, Simon G: Experimental investigation of impinging jet heat transfer and erosion effect using al2o3water nanofluid. 6^{th }IASME/WSEAS International Conference on Fluid Mechanics and Aerodynamics (FMA'08), Rhodes, Greece, August 2022, 2008

Gherasim I, Roy G, Nguyen CT, VoNgoc D: Heat transfer enhancement and pumping power in confined radial flows using nanoparticle suspensions (nanofluids).
International Journal of Thermal Sciences 2011, 50:369377. Publisher Full Text

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

Feng Y, Kleinstreuer C: Nanofluid convective heat transfer in a paralleldisk system.
International Journal of Heat and Mass Transfer 2010, 53:46194628. Publisher Full Text

Chakraborty S, Chakraborty A, Das S, Mukherjee T, Bhattacharjee D, Rayisij R: Application of water basedtio2 nanofluid for cooling of hot steel plate.
ISIJ International 2010, 50(1):124127. Publisher Full Text

Chang T, Syu S, Yang Y: Effects of particle volume fraction on spray heat transfer performance of Al2O3water nanofluid.
International Journal of Heat and Mass Transfer 2012, 55:10141021.

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

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.