Abstract
In this study a steady state twodimensional mixed convection problem in an air filled cavity is investigated. The effects of linearly heated and uniformly cooled walls on flow and heat transfer characteristics within the cavity are determined. The strength of the fluid circulation within the cavity is found for both heated and cooled walls for different Reynolds and Grashof numbers. The nonlinear coupled equations are solved numerically using the penaltyGalerkin finite element method. Stream function and isotherm results are obtained for different Reynolds and Grashof numbers. The results for the heat transfer rate are presented in terms of both the local and the average Nusselt number. In general, the strength of the circulation is stronger for the case of cooled walls, and the anticlockwise circulation is significantly stronger for cooled walls, while the clockwise circulation is only slightly stronger for the cooled walls and the difference in strength decreases both with increasing Reynolds and Grashof numbers. Both the local and average Nusselt numbers are generally higher for the case of cooled side walls than that for heated side walls.
MSC: 34B15, 65N30, 76M20.
1 Introduction
Over the last decade mixed convection in a square cavity has become an increasingly attractive field of study. The popularity of such studies is due in part to the many applications of such flows in industrial and natural settings. These applications include, for example, the thermalhydraulics of nuclear reactors, drying technologies, the dynamics of lakes, and food processing. Mixed convection in a liddriven cavity has been extensively studied. Torrance et al.[1] investigated the fluid motion in a liddriven cavity where the upper wall was maintained at a different temperature to that of the other walls in the cavity. They observed that an increase in the aspect ratio led to an increase in secondary circulations in the lower part of the cavity.
Kawaguti [2] was one of the first to investigate the effect of the Reynolds number on forced convection in a square cavity. Results were found using finite difference methods for values of the Reynolds number Re between 0 and 64 for different ratios of lengths of the cavity. An attempt was made to find results for , but a convergent solution was not found. By modifying Kawaguti’s method, Burggraf [3] was able to attain stable solutions for forced convection for a Reynolds number up to 1,000. Heinrich and Marshall [4] used the penalty finite element method to investigate liddriven cavity flow up to . For a low Reynolds number, their results were in excellent agreement with published results. As the Reynolds number increased, however, their results were less accurate in comparison to other results. Early studies of natural convection in a cavity were of an experimental nature. An example of an early experimental study is that of Eckert and Carlson [5]. They made an experimental investigation into the flow and heat transfer in an enclosed cavity. The top and bottom plates were heated to different temperatures. Their results showed that, contrary to previous beliefs, it was very unlikely that a core of uniform temperature could exist at large Rayleigh numbers unless the height to width aspect ratio was very small. Ostrach [6] gives an excellent review of early experimental studies of convection in a cavity. Recently, Prasad and Koseff [7] investigated mixed convection in a liddriven cavity through experimental means where the lower surface was heated and the upper surface cooled. The lid speed and values of the Reynolds and Grashof numbers were varied. Their results indicated that the overall heat transfer rate was a very weak function of the Grashof number for the range of Reynolds numbers examined.
Moallemi and Jang [8] investigated flow in a liddriven cavity with the bottom wall heated. They studied the effects of small Prandtl numbers on the flow and heat transfer for various values of the Richardson number. Sivakumar et al.[9] investigated mixed convection in a liddriven cavity with a cooled right wall and a heat source in the left wall. Three different lengths of the heat source were examined. The finite volume method was used to solve the resulting equations. It was found that for low values of the Richardson number, reducing the heating portion length had no effect on the fluid flow. However, on increasing the Richardson number, it was seen that the flow depended heavily on the length of the heat source. It was also found that reducing the length of the heated portion of the wall resulted in a better heat transfer rate. Nithyadevi et al.[10] studied the case where the side walls were partially heated or partially cooled while the rest of the cavity was well insulated. They used the finite volume method to solve the equations. They studied nine different cases where the position of the heating and cooling sources was moved between the top, bottom and middle of the wall of the cavity. They found that the heat transfer rate was enhanced when there was a cooling source near the top of the cavity. For a combination of a bottom top heating source, the heat transfer rate was found to be high, while for the combination of a top bottom heating source, the heat transfer rate was low. Related studies have been done by, among others, Paraconi and Corvaro [11] and Pesso and Piva [12], Sathiyamoorthy et al.[13].
Recently, Basak et al.[14,15] studied mixed convection in a liddriven cavity. The moving upper lid was insulated in both studies. In Basak et al.[14] they investigated cooled side walls with a heated (uniformly and nonuniformly) bottom wall. In Basak et al.[15] they looked at a uniformly heated bottom wall with linearly heated left side wall. The right side wall was heated either linearly or cooled uniformly. Both papers used the penalty finite element method to solve the governing equations of motion and heat transfer. Both studies found that the strength of convection increased with increasing Grashof numbers and that liddriven flow was dominant for . Their work showed that the heat transfer rate for a uniformly heated bottom wall was higher than that for the nonuniformly heated bottom wall and that the heat transfer rate was larger for the case of one cooled right wall than it was for the case when both walls were linearly heated.
Research has also been done into mixed convection with more than one moving wall in a cavity. Oztop and Dagtekin [16] used the finite volume method to investigate three cases of moving walls with adiabatic top and bottom wall and cooled left wall and heated right wall. They observed that the fluid flow and the heat transfer in the cavity depended on both the Richardson number Ri and the direction of motion of the walls.
Cheng and Liu [17] investigated four cases of mixed convection in a square cavity. In the first case, the side walls were well insulated while the top moving wall was heated and the bottom wall cooled. In the second case, the bottom wall was heated and the top wall cooled. In the third and fourth cases, the top and bottom walls were well insulated and one side wall was cooled while the other was heated. They found that for the first case, when the Richardson number was greater than 1, the heat transfer was mainly through conduction. When the Richardson number was less than 1, forced convection dominated the fluid flow as was found by Oztop and Dagtekin [16]. In contrast to the findings of Oztop and Dagtekin, in the second case, when the Richardson number was greater than 1, the lower half of the cavity showed natural convection while the upper part was dominated by forced convection. When the Richardson number was less than 1, the fluid flow was similar to that for case 1. For the third case when the Richardson number was greater than 1, heat transfer was shown to be by conduction near the side walls and by convection in the center. Again when the Richardson number was less than one, the fluid flow was similar to that for cases 1 and 2. For the fourth case, a large portion of the cavity was dominated by natural convection for all values of the Richardson number.
Corcione [18] investigated the effect of the Rayleigh number and the width to height aspect ratio of the cavity on steady laminar natural convection in an air filled cavity. The cavity was heated from below and cooled from above while six different cases of side wall heating and cooling were investigated. They showed that the heat transfer rate from a cooled or heated boundary increased as the Rayleigh number increased. In the case with insulated side walls, the heat transfer rate from the heated bottom wall or top cooled wall increased as each of the side walls was replaced by a cooled or heated side wall. Wong [19] investigated mixed convection in a square cavity. The left wall was maintained at a constant velocity. The other walls were kept stationary. The top and bottom walls were insulated and the left moving wall was heated while the right side wall was cooled. A numerical solution was found using the consistent splitting scheme and the finite element method. Wong fixed the Reynolds number at 100 and varied the Grashof number to measure the effects of the Richardson number. When the Richardson number was 0.001, the flow was found to be strongly influenced by forced convection. At , mixed convection was found to occur. At , the fluid flow was mainly dominated by natural convection. These findings are similar to those of Oztop and Dagtekin [16] even though only one side wall was moving in this case. Chinyoka [20] investigated nonisothermal twodimensional unsteady fluid flow problems in a rectangular channel.
In this study mixed convection in a twodimensional cavity is investigated using the penaltyGalerkin finite element method to solve the governing momentum and heat equations. Much work has been done on mixed convection with uniformly heated or cooled walls. Since little work has been done on sinusoidal and linearly heated walls and even less work on a combination of these, the aim of this investigation is to combine linear and sinusoidal heating of different walls of the cavity with a moving lid to study the resulting mixed convective flow.
2 Mathematical modeling
A twodimensional square cavity of length L is considered for the current investigation. The cavity lid is sinusoidally heated while moving at a constant speed from left to right, and the cavity is heated at a constant rate from the bottom. Two cases of side wall heating are investigated. In the first instance, the cavity side walls are linearly heated, while in case 2 the side walls are uniformly cooled.
The equations that describe the fluid flow and the heat transfer characteristics within the cavity, subject to the Boussinesq approximation are as follows:
subject to the boundary conditions
where u and v are the velocity components in the x and y directions respectively, α is the thermal diffusivity, β is the coefficient of thermal expansion, ν is the kinematic viscosity, ρ is the mass density, g is the acceleration due to gravity, is the highest temperature attained along the bottom wall and the middle of the top wall and is the lowest temperature attained at the top corners in case 1 and along the side walls in case 2.
We nondimensionalize the equations by using the following change of variables:
where is the velocity of the upper wall and θ is the nondimensional fluid temperature. This leads to the equations
where the important parameters are the Prandtl number Pr, the Reynolds number Re, and the Grashof number Gr given, respectively, by
The corresponding boundary conditions are:
The geometry of the problem together with the associated boundary conditions is shown in Figure 1.
Figure 1. Schematic sketch of the cavity problem with, Case 1, linearly heated side walls and, Case 2, uniformly cooled walls.
3 Numerical solution
Equations (8)(11) were solved using the penaltyGalerkin finite element method (PGFEM). This method has the advantage of eliminating pressure as a dependent variable while still satisfying the continuity equation.
We introduce test functions (which are piecewise, once differentiable functions) and multiply equations (8)(11) by the test functions to obtain the weak forms of the equations, namely
The continuity equation in its weighted form becomes
The PGFEM supposes that we may replace equation (18) with
where ϵ is an arbitrarily small parameter. Rearranging these equations, we find
where λ is the penalty parameter. Mass conservation is satisfied in the limit or . Here λ is chosen to be 10^{7} as this has been shown to give consistent results (Basak et al.[14]). If required, the pressure field can be calculated using equation (20). This generally requires a higher level of convergence than is necessary for the velocity field. In this investigation we have chosen to restrict our results to the stream function and temperature contours. As an alternative to the penalty function technique above, the semiimplicit method for pressure linked equations (SIMPLE) technique may be used; see, for instance, Chinyoka [2022].
Substituting equation (20) into momentum equations (15) and (16) leads to the penalized momentum equations:
To determine the finite element approximate solution of equations (17), (21) and (22), we discretise the domain into 400 rectangular elements with a total of 441 nodes. In principle, more elements may be used to ensure even greater accuracy of results, but this greatly increases the computation time. Initially, one hundred elements were used, but although the computation of such a formulation was efficient, the accuracy of the solutions was impaired. Four hundred elements were found to give sufficient accuracy without impairing the efficiency of the solution method. There are 361 nodes which are not on a Dirichlet boundary. Associated with each node in the domain is a basis function . For convenience, the same functions were used as both test and basis functions. These basis functions are chosen to be linear Lagrange elements. We expand U, V and θ as linear combinations of the basis functions over the nonDirichlet boundary nodes as
where , and are the nodal values of U, V and θ at the jth node.
By substituting these expressions for U, V and θ into equations (17), (21) and (22), we obtain (after manipulation and simplification)
The set of nonlinear equations is solved using the Gaussian quadrature and reduced integration to prevent locking, that is, the deterioration in performance of the numerical scheme as (see, for instance, Qi et al.[23] and the references therein).
3.1 The stream function and the Nusselt number
Results of the fluid motion are usually interpreted in terms of the stream function ψ, where
Expanding the stream function in the same way as we did for U, V and θ, we obtain
The heat transfer coefficient in terms of the Nusselt number is defined as
where n denotes the normal to the plane. The average Nusselt number is defined as
4 Results and discussion
Results have been found for quadratic elements with 361 inner nodes, (the value for air), , and . Convergence was found to be generally poor for .
Benchmark results were found for uniformly heated bottom and cooled side walls with an adiabatic stationary lid. Figure 2 shows the stream function and temperature contours for the benchmark results. The results are in good agreement with those in the literature, particularly Basak et al.[24]. We note that in the figures below, a negative stream function value denotes clockwise flow while a positive value denotes anticlockwise flow.
Figures 35 show the results for the stream function and temperature contours for a moving sinusoidally heated lid with linearly heated side walls and a uniformly heated bottom wall for Grashof numbers between 10^{3} and 10^{5}.
Figure 3. Stream function and temperature contours for linearly heated side walls when, (a) , (b), (c).
Figure 3 shows that for , , two counter rotating circulations are formed in the cavity. The circulation on the lefthand side of the cavity is counterclockwise, while the circulation on the righthand side of the cavity is clockwise. The clockwise flow is significantly stronger than the anticlockwise flow and thus it can be observed that the flow is dominated by forced convection due to the moving lid. As was expected, the isotherms for , reach a maximum along the bottom wall and at the center of the top wall, while reaching a minimum at the top corners of the cavity. At the bottom of the cavity, isotherms span the entire width of the cavity. When , the isotherms instead curve upwards towards the lid of the cavity. Hot isotherms at the top of the cavity are restricted to the center of the top wall.
As Gr increases to 10^{4}, the anticlockwise circulation increases in size and in strength so that two symmetric rotations are formed within the cavity. Natural convection is now as equally dominant as the forced convection regime. The isotherm pattern is similar to that of , although the heat from the bottom wall has risen slightly higher into the middle of the cavity.
At , the clockwise circulation in the cavity increases in size and dominates most of the cavity. The circulation is slanted towards the top right and bottom left corners. A small third anticlockwise circulation is formed on the bottom righthand side of the cavity. It can be seen that the isotherms are compressed in the top corners of the cavity and in the middle of the bottom wall. The isotherms are widely dispersed in the lefthand corner of the cavity. This is due to the shape and angle of the clockwise circulation. Cool fluid is brought down on the righthand side of the cavity, while hot fluid is transported upwards towards the middle of the cavity on the lefthand side of the bottom wall. The strength of the circulations increase with increasing Grashof numbers.
For and (see Figure 4), the flow in the cavity is clearly dominated by forced convection. Only a very small secondary circulation is found for this low Grashof number. The isotherms are no longer perfectly symmetrical. As Gr increases to 10^{4}, the secondary circulation on the lefthand side of the cavity grows in size and strength. The strength of the clockwise circulation is stronger than that of the anticlockwise circulation. The flow is still dominated by liddriven flow, but natural convection now has a more significant effect on the flow. An increase in the Grashof number causes the heat to rise further towards the center of the cavity. When , the secondary circulation increases in size, becomes slanted and dominates most of the cavity and a small third clockwise circulation is formed on the bottom righthand side of the cavity. The isotherms at are also similar in shape to those when .
Figure 4. Stream function and temperature contours for linearly heated side walls when, (a) , (b), (c).
In Figure 5 at , and , the flow is again dominated by forced convection. A single clockwise circulation is formed due to strong inertial effects at the top wall. There is little change in the flow as the Grashof number increases from 10^{3} to 10^{4}. The isotherms for both and are dispersed on the lefthand side of the cavity and are compressed on the righthand side of the cavity. This is due to the domination of the forced convection flow regime. At , , a large slanted clockwise circulation is formed which dominates more of the cavity than the circulations when and when . Again a small third anticlockwise circulation is formed at the bottom righthand corner of the cavity. The isotherms are similar in shape for and .
Figure 5. Stream function and temperature contours for linearly heated side walls when, (a) , (b), (c).
For (see Figure 6), the flow is dominated by forced convection for , and . For , a single clockwise circulation is formed within the cavity. As Gr increases to 10^{4}, a small anticlockwise circulation is formed on the bottom righthand side of the cavity. As Gr increases to 10^{5}, there is little change in the fluid flow. The isotherms for both and are dispersed on the bottom lefthand side of the cavity and compressed on the bottom righthand side of the cavity. The warm isotherms are dispersed towards the top righthand side of the cavity, causing the cool isotherms to be compressed into the top right corner. As Gr increases to 10^{5}, the hot isotherms stretch slightly higher into the center of the cavity. It can be seen that as the Reynolds number increases, the strength of the circulation decreases.
Figure 6. Stream function and temperature contours for linearly heated side walls when, (a) , (b), (c).
Figures 79 show the results for the stream function and temperature contours for a moving sinusoidally heated lid with cooled side walls and uniformly heated bottom wall for , and a range of Grashof numbers between 10^{3} and 10^{5}. For and (Figure 7), the flow is slightly dominated by forced convection. Two counter rotating circulations are formed in the cavity, with the clockwise circulation being only slightly stronger than the anticlockwise circulation. The temperature distribution across the cavity shows that the cavity is at a high temperature across the bottom of the cavity and towards the center of the top of the cavity. The cavity is cooler towards the side walls and top corners as was expected. The hot isotherms still span the width of the cavity, but are pushed downwards in the corners. For , the isotherms are horizontal across the cavity, but for , the isotherms stretch vertically across the cavity. The temperature distribution obviously differs from that of the linearly heated walls at the side walls. It can also be seen that the contour lines are more widely dispersed towards the bottom of the cavity for the linearly heated case.
Figure 7. Stream function and temperature contours for cooled side walls when, (a), (b), (c).
As Gr increases to 10^{4}, the clockwise circulation increases in size to form two counter rotating circulations of a similar size and the natural convection is almost equally dominant to that of forced convection. The centers of the circulations are now at the same height. The temperature contours become slightly more compressed towards the side walls as the Grashof number increases. The hotter contour lines from the top wall are seen to be wider than that of and are found to be narrower in the center towards the bottom wall. As Gr increases to 10^{5}, the strength of the stream function increases and the size and position of the two circulations are identical. It can be seen that the hot contour lines have widened even further along the top wall causing the cooler contour lines to become compressed towards the top of the cavity. Towards the bottom of the cavity, the cool isotherms (those where ) are widely dispersed towards the center of the cavity while the warmer isotherms are compressed towards the bottom wall. This is particularly the case towards the bottom corners of the cavity.
In Figure 8 it is clear that for , , the flow is dominated by forced convection. The majority of the cavity is dominated by a clockwise rotating circulation. A small anticlockwise circulation is formed on the bottom lefthand side of the cavity. The temperature contours are similar to the case where . On close inspection it can be seen that the contours are slightly more compressed towards the top right corner than when .
Figure 8. Stream function and temperature contours for cooled side walls when, (a), (b), (c).
As Gr increases to 10^{4}, the anticlockwise circulation grows in size although the clockwise rotating circulation is still larger in size and in strength indicating that the flow is still slightly dominated by forced convection. Again the contour lines are similar in shape to those for , and a slight compression towards the top righthand corner can be seen. As Gr increases further to 10^{5}, natural convection becomes equal to forced convection. There is no observable difference in the shape of the temperature contours from that of .
Figure 9 shows that for and and , the flow is dominated by forced convection. For , only one clockwise rotating circulation appears in the cavity. As Gr increases to 10^{4}, a small anticlockwise circulation is formed in the bottom lefthand corner of the cavity. As Gr increases to 10^{5}, the anticlockwise circulation increases in size and strength. The flow is still slightly dominated by forced convection. The temperature contours for show a marked difference to those for and . For , the hot isotherms (where ) are dispersed towards the bottom left and top right corners and compressed in the top left and bottom right corners. The cool isotherms are compressed in the top right and bottom left corners and dispersed towards the center of the cavity along the middle and bottom of the right wall and at the top of the left wall. The compression and dispersion is more significant along the right side wall than along the left side wall. As Gr increases to 10^{4}, the hot temperature contours from the bottom wall move slightly higher towards the center of the cavity, while the basic shape of the contours remains the same as those for . A further increase in Gr to 10^{5} causes the temperature contours to become compressed in the top two corners of the cavity and more dispersed in the bottom two corners. There is a significant compression of isotherms towards the top righthand corner. The hot temperature contours from the bottom wall become more elongated and extend further towards the center of the cavity.
Figure 9. Stream function and temperature contours for cooled side walls when, (a), (b), (c).
Figure 10 shows that for , and , the flow is dominated by forced convection with a single clockwise circulation being formed within the cavity. As Gr increases to 10^{5}, an anticlockwise circulation is formed on the bottom lefthand side of the cavity. The isotherms for and are similar. The hot isotherms are dispersed on the lefthand side of the cavity and compressed on the righthand side of the cavity. Warm isotherms () from the top of the cavity are stretched downwards towards the righthand side of the cavity. As Gr increases to 10^{5}, the warm isotherms stretch slightly higher into the center of the cavity than for and . The warm isotherms from the top of the cavity stretch slightly further downwards towards the center righthand side of the cavity.
Figure 10. Stream function and temperature contours for cooled side walls when, (a), (b), (c).
In general, the circulation is stronger for the case of cooled side walls than that for linearly heated walls. The anticlockwise circulation is significantly stronger for the cooled walls than for the linearly heated walls while the clockwise circulation is only slightly stronger for the cooled walls. For , the cooled walls produce two identical counter rotating circulations in contrast to the three tilted circulations resulting from the linearly heated walls. For , the stream function pattern for the cooled walls shows a slightly stronger clockwise circulation indicating a slight domination by forced convection, but again the flow pattern is different from that of the linearly heated side walls.
5 Heat transfer at the walls
Figure 11 shows the effect of the Grashof number on the heat transfer at each of the walls in the cavity for the case of linearly heated side walls.
Figure 11. Local Nusselt number for linearly heated side walls withand(+),(△) and(o).
For all values of the Grashof number, the local Nusselt number is equal to one at the edge of the bottom wall on both sides due to the linear heating of the side walls. For , the local Nusselt number falls below 1 and reaches a local minimum between and . There is not much variation in the Nusselt number across the cavity due to the relatively flat isotherms at the bottom of the cavity. When , the Nusselt number is of a similar value to that for between and . Beyond , the graph rises above that of reaching a maximum at .
When , the Nusselt number is greater than 1 for the majority of the cavity due to the compressed isotherms. The Nusselt number drops below 1 between and due to the widely dispersed isotherms on the lefthand side of the cavity. The Nusselt number is at a maximum at . The local Nusselt number at the edges of the top wall is equal to −1 for all values of the Grashof number due to the linearly heated walls. For , the Nusselt number rises and reaches a local maximum at and then gradually decreases. For , the Nusselt number follows a similar pattern to that for , although is lower than the graph of . At , the graph of the Nusselt number is close to being symmetric in shape. The Nusselt number drops to a local minimum at and then rises to a local maximum at the center of the cavity. The Nusselt number then drops again to a local minimum at . When , the local Nusselt number at the righthand wall gradually increases as y increases. For , a similar occurrence happens, although the Nusselt number is below that for until where the Nusselt number becomes greater than that of . For , the Nusselt number decreases below 0 attaining a local minimum at due to the compressed isotherms at the bottom corner of the right wall. The Nusselt number then increases reaching a local maximum at .
The local Nusselt number at the lefthand wall for and is very similar in shape to that for the righthand wall due to a near symmetrical isotherm pattern at the two walls. However at , the local Nusselt number at the lefthand wall is above 0 for the entire height of the cavity and has an oscillatory pattern. The Nusselt number initially rises attaining a local maximum at due to the dispersed isotherms along the bottom half of the lefthand wall. Towards the top of the cavity, the isotherms become more dispersed once again and the Nusselt number again rises to a local maximum at .
The graph at the lefthand wall is similar to that of Basak et al.[15] which also exhibits a slightly oscillatory shape. However, the graph in the present study has a local maximum near the top corner and then decreases in contrast to the graph in Basak et al.[15] which increases monotonically to the wall. This is due to the difference between the adiabatic and the heated lid.
Figure 12 shows the average Nusselt number at the four walls of the cavity. At the bottom wall, the Nusselt number is seen to increase as the Grashof number increases. The graph of the average Nusselt number at the bottom wall is almost identical to the results of Basak et al.[15]. At the top wall, the opposite is seen to happen. This is due to the fact that as the Grashof number increases, the isotherms become more compressed in the top corners. At the lefthand side wall, the average Nusselt number is lower for than it is for . This is due to the slight dip in the isotherms at the left wall. The average Nusselt number at the left wall then increases for . The average Nusselt number at the right side wall decreases as the Grashof number increases. The local Nusselt number at the bottom wall for cooled side walls is shown in Figure 13. For all values of the Grashof number, the Nusselt number decreases to where it attains a local minimum. The local Nusselt number when is lower than that when which, in turn, is lower than that when . In contrast to the case of linearly heated side walls, the local Nusselt number at the bottom wall for cooled walls is concave up due to the fact that the isotherms are equal to 0 at the corners of the bottom wall for the cooled walls.
Figure 12. The average Nusselt number for linearly heated side walls with.
Figure 13. The local Nusselt number for cooled side walls with,(+),(△) and(o).
The local Nusselt number decreases to a local minimum at for and then increases slightly towards the top of the cavity. For , the Nusselt number decreases to a local minimum near and then slowly increases towards the top of the cavity. When , the local Nusselt number decreases to a local minimum near due to the dispersion of isotherms near the bottom corner(s). The Nusselt number then increases above that of and at and , respectively, due to the more compressed isotherms in the top corners. The Nusselt number then increases to a local maximum near . There is no oscillatory pattern when as was observed in the case of the heated side walls.
6 Conclusion
A numerical study has been performed using the penaltyGalerkin finite element method to analyze mixed convective heat transfer and fluid flow in an air filled square cavity. The effects of linearly heated and uniformly cooled walls on flow and heat transfer characteristics within the cavity have been studied. It was observed that
– In general, the strength of circulation was stronger for the case of cooled walls, and that the anticlockwise circulation was significantly stronger for cooled walls, while the clockwise circulation was only slightly stronger for cooled walls. The difference in strength decreased with increasing both Reynolds and Grashof numbers.
– Both the local and average Nusselt numbers were generally higher for the case of cooled side walls than those for heated side walls.
– The local Nusselt numbers at the left and right walls were similar in the case of cooled side walls due to symmetric patterns in the temperature isotherms. However, particularly for , the heated side walls had different local Nusselt numbers at the two different walls.
– The local Nusselt number at the top wall was similar in shape for the two cases although the linearly heated case was lower in value than for the cooled side walls.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
DSD carried out the numerical computations. PS participated in the design of the study and helped to draft the manuscript.
Acknowledgements
The authors are grateful to the University of KwaZuluNatal for financial support.
References

