** Journal of Electronics Cooling and Thermal Control ** Vol. 3 No. 1 (2013) , Article ID: 29390 , 15 pages DOI:10.4236/jectc.2013.31002

Effect of Sinusoidal Heating on Natural Convection Coupled to Thermal Radiation in a Square Cavity Subjected to Cross Temperature Gradients

^{1}Laboratory of Flows and Transfers Modelling (LAMET), Physics Department, Faculty of Sciences and Technics, Sultan Moulay Slimane University, Béni-Mellal, Morocco

^{2}Laboratory of Fluid Mechanics and Energetics (LMFE), Unit Affiliated to CNRST (URAC 27), Physics Department, Faculty of Sciences Semlalia, Cadi Ayyad University, Marrakech, Morocco

Email: hasnaoui@ucam.ac.ma

Received November 30, 2012; revised January 3, 2013; accepted January 19, 2013

**Keywords:** Natural Convection; Thermal Radiation; Heatlines; Cross Gradients of Temperature; Periodic Heating; Resonant Heat Transfer; Numerical Study

ABSTRACT

Coupled natural convection and surface radiation within a square cavity, filled with air and submitted to discrete heating and cooling from all its walls, is studied numerically. The thermally active elements are centrally located on the walls of the cavity. Two heating modes, called SB and SV, are considered. They correspond to bottom and vertical left elements sinusoidally heated in time, respectively, while the top and vertical right ones are constantly cooled. The remaining portions of all the walls are considered adiabatic. The parameters governing the problem are the amplitude and the period of the temporally sinusoidal temperature, the emissivity of the walls, the relative lengths of the active elements and the Rayleigh number. The effect of such parameters on flow and thermal fields and the resulting heat transfer is examined. It is shown that, during a flow cycle, the flow structure can present complex behavior, depending on the emissivity and the amplitude and period of the exciting temperature. The rate of heat transfer is generally enhanced in the case of sinusoidal heating. Also, the resonance phenomenon existence, characterized by maximum fluctuations in flow intensity and heat transfer, is proved in this study.

1. Introduction

Natural convection coupled with surface radiation in closed cavities has been extensively studied using numerical simulations and experiments, owing to the practical importance of such configurations in many engineering applications (convective heat losses from solar collectors, thermal design of buildings, air conditioning and, recently, electronic cooling). The majority of the existing studies, which are of numerical nature, concerned with rectangular cavities where the temperature gradient is either horizontal or vertical, including different kinds of boundary conditions [1-6]. Results of these studies show that radiation affects the dynamical and thermal structures of the fluid, reduces natural convection heat transfer component, and contributes to increase the total amount of heat exchanged in the configurations considered. Most of the works conducted in the past on natural convection coupled with radiation inside rectangular enclosures have been substantially oriented to study unidirectional heat transfers resulting from imposed temperature gradients (due to heat fluxes or temperature differences) either parallel or normal to gravity. In some practical situations, much more complex boundary conditions may be encountered where horizontal and vertical temperature gradients could be simultaneously imposed across the cavity. This justifies the presence of published works where the rectangular cavities are heated from below and cooled from above and simultaneously submitted to various specified thermal boundary conditions at the sidewalls. In the absence of radiation effect, Corcione [7] investigated steady natural convection of airfilled rectangular enclosures heated from below, cooled from above and submitted to various thermal boundary conditions at the sidewalls. With respect to the basic thermal configuration wherein the sidewalls are insulated, the heat transfer rate resulting from the hot bottom wall increases when each adiabatic sidewall is replaced with a cooled one, while showing only slight decreases as such replacement is carried out through a heated sidewall. Cianfrini et al. [8] studied natural convection heat and momentum transfer in air-filled tilted square enclosures with differentially heated opposite walls. The results obtained for the average Nusselt number of the whole cavity were expressed through a semi-empirical dimensionless correlation. More recently, Deng [9] studied laminar natural convection in a two-dimensional square enclosure having two and three discrete heat source—sink pairs on the vertical sidewalls. The effects of the size and the arrangement of the sources and sinks on the fluid flow and heat transfer characteristics were investigated. El Ayachi et al. [10] analyzed the interaction of radiation and natural convection in a square cavity filled with air and discretely heated and cooled from the four walls. In this study, the bottom active wall is at higher temperature than that of the top active wall. As results of their study, compared to pure natural convection, heat transfer across the cavity was found to be enhanced notably with the increase of the active wall emissivity, while the effect of the insulated wall emissivity was insignificant.

