This article is part of the series Nanofluids.

Open Access Highly Accessed Nano Express

Numerical study of a confined slot impinging jet with nanofluids

Oronzio Manca*, Paolo Mesolella, Sergio Nardini and Daniele Ricci

Author affiliations

Dipartimento di Ingegneria Aerospaziale e Meccanica, Seconda Università degli Studi di Napoli, Via Roma 29 - 81031 Aversa, Italy

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2011, 6:188  doi:10.1186/1556-276X-6-188


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


Received:14 December 2010
Accepted:1 March 2011
Published:1 March 2011

© 2011 Manca et al; 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

Background

Heat transfer enhancement technology concerns with the aim of developing more efficient systems to satisfy the increasing demands of many applications in the fields of automotive, aerospace, electronic and process industry. A solution for obtaining efficient cooling systems is represented by the use of confined or unconfined impinging jets. Moreover, the possibility of increasing the thermal performances of the working fluids can be taken into account, and the introduction of nanoparticles in a base fluid can be considered.

Results

In this article, a numerical investigation on confined impinging slot jet working with a mixture of water and Al2O3 nanoparticles is described. The flow is turbulent and a constant temperature is applied on the impinging. A single-phase model approach has been adopted. Different geometric ratios, particle volume concentrations and Reynolds number have been considered to study the behavior of the system in terms of average and local Nusselt number, convective heat transfer coefficient and required pumping power profiles, temperature fields and stream function contours.

Conclusions

The dimensionless stream function contours show that the intensity and size of the vortex structures depend on the confining effects, given by H/W ratio, Reynolds number and particle concentrations. Furthermore, for increasing concentrations, nanofluids realize increasing fluid bulk temperature, as a result of the elevated thermal conductivity of mixtures. The local Nusselt number profiles show the highest values at the stagnation point, and the lowest at the end of the heated plate. The average Nusselt number increases for increasing particle concentrations and Reynolds numbers; moreover, the highest values are observed for H/W = 10, and a maximum increase of 18% is detected at a concentration equal to 6%. The required pumping power as well as Reynolds number increases and particle concentrations grow, which is almost 4.8 times greater than the values calculated in the case of base fluid.

List of symbols

Background

Heat transfer enhancement is very important in the industry, and several techniques are employed to realize this aim. Impinging jets, whether confined or unconfined, have been widely used for efficient cooling in industrial applications as a means of providing highly localized heat transfer coefficients, representing a possible solution. Depending on the application, flow conditions can range from laminar to highly turbulent ones. Applications of impinging jets include drying of textiles, film and paper, cooling of gas turbine components and the outer walls of combustors, freezing of tissue in cryosurgery and manufacturing, material processing and electronic cooling. There are numerous articles dealing with this problem both numerically and experimentally as reported in the literature reviews on the subject [1-6].

Several studies have been developed on impinging air jets [1,2]. Recently, a greater attention has been dedicated to the impinging liquid jet since orders of magnitude of the heat transfer rates are several times those of gas jets. Liquid jets have possible application to the cooling of heat engines [5,7], thermal control in electronic devices [8-10] and in the thermal treatment of metals and material processing [11-14].

In the application of jet impingements, circular or slot jets are the main jet configurations. For these two configurations, flow and heat transfer mechanics are significantly different. It seems that greater research activity on heat and mass transfer with circular impinging jets has been predominantly published [1-3,15,16]. However, investigations on heat and mass transfer with slot jet impingement have attracted more attention recently. In fact, slot jet impingements offer many more beneficial features, such as higher cooling effectiveness, greater uniformity and more controllability, as underlined in [17]. For example, these factors allow for fulfillment of the increasing heat flux and decreasing dimensions in electronics packages [17-24]. The common types of impinging jets are with or without confinement. Confined impinging jets have the advantages of smaller space design, while unconfined impinging jets have an advantage of simple design and easy fabrication. The two types of impinging jets have their own merits, and they are both commonly used as the cooling solutions, and the literature reviews on the subject have been provided in [2,3,6]. The effects of confinement on impinging jet heat transfer have been considered in [25-27]. Moreover, several studies show the importance of the subject and different cases have been investigated, such as confined slot-jet impingement on a moving plate [28], impinging jet on obliquely a flat surface [29], impinging jet on a porous medium [30] and slot jet impingement cooling on a semi-circular concave [31].

In order to obtain a heat transfer enhancement in jet impingement, different techniques have been employed, such as the insert of foams or fins [32]. These techniques determine a modification of the cooling system whereas the use of nanofluids in a coolant seems to be simpler in realizing a heat transfer enhancement [33]. However, nanofluids are to this day controversial in many areas such as inconsistencies in published data and disagreements on the heat transfer mechanisms, as observed by Gherasim et al. [34]. Various aspects of nanofluids have been covered in several reviews and some of these are given in [35-47].

