Next Article in Journal
Constitutive Explanations as a Methodological Framework for Integrating Thermodynamics and Economics
Next Article in Special Issue
Natural Convection and Entropy Generation in Nanofluid Filled Entrapped Trapezoidal Cavities under the Influence of Magnetic Field
Previous Article in Journal
Towards the Development of a Universal Expression for the Configurational Entropy of Mixing
Previous Article in Special Issue
Effect of Magnetic Field on Entropy Generation in a Microchannel Heat Sink with Offset Fan Shaped
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Entropy Generation in Natural Convection of Nanofluid inside a Square Cavity Having Hot Solid Block: Tiwari and Das’ Model

1
Department of Theoretical Mechanics, Faculty of Mechanics and Mathematics, Tomsk State University, 634050 Tomsk, Russia
2
Institute of Power Engineering, Tomsk Polytechnic University, 634050 Tomsk, Russia
3
Department of Mechanical Engineering, Technology Faculty, Fırat University, 23100 Elazig, Turkey
4
Department of Mechanical Engineering, Engineering Faculty, King Abdulaziz University, 21589 Jeddah, Saudi Arabia
5
Department of Mathematics, Babeş-Bolyai University, 400084 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
Submission received: 9 November 2015 / Revised: 17 December 2015 / Accepted: 21 December 2015 / Published: 31 December 2015
(This article belongs to the Special Issue Entropy in Nanofluids)

Abstract

:
A computational work has been performed in this study to investigate the effects of solid isothermal partition insertion in a nanofluid filled cavity that is cooled via corner isothermal cooler. Mathematical model formulated in dimensionless primitive variables has been solved by finite volume method. The study is performed for different geometrical ratio of solid inserted block and corner isothermal cooler, Rayleigh number and solid volume fraction parameter of nanoparticles. It is observed that an insertion of nanoparticles leads to enhancement of heat transfer and attenuation of convective flow inside the cavity.

1. Introduction

Solid blocks or bars can be used to control heat transfer as passive element in different shaped enclosures filled with nanofluids or pure fluids. This application can be observed in building design, electronic equipment, heat exchangers and some solar energy systems.
Buoyancy driven convection in a cavity filled with nanofluids with several pairs of heaters and coolers inside is analyzed computationally by Garoosi et al. [1]. They showed that the highest and the lowest impacts of design parameters, on the enhancement of heat transfer rate are caused by changing the position of heater and cooler and types of nanoparticles, respectively. Mahmoodi and Sebdani [2] solved the problem of free convection of Cu/water nanofluid in a cavity having adiabatic square bodies at its center with a numerical technique. They observed that for all Rayleigh numbers with the exception of Ra = 104, the average Nusselt number increases with increasing of nanoparticle volume fraction. Lee and Ha [3] considered a physical model with horizontal layer of fluid heated below and cold above with a conducting body placed at the center of the layer at different thermal conductivity ratio. They solved unsteady problem for different Rayleigh numbers. Multi-domain technique is used to handle a square-shaped conducting body. They presented the results for the case of dimensionless thermal conductivity of unity and also compared the results to those of pure natural convection. They also extended the similar geometry for different parameters [4]. Antar [5] calculated the heat transfer rate across building blocks in a cavity. He considered the conduction heat transfer in the block material and both natural convection and radiation in the cavity and found that increasing the number of cavities while keeping the block width unchanged decreases the heat flux significantly. Martyushev and Sheremet [6] numerically analyzed natural convection combined with thermal surface radiation in an enclosure with a heat-conducting solid walls of finite thickness and having a local heat source. They revealed that regardless of the considered solid–fluid interface, the average convective Nusselt number increases with the Rayleigh number and thermal conductivity ratio, and decreases with the surface emissivity and ratio of solid wall thickness to cavity spacing. The average radiative Nusselt number increases with the Rayleigh number, surface emissivity and thermal conductivity ratio and decreases with ratio of solid wall thickness to cavity spacing. Arefmanesh et al. [7] used the meshless local Petrov–Galerkin method to make simulation the buoyancy-driven flow and heat transfer in a differentially-heated enclosure having a baffle attached to its higher temperature side wall. They found that the average Nusselt number descends with increasing the baffle amplitude and the cold wall average Nusselt number increases as the baffle number of undulation augments. Bakkas et al. [8] made a computer solution on the two-dimensional laminar steady natural convection in a horizontal channel with the upper wall maintained at a constant temperature and the lower one provided with rectangular heating blocks, periodically distributed. In their work, the blocks are connected with adiabatic segments and their surfaces are assumed to release a uniform heat flux. They demonstrated that, depending on the length of the computational domain and the governing parameters, different flow structures can be obtained. Tsay et al. [9] made a numerical analysis to obtain characteristics of heat transfer from a block heat source module at different angles in two-dimensional cabinets. They observed that the maximum reduction in hot spot temperature is about 41% when two air vents are constructed on the cabinet wall. The variation of module angle results in the maximum difference of the hot spot temperature is 17% for closed cabinet, and 10% for ventilated cabinet. Deng and Tang [10] simulated the heatline in a square solid inserted cavity. They used finite volume and heatline formulation to make visualization. Visualization results by streamlines and heatlines directly exhibit the nature of fluid flow and heat transfer in macroscopical level, and hence, provide a more vigorous means to discuss the convective heat transfer accordingly. Shuja et al. [11] used two porous blocks located in a square cavity to see the effects of these blocks on heat and fluid flow. The flow conditions at the cavity inlet are kept the same for all the cases simulated by authors and found that the Nusselt number enhances with the increasing porosity and heat flux. Kuznetsov and Sheremet [12] studied conjugate natural convection in a rectangular enclosure having heat-conducting solid walls with a heat-generating solid block inside the cavity located on the bottom wall in conditions of nonuniform heat exchange with an environment. The authors found that an increase in the Grashof number leads to formation of the steady thermal plume, and also is reflected in reduction of the cooling degree of the cavity. It is revealed that variation of the heat-generating sizes and position of the heat source allows to control the fluid flow and heat transfer inside the cavity. Braga and Lemos [13] compared the heat transfer characteristics across a square cavity partially filled with a fixed amount of conducting solid material. The solid phase is shaped into two different geometries, namely square and cylindrical rods, which are horizontally displaced inside the cavity. When comparing the two geometries for the same Rayleigh number, the study showed that the average Nusselt number for cylindrical rods are slightly lower than of square rods.
Effects of inclination angle on natural convection flow in a square enclosure with a centered internal conducting square block are analyzed by Das and Reddy [14] by using finite volume technique. Ray and Chatterjee [15] made a work on solid object inserted cavity with corner heating and found similar results.
Calculation of entropy generation is an excellent tool to obtain the energy losses inside the system. The basic theory of entropy generation is widely presented by Bejan [16,17,18]. The application of entropy generation is also reviewed by Oztop and Al-Salem [19]. A numerical study on natural convection and entropy generation in a porous enclosure with heat sources is presented by Lami and Praka [20]. Again, combination of entropy generation with nanofluid filled cavity block insertion is given by Nayak et al. [21]. The entropy generation and Bejan number are evaluated to demonstrate the thermodynamic optimization of the mixed convection and they showed that the heat transfer rate increases remarkably by the addition of nanoparticles. Selimefendigil and Oztop [22] studied the natural convection and entropy generation of nanofluid filled cavity having different shaped obstacles with magnetic field effect. They observed that averaged heat transfer reduces by 21.35%, 32.85% and 34.64% for the cavity with circular, diamond and squared shaped obstacles, respectively, compared to cavity without obstacles at RaI = 106. It should be mentioned to this end the very good review paper on entropy generation in nanofluid flow by Mahian et al. [23].
The main aim of this study is to make a numerical analysis in a cavity filled with a nanofluid and a square heated object inserted. The cavity is cooled from the top wall and the left bottom corner. Thus, the study is a first attempt to control heat transfer via a passive element in a cavity with corner heater. Effects of entropy generation are studied to see the effective parameter on second law analysis of thermodynamics.