Torrance, K, Davis, R, Eike, K, Gill, P, Gutman, D, Hsui, A, Lyons, S, Zien, H: Cavity flows driven by buoyancy and shear. J. Fluid Mech.. 51, 221–231 (1972). PubMed Abstract  Publisher Full Text

Kawaguti, M: Numerical solution of the NavierStokes equations for the flow in a twodimensional cavity. J. Phys. Soc. Jpn.. 16, 2307–2315 (1961). Publisher Full Text

Burggraf, OR: Analytical and numerical studies of the structure of steady separated flows. J. Fluid Mech.. 24, 113–151 (1966). Publisher Full Text

Heinrich, JC, Marshall, RS: Viscous incompressible flow by a penalty function finite element method. Comput. Fluids. 9, 73–83 (1981). Publisher Full Text

Eckert, EG, Carlson, WO: Natural convection in an air layer enclosed between two vertical plates with different temperatures. Int. J. Heat Mass Transf.. 2, 106–110 (1961). Publisher Full Text

Ostrach, S: Natural convection in enclosures. Adv. Heat Transf.. 7, 161–227 (1972)

Prasad, AK, Koseff, JR: Combined forced and natural convection heat transfer in a deep liddriven cavity flow. Int. J. Heat Fluid Flow. 17, 460–467 (1996). PubMed Abstract  Publisher Full Text

