Novel Analytical and Numerical Methods in Heat Transfer Enhancement and Thermal Management
View this Special IssueResearch Article  Open Access
Chithra Devaraj, Eswaramurthi Muthuswamy, Sundararaj Kandasamy, "Numerical Investigation of Laminar Natural Convection in Inclined Square Enclosure with the Influence of Discrete Heat Source", Journal of Applied Mathematics, vol. 2015, Article ID 985218, 9 pages, 2015. https://doi.org/10.1155/2015/985218
Numerical Investigation of Laminar Natural Convection in Inclined Square Enclosure with the Influence of Discrete Heat Source
Abstract
Natural convection heat transfer in a twodimensional square enclosure at various angles of inclination is investigated numerically using a finite volume based computational procedure. The heat transfer is from a constant temperature heat source of finite length centred at one of the walls to the cold wall on the opposite side while the remaining walls are insulated. The effect of area ratio of the heat source from 0.2 to 1.0, Rayleigh number Ra from 10^{3} to 10^{7}, and angle of inclination of the enclosure varying from 0° to 360° on the flow field and heat transfer characteristics are investigated. Streamline and isothermal line patterns are found to be similar at low Rayleigh numbers whereas at high Rayleigh numbers the differences are significant due to the influence of the parameters considered. Average Nusselt number decreases drastically as the position of the heat source is moved above the horizontal centre line of the enclosure. Correlation of the average Nusselt number which depends on the parameters of interest is obtained in the general form . The correlation coefficients are determined by multiple regression analysis for the entire range of Rayleigh numbers analysed and the values found by correlation equations are in good agreement with the numerical results.
1. Introduction
Natural convection heat transfer in enclosures has been extensively studied by researchers, because of its practical significance in science and technology. Applications include heating and cooling of building spaces, solar energy collectors, heat exchangers, and effective cooling of electronic components and machinery. The fluid flow and heat transfer behaviour of such systems are analysed numerically and experimentally by a number of researchers, incorporating different geometrical and thermal boundary conditions, from low to high Rayleigh numbers. The studies conducted by different groups were extensively summarized by Wan et al. [1] and a set of benchmark quality data was also presented in this study for the natural convection range of . On the other hand, several investigations [2–4] were carried out by the researchers on high Rayleigh numbers above 10^{8} where flow becomes turbulent for which numerical solutions are obtained by applying suitable turbulence models. Mao and Zhang [5] have carried out numerical simulation of turbulent natural convection of compressible air in tall cavities to evaluate the accuracy of turbulent models for various temperature boundary conditions. This study reveals that turbulence model is good in predicting velocity distribution while LES model is good in predicting temperature distribution at high Rayleigh numbers. Unsteady natural convection in a rectangular cavity, numerically investigated by Patterson and Imberger [6], shows that the transient flows strongly depend on the Prandtl number Pr and ultimately converge to identical steady state solutions. The numerical analysis conducted on convection in a cavity with variable viscosity fluid by Kandaswamy et al. [7] confirms that the heat transfer rate increases with the increase in viscosity of the fluid. Analytical studies [8, 9] on natural convection in a cavity with volumetric heat generation exemplify the convergence behaviour of various computational methods. A finite method approach by Saha et al. [10] for natural convection in a square tilt open cavity reveals the strong influence of the inclination angle of the cavity on convection heat transfer in addition to the influence of thermal boundary conditions. Corcione [11] has applied SIMPLER algorithm to determine the effect of thermal boundary conditions at the sidewalls on natural convection in a rectangular enclosure heated from below and cooled from above for the range of Rayleigh numbers between 10^{3} and 10^{6}. A study on laminar natural convection in shallow air filled rectangular cavities was presented by Samy and Elsherbiny [12] which showed strong dependency of natural convection heat transfer on Rayleigh number, enclosure aspect ratio, and angle of inclination of the enclosure. Aminossadati and Ghasemi [13] investigated the effects of orientation of an inclined enclosure with two adjacent sidewalls at different temperatures on natural convection heat transfer and observed that at high Rayleigh numbers the flow pattern and the rate of heat transfer depend on the angle of inclination of the enclosure. Kandaswamy and Eswaramurthi [14, 15] have focused their study on the determination of the effect of density maximum of water around 4°C on natural convection in a square porous enclosure applying finite volume method. The influences of porosity, Grashof number, Darcy number, and internal heat generation on natural convection heat transfer were presented in detail in these studies.
Numerical and experimental studies were conducted by different researchers by introducing geometrical features such as baffles, partition walls, fins, and conducting blocks proposing various heat transfer applications [16–19]. In each study, the configuration is described to affect the flow pattern within the enclosure which in turn influences the natural convection heat transfer. A number of studies on square enclosures with time dependent boundary conditions were found in the literature. For illustration, Aswatha et al. [20] applied finite volume method to study the effect of different thermal boundary conditions at bottom wall such as uniform/sinusoidal/linearly varying temperatures on natural convection in enclosures. The results presented from the abovementioned study show that uniform temperature gives higher Nusselt number compared to the sinusoidal and linearly varying bottom wall temperature cases.
In many practical situations, there is a need to view the natural convection in inclined enclosures where the heating component which can be treated as heat source does not spread over the entire length of the inner wall. The directions of gravitational and buoyancy forces remain the same, irrespective of the changes in the angle of inclination of the enclosure which alters the position of the heat source. The combined effect of size of the heat source and the angle of inclination causes additional changes in the flow pattern and the heat transfer characteristics. In the present study, natural convection phenomena in the enclosure for steady laminar flow are investigated for different values of area ratios of the heat source in the range of 0.2 to 1.0 and the angle of inclination varied from 0° to 360°. Each configuration is analysed for various laminar ranges of Rayleigh numbers, .
2. Mathematical Formulation
The schematic diagram of the two dimensional square enclosure of size , consisting of a discrete heat source of length at a constant temperature at one of its inner walls, is depicted in Figure 1. The wall opposite to the heat source acts as heat sink at constant temperature . Except the heat source and the sink, the remaining parts of the enclosure walls are assumed to be adiabatic. The origin of the coordinate system is fixed at one corner of the wall containing heat source. Numerical analyses are carried out for different area ratio of heat source equal to 0.2, 0.4, 0.6, 0.8, and 1.0 assuming that the enclosure contains fluid which is incompressible and laminar. Except the density which is approximated by the Boussinesq model, the fluid properties are assumed to be constant as elaborated by Gray and Giorgini [21]. Heat is transferred from the heat source to the cold wall. Thermally induced density gradient develops a buoyant flow and the gravitational force always acts on the fluid particles in the vertically downward direction irrespective of the angle of inclination of the enclosure.
2.1. Governing Equations
Conservation of mass, momentum, and energy equations are the governing equations to be solved for natural convection flow:
2.2. Boundary Conditions
Bottom wall: (; ): . Top wall: (; ): Heat source (; ): . Left wall lower side to heat source: (; ): . Left wall upper side to heat source: (; ): . Right wall: (; ): .
The properties of the fluid and the acceleration due to gravity are assigned with the appropriate values to obtain the desired Rayleigh number which is given byThe average Nusselt number over the cold wall, which is a measure of convection heat transfer in the enclosure, is estimated by using the relation
3. Computational Method and Grid Independence Study
In the present study, the finite volume based computational procedure is applied to solve the governing equations. The SIMPLE algorithm suggested by Patankar [22] is used for pressure velocity coupling and the second order upwind scheme is used for solving momentum and energy equations. Velocity gradients of the fluid would be high near the inner wall due to the boundary layer formation in comparison with the far field region. Rectangular form grid structure is created as illustrated in Figure 2 with fine grids close to the inner wall to capture the boundary layer accurately and with coarse grids near the centre of the enclosure to reduce the computational time. Grid independent test is conducted for grid sizes , , , , , , and for a typical square enclosure with heat source forming over the complete left wall (), top and bottom adiabatic walls, and complete right wall (cold surface) acting as heat sink. The variations of local Nusselt number over the cold wall for different grid structures at Rayleigh numbers 10^{4} and 10^{6} are shown in Figure 3. Grid comparison indicates that grid independent results are obtained for grid sizes and above. In order to ensure better accuracy of results, grid structure is selected as the ideal one for all the analyses in the present study. The convergence criteria were set to 10^{−6} for all the relative (scaled) residuals.
4. Validation Study
The methodology applied in the present study is validated rigorously by conducting similar studies reported by Wan et al. [1] as benchmark data for a wide range of Rayleigh numbers. Contour plots and plots of various parameters are compared and found to produce similar results between the two studies. The comparison of nondimensional horizontal velocities () at the midwidth () is presented in Figure 4 to demonstrate the attainment of validation study. Pointwise measurements for the Rayleigh numbers 10^{3} to 10^{6} are perceived to be in good agreement between the two results. The lack of comparison in the pointwise measurement begins to appear for Rayleigh number equal to 10^{7} at few corresponding points of the cavity. For example there is a maximum difference of 10% for Rayleigh number 10^{7} at 0.1 unit distances from the top and bottom walls of the cavity due to high velocity gradients at these levels. However, the present study is focused on the laminar range of Rayleigh numbers not exceeding 10^{7}.
5. Results and Discussion
The results of numerical investigation focused on the influence of area ratio of the heat source and the angle of inclination of the enclosure for Rayleigh numbers 10^{3} to 10^{7} are presented in this section. Contour plots are used for qualitative description considering typical cavities with heat source area ratios 0.4 and 0.8 and angles of inclination 0° and 60°. The influences on performance parameters are described quantitatively using plots.
5.1. Influence of Area Ratio and Angle of Inclination on Flow Pattern
Figure 5 shows the results associated with flow field in terms of streamlines for enclosure angles 0° and 60° with area ratios equal to 0.4 and 0.8, respectively, for different Rayleigh numbers. At low Rayleigh numbers the viscous forces are more predominant over buoyant forces due to which the flow circulations are weak. The influences of enclosure angle of inclination and the area ratio on streamline pattern are less significant at low Rayleigh numbers up to 10^{4} and the streamlines obtained are concentric forming single cell flow pattern. The development of the single cell is in the clockwise direction due to the rise of the fluid near the hot wall and consequent falling of the fluid near the cold wall. As Rayleigh number increases, the higher buoyant forces alter the flow pattern and the streamlines are closely spaced near the walls. It is an indication of faster moving fluid particles intensifying natural convection heat transfer. Formation of the streamline at high Rayleigh numbers is found to additionally depend on the area ratio of the heat source and the angle of inclination of the enclosure. For instance, for 0° angle of inclination of the enclosure the flow field tends to evolve into two counterrotating cells at Rayleigh number equal to 10^{5} where such evolution was initiated only at Rayleigh number equal to 10^{6} for 60° angle of inclination of the enclosure. As the Rayleigh number increases, the streamline spacing decreases and hence the magnitude of fluid velocity increases. When the enclosure is not inclined, the pattern of the streamline spacing is not symmetrical on the upper and lower sides of the heat source. More closely packed streamlines are noticed on the upper side of the hot surface. For 60° angle of inclination of enclosure the streamline spacing pattern appears symmetrically on the upper and lower sides of the hot surface.
Isothermal patterns for the area ratios 0.4 and 0.8 and the angles of inclination 0° and 60° in Figure 6 show the influence of these two parameters on temperature variation which in turn affects the rate of heat transfer. Thermal gradients are high near the hot and cold surfaces at all Rayleigh numbers. Isothermal lines are parallel to the hot and cold walls in the middle portion of the enclosure at low Rayleigh numbers and as the Rayleigh number increases these lines turn out to be almost perpendicular. The effect of angle of inclination of enclosure is noticeable in terms of isothermal lines plotted for 60° angle of inclination. The isothermal lines are closely packed near the adiabatic wall which is on the downstream side of the heat source.
The effect of area ratio for the same Rayleigh number and the angle of inclination can also be observed in Figures 5 and 6. Though the trends of the curves remain the same as the area ratio increases, additional concentration of the streamlines and isothermal lines is noticed, which indicates the increase in heat transfer rate with area ratio.
5.2. Variation of Average Nusselt Number
The angle of inclination of the enclosure strongly influences the natural convection heat transfer particularly at high Rayleigh numbers, due to change in streamline and isothermal line patterns. For illustration, the quantitative variations of average Nusselt number with for different Rayleigh numbers with hot surface area ratio equal to 0.8 are shown in Figure 7. At low Rayleigh numbers 10^{3} and 10^{4} the flow circulation within the enclosure is weak and hence average Nusselt number is independent of its angle of inclination. At high Rayleigh numbers the average Nusselt number decreases as increases from 0° to 90° and attains its initial value corresponding to 0° when rises from 90° to 180°. When is increased beyond 180° and approaches 270°, the hot surface ascends above the centre of the enclosure and the cold surface moves towards the bottom level. The buoyancy force deviating away from the cold wall contributes additional drop in rate of heat transfer. The lowest value of heat transfer rate is obtained at equal to 270° in which the hot and cold surfaces are on the extreme upper and lower sides, respectively, and the direction of the buoyancy force is exactly towards the hot wall. Increasing from 270° to 360° changes the average Nusselt number to its initial value which corresponds to 0° angle of inclination. The same trend of variation of with is observed for all the Rayleigh numbers investigated with different area ratios of the heat source.
5.3. Heat Transfer Dimensionless Correlations
In steady state condition all the heat flux () from the heat source is absorbed by the cold wall and hence the over the cold wall is a measure of heat transfer capacity of the enclosure. Based on the numerical results obtained in this study a correlation equation can establish the dependence of average Nusselt number on Rayleigh number and area ratio of the heat source, for a particular angle of inclination of the enclosure, as follows:
The values of correlation coefficients , , and which are obtained by multiple regression analysis for all the investigated angles of inclination of the enclosure together with the maximum relative error and the average standard deviation are detailed in Table 1. The values of and are estimated by the following equations as suggested by Chapra and Canale [23]:

