Heat Mass Transfer DOI 10.1007/s00231-015-1583-6

ORIGINAL

Modelling of the mixed convection in a lid‑driven cavity with a constant heat flux boundary condition Przemysław Błasiak1 · Piotr Kolasin´ski1 

Received: 15 September 2014 / Accepted: 2 May 2015 © Springer-Verlag Berlin Heidelberg 2015

Abstract  In this paper steady state two-dimensional mixed convection heat transfer problem in a lid-driven cavity heated via an uniformly distributed heat flux on the bottom wall is investigated numerically. The lid moves with constant velocity and is kept at low constant temperature, between two ideally thermally insulated vertical walls. A wide range of Prandtl Pr and Richardson Ri number is examined to study their effects on heat transfer rate and fluid flow. Governing parameters are 0.001 ≤ Ri ≤ 1.0 and 0.71 ≤ Pr ≤ 56.00. Grashof number Gr is fixed at 104. The results are presented in the form of isotherms and streamlines plots. Also, local and mean Nusselt number are depicted on charts. Numerical values of the surface averaged Nusselt number are also presented. Results show that increase of Prandtl number strongly influences enhancement of heat transfer rate and that decreasing of Richardson number increases surface averaged Nusselt number. Mechanisms responsible for intensification of heat transfer are identified and physical explanation of this phenomenon are also given. List of symbols A Surface (m2) g Gravitational acceleration (m/s2) h Heat transfer coefficient (W/m2 K) k Thermal conductivity (W/m K) L Cavity side length (m) n Iteration loop counter

* Przemysław Błasiak [email protected] 1



Department of Thermodynamics, Theory of Machines and Thermal Systems, Wrocław University of Technology, Wybrzez˙e Wyspian´skiego 27, 50‑370 Wrocław, Poland

Nu Nusselt number Pr Prandtl number r Rate of grid stretching Re Reynolds number Ri Richardson number s Spatial location of the node T Temperature (K) T Average temperature of the wall (K) U, V Dimensionless velocities in x and y direction q′′ Heat flux (W/m2) x, y Cartesian coordinates (m) X, Y Dimensionless Cartesian coordinates Greek symbols α Thermal diffusivity (m2/s) β Thermal expansion coefficient (1/K) Δ Step size Θ Dimensionless temperature ν Kinematic viscosity (m2/s) Subscripts i, j Nodes numbers in the mesh T Parameters on the top (cold) wall B Parameters on the bottom (hot) wall

1 Introduction Many engineering problems, where mixed convection flow and heat transfer occurs, can be predicted with the use of a simple model known as lid-driven cavity flow. This model can be used for the description of heat transfer phenomenon in the following applications: electronic devices cooling, geothermal systems, lubrication technologies, chemical processes, solar energy systems, etc. Detailed application

13



of this model was presented in [1]. Additionally heat and fluid flow in a square cavity are a standard test case for code validation. This is well-known in the field of applied mathematics and numerical methods [2, 3]. Mixed convection heat transfer and fluid flow in a square cavity were studied by many scientists and it is well described in the literature. Modelling of this problem can be classified into two categories, based on applied boundary conditions, influencing buoyancy forces. The first group is related to Dirichlet boundary conditions i.e. constant temperature on the surface of the wall. The second group is related to Neumann boundary conditions i.e. constant heat flux on the surface of the wall. Most studies described in the literature are dedicated to the first group [4–9]. However this study is related to modelling heat and fluid flow in a square cavity, in which, constant heat flux boundary condition is employed on the bottom horizontal wall. Studies related to mixed convection heat and fluid flow where the Neumann boundary conditions are applied, are rarely described in the literature. Mixed convection heat transfer in a cavity, in the presence of thermal sources, were studied in series of papers by Papanicolaou and Jaluria [10–13]. Authors examined numerically mixed convection heat transfer from an isolated thermal source with uniform surface heat flux located in a rectangular cavity with a fluid outlet [10]. Analysed Richardson and Reynolds number were in the range 0–10 and 50–2000 respectively. The working fluid was air. Authors found that average Nusselt number increased with increasing Richardson number at fixed Reynolds number. More intensive increase in heat transfer rate with increasing Richardson number was also observed. Additionally influence of different configurations resulting from various locations of heat source and outlet was investigated. In [11] authors presented numerical study on transition from steady, laminar regime to a periodic regime in a cavity filled with air subjected to heating via a heat source. The heating source was located on one of the vertical walls. Influence of Grashof and Reynolds numbers on resulting heat transfer and fluid flow were investigated. It was found that flow became unstable beyond a critical value of Richardson number. Small fluid elements at the centre of the cavity were cooled and heated periodically. Consequently this resulted in circulation of bulk mass of the fluid which was caused by density change. In [12] authors numerically studied conjugate heat transfer in the case of an isolated, constant source of heat flux, located on one of the walls. The configuration of the cavity included an external flow which flowed into cavity through an inlet located on the left vertical wall and flowed out via outlet located on the right vertical wall. Influence of the following parameters on a heat transfer in the cavity were tested: Reynolds number, Richardson number, ratio of the wall material thermal conductivity to fluid thermal conductivity

