Next Article in Journal
Constant Slope Maps and the Vere-Jones Classification
Previous Article in Journal
Product Design Time Forecasting by Kernel-Based Regression with Gaussian Distribution Weights
Previous Article in Special Issue
Entropy Generation on MHD Eyring–Powell Nanofluid through a Permeable Stretching Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Buoyancy-Induced Flow and Entropy Generation of Nanofluid-Filled Open Cavities Having Adiabatic Diamond Shaped Obstacles

1
College of Engineering, Mechanical Engineering Department, Haïl University, Haïl City 2240, Saudi Arabia
2
Unité de Recherche de Métrologie et des Systèmes Energétique, Ecole Nationale d’Ingénieurs, University of Monastir, Monastir City 5000, Tunisia
3
Young Researchers and Elite Club, Mashhad Branch, Islamic Azad University, Mashhad 9187147578, Iran
4
Department of Mechanical Engineering, Technology Faculty, Fırat University, Elazig 23119, Turkey
5
Mechanical Engineering Department, Faculty of Engineering, King Abdulaziz University, P.O. Box 80204, Jeddah 21589, Saudi Arabia
6
Unité de Recherche Matériaux, Énergie et Énergies Renouvelables, Faculté des Sciences de Gafsa, Gafsa 2100, Tunisia
*
Author to whom correspondence should be addressed.
Submission received: 1 May 2016 / Revised: 14 June 2016 / Accepted: 15 June 2016 / Published: 21 June 2016
(This article belongs to the Special Issue Entropy in Nanofluids)

Abstract

:
A three dimensional computational solution has been obtained to investigate the natural convection and entropy generation of nanofluid-filled open cavities with an adiabatic diamond shaped obstacle. In the model, the finite volume technique was used to solve the governing equations. Based on the configuration, the cavity is heated from the left vertical wall and the diamond shape was chosen as adiabatic. Effects of nanoparticle volume fraction, Rayleigh number (103 ≤ Ra ≤ 106) and width of diamond shape were studied as governing parameters. It was found that the geometry of the partition is a control parameter for heat and fluid flow inside the open enclosure.

Graphical Abstract

1. Introduction