2. Basic Equations

Consider the free convection in a Cu-water nanofluid located in a square differentially heated cavity with hot centered block. A schematic geometry of the problem under investigation is shown in Figure 1, where x ¯ axis is measured in the horizontal direction along the bottom wall of the cavity and y ¯ axis is measured along the left vertical wall of the cavity and L is the size of the cavity. It is assumed that the top wall and left corner of the cavity are cooled with constant temperature T c , while the centered square block is heated with constant temperature T h . It is also assumed that other walls of the cavity are adiabatic. The walls of the cavity are assumed to be impermeable.
Figure 1. Physical model and coordinate system.
Figure 1. Physical model and coordinate system.
Entropy 18 00009 g001
Except for the density, the properties of the fluid are taken to be constant. It is further assumed that the effect of buoyancy is included through the Boussinesq approximation. Under the above assumptions, the conservation equations for mass, momentum and thermal energy can be written as follows
V = 0
ρ nf ( V , ) V = p + μ nf 2 V ( ρ β ) nf ( T T c ) g
( ρ C p ) nf ( V , T ) = k nf ( 2 T x ¯ 2 + 2 T y ¯ 2 )
The used physical properties of the nanofluid are described in detail earlier [24].
Equations (1)–(3) can be written in Cartesian coordinates as
u ¯ x ¯ + v ¯ y ¯ = 0
ρ nf ( u ¯ u ¯ x ¯ + v ¯ u ¯ y ¯ ) = p ¯ x ¯ + μ nf ( 2 u ¯ x ¯ 2 + 2 u ¯ y ¯ 2 )
ρ nf ( u ¯ v ¯ x ¯ + v ¯ v ¯ y ¯ ) = p ¯ y ¯ + μ nf ( 2 v ¯ x ¯ 2 + 2 v ¯ y ¯ 2 ) + ( ρ β ) nf g ( T T c )
T t + u ¯ T x ¯ + v ¯ T y ¯ = α nf ( 2 T x ¯ 2 + 2 T y ¯ 2 )
The existence of irreversibility source in the flow field such as the viscous dissipation effect and heat transfer causes the entropy generation. In Cartesian coordinates, the dimensional local entropy generation S ¯ g e n determined by Woods can be expressed as follows in the present study:
S ¯ g e n = k nf T 0 2 [ ( T x ¯ ) 2 + ( T y ¯ ) 2 ] + μ nf T 0 [ 2 ( u x ¯ ) 2 + 2 ( v y ¯ ) 2 + ( u y ¯ + v x ¯ ) 2 ]
where T 0 = 0.5 ( T c + T h ) .
Equation (8) constitutes two terms: the first is the local entropy generation due to the heat transfer ( S ¯ g e n , h t ) ; and the second is the dimensional local entropy generation due to fluid friction ( S ¯ g e n , f f ) .
Further, we introduce the following dimensionless variables
x = x ¯ / L ,   y = y ¯ / L ,   u = u ¯ / g β ( T h T c ) L ,   v = v ¯ / g β ( T h T c ) L , p = p ¯ / ρ n f g β ( T h T c ) L ,   θ = ( T T c ) / ( T h T c )
We are then left with the following equations
u x + v y = 0
u u x + v u y = p x + H 1 ( φ ) Pr R a ( 2 u x 2 + 2 u y 2 )
u v x + v v y = p y + H 1 ( φ ) Pr R a ( 2 v x 2 + 2 v y 2 ) + H 2 ( φ ) θ
u θ x + v θ y = H 3 ( φ ) R a Pr ( 2 θ x 2 + 2 θ y 2 )
The corresponding boundary conditions for these equations are given by
u = v = 0 on solid walls θ = 0 on cooled walls θ = 1 on heated walls of centered body θ / n = 0 on adiabatic walls
Here Pr = μ f ( ρ C p ) f / ( ρ f k f ) is the Prandtl number, R a = g ( ρ β ) f ( T h T c ) L 3 / ( α f μ f ) is the Rayleigh number, and the functions H 1 ( φ ) , H 2 ( φ ) and H 3 ( φ ) are given by
H 1 ( φ ) = 1 ( 1 φ ) 2.5 [ 1 φ + φ ρ p / ρ f ]
H 2 ( φ ) = 1 φ + φ ( ρ β ) p / ( ρ β ) f 1 φ + φ ρ p / ρ f
H 3 ( φ ) = k p + 2 k f 2 φ ( k f k p ) [ k p + 2 k f + φ ( k f k p ) ] [ 1 φ + φ ( ρ С p ) p / ( ρ С p ) f ]
These functions depend on the nanoparticles concentration ϕ and physical properties of the fluid and solid nanoparticles.
The dimensionless local entropy generation S g e n is obtained by using dimensionless parameters presented in Equation (9) given as:
S g e n = S ¯ g e n T 0 2 L 2 k f ( T h T c ) 2 = k nf k f [ ( θ x ) 2 + ( θ y ) 2 ] + + χ [ 2 ( u x ) 2 + 2 ( v y ) 2 + ( u y + v x ) 2 ] = S g e n , h t + S g e n , f f
In Equation (18), χ is the irreversibility factor. It is expressed as:
χ = μ nf T 0 k f g β L ( T h T c )
The integration of Equation (18) in the entire computational domain gives the dimensionless average entropy generation, S g e n , a v g , expressed as follows:
S g e n , a v g = 1 ϑ S g e n d ϑ = S g e n , h t , a v g + S g e n , f f , a v g
Further, the Bejan number Be is a parameter that shows the importance of heat transfer irreversibility in the domain and is defined as
B e = S g e n , h t S g e n , h t + S g e n , f f
The relative global dominance of heat transfer irreversibility is predicted by B e a v g (average Bejan number) which can be defined as
B e a v g = S g e n , h t , a v g S g e n , a v g
It may be noted that B e a v g > 0.5 shows that irreversibility due to heat transfer dominates in the flow. Fluid friction irreversibility dominates when B e a v g < 0.5 and when B e a v g = 0.5 , the heat transfer and fluid friction entropy generation are equal.
The physical quantities of interest are the local Nusselt number ( N u = k nf k f θ n ) along the heat source surface and the average Nusselt number, which is defined as
N u ¯ = k nf L 4 k f l ( 1 l / L 2 1 + l / L 2 θ y | y = 1 l / L 2 d x + 1 l / L 2 1 + l / L 2 θ x | x = 1 + l / L 2 d y + 1 l / L 2 1 + l / L 2 θ y | y = 1 + l / L 2 d x + 1 l / L 2 1 + l / L 2 θ x | x = 1 l / L 2 d y )

