The Scientific World Journal

Volume 2014 (2014), Article ID 617492, 11 pages

http://dx.doi.org/10.1155/2014/617492

Research Article

Natural Convection in a Differentially Heated Square Enclosure with a Solid Polygon

R. Roslan,¹ H. Saleh,² and I. Hashim2,3,4

¹Faculty of Science, Technology & Human Development, Universiti Tun Hussein Onn Malaysia, 86400 Parit Raja, Batu Pahat, Johor, Malaysia

²School of Mathematical Sciences, Faculty of Science & Technology, Universiti Kebangsaan Malaysia (UKM), 43600 Bangi Selangor, Malaysia

³Solar Energy Research Institute, Universiti Kebangsaan Malaysia (UKM), 43600 Bangi Selangor, Malaysia

⁴Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80257, Jeddah 21589, Saudi Arabia

Received 18 January 2014; Revised 27 April 2014; Accepted 20 May 2014; Published 2 June 2014

Academic Editor: Mohammed Hasnaoui

Copyright © 2014 R. Roslan 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.

Abstract.

The aim of the present numerical study is to analyze the conjugate natural convection heat transfer in a differentially heated square enclosure containing a conductive polygon object. The left wall is heated and the right wall is cooled, while the horizontal walls are kept adiabatic. The COMSOL Multiphysics software is applied to solve the dimensionless governing equations. The governing parameters considered are the polygon type, , the horizontal position, , the polygon size, , the thermal conductivity ratio, , and the Rayleigh number, . The critical size of the solid polygon was found exists at low conductivities. The heat transfer rate increases with the increase of the size of the solid polygon, until it reaches its maximum value. Here, the size of the solid polygon is reaches its critical value. Further, beyond this critical size of the solid polygon, will decrease the heat transfer rate.

1. Introduction

Natural convection heat transfer in differentially heated enclosures from side or below has received considerable attention over the past few decades, largely due to a wide variety of applications, which include solar collector technology, energy storage, nuclear reactor and technology. Another practical application of natural convection is encountered when an obstacle such as an inserted object is placed inside the enclosure. In applications involving building energy components, such as walls, windows, or an electronic enclosure, the inserted solid object may reduce the flow, thereby reducing the heat transfer rate across the enclosure. On the other hand, heat transfer may be enhanced if the solid object has a relatively high thermal conductivity.

House et al. [1] firstly examined the effect of a centered heat-conducting object on natural convection heat transfer in a vertical enclosure for various Rayleigh number, Prandtl number, object size, and ratio of thermal conductivities. They obtained that the heat transfer may be enhanced or reduced by a square object with a thermal conductivity ratio less or greater than unity. Later, Liu and Phan-Thien [2] included a radiation effect. Oh et al. [3], Deng and Tang [4], and Zhao et al. [5] visualized the problem by the heatlines concept. Merrikh and Mohamad [6] found that placing the solid objects near to the walls reduces the rate of heat transfer due to the blockage effects, but placing low conductor objects far from the boundary layer region may enhance the rate of heat transfer compared with enclosures without obstacles. Ha et al. [7] used a spectral multidomain methodology to handle a square object located at the center and concluded that the Rayleigh number for the fluid and temperatures fields to become nonsymmetrical and time dependent depends on the thermal boundary conditions of the the object. Bhave et al. [8] reported an existence of the critical size of the adiabatic object. The heat transfer increases with the increase of the size of the object, until it reaches a critical size, where the heat transfer is at maximum value. Beyond this critical size of the adiabatic object, will decrease the heat transfer rate. Das and Reddy [9] and Aminossadati and Ghasemi [10] studied conjugate natural convection in enclosures with a given inclination angle. They obtained the higher solid conductivity results in a better heat transfer and increasing the inclination angle improves the heat transfer performance.

Fu et al. [11] placed the isothermal cylinder object in the left half enclosure. They found a heat transfer enhancement by inserting a rotation object. Costa and Raimundo [12] considered a conductive rotating cylinder inserted in the center of a square enclosure. They concluded that the thermophysical properties of the cylinder object were important on the overall heat transfer process across the enclosure. Hussain and Hussein [13] analyzed effect of inserting a conductive rotating cylinder at different vertical locations inside a differentially heated square cavity. Besides, inserting the quadrilateral or cylinder objects, Shih et al. [14] inserted the adiabatic triangle object in the enclosure. They concluded that a bigger triangle size exhibits the highest heat transfer performance.