The work on open systems is vital in the design of electronic cooling, heating and cooling of building, ventilation, food industry and solar applications. Bilgen and Öztop [1] solved via a finite volume method the two-dimensional governing equations of 2D natural convection in partially open inclined cavities. Prakash et al. [2] studied different heat transfer areas using Fluent. A buoyancy-induced flow loss is observed to increase with increasing opening ratio. This increment in buoyancy-induced flow loss for different inclination angles is found to range between 30%–80% when the opening ratio is increased from 1/4 to 1/2 for all closed space shapes. Zamora and Kaiser [3] solved the convective flows in cubical open cavities. They investigated the influence of an internal wall (adiabatic or isothermal) as well as the slope of an external heated wall on flow structure and heat transfer.
The number of studies on buoyancy-induced flow in nanofluid-filled open cavities, has increased recently. In this context, buoyancy-induced flow of a SWCNT-nanofluid (single-walled carbon nanotube nanofluids) in an open enclosure was studied by Jafari et al. using a lattice Boltzmann simulation [4]. It is mainly a square closed space heated with an isothermal vertical wall. They observed that the aspect ratio of the closed space plays an important role in buoyancy-induced flow and an increase of this parameter leads to heat transfer reduction in the target problem. Sheremet et al. [5] applied the curvilinear boundary conditions on a nanofluid-filled wavy walled closed space and they found that both waviness and nanoparticle volume rate are effective parameters in buoyancy-induced flow and heat transfer.
Using the lattice Boltzmann approach, Kefayati [6] studied the entropy generation due to laminar buoyancy induced flow in a square closed space filled with non-Newtonian nanofluid. He observed that total entropy production rises with the Ra number and the increase in volume fraction enhances entropy generation due to heat and fluid flow in different power-law indexes. Later, Kefayati [7] employed the lattice Boltzmann method to investigate the free convection in an open cavity in the presence of MHD flow using alumina/water nanofluid. Mejri and Mahmoudi [8] applied a magnetic field on buoyancy-induced flow heat transfer in a nanofluid-filled open enclosure.
Öztop and Salem [9] performed a review to present on entropy-related work with energy systems including buoyancy-induced flow. There are few works on three-dimensional analysis of entropy generation. Öztop et al. [10] numerically studied the three-dimensional heat transfer and entropy generation due to combined buoyancy and thermocapillary convection in a cubic cavity. Their results elucidated that at low Rayleigh numbers, the Marangoni number influence on the total entropy generation rises. Entropy generation for buoyancy-induced flow in a partitioned closed space, with adiabatic horizontal and isothermally cooled vertical walls, is studied numerically by Famouri and Hooman [11]. They presented the effects of fluid friction on entropy production is weak and the heat transfer irreversibility increases almost linearly with the Nu number and the temperature difference.
Mehrez et al. [12] exploited the finite-volume method to study the MHD flow effect on the entropy generation due to the use of Cu–water nanofluids in an open cavity where the bottom side was heated. They concluded that the average Nusselt number and entropy generation were enhanced with particle loading. Furthermore, it was found that the amount of increase in the Nusselt number and entropy generation depends on the intensity of the magnetic field (Hartmann number), and the angle at which the magnetic field is applied.
Heat transport can be controlled by installing differently shaped passive elements in open or closed cavities. In this context, conjugate heat transfer in a partially open square closed space with a vertical heat source has been simulated by Koca [13]. In this case, the closed space had an opening on the top side where the position and size of opening was changed. It was found that ventilation location has a significant influence on the heat transfer rate. Nasrin [14] concluded that obstacles could be used as effective tools to control the heat transfer in cavities where MHD flow is considered in the problem.
Hussein et al. [15] solved the three-dimensional buoyancy-induced flow and entropy generation in a trapezoidal closed space with different slopes filled with air for Rayleigh numbers in the range of 103 and 105. The closed space was heated from the bottom wall while the left and right sides were kept at a lower temperature compared to the bottom. They found that at low Rayleigh numbers, the slope of the closed space has no effect on the total entropy generation.
This literature review indicates that most studies on buoyancy-induced flow in open cavities have been considered as 2D dimensional problems. The main aim of this paper was to perform a 3D analysis of heat transport of nanofluids in an open cavity with a diamond shaped obstacle in the middle of a closed space.

2. Definition of the Physical Model

Figure 1 displays a schematic of the present problem. As shown, the left side of the closed space has a higher temperature than the right side that is considered as the opening. Other walls are assumed to be adiabatic. A diamond shaped obstacle with a side size of Ld and insulated walls is inserted in the middle of the closed space. The space between the closed space and the obstacle is filled with Al2O3/water nanofluid.
The nanofluid is assumed to be incompressible and the flow is considered laminar. The base fluid and nanoparticles are assumed to be at thermal equilibrium. Table 1 presents the thermo-physical properties of the base fluid (water) and the nanoparticles.

2.1. Governing Equations

Equations governing the studied configuration were written using the 3D vorticity-vector potential formalism ( ψ ω ) . This formalism allows the elimination of the pressure which is delicate to treat. Vorticity and vector potential are defined respectively by the following two relations:
ω = × V
and
V = × ψ
Using the dimensionless variables the governing equations are written in the following form:
ω = 2 ψ
ω t + ( V ) ω = ( ω ) V + ν n f ν f Pr 2 ω β n f β f R a Pr × T g
T t + V T = α n f α f 2 T
The velocity ( V ), time ( t ), vector potential ( ψ ) and vorticity ( ω ), are put respectively in the dimensionless forms by α / l , l 2 / α , α and l 2 / α and the dimensionless temperature is written as:
T = ( T T c ) / ( T h T c )
In these equations, the dimensionless Pr and Ra numbers are respectively defined as:
Pr = ν α  and Ra = g β Δ T l 3 ν α
The effective density, the specific heat, the effective thermal conductivity (Maxwell–Garnetts model) and the effective dynamic viscosity (Brinkman model) of the nanofluid are given respectively as [16,17]:
ρ n f = ( 1 ϕ ) ρ f + ϕ ρ s
( ρ C p ) n f = ( 1 ϕ ) ( ρ C p ) f + ϕ ( ρ C p ) s
k n f k f = k s + 2 k f 2 ϕ ( k f k s ) k s + 2 k f + ϕ ( k f k s )
μ n f = μ f ( 1 ϕ ) 2.5