The employment of nanofluids in impinging jets has been investigated recently by some researchers and, to the best of our knowledge, their investigations have been reported in [34,48-60]. The numerical investigation on hydrodynamic and thermal fields of Al2O3/water nanofluid in a radial laminar flow cooling system carried out by Roy et al. [48] can be can be considered as the first article on an impinging jet. Those authors found that considerable heat transfer enhancement was observed up to 200% in the case of a nanofluid with 10% in nanoparticle volume concentration at a Reynolds number equal to 1200. However, a significant increase in wall shear stress was noticed increasing the nanoparticle volume concentration. The laminar-forced convection flow of nanofluids between two coaxial and parallel disks with central axial injection was investigated numerically considering temperature-dependent properties by Palm et al. [51]. Results indicated a heat transfer benefit by adopting Al2O3/water nanofluid with a volume fraction of nanoparticles of 4%. An increase of 25% was evaluated in terms of average wall heat transfer coefficient, when referred to the water. Moreover, the use of temperature-dependent properties determined for greater heat transfer predictions with corresponding decreases in wall shear stresses when compared to evaluations employing constant properties. A numerical study on steady, laminar radial flow of a nanofluid in a simplified axi-symmetric configuration with axial coolant injection was performed by Roy et al. [52] for electronic cooling applications. Also in this investigation increases in heat removal capabilities were detected with the use of nanofluids.

An experimental investigation in a confined and submerged impinging jet on a flat, horizontal and circular heated surface with nanofluid (Al2O3 dispersed in water) was carried out by Nguyen et al. [56]. Experimental results were obtained for both laminar and turbulent flow regimes and they showed that, depending on the combination of nozzle-to-heated surface distance and particle volume fraction, the use of a nanofluid can determine a heat transfer enhancement in some cases, but an adverse effect on the convective heat transfer coefficient may occur in other cases. A circular confined and submerged jet impinging on a horizontal hot plate was numerically simulated by Vaziei and Abouali [57]. Water and 36-nm Al2O3-water nanofluid with various particle volume fractions were considered as a working fluid for cooling the hot plate. Both laminar and turbulent impinging jets in various nozzle-to-plate distances and Reynolds numbers were simulated. The results showed that the use of Al2O3 nanoparticles in laminar jets enhanced the heat transfer but for the turbulent jets Al2O3-water nanofluid had a lower performance for heat removal compared with the base fluid. The heat transfer enhancement capabilities of Al2O3/water inside a confined impinging jet cooling device was numerically studied by Gherasim et al. [34]. Results highlighted those limitations in the use of this nanofluid type in a radial flow configuration, due to the significant increase in the associated pumping power. Steady laminar incompressible thermal alumina-water flow between parallel disks was simulated by Feng and Kleinstreuer [58]. The results indicated that the Nusselt number increases with higher nanoparticle volume fraction, smaller nanoparticle diameter, reduced disk-spacing and larger inlet Reynolds number. The laminar forced convective heat transfer features of Al2O3/water nanofluid in the confined radial flow were numerically investigated by Yang and Lai [59,60] with constant [59] and temperature-dependent properties [60]. Results showed the same trend given in the previous published works: the Nusselt number increases with the increases in Reynolds number and nanoparticle volume fraction, though the increase in pressure drop is more significant with the increase of particle concentration. Furthermore, temperature-dependent thermo-physical properties of nanofluids were found to have a marked bearing on the simulation results.

It seems that a slot-confined and submerged impinging jet on a flat surface with nanofluids has not been investigated in both laminar and turbulent flow regimes in spite of its importance in engineering applications such as electronic cooling and material processing.

In this article, a numerical investigation on turbulent flow on a slot-confined and submerged impinging jet on an isothermal flat surface is carried out. The results are given to evaluate the fluid dynamic and thermal features of the considered geometry with Al2O3/water as the working nanofluid adopting the single phase model.

Methods

Geometrical model

A computational thermo-fluid dynamic analysis of a two-dimensional model, Figure 1a, which regards the impinging jet on a heated wall with nanofluids, is considered in order to evaluate the thermal and fluid-dynamic performances, and study the velocity and temperature fields. The two-dimensional model has a length L equal to 310 mm while the height H ranges from 24.8 to 124 mm and the jet orifice width W is 6.2 mm. A constant temperature value of 343 K is applied on the impingement bottom surface. Different values of H/W ratio, equal to 4, 6, 8, 10, 15 and 20, are considered. The working fluid is water or a mixture of water and γ-Al2O3 nanoparticles with a diameter of 38 nm, at different volume fractions equal to 1, 4 and 6%.

thumbnailFigure 1. Sketch of the model. (a) Geometry and boundary conditions; (b) details on the adopted mesh near the impingement surface.

Physical properties of nanofluids

The working fluid is water or a mixture of water and γ-Al2O3 nanoparticles with a diameter of 38 nm, at different volume fractions equal to 1, 4 and 6%. In Table 1, the values of density, specific heat, dynamic viscosity and thermal conductivity, given by Rohsenow et al. [61], are reported for water and γ-Al2O3. The presence of nanoparticles and their concentrations influence the mixture properties. A single-phase model was adopted, and the following equations were used for computing the thermal and physical properties of the considered nanofluids [62-65], given in Table 2. Density was evaluated using the classical formula developed for conventional solid-liquid mixtures, while the specific heat values were obtained by assuming the thermal equilibrium between particles and surrounding fluid [62,63].

Table 1. Material properties at a temperature of 293 K

Table 2. Properties of nanofluids, single-phase model.

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

(1)

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

(2)

Nanofluids may be considered as Newtonian fluids for low volume fractions, e.g., up to 10%, and for small temperature increases. In this way, for the viscosity as well as for thermal conductivity, formulas given by [64,65] were adopted:

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

(3)

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

(4)

However, it is well known that the evaluation of these properties by various research groups differs from each other because of the numerical and experimental approaches and processes adopted [64,65].

Mathematical description and governing equations

Steady-state, turbulent, incompressible, single-phase, and constant properties flow conditions are considered in the present analysis. The governing equations of continuity, momentum and energy are solved in rectangular coordinates:

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

(5)

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

(6)

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