In the above-mentioned studies, thermal boundary conditions were assumed to be either steady isothermal or constant heat flux wall conditions. However, in many engineering applications, the energy provided to the system is variable in time and gives rise to unsteady natural-convection flow. Solar collectors and printed circuit boards are examples of such systems submitted to variable thermal boundary conditions. In addition, thermal and dynamical behaviors of a fluid subjected to time-dependent thermal conditions are impossible to predict on the basis of the results obtained with constant temperature or heat flux conditions. This justifies the presence of some works in the literature in which the variable aspect of the thermal boundary conditions was considered. In an earlier study, Lage and Bejan [11] studied theoretically the problem of natural convection in enclosures with one side heated with a pulsating heat flux. They showed that the buoyancy-induced flow resonates to a certain frequency of the pulsating heat input and the resonance phenomenon is characterized by maximum fluctuations observed in the heat transfer evolution with the period of the time-dependent temperature. The numerical results obtained by Lakhal et al. in the case of a square cavity totally [12] or partially [13] heated from below with periodic variable temperatures showed that the resulting flow structure and heat transfer were strongly dependent on the amplitude and the period of the variable temperature. Later, Lakhal et al. [14] studied transient natural convection in a square cavity isothermally cooled from above and partially heated from the side. The heating temperature was varied sinusoidally in time and in a pulsating manner. The pulsating temperature leads to better enhancement of the heat transfer in comparison with that engendered by the sinusoidal varying temperature. Antohe and Lage [15] investigated theoretically and numerically the transport of momentum and heat while considering clear fluid and fully saturated porous medium differentially heated enclosures with a time-periodic pulsating heat flux. The numerical simulations indicated that the natural-convection activity within the enclosure reaches several local maxima for certain values of the heating frequency, referred to as resonance frequency. Kwak and Hyun [16] and Kwak et al. [17] also studied such a resonant phenomenon of natural convection. It was clearly noted from their studies that the resonant frequency might be well predicted by the Brunt-Väisälä frequency, which is characterized by the thermal stratification of a system. The case of convection heat transfer was considered by Abourida et al. for a square cavity filled with air and submitted to variable thermal boundary conditions on its horizontal [18] and vertical [19] walls. Different modes of heating were considered. At high Rayleigh numbers, results reported in [18] showed that the periodic heating can be used to enhance slightly or to reduce notably the heat losses in comparison with the case of constant temperatures. For a cavity with vertical active walls [19], it was reported that the buoyancy-driven flow presents a tendency to enter in resonance with the periodic heating only when the cold temperature is maintained constant. The resonance phenomenon was not encountered when both hot and cold temperatures were varied simultaneously. Douamna et al. [20] considered two-dimensional transient natural convection in a horizontal channel periodically heated from below with a temperature varying sinusoidally with time. Results of the study showed that the three different routes leading to chaos are identified by progressively varying the amplitude and the period of the variable heating temperature. In a recent work, Zhao et al. [21] analyzed transient laminar natural convection induced by two discrete heating elements flush-mounted on one vertical wall of a square enclosure. The resonance frequencies were obtained numerically and predicted by theoretical analysis. More recently, El Ayachi et al. [22] investigated natural convection coupled with thermal radiation in a square cavity, differentially heated with a temperature varying periodically in time. The existence of the resonance phenolmenon was observed in this study and found to be more intensified by increasing the emissivity and the amplitude of the exciting temperature. The critical period of the exciting temperature leading to the resonance phenomenon was found to be independent on the amplitude the excitation and the emissivity of the walls.

The importance of natural convection coupled with thermal radiation in a cavity subjected to periodic boundary conditions in time, is justified by the relevance of such a transitional process to many technological applications. The power supply of electronic circuits by an alternating current, the collectors of solar energy, rooms housing and building hollow blocks, in which recirculation is periodically driven by daily solar heating, are concrete examples.

In our knowledge, works dealing with time periodic combined natural convection-radiation in rectangular cavities subjected to crossed thermal gradients are nonexistent. This work is, therefore, a numerical contribution to the study of natural convection and surface radiation within a square cavity filled with air and discretely heated and cooled, in a periodic manner, from the four walls. The main parameters governing the problem are the emis-sivity of the walls and the characteristics of the exciting temperature (amplitude and period). The effect of these parameters on heat transfer and fluid flow within the cavity is examined.

