Next Article in Journal
An Entropy-Based Neighborhood Rough Set and PSO-SVRM Model for Fatigue Life Prediction of Titanium Alloy Welded Joints
Next Article in Special Issue
Second Law Analysis of Unsteady MHD Viscous Flow over a Horizontal Stretching Sheet Heated Non-Uniformly in the Presence of Ohmic Heating: Utilization of Gear-Generalized Differential Quadrature Method
Previous Article in Journal
A Method for Detecting Dynamic Mutation of Complex Systems Using Improved Information Entropy
Previous Article in Special Issue
Doubly Nonnegative and Semidefinite Relaxations for the Densest k-Subgraph Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural Convection and Irreversibility Evaluation in a Cubic Cavity with Partial Opening in Both Top and Bottom Sides

by
Hakan F. Oztop
1,
Mohammed A. Almeshaal
2,
Lioua Kolsi
3,4,
Mohammed Mehdi Rashidi
5,6,* and
Mohamed E. Ali
7
1
Department of Mech. Engineering, Technology Faculty, Firat University, Elazig TR-23119, Turkey
2
Department of Mechanical Engineering, College of Engineering, Al Imam Mohammad Ibn Saud Islamic University, Riyadh 11432, Saudi Arabia
3
Department of Mech. Engineering, College of Engineering, Hail University, Hail 2240, Saudi Arabia
4
Laboratoire de Métrologie et des Systèmes Énergétiques, École Nationale d’Ingénieurs, University of Monastir, Monastir 5000, Tunisia
5
Shanghai Key Laboratory of Vehicle Aerodynamics and Vehicle Thermal Management Systems, Tongji University, Shanghai 200072, China
6
ENN-Tongji Clean Energy Institute of Advanced Studies, Tongji University, Shanghai 200072, China
7
Department of Mechanical Engineering, College of Engineering, King Saud University, Riyadh 11421, Saudi Arabia
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(2), 116; https://doi.org/10.3390/e21020116
Submission received: 3 December 2018 / Revised: 11 January 2019 / Accepted: 21 January 2019 / Published: 27 January 2019

Abstract

:
A numerical study on natural convection in a cubical cavity with partial top and bottom openings is performed in this paper. One of the vertical walls of the cavity has higher temperature than that of the opposite one; the remaining walls are insulated perfectly. Three-dimensional simulations of governing equations have been performed using a finite volume technique. The results are presented for different parameters such as opening length and Rayleigh number. It is observed that heat transfer rate and fluid flow can be controlled via opening ratio size and Rayleigh number.

1. Introduction