13

Heat Mass Transfer

(air) and location of the heat source. Results showed that with increasing heat flux applied on heat sources surfaces, heat transfer rate from these surfaces also increased. Moreover it was observed that circulation of the fluid in the cavity was stronger and temperatures of heat sources were higher. Different locations of the heat source strongly influenced heat transfer rate. Authors concluded that heat source location led to buoyancy-assisting mixed convection flows, resulted in lower temperatures on the surfaces and higher heat transfer rates. In [13] authors performed a numerical study of combined forced and natural convective cooling of heat-dissipating electronic components. Those components were located in a rectangular cavity and cooled by an external through-flow of air. Richardson number was varied from 0.01 to 10. Influence of heat sources number and location on heat transfer rate and fluid flow in the cavity was analysed. As it was shown location of the heat source on the right vertical wall is the most beneficial case in terms of cooling. Hsu and Wang [14] carried out a numerical study on mixed convection in a cavity. Discrete heat sources were embedded on a vertical board located on the bottom wall of the cavity. Air flow was directed into the cavity through an inlet located on the left vertical wall and flowed out through an outlet located on the opposite wall. Hsu and Wang studied influence of the same parameters as Papanicolaou and Jaluria [12] on heat transfer and fluid flow. It was shown that at fixed Reynolds number decrease of Richardson number resulted in decrease of average Nusselt number. Highest temperature gradients were located near heat sources. Also it was observed that with increasing Reynolds number heat transfer rate also increased. Guo and Sharif [15] numerically analysed mixed convection heat transfer in a rectangular cavity with constant heat flux located on the centre part of the bottom wall. The side walls were isothermal and moved in the vertical direction. Heat source length was changed in a range 0.2–0.8 L. Computations were done for Richardson number range 0.1–10. Authors studied impact of Richardson number, heat source length and location, and aspect ratio of the cavity on maximum temperature and Nusselt number. It was found that low Ri values influenced temperature in such a way that it was more uniformly distributed within the cavity. Also it was found relatively large region of the cavity is affected by the heat source. It was also noticed that the region affected by the heat source became larger with increasing heat source length. Ogut [16] examined mixed convection in an inclined lid-driven cavity filled with air. Two vertical walls were adiabatic, while the bottom wall was heatedup via constant heat flux. The moving lid was isothermal. Richardson number was in a range 0.1–100 and inclination angle varied from 0° to 180°. Author found that inclination angle did not have significant effect on heat transfer and

Heat Mass Transfer

fluid flow for Ri = 0.1. Richardson number significantly influenced heat transfer rate and with increasing Richardson number the average Nusselt number increased as well. Shahi et al. [17] numerically investigated mixed convection nano-fluid flows in a square cavity with inlet and outlet ports. Constant heat source was symmetrically located at the bottom wall. Richardson number was changed in the range of 0 ≤ Ri ≤ 10. It was found that increase in solid volume fraction, resulted in increase in average Nusselt number at the heat source surface and in decrease of average bulk temperature. Mansour et al. [18] examined mixed convection in a square cavity partially heated from below and filled with water-base nano-fluid. Authors considered impact of Reynolds number, solid volume fraction, different heat source lengths and locations on local and average Nusselt number along the heat source. The heat source length was varied by 0.2 ≤ L ≤ 0.8. It was found that for high Reynolds numbers (Re = 100) the lid-driven effect was dominant. For small values of Reynolds number (Re = 10) an increase in effect of heat source presence on heat and fluid flow was observed. Cho et al. [19] numerically investigated heat transfer in water-based nano-fluids confined within a lid-driven cavity with wavy left and right walls. On the left wall constant heat flux was applied and the right wall was maintained at a lower temperature. Both the upper and lower walls were thermally insulated and moved in the horizontal direction. Results showed that for studied range of Richardson number, mean Nusselt number increased, when volume fraction of nano-particles and Grashof number increased. Above described literature review shows that studies dedicated to mixed convection heat transfer in a lid-driven cavity with Neuman boundary conditions are limited. Additionally none of the researchers, except Ogut [16], did not examine the case, where heat flux was applied on the whole surface of the wall. However the main aim of his work was to analyse effect of an inclination angle on heat transfer rate. Thus, Prandtl number was kept constant. In present study mixed convection heat transfer in liddriven cavity is investigated. Fluid in the cavity is heated through the surface of the bottom wall via a constant heat flux and cooled via the moving lid, which is kept at constant, low temperature. Two vertical walls are thermally insulated. This simple configuration was not considered in the literature. Nevertheless, in our opinion, it is an important problem since the constant heat flux boundary condition is more physically realistic. Mathematical treatment of a heat transfer with use of constant temperature boundary condition is easier. However, most engineering issues do not fulfil isothermal conditions at boundaries. Modelling of heat transfer phenomena using constant heat flux boundary condition better reflects actual situation. As previously shown, literature on this topic is rare, and thorough study