2. Problem Formulation

The schematic view of the geometry considered in the present study is given in Figure 1. It consists of a square cavity discretely heated and cooled from the four walls. Two heating modes are considered. In the first one, named SB heating mode, the bottom active element is heated with a temperature varying sinusoidally in time, while the vertical left active element is maintained hot at a constant temperature. In the second mode, called SV heating mode, the vertical hot element temperature varies sinusoidally in time while that of the bottom active wall is considered constant. The inner surfaces, in contact with the fluid, are assumed to be gray, diffuse emitters and reflectors of radiation. The flow is conceived to be laminar, two-dimensional and incompressible with negligible viscous dissipation. All the thermophysical properties of the fluid are assumed constant except the density in the buoyancy term which is assumed to vary linearly with temperature (Boussinesq approximation); such a variation gives rise to the buoyancy forces.

Taking into account the above-mentioned assumptions, the non-dimensional governing equations, written in vorticity-stream function formulation, are as follows:

(1)

(2)

Figure 1. Studied configuration and imposed thermal excitations.

(3)

The dimensionless stream function and the vorticity are related to the non-dimensional velocity components by the following expressions:

(4)

2.1. Boundary Conditions

The dimensionless boundary conditions, associated to the problem are such:

on the cavity walls (5a)

on the cooled elements (5b)

SB Mode

(5c)

SV Mode

(5d)

(5e)

“n” being the normal direction to the considered adiabatic wall.

2.2. Coupling between Radiation and Natural Convection

The working fluid (air) is considered to be perfectly transparent to radiation. Thus air does not participate directly to the radiative heat transfer process (the energy equation remains unchanged) but indirectly via the solid surfaces that are implicated in the radiation exchange. The coupling of the thermal model is performed by computing the radiative exchanges. The calculation of the radiative heat exchange between the internal walls of the cavity is based on the radiosity method. The grid used for convection is retained in the presence of radiation and consists of 324 isothermal elementary segments. Each segment is sufficiently short to be assumed isothermal. The view factors between the isothermal elementary surfaces were determined by the Hottel’s [23] crossed string method. The summation rules were checked to be sure that the view factors summation equals unity for each surface. The inner surfaces of the enclosure are assumed to be opaque, diffuse and gray. The dimensional radiosity equation for the i^{th} element of the enclosure, in the case of a radiatively non-participating medium is:

is the outgoing radiation energy flux of the surface S_{i}. It consists of emission, E_{i}, plus the reflection of incoming energy,.

The non-dimensional radiosity equation may be written as:

(6)

The non-dimensional net radiative heat flux leaving a surface S_{i} is evaluated by:

2.3. Heat Transfer

(7)

At each time step, the mean Nusselt numbers, characterizing the contributions of natural convection and thermal radiation through the heated walls, are evaluated as:

- on the vertical heated wall

(8)

(9)

- on the horizontal heated wall

(11)

(12)

The instantaneous convective, radiative and total Nusselt numbers across the whole cavity are defined respectively as:

(12)

The mean Nusselt numbers, averaged in time over periods are calculated as:

(13)

(14)

where t_{cv} and t_{rd} are respectively the periods of the temporal variations of convective and radiative Nusselt numbers (they are identical in general).

2.4. Heatlines Formulation

Visualization of the paths followed by the heat flows through the enclosure requires the use of the heatlines concept (lines of constant heat function H). Such lines are defined through the first derivatives of the function H as follows [24]:

(15)

The dimensionless heat function equation can be derived easily from Equation (15) as:

(16)

Equation (16) is a conduction type problem with a source term becoming zero if the fluid flow subsides (fluid at rest). The solution of Equation (16) yields the values of the heat function in the nodes of the computational domain. The contour plots of the heat function values provide heatline patterns. The H field is defined through its first-order derivatives and only differences of H values are required instead of its intrinsic values. Thus, we have the freedom to state that (an arbitrary reference value for the heat function H). The boundary conditions for the dimensionless heat function, Equation (16), are obtained by integrating Equation (15) along the considered boundary. The following equations can be used to determine the values of heat function at the walls of the cavity:

(17)

3. Method of Solution