2.2. Boundary Conditions

For the considered configuration, the boundary conditions are:
Temperature:
T = 1  for  x = 0
T n = 0  on adiabatic walls
T i n = T c  if  n . V < 0
and
T n | o u t = 0  if  n . V 0  At open boundary
Vorticity:
ω x = 0 ,   ω y = V z x  and  ω z = V y x  at  x = 0
ω x = V z y ,   ω y = 0  and  ω z = V x y  at  y = 0  and  1
ω x = V y z ,   ω y = V x z  and  ω z = 0  at  z = 0  and  1
Vector potential:
ψ x x = ψ y = ψ z = 0  at  x = 0
ψ x = ψ y y = ψ z = 0  at  y = 0  and  1
ψ x = ψ y = ψ z z = 0  at  z = 0  and  1
Velocity:
V x = V y = V z = 0  on all walls
V x x = V y x = V z x = 0  at open boundary
The expression of the dimensional generated entropy ( S g e n ) is presented by Equation (22):
S g e n = { k n f T 0 2 [ ( T x ) 2 + ( T y ) 2 + ( T z ) 2 ] } + μ n f 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 }
Put in its dimensionless form, the local generated entropy ( N s ) is written as:
N s = k n f k f [ ( T x ) 2 + ( T y ) 2 + ( T z ) 2 ] + ϕ S μ n f μ f { 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 ] }
In this equation another dimensionless parameter appears, this number is the irreversibility coefficient defined by:
ϕ S = ( α l Δ T ) 2 T 0
The local generated entropy is the sum of two types of entropies; the first one is due to temperature gradients (Ns–th) and the second is due to viscous effects (Ns-fr).
The total generated entropy ( S t o t ) is given by:
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
The Bejan number (Be) is the ratio of heat transfer irreversibility to the total irreversibility due to the heat transfer and the fluid friction as:
B e = S t h S t h + S f r
The local Nusselt number (Nu) is defined as follows:
N u = ( k n f k f ) T x | x = 0 , 1
The average Nusselt number on the hot wall ( N u a v ), is expressed by:
N u a v = 0 1 0 1 N u d y d z
The FORTRAN language was used to write the numerical program allowing to solve the mathematical model described above. Governing Equations (2)–(4) and (23) were discretized using the control volume method based on the central-difference scheme and solved using the successive relaxation iteration scheme. Although grids are uniform in all directions, additional nodes are added on boundaries. For all numerical tests, the time step and spatial mesh are fixed respectively at (10−4) and (813) and the convergence criterion for each step of time is chosen as:
i 1 , 2 , 3 max | ψ i n ψ i n 1 | max | ψ i n | + max | T i n T i n 1 | 10 5

3. Validation of the Code and Grid Dimension

For 3D analysis, the findings of Wakashima and Saitoh [18] and Fusegi et al. [19] on buoyancy- induced flow are validated with Pr = 0.71 and shown in Table 2. The present results are in excellent agreement with the previously reported results. After a grid sensivity test, the 81 × 81 × 81 grid distribution with additional nodes in boundaries was retained. For nanofluids buoyancy-induced flow, the code was validated by comparing with the studies of Öztop and Abu-Nada [20] (Figure 2) and of Jahanshahi et al. [21] (Figure 3). It can be said from this comparison that there is good agreement.

4. Results and Discussion