3. Numerical Method and Validation

The governing Equations (10)–(13) with corresponding boundary conditions Equation (14) have been solved on the basis of finite volume method [25,26]. The convective terms are discretized applying the power scheme and the diffusive terms are discretized with the central scheme. The SIMPLE algorithm is used to resolve the pressure–velocity coupling. The resulting algebraic equations systems for momentum equations were solved applying the line by line method with the implicit scheme of altering directions. The resulting algebraic equations system for energy equation was solved applying the explicit Buleev method for three-dimensional seven-point equations [27]. If values of mass balance for each control volume as well as the residual values of the different equations are sufficiently low overall convergence is obtained (typically 10−6). The above convergence criterion assures an acceptable solution [28,29]. The computer code was developed by authors using the C++ programming language.
The performance of part concerning natural convection inside the square cavity with a heated centered body was tested against the results of Asan [28]. Figure 2 and Figure 3 show a good agreement between the obtained streamlines and isotherms for different values of the Rayleigh number and solid block size and the numerical data of Asan [28]. The average Nusselt numbers of internal hot square block, which are also obtained by the present code, are compared with results of Asan [28] in Table 1.
Table 1. Comparison of present calculations with those of Asan [28].
Table 1. Comparison of present calculations with those of Asan [28].
AuthorsRa = 103Ra = 104Ra = 105Ra = 106
Asan [28]5.8135.9956.18810.425
Present results5.9135.9246.00610.598
Figure 2. Comparison of streamlines ψ and isotherms θ for Ra = 104, l/L = 0.3: numerical data of Asan [28] (a) and the present results (b).
Figure 2. Comparison of streamlines ψ and isotherms θ for Ra = 104, l/L = 0.3: numerical data of Asan [28] (a) and the present results (b).
Entropy 18 00009 g002
The performance of entropy generation part of the model was tested against the results of Ilis et al. [29] and Bhardwaj et al. [30] for steady-state natural convection in a differentially heated square cavity filled with the regular fluid for Prandtl number 0.7. Figure 4 and Figure 5 show a good agreement between the obtained fields of local entropy generation due to heat transfer and fluid friction with χ = 10 4 for different Rayleigh numbers and the numerical data of Ilis et al. [29] and Bhardwaj et al. [30].
Figure 3. Comparison of streamlines ψ and isotherms θ for Ra = 105, l/L = 0.6: numerical data of Asan [28] (a) and the present results (b).
Figure 3. Comparison of streamlines ψ and isotherms θ for Ra = 105, l/L = 0.6: numerical data of Asan [28] (a) and the present results (b).
Entropy 18 00009 g003
Figure 4. Comparison of local entropy generation due to heat transfer S g e n , h t and fluid friction S g e n , f f for Ra = 103: numerical data of Ilis et al. [29] (a); numerical data of Bhardwaj et al. [30] (b); and the present results (c).
Figure 4. Comparison of local entropy generation due to heat transfer S g e n , h t and fluid friction S g e n , f f for Ra = 103: numerical data of Ilis et al. [29] (a); numerical data of Bhardwaj et al. [30] (b); and the present results (c).
Entropy 18 00009 g004
For the purpose of obtaining grid independent solution, a grid sensitivity analysis is performed. The grid independent solution was performed by preparing the solution for steady free convection in a square cavity with centered heated filled with a Cu-water nanofluid at Ra = 105, Pr = 6.82, ϕ = 0.02, l / L = 0.3 , h / L = 0.5 . Three cases of the uniform grid are tested: a grid of 80 × 80 points, a grid of 100 × 100 points, a grid of 120 × 120 points, and a grid of 140 × 140 points. Figure 6 shows an effect of the mesh parameters on the temperature profiles at middle cross-section x = 0.5.
On the basis of the conducted verifications, the uniform grid of 120 × 120 points has been selected for the following analysis.
Figure 5. Comparison of local entropy generation due to heat transfer S g e n , h t and fluid friction S g e n , f f for Ra = 105: numerical data of Ilis et al. [29] (a); numerical data of Bhardwaj et al. [30] (b); and the present results (c).
Figure 5. Comparison of local entropy generation due to heat transfer S g e n , h t and fluid friction S g e n , f f for Ra = 105: numerical data of Ilis et al. [29] (a); numerical data of Bhardwaj et al. [30] (b); and the present results (c).
Entropy 18 00009 g005
Figure 6. Variation of the temperature profiles at x = 0.5 versus the mesh parameters.
Figure 6. Variation of the temperature profiles at x = 0.5 versus the mesh parameters.
Entropy 18 00009 g006