The non linear partial differential governing equations, Equations (1)-(3), were discretized using a finite difference technique. First and second spacial derivatives were approximated by means of central differences and the time derivatives by a forward difference. The integration of Equations (1) and (2) was ensured by the Alternating Direction Implicit method (ADI). At each time step, the Poisson equation, Equation (3), was treated by using the Point Successive Over-Relaxation method (PSOR) which with an optimum over-relaxation coefficient equal to 1.92 for the grid adopted in the present study. The set of Equations (6), representing the radiative heat transfer between the different elementary surfaces of the cavity, was solved by using the Gauss-Seidel method. The numerical code was validated against the results of Akiyama and Chong [2] obtained in the case of a square cavity differentially heated. Results of comparisons, made in terms of convective Nusselt numbers, evaluated at the heated wall and presented in Table 1, show a fairly good agreement with relative maximum deviations limited to for for Ra varying in the range. Results of a qualitative comparison in terms of isotherms, streamlines and heatlines, against those obtained by Deng [9] are presented in Figure 2 for Ra = 10^{5} and Pr = 0.71 in the case of a square cavity discretely heated and cooled on the four walls without radiation effect; a good qualitative agreement can be seen by comparing these figures. Details relative to the sensitivity of the results vis-à-vis of the grid are given in Table 2. The results presented in this table justify the selected grid of 81 ´ 81 as a reasonable compromise between computational effort and required accuracy. More precisely, the maximum deviations observed in terms of and remain within 0.61%, 0.90%, 1.28% and 0.11%, respectively, when the grids of 61 ´ 61 and 101 ´ 101 are considered.

4. Results and Discussion

The main parameters governing the problem are the amplitude of the exciting temperature, it’s period, the emissivity of the walls

, the Prandtl number, Pr, the Rayleigh number, Ra, and the relative length of active elements, B. To highlight the influence of a, t and e, the values of B, Pr and Ra are fixed to 0.5, 0.72 (air) and 10^{6} respectively. In the following, the results obtained are presented in terms of temporal variations of maximum and minimum stream functions and total Nusselt number Nu (Nu being the sum of the convective and radiative Nusselt numbers). The mean values of these quantities, averaged in time over one flow cycle are also discussed to point out the effect of the variable temperature on the coupling between natural convection and surface radiation for both cases of heating considered (SB and SV modes).

Table 1. Effect of Ra and e on, evaluated on the heating wall of a square cavity for and.

Table 2. Grid sensitivity analysis for Ra = 10^{6} and various values of e.

(a)(b)

Figure 2. Numerical code validation in terms of heatlines in the case of a differentially heated square cavity for Ra = 10^{5} and Pr = 0.7: (a) Deng [9] and (b) Our code.

4.1. Temporal Variation of Ψ_{min}, Ψ_{max} and Nu

For SB and SV heating modes, the evolution with time of Ψ_{min}, Ψ_{max}, and Nu, are presented in Figures 3-5 for (absence of radiation), various values of a, (for SB mode) and (for SV mode). The choice of these two values of t is guided by a series of tests showing significant effects on fluid flow and heat transfer, which are discussed later. The reference case corresponding to a constant heating temperature (a = 0) is represented in these figures by a continuous horizontal line. Figures 3 and 4 show that, globally the evolution in time of Ψ_{min} and Ψ_{max} is periodic, the amplitudes of the oscillations are visibly enhanced when the amplitude a of the exciting temperature is increased. It can be seen from Figure 3 that, in the case of the SB heating mode, the main flow intensity is higher in comparison to the case of SV heating mode. It can also be noted that the monocellular flow structure, observed for a = 0, disappears in favor of a secondary flow whose size and intensity increase with a (Figure 4). In terms of heat transfer, the increase of the amplitude a is accompanied by a proportional increase in the amplitude of the sinusoidal evolution of Nu (Figure 5).

The effect of radiation on the temporal variation of Ψ_{min}, Ψ_{max} and Nu is shown in Figures 6-8 in the case of highly emissive walls.

The amplitude of Ψ_{min} oscillations is clearly reduced for both heating configuration considered (Figure 6),

Figure 3. Effect of the amplitude a on for ε = 0.

while that of Ψ_{max} is slightly affected, negatively in the case of the SB heating mode and positively in the SV one (Figure 7). It should be noted that for the two heating modes, SB and SV, the clockwise cell, constituting the

Figure 4. Effect of the amplitude a on Ψ_{min}(t) for ε = 0.

Figure 5. Effect of the amplitude a on Nu(t) for ε = 0.

