The Scientific World Journal

Volume 2014 (2014), Article ID 786102, 8 pages

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

Research Article

Conjugate Heat Transfer in Rayleigh-Bénard Convection in a Square Enclosure

Habibis Saleh¹ and Ishak Hashim1,2,3

¹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 2 January 2014; Accepted 3 March 2014; Published 28 May 2014

Academic Editors: K. Hooman and S. Nardini

Copyright © 2014 Habibis Saleh and Ishak Hashim. 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.

Conjugate natural convection-conduction heat transfer in a square enclosure with a finite wall thickness is studied numerically in the present paper. The governing parameters considered are the Rayleigh number , the wall-to-fluid thermal conductivity ratio , and the ratio of wall thickness to its height . The staggered grid arrangement together with MAC method was employed to solve the governing equations. It is found that the fluid flow and the heat transfer can be controlled by the thickness of the bottom wall, the thermal conductivity ratio, and the Rayleigh number.

1. Introduction

The beautiful hexagonal patterns can be obtained by a simple experiment first conducted by Henri Bénard, a French physicist, in 1900. Later, in 1916, Lord Rayleigh proposed his theory of a feedback coupling resting on buoyancy: a fluid particle hotter than its environment encounters ever colder fluid as it rises, which leads to the instability. The critical Rayleigh numbers, which signal the onset of natural convection in rectangular boxes heated from below and cooled from above, have been obtained theoretically by Davis [1] and Catton [2]. Corcione [3] proposed the heat transfer correlations in terms of the average Nusselt number with Rayleigh number and the aspect ratio of the rectangular enclosure heated from below. This convection problem has attracted a great deal of attention from researchers because of its presence in both nature and industrial applications. In nature, the convection cells formed from air rising above sunlight-warmed land or water are a major feature of all weather systems. Convection is also seen in the rising plume of hot air from fire, oceanic currents, and sea-wind formation. A very common industrial application of natural convection is free air cooling without the aid of fans: this can happen on small scales (computer chips) to large scale process equipment.

The literatures indicated that most of the reported works on Rayleigh-Bénard convection do not study the effect of a conductive bottom wall. Applications of conductive bottom wall can be found, for example, in high performance insulation for buildings. When the conductivities of the wall and fluid are comparable and the wall thickness is finite, conduction-convection analysis is necessary. This coupled conduction-convection problem is known as conjugate convection. Conjugate natural convection in a rectangular enclosure surrounded by walls was firstly examined by Kim and Viskanta [4, 5]. Their results show that wall conduction effects reduce the average temperature differences across the cavity, partially stabilize the flow, and decrease the heat transfer rate. A differentially heated vertical square enclosure with two finite thickness horizontal walls was investigated by Mobedi [6]. Kaminski and Prakash [7] and Misra and Sarkar [8] performed a numerical study on conjugate convection in a square cavity with a thick conducting wall on one of its vertical sides. The influence of wall conduction on natural convection in an inclined square cavity was researched by Acharya and Tsang [9]. The effects of the Rayleigh number, dimensionless conductivity ratio, dimensionless wall width, and inclination angle on the natural convection in an inclined enclosure bounded by a solid wall were investigated by Ben Yedder and Bilgen [10]. Later, Nouanegue et al. [11] extended the work of Ben Yedder and Bilgen [10] to include a radiation effect. Zhang et al. [12] studied conjugate convection in an enclosure with time-periodic sidewall temperature and inclination. Recently, Aminossadati and Ghasemi [13] investigated conjugate convection in an inclined nanofluid-filled enclosure.

The investigation of the effect of the conductive bottom wall on convective flows in a square enclosure has not received much attention. Recent works include those of Varol et al. [14] and Saleh et al. [15] for Darcy-Bénard convection in a porous enclosure. The aim of this work is to examine the effect of the conductive bottom wall on Rayleigh-Bénard convection in a square enclosure. This effect on the flow development, temperature distribution, and heat transfer rate in the wall and fluid will be presented graphically.