A computational work has been done to solve the equations of buoyancy-induced flow in a 3D open cavity with a diamond shaped obstacle. The working fluid is chosen as a nanofluid and entropy generation is also analyzed for different governing parameters such as nanoparticle volume fraction, Ra number and length of the diamond obstacle.
Streamlines and velocity vector projection at z = 0.5 plan for different Ra numbers and nanoparticle fractions are presented in Figure 4 and Figure 5, respectively, at Ld = 0.3. The flow enters the closed space from the bottom side and goes to the top side by surrounding the partition. Also, a small part of the fluid comes into the closed space and impinges on the inserted body and goes back. It shows small differences between nanofluid and pure fluid. Vena contracta shaped flow is formed around the inserted body. With increasing flow velocity and Ra number, a circulation cell is formed on the corner of the inserted body. This circulation cell becomes bigger with nanoparticle addition. Its dimension also increases with increasing Ra number.
Figure 6, presents the iso-surfaces of temperature for Ld = 0.3 at different Ra numbers for both nanofluid and water. For higher values of Ra number, the thermal boundary layer becomes thinner. Thus, insertion of the body becomes more effective for the lower values of Ra number. It is noticed that the inserted body has adiabatic boundary conditions. Thus, heat is captured under the left bottom side of the partition.
Figure 7 illustrates the isotherms for different Ra numbers and nanoparticle volume fraction at z = 0.5. As seen from the figure, distances among isotherms are almost same due to the domination of the conduction mode of heat transfer at low Ra numbers. Thus, isotherms are distributed diagonally due to the inserted insulated material. The effects of the inserted body become less due to thinner thermal boundary layer at higher Ra number. Local thermal entropy generation is displayed for different Ra number. In Figure 8, gray lines show values for pure fluid and black ones are given for nanofluids. Sth value becomes effective at the corners of the inserted body due to increases in kinetic energy. It is also cumulated mostly at the bottom side of the left vertical wall due to heat transport way. The inserted body behaves as an insulated wall. When values of Ra number increase the thermal boundary layer becomes thinner and the inserted body becomes insignificant on Sth value. Figure 9 illustrates the entropy generation due to fluid friction. Again, the edge of the inserted body becomes significant in entropy generation. Also, boundaries are effective for entropy generation. Especially, the top outlet edge becomes effective in entropy generation due to fluid friction, because friction is increased at that part of the flow.
Total entropy generation is presented in Figure 10 for the same parameters as Figure 8 and Figure 9. Figure 11 shows the variation of average Nusselt number along the hot wall as a function of nanoparticle volume fraction at different Ra numbers. For Ra = 103 and 104, heat transport becomes constant with volume fraction, due to the domination of conduction mode. For higher Ra, heat transfer increases by increasing the nanoparticle volume fraction.
Geometric effects of the inserted body are presented in Figure 12 for the highest value of Ra number. As seen from the figure, the dimension of the partition is the main effective parameter in the flow structure inside the closed space. In any case, a circulation cell occurs due to flow separation of the corner of the partition. Also, another mini-circulation cell is formed under this circulation cell. For the biggest partition, the flow is obstructed from the outside. Nanoparticle addition has less effect on flow structure.
In the same manner, isotherms are given in Figure 13 for different nanoparticle volume fraction at z = 0.5 for Ra = 106. The small partition size becomes insignificant in the temperature distribution. However, temperature distribution is affected by the partition and the high temperature penetrates into the cavity.
Average Nusselt number is illustrated for the considered parameters in Figure 14. As seen from the figure, the general trend with geometrical parameter is a small decrease. However, higher heat transfer is obtained with the addition of nanoparticles into the base fluid. Variation of Sth value with geometrical parameter is plotted in Figure 15 at different Ra numbers. As given in the figure, entropy generation becomes almost constant with geometrical parameters for lower Ra numbers. However, it is increased with the increase of geometrical parameter at higher values of the Ra number due to increase in heat transfer. Entropy generation due to friction as a function of geometrical parameter for different nanoparticle volume fraction is shown in Figure 16(a) for lower Ra numbers and Figure 16b for higher Ra numbers. It is an interesting result that heat transfer is decreased for Ld > 0.3 due to the obstruction effect of the inserted body. Entropy generation due to fluid friction becomes constant for different geometrical parameters at Ra = 105 for both nanofluid and pure fluid. It is increased with geometrical parameter for Ra = 106 and higher heat transfer occurs for nanofluids.
Figure 17 illustrates the total entropy generation for the chosen parameters. Total entropy generation becomes almost constant for Ra = 105 with Ld values and higher values are formed for nanofluids. However, entropy generation is increased almost linearly with Ld values for Ra = 106. On the contrary, entropy generation decreases with Ld values due to obstruction of the inserted body. Finally, Figure 18 presents the variation of Bejan number with Ld values at different Ra number and nanoparticle volume fraction. This value becomes almost constant for the highest value of Ra number for all values of geometrical parameters. On the contrary, Bejan number is decreased with increasing the geometrical parameter but in any cases, higher values are formed for the case of nanofluids. This decreasing trend is also valid for Ra = 104.