is needed. Furthermore all of the papers were aimed at calculating local and mean Nusselt number along the surface where the heat flux was employed. In such situations the method of Nusselt number calculation is well known from e.g. [15, 16, 18]. Thus, no results for Nusselt number on other surfaces are available in the literature. In the present study the main goal is to investigate variation of local Nusselt number on the moving lid which is maintained at a constant temperature. In order to calculate local and mean Nusselt number at this surface, approach similar to this presented in [20] is used. In this work we decided to focus on following problems: the first is to examine influence of a wide range of Prandtl number and Richardson number on heat transfer and fluid flow in the cavity. The second objective is to explain mechanisms of scraping of the boundary layer phenomenon in the cavity. The third goal is to provide a benchmark solution for this simple configuration with use very fine resolution grid (1024 × 1024).

2 The problem description The scheme of the model and the boundary conditions are depicted in Fig. 1 in form of a two-dimensional square cavity of length side L. The top lid surface moves with constant velocity UT from left side of the cavity to the right, while the other walls are stationary. The top lid surface is maintained at constant temperature TT. A working fluid is heated through the bottom wall via a constant heat flux q′′. Vertical walls are adiabatic. Acceleration of gravity acts in the negative vertical direction. In order to examine influence of Prandtl number on heat transfer, following working fluids were chosen for

Fig. 1  Numerical domain and the boundary conditions

13



Heat Mass Transfer

simulation: air, water and industrial aniline. Prandtl numbers of these fluids are 0.71, 9.54 and 56.00 respectively. Grashof number was fixed at 104.

3 Mathematical model Two-dimensional, steady state, laminar fluid flow and heat transfer in a square cavity are modelled. Fluid is treated as Newtonian and incompressible. Thermophysical properties at reference temperature TT are assumed to be constant, except density in the buoyancy term which follows the Boussinesq approximation [21]. It is also assumed that radiation heat transfer is negligible compared to convective and conductive mechanisms of heat exchange in the fluid. Additionally viscous dissipation is also neglected. The governing equations for the problem are the continuity, momentum and energy equation. After applying these assumptions, these equations in non-dimensional form can be written as follows

∂V ∂U + =0 ∂X ∂Y

(1)

  ∂U ∂P 1 ∂ 2U ∂ 2U ∂U +V =− + + U ∂X ∂Y ∂X Re ∂X 2 ∂Y 2

(2)

U

  ∂V ∂P 1 ∂ 2V ∂ 2V Gr ∂V +V =− + + + 2 θ (3) 2 2 ∂X ∂Y ∂Y Re ∂X ∂Y Re

∂Θ 1 ∂ 2Θ ∂ 2Θ ∂Θ +V = + U ∂X ∂Y PrRe ∂X 2 ∂Y 2 



y u v p T − TT x ,Y = ,U = ,V = ,P = ,Θ = 2 L L UT UT �T ρUT

Re =

gβL 3 �T Gr ν UT L , Gr = , Pr = , Ri = 2 2 ν ν α Re

(4)

where T¯ B is an average temperature at the bottom surface wall (AB) and is defined as:



AB

TB dA AB

(10)

4 Numerical details

(5)

(6)

All distances are normalized by L, all velocities are normalized by UT, and pressure is normalized by ρU2T, where ρ is fluid density. Temperature is normalized as (T − TT)/ΔT, where ΔT is a temperature difference defined as follows: q′′ L T = k (7)

13

For each wall of the cavity the no-slip and the no-permeable wall boundary conditions were applied. Once velocity and temperature distributions in the cavity are obtained, local and area-averaged Nusselt numbers are calculated as follows:   ∂T  ∂y Y =1 L hL 1 L (9) =− NudX , Nu = Nu = k L 0 TB − T T

T¯ B =

where U and V are dimensionless components of the velocity vector in the X and Y direction respectively, P is dimensionless pressure and Θ is dimensionless temperature. Non-dimensional parameters are Reynolds Re, Grashof Gr and Prandtl Pr numbers. Non-dimensional variables and non-dimensional parameters are defined as follows: X=