Moallemi, M, Jang, KS: Prandtl number effects on laminar mixed convection heat transfer in a lid driven cavity. Int. J. Heat Mass Transf.. 35, 1881–1892 (1992). PubMed Abstract  Publisher Full Text

Sivakumar, V, Sivasankaran, S, Prakash, P, Lee, J: Effect of heating location and size on mixed convection in liddriven cavities. Comput. Math. Appl.. 59, 3053–3065 (2010). Publisher Full Text

Nithyadevi, N, Kandaswamy, P, Lee, J: Natural convection in a rectangular cavity with partially active side walls. Int. J. Heat Mass Transf.. 50, 4688–4697 (2007). PubMed Abstract  Publisher Full Text

Paraconi, M, Corvaro, F: Natural convection in a square enclosure with a hot source. Int. J. Therm. Sci.. 48, 1683–1695 (2009). Publisher Full Text

Pesso, T, Piva, S: Laminar natural convection in a square cavity: low Prandtl numbers and large density differences. Int. J. Heat Mass Transf.. 52, 1036–1043 (2009). Publisher Full Text

Sathiyamoorthy, M, Basak, T, Roy, S, Mahanti, NC: Effect of the temperature difference aspect ratio on natural convection in a square cavity for nonuniform thermal boundary conditions. J. Heat Transf.. 129, 1723–1728 (2007). Publisher Full Text

