OPEN ACCESS
Numerical simulations for 2D unsteady, incompressible flow with heat transfer in a foursided liddriven rectangular domain are reported in the present study. For the foursided liddriven rectangular domain, the lower wall is moved to the left, the upper wall is moved to the right, while the right wall is moved upwards and the left wall is moved downwards. All four walls move with equal speed. Different constant temperatures are applied to the left and right moving walls, and thermal insulation is applied to the upper and bottom moving walls. The governing equations are discretized using the QUICK scheme of finite volume methods. The SIMPLE algorithm is adopted to compute the numerical solutions of the flow variables, u velocity, v  velocity, P , and q as well as local and average Nusselt numbers for 50 £ Re £ 1, 500 and Pr = 6.63 . Due to the force generated by moving fluid, the direction of moving walls and the Reynolds number affect fluid flow in the rectangular domain in addition, at different Reynolds numbers along the cold wall of the domain, the variation in average and local Nusselt numbers reveals that overall heat transfer increases isotherms showed that as Reynolds numbers increase, the horizontal temperature gradient near the vertical walls decreases, because of which heat transfer decreases.
heat transfer, isotherms, nusselt number, velocity, streamlines
Fluid flow and heat transfer in a foursided liddriven rectangular domain has been the subject of intensive research in recent years. This is due to its significant applications such as cooling of electronic devices, furnaces, heat exchangers, boiler tubes, cooling of cylinder heads in I.C. engines, heating and cooling of buildings, heating of electric irons, heat treatment of engineering components, quenching of ingots, freezing of foods, etc.
Extensive literature studies have focused on heat transfer and fluid flow in rectangular or square cavities. These studies fall into two categories. The first deals with horizontal top [110] or bottom [11] wall sliding liddriven two dimensional cavities, in which the top wall has a constant velocity [13] or oscillates [4, 5], and behaves similarly in three dimensional cavities [6, 7, 9, 10, 16]. Various boundary conditions are applied to other solid walls in such cases. The second one is concerned with sidedriven, differentially heated cavities in these cases, left or right vertical wall or both vertical walls move with a constant velocity in their planes [12, 13, 14] in these studies, to create a temperature gradient in the cavity usually the liddriven side and the one opposing are heated differentially.
Lorenzini et at. [12] have investigated the constructal design of rectangular fin intruded into mixed convective liddriven cavity flows. Aydn [13] numerically studied mechanisms of aiding and opposing forces in a shear and buoyancydriven cavity. The square cavity had one vertical hot wall moving upwards or downwards, the opposite cold wall fixed, and both horizontal walls adiabatic. Oztop and Dagtekin [14] examined mixed convection in a twosided, liddriven differentially heated square cavity.
Kuhlmann et al. [15] conducted a numerical and experimental study on steady flow in rectangular twosided liddriven cavities. They found that the basic twodimensional flow was not always unique. For low Reynolds numbers, it consists of two separate corotating vortices adjacent to the moving walls.
Blohm and Kuhlmann [16] studied experimentally incompressible fluid flow in a rectangular container driven by two facing side walls which move steadily in antiparallel for Reynolds numbers up to 1,200. Two rotating cylinders of large radii close the cavity tightly, and create the moving side walls. Beyond a first threshold, threedimensional cells bifurcate supercritically out of the basic flow state. When both side walls move at the same rate (driven symmetrically), oscillatory stability was tricritical.
Xu et al. [17] have investigated the unsteady flow with heat transfer adjacent to the finned side wall of a differentially heated cavity with conducting adiabatic fin. Basak et at. [18] have investigated the effects of thermal boundary conditions on natural convection flows within a square cavity. Ilis et al.[19] have studied the effect of aspect ratio on entropy generation in a rectangular cavity with differentially heated vertical walls. Study of combined free convection and surface ratiation in closed cavities partially heated from below studied by Singh and Singh [20]. Wahba [21] has investigated the multiplicity of states for twosided and foursided lid driven cavity flows.
The abovementioned literature survey pertinent to the present problem under consideration revealed that liddriven cavities have interesting applications in various fields. However, no studies in the literature considered the case of a foursided liddriven rectangular domain with fluid flow and heat transfer. Specifically, in the present paper discussing the foursided liddriven rectangular domain, the lower wall is moved to the left, the upper wall is moved to the right, while the right wall is moved upwards and the left wall is moved downwards. All four walls move with equal speed.
What motivated us is the enormous scope of applications of unsteady incompressible flow with heat transfer as discussed earlier. Literature survey also revealed that the problem of fluid flow and heat transfer in a foursided lid driven rectangular domain, along with slip wall, and temperature boundary conditions, has not been studied numerically. Furthermore, to investigate the importance of the applications enumerated earlier, there is a need to determine numerical solutions of the unknown flow variables to fulfil this requirement, we present numerical simulations of the problem of fluid flow and heat transfer in a foursided lid driven rectangular domain, along with slip wall and temperature boundary conditions, using the SIMPLE algorithm.
Our main target of this work is to numerically investigate fluid flow and heat transfer in a foursided liddriven rectangular domain. We are used the QUICK scheme of finite volume methods to discretize the governing equations.SIMPLE algorithm is adopted to compute the numerical solutions of the flow variables, velocity, velocity, ,and as well as local and average Nusselt numbers for$50 \leq \operatorname{Re} \leq 1,500$ and $\mathrm{Pr}=6.63$ at time $t=0.001 \mathrm{s}$ .
The summary of the layout of the current work is as follows: Section 2 describes mathematical formulation that includes the physical description of the problem, governing equations, and initial and boundary conditions. Section 3 describes the numerical solution of the governing equations along with validity of results obtain with the benchmark solutions. Section 4 discusses the numerical results. Section 5illustrates the conclusions of this study.
2.1 Physical description
Geometry of the problem considered in this work along with the boundary conditions is depicted in the Figure 1. A foursided liddriven rectangular domain around the point (1,0.5) in which laminar unsteady incompressible flow is considered. The lower wall is moved to the left, the upper wall is moved to the right, while the right wall is moved upwards and the left wall is moved downwards. All four walls move with equal speed. The vertical lids have different constant temperatures. The horizontal walls are adiabatic.The left wall considered as the hot wall and the right is as the cold wall.
Figure 1. Geometry of the foursided lid driven rectangular domain problem.
We are assuming that, the values of $T_{0}$ and $\Delta T$ are chosen suitably so that the temperature defined on the left wall is greater than that of right wall. We are assumed that, at all four corner points of the computational domain, velocity Components $(u, v)$ vanish. It may be noted here regarding specifying the boundary conditions for pressure, the convention followed is that either the pressure at boundary is given or velocity components normal to the boundary are specified.
2.2 Governing equations
The flow is assumed to be twodimensional, unsteady state,laminar, and the fluid is incompressible. The dimensionless forms of the governing equations are the continuity, and the $x$ and $y$ components of the NavierStokes and the energy equations, assuming negligible dissipation and constant thermophysical properties, as given below:
$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$ (1)
$\frac{\partial u}{\partial t}+u \frac{\partial u}{\partial x}+v \frac{\partial u}{\partial y}=\frac{\partial P}{\partial x}+\left(\frac{1}{R e}\right)\left(\frac{\partial^{2} u}{\partial x^{2}}+\frac{\partial^{2} u}{\partial y^{2}}\right)$ (2)
$\frac{\partial v}{\partial t}+u \frac{\partial v}{\partial x}+v \frac{\partial v}{\partial y}=\frac{\partial P}{\partial y}+\left(\frac{1}{R e}\right)\left(\frac{\partial^{2} v}{\partial x^{2}}+\frac{\partial^{2} v}{\partial y^{2}}\right)$ (3)
$\frac{\partial \theta}{\partial t}+u \frac{\partial \theta}{\partial x}+v \frac{\partial \theta}{\partial y}=\left(\frac{1}{P r}\right)\left(\frac{\partial^{2} \theta}{\partial x^{2}}+\frac{\partial^{2} \theta}{\partial y^{2}}\right)$ (4)
where $u, v, P, \theta, R e$, and $P r$ are the dimensionless velocity components in $x$ and $y$ directions, the dimensionless pressure, the dimensionless temperature, the Reynolds number, and the Prandtl number respectively. We define the following nondimensional variables
$\left.\begin{array}{c}{x=\frac{x^{\prime}}{L}, y=\frac{y^{\prime}}{L}, u=\frac{u^{\prime}}{V_{P}}, v=\frac{v^{\prime}}{V_{p}}} \\ {P=\frac{P^{\prime}}{\rho V_{P}^{2}}, \theta=\frac{TT_{c}}{T_{h}T_{c}}, R e=\frac{V_{P} L}{v}, \operatorname{Pr}=\frac{v}{\alpha}}\end{array}\right\}$ (5)
The initial, noslip and slip wall boundary conditions are given by:
for $t=0$
$u(x, y, 0)=0, v(x, y, 0)=0, \theta(x, y, 0)=295.55$ (6)
$\left.\begin{array}{llll}{\text { for } t>0,} & {} & {} & {} \\ {x=0,} & {0<y<1,} & {u=0,} & {v=1,} & {\theta_{h}=\theta_{0}+\Delta \theta} \\ {x=2,} & {0<y<1,} & {u=0,} & {v=1,} & {\theta_{c}=\theta_{0}\Delta \theta} \\ {y=0,} & {0<x<2,} & {u=1,} & {v=0,} & {\frac{\partial \theta}{\partial y}=0} \\ {y=1,} & {0<x<2, u=1,} & {v=0,} & {\frac{\partial \theta}{\partial y}=0}\end{array}\right\}$ (7)
Local Nusselt number is
$N u_{y}=\frac{1}{A} \frac{\left(\frac{\partial \theta}{\partial x}\right)_{w}}{\left(\theta_{h}\theta_{c}\right)}$ (8)
Integration of the local Nusselt number along the wall is used to calculate the average Nusselt number.
$\overline{N u}=\frac{1}{A} \int_{0}^{1} N u_{y} d x$ (9)
The stream function is calculated from the definition
$u=\frac{\partial \psi}{\partial y}, v=\frac{\partial \psi}{\partial x}$ (10)
It is taken $\psi=0$ at the solid boundaries.
We are used a finite volume method to solve the governing equations in the discretizing controlling equations, the second order Quick scheme is selected. SIMPLE algorithm is adopted in solving the conservation equations. Discretization of the governing equations using Quick scheme is being skipped here as it is available in the literature [22].
All the results obtained in this study converged to a maximum residual of $10^{8}$. Furthermore, we are considered three different grid systems:$(100 \times 50),(150 \times 75)$, and $(200 \times 100)$ to ensure gridindependent results and $(150 \times 75)$ grid points were considered to obtain accurate solution in the entire computation of this study. We are also noted that by increasing the grid density by 77%, from $(150 \times 75)$ to $(200 \times 100)$ , the relative change in the value of the average Nusselt number is less than 1% (0.43%), which confirm that computed results on the $(200 \times 100)$ grid are indeed gridindependent.
The validation of current simulation has been verified with Iwatsu et al. [2] and Oztop and Dagtekin [14]. There is a good agreement for the average Nusselt numbers in the current study when compared to those of [2] and [14] as shown in table 1.
Table 1. Comparison of the average Nusselt number obtained in the present study with Iwatsu et al. [2] and Oztop and Dagtekin [14].
number 
Iwatsu et al. [2] Nu 
Oztop and Dagtekin [14] Nu 
This study Nu 
100 
1.34 
1.33 
1.30 
400 
3.62 
3.60 
3.59 
1,000 
6.29 
6.21 
6.21 
Fluid flow and temperature fields in a foursided liddriven rectangular domain are examined. The numerical simulations are performed with the same Reynolds numbers on both sides of the domain. Simulations in the rectangular domain are performed for the Reynolds number range from 50 to 1,500.The Foursided liddriven rectangular domain is analysed according to the direction of the moving plate shown in Figure 1.
Figure 2. Streamlines (a  c) and isotherms (b  d) for: Re = 50 and Re = 100.
Streamlines (a  c) and isotherms (b  d) for $R e=$ 50 and 100 are shown in Figure 2. For$R e=$ 50 and $R e=$ 100 in Figure 2(a) and (c), we observe, there is a primary cell with two secondary cells: the primary cell is at the center of the domain, however, the cell is not quite centered on the symmetry lines. The two secondary cells are weak and formed near the moving walls on both sides. For Figure 2(a) and (c),
Figure 3. (a) Streamlines and (b) isotherms for: Re = 500.
Streamlines and isotherms for$R e=$ 500 shown in Figure 3.For $R e=$500 in Figure 3(a), we observe that there are two secondary cells only. Both weaker cells are near the moving left and right walls. We also observe that streamlines become dense near the moving left and right walls.
For Figure 3(a) and (c),$\left\psi_{e x t}\right=1.68 \times 10^{3}, x=1.01, y=0.51$
Figure 4. Streamlines (a  c) and isotherms (b  d) for: Re = 1,000 and Re = 1,500.
Streamlines (a  c) and isotherms (b  d) for $R e=$ 1,000 and 1,500 are shown in Figure 4. For $R e=$ 1,000 and $R e=$1,500 in Figure 4(a) and (c), we observe, there is a primary cell with two secondary cells: the primary cell is at the center of the domain, however, the cell is not quite centered on the symmetry lines. The two secondary cells are weak and formed near the moving walls on both sides.
For Figure 4(a) and (c),
From Figure 2(a) and (c), Figure 4(a) and (c), we observe,there is a primary cell with two secondary cells: the primary cell is at the center of the domain; however, the cell is not quite centered on the symmetry lines. The two secondary cells are weak and formed near the moving walls on both sides. From Figure 3(a), we observe that there are two secondary cells only. Both weaker cells are near the moving left and right walls. We also observe that streamlines become dense near the moving left and right walls.
This happens because of the fluid rises along the right cold wall and sinks on the left hot wall due to forces generated by the moving fluid. Formation of two rotating cells at each side and a rotating cell at the center makes the heat transfer from left to right possible. The same phenomena have not been observed with one and two vertical sided liddriven cavities in the literature isotherms in Figure 2(b), (d), Figure 3(b), and Figure 4(b) and (d) show that as the Reynolds numbers increase the horizontal temperature gradient near left and right vertical walls decreases, because of which heat transfer decreases.
Based on the numerical solutions for uvelocity, Figure 5 illustrates the variation of velocity along the vertical line through the geometric center of the rectangular domain at Reynolds numbers $R e=$ 50, 100, 500, 1,000, and 1,500. We can see that velocity increases from the bottom wall to the top wall of the rectangular domain.
Figure 5. velocity along the vertical line through the geometric center of the domain.
Figure 6. velocity along the horizontal line through the geometric center of the domain.
Based on the numerical solutions for vvelocity, Figure 6 illustrates the variation of velocity along the horizontal line through the geometric center of the rectangular domain at Reynolds numbers $R e=$ 50, 100, 500, 1,000, and 1,500. We can see that as $R e=$ is increased from 50 to 1,500 the  velocity profile looks more negative which is not the case for individual Reynolds numbers. We found that, velocity becomes almost similar in shape with increase Reynolds number. So, an oscillatory flow of velocity profiles has been observed in the domain.
Figure 7. Local Nusselt number along the cold wall.
Based on the numerical solutions for the local Nusselt number, Figure 7. illustrates the variation of the local Nusselt number along the cold wall of the rectangular domain at Reynolds numbers $R e=$ 50, 100, 500, 1,000, and 1,500. We observed that, for each Reynolds number, the local Nusselt number decreases and thereby there is an increase in decay of heat.
Figure 8. Average Nusselt number as a function of Reynolds number.
Overall heat transfer: Based on the numerical solutions for the average Nusselt number, Figure 8 illustrates the variation of the average Nusselt number as a function of Reynolds number. We can see that, as is increased from 50 to 1,500 the average Nusselt number is also increased. We observe that overall heat increases as increase of Reynolds numbers.
This paper presents a numerical study to examine the effects of moving walls on the fluid flow and heat transfer in a foursided liddriven rectangular domain insulated boundary conditions are imposed on horizontal side walls. The governing equations are solved using the SIMPLE finite volume method. We compared previously published work on special cases of the problem and found good agreement. We presented and discussed graphical results for various parametric conditions. Heat transfer mechanisms and flow characteristics inside a rectangular domain are strongly dependent on moving walls and the Reynolds number. As the Reynolds numbers increase, so do both the local and average Nusselt number.
Numerical solutions for velocity illustrate the variation of velocity along the vertical line through the geometric center of the rectangular domain at Reynolds umbers $R e=$ 50, 100, 500, 1,000, and 1,500. We observed that velocity increases with increasing values of Reynolds from the bottom wall to the top wall of the rectangular domain. The numerical solutions for velocity illustrate the variation of velocity along the horizontal line through the geometric center of the rectangular domain at Reynolds numbers $R e=$ 50, 100, 500,1,000, and 1,500. We can see that as is increased from 50 to 1,500 the velocity profile looks more negative which is not the case for individual Reynolds numbers. We found that,velocity becomes almost similar in shape with increase Reynolds number. So, an oscillatory flow of velocity profiles has been observed in the domain.
From streamlines plots for Reynolds numbers 50, 100,1,000, and 1,500 of the rectangular domain, we observe, there is a primary cell with two secondary cells: the primary cell is at the center of the domain; however, the cell is not quite centered on the symmetry lines. The two secondary cells are weak and formed near the moving walls on both sides. While for Reynolds number 500, we observed that there are two secondary cells only. Both weaker cells are near the moving left and right walls. We also observe that streamlines become dense near the moving left and right walls. This happens due to the rises of the fluid along the right cold wall and sinks on the left hot wall due to forces generated by the moving fluid.Formation of two rotating cells at each side and a rotating cell at the center makes the heat transfer from left to right possible.The same phenomena have not been observed with one and two verticalsided liddriven cavities in the literature isotherms showed that, as Reynolds numbers increase the horizontal temperature gradient near the vertical walls decreases, because of which heat transfer decreases.
Numerical solutions for the local Nusselt number illustrates the variation of local Nusselt numbers along the cold wall of the rectangular domain at Reynolds numbers 50, 100, 500,1,000, and 1,500. We observed that, for each Reynolds number, the local Nusselt number decreases and thereby there is an increase in decay of heat. The numerical solutions for the average Nusselt number illustrates the variation of average Nusselt numbers as a function of Reynolds number.We observe that the average Nusselt number, or overall heat,increases with increasing Reynolds numbers.
The first author acknowledge the support from research council, University of Delhi for providing research and development grand 2015–16 vide letter No. RC/20159677 to carry out this work.
Δt 
time spacing 
$\frac{\partial}{\partial n}$ 
differentiation along the normal to the boundary 
A 
aspect ratio 
p’ 
pressure, pa 
P 
dimensionless pressure 
Nu Nu 
local Nusselt number 
H 
cavity height, m 
L 
cavity height, m 
T 
temperature, K 
Vp 
liddriven plate velocity, m.s1 
x’, y’ 
cartesian coordinates 
x, y 
dimensionless Cartesian coordinates 
u’, v’ 
velocities components in x,y direction, 