where q′′ is a heat flux applied on the bottom cavity wall. Thermophysical properties of the fluid appearing in Eq. 6 are kinematic viscosity (ν), thermal expansion coefficient (β), thermal conductivity (k) and thermal diffusivity (α). Non-dimensional ratio Gr/Re2 in Eq. 3 is Richardson number. It represents the magnitude of natural convection relative to the forced convection. Ri number is commonly used parameter for assessment the convection flow regime. If Ri ≫ 1 than natural convection regime occurs, whereas forced convection is dominant for Ri ≪ 1. In the range of Ri = 0.1–1.0 both natural and forced convection play significant role in heat transfer and neither can be neglected. Boundary conditions are depicted in Fig. 1 and are specified as follows:  Y = 0, 0 < X < 1, U = 0, V = 0 ∂Θ/∂Y = −1    Y = 1, 0 < X < 1, U = 1, V = 0 Θ = 0 (8) X = 0, 0 < Y < 1, U = 0, V = 0 ∂Θ/∂X = 0    X = 1, 0 < Y < 1, U = 0, V = 0 ∂Θ/∂X = 0

The set of governing equations was discretized and solved with the use of ANSYS-CFX software [22]. The discretization method is an element-based finite volume approach [23] using non-staggered grid layout such that the control volumes are the same for all transport equations. In order to avoid a checkerboard pressure field, the Rhie-Chow algorithm is used [24]. ANSYS-CFX is a fully implicit coupled solver, which solves continuity equation Eq. 1 and momentum equations Eqs. 2 and 3 as a single system. Thus the pressure–velocity coupling is inherent in the solution procedure (pressure correction equation is not needed as opposed to the segregated solvers). In the simulations very fine meshes with resolution of 1024 × 1024 are used in order to accurately resolve boundary layers near the cavity walls. Two mesh-laws for nodes distribution are applied: uniform and geometric. For the geometric mesh type the arrangement of the nodes in the mesh is given by [25]:

Heat Mass Transfer

(11)

si+1 = si + r

where si is the spatial location of the node, Δ is the step size and r is the rate of grid stretching. The uniform mesh is used when the working fluid is air. In other cases geometric mesh is used. Use of different meshes was caused by convergence problems. For fluids with higher Prandtl numbers (compared to air), it was impossible to obtain convergent solution on the uniform mesh. Additionally to facilitate convergence of the solution each simulation was calculated twice, first on the coarse mesh and then on the 1024 × 1024 mesh. Solution from the coarse mesh was mapped on the dense mesh and then used as initial conditions. Convergence criteria used for each set of parameters were similar as those adopted in [4]  n+1  Ψ (i, j) − Ψ n (i, j)   < 10−6 MaxΨ n+1 (i, j)

   NuT  < 0.05 and 1.0 − NuB 

(12)

where i and j denote coordinates of the nodes in the numerical mesh, n is the iteration loop counter, Ψ is a dependent variable (u, v or T) and NuT and NuB are area-averaged

Nusselt numbers on the top and bottom walls respectively. Calculations were performed in a parallel mode on Supernova supercomputer of Wroclaw Centre for Networking and Supercomputing, being a part of Wroclaw University of Technology. Computational time was lengthened with increasing Prandtl number and decreasing Richardson number. For example to fulfil convergence criteria Eq. 12 in the case of Ri = 1 and Pr = 0.71 calculations lasted few hours. In the case of Ri = 0.001 and Pr = 56.0 calculations lasted 10 days. To validate described earlier study, two different testcases were considered. The first was standard flow without heat transfer in a lid-driven cavity of length side L. Calculations were carried out on the uniform mesh of resolution 1024 × 1024 for Reynolds number Re  = 1000 and compared with the very accurate results presented in [2]. Excellent agreement was obtained and the comparison of plots for streamlines, vorticity and pressure are presented in Fig. 2. In Tables 1 and 2 numerical values of non-dimensional velocity, pressure and vorticity are presented. The second test-case was the same as this treated by Ogut [16]. Validation was performed on uniform mesh

Fig. 2  Comparison of the steady state results for Re = 1000 (mesh size 1024 × 1024). From left to right stream function, vorticity and pressure. First row results presented in [2], second row present results

13

Table 1  Horizontal nondimensional velocity, pressure and vorticity through the vertical centerline of the cavity at Re = 1000 on grid 1024 × 1024

Table 2  Vertical nondimensional velocity, pressure and vorticity through the vertical centerline of the cavity at Re = 1000 on grid 1024 × 1024

Heat Mass Transfer Y

U [2]

U [present]

P [2]

P [present]

Ω [2]

Ω [present]

1.0000 0.9688 0.9531 0.7344 0.5000 0.2813 0.1016 0.0625

−1.00000 −0.58031 −0.47239 −0.18861 0.06205 0.28040 0.30029 0.20227

−1.00000 −0.58072 −0.47224 −0.18862 0.06205 0.28033 0.30035 0.20223

0.052971 0.051493 0.050314 0.012113 0.000000 0.040381 0.104416 0.10916

0.052899 0.051484 0.050356 0.0121126 0.000000 0.040410 0.104426 0.10911

0.0000

0.00000

0.00000

0.11056

0.11056

14.792 9.4781 4.8628 2.0909 2.0669 2.2678 −1.6352 −2.3174

−4.1554

14.884 9.4941 4.8573 2.0908 2.0668 2.2673 −1.6344 −2.3172