4. Results and Discussion

Numerical analysis has been conducted at the following values of the governing parameters: Rayleigh number (Ra = 103–106), solid volume fraction parameter of nanoparticles (φ = 0.0–0.05), dimensionless centered solid block size ( l / L = 0.2 0.6 ) and dimensionless cooler length in x and y directions ( h / L = 0.25 0.75 ) . Particular efforts have been focused on the effects of these governing parameters on the fluid flow and heat transfer inside the cavity. Streamlines, isotherms, entropy generation profiles, average Nusselt numbers, average total entropy generation and average Bejan number for different values of key parameters mentioned above are illustrated in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12.
Figure 7 presents distributions of streamlines, isotherms, local entropy generation due to heat transfer and fluid friction inside the cavity for different values of the Rayleigh number at l / L = 0.2 , h / L = 0.25 , ϕ = 0.03. Regardless of the Rayleigh number values, two convective cells are formed inside the cavity astride a heat source. It should be noted that one can find ascending flows close to the centered heater and descending flows near the vertical adiabatic wall. These two vortices characterize the presence of major motion from the left side of the heat source and minor circulation from the right side of the heater. An appearance of such division is due to the presence of the cooler in the left bottom corner that intensifies the left convective circulation: | ψ | max R a = 10 3 = 0.00173 < | ψ | max R a = 10 4 = 0.0042 < | ψ | max R a = 10 5 = 0.0116 . An increase in Ra leads to an intensification of convective flow with a displacement of the convective cores along the vertical direction. One can find that for Ra = 104 there are two convective cores inside the left vortex that can be explained by the effect of the cooling from both the top wall and the left bottom corner. At the same time, an increase in Ra leads to an intensification of the right convective cell with some deformation of this vortex in the bottom part from the main circulation.
Figure 7. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for l / L = 0.2 , h / L = 0.25 , ϕ = 0.03: Ra = 103 (a); Ra = 104 (b); Ra = 105 (c); and Ra = 106 (d).
Figure 7. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for l / L = 0.2 , h / L = 0.25 , ϕ = 0.03: Ra = 103 (a); Ra = 104 (b); Ra = 105 (c); and Ra = 106 (d).
Entropy 18 00009 g007
The temperature fields are also changed with Ra. For small values of the Rayleigh number (Ra = 103 and 104, Figure 7a,b), the main heat transfer regime is a heat conduction. This regime is characterized by isotherms that are parallel to the top horizontal isothermal wall. At the same time, one can find that heat fluxes from the top and corner coolers at Ra = 103 deform the heat flux from the centered heater. Therefore the main heating occurs in the right corner where the walls are adiabatic. An increase in Ra (= 104) leads to less intensive heating of the right corner of the cavity due to an appearance of convective heat flux and distributions of isotherms also along the vertical axis. The latter illustrates weak counteraction to the penetration of cold temperature from the top wall inside the cavity. In the case of Ra = 105 (Figure 7c) a thermal plume over the heat source is formed characterizing an amplification of the buoyancy force. Such development of the heat flux over the heat source prevents to the cooling of the cavity from the top wall, while cooling from the corner occurs in the bottom part of the cavity along the horizontal axis. It is worth noting that, for this value of the Rayleigh number, the thermal plume is a symmetric distribution of isotherms over the heat source. Further increase in Ra leads to an intensification of heat removal from the central heater along the vertical axis with asymmetric distributions of isotherms in this plume. At the same time, one can find more essential cooling of the bottom part of the cavity at Ra = 106 (Figure 7d).
For the analyzed problem, it is interesting to define the entropy generation areas inside the cavity. Figure 7 shows also the local entropy generation due to heat transfer Sgen,ht and fluid friction Sgen,ff. For the small values of the Rayleigh number (Ra = 103 and 104) distributions of Sgen,ht and Sgen,ff are uniform in general with the exception of the heater and cooler end points. For this Rayleigh number value the intensive entropy generation places due to heat transfer are the end points of the heater and corner cooler. Such distributions can be explained by the presence of some singularities in these points such as salient points for the heater or border points between isothermal and adiabatic conditions. At the same time, more intensive entropy generation due to fluid friction occurs close to the cooler and vertical walls of the heater that can be explained by the presence of essential velocity gradients in these parts of the cavity. An increase in the Rayleigh number leads to an intensification of the entropy generation inside the cavity. Formation of the thermal plume over the heater and a thin thermal boundary layer around the heater leads to an intensification of the entropy generation due to fluid friction in these zones. Moreover it is possible to conclude that for high values of the Rayleigh number distributions of local entropy generation are the similar to the distributions of isotherms inside the thermal plume.
Figure 8. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and Rayleigh number for l / L = 0.2 , h / L = 0.25 .
Figure 8. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and Rayleigh number for l / L = 0.2 , h / L = 0.25 .
Entropy 18 00009 g008
An intensification of the convective heat transfer rate around the heater Equation (22) with buoyancy force is presented in Figure 8. It should be noted that the average Bejan number reflects the ratio between the average entropy generation due to heat transfer and average total entropy generation. For low values of Ra when convective flow is characterized by small velocities and the entropy generation occurs due to heat transfer, therefore we have B e a v g 1 . An increase in Ra greater than 105 leads to an essential decrease in the average Bejan number (Figure 9b) with significant increase in the average total entropy generation (Figure 9c). Such behavior is due to an intensification of convective flow with an increase in the Rayleigh dissipative function. It should be noted that an increase in N u ¯ with Ra is similar qualitatively to an increase in S g e n , a v g with the Rayleigh number.
An influence of the nanoparticles volume fraction on the isolines of stream function, temperature, local entropy generation due to heat transfer and fluid friction is presented in Figure 9 for Ra = 105, l / L = 0.2 , h / L = 0.25 . An increase in ϕ leads to inessential changes in the local entropy generation while streamlines and isotherms have some changes. It should be noted that an increase in the concentration of nanoparticles leads to an attenuation of the convective flow taking into account the values of | ψ | max : | ψ | max φ = 0 = 0.0118 > | ψ | max φ = 0.05 = 0.0113 . At the same time, one can find less intensive cooling of the bottom part with ϕ taking into account the position of isotherm θ = 0.3.
Figure 9. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for Ra = 105, l / L = 0.2 , h / L = 0.25 : ϕ = 0.0 (a) and ϕ = 0.05 (b).
Figure 9. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for Ra = 105, l / L = 0.2 , h / L = 0.25 : ϕ = 0.0 (a) and ϕ = 0.05 (b).
Entropy 18 00009 g009
The effect of the concentration of nanoparticles on the average Nusselt and Bejan numbers and total entropy generation is presented in Figure 8. It is worth noting that an increase in ϕ leads to an increase in all these parameters. More intensive increase in N u ¯ and S g e n , a v g with nanoparticles volume fraction occurs for high values of Ra. As has been mentioned above, for high values of Ra, B e a v g 0.5 , therefore fluid friction irreversibility dominates.
The effect of the dimensionless centered solid block size on the distributions of local parameters is demonstrated in Figure 10. An increase in l / L leads to essential changes in all parameters. It is worth noting that an increase in l / L leads to a decrease in the distance between the top cold wall and the hot surface of the centered block therefore we have a local Rayleigh–Benard problem for narrow channel. As a result such narrowing leads to a formation of several convective cells of low intensity in this nanofluid gap [28]. Moreover, all these recirculations illustrate an appearance of ascending and descending thermal plumes over the heat source while under the heat source one can find heat conduction as a dominating heat transfer regime.
Figure 10. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for Ra = 106, h / L = 0.25 , ϕ = 0.03: l / L = 0.4 (a) and l / L = 0.6 (b).
Figure 10. Streamlines ψ, isotherms θ, local entropy generation due to heat transfer Sgen,ht, local entropy generation due to fluid friction Sgen,ff for Ra = 106, h / L = 0.25 , ϕ = 0.03: l / L = 0.4 (a) and l / L = 0.6 (b).
Entropy 18 00009 g010
An increase in the dimensionless centered solid block size leads to an essential decrease in the average Nusselt number (Figure 11a) and an increase in the average Bejan number (Figure 11b) and average total entropy generation (Figure 11c). A reduction of N u ¯ can be explained by an attenuation of the convective flow and also by narrowing of the nanofluid cavity. An increase in B e a v g and S g e n , a v g with l/L is due to an increase in the average entropy generation due to heat transfer inside the cavity where one can find an interaction of boundary layers in the case of narrow nanofluid zones.
Figure 11. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and dimensionless centered solid block size for Ra = 106, h / L = 0.25 .
Figure 11. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and dimensionless centered solid block size for Ra = 106, h / L = 0.25 .
Entropy 18 00009 g011
Figure 12 illustrates an influence of the dimensionless cooler length in x and y directions on the heat transfer rate, average Bejan number and average total entropy generation. An increase in h / L leads to an enhancement of the heat transfer rate, average Bejan number and average total entropy generation due to an increase in the temperature gradient and as a result an increase in the convective velocity. The latter leads to an increase in the entropy generation due to heat transfer.
Figure 12. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and dimensionless cooler length for Ra = 105, l / L = 0.2 .
Figure 12. Variation of the average Nusselt number (a); the average Bejan number (b); and the average entropy generation due to heat transfer (c) versus the nanoparticle volume fraction and dimensionless cooler length for Ra = 105, l / L = 0.2 .
Entropy 18 00009 g012

