Print version ISSN 0104-6632
Braz. J. Chem. Eng. vol.24 no.3 São Paulo July/Sept. 2007
FLUID DYNAMICS; HEAT AND MASS TRANSFER; AND OTHER TOPICS
V. C. MarianiI, *; L. S. CoelhoII
IPontifical Catholic University of Paraná, Graduate Program in Mechanical Engineering, PUCPR/CCET/PPGEM, Imaculada Conceição 1155, CEP: 80215-901, Curitiba PR, Brazil. E-mail: email@example.com
IIPontifical Catholic University of Paraná, Graduate Program in Industrial Systems and Engineering,Automation and Systems Laboratory, PUCPR/CCET/PPGEPS, Imaculada Conceição 1155, CEP: 80215-901, Curitiba - PR, Brazil. E-mail: firstname.lastname@example.org
A numerical study was conducted to investigate steady heat transfer and flow phenomena of natural convection of air in enclosures, with three aspect ratios (H/W = 1, 2, and 4), within which there is a local heat source on the bottom wall at three different positions, Wh. This heat source occupies 1% of the total volume of the enclosure. The vertical walls in the enclosures are insulated and there is an opening on the right wall. The natural convection is influenced by the difference in temperature between the left and right walls, represented by a Rayleigh number (Rae), and by local heat source, represented by a Rayleigh number (Rai). Numerical simulations were performed for several values of the Rayleigh number ranging between 103 and 106, while the intensity of the two effects the difference in temperature on the vertical walls and the local heat source was evaluated based on the Rai/Rae ratio in the range between 0 and 2500. The analysis proceeds by observing variations in the streamlines and isotherms with respect to the different Rae, R ratios, aspect ratios, of the radius and positions of the local heat source. The average Nusselt numbers on the hot and cold walls are influenced by different values of the parameters R, Rae, Wh, and H/W. Results show the presence of different flow patterns in the enclosures studied. Thus, the flow and heat transfer can be controlled by external heating, and local heat source.
Keywords: Natural convection; Average Nusselt number; Rectangular enclosures; Numerical study; heat source.
In recent years, numerical modeling of the convective heat transfer problem has been an area of great interest due to its broad applications in engineering. Compared to the experimental method, numerical analysis provides a more direct way to enhance/reduce heat transfer effectively so as to improve the performance or to optimize the structure of a thermal device.
Natural convection in enclosures has been studied both experimentally and numerically, due to the considerable interest in its many engineering applications, such as building insulation, solar energy collection, cooling of heat-generating components in the electrical and nuclear industries, and flows in rooms due to thermal energy sources (Yang, 1987).
Numerical studies of natural convection heat transfer and flow in closed enclosures without a local heat source are reported in the literature; we can cite the work of Davis (1983), Hortmann et al. (1990), Le Quéré (1991), Mohamad (1998), Corcione (2003), and Ben-Nakhi and Chamkha (2006).
Other authors have studied the natural convection caused by a heat-generating conducting body located inside an enclosure: Chu and Churchill, 1976; Khalilollahi and Sammakia, 1986; Keyhani et al., 1988; Farouk, 1988; Ho and Chang, 1994; Ha et al., 1999; Deng and Tang, 2002; Oztop et al., 2004; Bazylak et al., 2006.
Numerous studies on natural convection caused only by external heating in partially open enclosures have been conducted by Chan and Tien, 1985; Angirasa et al., 1995; Polat and Bilgen, 2002; Bilgen and Oztop, 2005; Lauriat and Desrayaud, 2006. However, few results have been reported for natural convection caused simultaneously by both external heating in partially open enclosures and an internal local heat sources although problems of this type are frequently important and their study is necessary for understanding the performance of complex natural convection flow and heat transfer.
Indirectly related to the present study, Xia and Zhou (1992) studied a square and partially open enclosure with an internal heat source. These authors change the position on the bottom wall or left vertical wall for only three R ratios. They found that the opening was advantageous to the flow and heat transfer in the cavity. In this case, the characteristics of flow and heat transfer changed with heat source location, external and internal Rayleigh number, and opening size. Reinehr et al. (2002) examined natural convection using the aspect ratio H/W = 2, with an internal heat source whose position was varied only on the bottom wall. In that work, no heat transfer results were reported and a limited number of Rae and R ratios were also studied.
The present work is a numerical study of natural convection due to the temperature difference between left and right walls and an internal local heat source in three partially open enclosures, for which few results have been reported in the literature. The enclosures have an opening in the cooled right vertical wall, while the left vertical wall is heated and the upper and lower walls are adiabatic. Natural convection is induced by the difference in temperature between the vertical walls, and it is represented by the Rayleigh number (Rae) and by an internal local heat source represented by the Rayleigh number (Rai), occupying approximately 1% of the enclosure volume.
The study is conducted numerically under the assumption of steady laminar flow for three different values of both the height-to-width aspect ratio of the enclosure of 1, 2, and 4 and the Rayleigh number based on enclosure height in the range between 103 and 106. The Rai/Rae ratio in the range between 0 and 2500 and the internal local heat source position at Wh = 0.25, 0.5, and 0.75 on the bottom wall are evaluated. In this context, the influence on flow patterns, temperature distributions and heat transfer rates is analyzed and discussed.
To model the flow under study, we use the conservation equations for mass, momentum, and energy for the two-dimensional, steady, and laminar flow. For the moderate temperature difference considered in this work, all the physical properties of the fluid, µ , k, and cp, are considered constant except density, in the buoyancy term, which obeys the Boussinesq approximation. In the energy conservation equation, we neglect the effects of compressibility and viscous dissipation. Thus, the dimensionless equations that govern the flow are
Definitions of the dimensionless parameters are listed in the Nomenclature section. The fluid in the interior environment is atmospheric air with the Prandtl number, Pr = 0.71 (air is the working fluid). The Rayleigh number (Rae) is represented by the difference in temperature between the vertical walls, 103 < Rae < 106. The intensity of heat produced by the source is represented by the Rayleigh number (Rai), which is based on the volumetric heat generation rate. The influence of the intensity of the two Rayleigh numbers is evaluated by means of the equation,
where 0 < R < 2500.
The stream function is defined by the following equations
where the numerical results illustrating the fluids behavior are presented in terms of the dimensionless stream function,
Other physical quantities of interest in the present study are the average Nusselt numbers for the hot and cold walls; these variables are defined respectively as
The air is studied in three rectangular enclosures with width W and height H heated on the left-hand vertical wall and cooled on the upper half of the right-hand vertical wall. These two walls have a prescribed temperature, while the horizontal walls are adiabatic. Schematic representations of the configurations studied are presented in Fig. 1.
The position Wh = 0.25, 0.5, and 0.75, illustrated in Fig. 1 is the distance from the left lateral wall to the center of the heat source. The internal heat generation source is located on the adiabatic bottom wall at different places and occupies 1% of the total volume of the enclosure. The lower half of the right wall is open and is in contact with the air outside the enclosure. The boundary conditions used are
where K is the ratio of the thermal conductivities between the heated source and the fluid, qs is the dimensionless temperature in the heated source, and qf is the dimensionless temperature in the fluid. The velocity and temperature profiles for the opening are assumed to be,
The differential equations, represented by equations (1) to (4), together with respective boundary conditions, equations (11) and (12), are solved using the finite volume method (MVF) described in Patankar (1980). In this method the solution domain is divided into small finite control volumes. The differential equations are integrated into each of those control volumes. From this integration there were algebraic equations which, when solved simultaneously or separately, supplied pressure and velocity components. A power law scheme is adopted for the convection-diffusion formulation. For the pressure-velocity coupling the SIMPLEC algorithm (semi implicit method for pressure linked equations consistency) is used (Patankar, 1980).
The discretized equations are solved iteratively, using the line-by-line method known as the Thomas algorithm or TDMA (tridiagonal matrix algorithm). An underrelaxation parameter of 0.5 was used in order to obtain a stable convergence for the solution of momentum and energy equations, while there was no need for such a parameter in the solution of the pressure equation.
Validation of the computer code for this work was verified for the natural convection problem in a closed enclosure without a local heat source. The results presented here are for two Rayleigh numbers, Rae = 105 and Rae = 106. Hortman et al. (1990) found the average Nusselt numbers 4.616 and 4.525 for Rae = 105 and grids 42 × 42 and 82 × 82, respectively, while in our study 4.604 and 4.535, respectively, were found. For the Rae = 106 they found the average Nusselt numbers 9.422 and 8.977 for the same grids, while in our study 9.487 and 8.975, respectively, were found.
Grid-independence tests were conducted for all the configurations studied in this work. Three different grid sizes (22 × 22, 42 × 42 and 82 × 82) were used and, for example, for H/W = 1, Rae = 2500, and R = 1000 the average Nusselt numbers for the hot wall obtained were -3.87, -4.01, and -4.06 and for the cold wall 4.45, 4.40, and 4.38, respectively, were found for the grids. Because of the small differences for the 42 × 42 and 82 × 82 grids, a 42 × 42 uniform grid was chosen for all the simulations presented in this work. Staggered storage of the variables was used. The numerical solution is considered to be converged when the maximum absolute value of the mass conservation was smaller than 10-10.
NUMERICAL RESULTS AND DISCUSSION
In order to compare the numerical code specifically developed for the present study, some solutions obtained in the square cavity, H/W = 1, when the internal heat source was centrally located at the bottom horizontal wall (at position Wh = 0.5) were compared with results of Xia and Zhou (1992) and Reinehr et al. (2002), showing good agreement. These results are presented in Table 1, where different flow patterns occur with the change in Rae or R. A comparative analysis between the values in Table 1 shows relative deviation (RD),
where qmax is proposed in Xia and Zhou (1992) and Reinehr et al. (2002) and Dqmax is the difference between the results of this work and those of the two studies mentioned. The results of Xia and Zhou (1992) are situated in the interval of 3.11% < qmax < 8.3%, while for Reinehr et al. (2002) the interval is 0.2% < qmax < 7%, showing good agreement between the results.
According to Table 2 the maximum dimensionless temperature increases considerably with the increasing R ratio, keeping Rae constant, independently of the position, heat sources and aspect ratio H/W. The maximum dimensionless temperature decreases with the increasing Rae when the R ratio is kept constant. A comparison of the maximum dimensionless temperatures at the same values of R and Rae, with the heat source at different positions inside the enclosures indicates that the differences are not significant. When the height of the enclosures increases, the dimensionless temperature also increases, as shown in Table 2.
In Table 3 the dependence of on R, Rae, H/W, and Wh is shown. In the range studied, undergoes a process from positive to negative with the increasing R for all enclosures studied. When the effect of external heating (Rae) on heat transfer in the partially open enclosures is larger than that of the local heat source, is positive and decreases with increasing R at a given Rae, and the hot wall plays a role in heating the fluid in the enclosures, but when Rae is smaller, is negative and || increases with R and the hot wall plays a role in cooling the fluid in the enclosures. This behavior was also noted by Xia and Zhou (1992). The changes with the different positions of local heat source, generally increasing with Wh. The || value increases with H/W. It can be observed that the value of Nuh generally doubles with the increasing of enclosures height from H/W = 1 to H/W = 2 and H/W = 2 to H/W = 4. In Table 4 the dependence of on R, Rae, H/W, and Wh is shown. In the range studied, increases with R and Rae. is generally larger when the local heat source is closer to the cold wall. When the enclosure height increases, increases significantly.
Streamlines and isotherms for R = 2500 are shown in Figs. 2 to 7 for all enclosures and Wh = 0.25, 0.5, and 0.75. An analysis of these figures indicates that when the heating source is located at the center, Wh = 0.5, or close to the opening of the enclosure, Wh = 0.75, various differences occur in the flow pattern according to the change in Rae or R. When the local heat source is located at the position Wh = 0.25 close to the hot wall, its influence on the flow patterns and on the heat transfer is negligible.
Verifying the position Wh = 0.5 in Figs. 2 and 3, it can be seen that for Rae = 103 two eddies appear inside the enclosures with H/W = 1 and 2 (see Fig. 2a and 2b) and three eddies appear inside the enclosure with H/W = 4 (see Fig. 2c). They circulate counterclockwise at the left of the enclosure and clockwise at the bottom right, where the intensity of both eddies is important for H/W = 1, i.e., flow and heat transfer are controlled by the internal local heat source and the difference in temperature on the vertical walls has a negligible influence. The left eddy close to the hot wall for H/W = 2 and 4 is less intense than the main flow because the local heat source becomes less important than the difference in temperature between the vertical walls.
In Figs. 4 and 5, for the position Wh = 0.75 it can be seen that at Rae = 104 two main eddies emerge in the enclosure, one circulating counterclockwise close to the hot wall and the other circulating clockwise close to the cold wall. In Figs. 4a and 5a the flow and heat transfer are controlled by the internal local heat source because the intensity of the counterclockwise eddy is much greater than that of the clockwise eddy, while in Figs. 4b, 4c, 5b, and 5c, the flow and heat transfer are controlled by the difference in temperature on the vertical walls.
In Figs. 6 and 7, it can be observed that at Wh = 0.25 for Rae = 105 the enclosures have only a very weak counterclockwise eddy at the backwind side of the local heat source and the larger eddy is a clockwise flow moving upwards along the hot wall and downwards along the cold wall, through the inflow and outflow openings. In these enclosures the flow and heat transfer are controlled by the difference in temperature of the vertical walls.
In Fig. 8 the fluid dynamic behavior of the air in the square enclosure (H/W = 1) is shown. Note that the increase in temperature differences between the vertical walls (Rae number) affects the fluid dynamic behavior, increasing the intensity of the flow in the enclosure. This behavior was observed for the various values of R employed, but is illustrated here only for R = 2500 and for the local heat source at the central position on the bottom wall. The main eddy increases with the Rae number, although it reaches the full height of the enclosure. In Fig. 9 some of the isotherms obtained for differents Rae are shown. Note that the thermal behavior changes as a function of the increasing Rae number.
The on the hot and cold walls are shown in Figs. 10a, 10b, and 10c, respectively, for H/W = 1, 2, and 4 at Wh = 0.5 and R = 2500, as a function of Rae. As shown in the isotherms in Fig. 3 and by the dimensionless temperature in Table 2, the fluid temperature decreases with increasing Rae for the same R in all enclosures studied in this work. The temperature gradient on the hot wall increases with Rae, as shown in Table 3. In this table for approximately R > 400, the has negative values, meaning that the fluid temperature next to the hot wall is higher than the hot wall temperature and the heat flow changes direction from the enclosure to the hot wall, while Nuc increases with Rae, meaning that the heat flow from the enclosure to the cold wall increases with this parameter. For example, and for Rae = 105 are larger than those for Rae = 104 for the same R and all enclosures, due to the increasing convective heat transfer with increasing Rae. The dimensionless stream functions ymax and |ymin| increase with Rae.
In Fig. 11 the variation in , , qmax, dimensionless ymax, and ymin with the R ratio is shown. In the range studied, the curves are similar, independent of the H/W value; nevertheless the values obtained are different. The increases with R, meaning that the heat flow from the enclosure to the cold wall increases with this parameter, and on the hot wall, undergoes a process from positive to negative with increasing R; this effect was shown in the results in Table 3. The dimensionless stream functions ymax and |ymin| increase with R. The other variables increase with R for fixed Rae. If the effect of external heating on heat transfer in the partially open enclosure is larger than that of the local heat source, then is positive and decreases with the increasing R at a given Rae.
Fig. 12 illustrates the variation in , , qmax, dimensionless ymax, and ymin with the location of local heat source for Rae = 104 and R = 1000. For this configuration, when the local heat source moves on the adiabatic wall (bottom wall) small changes in the average Nusselt number and in the other variables occur. When the heat source is located near the opening of the enclosures, the, (generally) and ymax values become larger in all enclosures studied. It can be seen in Fig. 12 that the values of all variables are modified smoothly with the change in aspect ratio of all enclosures.
This work was concerned with the numerical modeling of natural convection in three enclosures, within which there is a local heat source, occupying 1% of the total volume of the enclosures, located on the bottom wall at three different positions. The governing parameters were the Rayleigh number and the R ratio characterizing the heat transfer regime in natural convection. In view of the results, the findings may be summarized as follows:
i) with the increase in Rayleigh number (Rae), i. e., in the difference in temperature between the vertical walls, the maximum dimensionless temperature in the internal enclosures decreases, maintaining the R ratio constant, while and increase.
ii) with the increase in the R ratio, the maximum dimensionless temperature in the enclosure increases, maintaining the Rayleigh number constant for all positions of the local heat source analyzed, while decreases and increases.
iii) the position of the local heat source, Wh, influences the fluid dynamics of the air as well as the heat transfer rate in the enclosures.
iv) the aspect ratio H/W influences meaningfully the results obtained (for example, for maximum dimensionless temperature, small enclosures have lower temperatures, however large enclosures have higher temperatures for the same configurations) for , when the value is positive then it increases and when is negative it decreases with increasing height of the enclosure, while increases with increasing height of the enclosure.
v) The opening is advantageous for the flow and heat transfer in the enclosures, and its characteristics are complicated and change with location of internal heat source and external and internal Rayleigh numbers.
|cp||heat capacity||(J/kg K)|
|g||acceleration of gravity||(m/s2)|
|k||thermal conductivity||(W/m K)|
|K||ratio of the thermal conductivities between the heated source and the fluid||(-)|
|average Nusselt number||(-)|
|P||dimensionless pressure||(=(p + rgy)H3 / ra2)|
|Pr||Prandtl number||(= n/a)|
|q||rate of local heat generation by the heated protrusion||(W/m2)|
|R||Rayleigh number ratio||(-)|
|Rae||external Rayleigh number||(=(gb(Th - Tc)H3 / na))|
|Rai||internal Rayleigh number||(=(gbqH5 / naK))|
|u||velocity in x direction||(m/s)|
|v||velocity in y direction||(m/s)|
|U||dimensionless velocity in x direction||(= uH/a)|
|V||dimensionless velocity in y direction||(= vH/a)|
|Wh||heat generation source location||(m)|
|x, y||Cartesian coordinates||(m)|
|X, Y||dimensionless Cartesian coordinates||(= x/H, = y/H)|
|a||thermal diffusivity||(= k/rcp) (m2/s)|
|b||volumetric coefficient of thermal expansion||(K-1)|
|q||dimensionless temperature||(=(T-Tc)/(Th- Tc))|
|qs||dimensionless temperature in the heated source||(-)|
|qf||dimensionless temperature in the fluid||(-)|
|m||dynamic viscosity||(Pa s)|
Angirasa, D., Eggels, J. G. M. and Nieuwstadt, F. T. M., Numerical simulation of transient natural convection from an isothermal cavity open on a side, Numerical Heat Transfer, Part A: Applications, Vol. 28, No. 6, pp. 755-768 (1995). [ Links ]
Bazylak, A., Djilali, N. and Sinton, D., Natural convection in an enclosure with distributed heat sources, Numerical Heat Transfer, Part A: Applications, Vol. 49, No. 7, pp. 655-667 (2006). [ Links ]
Ben-Nakhi, A. and Chamkha, A. J., Effect of length and inclination of a thin fin on natural convection in a square enclosure, Numerical Heat Transfer, Part A: Applications, Vol. 50, No. 4, pp. 389-407 (2006). [ Links ]
Bilgen, E. and Oztop, H., Natural convection heat transfer in partially open inclined square cavities, International Journal of Heat and Mass Transfer, Vol. 48, No. 8, pp.1470-1479 (2005). [ Links ]
Chan, Y. L. and Tien, C. L., A numerical study of two-dimensional laminar natural convection in shallow open cavities, International Journal of Heat Mass Transfer, Vol. 28, No. 3, pp. 603-612 (1985). [ Links ]
Chu, H. H. and Churchill, S. W., The effect of heater size, location, aspect-ratio and boundary conditions on two-dimensional laminar natural convection in rectangular channels, Journal of Heat Transfer, Vol. 98, No. 2, pp. 194-201 (1976). [ Links ]
Corcione, M., Effects of the thermal boundary conditions at the sidewalls upon natural convection in rectangular enclosures heated from below and cooled from above, International Journal of Thermal Sciences, Vol. 42, No. 2, pp. 199-208 (2003). [ Links ]
Davis, G. de V., Natural convection of air in a square cavity: A benchmark solution, International Journal for Numerical Methods in Fluids, Vol. 3, No. 3, pp. 249-264 (1983). [ Links ]
Deng, Q. and Tang, G., Numerical visualization of mass and heat transport for conjugate natural convection/heat conduction by streamline and heatline, International Journal of Heat and Mass Transfer, Vol. 45, No. 11, pp. 2373-2385 (2002). [ Links ]
Farouk, B., Turbulent thermal convection in an enclosure with internal heat generation, Journal of Heat Transfer, Vol. 110, No. 1, pp. 126-132 (1988). [ Links ]
Ha, M. Y., Jung, M. J. and Kim, Y. S., Numerical study on transient heat transfer and fluid flow of natural convection in an enclosure with a heat-generating conducting body, Numerical Heat Transfer, Part A: Applications, Vol. 35, No. 4, pp. 415-433 (1999). [ Links ]
Ho, C. J. and Chang, J. Y., A study of natural convection heat transfer in a vertical rectangular enclosure with two-dimensional discrete heating: Effect of aspect ratio, International Journal of Heat and Mass Transfer, Vol. 37, No. 6, pp. 917-925 (1994). [ Links ]
Hortmann, M., Peric, M. and Scheuerer, G., Finite volume multigrid prediction of laminar natural convection: Benchmark solutions, International Journal for Numerical Methods in Fluids, Vol. 11, No. 2, pp. 189-207 (1990). [ Links ]
Keyhani, M., Prasad, V. and Cox, R., An experimental study of natural convection in a vertical cavity with discrete heat sources, Journal of Heat Transfer, Vol. 100, No. 3, pp. 616-624 (1988). [ Links ]
Khalilollahi, A. and Sammakia, B., Unsteady natural convection generated by a heated surface within an enclosure, Numerical Heat Transfer, Part A: Applications, Vol. 9, No. 6, pp. 715-730 (1986). [ Links ]
Lauriat, G. and Desrayaud, G., Effect of surface radiation on conjugate natural convection in partially open enclosures, International Journal of Thermal Sciences, Vol. 45, No. 4, pp. 335-346 (2006). [ Links ]
Le Quéré, P., Accurate solutions to the square thermally driven cavity at high Rayleigh number, Computers Fluids, Vol. 20, No. 1, pp. 29-41 (1991). [ Links ]
Mohamad, A. A., Benchmark solution for unsteady state CFD problems, Numerical Heat Transfer, Part A: Applications, Vol. 34, No. 6, pp. 653-672 (1998). [ Links ]
Oztop, H. F., Dagtekin, I. and Bahloul, A., Comparison of position of a heated thin plate located in a cavity for natural convection, International Comm. Heat Mass Transfer, Vol. 31, No. 1, pp. 121-132 (2004). [ Links ]
Patankar, S. V., Numerical heat transfer and fluid flow, Hemisphere Washington, D.C. (1980). [ Links ]
Polat, O. and Bilgen, E., Laminar natural convection in inclined open shallow cavities, International Journal of Thermal Sciences, Vol. 41, No. 4, pp. 360-368 (2002). [ Links ]
Reinehr, E. L., Souza, A. A. U. and Souza, S. M. A., Fluid dynamic behavior of air with natural convection and heat generation source in confined environment, In Proceedings of XIV Brazilian Congress of Chemical Engineering, Natal, Rio Grande do Norte, Brazil, pp. 1-8 (2002) (in Portuguese). [ Links ]
Xia, J. L. and Zhou, Z. W., Natural convection in an externally heated partially open cavity with a heated protrusion. FED-vol. 143/ [ Links ]HTD, Vol. 232, [ Links ]Measurement and Modeling of Environmental Flows - ASME, Vol. 232, pp. 201-208 (1992). [ Links ]
Yang, K. T., Natural convection in enclosures, in Handbook of Single Phase Convection Heat Transfer, Wiley, New York, NY, USA (1987). [ Links ]
(Received: May 30, 2005 ; Accepted: May 25, 2007)
* To whom correspondence should be addressed