−4.1248

X

V [2]

V [present]

P [2]

P [present]

Ω [2]

Ω [present]

0.0000 0.0391 0.0547

0.00000 − 0.29330 − 0.41018

0.00000 −0.29355 −0.41023

0.077429 0.078658 0.077128

0.077432 0.078687 0.077111

−0.42639 0.02578 0.33391 0.33298 0.29622

0.049004 0.000000 0.047259 0.084369 0.087625

0.049120 0.000000 0.047256 0.084321 0.087622

−5.4967 − 8.2462 − 6.5097

− 5.6046 − 8.2436 − 6.5077

0.00000

0.00000

0.090448

0.090477

−7.6333

−7.5486

0.1406 0.5000 0.7734 0.9062 0.9297 1.0000

− 0.42634 0.02580 0.33398 0.33290 0.29622

3.4294 2.0669 2.0010 −0.8251 −1.5025

3.4292 2.0668 2.0013 −0.8237 −1.5026

Fig. 3  Comparison of the steady state results for Ri = 0.01 (Re = 1000; mesh size 1024 × 1024). Left column stream function, right column isotherms of non-dimensional temperature Θ. First row results presented in [16], second row present results

13

Heat Mass Transfer Fig. 4  Distribution of streamlines (left) and isotherms (right) in the cavity for Ri = 1.0 (Re = 100) and different Prandtl numbers. From top to the bottom Pr = 0.71, 9.47 and 56.00 respectively

Table 3  The set of parameters used in the simulations Fluid Air Water Aniline

Pr

Gr

Ri

0.71 9.47

104 104

0.001–1.0 0.001–1.0

56.00

104

0.001–1.0

with resolution 1024 × 1024 for Richardson number Ri  = 0.1 and Prandtl number Pr  = 0.71. Comparison of plots for streamlines and isotherms is depicted in Fig.  3. Additionally difference between area-averaged Nusselt number from this study and Ogut results was <2 %.

13



Heat Mass Transfer

Fig. 5  Distribution of streamlines (left) and isotherms (right) in the cavity for Ri = 0.1 (Re = 316) and different Prandtl numbers. From top to the bottom Pr = 0.71, 9.47 and 56.00 respectively

5 Results and discussion In this study the fluid flow and heat transfer phenomena are investigated for a wide range of Richardson numbers and Prandtl numbers. Sixteen numerical simulations were performed. Prandtl numbers were Pr = 0.71, 9.47 and 56.00 respectively. Ri values range was 0.001–1.0. Grashoff number was fixed at 104. Correspondingly the range of Reynolds numbers was from Re = 100 to Re = 3162. The set of parameters used during simulations is presented in Table 3.

13

The calculated streamlines and isotherms fields for Ri  = 1.0 (Re = 100) and for three different values of Pr are shown in Fig. 4. Movement of the fluid is driven by the shear induced by the moving lid. This results in fluid circulation in the clockwise direction. In the lower right cavity corner secondary counter-clockwise vortex appears. It can be seen that Pr does not have an influence on the pattern of the flow field. On the other hand Pr strongly affects temperature distribution. It is clear that with increase of Pr number temperature gradient also increases. The temperature field

Heat Mass Transfer Fig. 6  Distribution of streamlines (left) and isotherms (right) in the cavity for Ri = 0.01 (Re = 1000) and different Prandtl numbers. From top to the bottom Pr = 0.71, 9.47 and 56.00 respectively

is formed due to combined effects of natural convection caused by heat flux on the bottom wall and forced convection caused by the driving lid moving in the right direction. Cold fluid particles localized near the driving lid are moved to the right vertical wall due to shear forces. Then they are heated on the bottom wall. At the same time hot fluid particles localized at the bottom cavity part move towards left vertical wall. This is caused by buoyancy forces. After reaching the top of cavity area fluid particles are cooled by the low temperature fixed at the lid. These two phenomena

result in circulating fluid flow. Due to fluid movement isotherms also rotate in the clockwise direction. In the limit for the highest values of Pr, isotherms are moved to the left and to the bottom walls. Thus temperature boundary layer thickness decreases. It results in increase of temperature gradient near the walls. Consequently for fluids with large Pr, higher values of Nu are expected (Table 4). Figure 5 shows distribution of streamlines and isotherms in the cavity for Ri = 0.1 (Re = 316) and for three different values of Pr. Due to higher velocity of the moving lid

13



Heat Mass Transfer

Fig. 7  Distribution of streamlines (left) and isotherms (right) in the cavity for Ri = 0.001 (Re = 3162) and different Prandtl numbers. From top to the bottom Pr = 0.71, 9.47 and 56.00 respectively