m.s1 
u, v 
dimensionless velocities 
Re 
Reynolds number 
Pr 
Prandtl number 
Greek symbols 

a 
thermal diffusivity, m2.s1 
υ 
kinematic viscosity, m2.s1 
ψ 
stream function 
ρ 
fluid density, kg.m3 
µ θ 
dynamic viscosity, kg. m1.s1 dimensionless temperature 
Subscripts 

c 
cold wall 
h 
hot wall 
p w 
plate wall 
[1] Torrance K., Davis R., Eike K., Gill D., Gutman D.,Hsui A., Lyons S., Zien H. (1972). Cavity flows driven by buoyancy and shear, J. Fluid Mech, Vol. 51, No. 2,pp. 221231. DOI: 10.1017/S0022112072001181
[2] Iwatsu R., Hyun J.M., Kuwahara K. (1993). Mixed convection in a driven cavity with a Stable vertical temperature gradient, Int. J. Heat Mass Transfer, Vol.36, No. 6, pp. 16011608. DOI: 10.1016/S00179310(05)800699
[3] Ramanan N., Homsy G.M. (1994). Linear stability of liddriven cavity flow, Phys. Fluids, Vol. 6, No. 8, pp.26902701. DOI: 10.1063/1.868158
[4] Iwatsu R., Hyun J.M., Kuwahara K. (1992).Convection in a differentiallyheated square cavity with a torsionallyoscillating lid, Int. J. Heat Mass Transfer, Vol. 35, No. 5, pp. 10691076. DOI:10.1016/00179310(92)90167Q
[5] Iwatsu R., Hyun J.M., Kuwahara K. (1992). Numerical simulation of flows driven by a torsionally oscillating lid in a square cavity, J. Fluids Eng., Vol. 114, No. 2,pp. 143149. DOI: 10.1115/1.2910008
[6] Freitas C.F., Street R.L. (1988). Nonlinear transient phenomena in a complex recirculating flow: anumerical investigation, Int. J. Numer. Meth. Fluids,Vol. 8, No. 7, pp. 769802. DOI:10.1002/fld.1650080704
[7] Iwatsu R., Hyun J.M. (1995). Threedimensional drivencavity flows with a vertical temperature gradient, Int. J. Heat Mass Transfer, Vol. 38, No. 18,pp. 33193328. DOI: 10.1016/00179310(95)00080S
[8] Mohamad A.A., Viskanta R. (1991). Transient low Prandtl number fluid convection in a liddriven cavity,277Numer. Heat Transfer A, Vol. 19, No. 2, pp. 187205.DOI: 10.1080/10407789108944845
[9] Prasad A.K., Koseff J.R. (1996). Combined forced and natural convection heat transfer in a deep liddriven cavity flow, Int. J. Heat Fluid Flow, Vol. 17, No. 5, pp.460467. DOI: 10.1016/0142727X(96)000549
[10] Mohamad A.A., Viskanta R. (1995). Flow and heat transfer in a liddriven cavity with stably stratified fluid, Appl. Math. Model, Vol. 19, No. 8, pp. 465472.DOI: 10.1016/0307904X(95)00030N
[11] Chen C.J., NassariNeshat H., Ho K.S. (1981). Finiteanalytical numerical solution of heat transfer in twodimensional cavity flow, Numer. Heat Transfer, Vol. 4,No. 2, pp. 179197. DOI:10.1080/01495728108961786
[12] Lorenzini G., Machado B.S., Isoldi L.A., dos Santos E.D., Rocha L.A.O. (2016). Constructal design of rectangular fin intruded into mixed convective liddriven cavity flows, ASME Journal of Heat Transfer,Vol. 138, No. 10, pp. 1025011/10250112. DOI:10.1115/1.4033378
[13] O. Aydın (1999). Aiding and opposing mechanisms of mixed convection in a shear and buoyancydriven cavity, Int. Comm. Heat Mass Transfer, Vol. 26, No. 7,pp. 10191028. DOI: 10.1016/S07351933(99)000913
[14] Oztop H.F., Dagtekin I. (2004). Mixed convection in twosided liddriven square cavity, Int. J. Heat Mass Transfer, Vol. 47, No. 89, pp. 17611769. DOI:10.1016/j.ijheatmasstransfer.2003.10.016
[15] Kuhlmann H.C., Wanschura M., Rath H.J. (1997).Flow in twosided liddriven cavities: nonuniqueness,instabilities and cellular structures, J. Fluid Mech, Vol.336, pp. 267299. DOI: 10.1017/S0022112096004727
[16] Blohm C.H., Kuhlmann H.C. (2002). The twosided liddriven cavity: experiments on stationary and timedependent flows, J. Fluid Mech, Vol. 450, pp. 6795.DOI: 10.1017/S0022112001006267
[17] Basak T., Roy S., Balakrishnan A.R. (2006). Effects of thermal boundary conditions on natural convection flows within a square cavity, Int. J. Heat Mass Transfer, Vol. 49, No. 2324, pp. 45254535. DOI:10.1016/j.ijheatmasstransfer.2006.05.015
[18] Xu F., Patterson J.C., Lei C. (2008). An experimental study of the unsteady thermal flow around a thin fin on a sidewall of a differentially heated cavity, Int. J. Heat Fluid Flow, Vol. 29, No. 4, pp. 11391153. DOI:10.1016/j.ijheatfluidflow.2008.01.001
[19] Ilis G.G., Mobedi M., Sunden B. (2008). Effect of aspect ratio on entropy generation in a rectangular cavity with differentially heated vertical walls, Int. Comm. Heat Mass Transfer, Vol. 35, No. 6, pp. 696703. DOI: 10.1016/j.icheatmasstransfer.2008.02.002
[20] Singh S.N., Singh D.K. (2015). Study of combined free convection and surface ratiation in closed cavities partially heated from below, Int.J. Heat and Tech, Vol.33, No. 2, pp. 18. DOI: 10.18280/ijht.330201
[21] Wahba E.M. (2009). Multiplicity of states for twosided and foursided lid driven cavity flows,Computers & Fluids, Vol. 38, No. 2. pp. 247253,.DOI: 10.1016/j.compfluid.2008.02.001
[22] Farhangnia M., Biringen S., Peltier L.J. (1996).Numerical simulation of twodimensional buoyancydriven turbulence in a tall rectangular cavity, Int. J. Numer. Methods Fluids, Vol. 23, No. 12, pp. 13111326. DOI: 10.1002/(SICI)10970363(19961230)23:12<1311::AIDFLD447>3.0.CO;28