5. Conclusions

A computational study has been done in this work to obtain the heat transfer, fluid flow and entropy generation for different parameters such as nanofluids, geometrical parameter and Ra numbers. An inserted body can be used to control heat transfer, fluid flow and energy efficiency inside the closed space. Thus, Bejan number becomes almost fixed for the highest value of Ra number and it decreases with increasing geometrical value, but it decreases for lower values of Ra numbers. In any case both heat transfer and entropy generation augment with addition of nanoparticle. The increase of geometrical parameter enhances the total entropy generation due to the enhancement of entropy generation with heat transfer. Circulation cells inside the enclosure are a function of the geometrical parameters of the inserted body.

Author Contributions

In this paper, Hakan F. Öztop, Lioua Kolsi and Omid Mahian Conceived and designed the numerical experiments. Lioua Kolsi and Walid Aich performed the experiments. Nidal Abu-Hamdeh, Mohamed Naceur Borjini and Habib Ben Aissia analyzed the data. Hakan F. Öztop, Omid Mahian and Lioua Kolsi wrote the paper.

Conflicts of Interest

The authors declare that there is no conflict of interest.

Nomenclature

BeBejan number
CpSpecific heat at constant pressure (J/kg·K)
GGravitational acceleration (m/s2)
kThermal conductivity (W/m·K)
lEnclosure width and height (m)
LdDiamond width
NUnit vector normal to the wall
N s Dimensionless local generated entropy
NuLocal Nusselt number
PrPrandtl number
RaRayleigh number
S g e n Generated entropy (kJ/kg·K)
tDimensionless time ( t α / l 2 )
TDimensionless temperature [ ( T T c ) / ( T h T c ) ]
T c Cold temperature (K)
T h Hot temperature (K)
ToBulk temperature [To = ( T c+ T h)/2] (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)
ϕ Nanoparticles volume fraction
ϕ S Irreversibility coefficient
ψ Dimensionless vector potential ( ψ / α )
ω Dimensionless vorticity ( ω α / l 2 )
Δ T Dimensionless temperature difference

Subscripts

AvAverage
x, y, zCartesian coordinates
frFriction
fFluid
nfNanofluid
sSolid
thThermal
totTotal

Superscript

Dimensional variable