The values of and vary between 0.0143 and 0.1523 and between 0.2216 and 3.567, respectively, which prove good agreement between the numerical results and the values obtained using correlation equation.
6. Conclusions
Natural convection heat transfer in a square enclosure heated from the discrete heat source and cooled from the opposite wall with the remaining walls insulated has been investigated. The studies are conducted for several area ratios of the heat source from 0.2 to 1.0, for the angles of inclination of the enclosure 0° to 360°, and for values of Rayleigh number in the range between 10^{3} and 10^{7}. Streamline and isothermal line patterns are found to be similar for all the cases at low Rayleigh numbers whereas considerable differences are observed at high Rayleigh numbers. Studies conducted on the influence of angle of inclination reveal that the average Nusselt number decreases drastically as the position of the heat source is moved above the horizontal centre line of the enclosure. Heat transfer correlation equations have been developed for each angle of inclination of enclosure investigated which are in good agreement with the numerical results.
Nomenclature
:  Area ratio of heat source () 
:  Acceleration due to gravity 
:  Length of the enclosure wall 
Nu:  Nusselt number 
:  Average Nusselt number 
:  Fluid static pressure 
Pr:  Prandtl number 
:  Heat flux 
Ra:  Rayleigh number 
:  Length of the heat source 
:  Temperature 
:  Temperature of cold wall 
:  Temperature of the heat source 
:  Velocity components in and directions 
:  Cartesian coordinates 
:  Correlation coefficients. 
:  Thermal diffusivity 
:  Volume expansion coefficient 
:  Angle of inclination of the enclosure 
:  Kinematic viscosity 
:  Fluid density 
:  Maximum relative error 
:  Average standard deviation. 
:  Cold wall 
:  Hot wall 
avg:  Average value 
corr:  Correlation. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 D. C. Wan, B. S. V. Patnaik, and G. W. Wei, “A new benchmark quality solution for the buoyancydriven cavity by discrete singular convolution,” Numerical Heat Transfer, Part B: Fundamentals, vol. 40, no. 3, pp. 199–228, 2001. View at: Publisher Site  Google Scholar
 N. C. Markatos and K. A. Pericleous, “Laminar and turbulent natural convection in an enclosed cavity,” International Journal of Heat and Mass Transfer, vol. 27, no. 5, pp. 755–772, 1984. View at: Publisher Site  Google Scholar
 P. Le Quere, “Accurate solutions to the square thermally driven cavity at high Rayleigh number,” Computers and Fluids, vol. 20, no. 1, pp. 29–41, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. Barakos, E. Mitsoulis, and D. Assimacopoulos, “Natural convection flow in a square cavity revisited: laminar and turbulent models with wall functions,” International Journal for Numerical Methods in Fluids, vol. 18, no. 7, pp. 695–719, 1994. View at: Publisher Site  Google Scholar
 Y. Mao and Y. Zhang, “Evaluation of turbulent models for natural convection of compressible air in a tall cavity,” Numerical Heat Transfer, Part B: Fundamentals, vol. 64, no. 5, pp. 351–364, 2013. View at: Publisher Site  Google Scholar
 J. Patterson and J. Imberger, “Unsteady natural convection in a rectangular cavity,” Journal of Fluid Mechanics, vol. 100, no. 1, pp. 65–86, 1980. View at: Publisher Site  Google Scholar
 P. Kandaswamy, S. Sivasankaran, M. Bhavanilakshmi, and R. S. Ramaswamy, “Convection in a cavity with variable viscosity fluid,” Indian Journal of Pure and Applied Mathematics, vol. 36, no. 1, pp. 641–650, 2005. View at: Google Scholar
 M. V. Joshi, U. N. Gaitonde, and S. K. Mitra, “Analytical study of natural convection in a cavity with volumetric heat generation,” Journal of Heat Transfer, vol. 128, no. 2, pp. 176–182, 2006. View at: Publisher Site  Google Scholar
 C. An, C. B. Vieira, and J. Su, “Integral transform solution of natural convection in a square cavity with volumetric heat generation,” Brazilian Journal of Chemical Engineering, vol. 30, no. 4, pp. 883–896, 2013. View at: Publisher Site  Google Scholar
 G. Saha, S. Saha, and A. H. Mamun, “A finite element method for steady state natural convection in a square tilt open cavity,” ARPN Journal of Engineering and Applied Sciences, vol. 2, no. 2, pp. 41–49, 2007. View at: Google Scholar
 M. Corcione, “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. View at: Publisher Site  Google Scholar
 M. Samy and Elsherbiny, “Laminar natural convection in shallow air rectangular cavities,” Alexandria Engineering Journal, vol. 43, no. 5, pp. 603–613, 2004. View at: Google Scholar
 S. M. Aminossadati and B. Ghasemi, “The effects of orientation of an inclined enclosure on laminar natural convection,” International Journal of Heat and Technology, vol. 23, no. 2, pp. 43–49, 2005. View at: Google Scholar
 P. Kandaswamy and M. Eswaramurthi, “Density maximum effect on buoyancydriven convection of water in a porous cavity with variable side wall temperatures,” International Journal of Heat and Mass Transfer, vol. 51, no. 78, pp. 1955–1961, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Eswaramurthi and P. Kandaswamy, “Transient doublediffusive convection of water around 4°C in a porous cavity,” Journal of Heat Transfer, vol. 131, no. 5, Article ID 052601, 7 pages, 2009. View at: Publisher Site  Google Scholar
 A. Himsar, K. Kouki, D. Mashasi, S. Takeo, T. Hiroshi, and S. Jun, “Laminar natural convection heat transfer in an air filled square cavity with two insulated baffles attached to its horizontal walls,” Thermal Science and Engineering, vol. 14, no. 3, pp. 35–46, 2006. View at: Google Scholar
 K. Kahveci, “Natural convection in a partitioned vertical enclosure heated with a uniform heat flux,” Journal of Heat Transfer, vol. 129, no. 6, pp. 717–726, 2007. View at: Publisher Site  Google Scholar
 X. Shi and J. M. Khodadadi, “Laminar natural convection heat transfer in a differentially heated square cavity due to a thin fin on the hot wall,” Journal of Heat Transfer, vol. 125, no. 4, pp. 624–634, 2003. View at: Publisher Site  Google Scholar
 J. F. Zou, Y. Gao, and W. K. Chow, “Numerical simulations on laminar natural convection in a square cavity with a conducting circular block,” International Journal of Heat and Technology, vol. 28, no. 1, pp. 1–8, 2010. View at: Google Scholar
 Aswatha, C. J. Gangadhara Gowda, S. N. Sridhara, and K. N. Seetharamu, “Effect of different thermal boundary conditions at bottom wall on natural convection in cavities,” Journal of Engineering Science and Technology, vol. 6, no. 1, pp. 109–130, 2011. View at: Google Scholar
 D. D. Gray and A. Giorgini, “The validity of the boussinesq approximation for liquids and gases,” International Journal of Heat and Mass Transfer, vol. 19, no. 5, pp. 545–551, 1976. View at: Publisher Site  Google Scholar
 S. V. Patankar, Numerical Heat Transfer and Fluid Flow, Hemisphere Publications, Washington, DC, USA, 1980.
 S. C. Chapra and R. P. Canale, Numerical Methods for Engineers, McGrawHill, Singapore, 3rd edition, 2008.
Copyright
Copyright © 2015 Chithra Devaraj et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.