5. Conclusions

The steady problem of entropy generation and natural convection of Cu/water nanofluid within a differentially heated cavity with a centered hot square block and corner cooler was formulated in dimensionless primitive variables using the nanofluid model proposed by Tiwari and Das and solved numerically on the basis of finite volume method of a second-order accuracy. Profiles of stream function, temperature, local entropy generation due to heat transfer and local entropy generation due to fluid friction in a wide range of governing parameters have been obtained. The main findings can be listed as:
(1)
An increase in the Rayleigh number leads to an intensification of convective flow and heat transfer, and a reduction of the average Bejan number due to an enhancement of fluid friction. An increase in N u ¯ with Ra is similar qualitatively to an increase in S g e n , a v g with the Rayleigh number.
(2)
An increase in the nanoparticles volume fraction leads to an attenuation of convective flow and less intensive cooling of the bottom part of the cavity. Average Nusselt and Bejan numbers and average total entropy generation are increasing functions of ϕ.
(3)
An increase in the dimensionless centered solid block size leads to a formation of several convective cells of low intensity in nanofluid gap between top cold wall and top surface of the heater. In addition, the heat transfer rate decreases with l / L while average Bejan number and average total entropy generation are an increasing functions of l / L .
(4)
All considered average parameters are increasing functions of the dimensionless cooler length due to a widening of the heat sink surface.