References

  1. Bilgen, E.; Öztop, H. Natural convection heat transfer in partially open inclined square cavities. Int. J. Heat Mass Transf. 2005, 48, 1470–1479. [Google Scholar] [CrossRef]
  2. Prakash, M.; Kedare, S.B.; Nayak, J.K. Numerical study of natural convection loss from open cavities. Int. J. Ther. Sci. 2012, 51, 23–30. [Google Scholar] [CrossRef]
  3. Zamora, B.; Kaiser, A.S. 3D effects in numerical simulations of convective flows in cubical open cavities. Int. J. Ther. Sci. 2014, 77, 172–185. [Google Scholar] [CrossRef]
  4. Jafari, M.; Farhadi, M.; Akbarzade, S.; Ebrahimi, M. Lattice Boltzmann simulation of natural convection heat transfer of SWCNT-nanofluid in an open enclosure. Ain. Shams Eng. J. 2015, 6, 913–927. [Google Scholar] [CrossRef]
  5. Sheremet, M.A.; Pop, I.; Shenoy, A. Unsteady free convection in a porous open wavy cavity filled with a nanofluid using Buongiorno’s mathematical model. Int. Commun. Heat Mass Transf. 2015, 67, 66–72. [Google Scholar] [CrossRef]
  6. Kefayati, G.H.R. FDLBM simulation of entropy generation due to natural convection in an enclosure filled with non-Newtonian nanofluid. Powder Technol. 2015, 273, 176–190. [Google Scholar] [CrossRef]
  7. Kefayati, G.H.R. Effect of a magnetic field on natural convection in an open cavity subjugated to water/alumina nanofluid using Lattice Boltzmann method. Int. Commun. Heat Mass Transf. 2013, 40, 67–77. [Google Scholar] [CrossRef]
  8. Mejri, I.; Mahmoudi, A. MHD natural convection in a nanofluid-filled open enclosure with a sinusoidal boundary condition. Chem. Eng. Res. Design 2015, 98, 1–16. [Google Scholar] [CrossRef]
  9. Öztop, H.F.; Al-Salem, K. A review on entropy generation in natural and mixed convection heat transfer for energy systems. Renew. Sustain. Energy Rev. 2012, 16, 911–920. [Google Scholar] [CrossRef]
  10. Öztop, H.F.; Kolsi, L.; Al-Salem, K. Numerical study of three-dimensional combined buoyancy and thermocapillary convection and evaluation of entropy generation. Int. J. Numer. Methods Heat Fluid Flow 2014, 24, 148–168. [Google Scholar] [CrossRef]
  11. Famouri, M.; Hooman, K. Entropy generation for natural convection by heated partitions in a cavity. Int. Commun. Heat Mass Transf. 2008, 35, 492–502. [Google Scholar] [CrossRef]
  12. Mehrez, Z.; El Gafsi, A.; Belghith, A.; Le Quéré, P. MHD effects on heat transfer and entropy generation of nanofluid flow in an open cavity. J. Magn. Magn. Mater. 2015, 374, 214–224. [Google Scholar] [CrossRef]
  13. Koca, A. Numerical analysis of conjugate heat transfer in a partially open square cavity with a vertical heat source. Int. Commun. Heat Mass Transf. 2008, 35, 1385–1395. [Google Scholar] [CrossRef]
  14. Nasrin, R. Finite element simulation of hydromagnetic convective flow in an obstructed cavity. Int. Commun. Heat Mass Transf. 2011, 38, 625–632. [Google Scholar] [CrossRef]
  15. Hussein, A.K.; Kolsi, L.; Chand, R.; Sivasankaran, S.; Nikbakhti, R.; Li, D.; Borjini, M.N.; Habib, B.A. Three-dimensional unsteady natural convection and entropy generation in an inclined cubical trapezoidal cavity with an isothermal bottom wall. Alexandria Eng. J. 2016, 55, 203–214. [Google Scholar] [CrossRef]
  16. Mahian, O.; Kianifar, A.; Kleinstreuer, C.; Moh’d, A.A.; Pop, I.; Sahin, A.Z.; Wongwises, S. A review of entropy generation in nanofluid flow. Int. J. Heat Mass Transf. 2013, 65, 514–532. [Google Scholar] [CrossRef]
  17. Costa, V.A.F.; Oliveira, M.S.A.; Sousa, A.C.M. Control of laminar natural convection in differentially heated square enclosures using solid inserts at the corners. Int. J. Heat Mass Transf. 2013, 46, 3529–3537. [Google Scholar] [CrossRef]
  18. 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. 2004, 47, 853–864. [Google Scholar] [CrossRef]
  19. Fusegi, T.; Hyun, J.M.; Kuwahara, K.; Farouk, B. A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure. Int. J. Heat Mass Transf. 1991, 34, 1543–1557. [Google Scholar] [CrossRef]
  20. Öztop, H.F.; Abu-Nada, E. Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  21. Jahanshahi, M.; Hosseinizadeh, S.F.; Alipanah, M.; Dehghani, A.; Vakilinejad, G.R. Numerical simulation of free convection based on experimental measured conductivity in a square cavity using Water/SiO2 nanofluid. Int. Commun. Heat Mass Transf. 2010, 37, 687–694. [Google Scholar] [CrossRef]