circulation increases, and secondary counter-clockwise vortex in the right cavity corner enlarges. Moreover the core of the primary vortex is shifted slightly to the left and to the bottom. Curvature of the streamlines is more deformed near the left top cavity corner compared to previous case (with Ri = 1.0). This is primarily due to shearing effect of the driving lid. With increasing Pr one can see that it does not influence the pattern of the flow field. Regarding temperature field the situation is contrary. Circulation of isotherms is more visible compared to Ri = 1.0. In the case

13

of Pr = 0.71, in the centre and in the right cavity part the zone with uniform fluid temperature expands. An increasing Pr causes that zone expands and for Pr = 56.00 it occupies about 85 % of the cavity. Thickness of temperature boundary layer decreases with increasing Pr and layer is shifted to the left and bottom wall. The highest temperature gradients exist in thin regions on the left and bottom wall where the thinnest temperature boundary layer is localized. In the lower right corner, where the secondary counter-clockwise vortex appears, temperature boundary layer is thicker.

Heat Mass Transfer

Strong vorticity field in this corner results in increasing temperature gradient, but only in the thin layer. Overall thickness of the boundary layer in this zone is large compared to the middle zone of bottom wall and it constitutes thermal resistance. However with increasing Pr isotherms are moved towards the corner and temperature gradients increase. Streamlines and isotherms for Ri = 0.01 and for three different Pr values are presented in Fig. 6. In this case Re  = 1000 and fluid circulation is strong. The third vortex occurs in the bottom left cavity corner. This vortex also rotates counter-clockwise direction, same as secondary vortex in the bottom right corner. The vortex-occupied area in the right corner slightly extends. The primary vortex core displacement is significant. Its position is close to the cavity centre, this results from increasing inertia forces. For each Pr value the flow field is almost the same, this shows that it does not affect streamlines. In contrast Pr influences on heat transfer. It is evident that with increasing Pr the zone with uniform temperature distribution also increases. The

temperature boundary layer becomes thinner. This results in enhancement of heat transfer (due to higher temperature gradients, which exist predominantly on the left and bottom wall and in areas occupied by secondary vortices). For Pr = 56.00 one can see that secondary vortices contribute to thicken the temperature boundary layer. Consequently the width of the temperature boundary layer in these places is almost constant, despite temperature gradient increases. Fully developed vortex flow in the cavity for Ri = 0.001 (Re  = 3162) is depicted in Fig. 7. In this case degeneration of the streamlines in the left top corner (caused by the shearing effect), is so strong that the fourth counter-clockwise vortex appears. Circulation of the fluid is very strong. This results in increase of the zones occupied by the left and right bottom secondary vortex. Additionally the centre of the primary vortex core almost coincides with the cavity centre. Similarly to previous cases (with lower values of Ri) also in this case Pr does not influence streamlines in the cavity. The resulting streamlines caused by shearing forces

Fig. 8  Distribution of the Nusselt number along the lid for Prandtl number Pr = 0.71 and different Richardson numbers. From left to right and from top to the bottom Ri = 1, 0.1, 0.01 and 0.001 respectively

13



driven by the moving lid and viscosity of the fluid. On the contrary impact of Pr on temperature distribution is apparent, in terms of isotherms distribution. Even for low Prandtl number (Pr = 0.71), zone with uniform temperature field is very large and it covers about 75 % of the cavity. With increasing Pr the thickness of temperature boundary layer is very thin. However in the bottom corners, where secondary vortices occur, the temperature boundary layer width remains unchanged. Nevertheless temperature gradients increase in these areas with increasing Pr. The opposite effect occurs on the top wall, near the left top corner. Highest temperature gradients exist here due to the thinnest temperature boundary layer. In the vicinity of the left top corner highest values for local Nu are expected. This is confirmed by simulation results, treated in the following. Distributions of the local Nusselt number along the moving wall for Pr = 0.71 and four different values of Ri are depicted in Fig. 8. The expectation that the highest values of Nu will be located in the left top cavity corner is clearly confirmed for each value of Ri. The maximum Nu occurs

Heat Mass Transfer

for X = 0 and decreases with increasing X. This phenomenon can be explained as follows. Let it be assumed that the lid is stationary and infinitely long, and vertical walls are blades which move from right to left direction (opposite to the direction of lid motion in the actual case). This manipulation is fully allowable and does not influence the solution, due to governing equations for actual and imaginary problem are identical. Movement of the blades causes scraping of the temperature boundary layer. Just behind the blade temperature boundary layer is very thin and grows along the lid. Thinner temperature boundary layer results in smaller thermal resistance. Theoretically with increasing Re (decreasing Ri) in the left top cavity corner thickness of the temperature boundary layer tends to zero and Nu tends to infinity. This can be observed in the presented Nu plots (Figs. 8, 9, 10). Effect of temperature boundary layer scraping by the walls (blades) is larger with decreasing Ri i.e. increasing Re. However developing of the temperature boundary layer is fast and Nu decreases abruptly only in a range of X = 0–0.1. With decrease of Ri local Nu curve