main flow, remains dominant during the evolution of the cycle, while the secondary counterclockwise flow, plays a less important role. For the heat transfer, Figure 8, shows that when radiation of the walls is considered, the amplitude of Nu oscillations, is significantly enhanced; indicating an improvement of the total heat transfer when

Figure 6. Effect of the amplitude a on Ψ_{min}(t) for ε = 1.

Figure 7. Effect of the amplitude a on Ψ_{min}(t) for ε = 1.

radiation is taken into account.

4.2. Streamlines, Isotherms and Heat-Lines over a Flow Cycle

From the temporal evolution of Ψ_{max}, discussed above, it

Figure 8. Effect of the amplitude a on Nu(t) for ε = 1.

appears that the secondary flow intensity is less important in the case of SV mode than in the SB one. That is why in this subparagraph, the study will be limited to this later configuration. Figures 9 and 10 correspond to the case of constant thermal excitation and are used as reference, respectively for e = 0 and e = 1.

On the other hand, in Figures 11 and 12 are presented the streamlines, isotherms and heatlines, obtained during a flow cycle in the case of sinusoidal exciting temperature, respectively for e = 0 (non-radiative walls) and e = 1 (strongly radiative walls). Also, in Figures 11 and 12, the letters (a), (b), (c), (d) and (e) correspond to the selected times in Figures 4 and 7. As can be seen from Figure 9, the flow is dominated mainly by a large cell zone surrounding three small others located in the central part of the cavity. Once the modulated temperature is applied (a = 1), a fundamental change in the flow structure is observed (Figure 11). Thus, the streamlines show that during the evolution of the cycle, the flow passes through various states characterized by the dominance, at the instant (a), of large clockwise dominant cell and the appearance, at instants (b), (c) and (d), of a secondary counterclockwise cell, small and of

(a)(b)(c)

Figure 9. (a) Streamlines; (b) Isotherms and (c) Heatlines in the case of constant heating (a = 0) for ε = 0.

(a)(b)(c)

Figure 10. (a) Streamlines; (b) Isotherms and (c) Heatlines in the case of constant heating (a = 0) for ε = 1.

(a)(b)(c)(d)(e)

Figure 11. Streamlines (left), isotherms (center) and heatlines (right) over a flow cycle for the SB mode in the case where e = 0. (a) Y_{min} = −75.63, Y_{max} = 0.124; (b) Y_{min} = −64.52, Y_{max} = 8.57; (c) Y_{min} = −59.40, Y_{max} = 17.49; (d) Y_{min} = −56.9 , Y_{max} = 8.48; (e) Y_{min} = −71.75, Y_{max} = 0.03.

(a)(b)(c)(d)(e)

Figure 12. Streamlines (left), isotherms (center) and heatlines (right) over a flow cycle for the SB mode in the case where ε = 1. (a) Ψ_{min} = −79.4, Ψ_{max} = 2.2; (b) Ψ_{min} = −67.63, Ψ_{max} = 6.46; (c) Ψ_{min} = −64.53, Ψ_{max} = 12.83; (d) Ψ_{min} = −65.41, Ψ_{max} = 6.43; (e) Ψ_{min} = −80.05, Ψ_{max} = 1.83.

low intensity in the lower left corner of the cavity. The corresponding isotherms are tightened at the active elements, indicating that the convective heat exchange is good in this region of the cavity. Lines corresponding to heat flux (heatlines) show a similarity with the streamlines and demonstrate the dominance of convective heat transfer (high values of the stream function).

Furthermore, for e = 1 and a = 0, Figure 10 shows a more complex flow structure consisting of a large and intense main cell zone and four small secondary counterclockwise cells taking place in the four corners of the cavity. By introducing a sinusoidal thermal excitation, a specific examination of the streamlines of Figure 12 reveals that the flow structure, during one cycle of evolution, is dominated by a large and intense clockwise cell, surrounded, at the four corners of the cavity, by secondary counterclockwise cells more or less important according to the considered instants. Because of this complexity, the isotherms show a consistent temperature distribution, generally characterized by large thermal gradients in the vicinity of the active elements. Meanwhile, the heatlines show a similar shape to that of the streamlines, reflecting the dominance of convective heat transfer in such a situation.

4.3. Effect of the Period τ on and