Figure 1. Physical model with boundary conditions.
Figure 1. Physical model with boundary conditions.
Entropy 18 00232 g001
Figure 2. Nusselt number versus ϕ for Ra = 104; comparison with results of Öztop and Abu-Nada [20].
Figure 2. Nusselt number versus ϕ for Ra = 104; comparison with results of Öztop and Abu-Nada [20].
Entropy 18 00232 g002
Figure 3. Comparison of the temperature on axial midline between the present results (cubic closed space; z = 0.5 plan) and results of Jahanshahi et al. [21] (Pr = 6.2, ϕ = 0.1 and Ra = 6.2 × 104).
Figure 3. Comparison of the temperature on axial midline between the present results (cubic closed space; z = 0.5 plan) and results of Jahanshahi et al. [21] (Pr = 6.2, ϕ = 0.1 and Ra = 6.2 × 104).
Entropy 18 00232 g003
Figure 4. Streamlines for Ld = 0.3 and different Ra; (Left) ϕ = 0; and (Right) ϕ = 0.15.
Figure 4. Streamlines for Ld = 0.3 and different Ra; (Left) ϕ = 0; and (Right) ϕ = 0.15.
Entropy 18 00232 g004
Figure 5. Velocity vector projection at z = 0.5 plan for Ld = 0.3 and different Ra; (Left) ϕ = 0; and (Right) ϕ = 0.15.
Figure 5. Velocity vector projection at z = 0.5 plan for Ld = 0.3 and different Ra; (Left) ϕ = 0; and (Right) ϕ = 0.15.
Entropy 18 00232 g005
Figure 6. Iso-surfaces of temperature for Ld = 0.3 and different Ra; Dark (ϕ = 0.15) clear (ϕ = 0).
Figure 6. Iso-surfaces of temperature for Ld = 0.3 and different Ra; Dark (ϕ = 0.15) clear (ϕ = 0).
Entropy 18 00232 g006
Figure 7. Isotherms at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Figure 7. Isotherms at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Entropy 18 00232 g007
Figure 8. Local thermal entropy generation (Sth) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Figure 8. Local thermal entropy generation (Sth) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Entropy 18 00232 g008
Figure 9. Local entropy generation due to friction (Sfr) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Figure 9. Local entropy generation due to friction (Sfr) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Entropy 18 00232 g009
Figure 10. Local total entropy generation (Stot) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Figure 10. Local total entropy generation (Stot) at z = 0.5 plan for Ld = 0.3 and different Ra; (gray) ϕ = 0; (black) ϕ = 0.15.
Entropy 18 00232 g010
Figure 11. Average Nusselt number at hot wall as function of ϕ for Ld = 0.3 and different Ra.
Figure 11. Average Nusselt number at hot wall as function of ϕ for Ld = 0.3 and different Ra.
Entropy 18 00232 g011
Figure 12. Velocity vector projection at z = 0.5 plan for Ra = 106 and different Ld = 0.3.
Figure 12. Velocity vector projection at z = 0.5 plan for Ra = 106 and different Ld = 0.3.
Entropy 18 00232 g012
Figure 13. Isotherms at z = 0.5 plan for Ra = 106 and different Ld (gray) ϕ = 0; (black) ϕ = 0.15.
Figure 13. Isotherms at z = 0.5 plan for Ra = 106 and different Ld (gray) ϕ = 0; (black) ϕ = 0.15.
Entropy 18 00232 g013
Figure 14. Average Nusselt number at hot wall for ϕ = 0 and ϕ = 0.15 and different Ra.
Figure 14. Average Nusselt number at hot wall for ϕ = 0 and ϕ = 0.15 and different Ra.
Entropy 18 00232 g014
Figure 15. Thermal entropy generation as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra.
Figure 15. Thermal entropy generation as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra.
Entropy 18 00232 g015
Figure 16. Entropy generation due to friction as a function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra; (a) Ra = 103 and 104; (b) Ra = 105 and 106.
Figure 16. Entropy generation due to friction as a function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra; (a) Ra = 103 and 104; (b) Ra = 105 and 106.
Entropy 18 00232 g016aEntropy 18 00232 g016b
Figure 17. Total entropy generation as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra; (a) Ra = 103 and 104; (b) Ra = 105 and 106.
Figure 17. Total entropy generation as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra; (a) Ra = 103 and 104; (b) Ra = 105 and 106.
Entropy 18 00232 g017
Figure 18. Bejan number as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra.
Figure 18. Bejan number as function of Ld for ϕ = 0 and ϕ = 0.15 and different Ra.
Entropy 18 00232 g018
Table 1. Thermo-physical properties of nanoparticles and base fluid.
Table 1. Thermo-physical properties of nanoparticles and base fluid.
PropertiesAl2O3Water
Cp (J·kg−1·°C−1)7654179
ρ (kg·m−3)3970997.1
k (W·m−1·°C−1)400.613
α × 107 (m2·s−1)131.71.47
β × 10−5 (K−1)0.8521
Table 2. Comparison of present results with the 3D results of (Wakashima and Saitho [18]) and (Fusegi et al. [19]) for differentially heated cubic closed space and Pr = 0.71.
Table 2. Comparison of present results with the 3D results of (Wakashima and Saitho [18]) and (Fusegi et al. [19]) for differentially heated cubic closed space and Pr = 0.71.
RaAuthors ψ z (Center) ω z (Center) V x max (y) V y max (x) N u a v
104Present work0.055281.10630.199 (0.826)0.221 (0.112)2.062
Wakashima and Saitho [18]0.054921.10180.198 (0.825)0.222 (0.117)2.062
Fusegi et al. [19]--0.201 (0.817)0.225 (0.117)2.1
105Present work0.0340.2620.143 (0.847)0.245 (0.064)4.378
Wakashima and Saitoh [18]0.034030.25730.147 (0.85)0.246 (0.068)4.366
Fusegi et al. [19]--0.147 (0.855)0.247 (0.065)4.361
106Present work0.019720.12840.0832 (0.847)0.254 (0.032)8.618
Wakashima and Saitho [18]0.019760.13660.0811 (0.86)0.2583 (0.032)8.6097
Fusegi et al. [19]--0.0841 (0.856)0.259 (0.033)8.77