Acknowledgments

This work of Mikhail A. Sheremet was conducted as a government task of the Ministry of Education and Science of the Russian Federation, Project Number 13.1919.2014/K. The authors wish to express their very sincere thanks to the very competent Reviewers for the very good comments and suggestions.

Author Contributions

Mikhail A. Sheremet conducted the simulations for the paper and wrote some sections of the manuscript. Hakan F. Oztop, Ioan Pop and Nidal Abu-Hamdeh prepared the Introduction part, wrote some other sections and analyzed the obtained numerical data. All authors contributed to the conclusions and assisted with reviewing and revising the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Belocal Bejan number
Beavgaverage Bejan number
Cpspecific heat at constant pressure
ggravitational acceleration vector
H 1 ( φ ) ,   H 2 ( φ ) ,   H 3 ( φ ) special functions
kthermal conductivity
Lsize of the cavity
Nulocal Nusselt number
N u ¯ average Nusselt number
pdimensionless pressure
p ¯ dimensional pressure
PrPrandtl number
RaRayleigh number
S g e n ¯ dimensional local entropy generation
S g e n dimensionless local entropy generation
S g e n , a v g dimensionless average entropy generation
S g e n , f f dimensionless local entropy generation due to fluid friction
S g e n , f f , a v g dimensionless average entropy generation due to fluid friction
S g e n , h t dimensionless local entropy generation due to the heat transfer
S g e n , h t , a v g dimensionless average entropy generation due to heat transfer
Tdimensional nanofluid temperature
T c temperature of the cold walls
T h temperature of the heat source
Vvelocity vector
u, vdimensionless velocity components
u ¯ , v ¯ dimensional velocity components
x, ydimensionless Cartesian coordinates
x ¯ , y ¯ dimensional Cartesian coordinates
Greek symbols
αthermal diffusivity
βthermal expansion coefficient
θdimensionless temperature
ϑ volume occupied by the nanofluid
μdynamic viscosity
ρdensity
ρ C p heat capacitance
ρ β buoyancy coefficient
φuniform concentration of the nanoparticles
χirreversibility factor
ψdimensionless stream function
Subscript
ccold
ffluid
nfnanofluid