Convective heat transfer represents an important issue for various practical and engineering fields such as electronic cooling, solar heaters or building design. To make an efficient design, heat transfer and fluid flow are highly important parameters for these kinds of thermal systems.
Abib and Jaluria [1] computationally studied the free convection in cavities with partial openings. They solved the governing equation using the stream-function vorticity formulation based on the Boussinesq approximation. They noted that an increase of Ra causes a decrease in size of the recirculation cell that moves toward the vertical wall. Based on the Finite Volume Method (FVM), Polat and Bilgen [2] solved governing equation of natural convection in an opened tilted shallow cavitiy. Bilgen and Oztop [3] studied numerically the natural convection in an inclined cavity having partial opening. They observed that the inclination angle can be a means of control of heat transfer and flow structure inside the cavity. Bondareva et al. [4] used the heat-line visualization to analyze the flow structure and heat transfer in an open cavity filled with nanofluid and equipped with thick walls. Specifically, they solved a conjugate problem. Koufi et al. [5] studied numerically the mixed convection in open cavities. The authors considered the turbulent regime and tested the effects of the openings. Malekshah and Salari [6] performed both numerical and experimental analyses of free convection in cuboid geometries containing two immiscible fluids. Singh and Singh [7] investigated the combined radiation natural convection in inclined open cavities using the FVM. They concluded that the inclination is strongly determinant on heat transfer. Based on the Lattice Boltzmann Method (LBM), Sheikholeslami [8] investigated the hydrothermal behavior of the Magnetohydrodynamic (MHD) convection of nanofluid in an open cavity filled with porous media. He concluded that the rate of fluid exiting the cavity through the openings increases with Da and Ra. Bondareva et al. [9] used heat-line visualization to study the MHD flow in a wavy open cavity. Oztop et al. [10] considered three-dimensional partially open enclosures and performed a numerical investigation on three-dimensional natural convection with an evaluation of the different kinds of entropy generations.
Oztop et al. [11] studied the 2D natural convection in partially open cavities filled with porous media. They found that the increase of Grashof number increases the rate of heat transfer. Hinojosa et al. [12] analyzed the irreversibility generated by the natural convective flow and surface thermal radiation in a square open cavity. The surface thermal radiation was found to cause an important increase of the overall entropy generation. Hussain and Mustafa [13] studied the natural convection of a nanofluid in a locally heated parallelogrammic cavity having openings in its walls.
Singh and Singh [14] made a numerical solution on 2D natural convection in cavities with opened top wall considering the surface radiation. The effect of the volumetric heat generating source location on temperature field was studied. Correlations allowing the estimation of the maximum dimensionless temperature were proposed.
Bilgen and Muftuoglu [15] imposed a uniform heat flux on bouandries of open cavities. They found that heat transfer and the volume flow rates increases with Ra. Also, other studies related with open cavities and natural convection can be found in the literature, e.g., Prakash et al. [16], Gonzalez et al. [17] and Mohamad et al. [18]. Entropy generation can be calculated from the obtained data of temperature and velocities. It gives an opportunity to obtain energy losses inside the system. A numerical analysis of entropy production in an open cavity with pulsating flow in a horizontal channel has been presented by Zamzari et al. [19]. In a similar work, Mehrez et al. [20] performed work on entropy generation of nanofluids flow in an open cavity. Other works related to the subject can be found in [21,22,23,24,25,26,27,28,29,30,31,32,33].
In this work, a computational analysis on heat transfer, entropy generation and fluid flow due to natural convection in a cavity with partial opening from top and bottom sides is performed. Based on authors’ knowledge and the above literature, this work is a first step for the understanding of such 3D configuration.

2. Physical Model

The studied configuration is presented with coordinates in Figure 1. The model consists of a cubical cavity having two openings, i.e., one each at the bottom and ceiling. The right and left walls are differentially heated, and all remaining walls are considered adiabatic. The fluid (air) is considered incompressible, and the Boussinesq approximation is considered.

3. Mathematical Formulation

The 3D formalism ( ψ ω ) is used to simplify the treatment and the resolution of the equations governing the studied configuration by suppressing the pressure gradients terms. The vector potential ( ψ ) and vorticity ( ω ) in 3D geometries are expressed by:
ω = × V   and   V = × ψ
In the dimensionless form, equations governing the phenomenon are:
ω = 2 ψ
ω t + ( V × ) ω ( ω × ) V = Δ ω + R a × P r × [ T z ; 0 ; T x ]
T t + V × T = Δ T
The dimensionless numbers in the above equations can be expressed as:
P r = ν α   and   R a = g × β × Δ T × l 3 ν × α
Boundary conditions
The considered boundary conditions are as follows:
Temperature:
T = 1   at   x = 0
T = 0   at   x = 1
T n = 0   on   all   other   walls
At   open   boundary :   T i n = T c   if   n × V < 0   and   T n | o u t = 0   if   n × V 0
Velocities
V x = V y = V z = 0   on   all   walls
At   open   boundaries :   V x y = V y y = V z y = 0
Vorticities
ω x = 0 ,   ω y = V z x ,   ω z = V y x   at   x = 0   and   1
ω x = V z y , ω y = 0 ,   ω z = V x y   at   y = 0   and   1
ω x = V y z ,   ω y = V x z ,   ω z = 0   at   z = 0   and   1
Vector potential:
ψ x x = ψ y = ψ z = 0   at   x = 0   and   1
ψ x = ψ y y = ψ z = 0   at   y = 0   and   1
ψ x = ψ y = ψ z z = 0   at   z = 0   and   1
The generated entropy is expressed by:
S g e n = 1 T 2 × q × T + μ T × φ
With:
ϕ = 2 [ ( V x x ) 2 + ( V y y ) 2 + ( V z z ) 2 ] + ( V y x + V x y ) 2 + ( V z y + V y z ) 2 + ( V x z + V z x ) 2
Thus:
S g e n = k T 0 2 [ ( T x ) 2 + ( T y ) 2 + ( T z ) 2 ] + μ T 0 { 2 × [ ( V x x ) 2 + ( V y y ) 2 + ( V z z ) 2 ] + ( V y x + V x y ) 2 + ( V z y + V y z ) 2 + ( V x z + V z x ) 2 }
The dimensionless local entropy generation is expressed as:
N s = S g e n 1 k ( l T 0 Δ T ) 2
N s = [ ( T x ) 2 + ( T y ) 2 + ( T z ) 2 ] + ϕ × { 2 [ ( V x x ) 2 + ( V y y ) 2 + ( V z z ) 2 ] + [ ( V y x + V x y ) 2 + ( V z y + V y z ) 2 + ( V x z + V z x ) 2 ] }
with ϕ = μ α 2 T m l 2 k Δ T 2 is the irreversibility coefficient.
The total dimensionless entropy generation is:
S t o t = v N s d v = v ( N s t h + N s f r ) d v = S t h + S f r
Average Bejan number (Be) is evaluated using the following expression:
B e = S t h S t h + S f r + S J
Local and average Nusselt numbers are expressed respectively using:
N u = T x | x = 0   and   N u a v = 0 1 0 1 N u d y d z
The solutions of the above described governing equations were obtained using a code developed in the FORTRAN language. The FVM is used to develop the governing equations, the convective terms and temporal derivatives are discretized via the central-difference scheme and the fully implicit procedure, respectively. The solution is considered satisfactory if:
i 1 , 2 , 3 max | ψ i n ψ i n 1 | max | ψ i n | + max | T i n T i n 1 | 10 4