Fig. 9  Distribution of the Nusselt number along the lid for Prandtl number Pr = 9.47 and different Richardson numbers. From left to right and from top to the bottom Ri = 1, 0.1, 0.01 and 0.001 respectively

13

Heat Mass Transfer

Fig. 10  Distribution of the Nusselt number along the lid for Prandtl number Pr = 56 and different Richardson numbers. From left to right and from top to the bottom Ri = 1, 0.1, 0.01 and 0.001 respectively

flattens and curve slope near the left vertical wall increases. One can observe that scraping effect occurs and increases 10 times the Nu maximum (for Ri = 0.001) in comparison with the case where Ri = 1.0. Figure 9 shows local variations of Nu along the top wall for Pr = 9.47 and Ri ranges from 1.0 to 0.001. It is clear that Nu increases with increasing Pr. The curves of local Nu distribution are very similar in shape. In comparison with case of Pr = 0.71, local Nu plots are more flattened and Nu values are higher. Sharp peak of Nu near the left top cavity corner is much more visible compared to case of Pr = 0.71. This is caused by superposition of the scraping effect and fluid impingement on the top wall near the left top corner. The second phenomena is triggered out by buoyancy forces and fluid circulation. In the cases where Pr is higher, buoyancy forces are stronger and superimposing of mentioned effects results in significant increase of local Nu maximum. For Ri values smaller than 0.01 fluid circulation is very strong and forced convection dominates (even for Pr = 56.00). With decreasing Ri local Nu values

increase. Irrespective of decreasing Ri abrupt increasing of Nu always occurs in the range of X = 0–0.1. After reaching X = 0.1 the curve slope monotonically decreases. This means that the area size, where the highest Nu values occur, does not depend on Ri. The same conclusions can be adopted in case of Pr = 56.00. The results for this case are depicted in Fig. 10. Consequently the area size, where the highest Nu values occur, does not depend on Pr as well. This is very important conclusion because it shows that mechanism of temperature boundary layer developing is similar in each studied cases and this is governed by the universal law. Additionally increasing Pr does not change the character of local Nu. The only result is substantial increase of Nu and area averaged Nusselt number Nu. In Table 4 average values of Nusselt numbers over the surface of the top lid are presented. With increasing Pr and decreasing Ri rate of the heat transfer increases. This is confirmed by increasing Nu. For the smallest values of Ri heat transfer rate is almost 5 times larger for Pr = 56.00 case compared to Pr = 0.71 case. With increasing Pr and decreasing Ri

13



Heat Mass Transfer

Despite the simplicity of this case it was not considered in the literature. Calculations involved a wide range of Richardson number and Prandtl number, while the Grashof number was fixed. The main goal was to identify the effects of Richardson number and Prandtl number on the area averaged Nusselt number. Numerical simulations were carried out on very fine meshes (resolution 1024 × 1024) to provide benchmark solution. Based on the obtained results the following findings can be summarized:

Fig. 11  Variation of the average Nusselt number with respect to Richardson number calculated at the lid

Table 4  Average Nusselt numbers Nu on the top and bottom wall Pr

0.71

9.47

56.00

Ri

Error (%)

Nu

Top

Bottom

1.0 0.1 0.01 0.001 1.0 0.1 0.01 0.001

2.5789 3.7763 6.2027 10.2927 5.3224 7.7682 13.9427 26.3319

2.5808 3.7788 6.2132 10.3495 5.4943 8.0143 14.3995 27.3229

0.07 0.07 0.17 0.55 3.13 3.07 3.17 3.63

1.0 0.1 0.01

8.8690 14.0100 26.5014

9.2233 14.4181 27.3266

3.84 2.83 3.02

0.001

50.3581

52.5920

4.25

time needed for simulation performing also lengthened. Furthermore problems with fulfilment of the condition defined by Eq. 12 in terms of Nu appeared, and error began to rise. For the limiting case in this study, i.e. for Pr = 56.00 and Ri = 0.001 the error was maximal but <5 %. In Fig. 11 the calculated values of Nu at the top wall, in terms of Ri, are depicted. With decreasing Ri heat transfer rate increases and the curve slope increases with increasing Pr.

6 Conclusions In this study numerical modelling of mixed convection heat transfer and fluid flow in a lid-driven cavity with constant heat flux applied on the bottom wall was performed.

13