Present work aims to investigate the fluid flow and heat transfer characteristics for various polygon types and thermal properties placed inside the center of the square enclosure. Complete two-dimensional numerical simulation and systematical generalization of the conjugate heat transfer behavior occurring in the enclosure by varying the obstacles’ shape is carried out.

2. Mathematical Formulation

Figure 1(a) presents the coordinate systems and a square enclosure having a conductive regular polygon placed at (). It is assumed that the dimension in the -direction is large enough and the end effects on the flow are negligible; that is, fluid flow and heat transfer are two-dimensional. The vertical walls of both enclosures are maintained constant and at uniform different levels of temperature, and . The horizontal walls are kept insulated. Free space between polygons and walls enclosure are filled with the Newtonian fluid which is incompressible but expands or contracts with temperature changes. This assumption leads to the use of the Boussinesq approximation. Under the above assumptions, thermal radiation is neglected; the governing equations for steady natural convection flow using conservation of mass, momentum, and energy can be written as

and the energy equation for the solid polygon is

where the subscripts and stand for the fluid and the solid, respectively. No-slip condition is assumed at all the solid-fluid interfaces. Using the following nondimensional variables:

The resulting nondimensional forms of (1)-(2) are

The values of the nondimensional velocity are zero in the solid region and on the solid-fluid interfaces. The boundary conditions for the nondimensional temperatures are:

where is the thermal conductivity ratio. All sides of the solid polygon are equal in length (regular polygon). Number of sides is denoted by , where for the shape is triangle, for the shape is square, for the shape is pentagon, for the shape is hexagon, for the shape is heptagon, for the shape is octagon, and for the shape becomes cylinder. The polygon to enclosure area ratio is defined by

The fluid motion is displayed using the stream function obtained from velocity components and . The relationships between the stream function and the velocity components are and , which yield a single equation as

(a)

(b)

Figure 1: (a) Schematic representation of the model. (b) Mesh distribution.

The physical quantity of interest in this problem is the heat transfer rate. The heat transfer rate distribution is obtained by applying Fourier’s law at the hot wall; that is,

which becomes in terms of the dimensionless variables

where is local Nusselt number based on . The average Nusselt number on the heated left wall is then found by integrating the local distributions; that is, by using

3. Computational Methodology

The governing equations along with the boundary conditions are solved numerically by the CFD software package COMSOL Multiphysics. COMSOL Multiphysics (formerly FEMLAB) is a finite element analysis, solver, and simulation software package for various physics and engineering applications. We consider the following application modes in COMSOL Multiphysics. The Incompressible Navier-Stokes application mode, (ns), the Convection-Conduction application mode (cc) and the Diffusion application mode (di) are used for equations (4), (5), and (6), respectively. P2-P1 Lagrange elements and the Galerkin least-square method are used to assure stability. A parallel direct solver (PARDISO) and the damped Newton method are implemented to solve the discretized equations. Here the convergence criterion is set to .

In this study, mesh generation on square enclosure containing polygon object is made by using triangles. The triangular mesh distribution is shown in Figure 1(b). Several grid sensitivity tests were conducted to determine whether the mesh scheme is sufficient and to ensure that the results are grid independent. We use the COMSOL default settings for predefined mesh sizes, that is, extremely coarse, extra coarse, coarser, coarse, normal, fine, finer, extra fine, and extremely fine. In the tests, we consider the parameters ,  , , , and as shown in Table 1. Considering the accuracy and CPU time, a finer mesh size was selected for all the computations done in this paper.

Table 1: Grid sensitivity check at , , , , and .

Predefined mesh size
Mesh elements

CPU time (s)

Extremely coarse
261
4.4311
2

Extra coarse
459
4.4860
2

Coarser
677
4.5051
2

Coarse
1309
4.5101
2

Normal
1973
4.5120
2

Fine
2969
4.5162
2

Finer
8871
4.5193
3

Extra fine
24827
4.5195
6

Extremely fine
32411
4.5201
7