(7)

where E is the total energy, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/188/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/188/mathml/M9">View MathML</a> and (τij)eff is the deviatoric stress tensor, defined as

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

(8)

The k-ε standard model with enhanced wall treatment is assumed. The transport equations are as follows [66]:

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

(9)

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

(10)

where Gk is the production of turbulent kinetic energy due to mean velocity gradients, Gb represents the generation of the turbulent kinetic energy due to buoyancy while YM is referred to the fluctuation rates related to the overall dissipated turbulent thermal energy. In particular, Gk may be expressed by

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

(11)

where C1ε and C2ε are constants; while the term <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/188/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/188/mathml/M14">View MathML</a> defines the dependence rate of ε on buoyancy; σk and σe represent the turbulent Prandtl numbers based on k and ε, respectively; while Sk and Sε are further generation terms. The turbulent viscosity is defined by

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

(12)

where Cμ is a constant. The model constant values are the following:

C1ε = 1.44, C2ε = 1.92, Cμ = 0.09, σk = 1.0 and σε = 1.3.

The enhanced wall treatment approach has been considered. The assigned boundary conditions are

- Inlet jet section: uniform velocity and temperature profile;

- Outlet section: pressure outlet;

- Bottom wall: velocity components equal to zero and constant temperature;

- Upper wall: velocity components equal to zero and adiabatic condition.

The dimensionless parameters considered here are

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

(13)

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

(14)

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

(15)

where uj is the jet velocity, W is the jet width, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/188/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/188/mathml/M19">View MathML</a> is the impingement surface heat flux, TH and TJ represent the temperature of the impingement surface and the jet temperature, respectively.

Numerical procedure

The governing equations of continuity, momentum and energy, reported in the previous section, are solved by the finite volume method by means of FLUENT code [67]. A steady-state solution and a segregated method are chosen to solve the governing equations, which are linearized implicitly with respect to dependent variables of the equation. A second-order upwind scheme is chosen for energy and momentum equations. The SIMPLE coupling is chosen as scheme to couple pressure and velocity. The convergence criteria of 10-5 for the residuals of the velocity components and of 10-8 for the residuals of the energy are assumed. It is assumed that the incoming flow is turbulent at ambient temperature and pressure. Different inlet uniform velocities, corresponding to Reynolds numbers ranging from 5000 to 20000, were considered and they are reported in Table 3. Furthermore, the inlet turbulence intensity value is set to 2%.

Table 3. Inlet velocities (m/s).

The enhanced wall treatment functions are activated to increase the model accuracy in the near-wall region. It is a two-layer method with enhanced functions. The domain is divided into two regions, the near-wall region and the core ones, according to the turbulent Reynolds number, based on the distance-to-wall term y.

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

(16)

The core region, for Rey > 200, is solved by means of the standard k-ε model, while in the other region the Wolfstein model is applied [68].

Along the solid walls, no slip condition is employed, whereas a velocity inlet is given for the jet orifice and pressure conditions are set for the outlet sections.

Four different grid distributions are tested on the model with H/W ratio equal to 6 at Re = 5000, with water (ϕ = 0%) as working fluid, to ensure that the calculated results are grid independent. The four grids have 4950 (90 × 55), 19800 (180 × 110), 79200 (360 × 220), and 316800 (720 × 440) nodes, respectively. The grid mesh is structured in each case with grid adoption for y+ = 1 at adjacent wall region and a sketch is shown in Figure 1b. For the adiabatic wall and the bottom surface, nodes are distributed by means of an exponential relation (n = 0.9), to have a fine mesh near the impingement region, where an equi-spatial distribution is chosen. On the vertical ones, a bi-exponential (n = 0.8) distribution is considered.

Comparing the third- and fourth-mesh configurations, in terms of average and stagnation point Nusselt number, results are very close, and the relative errors are very little, as reported in Figure 2. As a result, the third grid case has been adopted because it ensured a good compatibility between the machine computational time and the accuracy requirements.

thumbnailFigure 2. Analysis of grid independence in terms of Nusselt number, H/W = 6, R = 5000 and ϕ= 0%.

Results are validated by comparing the obtained numerical data with the experimental and numerical ones, given in [28,69,70]. Figure 3 presents the comparison in terms of average Nusselt number profiles, for the cases, characterized by Re = 11000, H/W = 6, TJ = 373 K and TH = 338 K. It is observed that the numerical results, obtained in this work, fit very well with the experimental ones given in [5,6] both near the stagnation point region and at the side one.

thumbnailFigure 3. Validation of numerical results in terms of Nusselt number.

Results and discussion

A computational thermo-fluid dynamic analysis of a two-dimensional model, regarding a confined impinging jet on a heated wall with nanofluids, is considered to evaluate the thermal and fluid-dynamic performances and study the velocity and temperature fields. Different inlet velocities are considered to ensure a turbulent regime, and the working fluids are water and mixtures of water and γ-Al2O3 at different volume fractions, treated by a single-phase model approach. The range of Reynolds numbers, geometric ratio and volume fractions are given below:

• Reynolds number, Re: 5000, 10000, 15000 and 20000;

H/W ratio: 4, 6, 8, 10, 15 and 20;

• particle concentrations, ϕ: 0, 1, 4 and 6%.

Results are presented in terms of average and local Nusselt number profiles, as a function of Reynolds number, H/W ratio and particle concentrations; moreover, dimensionless temperature fields and stream function contours are provided.