1. In the studied range of Richardson number, Prandtl number does not have an impact on flow field pattern but strongly influences temperature field and heat transfer phenomena. 2. In the case of constant Grashof number heat transfer rate increases with increasing Prandtl number and decreasing Richardson number. 3. For Richardson number values smaller than 0.01 forced convection dominates. 4. The highest values of Nusselt number are observed in the left top cavity corner. This is caused by combined effect of two mechanisms: scraping of temperature boundary layer phenomenon and fluid impingement with the top wall triggered out by circulation fluid movement and buoyancy forces. 5. The size of the area where the highest values of Nusselt numbers occurs does not depend on Richardson number and Prandtl number, which implies a universal law of rebuilding of the temperature boundary layer. 6. Increasing Prandtl number results in increase of temperature gradient and in reducing the thickness of the temperature boundary layer. 7. For low Ri values the temperature in the cavity centre is uniformly distributed. The same effect was noticed by Guo and Sharif [15]. Acknowledgments  Calculations have been carried out using resources provided by Wroclaw Centre for Networking and Supercomputing (http://wcss.pl, the division of Wroclaw University of Technology). This work is co-financed by the European Union as part of the European Social Fund.

References 1. Shankar P, Despande DE (2000) Fluid mechanics in the driven cavity. Annu Rev Fluid Mech 32:93–136 2. Bruneau CH, Saad M (2006) The 2d lid-driven cavity problem revisited. Comput Fluids 35:326–348 3. Vierendeels J, Merci B, Dick E (2003) Benchmark solutions for the natural convective heat transfer problem in a square cavity with large horizontal temperature differences. Int J Numer. Methods Heat Fluid Flow 13:1057–1078 4. Moallemi MK, Jang KS (1992) Prandtl number effects on laminar mixed convection heat transfer in a lid-driven cavity. Int J Heat Mass Transf 35:1881–1892

Heat Mass Transfer 5. Aydin O, Yang WJ (2000) Mixed convection in cavities with a locally heated lower wall and moving sidewalls. Numer Heat Transf Part A 37:695–710 6. Prasad YS, Das MK (2007) Hopf bifurcation in mixed convection flow inside a rectangular cavity. Int J Heat Mass Transf 50:3583–3598 7. Basak T, Roy S, Sharma PK, Pop I (2009) Analysis of mixed convection flows within a square cavity with uniform and nonuniform heating of bottom wall. Int J Therm Sci 48:891–912 8. Oztop HF, Dagtekin I (2004) Mixed convection in two-sided liddriven differentially heated square cavity. Int J Heat Mass Transf 47:1761–1769 9. Santos ED, Piccoli GL, Franca FHR, Petry AP (2011) Analysis of mixed convection in transient laminar and turbulent flows in driven cavities. Int J Heat Mass Transf 54:4585–4595 10. Papanicolaou E, Jaluria Y (1991) Mixed convection from an isolated heat source in a rectangular enclosure. Numer Heat Transf Part A 18:427–461 11. Papanicolaou E, Jaluria Y (1992) Transition to a periodic regime in mixed convection in a square cavity. J Fluid Mech 239:489–509 12. Papanicolaou E, Jaluria Y (1993) Mixed convection from a localized heat source in a cavity with conducting walls: a numerical study. Numer Heat Transf Part A 23:463–484 13. Papanicolaou E, Jaluria Y (1994) Mixed convection from simulated electronic components at varying relative positions in a cavity. J Heat Transf 116:960–970 14. Hsu TH, Wang SG (2000) Mixed convection in a rectangular enclosure with discrete heat sources. Numer Heat Transf Part A 38:627–652 15. Guo G, Shariff MAR (2004) Mixed convection in rectangular cavities at various aspect ratios with moving isothermal sidewalls

and constant flux heat source on the bottom wall. Int J Therm Sci 43:465–475 16. Ogut EB (2010) Mixed convection in an inclined lid-driven enclosure with a constant flux heater using differential quadrature (dq) method. Int J Phys Sci 5:2287–2303 17. Shahi M, Mahmoudi AH, Talebi F (2010) Numerical study of mixed convective cooling in a square cavity ventilated and partially heated from the below utilizing nanofluid. Int Commun Heat Mass Transf 37:201–213 18. Mansour MA, Mohamed RA, Abd-Elaziz MM, Ahmed SE (2010) Numerical simulation of mixed convection flows in a square lid-driven cavity partially heated from below using nanofluid. Int Commun Heat Mass Transf 37:1504–1512 19. Cho CC, Chen CL, Chen CK (2013) Mixed convection heat transfer performance of water-based nanofluids in lid-driven cavity with wavy surfaces. Int J Therm Sci 68:181–190 20. Ganzarolli MM, Milanez LF (1995) Natural convection in rectangular enclosures heated from below and symmetrically cooled from the sides. Int J Heat Mass Transf 38:1063–1073 21. Ferziger JH, Peric M (2002) Computational methods for fluid dynamics, 3rd edn. Springer, Berlin 22. ANSYS® Academic Research, Release 14.5, Help System, Solver Theory Guide, ANSYS Inc 23. Patankar SV (1980) Numerical heat transfer and fluid flow, 1st edn. McGraw-Hill Book Company, New York 24. Rhie C, Chow W (1982) A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation. AIAA-820998 Paper, vol 6, pp 1982–1998 25. Lage JL, Bejan A (1991) The Ra-Pr domain of laminar natural convection in an enclosure heated from the side. Numer Heat Transf Part A 19:21–41

13