References

  1. Garoosi, F.; Bagheri, G.; Talebi, F. Numerical simulation of natural convection of nanofluids in a square cavity with several pairs of heaters and coolers (HACs) inside. Int. J. Heat Mass Transf. 2013, 67, 362–376. [Google Scholar] [CrossRef]
  2. Mahmoodi, M.; Sebdani, S.M. Natural convection in a square cavity containing a nanofluid and an adiabatic square block at the center. Superlattices Microstruct. 2012, 52, 261–275. [Google Scholar] [CrossRef]
  3. Lee, J.R.; Ha, M.Y. A numerical study of natural convection in a horizontal enclosure with a conducting body. Int. J. Heat Mass Transf. 2005, 48, 3308–3318. [Google Scholar] [CrossRef]
  4. Lee, J.R.; Ha, M.Y. Numerical simulation of natural convection in a horizontal enclosure with a heat-generating conducting body. Int. J. Heat Mass Transf. 2006, 49, 2684–2702. [Google Scholar] [CrossRef]
  5. Antar, M.A. Thermal radiation role in conjugate heat transfer across a multiple-cavity building block. Energy 2010, 35, 3508–3516. [Google Scholar] [CrossRef]
  6. Martyushev, S.G.; Sheremet, M.A. Conjugate natural convection combined with surface thermal radiation in an air filled cavity with internal heat source. Int. J. Ther. Sci. 2014, 76, 51–67. [Google Scholar] [CrossRef]
  7. Arefmanesh, A.; Najafi, M.; Musavi, S.H. Buoyancy-driven fluid flow and heat transfer in a square cavity with a wavy baffle—Meshless numerical analysis. Eng. Anal. Bound. Elem. 2013, 37, 366–382. [Google Scholar] [CrossRef]
  8. Bakkas, M.; Amahmid, A.; Hasnaoui, M. Numerical study of natural convection heat transfer in a horizontal channel provided with rectangular blocks releasing uniform heat flux and mounted on its lower wall. Energy Convers. Manag. 2008, 49, 2757–2766. [Google Scholar] [CrossRef]
  9. Tsay, Y.L.; Cheng, J.C.; Zhuang, Y.L. Natural convective characteristics of a heat source module at different angles in the closed and ventilated cabinets. Int. Commun. Heat Mass Transf. 2010, 37, 1048–1056. [Google Scholar] [CrossRef]
  10. Deng, Q.H.; Tang, G.F. Numerical visualization of mass and heat transport for conjugate natural convection/heat conduction by streamline and heatline. Int. J. Heat Mass Transf. 2002, 45, 2373–2385. [Google Scholar] [CrossRef]
  11. Shuja, S.Z.; Yilbas, B.S.; Kassas, M. Flow over porous blocks in a square cavity: Influence of heat flux and porosity on heat transfer rates. Int. J. Ther. Sci. 2009, 48, 1564–1573. [Google Scholar] [CrossRef]
  12. Kuznetsov, G.V.; Sheremet, M.A. Conjugate natural convection in an enclosure with a heat source of constant heat transfer rate. Int. J. Heat Mass Transf. 2011, 54, 260–268. [Google Scholar] [CrossRef]
  13. Braga, E.J.; Lemos, M.J.S. Laminar natural convection in cavities filled with circular and square rods. Int. Commun. Heat Mass Transf. 2005, 32, 1289–1297. [Google Scholar] [CrossRef]
  14. Das, M.K.; Reddy, K.S.K. Conjugate natural convection heat transfer in an inclined square cavity containing a conducting block. Int. J. Heat Mass Transf. 2006, 49, 4987–5000. [Google Scholar] [CrossRef]
  15. Ray, S.; Chatterjee, D. MHD mixed convection in a lid-driven cavity including heat conducting solid object and corner heaters with Joule heating. Numer. Heat Transf. A 2014, 66, 530–550. [Google Scholar] [CrossRef]
  16. Bejan, A. Second law analysis in heat transfer. Energy Int. J. 1980, 5, 721–732. [Google Scholar] [CrossRef]
  17. Bejan, A. Entropy Generation through Heat and Fluid Flow; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1982. [Google Scholar]
  18. Bejan, A. Entropy Generation Minimization: The Method of Thermodynamic Optimization of Finite-Size Systems and Finite-Time Processes; CRC Press: New York, NY, USA, 1995. [Google Scholar]
  19. Oztop, 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]
  20. Lami, P.A.K.; Praka, K.A. A numerical study on natural convection and entropy generation in a porous enclosure with heat sources. Int. J. Heat Mass Transf. 2014, 69, 390–407. [Google Scholar] [CrossRef]
  21. Nayak, R.K.; Bhattacharyya, S.; Pop, I. Numerical study on mixed convection and entropy generation of Cu–water nanofluid in a differentially heated skewed enclosure. Int. J. Heat Mass Transf. 2015, 85, 620–634. [Google Scholar] [CrossRef]
  22. Selimefendigil, F.; Öztop, H.F. Natural convection and entropy generation of nanofluid filled cavity having different shaped obstacles under the influence of magnetic field and internal heat generation. J. Taiwan Inst. Chem. Eng. 2015, 56, 42–56. [Google Scholar] [CrossRef]
  23. Mahian, O.; Kianifar, A.; Kleinstreuer, C.; Al-Nimr, M.A.; Pop, I.; Sahin, A.Z.; Wongwises, S. A review on entropy generation in nanofluid flow. Int. J. Heat Mass Transf. 2013, 65, 514–532. [Google Scholar] [CrossRef]
  24. Sheremet, M.A.; Grosan, T.; Pop, I. Free convection in a square cavity filled with a porous medium saturated by nanofluid using Tiwari and Das’ nanofluid model. Transp. Porous Media 2015, 106, 595–610. [Google Scholar] [CrossRef]
  25. Ferziger, J.H.; Peric, M. Computational Methods for Fluid Dynamics; Springer: Berlin, Germany, 2002. [Google Scholar]
  26. Versteeg, H.K.; Malalasekera, W. An Introduction to Computational Fluid Dynamics. The Finite Volume Method; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1995. [Google Scholar]
  27. Buleev, N.I. The method of incomplete factorization for solving two-dimensional and three-dimensional equations of the diffusion type. USSR Comput. Math. Math. Phys. 1970, 10, 307–310. [Google Scholar] [CrossRef]
  28. Asan, H. Natural convection in an annulus between two isothermal concentric square ducts. Int. Commun. Heat Mass Transf. 2000, 27, 367–376. [Google Scholar] [CrossRef]
  29. Ilis, G.G.; Mobedi, M.; Sunden, B. Effect of aspect ratio on entropy generation in a rectangular cavity with differentially heated vertical walls. Int. Commun. Heat Mass Transf. 2008, 35, 696–703. [Google Scholar] [CrossRef]
  30. Bhardwaj, S.; Dalal, A.; Pati, S. Influence of wavy wall and non-uniform heating on natural convection heat transfer and entropy generation inside porous complex enclosure. Energy 2015, 79, 467–481. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Sheremet, M.A.; Oztop, H.F.; Pop, I.; Abu-Hamdeh, N. Analysis of Entropy Generation in Natural Convection of Nanofluid inside a Square Cavity Having Hot Solid Block: Tiwari and Das’ Model. Entropy 2016, 18, 9. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010009

AMA Style

Sheremet MA, Oztop HF, Pop I, Abu-Hamdeh N. Analysis of Entropy Generation in Natural Convection of Nanofluid inside a Square Cavity Having Hot Solid Block: Tiwari and Das’ Model. Entropy. 2016; 18(1):9. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010009

Chicago/Turabian Style

Sheremet, Mikhail A., Hakan F. Oztop, Ioan Pop, and Nidal Abu-Hamdeh. 2016. "Analysis of Entropy Generation in Natural Convection of Nanofluid inside a Square Cavity Having Hot Solid Block: Tiwari and Das’ Model" Entropy 18, no. 1: 9. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010009

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