Share and Cite

MDPI and ACS Style

Kolsi, L.; Mahian, O.; Öztop, H.F.; Aich, W.; Borjini, M.N.; Abu-Hamdeh, N.; Aissia, H.B. 3D Buoyancy-Induced Flow and Entropy Generation of Nanofluid-Filled Open Cavities Having Adiabatic Diamond Shaped Obstacles. Entropy 2016, 18, 232. https://0-doi-org.brum.beds.ac.uk/10.3390/e18060232

AMA Style

Kolsi L, Mahian O, Öztop HF, Aich W, Borjini MN, Abu-Hamdeh N, Aissia HB. 3D Buoyancy-Induced Flow and Entropy Generation of Nanofluid-Filled Open Cavities Having Adiabatic Diamond Shaped Obstacles. Entropy. 2016; 18(6):232. https://0-doi-org.brum.beds.ac.uk/10.3390/e18060232

Chicago/Turabian Style

Kolsi, Lioua, Omid Mahian, Hakan F. Öztop, Walid Aich, Mohamed Naceur Borjini, Nidal Abu-Hamdeh, and Habib Ben Aissia. 2016. "3D Buoyancy-Induced Flow and Entropy Generation of Nanofluid-Filled Open Cavities Having Adiabatic Diamond Shaped Obstacles" Entropy 18, no. 6: 232. https://0-doi-org.brum.beds.ac.uk/10.3390/e18060232

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