Based on a set of preliminary tests, the time-averaged values of the stream functions and the Nusselt number (Y_{max}, Y_{min}, and Nu), obtained for various periods, let suppose the existence of the resonance phenomenon within the cavity. This phenomenon is characterized by maximum fluctuations in the response of the system to an external excitation with well-defined period. This effect was reported in the past by several authors, both in the absence of radiation effect [11-21] and in the presence of radiation [22]. Hence, variations of and are presented versus t for both considered heating modes and various values of e. The values corresponding to the case of the isothermal heating sources are also presented in the same figures as references.

In the absence of radiation, variations of and are presented in Figures 13-15, respectively.

Variations of the main flow intensity, , with t (Figure 13) show that increases with τ to a peak, characterizing a resonance phenomenon, for a critical value of τ, which is of 0008 for the SB mode and 0.00825 for the SV one. Such a difference between these critical periods can be attributed to the difference between the flow intensities for the two heating modes. Higher intensities of the main flow would lead to lower values of the

Figure 13. Variations of with the period τ for ε = 0 and different values of a.

Figure 14. Variations of with the period τ for ε = 0 and different values of a.

resonant periods. Above this threshold, a further increase of τ leads to a continuous weakening of the flow intensity toward asymptotic values greater than those corresponding to the constant heating temperature, reached at high values of τ. By increasing the amplitude a of the variable temperature, the peak becomes more important and the resonance phenomenon becomes more pronounced,

Figure 15. Variations of with the period τ for ε = 0 and different values of a.

but without changing the critical periods.

It was reported in previous studies that the resonance occurs when the time duration, necessary for the flow to rotate through the full cavity (period of a fluid wheel [11]), is equal to that of the cyclic modulated temperature. The rotational dimensional period required for a quantity of fluid to rotate along the vertical and horizontal walls in an enclosure can be expressed as [11]:

where and are the length and height of the enclosure, respectively. The scale of the characteristic velocity is determined by Lage and Bejan [11].

For the secondary flow, Figure 14 also shows a resonance phenomenon in the evolution of which is well marked and premature for the SB mode with almost the same critical period observed in the evolution of, while for the SV mode, the occurrence of this resonance is delayed and depends on the amplitude a. Beyond the critical period, decreases and stabilizes at values depending on the heating mode considered.

Regarding variations of with t (Figure 15), the resonance phenomenon is also observed and appears to be more important in the case of SB mode, with critical periods identical to those on (0.008 for the SB mode and 0.00825 for the SV one). The later values are not affected by the increase of the oscillating amplitude a, a result in accordance with previous numerical findings [22]. Also, it is important to note that beyond the critical value of t, decreases to a minimum value reached for t_{0} which depends on the heating mode imposed. In the case of SB mode, , while in the case of the SV one, t_{0} is substantially less and is 0.01, values beyond which increases until an asymptotic trend, depending on a.

In the case of highly emissive walls, Figures 16-18 show that the fluid circulation and the total heat transfer are more enhanced in comparison with the case of e = 0. When e passes from 0 to 1, the enhancement of

Figure 16. Variations of with the period τ for ε = 1 and different values of a.

Figure 17. Variations of with the period τ for ε = 1 and different values of a.

the time averaged heat transfer, , is of about 142% at the resonance.

The evolutions with τ are similar to those obtained without radiation effect. The resonance phenomenon is obtained for the same critical periods for both considered heating modes.

The contribution of radiation to total heat transfer

Figure 18. Variations of with the period τ for ε = 1 and different values of a.

through the cavity is also examined in this study for the two considered configurations, SB and SV heating modes. The obtained results show that for a moderate value of the emissivity, the contribution of radiation is not negligible and it is close to that obtained for the steady state case (a = 0). Note that in the case of the SV mode, an increase of a leads to a decrease in the contribution of radiation. In the SB heating mode, this decrease occurs specifically for high values of t. For e = 1, a higher contribution of radiation to total heat transfer is obtained (between 59% and 62%) when compared to that of natural convection.

5. Conclusion

The problem of natural convection coupled with thermal radiation inside a square cavity, submitted to cross gradients of temperature, is studied numerically. Two different heating modes (SB and SV), with a periodic variation in time, are considered. In such a situation, the solutions obtained with and without thermal radiation, are periodic in time with periods identical to that of the exciting temperature. A resonance phenomenon, characterized by maximum fluctuations of the flow intensity and heat transfer rate, is observed. The resonance period of SB mode is slightly lower than that of SV one. Both these periods have the particularity of being independent on the excitation amplitude and the walls emissivity. Except for low values of the excitation period, time variable heating improves, generally, heat transfer than constant heating, which makes the former interesting for the heat removal through the cavity. Moreover, the control of the coupling convection-radiation, in the unsteady regime, appears to be a good means to control heat transfer for very specific values of the governing parameters.