Figures 4 and 5 depict the stream lines contours and the temperature fields, respectively, for the representative cases with H/W = 4 and 10, at Re = 10000 and 20000 and ϕ = 0 and 4%. According to Figure 4, two counter-rotating vortex structures are generated as the jet impinges on the bottom surface and only one stagnation point, where velocity and temperature gradients are very high, is observed. This is due to the jet entrainment and confining effects of the upper adiabatic surfaces. Vortex intensity and size depend on H/W ratio, factors such as the confining effects, Reynolds number, and particle concentrations. It can be seen in Figure 4a, b, at Re = 10000, H/W = 10 and ϕ = 0 and 4%, the introduction of particles leads to a little smoother eddies with a low intensity increase, because the nanofluid viscosity is higher than water. As H/W ratio decreases from 10 to 4, at Re = 10000 and ϕ = 4%, vortices are less strong and smaller as they extinguish at x/W values equal to about -30 and 30, as pointed out in Figure 4b, c. As Re increases, the separation area near the inlet section becomes larger while the fluid stream results to be more compressed towards the impingement surface, as observed in Figure 4d.

thumbnailFigure 4. Stream functions contours. (a) H/W = 10, Re = 10000 and ϕ = 0%; (b) H/W = 10, Re = 10000 and ϕ = 4%; (c) H/W = 4, Re = 10000 and ϕ = 4%; (d) H/W = 10, Re = 20000 and ϕ = 4%.

thumbnailFigure 5. Temperature fields. (a) H/W = 10, Re = 10000 and ϕ = 0%; (b) H/W = 10, Re = 10000 and ϕ = 4%; (c) H/W = 4, Re = 10000 and ϕ = 4%; (d) H/W = 10, Re = 20000 and ϕ = 4%.

The temperature fields, depicted in Figure 5, follow the stream line patterns. For increasing concentrations, nanoparticles produce an increase of fluid bulk temperature, because of the elevated thermal conductivity of mixtures. Near the impingement surface, temperature grows and tends to decrease for increasing x/W values. For larger Reynolds numbers, the efficiency of heat transfer increases.

The variation of local Nusselt number along the impingement plate for Re = 20000, H/W = 4 and ϕ = 6% and for Re = 5000, H/W = 6 and different concentrations, is shown in Figure 6a, b, respectively. It is observed that the highest values of Nux are evaluated at the stagnation point for all the considered cases; their values are 214 and 239 for H/W = 4 and H/W = 10, respectively. For low H/W values, local Nusselt number decreases more quickly than high H/W ratios. At the end of the plate, for any considered H/W, Nux reaches similar values equal to about 25, as observed in Figure 6a. In Figure 6b, it is shown how the variation of nanofluid concentration affects the heat transfer. Higher heat transfer enhancements are observed for ϕ = 4, 6%, especially, near the impingement location. This does not happen only for H/W = 4 as can be understood from the average Nusselt number value trends, reported later, in comparison with other H/W ratios.

thumbnailFigure 6. Local Nusselt number profiles along x/W: (a) H/W = 4 and 10, Re = 20000, ϕ = 6%; (b) H/W = 6, Re = 5000, ϕ = 0, 1, 4 and 6%.

In Figure 7, the variation of local qw/q0w ratio is shown. The qw/q0w value represents the local ratio between the local total heat flux and total heat flux at stagnation point for any case. The maximum value is reached at the stagnation point of any considered case. As Re increases, qw/q0w ratio increases. Difference in terms of qw/q0w is more significant passing from Re = 5000 to 10000 than the other considered Re. In fact, at x/W = 4, there is a difference of 0.12 in terms of qw/q0w while in the other cases, the largest difference is 0.9. The heat transfer augmentation is more significant near the stagnation point than in correspondence with the end of the impinged plate. In Figure 7b, it is observed as the nanofluid concentration has very little influence on qw/q0w. The effects of H/W are underlined in Figure 7c: near the stagnation point, qw/q0w ratio has almost the same value for all H/W. From x/W = 4 curves spread out and qw/q0w increases as H/W increases. This affects the results in terms of average Nusselt number, calculated at different H/W ratios.

thumbnailFigure 7. Profiles of qw/q0w ratio along x/W: (a) H/W = 4, ϕ = 0%, Re = 5000, 10000, 15000 and 20000; (b) H/W = 4, Re = 5000, ϕ = 0, 1, 4 and 6%; (c) H/W = 4, 6, 8 and 10, ϕ = 0% and Re = 5000.

The average Nusselt number profiles as function of Re are depicted in Figure 8 for H/W = 4, 6, 8, and 10. Profiles increase as Re increases for all the considered cases. It is observed that as ϕ increases Nuavg becomes higher for a fixed value of Re. Passing from ϕ = 0% to ϕ = 1%, a significant increase of Nuavg, only for H/W = 4 configuration is noted, where it passes from 35 to 37 at Re = 15000 or 65 to 69 for Re = 20000. For the other cases, a significant heat transfer enhancement is found for the highest ϕ values; in fact in these cases, passing from ϕ = 0% to ϕ = 1%, the maximum enhancement is found to be equal to 1.22 times for H/W = 10 at Re = 20000.

thumbnailFigure 8. Average Nusselt number profiles as function of Re, ϕ= 0, 1, 4 and 6%: (a) H/W = 4; (b) H/W = 6; (c) H/W = 8; (d) H/W = 10.