4. Verification and Grid Sensitive Study

As first verification the results of the present code have been compared with the 2D results of Bilgen and Oztop [3]. As shown in Figure 2, a good concordance in the flow structure and temperature field is encountered. A second verification based on 3D works of Wakashima and Saitho [34] for air filled cubic cavity is conducted and presented in Table 1. It can be concluded from the table that the code gives satisfactory results compared with those presented in the literature.
The grid sensitivity test has been performed for Pr = 0.7, Ra = 105 and d = 0.5. The average Nusselt number and the maximum of x-component velocity were used as testing parameters. Four spatial meshes of 613, 713, 813 and 913 were compared. The results of the grid dependency test are presented in Table 2. The increases from the grid 813 to 913 are 0.149% for Nuav and 0.734% for V x max . Thus, for computational economy and accuracy, spatial mesh size of 813 and time step of 10−4 are chosen to perform all simulations in the present work.

5. Results and Discussion

A numerical study is performed to investigate the natural convection heat transfer and fluid flow in a cavity with partially opened top and bottom sides. The effects of Ra and the dimension of the opening are highlighted.
Figure 3 illustrates the velocity vector projections and magnitude of velocity at z = 0.5 for Ra = 105 and different opening ratios. For d = 0.2, a huge vortex is located on the right side. The size of this clockwise circulating cell is reduced for d = 0.4. It becomes smaller and is pushed to right bottom side for d = 0.6. For d = 0.8 the cell. The flow goes from the bottom to the top opening, and lines are almost parallel to vertical walls with the very small circulating cell becoming very small. For d = 1, the circulation cell disappears and the flow passes directly from the bottom to the top of the cavity, except at the top right region, where the fluid enters and leaves from the top opening. It is interesting to note that flow comes into the cavity from the top side and dissipates suddenly for fully opened cases. In all cases, no symmetric structure was encountered.
In order to show the 3D structure of the flow, Figure 4 illustrates some particles trajectories for different opening widths at Ra = 105. The obtained results are compared with a fully opened case. As seen from the figure, the opening ratio directly effects the circulation cell dimension and location. For a fully opened case, no circulation cell is formed inside the cavity. This is an interesting result; flow does not have any circulation near the hot side for all cases due to a higher flow velocity. Due to the existence of circulation vortexes (for d < 1), the fluid passes transversally through the enclosure, which represents a purely 3D characteristic. For all presented cases, the circulation vortexes are convergent from the front and back walls to the center of the cavity.
As seen from the results of the temperature fields presented in Figure 5, iso-surfaces of temperature are pushed to the heated wall with increases in the opening ratio. This result is due to the increase of the flow velocity created by the presence of the openings. For d = 0.2, vertical stratification exists in the central region of the cavity due to the great size of the circulation vortex. By increasing the openings size, the vertical stratification disappears and the iso-surfaces of temperature become more piled near the hot wall due to the disappearance of the recirculation vortexes caused by the direct flow between the bottom and top openings.
Local entropy-generated contours at z = 0.5 plan are shown in Figure 6, Figure 7 and Figure 8 for different opening sizes and different Rayleigh numbers. For low Rayleigh numbers and small opening ratios, thermal entropy generation (Sth) is distributed across the entire domain. By increasing the Rayleigh number, Sth becomes concentrated near the hot wall, and by increasing d, it becomes concentrated near the bottom right corner due to the increase of temperature gradients in this region, as shown in Figure 5. In contrast to thermal entropy generation, frictional irreversibility becomes distributed across the entire cavity by increasing the opening ratio due to the enlargement of the domain from where the flow passes from the bottom to the top opening, causing an intensification of the flow, and consequently, an increase of the fluid-fluid friction. In all cases, the total entropy profile is similar to that of the thermal one, showing a dominance of Sth compared to Sfr, especially for low Rayleigh numbers, where Sth dominates in the domain. For higher Ra, the dominance is more pronounced on the right side of the cavity. This ascertainment is boosted by the Bejan number profile.
Figure 9a–e present the variations of average Nusselt number, thermal, viscous and total entropy generations and Bejan number with opening ratio. It is observed that higher heat transfer occurs for higher values of Ra, and that heat transfer increases almost linearly with opening ratio. This increment is clear for the lowest value of the Rayleigh number. Similarly, thermal entropy generation increases with both Rayleigh numbers and opening ratio almost linearly. But there is a huge increment for fully opened cavities.
The opening ratio becomes insignificant on viscous entropy generation for lower values of Rayleigh number. In contrast, a wavy distribution is formed for entropy generation due to fluid friction versus opening ratio. It is an interesting that there is a maximum value around d = 0.2. This result can serve for optimization purposes. The variation of total entropy generation shows a similar trend with thermal entropy generation except for high Ra and low opening sizes. The Bejan number is almost constant with opening ratio for lower Rayleigh numbers, and is around of 1, showing the dominance of thermal entropy generation. In contrast, a linear increase occurs for high Rayleigh numbers.