Basak, T, Roy, S, Sharma, PK, Pop, I: Analysis of mixed convection flows within a square cavity with uniform and nonuniform heating of bottom wall. Int. J. Therm. Sci.. 48, 891–912 (2009). Publisher Full Text

Basak, T, Roy, S, Sharma, PK, Pop, I: Analysis of mixed convection flows within a square cavity with linearly heated side wall(s). Int. J. Heat Mass Transf.. 52, 2224–2242 (2009). Publisher Full Text

Oztop, H, Dagtekin, I: Mixed convection in twosided liddriven differentially heated square cavity. Int. J. Heat Mass Transf.. 47, 1761–1769 (2004). Publisher Full Text

Cheng, TS, Liu, WH: Effect of temperature gradient orientation on the characteristics of mixed convection flow in a liddriven square cavity. Comput. Fluids. 39, 965–978 (2010). Publisher Full Text

Corcione, M: Effects of the thermal boundary conditions at the sidewalls upon natural convection in rectangular enclosures heated from below and cooled from above. Int. J. Therm. Sci.. 42, 199–208 (2003). Publisher Full Text

Wong, JC: Numerical simulation of twodimensional laminar mixedconvection in a liddriven cavity using the mixed finite element consistent splitting scheme. Int. J. Numer. Methods Heat Fluid Flow. 17, 46–93 (2007). Publisher Full Text

Chinyoka, T: Twodimensional flow of chemically reactive viscoelastic fluids with or without the influence of thermal convection. Commun. Nonlinear Sci. Numer. Simul.. 16, 1387–1395 (2011). Publisher Full Text

Chinyoka, T: Viscoelastic effects in doublepipe singlepass counter flow heat exchangers. Int. J. Numer. Methods Fluids. 59, 677–690 (2009). Publisher Full Text

Chinyoka, T: Modelling of crossflow heat exchangers with viscoelastic fluids. Nonlinear Anal., Real World Appl.. 10, 3353–3359 (2009). Publisher Full Text

Qi, H, Wang, LH, Zheng, WY: On lockingfree finite element schemes for threedimensional elasticity. J. Comput. Math.. 23, 101–112 (2005)

Basak, T, Roy, S, Balakrishnan, AR: Effects of thermal boundary conditions on natural convection flows within a square cavity. Int. J. Heat Mass Transf.. 49, 4525–4535 (2006). Publisher Full Text