The heat transfer enhancement is evident, also observing the average heat transfer coefficient profiles, described in Figure 9. Results are given for different Re, H/W ratios and concentrations. The maximum values of havg are calculated for the highest values of Re, H/W and concentrations considered. In fact, for H/W = 10 and Re = 20000, it results that havg is equal to about 7600, 8000, 9400 and 10500 W/m2K, as depicted in Figure 9d, while, at H/W = 4 and Re = 20000, havg are equal to about 6200, 6800, 7700, and 8600 W/m2K, for ϕ = 0, 1, 4, and 6%, respectively.

thumbnailFigure 9. Average convective heat transfer coefficient profiles as function of Re, ϕ= 0, 1, 4 and 6%: (a) H/W = 4; (b) H/W = 6; (c) H/W = 8; (d) H/W = 10.

Figure 10 shows the average Nusselt number profiles, referred to the values calculated for the base fluid, as a function of Reynolds number for particle concentrations equal to 1, 4 and 6% at H/W ratio of 4. It is observed that the ratio Nuavg/Nuavg, bf is greater than one for all the configurations analyzed and rises slightly for increasing Reynolds numbers and concentrations; in fact, the highest value of 1.18 is detected at Re = 20000 and ϕ = 6%.

thumbnailFigure 10. Profiles of Nuavg/Nuavg, bf ratio as a function of Re for different values of particle concentrations, H/W = 4.

The results in terms of local Nusselt numbers, calculated for the stagnation point, are depicted in Figure 11. They are provided as a function of Reynolds numbers and given for different concentrations for different H/W ratios, equal to 4, 6, 8, and 10. It is shown that profiles increase almost linearly with increasing Reynolds numbers for all the considered concentrations and H/W ratios. Moreover, the Nu0 values are the highest for ϕ = 6% for all the considered Reynolds numbers. For example, comparing the results for ϕ = 1, 4, and 6%, with the base fluid ones, an increase in values of 2.7, 10.8, and 16.2% are detected for H/W = 4 at Re = 20000, respectively. Moreover, Nu0 values rises as H/W increases for Re > 10000, as observed in Figure 11b, c, d.

thumbnailFigure 11. Stagnation point values of local Nusselt number. Values of local Nusselt number in correspondence with the stagnation point, for different Re and concentrations: (a) H/W = 4; (b) H/W = 6; (c) H/W = 8; (d) H/W = 10.

In fact, Figure 12 shows that Nu0 is maximum in correspondence with H/W = 4 for Re < 10000 and H/W = 10 for higher Reynolds numbers for all the concentrations. For ϕ = 0%, at Re = 5000 Nu0 values are about 70, 81, 86, and 87, while at Re = 20000, Nu0 = 195, 197, 200, and 205, for H/W = 4, 6, 8, and 10, respectively. The results for ϕ = 6% are depicted in Figure 12b; it is shown that at Re = 5000 the maximum value of the stagnation point Nusselt number is about 102, 100, 93, and 82, for H/W = 4, 6, 8, and 10, respectively. For the same geometrical configurations, at Re = 20000, Nu0 values are equal to 215, 225, 235, and 240.

thumbnailFigure 12. Stagnation point Nusselt number values as a function of Re. Nusselt number values of stagnation point as a function of Re, for different H/W ratios: (a) ϕ = 0%; (b) ϕ = 6%.

Results in terms of average Nusselt numbers are shown in Figure 13, for different H/W ratios and ϕ = 0, 6%. The profiles increase linearly as Re increases as well as H/W ratio. In fact, the highest values of Nuavg are detected for H/W = 10 while the minimum ones for H/W = 4. Moreover, average Nusselt numbers increase as ϕ increases; thus, Nuavg values are equal to 42 and 79 for water, as depicted in Figure 13a, while for ϕ = 6%, they are equal to 48 and 92, as pointed out by Figure 12b, at Re = 10000 and 20000, respectively.

thumbnailFigure 13. Average Nusselt number profiles as a function of Re for different H/W ratios: (a) ϕ = 0%; (b) ϕ = 6%.

Figure 14 confirms that the configurations with H/W = 10 exhibit the maximum values of the average Nusselt numbers for all the considered Reynolds numbers and concentrations. In fact, at Re = 5000 and 20000, the profiles increase as H/W rises until H/W = 10, and then they decrease for H/W = 15 and 20.

thumbnailFigure 14. Average Nusselt number profiles as a function of H/W for Re = 5000 and 20000, ϕ= 0 and 6%.

The pumping power is defined as PP = VΔP, and its profiles are shown in Figure 15, for all the considered H/W values, concentrations and as a function of Reynolds number. The required power has a square dependence on Re. It increases as H/W and particle concentration increase. For example, as observed in Figure 15a, at H/W = 4, for water PP = 15 and 90 W at Re = 10000 and 20000, respectively, while for ϕ = 6% PP = 50 and 410 W. At the same Re, for H/W = 10, PP is equal to 18 W, as underlined in Figure 15d, and 98 W for water, and 58 and 470 W for ϕ = 6%, respectively.

thumbnailFigure 15. Profiles of the required pumping power as a function of Re, for ϕ= 0, 1, 4 and 6%: (a) H/W = 4; (b) H/W = 6; (c) H/W = 8; (d) H/W = 10.

The pumping power ratio, referred to the base fluid values, is described in Figure 16. It is observed that the ratio does not seem to be dependent on Re, and PP/PPbf ratio increases as concentration increases, as expected. In fact, at Re = 15000, the required pumping power is 1.2, 2.6 and 4.8 times greater than the values calculated in case of water.

thumbnailFigure 16. Pumping power profile, referred to the base fluid values as a function of Re, ϕ = 1, 4 and 6%, H/W = 4.