6. Conclusions

A numerical study has been performed to examine the heat transfer and fluid flow in a cubical cavity with top and bottom openings. The main findings can be listed as
  • Heat transfer enhances with increasing the opening ratios and Rayleigh numbers.
  • The flow field can be controlled in a cavity with the width of opening part. Both the dimension and location of circulation cells can be controlled.
  • The novelty of this work was to undertake an analysis of the natural cooling and natural ventilation problem of a model room. Thus, the results can be used for building ventilation.
  • As expected, the flow is more pronounced near of heated part compared to other regions.
  • The obtained results can be used for some kinds of filters or heating and cooling systems.
  • Increasing the opening part increases the entropy generation almost linearly for lower values of the Rayleigh number.

Author Contributions

Data curation, M.M.R.; Investigation, H.F.O.; Methodology, L.K.; Software, M.A.A.; Validation, M.E.A.

Funding

This research received no external funding.

Acknowledgments

First and last authors extend their appreciation to the International Scientific Partnership Program ISPP at King Saud University for funding this research work through ISPP#131.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CpSpecific heat at constant pressure (J/kg·K)
dDimensionless opening width
gGravitational acceleration (m/s2)
hDimensionless opening length, h’/l
kThermal conductivity (W/m·K)
lDimensionless cavity width
NuLocal Nusselt number
NuavAverage Nusselt number
PrPrandtl number
RaRayleigh number
tDimensionless time ( t × α / l 2 )
TDimensionless temperature [ ( T T c ) / ( T h T c ) ]
T c Cold temperature (K)
T h Hot temperature (K)
V Dimensionless velocity vector ( V · l / α )
x, y, zDimensionless Cartesian coordinates ( x / l , y / l , z / l )
Greek Symbols
α Thermal diffusivity (m2/s)
β Thermal expansion coefficient (1/K)
ρ Density (kg/m3)
μ Dynamic viscosity (kg/m·s)
ν Kinematic viscosity (m2/s)
ψ Dimensionless vector potential ( ψ / α )
ω Dimensionless vorticity ( ω × α / l 2 )
Δ T Dimensionless temperature difference
Subscripts
avAverage
x, y, zCartesian coordinates

