Transient natural convection in a square cavity filled with a porous medium is studied numerically. The cavity is assumed heated from one vertical wall and cooled at the top, while the other walls are kept adiabatic. The governing equations are solved numerically by a finite difference method. The effects of Rayleigh number on the initial transient state up to the steady state are investigated for Rayleigh number ranging from 10 to 2×102. The evolutions of flow patterns and temperature distributions were presented for Rayleigh numbers, Ra=102 and 103. It is observed that the time taken to reach the steady state is longer for low Rayleigh number and shorter for high Rayleigh number.
1. Introduction
Convective heat transfer has attracted significant attention due to wide applications in engineering such as operation of solar collectors, cooling systems in electronics equipments, and insulations of buildings. Many studies with application to the previous research areas may be found in the books by Nield and Bejan [1], Ingham and Pop [2], and Vafai [3].
In recent years, most researchers devoted their studies on natural convection related to either a vertically or horizontally imposed heat flux or temperature difference; see, for example, Bejan [4], Rasoul and Prinos [5], Goyeau et al. [6], Saeid and Pop [7], Baytas and Pop [8], and Zeng et al. [9]. Works on natural convection with differentially heated neighbouring walls include those of Ganzarolli and Milanez [10], Ameziani et al. [11], Ishihara et al. [12], Rahman and Sharif [13], Aydin and Yang [14], Ece and Büyük [15], Frederick and Berbakow [16], and Dalal and Das [17].
The heat transfer characteristics of natural convection in a square and rectangular cavity heated from one side and cooled from the ceiling have been studied by Aydin [18] and Aydin et al. [19, 20]. Aydin and Ünal [21] conducted a numerical study on transient buoyant flow inside an air-filled 45°-inclined enclosure. Varol et al. [22] studied the natural convection in a porous enclosure divided by a triangular massive partition and heated and cooled on adjacent walls. In addition, Revnic et al. [23] investigated the natural convection in an inclined square cavity heated and cooled on adjacent walls and filled with a porous medium. Finally, we mention that Yesiloz and Aydin [24] studied experimentally and numerically the natural convection in an inclined quadrantal cavity heated from below and cooled on one side.
In this work, transient natural convection in a square cavity filled with a porous medium is studied numerically. The square enclosure is heated from one vertical side and cooled from the ceiling. The vertical wall is assumed heated to a uniform temperature or heat flux while the ceiling is cooled at a uniform temperature. The other vertical wall and the floor are adiabatic.
2. Mathematical Formulation
The schematic diagram of two-dimensional square cavity is shown in Figure 1. It is assumed that the left vertical wall is suddenly heated to the constant temperature or heat flux Th and the top wall is suddenly cooled to the constant temperature Tc where Th>Tc by equal amount relative to an initially uniform temperature distribution. The right vertical and bottom walls are adiabatic.
Schematic diagram of the physical model and coordinate system.
Applying Darcy’s flow model and the Boussinesq approximation, the governing equations are as follows [7]: ∂u∂x+∂v∂y=0,∂u∂y-∂v∂x=-gβKν∂T∂x,σ∂T∂t+u∂T∂x+v∂T∂y=α(∂2T∂x2+∂2T∂y2),
where u and v are the velocity components along the x- and y-axes, respectively, T is the fluid temperature, and the definitions of the other symbols are mentioned in the Nomenclature. Equations (2) are subject to the following boundary conditions:u(x,y,0)=v(x,y,0)=0,T(x,y,0)=T0,u(0,y,t)=v(0,y,t)=0,T(0,y,t)=Th(y),u(L,y,t)=v(L,y,t)=0,∂T(L,y,t)∂x=0,u(x,0,t)=v(x,0,t)=0,∂T(x,0,t)∂y=0,u(x,L,t)=v(x,L,t)=0,T(x,L,t)=Tc(x).
Using the stream functions ψ defined by u=∂ψ/∂y, v=-∂ψ/∂x, and nondimensional variables:X=xL,Y=yL,θ=T-T0Th-Tc,Ψ=ψα,τ=αtσL2,
where T0=(Th+Tc)/2, the governing equations (2) can be written in dimensionless forms:∂2Ψ∂X2+∂2Ψ∂Y2=-Ra∂θ∂X,∂θ∂τ+∂Ψ∂Y∂θ∂X-∂Ψ∂X∂θ∂Y=∂2θ∂X2+∂2θ∂Y2,
where Ra is the Rayleigh number defined asRa=gβK(Th¯-Tc)Lνα.
The boundary conditions (2.2) can be written asΨ(X,Y,0)=0,θ(X,Y,0)=0,Ψ(0,Y,τ)=0,θ(0,Y,τ)=0.5,Ψ(1,Y,τ)=0,∂θ(1,Y,τ)∂X=0,Ψ(X,0,τ)=0,∂θ(X,0,τ)∂Y=0,Ψ(X,1,τ)=0,θ(X,1,τ)=-0.5.
The local Nusselt numbers along the hot and cold wall are given, respectively, by [23]Nuh(Y)=-2∂θ∂X(0,Y),Nuc(X)=-2∂θ∂Y(X,1),
and the average Nusselt numbers from the heated wall Nuh(Y) and cold wall Nuc(X) are, respectively, defined as
Nuh¯=∫01Nuh(Y)dY,Nuc¯=∫01Nuc(X)dX.
3. Numerical Scheme
The coupled system of (2.4) and (2.5) subject to boundary conditions (2.7)–(2.11) was solved numerically using a finite difference method. The alternating direction implicit method was applied for discretizing the equations. The unknowns Ψ and θ were calculated until the following convergence criterion was fulfilled:∑i,j|ζi,jn+1-ζi,jn|∑i,j|ζi,jn+1|≤ϵ,
where ζ is either Ψ or θ, n represents the iteration number, and ϵ is the convergence criterium. In this study, the convergence criterion was set at ϵ=10-5. In order to validate the computation coding, accuracy test had been performed for Ra=102 and Ra=103 for several grid sizes. Due to lack of suitable results from literature, we only compared our results with those from Revnic et al.[23] as shown in Table 1. There has been found a reasonably good agreement between the results obtained using 101×101 and results using Richardson extrapolation and grid size 101×101 as reported by Revnic et al. [23]. We noticed that further refinement in grid size would increase the average Nusselt numbers. For Ra=102,Ra=103 when the grid refined from 101×101 to 141×141 the increase in Nusselt number was 3.9 percent and 1.6 percent, respectively. However, in view of maximum value of stream function, Ψmax, an increase from 101×101 to 141×141 nodal points would increase the Ψmax less than 0.3 percent for both Ra. Therefore, in a way to compromise the accuracy and cost in calculations, the grid size 101×101 was used in all the calculations.
Comparison of the average Nusselt number, Nu¯, with results from [23] for heated and cold walls for Ra=102.
Reference
Mesh size
Nuh¯
Nuc¯
Revnic et al. [23]
26×26
8.2636
8.2240
51×51
9.2162
9.2008
101×101
10.1319
10.1247
Richardson extrapolation
10.4371
10.4326
Present results
26×26
8.5553
8.2786
51×51
9.3885
9.2271
101×101
10.2297
10.1383
4. Results and Discussion
The main objective of the present work is to investigate the effects of Rayleigh number, Ra, on natural convection in an enclosure heated and cooled on adjacent walls and filled with a porous medium. The computations were carried out for Ra=10 to 2×103. At first, the left and top walls were cooled and the fluid was at a uniform temperature and motionless in the enclosure. After sudden heating on the left wall, the fluid near the hot wall began to rise and flow along the cooled top wall. Then, the flow reached the other side of the wall and moved downward to the bottom wall. Finally, after passing the adiabatic walls, the flow reached the hot wall and completed the flow cycle. Due to buoyancy force, the fluid moved from the left region to the right region of the cavity, yielding a clockwise circulation cell inside the enclosure.
Figure 2 shows the transient result of the flow and temperature field for Ra=102. Initially, at t=0, there was no motion in the enclosure. At t=0.001, a single clockwise rotating cell appeared close to the hot wall. At this moment, the convective motion was not fully developed. Thus the isotherms crowded and formed a boundary layer parallel to the isothermal walls, implying that the conduction was the dominant heat transfer mechanism. After a short time, the flow grew in magnitude to their steady-state condition. As time increased (t=0.1), temperature had been distributed from the left to the right wall; therefore the streamlines intensity increased and the center of cell moved to the center of the enclosure. Meanwhile the isotherms formed thermal boundary layer along the hot and cold wall. This indicated that convection occurred and dominated the heat transfer in the enclosure. As time progressed longer (t=0.8), no significant difference was observed in the flow and temperature fields, which indicated that the flow had reached the steady-state regime.
Time history for streamlines (left) and isotherms (right) for Ra=100.
t=0.001
t=0.01
t=0.05
t=0.1
t=0.5
t=0.8
The transient development of flow structure for Ra=103 is illustrated in Figure 3. At the beginning, the flow structure was almost similar to that of the case Ra=102. However, after a short time (t=0.05), the streamline had been extended throughout the cavity and convection became more important. In addition, for t>0.1, the flow was then going to attain the steady-state condition. The thermal boundary layer was more pronounced and grew faster than the lower Ra case.
Time history for streamlines (left) and isotherms (right) for Ra=1000.
t=0.001
t=0.005
t=0.01
t=0.05
t=0.1
t=0.2
The steady-state condition for the cavity with different values of Ra is shown in Figure 4. For Ra=10, it was observed that flow exhibited a unicellular recirculation pattern and the corresponding isotherm formed a diagonally symmetric structure indicating that the heat transfer mechanism was still under influence of conduction. As the convective motion increased with the increasing in Ra, the colder flow tended to occupy the lower part of the enclosure and the temperature gradient was more severe near the heated vertical wall and cooled ceiling. Therefore, a thermal boundary layers were observed near the heated and cooled walls.
Streamlines (top) and isotherms (bottom) at steady-state for difference Ra.
Ra=10
Ra=100
Ra=1000
Ra=2000
Figure 5 depicts the transient responses of the average Nusselt number, Nuh¯, on the hot wall for different values of the Rayleigh number. Initially, the value of Nuh¯ was large because of the high-temperature gradient close to the hot wall. For Ra=10, the Nuh¯ decreased monotonously with time until it converged to a steady-state value. This behavior was related to conduction as the dominant heat transfer mechanism for low Ra value. However, for Ra=102, Nuh¯ dropped to a minimum value first and then increased continuously to a steady-state value. Similar behavior could be observed for Ra>102, dropping to a minimum and increasing with some oscillations to a steady-state value. As we observed, the steady-state value was reached faster with the increasing in Ra. This was a consequence of increased convective motion with the increasing on Raleigh number.
Variations of the average Nusselt number on the hot wall with time, t=0 to t=1.0.
Figure 6 illustrates the variation of the local Nusselt number Nuh along the hot wall (left) and Nuc for cold wall (right) at the steady-state regime for different Ra. It was observed that, for all Ra considered, except for Ra=10, the Nuh decreased to a minimum value first, followed by a sharp increase near the top wall. For all values of Ra, the Nuc decreased continuously along the ceiling to the adiabatic right wall.
Variations of the local Nuh and local Nuc on the hot and cold wall, respectively, for different Ra at steady-state condition.
5. Conclusion
The present study considers a transient natural convection in a two-dimensional square cavity filled with a porous medium heated from the left vertical wall and cooled from the top wall, while the other walls are kept adiabatic. The effect of Rayleigh number on the transient thermal behavior was investigated. Initially, the heat transfer process was characterized by pure conduction. However, after a short time and for Rayleigh number Ra>10, the convection dominated the flow motion in the enclosure. Increasing the Rayleigh number could lead to more convective motion and shorter time to reach a steady-state condition.
Nomenclatureg:
Magnitude of gravitational acceleration
K:
Permeability of the porous medium
L:
Cavity height/width
Nu:
Nusselt number
Ra:
Rayleigh number for porous medium
t:
Time
T:
Fluid temperature
u,v:
Velocity components along x- and y-axes, respectively
U,V:
Nondimensional along X- and Y-axes, respectively
x,y:
Cartesian coordinates
X,Y:
nondimensional Cartesian coordinates.
Greek Symbolsα:
Effective thermal diffusivity
β:
Coefficient of thermal expansion
θ:
Nondimensional temperature
κ:
Wave number
ν:
Kinematic viscosity of the fluid
σ:
Ratio of composite material heat capacity to convective fluid heat capacity
τ:
Non-dimensional time
ψ:
Stream function
Ψ:
Nondimensional stream function.
Subscripts and Superscriptsc:
Cold wall
h:
Hot wall.
Acknowledgments
The first author gratefully acknowledges the financial support from the Ministry of Higher Education of Malaysia and Universiti Teknologi MARA.
NieldD. A.BejanA.20062ndNew York, NY, USASpringerZBL1189.76205InghamD. B.PopI.2005Oxford, UKElsevier ScienceZBL1189.80022VafaiK.20052ndBoca Raton, Fla, USATaylor & Francis10.1201/9780415876384ZBL1189.76121BejanA.On the boundary layer regime in a vertical enclosure filled with a porous medium197962931022-s2.0-001844349510.1016/0094-4548(79)90001-8RasoulJ.PrinosP.Natural convection in an inclined enclosure1997754384782-s2.0-003144981310.1108/09615539710187783ZBL1064.76613GoyeauB.SongbeJ. P.GobinD.Numerical study of double-diffusive natural convection in a porous cavity using the Darcy-Brinkman formulation1996397136313782-s2.0-003013150810.1016/0017-9310(95)00225-1ZBL0963.76585SaeidN. H.PopI.Transient free convection in a square cavity filled with a porous medium2004478-9191719242-s2.0-164252102210.1016/j.ijheatmasstransfer.2003.10.014ZBL1106.76457BaytasA. C.PopI.Free convection in a square porous cavity using a thermal nonequilibrium model20024198618702-s2.0-003633839710.1016/S1290-0729(02)01379-0ZengM.WangQ. W.HuangZ. P.WangG.OzoeH.Numerical investigation of natural convection in an enclosure filled with porous medium under magnetic field200752109599712-s2.0-3464881696610.1080/10407780701364346ZBL1113.80016GanzarolliM. M.MilanezL. F.Natural convection in rectangular enclosures heated from below and symmetrically cooled from the sides1995386106310732-s2.0-002927842510.1016/0017-9310(94)00217-JZBL0925.76695AmezianiD. E.BennacerR.BouhadefK.AzziA.Effect of the days scrolling on the natural convection in an open ended storage silo20094812225522632-s2.0-7034978191210.1016/j.ijthermalsci.2009.04.005IshiharaI.FukuiT.MatsumotoR.Natural convection in a vertical rectangular enclosure with symmetrically localized heating and cooling zones20022333663722-s2.0-003660266310.1016/S0142-727X(02)00184-4RahmanM.SharifM. A. R.Numerical study of laminar natural convection in inclined rectangular enclosures of various aspect ratios20034443553732-s2.0-024241513710.1080/713838233AydinO.YangW. J.Natural convection in enclosures with localized heating from below and symmetrical cooling from sides20001055185292-s2.0-003369742510.1108/09615530010338196ZBL1171.76444EceM. C.BüyükE.Natural-convection flow under a magnetic field in an inclined rectangular enclosure heated and cooled on adjacent walls20063885645902-s2.0-3374566845710.1016/j.fluiddyn.2006.04.002ZBL1178.76322FrederickR. L.BerbakowO.Natural convection in cubical enclosures with thermal sources on adjacent vertical walls200241333134010.1080/10407780252780199DalalA.DasM. K.Natural convection in a rectangular cavity heated from below and uniformly cooled from the top and both sides20064933013222-s2.0-3034448488210.1080/10407780500343749AydinO.Transient natural convection in rectangular enclosures heated from one side and cooled from above19992611351442-s2.0-003298206810.1016/S0735-1933(98)00129-8AydinO.ÜnalA.AyhanT.Natural convection in rectangular enclosures heated from one side and cooled from the ceiling19994213234523552-s2.0-034424199010.1016/S0017-9310(98)00319-6ZBL1042.76557AydinO.ÜnalA.AyhanT.Numerical study on buoyancy-driven flow in an inclined square enclosure heated and cooled on adjacent walls19993665855992-s2.0-034396225410.1080/104077899274589AydinO.ÜnalA.Transient buoyant convection in a 45°-inclined enclosure heated and cooled on adjacent walls20061644514692-s2.0-3374482869910.1108/09615530610653082VarolY.OztopH. F.PopI.Natural convection flow in porous enclosures with heating and cooling on adjacent walls and divided by a triangular massive partition20083544764912-s2.0-3914914185210.1016/j.icheatmasstransfer.2007.09.010RevnicC.PopI.GrosanT.InghamD. B.Natural convection in an inclined square cavity with heated and cooled adjacent walls and filled with a porous medium20091202910.2174/1877729500901010020YesilozG.AydinO.oaydin@ktu.edu.trNatural convection in an inclined quadrantal cavity heated and cooled on adjacent walls20113561169117610.1016/j.expthermflusci.2011.04.002