Conclusions

A numerical analysis of a two-dimensional model on a confined impinging jet with nanofluids has been carried out to evaluate the thermal and fluid-dynamic performances and study the velocity and temperature fields. The bottom impinged wall is heated at a constant temperature and different fluid velocities are considered in the range 5000-20000. The base fluid is water and different volume concentrations of Al2O3 nanoparticles are taken into account by adopting a single-phase model approach. Furthermore, different H/W ratios have been studied. The dimensionless stream function contours showed that the vortex intensity and size depend on H/W ratio, such as on the confining effects, Reynolds number and particle concentrations. Furthermore, for increasing concentrations, nanofluids produce an increase of fluid bulk temperature, because of the elevated thermal conductivity of mixtures. The local Nusselt number profiles present the highest values at the stagnation point and the lowest at the end of the heated plate. The highest values of the average Nusselt numbers increase as the particle concentrations and Reynolds numbers increase and the highest values are observed for H/W = 10. A maximum increase of 18% is detected at ϕ = 6%. The required pumping power increases as well as Reynolds number, and particle concentrations grow, which is almost 4.8 times greater than the values calculated in the case of water. For list of symbols please see Table 4

Table 4. List of symbols

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All the authors have made substantial contributions in order to write this work. PM and DR developed the numerical model, ran the simulation and acquired data. The analysis and the interpretation of data have been carried out together with OM and SN. All the authors have been involved in drafting the manuscript and revising it critically and OM and SN have given final approval of the version to be published.

Acknowledgements

This study was supported by SUN Grant Initiative of 2008 grant and by MIUR with Articolo D.M. 593/2000 Grandi Laboratori "EliosLab".