REFERENCES

- C. Balaji and S. P. Venkateshan, “Interaction of Surface Radiation with Free Convection in a Square Cavity,” International Journal of Heat Fluid Flow, Vol. 14, No. 3, 1993, pp. 260-267. doi:10.1016/0142-727X(93)90057-T
- M. Akiyama and Q. P. Chong, “Numerical Analysis of Natural Convection with Surface Radiation in a Square Cavity,” Numerical Heat Transfer Part A, Vol. 31, No. 4, 1997, pp. 419-433. doi:10.1080/10407789708913899
- G. Colomer, M. Costa, R. Consul and A. Oliva, “Three Dimensional Numerical Simulation of Convection and Radiation in a Differentially Heated Cavity Using the Discrete Ordinates Method,” International Journal of Heat Mass Transfer, Vol. 47, No. 2, 2004, pp. 257-269. doi:10.1016/S0017-9310(03)00387-9
- C. Y. Han and S. W. Baek, “The Effects of Radiation on Natural Convection in a Rectangular Enclosure Divided by Two Partitions,” Numerical Heat Transfer Part A, Vol. 37, No. 3, 2000, pp. 249-270. doi:10.1080/104077800274280
- E. H. Ridouane, M. Hasnaoui, A. Amahmid and A. Raji, “Interaction between Natural Convection and Radiation in a Square Cavity Heated from Below,” Numerical Heat Transfer Part A, Vol. 45, No. 3, 2004, pp. 289-311. doi:10.1080/10407780490250373
- R. El Ayachi, A. Raji, M. Hasnaoui and A. Bahlaoui, “Combined Effect of Radiation and Natural Convection in a Square Cavity Differentially Heated with a Periodic Temperature,” Numerical Heat Transfer Part A, Vol. 53, No. 12, 2008, pp. 1339-1356. doi:10.1080/10407780801960043
- 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, 2003, pp. 199-208. doi:10.1016/S1290-0729(02)00019-4
- C. Cianfrini, M. Corcione and P. P. Dell’Omo, “Natural Convection in Tilted Square Cavities with Differentially Heated Opposite Walls,” International Journal of Thermal Sciences, Vol. 44, No. 5, 2005, pp. 441-451. doi:10.1016/j.ijthermalsci.2004.11.007
- Q.-H. Deng, “Fluid Flow and Heat Transfer Characteristics of Natural Convection in Square Cavities Due to Discrete Source—Sink Pairs,” International Journal of Heat Mass Transfer, Vol. 51, No. 25-26, 2008, pp. 5949-5957. doi:10.1016/j.ijheatmasstransfer.2008.04.062
- R. El Ayachi, A. Raji, M. Hasnaoui, A. Abdelbaki and M. Naïmi, “Combined Effects of Radiation and Natural Convection in a Square Cavity Submitted to Cross Gradients of Temperature: Case of Partial Heating and Cooling,” Computational Thermal Sciences, Vol. 3, No. 1, 2011, pp. 73-78. doi:10.1615/ComputThermalScien.v3.i1.60
- J. L. Lage and A. Bejan, “The Resonance of Natural Convection in a Horizontal Enclosure Heated Periodically from the Side,” International Journal of Heat Mass Transfer, Vol. 36, No. 8, 1993, pp. 2027-2038. doi:10.1016/S0017-9310(05)80134-6
- E. K. Lakhal, M. Hasnaoui, P. Vasseur and E. Bilgen, “Convection Naturelle Dans une Cellule Carrée Chauffée Périodiquement par le Bas: Étude Numérique,” Revue Generale de Thermique, Vol. 392, No. 392-393, 1994, pp. 480-485.
- E. K. Lakhal, M. Hasnaoui, P. Vasseur and E. Bilgen, “Natural Convection in a Square Enclosure Heated Periodically from Part of the Bottom Wall,” Numerical Heat Transfer Part A, Vol. 27, No. 3, 1995, pp. 319-333. doi:10.1080/10407789508913703
- E. K. Lakhal, M. Hasnaoui and P. Vasseur, “Numerical Study of Transient Natural Convection in a Cavity Heated Periodically with Different Types of Excitations,” International Journal of Heat Mass Transfer, Vol. 42, No. 21, 1999, pp. 3927-3941. doi:10.1016/S0017-9310(99)00070-8
- B. V. Antohe and J. L. Lage, “Amplitude Effect on Convection Induced by Time-Periodic Horizontal Heating,” International Journal of Heat Mass Transfer, Vol. 39, No. 6, 1996, pp. 1121-1133. doi:10.1016/0017-9310(95)00207-3
- H. S. Kwak and J. M. Hyun, “Natural Convection in an Enclosure Having a Vertical Sidewall with Time Varying Temperature,” Journal of Fluid Mechanics, Vol. 329, 1996, pp. 65-88. doi:10.1017/S0022112096008828
- H. S. Kwak, K. Kuwahara and J. M. Hyun, “Resonant Enhancement of Natural Convection Heat Transfer in a Square Enclosure,” International Journal of Heat Mass Transfer, Vol. 41, No. 18, 1998, pp. 2837-2846. doi:10.1016/S0017-9310(98)00018-0
- B. Abourida, M. Hasnaoui and S. Douamna, “Transient Natural Convection in a Square Enclosure with Horizontal Walls Submitted to Periodic Temperatures,” Numerical Heat Transfer Part A, Vol. 36, No. 7, 1999, pp. 737-750. doi:10.1080/104077899274543
- B. Abourida, M. Hasnaoui and S. Douamna, “Convection Naturelle Dans une Cavité Carrée Avec des Parois Verticales Soumises à des Températures Périodiques,” Revue Generale de Thermique, Vol. 37, No. 9, 1998, pp. 788- 800.
- S. Douamna, M. Hasnaoui and B. Abourida, “Two-Dimensional Transient Natural Convection in a Repetitive Geometry Submitted to Variable Heating from Below: Numerical Identification of Routes Leading to Chaos,” Numerical Heat Transfer Part A, Vol. 37, No. 7, 2000, pp. 779-799. doi:10.1080/104077800274073
- F. Y. Zhao, D. Liu and G. F. Tang, “Resonant Response of Fluid Flow Subjected to Discrete Heating Elements,” Energy Conversion Management, Vol. 48, No. 9, 2007, pp. 2461-2472. doi:10.1016/j.enconman.2007.04.008
- R. El Ayachi, A. Raji, M. Hasnaoui and A. Bahlaoui, “Combined Effect of Radiation and Natural Convection in a Square Cavity Differentially Heated with a Periodic Temperature,” Numerical Heat Transfer Part A, Vol. 53, No. 12, 2008, pp. 1339-1356. doi:10.1080/10407780801960043
- H. C. Hottel and A. F. Saroffim, “Radiative Heat Transfer,” McGraw-Hill, New York, 1967.
- S. Kimura and A. Bejan, “The Heatline Visualization of Convective Heat Transfer,” ASME Journal of Heat Transfer, Vol. 105, No. 4, 1983, pp. 916-919. doi:10.1115/1.3245684