To validate the computational code, the previously published problems on conjugate natural convection in a square enclosure with a conducting quadrilateral and cylinder objects were solved. Figure 2 shows the comparison of the present computed heatlines (a, c) against that of Zhao et al. [5] for the quadrilateral, (b) and Costa and Raimundo [12] for the cylinder, (d) at , , and . The comparison reveals good agreement with the reported studies. These comprehensive verification efforts demonstrated the robustness and accuracy of the present computation.

(a)

(b)

(c)

(d)

Figure 2: Comparison of the present computed heatlines (a, c) against that of Zhao et al. [5] for quadrilateral, (b) and Costa and Raimundo [12] for cylinder, (d) at , , and .

4. Results and Discussion

The analyses in the undergoing numerical investigation are performed in the following range of the associated dimensionless groups: the solid polygon to enclosure area ratio, , the horizontal position, , the thermal conductivity ratio, , and the Rayleigh number, . The polygon shapes vary from triangle to cylinder, where cylinder is categorized as a special polygon having infinity sides.

Figure 3 illustrates the streamlines for various types of solid polygon, where the solid area () attains the value , , and , respectively. The thermal conductivity ratio is fixed at and and the Rayleigh number at . The fluid temperature adjoining the hot surface rises and moves from the left to the right, falling along the cold surface, then rising again at the hot surface. This movement creates a clockwise circulation cell in free space between the polygon and walls enclosure. The cell shape adjacent to the outer polygon surface was bent by the presence triangle, quadrilateral and pentagon. The bent takes higher position as the polygon size is made bigger. The strength of the flow circulations decreases by increasing the solid area for the same polygon type. It is obvious that increasing the leads to smaller space for the flow to circulate. At fixed , increasing increases the flow strength, as indicated from values. This increasing occurs except at small cylinder, , where the flow strength is weaker than the small solid pentagon inserted inside in the enclosure.

(a)

(b)

(c)

(d)

Figure 3: Streamlines for different polygon size and number of polygon sides at , , and .

The corresponding heatlines, presented in Figure 4, show that open lines starting from the left hot surface and reaching the right cold surface are responsible for the heat exchange between the system and the environment. The heatlines crossing the hot surface are more crowded near the bottom side than those near the topside. Many of the heatlines follow the fluid flow. The heatlines indicate that heat is conducted vertically through the solid polygon. The closed-loop occurs at the various sizes and polygon types. We observed dual closed loop occurs in enclosure containing a small pentagon or a cylinder. The negative value in the closed-loop cells correspond to a passive zone in which heat rotates with small contribution to the heat transfer between the active walls. The bigger polygon sizes are observed to create greater suppression of the inner heat circulation as indicated from values. remains constant by increasing at relative small solid area. Increasing from 3 to 4 leads to a decreasing of the at medium size solid polygon. At relative large solid area, the values increases by increasing . The values remain constant by adjusting the and . The is located at the upper part of the enclosure; this refers to the total amount of the system heat transfer.

(a)

(b)

(c)

(d)

Figure 4: Heatlines for different polygons to enclosure area ratio and number of polygon sides at , , and .

Figure 5 presents the streamlines, isotherms, and heatlines to show the effects of thermal conductivity ratio on flow fields, thermal and heat path characteristics for solid pentagon, , , and . The streamlines shown in the left of Figure 5 are all circulated as vortices, since there is no mass or fluid exchange between the system of the enclosure and its environment. The flow strength measured by increases with increasing thermal conductivity ratio . The isotherms in the solid pentagon are sparser, while the isotherms in the fluid remain unchanged by increasing the solid conductivity. The inner heat circulation is enhanced when the solid conductivity is strengthened. The decreases as the solid conductivity is strengthened.

(a)

(b)

(c)

(d)

Figure 5: Streamlines (left), isotherms (middle), and heatlines (right) evolutions by varying thermal conductivity ratio for , , , and .

Figure 6 presents the streamlines, isotherms, and heatlines to show the effects of Rayleigh number on flow fields, thermal and heat path characteristics for solid pentagon, , , and . At low Rayleigh number, the streamlines present a circular shape and are distributed uniformly. This is a signal of very weak convection. Consequently, the heat transfer across the enclosure is dominantly due to conduction process and because of that the heatlines and the isotherms both present pseudoconduction pattern. The flow strength measured by increases with increasing the Rayleigh number. The temperature stratification is step by step developed. At , the nonuniform distribution of streamlines is observed. The circular cells are elongated and indicate that convection is strong in the outer region and weak in the core. The isotherms in the solid pentagon are distributed horizontally, while the isotherms in the fluid at the isothermal walls are denser than before. The sparse heatlines in the solid pentagon indicate that the role of conduction is markedly weakened in the process of heat transfer. At strong buoyancy force, , the streamlines are grouping in the thin outer region and sparse in the inner region. Crowded vertical thermal boundary layers are formed at isothermal walls. The flow circulation, inner heat circulation, and maximum heatlines values as presented in Figure 6 are enhanced as the buoyancy force is taken higher. The vortex of clockwise closed-loop cells moves from lower-middle region to upper-left region by increasing .

(a)

(b)

(c)

(d)

Figure 6: Streamlines (left), isotherms (middle), and heatlines (right) evolutions by varying Rayleigh number for , , , and .

Variations of the average Nusselt number with the horizontal position are shown in Figure 7 for , a polygon size, , and different values of the Rayleigh number. The Nusselt number is seen to increase with increasing the Rayleigh number due to the increasing effect of convection. In general, when the inner polygon is displaced to the center of the cavity, the increases due to the weakening of the interaction between the active walls and the solid polygon which becomes located outside the thermal boundary layers. For , the heat transfer is dominated by conduction and the Nusselt number remains unchanged by changing the polygon’s position. For , presents a parabolic profile with a maximum value at . At , exhibits a flat profile in the range . This flat profile becomes wider for , as a consequence of the development of the thermal boundary regime for this high value of .

Figure 7: Variation of with for different values of at and .

Variations of the average Nusselt number with the number of polygon sides are shown in Figure 8 for different values of the Rayleigh number at and and polygon size, . We applied the inverse of formula (8) to maintain the area of each polygon equivalent. Applying formula (8) to an input and gives a result of ; then applying the inverse formula (8) to gives the result and . and are inversely proportional, a longer will be used for a lower or vice versa. Figure 8 shows that the average Nusselt number is unchanged by varying the polygon types at the considered values of . Careful investigation reveals that the values shown in Figure 6 are equivalent with the values in this figure. This is a proof that the overall heat transfer performance can be read directly from the positive heatline value at the upper part of the enclosure.

Figure 8: Variation of with for different values of at , , and .

Variations of the average Nusselt number with the number of polygon sides are shown in Figure 9 for different thermal conductivity ratios at , , and . Increasing from 3 to 4 leads to a decreasing of the average Nusselt number. This means that the square object leads to greater suppression of the inner heat circulation than the triangle object as indicated from the (see Figure 4). The average Nusselt number was shown to be stable for . Increasing thermal conductivity ratio leads to a decrease of the Nusselt number for a fixed value of .

Figure 9: Variation of with for different values of at , , and .

Variations of the average Nusselt number with the polygon sizes are shown in Figure 10 for different thermal conductivity ratio at , , and . Figure 10 exhibits the existence of a critical size of the solid pentagon at low conductivities, , ; below which, the increase of the size increases the and above which the increase of the size decreases the . The is observed to decrease by increasing the solid size having high conductivities, .

Figure 10: Variation of with for different values of at , , and .

5. Conclusions

The present numerical simulations study the effects of various solid polygons properties on natural convection inside a square enclosure. The dimensionless forms of the governing equations were solved using the COMSOL Multiphysics software. Detailed computational results for fluid flow, temperature distribution, heat path, and heat transfer characteristics in the enclosure have been presented in graphical forms. The main conclusions of the present analysis are as follows.

(1)The strength of the flow and inner heat circulations is much higher for greater . The maximum heatlines value is not sensitive by changing the .(2)The heat transfer rate remains stable for and polygons located at the center of the enclosure will assure the maximum heat transfer rate.(3)The critical size of the solid polygon was found exists at low conductivities. The heat transfer rate increases with the increase of the size of the solid polygon, until it reaches a maximum value. At this moment, the size of the solid polygon is reaches its critical value. Further, beyond this critical size of the solid polygon, will decrease the heat transfer rate.

The theoretical prediction in this paper is hoped to be a useful guide for the experimentalists to study the various combinations of the polygon shaped and its thermal conductivity properties to control the fluid flow and thermal performance of enclosure at different sizes. The factors of polygon location, orientation, and rotation with different angular speed will be the focus of our research undertaking.

Nomenclature

: | Polygon to enclosure area ratio
: | Gravitational acceleration
: | Thermal conductivity ratio
: | Thermal conductivity
: | Width and height of enclosure
: | Number of polygon sides
: | Average Nusselt number
: | Radius of polygon & dimensionless radius of polygon
: | Rayleigh number
: | Temperature
, : | Velocity components in the - and -directions
, , : | Space coordinates & dimensionless space coordinates.

Greek Symbols

: | Thermal diffusivity
: | Thermal expansion coefficient
: | Dimensionless temperature
: | Kinematic viscosity.

Subscript

: | Cold
: | Fluid
: | Hot
: | Solid.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

J. M. House, C. Beckermann, and T. F. Smith, “Effect of a centered conducting body on natural convection heat transfer in an enclosure,” Numerical Heat Transfer A: Applications, vol. 18, no. 2, pp. 213–225, 1990.

Y. Liu and N. Phan-Thien, “Complete conjugate conduction convection and radiation problem for a heated block in a vertical differentially heated square enclosure,” Computational Mechanics, vol. 24, no. 3, pp. 175–186, 1999.

J. Y. Oh, M. Y. Ha, and K. C. Kim, “Numerical study of heat transfer and flow of natural convection in an enclosure with a heat-generating conducting body,” Numerical Heat Transfer A: Applications, vol. 31, no. 3, pp. 289–303, 1997.

Q.-H. Deng and G.-F. Tang, “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.

F.-Y. Zhao, D. Liu, and G.-F. Tang, “Conjugate heat transfer in square enclosures,” Heat and Mass Transfer, vol. 43, no. 9, pp. 907–922, 2007.

A. A. Merrikh and A. A. Mohamad, “Blockage effects in natural convection in differentially heated enclosures,” Journal of Enhanced Heat Transfer, vol. 8, no. 1, pp. 55–72, 2001.

M. Y. Ha, I.-K. Kim, H. S. Yoon et al., “Two-dimensional and unsteady natural convection in a horizontal enclosure with a square body,” Numerical Heat Transfer A: Applications, vol. 41, no. 2, pp. 183–210, 2002.

P. Bhave, A. Narasimhan, and D. A. S. Rees, “Natural convection heat transfer enhancement using adiabatic block: optimal block size and Prandtl number effect,” International Journal of Heat and Mass Transfer, vol. 49, no. 21-22, pp. 3807–3818, 2006.

M. K. Das and K. S. K. Reddy, “Conjugate natural convection heat transfer in an inclined square cavity containing a conducting block,” International Journal of Heat and Mass Transfer, vol. 49, no. 25-26, pp. 4987–5000, 2006.

S. M. Aminossadati and B. Ghasemi, “Conjugate natural convection in an inclined nanofluid-filled enclosure,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 22, no. 4, pp. 403–423, 2012.

W.-S. Fu, C.-S. Cheng, and W.-J. Shieh, “Enhancement of natural convection heat transfer of an enclosure by a rotating circular cylinder,” International Journal of Heat and Mass Transfer, vol. 37, no. 13, pp. 1885–1897, 1994.

V. A. F. Costa and A. M. Raimundo, “Steady mixed convection in a differentially heated square enclosure with an active rotating circular cylinder,” International Journal of Heat and Mass Transfer, vol. 53, no. 5-6, pp. 1208–1219, 2010.

S. H. Hussain and A. K. Hussein, “Mixed convection heat transfer in a differentially heated square enclosure with a conductive rotating circular cylinder at different vertical locations,” International Communications in Heat and Mass Transfer, vol. 38, no. 2, pp. 263–274, 2011.

Y.-C. Shih, J. M. Khodadadi, K.-H. Weng, and A. Ahmed, “Periodic fluid flow and heat transfer in a square cavity due to an insulated or isothermal rotating cylinder,” Journal of Heat Transfer, vol. 131, no. 11, pp. 1–11, 2009.