References

  1. Martin H: Heat and mass transfer between impinging gas jets and solid surface.

    Adv Heat Transfer 1977, 13:1-60. Publisher Full Text OpenURL

  2. Downs SJ, James EH: Jet impingement heat transfer-a literature survey.

    Tech. Rep. 1987-HT-35, ASME OpenURL

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

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

  4. Viskanta R: Heat transfer to impinging isothermal gas and flame jets.

    Exp Thermal Fluid Sci 1993, 6:111-134. Publisher Full Text OpenURL

  5. Webb B, Ma CF: Single-phase liquid jet impingement heat transfer.

    Adv Heat Transfer 1995, 26:105-117. Publisher Full Text OpenURL

  6. Tesar V, Travnicek Z: Increasing heat and/or mass transfer rates in impinging jets.

    J Vis 2005, 8:91-98. Publisher Full Text OpenURL

  7. Ma CF, Bergles AE: Convective heat transfer on a small vertical heated surface in an impingement circular liquid jet. In Heat Transfer Science and Technology. Edited by Wang BX. New York: Hemisphere; 1990:193-200. OpenURL

  8. Ma CF: Fundamental research on extremely small size liquid jet impingement heat transfer.

    Proceedings of the 3rd International Thermal Energy Congress: 28 June-1 August 1997; Kitakyushu, Japan 1997, 195-202. OpenURL

  9. Jiji LM, Dagan Z: Experiment investigation of single phase multi-jet impingement cooling of an array of microelectronic heat sources. In Cooling Technology for Electronic Equipment. Edited by Aung W. New York: Hemisphere; 1988:332-352. OpenURL

  10. Punch J, Walsh E, Grimes R, Jeffers N, Kearney D: Jets and rotary flows for single-phase liquid cooling: an overview of some recent experimental findings.

    Proceedings of the 11th International Conference on Thermal, Mechanical and Multi-Physics Simulation, and Experiments in Microelectronics and Microsystems, EuroSimE 2010: 26-28 April 2010; Bordeaux 2010, 5464505. OpenURL

  11. Kohing FC: Waterwall: water cooling system.

    Iron-Steel Eng 1985, 62:30-36. OpenURL

  12. Ma CF, Yu J, Lei DH, Gan YP, Tsou FK, Auracher H: Transient jet impingement boiling heat transfer on hot surfaces. In Multiphase Flow and Heat Transfer-Second International Symposium. Edited by Chen XJ et al. New York: Hemisphere; 1990:349-357. OpenURL

  13. Viskanta R: Heat transfer in material processing. In Handbook of Heat Transfer. 3rd edition. Edited by Rohsenow WM, Hartnett JP, Cho YI. New York: McGraw-Hill; 1998. OpenURL

  14. Schuettenberg S, Krause F, Hunkel M, Zoch HW, Fritsching U: Quenching with fluid jets.

    Mater Werkst 2010, 40:408-413. Publisher Full Text OpenURL

  15. Liu Z, Lienhard VJH, Lombara JS: Convective heat transfer by impingement of circular liquid jets.

    J Heat Transfer 1991, 113:571-582. Publisher Full Text OpenURL

  16. Ma CF, Zhao YH, Masuoka T, Gomi T: Analytical study on impingement heat transfer with single-phase free-surface circular liquid jets.

    J Thermal Sci 1996, 5:271-277. Publisher Full Text OpenURL

  17. Chen YC, Ma CF, Qin M, Li YX: Theoretical study on impingement heat transfer with single-phase free-surface slot jets.

    Int J Heat Mass Transfer 2005, 48:3381-3386. Publisher Full Text OpenURL

  18. Zhuang Y, Ma CF, Qin M: Experimental study on local heat transfer with liquid impingement flow in two-dimensional micro-channels.

    Int J Heat Mass Transfer 1997, 40:4055-4059. Publisher Full Text OpenURL

  19. Lin ZH, Chou YJ, Hung YH: Heat transfer behaviors of a confined slot jet impingement.

    Int J Heat Mass Transfer 1997, 40:1095-1107. Publisher Full Text OpenURL

  20. Chiriac VA, Ortega A: A numerical study of the unsteady flow and heat transfer in a transitional confined slot jet impinging on an isothermal surface.

    Int J Heat Mass Transfer 2002, 45:1237-1248. Publisher Full Text OpenURL

  21. Park TH, Choi HG, Yoo JY, Kim SJ: Streamline upwind numerical simulation of two-dimensional confined impinging slot jets.

    Int J Heat Mass Transfer 2003, 46:251-262. Publisher Full Text OpenURL

  22. Sahoo D, Sharif MAR: Numerical modeling of slot-jet impingement cooling of a constant heat flux surface confined by a parallel wall.

    Int J Thermal Sci 2004, 43:877-887. Publisher Full Text OpenURL

  23. Lee HG, Yoon HS, Ha MY: A numerical investigation on the fluid flow and heat transfer in the confined impinging slot jet in the low Reynolds number region for different channel heights.

    Int J Heat Mass Transfer 2008, 51:4055-4059. Publisher Full Text OpenURL

  24. Sivasamy A, Selladurai V, Kanna PR: Jet impingement cooling of a constant heat flux horizontal surface in a confined porous medium: mixed convection regime.

    Int J Heat Mass Transfer 2010, 53:5847-5855. Publisher Full Text OpenURL

  25. Lytle D, Webb B: Air jet impingement heat transfer at low nozzle-plate spacings.

    Int J Heat Mass Transfer 1994, 37:1687-1697. Publisher Full Text OpenURL

  26. Behnia M, Parneix S, Shabany Y, Durbin PA: Numerical study of turbulent heat transfer in confined and unconfined impinging jets.

    Int J Heat Mass Transfer 1999, 20:1-9. OpenURL

  27. Choo KS, Kim SJ: Comparison of thermal characteristics of confined and unconfined impinging jets.

    Int J Heat Mass Transfer 2010, 53:3366-3371. Publisher Full Text OpenURL

  28. Sharif MAR, Banerjee A: Numerical analysis of heat transfer due to confined slot-jet impingement on a moving plate.

    Appl Thermal Eng 2009, 29:532-540. Publisher Full Text OpenURL

  29. Ibuki K, Umeda T, Fujimoto H, Takuda H: Heat transfer characteristics of a planar water jet impinging normally or obliquely on a flat surface at relatively low Reynolds numbers.

    Exp Thermal Fluid Sci 2009, 33:1226-1234. Publisher Full Text OpenURL

  30. Dórea FT, de Lemos MJS: Simulation of laminar impinging jet on a porous medium with a thermal non-equilibrium model.

    Int J Heat Mass Transfer 2010, 53:5089-5101. OpenURL

  31. Yang YT, Wei TC, Wang YH: Numerical study of turbulent slot jet impingement cooling on a semi-circular concave.

    Int J Heat Mass Transfer 2011, 54:482-489. Publisher Full Text OpenURL

  32. Naphon P, Wongwises S: Investigation on the jet liquid impingement heat transfer for the central processing unit of personal computers.

    Int Commun Heat Mass Transfer 2010, 37:822-826. Publisher Full Text OpenURL

  33. Choi SUS: Enhancing thermal conductivity of fluids with nanoparticles, developments and applications of non-Newtonian flows.

    ASME FED 1995, 231:99-105. OpenURL

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

    Int J Thermal Sci 2011, 50:369-377. Publisher Full Text OpenURL

  35. Das SK, Choi SUS, Patel HE: Heat transfer in nanofluids-a review.

    Heat Transfer Eng 2006, 27:3-19. Publisher Full Text OpenURL

  36. Buongiorno J: Convective transport in nanofluids.

    J Heat Transfer 2006, 128:240-250. Publisher Full Text OpenURL

  37. Daungthongsuk W, Wongwises S: A critical review of convective heat transfer of nanofluids.

    Renew Sustain Energy Rev 2007, 11:797-817. Publisher Full Text OpenURL

  38. Trisaksri V, Wongwises S: Critical review of heat transfer characteristics of nanofluids.

    Renew Sustain Energy Rev 2007, 11:512-523. Publisher Full Text OpenURL

  39. Yu W, France DM, Routbort JL, Choi SUS: Review and comparison of nanofluid thermal conductivity and heat Transfer enhancements.

    Heat Transfer Eng 2008, 29:432-460. Publisher Full Text OpenURL

  40. Keblinski P, Prasher R, Eapen J: Thermal conductance of nanofluids: is the controversy over?

    J Nanoparticle Res 2008, 10:1089-1097. Publisher Full Text OpenURL

  41. Kakaç S, Pramuanjaroenkij A: Review of convective heat transfer enhancement with nanofluids.

    Int J Heat Mass Transfer 2009, 52:3187-3196. OpenURL

  42. Choi SUS: Nanofluids: from vision to reality through research.

    J Heat Transfer 2009, 131:1-9. Publisher Full Text OpenURL

  43. Buongiorno J, et al.: A benchmark study on the thermal conductivity of nanofluids.

    J Appl Phys 2009, 106:1-9. Publisher Full Text OpenURL

  44. Özerinç S, Kakaç S, Yazıcıoğlu AG: Enhanced thermal conductivity of nanofluids: a state-of-the-art review.

    Microfluid Nanofluid 2010, 8:145-170. OpenURL

  45. Wang L, Fan J: Nanofluids research: key issues.

    Nanoscale Res Lett 2010, 5:1241-1252. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Godson L, Raja B, Lal DM, Wongwises S: Enhancement of heat transfer using nanofluids-an overview.

    Renew Sustain Energy Rev 2010, 14:629-641. Publisher Full Text OpenURL

  47. Venerus D, et al.: Viscosity measurements on colloidal dispersions (nanofluids) for heat transfer applications.

    Appl Rheol 2010, 20:445-462. OpenURL

  48. Roy G, Nguyen CT, Lajoie P: Numerical investigation of laminar flow and heat transfer in a radial flow cooling system with the use of nanofluids.

    Superlattices Microstruct 2004, 35:497-511. Publisher Full Text OpenURL

  49. Maiga S, Palm SJ, Nguyen CT, Roy G, Galanis N: Heat transfer enhancement by using nanofluids in forced convection flows.

    Int J Heat Fluid Flow 2005, 26:530-546. Publisher Full Text OpenURL

  50. Roy G, Palm SJ, Nguyen CT: Heat transfer and fluid flow of nanofluids in laminar radial flow cooling systems.

    J Thermal Sci 2005, 14:362-367. Publisher Full Text OpenURL

  51. Palm SJ, Roy G, Nguyen CT: Heat transfer enhancement with the use of nanofluids in radial flow cooling systems considering temperature-dependent properties.

    Appl Thermal Eng 2006, 26:2209-2218. Publisher Full Text OpenURL

  52. Roy G, Nguyen CT, Comeau M: Numerical investigation of electronic component cooling enhancement using nanofluids in a radial flow cooling system.

    J Enhanc Heat Transfer 2006, 13:101-115. Publisher Full Text OpenURL

  53. Schoeppler M: Correspondence address diverging ink-jet technologies and applications.

    Proceedings of the International Conference on Digital Printing Technologies: 16-17 September 2006; Denver 2006, 1-2. OpenURL

  54. Liu ZH, Qiu YH: Boiling heat transfer characteristics of nanofluids jet impingement on a plate surface.

    Heat Mass Transfer 2007, 43:699-706. Publisher Full Text OpenURL

  55. Liu ZH, Qiu YH: The boiling heat transfer of water based nanofluid jet impingement on a plate surface.

    J Shanghai Jiaotong Univ 2007, 41:1658-1661. OpenURL

  56. Nguyen CT, Galanis N, Polidori G, Fohanno S, Popa CV, Le Bechec A: An experimental study of a confined and submerged impinging jet heat transfer using Al2O3-water nanofluid.

    Int J Thermal Sci 2009, 48:401-411. Publisher Full Text OpenURL

  57. Vaziei P, Abouali O: Numerical study of fluid flow and heat transfer for Al2O3-water nanofluid impinging jet.

    Proceedings of the 7th International Conference on Nanochannels, Microchannels and Minichannels 2009: 22-24 June 2009; Pohang 2009, 977-984. OpenURL

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

    Int J Heat Mass Transfer 2010, 53:4619-4628. Publisher Full Text OpenURL

  59. Yang YT, Lai FH: Numerical study of heat transfer enhancement with the use of nanofluids in radial flow cooling system.

    Int J Heat Mass Transfer 2010, 53:5895-5904. Publisher Full Text OpenURL

  60. Yang YT, Lai FH: Numerical investigation of cooling performance with the use of Al2O3/water nanofluids in a radial flow system.

    Int J Thermal Sci 2011, 50:61-72. Publisher Full Text OpenURL

  61. Rohsenow WM, Hartnett JP, Cho YI: Handbook of Heat Transfer. 3rd edition. New York: McGraw-Hill; 1998. OpenURL

  62. Pak BC, Cho YI: Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles.

    Exp Heat Transfer 1998, 11:151-170. Publisher Full Text OpenURL

  63. Maiga SEB, Nguyen CT, Galanis N, Roy G: Heat transfer behaviours of nanofluids in a uniformly heated tube.

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

  64. Maiga SEB, Cong Tam N, Galanis N, Roy G, Mare T, Coqueux M: Heat transfer enhancement in turbulent tube flow using Al2O3 nanoparticle suspension.

    Int J Numer Methods Heat Fluid Flow 2006, 16:275-292. Publisher Full Text OpenURL

  65. Palm SJ, Roy G, Nguyen CT: Heat transfer enhancement with the use of nanofluids in radial flow cooling systems considering temperature dependent properties.

    Appl Thermal Eng 2006, 26:2209-2218. Publisher Full Text OpenURL

  66. Launder BE, Spalding DB: The numerical computation of turbulent flows.

    Comput Methods Appl Mech Eng 1974, 3:269-289. Publisher Full Text OpenURL

  67. FLUENT Computational Fluid Dynamic Code Version 6.3 User Guide, Fluent, Inc[http://www.fluent.com] webcite

  68. Wolfstein M: The velocity and temperature distribution of one-dimensional flow with turbulence augmentation and pressure gradient.

    Int J Heat Mass Transfer 1969, 12:301-318. Publisher Full Text OpenURL

  69. Cadek FFA: Fundamental investigation of jet Impingement heat transfer. In Ph.D Thesis. University of Cincinnati; 1968. OpenURL

  70. Gordon R, Akfirat JC: Heat transfer characteristics of impinging two-dimensional air jets.

    J Heat Transfer 1966, 88:101-108. OpenURL