References

  1. Abib, A.H.; Jaluria, Y. Numerical Simulation of the buoyancy-induced flow in a partially open enclosure. Numer. Heat Transf. 1988, 14, 235–254. [Google Scholar] [CrossRef]
  2. Polat, O.; Bilgen, E. Laminar natural convection in inclined open shallow cavities. Int. J. Therm. Sci. 2002, 41, 360–368. [Google Scholar] [CrossRef]
  3. Bilgen, E.; Oztop, H. Natural convection heat transfer in partially open inclined square cavities. Int. J. Heat Mass Transf. 2005, 48, 1470–1479. [Google Scholar] [CrossRef]
  4. Bondareva, N.S.; Sheremet, M.A.; Oztop, H.F.; Abu-Hamdeh, N. Heatline visualization of natural convection in a thick walled open cavity filled with a nanofluid. Int. J. Heat Mass Transf. 2017, 109, 175–186. [Google Scholar] [CrossRef]
  5. Koufi, L.; Younsi, Z.; Cherif, Y.; Naji, H. Numerical investigation of turbulent mixed convection in an open cavity: Effect of inlet and outlet openings. Int. J. Therm. Sci. 2017, 116, 103–117. [Google Scholar] [CrossRef]
  6. Malekshah, E.H.; Salari, M. Experimental and numerical investigation of natural convection in a rectangular cuboid filled by two immiscible fluids. Exp. Therm. Fluid Sci. 2017, 85, 388–398. [Google Scholar] [CrossRef]
  7. Singh, D.K.; Singh, S.N. Combined free convection and surface radiation in tilted open cavity. Int. J. Therm. Sci. 2016, 107, 111–120. [Google Scholar] [CrossRef]
  8. Sheikholeslami, M. Influence of magnetic field on nanofluid free convection in an open porous cavity by means of Lattice Boltzmann method. J. Mol. Liq. 2017, 234, 364–374. [Google Scholar] [CrossRef]
  9. Bondareva, N.S.; Sheremet, M.A.; Oztop, H.F.; Abu-Hamdeh, N. Heatline visualization of MHD natural convection in an inclined wavy open porous cavity filled with a nanofluid with a local heater. Int. J. Heat Mass Transf. 2016, 99, 872–881. [Google Scholar] [CrossRef]
  10. Oztop, H.F.; Kolsi, L.; Alghamdi, A.; Abu-Hamdeh, N.; Borjini, M.N.; Aissia, H.B. Numerical analysis of entropy generation due to natural convection in three-dimensional partially open enclosures. J. Taiwan Inst. Chem. Eng. 2017, 75, 131–140. [Google Scholar] [CrossRef]
  11. Oztop, H.F.; Al-Salem, K.; Varol, Y.; Pop, I. Natural convection heat transfer in a partially opened cavity filled with porous media. Int. J. Heat Mass Transf. 2011, 54, 2253–2261. [Google Scholar] [CrossRef]
  12. Hinojosa, J.F.; Buentello, D.; Xamán, J.; Tello, M.P. The effect of surface thermal radiation on entropy generation in an open cavity with natural convection. Int. Commun. Heat Mass Transf. 2017, 81, 164–174. [Google Scholar] [CrossRef]
  13. Hussain, A.K.; Mustafa, A.W. Natural convection in fully open parallelogrammic cavity filled with Cu–water nanofluid and heated locally from its bottom wall. Therm. Sci. Eng. Prog. 2017, 1, 66–77. [Google Scholar] [CrossRef]
  14. Singh, D.K.; Singh, S.N. Conjugate free convection with surface radiation in open top cavity. Int. J. Heat Mass Transf. 2015, 89, 444–453. [Google Scholar] [CrossRef]
  15. Bilgen, E.; Muftuoglu, A. Natural convection in an open square cavity with slots. Int. Commun. Heat Mass Transf. 2008, 35, 896–900. [Google Scholar] [CrossRef]
  16. Prakash, M.; Kedare, S.B.; Nayak, J.K. Numerical study of natural convection loss from open cavities. Int. J. Therm. Sci. 2012, 51, 23–30. [Google Scholar] [CrossRef]
  17. González, M.M.; Hinojosa, J.F.; Villafán-Vidales, H.I.; Bautista-Orozco, A.; Estrada, C.A. Theoretical and experimental study of natural convection with surface thermal radiation in a side open cavity. Appl. Therm. Eng. 2015, 75, 1176–1186. [Google Scholar] [CrossRef]
  18. Mohamad, A.A.; El-Ganaoui, M.; Bennacer, R. Lattice Boltzmann simulation of natural convection in an open ended cavity. Int. J. Therm. Sci. 2009, 48, 1870–1875. [Google Scholar] [CrossRef]
  19. Zamzari, F.; Mehrez, Z.; Cafsi, A.E.; Belghith, A.; Quéré, P.L. Numerical investigation of entropy generation and heat transfer of pulsating flow in a horizontal channel with an open cavity. J. Hydrodyn. 2017, 29, 632–646. [Google Scholar] [CrossRef]
  20. Mehrez, Z.; Bouterra, M.; Cafsi, A.E.; Belghith, A. Heat transfer and entropy generation analysis of nanofluids flow in an open cavity. Comput. Fluids 2015, 374, 214–224. [Google Scholar] [CrossRef]
  21. Rashidi, M.M.; Mohammadi, F.; Abbasbandy, S.; Alhuthali, M.S. Entropy generation analysis for stagnation point flow in a porous medium over a permeable stretching surface. J. Appl. Fluid Mech. 2015, 8, 753–765. [Google Scholar] [CrossRef]
  22. Rashidi, M.M.; Bhatti, M.M.; Abbas, M.A.; Ali, M.E.-S. Entropy generation on MHD blood flow of nanofluid due to peristalticwaves. Entropy 2016, 18, 117. [Google Scholar] [CrossRef]
  23. Abbas, M.A.; Bai, Y.; Rashidi, M.M.; Bhatti, M.M. Analysis of entropy generation in the flow of peristaltic nanofluids in channels with compliant walls. Entropy 2016, 18, 90. [Google Scholar] [CrossRef]
  24. Rashidi, M.M.; Nasiri, M.; Shadloo, M.S.; Yang, Z. Entropy Generation in a Circular Tube Heat Exchanger Using Nanofluids: Effects of Different Modeling Approaches. Heat Transf. Eng. 2017, 38, 853–866. [Google Scholar] [CrossRef]
  25. Bhatti, M.M.; Abbas, T.; Rashidi, M.M. Numerical study of entropy generation with nonlinear thermal radiation on magnetohydrodynamics non-newtonian nanofluid through a porous shrinking sheet. J. Magn. 2016, 21, 468–475. [Google Scholar] [CrossRef]
  26. Abbas, T.; Ayub, M.; Bhatti, M.M.; Rashidi, M.M.; Ali, M.-S. Entropy generation on nanofluid flow through a horizontal Riga plate. Entropy 2016, 18, 223. [Google Scholar] [CrossRef]
  27. Wang, Z.; Wei, Y.; Qian, Y. Numerical study on entropy generation in thermal convection with differentially discrete heat boundary conditions. Entropy 2018, 20, 351. [Google Scholar] [CrossRef]
  28. Wei, Y.; Wang, Z.; Qian, Y. A Numerical Study on Entropy Generation in Two-Dimensional Rayleigh-B´enard Convection at Different Prandtl Number. Entropy 2017, 19, 433. [Google Scholar] [CrossRef]
  29. Wei, Y.; Wang, Z.; Qian, Y.; Guo, W. Study on Bifurcation and Dual Solutions in Natural Convection in a Horizontal Annulus with Rotating Inner Cylinder Using Thermal Immersed Boundary-Lattice Boltzmann Method. Entropy 2018, 20, 733. [Google Scholar] [CrossRef]
  30. Wei, Y.; Yang, H.; Lin, Z.; Wang, Z.; Qian, Y. A novel two-dimensional coupled lattice Boltzmann model for thermal incompressible flows. Appl. Math. Comput. 2018, 339, 556–567. [Google Scholar] [CrossRef]
  31. Kumar, P.; Eswaran, V. A Numerical Simulation of Combined Radiation and Natural Convection in a Differential Heated Cubic Cavity. J. Heat Transf. 2009, 132, 023501. [Google Scholar] [CrossRef]
  32. Wang, P.; Zhang, Y.; Guo, Z. Numerical study of three-dimensional natural convection in a cubical cavity at high Rayleigh numbers. Int. J. Heat Mass Transf. 2017, 113, 217–228. [Google Scholar] [CrossRef] [Green Version]
  33. Wen, B.; Chini, G.P. Inclined porous medium convection at large Rayleigh number. J. Fluid Mech. 2018, 837, 670–702. [Google Scholar] [CrossRef]
  34. Wakashima, S.; Saitoh, T.S. Benchmark solutions for natural convection in a cubic cavity using the high-order time–space method. Int. J. Heat Mass Transf. 2003, 47, 853–864. [Google Scholar] [CrossRef]