Nomenclature

A_{r}: aspect ratio of the cavity

B: relative lengths of the active portions,

cv: convection F_{ij}: view factor from S_{i} surface to S_{j} one_{}

g: acceleration due to gravity,

: length of the vertical active portions, m

: height of the cavity, m J_{i}: dimensionless radiosity,

: width of the horizontal active elements, m

: length of the cavity, m N_{r}: convection-radiation interaction parameter,

Nu: average Nusselt number Pr: Prandtl number,

Q_{r}: dimensionless radiative heat flux,

Ra: Rayleigh number,

rd: radiation t: dimensionless time,

: dimensional fluid temperature, K T: dimensionless fluid temperature,

: temperature of the cold elements, K

: temperature of the heated elements, K

: temperature difference,

T_{r}: dimensionless reference temperature,

: dimensionless horizontal and vertical velocities,

: dimensionless coordinates,

Greek Symbols

: thermal diffusivity of fluid,

: thermal expansion coefficient of fluid,

: emissivity of the walls

: thermal conductivity of fluid,

n: kinematic viscosity of fluid,

t: dimensionless period of the heating temperature,

: dimensionless vorticity,

: dimensionless stream function,

: Stéfan-Boltzman constant,

Subscripts

C: cooled surface H: heated surface

Superscripts

¢: dimensional variable