2. Mathematical Formulation

A schematic diagram of a square enclosure with a finite wall thickness is shown in Figure 1. The bottom surface of the impermeable wall is heated to a constant temperature , and the top surface of the enclosure is cooled to a constant temperature , while the vertical walls are kept adiabatic.

Figure 1: Schematic representation of the model.

Thermophysical properties of the fluid in the flow field are assumed to be constant except the density variations causing a body force term in the momentum equation. The Boussinesq approximation is invoked for the fluid properties to relate density changes to temperature changes and to couple in this way the temperature field to the flow field. Under the above assumptions, 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 impermeable wall is

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

the resulting nondimensional forms of (1)-(2) are the following:

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

where is the thermal conductivity ratio. The physical quantities of interest in this problem are the average Nusselt number defined by

where represents the dimensionless heat transfer through the walls.

3. Computational Methodology

Staggered grid arrangement together with the Marker and Cell (MAC) method [16] is adopted to solve the governing equations (4) subject to the boundary conditions (5). Due to the lack of boundary conditions for pressure, the use of the staggered grid and MAC formulation provides an advantage. That is, one may locate the secondary grid along the boundaries of the domain where only specification of velocity boundary conditions is required but not of the pressure. The fictitious values of velocity outside the domain are obtained by extrapolation of the interior points as given by Hoffmann and Chiang [17]. The temperature conditions at the interface boundary are

where the subscript is used to represent the coordinate and is the number of nodal points in the -axis in the wall.

Regular and uniform grid distribution is used for the whole enclosure. The effect of grid resolution was examined in order to select the appropriate grid density as demonstrated in Figure 2. The results indicate that 120 × 120 grid can be used in the final computations. To validate the computational code, the previously published problems on conjugate natural convection in a square cavity with a conducting side wall were solved. Table 1 shows the comparison of the average Nusselt number between the present results and the available results found in the literature for various Rayleigh numbers and thermal conductivity ratios. The comparison was in excellent agreement with the results reported by [7, 18, 19]. These comprehensive verification efforts demonstrated the robustness and accuracy of the present computation.

Table 1: Comparison of the average Nusselt number between the present results and other works for various Rayleigh numbers and wall-to-fluid thermal conductivity ratios.

Kr

1
5
10
1
5
10
1
5
10

Present
0.868
1.019
1.041
1.349
1.830
1.916
2.089
3.424
3.725

[7]
0.870
1.020
1.040
—
—
—
2.080
3.420
3.720

[18]
0.870
1.020
1.040
1.350
1.830
1.920
2.080
3.420
3.720

[19]
0.850
1.030
1.040
—
—
—
2.040
3.300
3.600

Figure 2: Grid independency study: versus number of grid points 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 wall thickness, , the thermal conductivity ratio, , and the Rayleigh number, .

Figure 3 illustrates the effects of the wall thickness parameter for and on the thermal fields and flow fields in the fluid and in the bottom solid wall. As can be seen, the parameter affects the fluid and the solid temperatures as well as the flow characteristics. The strength of the flow circulation of the fluid is much higher for a thin solid bottom wall. A circular main cell is formed at a thin bottom; then, the main cell shape becomes elliptical and finally breaks up into dual contrarotative cells at . This is because the fluid adjacent to the hotter wall has lower density than the fluid at the middle plane. As a result, the fluid moves upward due to the Archimedes force from the middle portion of the top wall. When the fluid reaches the upper part of the enclosure, it is cooled, so its density increases; then, the fluid flows downward at the left and the right planes of the enclosure. This creates a successive cell that is well known as Bénard cells. It is important to note that the Rayleigh number in the present study is based on the total height of the enclosure.

(a)

(b)

(c)

Figure 3: Isotherms (left) and streamlines (right) at and .

To show the effect of the thermal conductivity ratio on the thermal fields and the circulation of the fluid in the enclosure, the isotherms and streamlines are presented in Figure 4 for and . Three different conductivity ratios are selected: , , and . It is observed that two contrarotative cells are formed as shown in Figures 4(a)–4(c). The clockwise (“−” sign) circulation cell refers to natural circulation; that is, the main cell and the counterclockwise (“+” sign) circulation cell refer to the secondary cell. As the conductivity ratio increases, the magnitude of the main cell increases, while the magnitude of the secondary cells decreases and shrinks. This phenomenon is due to the temperature gradient near the wall that increases with the increase of the parameter . Thus, much heat transfer from the bottom solid wall to the fluid is obtained for higher values of (good conductive solid wall). It is also observed that convection effects in the fluid become stronger for higher values of .

(a)

(b)

(c)

Figure 4: Isotherms (left) and streamlines (right) at and .

Figures 5(a)–5(c) show the effects of on the thermal fields and flow fields in the porous enclosure and in the bottom solid wall with constant values of and . As can be seen in Figure 5(a), two contrarotative cells were formed with the same size and strength. When takes higher values as depicted in Figures 5(b) and 5(c), the main cell circulation strengthens while the secondary (left) one weakens and shrinks. Thermal fields show that the temperature distribution is almost uniform in the solid bottom wall for all values of investigated. The thermal fields in the fluid are modified strongly by increasing as shown in Figures 5(b) and 5(c). This refers to the strength of convection current related to the values.

(a)

(b)

(c)

Figure 5: Isotherms (left) and streamlines (right) at and .

Variations of the average Nusselt number with the Rayleigh number are shown in Figure 6(a) for different values of the wall thickness and thermal conductivity ratio, . The result presented in Figure 6(a) shows that, for a thin solid wall, the heat transfer from the fluid increases with increasing . This is due to the increasing of domination of convection heat transfer by increasing the buoyancy force inside the fluid. Figure 6(a) also shows that becomes constant for the highest values of the thickness parameter of the solid bottom wall. Variations of the average Nusselt number with the Rayleigh number are shown in Figure 6(b) for different values of with constant . Obviously, the heat transfer increases by increasing . The heat transfer enhancement by increasing is more pronounced at higher values of thermal conductivity ratio as shown in Figure 6(b). This is due to the temperature gradient near the solid wall that increases with increasing as shown in Figure 4. Variations of the average Nusselt number with the wall thickness are presented in Figure 7 for different values of with constant . This figure shows that the heat transfer decreases by increasing the solid wall thickness . There is a considerable difference between the heat transfer for small and large values of . For , the heat transfer is almost constant. When takes higher values, heat transfer drops sharply by increasing and becomes a conduction mode. This is due to the bottom solid wall that behaves as an insulated material in this case. Finally, the average Nusselt number can be correlated pretty well with the wall thickness, the thermal conductivity ratio, and the Rayleigh number as follows:

(a)

(b)

Figure 6: Variations of with for different (a) and (b) .

Figure 7: Variations of with for different .

This correlation is valid for , , and . Note that the maximum error is less than for the valid ranges.

5. Conclusions

The present numerical simulations study the effects of conduction in bottom wall on Rayleigh-Bénard convection in a square enclosure. The dimensionless forms of the governing equations were solved using the finite difference method. Detailed computational results for flow and temperature fields and the heat transfer rates 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 circulation of the fluid is much higher with thin walls and/or higher value of the solid to fluid thermal conductivity ratio.(2)The number of contrarotative cells and the strength of the circulation of each cell can be controlled by the thickness of the bottom wall, the thermal conductivity ratio, and the Rayleigh number.(3)The average Nusselt number increases by increasing either the Rayleigh number and/or thermal conductivity ratio, but the average Nusselt number decreases by increasing the wall thickness.

Nomenclature

: | Wall thickness and dimensionless wall thickness
: | Gravitational acceleration
: | Thermal conductivity ratio
: | Thermal conductivity
: | Width and height of enclosure
: | Average Nusselt number
: | Rayleigh number
: | Temperature
: | Velocity components in the - and -directions
: | Space coordinates and dimensionless space coordinates.

Greek Symbols

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

Subscript

: | Cold
: | Fluid
: | Hot
: | Wall.

Conflict of Interests

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

References

S. H. Davis, “Convection in a box: linear theory,” Journal of Fluid Mechanics, vol. 30, pp. 465–478, 1967.

I. Catton, “The effect of insulating vertical walls on the onset of motion in a fluid heated from below,” International Journal of Heat and Mass Transfer, vol. 15, no. 4, pp. 665–672, 1972.

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.

D. M. Kim and R. Viskanta, “Study of the effects of wall conductance on natural convection in differently oriented square cavities,” Journal of Fluid Mechanics, vol. 144, pp. 153–176, 1984.

D. M. Kim and R. Viskanta, “Effect of wall heat conduction on natural convection heat transfer in a square enclosure,” Journal of Heat Transfer, vol. 107, no. 1, pp. 139–146, 1985.

M. Mobedi, “Conjugate natural convection in a square cavity with finite thickness horizontal walls,” International Communications in Heat and Mass Transfer, vol. 35, no. 4, pp. 503–513, 2008.

D. A. Kaminski and C. Prakash, “Conjugate natural convection in a square enclosure: effect of conduction in one of the vertical walls,” International Journal of Heat and Mass Transfer, vol. 29, no. 12, pp. 1979–1988, 1986.

D. Misra and A. Sarkar, “Finite element analysis of conjugate natural convection in a square enclosure with a conducting vertical wall,” Computer Methods in Applied Mechanics and Engineering, vol. 141, no. 3-4, pp. 205–219, 1997.

S. Acharya and C. H. Tsang, “Influence of wall conduction on natural convection in an inclined square enclosure,” Heat and Mass Transfer, vol. 21, no. 1, pp. 19–30, 1987.

R. Ben Yedder and E. Bilgen, “Laminar natural convection in inclined enclosures bounded by a solid wall,” Heat and Mass Transfer, vol. 32, no. 6, pp. 455–462, 1997.

H. F. Nouanegue, A. Muftuoglu, and E. Bilgen, “Heat transfer by natural convection, conduction and radiation in an inclined square enclosure bounded with a solid wall,” International Journal of Thermal Sciences, vol. 48, no. 5, pp. 871–880, 2009.

W. Zhang, C. Zhang, and G. Xi, “Conjugate conduction-natural convection in an enclosure with time-periodic sidewall temperature and inclination,” International Journal of Heat and Fluid Flow, vol. 32, no. 1, pp. 52–64, 2011.

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.

Y. Varol, H. F. Oztop, and I. Pop, “Conjugate heat transfer in porous triangular enclosures with thick bottom wall,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 19, no. 5, pp. 650–664, 2009.

H. Saleh, N. H. Saeid, I. Hashim, and Z. Mustafa, “Effect of conduction in bottom wall on Darcy-Bénard convection in a porous enclosure,” Transport in Porous Media, vol. 88, no. 3, pp. 357–368, 2011.

F. H. Harlow and J. E. Welch, “Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface,” Physics of Fluids, vol. 8, no. 12, pp. 2182–2189, 1965.

K. A. Hoffmann and S. T. Chiang, Computational Fluid Dynamics. Volume I, Engineering Education System, Wichita, Kan, USA, 2000.

M. Hriberšek and G. Kuhn, “Conjugate heat transfer by boundary-domain integral method,” Engineering Analysis with Boundary Elements, vol. 24, no. 4, pp. 297–305, 2000.

N. Wansophark, A. Malatip, and P. Dechaumphai, “Streamline upwind finite element method for conjugate heat transfer problems,” Acta Mechanica Sinica, vol. 21, no. 5, pp. 436–443, 2005.