Figure 1. Physical model with boundary conditions and coordinates.
Figure 1. Physical model with boundary conditions and coordinates.
Entropy 21 00116 g001
Figure 2. Comparison with results of Bilgen and Oztop [3].
Figure 2. Comparison with results of Bilgen and Oztop [3].
Entropy 21 00116 g002
Figure 3. Velocity vectors projections and magnitude of velocity at z = 0.5 plan for different opening widths and Ra = 105.
Figure 3. Velocity vectors projections and magnitude of velocity at z = 0.5 plan for different opening widths and Ra = 105.
Entropy 21 00116 g003aEntropy 21 00116 g003b
Figure 4. Particles trajectories for different opening widths and Ra = 105.
Figure 4. Particles trajectories for different opening widths and Ra = 105.
Entropy 21 00116 g004
Figure 5. Temperature field for different opening widths and Ra = 105.
Figure 5. Temperature field for different opening widths and Ra = 105.
Entropy 21 00116 g005
Figure 6. Local entropy generations at (z = 0.5) plan, for d = 0.1.
Figure 6. Local entropy generations at (z = 0.5) plan, for d = 0.1.
Entropy 21 00116 g006aEntropy 21 00116 g006b
Figure 7. Local entropy generations at (z = 0.5) plan, for d = 0.5.
Figure 7. Local entropy generations at (z = 0.5) plan, for d = 0.5.
Entropy 21 00116 g007aEntropy 21 00116 g007b
Figure 8. Local entropy generations at (z = 0.5) plan, for d = 0.9.
Figure 8. Local entropy generations at (z = 0.5) plan, for d = 0.9.
Entropy 21 00116 g008aEntropy 21 00116 g008b
Figure 9. Effect of opening size on (a) average Nusselt number, (b) Thermal entropy generation, (c) Viscous entropy generation, (d) Total entropy generation, (e) Bejan number.
Figure 9. Effect of opening size on (a) average Nusselt number, (b) Thermal entropy generation, (c) Viscous entropy generation, (d) Total entropy generation, (e) Bejan number.
Entropy 21 00116 g009aEntropy 21 00116 g009b
Table 1. Comparison with results of Wakashima and Saitho [34] (presented between parentheses).
Table 1. Comparison with results of Wakashima and Saitho [34] (presented between parentheses).
Ra ψ z
(center)
ω z
(center)
N u a v
1040.05528 (0.05492)1.1063 (1.1018)2.062 (2.062)
1050.034 (0.03403)0.262 (0.2573)4.378 (4.366)
1060.01972 (0.01976)0.1284 (0.1366)8.618 (8.6097)
Table 2. Grid sensitivity test for Pr = 0.7, Ra = 105, and d = 0.5.
Table 2. Grid sensitivity test for Pr = 0.7, Ra = 105, and d = 0.5.
Grid Nuavg V x max
6137.62135.442
7137.6571 137.565
8137.71246 141.871
9137.724 142.913

Share and Cite

MDPI and ACS Style

Oztop, H.F.; A. Almeshaal, M.; Kolsi, L.; Rashidi, M.M.; E. Ali, M. Natural Convection and Irreversibility Evaluation in a Cubic Cavity with Partial Opening in Both Top and Bottom Sides. Entropy 2019, 21, 116. https://doi.org/10.3390/e21020116

AMA Style

Oztop HF, A. Almeshaal M, Kolsi L, Rashidi MM, E. Ali M. Natural Convection and Irreversibility Evaluation in a Cubic Cavity with Partial Opening in Both Top and Bottom Sides. Entropy. 2019; 21(2):116. https://doi.org/10.3390/e21020116

Chicago/Turabian Style

Oztop, Hakan F., Mohammed A. Almeshaal, Lioua Kolsi, Mohammed Mehdi Rashidi, and Mohamed E. Ali. 2019. "Natural Convection and Irreversibility Evaluation in a Cubic Cavity with Partial Opening in Both Top and Bottom Sides" Entropy 21, no. 2: 116. https://doi.org/10.3390/e21020116

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop