Next Article in Journal / Special Issue
Sensitivity Analysis of Entropy Generation in Nanofluid Flow inside a Channel by Response Surface Methodology
Previous Article in Journal
Classification Active Learning Based on Mutual Information
Previous Article in Special Issue
Natural Convection and Entropy Generation in Nanofluid Filled Entrapped Trapezoidal Cavities under the Influence of Magnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Generation and Natural Convection of CuO-Water Nanofluid in C-Shaped Cavity under Magnetic Field

1
Mechanical Engineering Department, Prince Mohammad Bin Fahd University, Al-Khobar 31952, Kingdom of Saudi Arabia
2
Mechanical Engineering Department, College of Engineering, University of Basrah, Basrah 61004, Iraq
3
Department of Engineering, Mahdishahr Branch, Islamic Azad University, Mahdishahr 75915-35618, Iran
*
Author to whom correspondence should be addressed.
Entropy 2016, 18(2), 50; https://doi.org/10.3390/e18020050
Submission received: 28 November 2015 / Revised: 10 January 2016 / Accepted: 27 January 2016 / Published: 5 February 2016
(This article belongs to the Special Issue Entropy in Nanofluids)

Abstract

:
This paper investigates the entropy generation and natural convection inside a C-shaped cavity filled with CuO-water nanofluid and subjected to a uniform magnetic field. The Brownian motion effect is considered in predicting the nanofluid properties. The governing equations are solved using the finite volume method with the SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm. The studied parameters are the Rayleigh number (1000 ≤ Ra ≤ 15,000), Hartman number (0 ≤ Ha ≤ 45), nanofluid volume fraction (0 ≤ φ ≤ 0.06), and the cavity aspect ratio (0.1 ≤ AR ≤ 0.7). The results have shown that the nanoparticles volume fraction enhances the natural convection but undesirably increases the entropy generation rate. It is also found that the applied magnetic field can suppress both the natural convection and the entropy generation rate, where for Ra = 1000 and φ = 0.04, the percentage reductions in total entropy generation decreases from 96.27% to 48.17% for Ha = 45 compared to zero magnetic field when the aspect ratio is increased from 0.1 to 0.7. The results of performance criterion have shown that the nanoparticles addition can be useful if a compromised magnetic field value represented by a Hartman number of 30 is applied.

1. Introduction

Nowadays, energy saving and environmental agencies strictly recommend energy conservation in order to save what remains of energy sources of the world and to reduce pollution. Hence, optimization of thermal systems has received unique attention. Recently, second-law-based investigations have gained attention for studying thermal systems. Entropy generation has been used as a gauge to evaluate the performance of a thermal system. Bejan [1,2,3] focused on the different reasons behind entropy generation in applied thermal engineering. Generation of entropy destroys available work of a system. Therefore it makes good engineering sense to focus on irreversibility of heat transfer and fluid friction process. Second-law analysis in ordinary fluids filling clear cavities can be summarized as follow. Mahmoud and Fraser [4] applied the first and second laws of thermodynamics to forced convection inside a cylindrical annular space with isoflux boundary conditions. Mahmud and Sadrul Islam [5] applied the second law of thermodynamics to predict the nature of irreversibility in inclined wavy wall cavity. Bouabid et al. [6] studied the transient state study of entropy generation in inclined cavity with various aspect ratios. Cheng and Liang [7] concluded that the applicability of the entropy generation minimization to optimize thermodynamic systems is conditional if the maximum output power is the objective.
Analyses of the second law of thermodynamics are also focused on saturated porous cavities, as in Baytas [8], Al-Hadhrami et al. [9], Mahmud and Fraser [10], Hooman and Ejlali [11], Zahmatkesh [12], Varol et al. [13], and Basak et al. [14].
The low thermal conductivity of conventional heat transfer fluids, common water, oil, and ethylene glycol, has restricted the improvement of thermal systems performance. Fluids containing nano-sized solid particles offer a possible solution to conquer this problem. The nanofluid has greater effective thermal conductivity than a pure base fluid. Nanofluids, a name conceived by Choi [15], are fluids consisting of solid nanoparticles with size less than 100 nm suspended with a specified solid volume fraction. Nanofluids can be used to improve thermal management systems in many engineering applications such as transportation, micromechanics instrument and cooling devices. Due to the encouraged enhanced properties associated with nanofluids properties, huge works have been published since about two decades ago (Aminossadati and Ghasemi [16], Ghasemi and Aminossadati [17], Abu-Nada and Chamkha [18], Nemati et al. [19], Mahmoudi et al. [20], Matin and Pop [21], Sheikhzadeh et al. [22], and Sheikholeslami et al. [23,24]).
Recently, the analysis of the second law of thermodynamics is directed at the handling of nanofluids, as in Kashani et al. [25] who investigated the entropy generation due to laminar natural convection of Cu-water nanofluid near the density maximum of water in a two-dimensional enclosure with various patterns of vertical wavy walls. Cho et al. [26] studied the heat transfer performance and entropy generation of natural convection in a nanofluid-filled U-shaped cavity. Ting et al. [27] studied the entropy generation of viscous dissipative nanofluid flow in thermal non-equilibrium porous media embedded in micro channels. Mahmoudi at al. [28] examined the natural convection in a square enclosure filled with a water-Al2O3 nanofluid in the presence of magnetic field and uniform heat generation/absorption. Aghaei et al. [29] evaluated the effects of magnetic field on heat transfer and entropy generation on mixed convection of nanofluids with variable properties in a trapezoidal enclosure. Kefayati [30] analyzed the heat transfer and entropy generation on laminar natural convection of non-Newtonian nanofluids in the presence of an external horizontal magnetic field in a square cavity. Ellahi et al. [31] analyzed a mathematical model in order to study the shape of nanosize particles on entropy generation and natural convection boundary layer flow along an inverted cone. More recently, Ismael et al. [32] studied the entropy generation due to conjugate natural convection conduction heat transfer in a square domain under steady-state condition. They proposed a new criterion for assessment of the thermal performance. As the applications of natural convection in engineering systems have developed, it became necessary to investigate non rectangular/square cavity shape. Cooling of electronic chips/system, micro-electro-mechanical systems (MEMS), solar collectors, and heat exchangers are some of the many important applications that handle different cavity shapes. Biserani et al. [33] used Bejan’s constructal theory to optimize the geometry of H-shaped cavity that intrudes into a solid conducted wall. They optimized other cavities namely, C-shaped and T-shaped cavities and found that H-shaped is superior in thermal performance. Mahmoodi [34] studied free convection in L-shaped cavity filled with Cu-water nanofluid. Mahmood and Hashemi [35] studied the C-shaped cavity filled with a nanofluid. Cho et al. [36] investigated the natural convection enhancement of Al2O3-water nanofluid in a U-shaped cavity. Mansour et al. [37] also investigated the natural convection inside U-shaped cavity filled with Cu-water nanofluid but they termed their cavity as C-shaped. Mojumder et al. [38] studied the natural convection in C-shaped cavity filled with Cobalt-kerosene ferrofluid under the effect of externally applied magnetic field. Kasaeipoor et al. [39] studied the convection of Cu-Water nanofluid in a vented T-shaped cavity in the presence of an externally applied magnetic field. Al-Zamily [40] investigated numerically the effect of constant magnetic field on natural convection in a semi-circular enclosure filled with Cu-water nanofluid with the present of heat flux.
The literature survey above led us to be sure that the entropy generation in C-shaped cavity filled with a nanofluid and subjected to a magnetic field has not been investigated yet. Therefore, this topic will be the material of the present paper. It is believed that this study will contribute to improving the thermal performance of electronic chip cooling systems.

2. Problem Description and Mathematical Modeling

The geometry of the present study is a C-shaped cavity filled with a CuO-water nanofluid and subjected to a horizontal uniform magnetic field (B0), which is shown in Figure 1. The two-dimensional cavity has equal length and height of L. The internal walls of cavity with the length H are maintained at a relatively low temperature Tc. The external right walls are insulated and the other walls are maintained at a relatively high temperature Th. The flow is assumed to be laminar, steady, and incompressible. The nanofluid is assumed to be Newtonian, incompressible and has a low electrical conductivity. Moreover, the water and CuO nanoparticles are in thermal equilibrium. The Joule heating are assumed to be negligible compared to the applied magnetic field.
Figure 1. Schematic of C-shaped cavity under magnetic field.
Figure 1. Schematic of C-shaped cavity under magnetic field.
Entropy 18 00050 g001
The equations that govern the conservation of mass, momentum and energy can be written as
u x + v y = 0 u u x + v u y = 1 ρ n f [ p x + μ n f ( 2 u x 2 + 2 u y 2 ) ] u v x + v v y = 1 ρ n f [ p y + μ n f ( 2 v x 2 + 2 v y 2 ) ] + ( ρ β ) n f g ( T T c ) σ n f   B 0 2 v u T x + v T y = α n f [ 2 T x + 2 T y 2 ]
where u, v and T are the velocity components in the x direction, y direction and temperature, respectively. The local entropy generation equation given in [1] can be adopted for nanofluid as:
S c = k n f T 0 2 [ ( T x ) 2 + ( T y ) 2 ] + μ n f T 0 [ 2 { ( u x ) 2 + ( v y ) 2 } + ( u y + v x ) 2 ]
The dimensionless form of the governing equations can be obtained by introducing the dimensionless variables. These are given as follows:
X = x L ,   Y = y L ,   U = u L α f ,   V = v L α f ,   A R = H L ,   P = p L 2 ρ n f α f 2 ,   θ = T   T c T h T c
R a = g   β f   L 3   ( T h T c ) α f   ν f ,   H a = B 0 L σ f μ f ,   P r = ν f α f
Based on these dimensionless variables, the continuity, momentum and energy equations for the laminar and steady-state mixed convection of nanofluid in the presence of magnetic field can be written in non-dimensional form:
U X + V Y = 0 U U X + V U Y = P X + μ n f ρ n f   α f [ 2 U X 2 + 2 U Y 2 ] U V X + V V Y = P Y + μ n f ρ n f   α f [ 2 V X 2 + 2 V Y 2 ] ρ f   σ n f   ρ n f   σ f H a 2  Pr  V + ( ρ β ) n f ρ n f   β f R a Pr θ   U θ X + V θ Y = α n f α f [ 2 θ X 2 + 2 θ Y 2 ]
In dimensionless forms, local entropy generation can be expressed as following:
S c * = [ ( θ X ) 2 + ( θ Y ) 2 ] + μ n f   k f μ f   k n f B r U 2 [ 2 × { ( U X ) 2 + ( V Y ) 2 } + ( U Y + V X ) 2 ]
S c * = S c S G , c ,  and  S G , c = k ( Δ T ) 2 L 2 T 0 2 , B r = μ u 2 k Δ T = E c Pr , Ω = Δ T T 0
where Ω , S c * , and Br are dimensionless temperature difference ratio, dimensionless local entropy generation and Brinkman number, respectively. Brinkman number determines the relative importance between dissipation effects and fluid conduction effects [4]. In the present study, the ratio B r / Ω is fixed at 1.0.
As the distribution of volumetric entropy generation is obtained, it would be integrated over the whole domain to yield the global entropy generation (GEG) rate:
S m = S c * d
The non-dimensional boundary conditions along the cavity walls are described as follows: U, V = 0 on all the walls; θ = 1 on the heat source; θ = 0 on the cool walls; and ∂θ/∂X = 0 for adiabatic walls.
The thermo-physical properties of pure water (base fluid) and the CuO nanoparticles, at a constant temperature of 25 °C, are listed in Table 1 [41,42].
Table 1. Thermo physical properties of water and CuO nanoparticles [41,42].
Table 1. Thermo physical properties of water and CuO nanoparticles [41,42].
Physical propertiesFluid phaseCuO
Cp (j/kg·K)4179540
ρ (kg/m3)997.16500
Pr6.2
k ( w / m · K ) 0.61318
β ( k 1 ) 21 × 10 5 5.1 × 10 5
σ ( μ S / c m )0.05 2.7 × 10 8
The effective density of the nanofluid is defined as follows:
  ρ n f = ( 1 ϕ ) ρ f + ϕ ρ s
where ϕ is the solid volume fraction. The thermal diffusivity of the nanofluid is defined as follows:
α n f = k n f ( ρ c p ) n f
The heat capacitance and the thermal expansion coefficient of the nanofluid are given as follows [18]:
( ρ c p ) n f = ( 1 ϕ ) ( ρ c p ) f + ϕ ( ρ c p ) s   ,   ( ρ β ) n f = ( 1 ϕ )   ( ρ β ) f   +   ϕ ( ρ β ) s
The diameter of nanoparticles (ds) is 33 nm. The effective viscosity of CuO-water nanofluid, its Brownian-motion velocity and the effect of water temperature on this motion, can be expressed as follows [43]:
μ nf μ f = 1 + ρ s V b r d s 2 72   N   δ , V b r = 1 d s ( 18   k b   T π   ρ s   d s ) 1 / 2
where ρ s is the nanoparticles density. In addition, δ denotes the nanoparticles average central distance and N is a parameter for empirical data fitting. They are written as follows:
δ = ( π 6 ϕ ) 1 / 3 × d s   ,   N = ( c 1 ϕ + c 2 ) d s + ( c 3 ϕ + c 4 )
where c1 = −1.133 × 10−6, c2 = −2.771 × 10−6, c3 = 9 × 10−8 and c4 = −3.93 × 10−7.
Chon et al. [44] proposed the following empirical relation for the thermal conductivity of Al2O3-water nanofluid with spherical particles:
k nf k f = 1 + 64.7 ϕ 0.74 ( d f d s ) 0.369 × ( k s k f ) 0.747 Pr 0.9955 Re s 1.2321
In Equation (13),
Re s = ρ f   k b   T 3 π   [ A 1 × 10 A 2 / ( T A 3 ) ] 2   l f , A 1 = 2.414 × 10 5 , A 2 = 247.8 , A 3 = 140
Res is defined as the nanoparticles Reynolds number. In addition, l f = 17   nm is mean free path of water. The accuracy of the above-mentioned equation was also confirmed by Mintsa et al. [45] for CuO and adopted afterwards in many studies, such as by Popa et al. [46].
The effective electrical conductivity of nanofluid was presented by Maxwell [47] as follows:
σ n f σ f = 1 + 3 ( γ 1 ) ϕ ( γ + 2 ) ( γ 1 ) ϕ   ,   γ = σ s σ f
The local Nusselt number over the heat transfer walls is calculated by:
N u = k n f k f ( θ n ) n  with  n = ( X , Y )
The mean Nusselt number over the heat transfer walls is evaluated as:
N u m = 1 3 [ A B N u   d X + B C N u   d Y + C D N u   d X ]

3. Numerical Solution and Grid Dependency Test

Finite volume method using a collocated grid system have been used to solve the governing Equation (4). The computational domain has been uniformly discretized using displaced network method. Central difference scheme was used to discretize diffusion terms, whereas the upwind scheme was adopted for convection terms. The SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm [48] was followed for pressure velocity linkage. All numerical efforts were achieved by an in house-built FORTRAN computer code (FORTRAN 90). Given that the basis of the solving algorithm lays on an iterative technique, the following convergence criterion is used:
j i [ δ n + 1 δ n δ n + 1 | i , j ] 2 10 8
where   δ may be any variable of (U, V, θ ). Once the variables U, V, and θ are obtained, the local entropy generation is calculated using Equation (5).
It is necessary to examine the dependency of the simulation results to the grid structure. For this purpose, the effect of the grid points number on the average Nusselt number has been evaluated. The results for Ha = 30, ϕ = 0.04 , AR = 0.1, 0.7 and Ra = 1000, 15,000 are listed in Table 2. As shown, it is clear that for any grid structure that is finer than 100 × 100, the average Nusselt number remains unchanged. Thus, a uniform grid with the size of 100 × 100 for this study has been selected (Figure 2). In order to verify the written code with others, a comparison has been made with selected results reported in [35]. The average Nusselt number was investigated in a C-shaped cavity with three internal cold walls and three external hot walls, the results of which are shown in Figure 3. As shown, difference between the present study and results reported in [35] is acceptable. Extra comparison with Pirmohammadi and Ghasemi [49] was achieved at zero inclination angle and Ra = 105 for different Hartman number. The comparison results are tabulated in Table 3, and they further guarantee the validity of the present numerical code.
Table 2. Grid independency study ( ϕ = 0.04 and Ha = 30).
Table 2. Grid independency study ( ϕ = 0.04 and Ha = 30).
Grid points 40 × 40 60 × 60 80 × 80 100 × 100 120 × 120
AR = 0.1Ra = 10000.68840.67710.63330.63280.6328
Ra = 15,0000.68510.66740.66170.66050.6604
AR = 0.7Ra = 10005.22155.20695.20185.20085.2007
Ra = 15,0005.24335.21635.20525.20265.2024
Figure 2. Schematic presentation of the grid.
Figure 2. Schematic presentation of the grid.
Entropy 18 00050 g002
Figure 3. Validation of the present code against Mahmoodi and Hashemi [35] for a Cu-water nanofluid natural convection in cavity.
Figure 3. Validation of the present code against Mahmoodi and Hashemi [35] for a Cu-water nanofluid natural convection in cavity.
Entropy 18 00050 g003
Table 3. Comparison of mean Nusselt number with [49] for Ra = 105.
Table 3. Comparison of mean Nusselt number with [49] for Ra = 105.
Hartman NumberPresentPirmoammadi and Ghasemi [49]Percentage error
501.01041.0188−0.82%
701.81211.8372−1.36%

4. Results and Discussion

The results are gathered by inspecting the effects of Rayliegh number 1000 ≤ Ra ≤ 15,000, Hartman number 0 ≤ Ha ≤ 45, nanoparticles volume fraction 0 ≤ φ ≤ 0.06 and cavity aspect ratio 0.1 ≤ AR ≤ 0.7 on the flow, thermal, and entropy generation fields and are presented in the following subsections. The effect of nanoparticles volume fraction is included in each subsection.

4.1. Effect of Raylirgh Number

Contour maps of streamlines, isotherms, and isentropic lines for different Ra values at Ha = 30, φ = 0.04, and AR = 0.3 are presented in Figure 4. The streamlines manifests a single vortex behavior. Due to the imposed thermal boundary conditions, the vortex rotates clockwise with its core localized within the main cavity. Increasing Ra from 1000 to 15,000 leads to strengthening the vortex rotation. The nanofluid fluid is intensified close to the vertical cold wall of the rib where relatively denes streamlines are noticed. Undisturbed isotherms with C-shaped behavior is noticed up to Ra = 10,000. Beyond this value, a small perturbation takes place on the isotherms and crowded isotherms close to the rib upper edge are seen. However, evenly disturbed isotherms within the lower and upper grooves are still present with all Ra values. The isentropic contours show that the entropy generation is concentrated along the cold vertical wall of the rib and along a substantial part of the hot vertical wall. The remainder of the cavity behaves as an “idle” region from the entropy generation point of view. The streamlines and the isotherms above lead to attribute the entropy concentrators mainly to the nanofluid friction irreversibility.
Figure 4. The effects of Rayleigh on the Streamlines, Isotherms and Isentropic lines for nanofluid ϕ = 0.04 ) at Ha = 30 and AR = 0.3.
Figure 4. The effects of Rayleigh on the Streamlines, Isotherms and Isentropic lines for nanofluid ϕ = 0.04 ) at Ha = 30 and AR = 0.3.
Entropy 18 00050 g004
The parameters of Figure 4 above are adopted in depicting the distribution of local Nusselt number Nu along the hot walls. The results are presented in Figure 5, which shows that most of the convective heat transfer takes place along the horizontals hot walls. The local Nusselt number along the upper hot wall is an increasing function of Ra number, while it is a decreasing function of Ra along the lower hot wall. Relatively little variation of Nu takes place along the vertical hot wall (BC). Globally, the overall heat transfer is, as expected, an increasing function of Ra due to the dominance of inertia force over the viscous one. This can be drawn from the variation of the average Nusselt number over the hot walls, which is presented in Figure 6. However, a linear increase of the Num with the nanoparticles volume fraction φ is recorded with all studied Ra values. This can be surely attributed to the enhanced thermal energy transport due to the enhanced nanofluid thermal conductivity.
Figure 5. Variation of local Nusselt number for various Rayleigh numbers at Ha = 30, AR = 0.3 and ϕ = 0.04 .
Figure 5. Variation of local Nusselt number for various Rayleigh numbers at Ha = 30, AR = 0.3 and ϕ = 0.04 .
Entropy 18 00050 g005
Figure 6. Variation of average Nusselt number for various Rayleigh numbers with ϕ for AR = 0.3 and Ha = 30.
Figure 6. Variation of average Nusselt number for various Rayleigh numbers with ϕ for AR = 0.3 and Ha = 30.
Entropy 18 00050 g006
The normalized total entropy generation is illustrated in Figure 7 for the same parameters above. An increase of the nanoparticles volume fraction and Rayleigh number lead to an increase in the thermal energy and intensify the nanofluid flow. As a result, a noticeable increase of the total entropy generation with both φ and Ra is recorded in Figure 7. For Ra = 1000, a linear increase of the entropy generation rate is seen. For higher Ra values, a faster increasing function is recorded.
Figure 7. Variation of global entropy generation ratio ( S m , * = S m / S m , φ = 0 ) for various Rayleigh numbers with ϕ for AR = 0.3 and Ha = 30.
Figure 7. Variation of global entropy generation ratio ( S m , * = S m / S m , φ = 0 ) for various Rayleigh numbers with ϕ for AR = 0.3 and Ha = 30.
Entropy 18 00050 g007

4.2. Effect of Hartman Number

The effect of the applied magnetic field is studied by varying Hartman number (Ha = 0–45) at Ra = 10,000, φ = 0.04 and AR = 0.3. The streamlines patterns, which are presented in Figure 8, show a continuous reduction in vortex strength with increasing Ha. When Ha = 0, the vortex rotates strongly under the effect of natural convection. This limits it to within the main cavity. Increasing Ha increases the suppression force, which opposes the convection currents and this enables the weakened vortex to penetrate symmetrically through the cavity grooves. A threshold value of Ha = 30 can be characterized at which the penetration effect starts. Before this threshold value, the isotherms manifest a disturbed pattern, which is an indication of the convection dominance. Beyond the threshold value, evenly distributed isotherms are noticed, which indicates a weakening of the convection. In the absence of the magnetic field (Ha = 0), the entropy generation concentrators are not only localized along the vertical hot and cold walls, but the entropy is also generated from a region localized in the middle cavity. However, the later region diminishes with increasing Ha number. The local Nusselt number shown in Figure 9 indicates a convection effect at Ha = 0, where asymmetric distribution along the hot walls is noticed. The symmetrical behavior is mostly attained with increasing Ha. It is worth mentioning that the local Nusselt number increases with Ha along the hot upper wall (AB) only. This might be attributed to the magnetic field effect in carrying the hot nanofluid close to the wall into the cavity. The average Nusselt number normalized by that of pure fluid (Num,φ = 0) is presented in Figure 10. The noticed increase of the normalized Nusselt number with Ha reflects that the drag action of the magnetic field on the nanofluid is less than that on pure fluid.
Figure 8. The effects of Hartmann numbers on the Streamlines, Isotherm and Entropy generation for pure-water (—) and nanofluid ϕ = 0.04 (- - -) at Ra = 10,000 and AR = 0.3.
Figure 8. The effects of Hartmann numbers on the Streamlines, Isotherm and Entropy generation for pure-water (—) and nanofluid ϕ = 0.04 (- - -) at Ra = 10,000 and AR = 0.3.
Entropy 18 00050 g008
Figure 9. Variation of local Nusselt number for various Hartman numbers at Ra = 10,000, AR = 0.3 and ϕ = 0.04.
Figure 9. Variation of local Nusselt number for various Hartman numbers at Ra = 10,000, AR = 0.3 and ϕ = 0.04.
Entropy 18 00050 g009
Figure 10. Variation of average Nusselt number ratio ( Nu m , * = Nu m / Nu m , φ = 0 ) for various Hartman number with ϕ at Ra = 10,000 and AR = 0.3.
Figure 10. Variation of average Nusselt number ratio ( Nu m , * = Nu m / Nu m , φ = 0 ) for various Hartman number with ϕ at Ra = 10,000 and AR = 0.3.
Entropy 18 00050 g010
The entropy generation rate with Ha is shown in Figure 11. Since the natural convection is already weak at Ra = 1000, there is no pronounced action of the magnetic field. However, there is a steep reduction of the entropy generation rate with Ha within Ra > 1000. This steep reduction is up to the threshold Hartman number value. To evaluate the thermal performance of the present geometry with its physical constraints, we adopted the ε criterion proposed by Ismael et al. [32]. It is the ratio of the total entropy generation to the average Nusselt number. According to this criterion, which is depicted in Figure 12, it seems that applying the magnetic field at low Rayleigh number (Ra = 1000) is a useless process. Meanwhile, for higher Ra values (Ra > 1000), the magnetic field improves the thermal performance. A value of Ha = 30 can be accepted as the best if one considers the power consumed in the magnet generating the magnetic field.
Figure 11. Variation of global entropy generation for various Hartman numbers with Ra for AR = 0.3, and ϕ = 0.04 .
Figure 11. Variation of global entropy generation for various Hartman numbers with Ra for AR = 0.3, and ϕ = 0.04 .
Entropy 18 00050 g011
Figure 12. Variation of ε ( ε = S m / Nu m ) for various Rayleigh numbers and Hartman numbers for AR = 0.3 and ϕ = 0.04 .
Figure 12. Variation of ε ( ε = S m / Nu m ) for various Rayleigh numbers and Hartman numbers for AR = 0.3 and ϕ = 0.04 .
Entropy 18 00050 g012

4.3. Effect of the Aspect Ratio

Figure 13 shows the contour maps for different aspect ratio at Ra = 10,000, φ = 0.0, 0.04 and Ha = 30. For very small aspect ratio (AR = 0.1), the available free space in the square-like cavity enables a clockwise vortex with bit concavity close to the cold rib. Increasing the aspect ratio limits the space available for vertex rotation, which leads to a weakened strength of rotation. When AR is further increased to AR = 0.7, the cavity grooves look like stagnant regions. The lower aspect ratio cavity (AR = 0.1) limits the heat outfall, which in turn limits the thermal boundary layer around the rib. The thermal boundary layer is developed when AR increases. A pure conduction is noticed within the cavity grooves. This is because the little vertical distance between the horizontal hot and cold walls. The isentropic lines show that the entropy generation rate is significantly increased with increasing the aspect ratio. However, the dashed contours of φ = 0.4 indicate to little effect of nanoparticles on the nanofluid flow, heat transfer, and entropy generation field.
Figure 13. The effects of AR on the Streamlines, Isotherms and Isentropic lines for nanofluid at φ = 0, 0.04, Ra = 10,000 and Ha = 30.
Figure 13. The effects of AR on the Streamlines, Isotherms and Isentropic lines for nanofluid at φ = 0, 0.04, Ra = 10,000 and Ha = 30.
Entropy 18 00050 g013aEntropy 18 00050 g013b
The local Nusselt number distributions displayed in Figure 14 shows that for low aspect ratio, the hot walls contribute in heat transfer with lesser discrepancy. Increasing the aspect ratio leads to significant contribution weights among the hot walls in heat transfer. Figure 15 clearly shows that the aspect ratio is significantly increased the average Nusselt number. Moreover, enhancement percentages in the average Nusselt number are 19.5% and 17.6% for ϕ = 0.06 compared to pure fluid ( ϕ = 0.0) when the aspect ratio increased from 0.1 to 0.7. This result can be attributed to the limited space available, which reduces the strength of vortex rotation and then limits the utilization of nanoparticles addition. The total entropy generation increases with increasing the aspect ratio, as shown in Figure 16. Recalling Figure 15, it can be concluded that the heat transfer irreversibility has a great role in generating the entropy. It is worthwhile to show that the entropy generation rate decreases with aspect ratio AR = 0.1−0.5 at the same rate, while at AR = 0.7, the decreasing trend is less steep. The percentage reductions in total entropy generation are 96.27% and 48.17% for Ha = 45 compared to zero magnetic field when the aspect ratio increased from 0.1 to 0.7. Eventually, and according to the ε total entropy/Nusselt number criterion, the applied magnetic field serves in improving the performance (see Figure 17). A threshold value of Ha = 30 may be economically the best value for the applied magnetic field.
Figure 14. Variation of local Nusselt number for various AR at ϕ = 0.04 , Ra = 10,000 and Ha = 30.
Figure 14. Variation of local Nusselt number for various AR at ϕ = 0.04 , Ra = 10,000 and Ha = 30.
Entropy 18 00050 g014
Figure 15. Variation of average Nusselt number for various AR at Ra = 10,000 and Ha = 30.
Figure 15. Variation of average Nusselt number for various AR at Ra = 10,000 and Ha = 30.
Entropy 18 00050 g015
Figure 16. Variation of average Entropy generation for various AR with Hartman numbers for Ra = 10,000 and ϕ = 0.04 .
Figure 16. Variation of average Entropy generation for various AR with Hartman numbers for Ra = 10,000 and ϕ = 0.04 .
Entropy 18 00050 g016
Figure 17. Variation of ε ( ε = S m / Nu m ) for various AR with Hartman numbers for Ra = 10,000 and ϕ = 0.04 .
Figure 17. Variation of ε ( ε = S m / Nu m ) for various AR with Hartman numbers for Ra = 10,000 and ϕ = 0.04 .
Entropy 18 00050 g017

5. Conclusions

Entropy generation and natural convection of CuO-water nanofluid in a C-shaped cavity subjected to a horizontal uniform magnetic field is numerically studied. The effects of nanoparticles volume fraction, Rayleigh number, Hartman number, and the aspect ratio are studied. The effect of Brownian motion is considered in the nanofluid properties. The results have led to the following concluding remarks.
(1)
The addition of nanoparticles enhances the convective heat transfer inside the C-shaped cavity at all Rayleigh numbers, whereas the entropy generation increases with increasing the volume fraction of the nanoparticles. This increase becomes fast at higher Rayleigh number.
(2)
The average Nusselt number increases considerably when the hot and cold walls become narrower, i.e., at higher aspect ratio.
(3)
The nanofluid utilization becomes more pronounced at lower aspect ratio.
(4)
The applied magnetic field is an inactive process at lower Rayleigh number.
(5)
The entropy generation rate decreases rapidly with the applied magnetic field.
(6)
A threshold value of Hartman number equal to 30 can give the best thermal performance in the C-shaped cavity.

Acknowledgments

The authors wish to express their sincere thanks to the competent Reviewers for the valuable comments and suggestions.

Author Contributions

Ali Chamkha has directed to achieve this physical problem and proposed the formulation of this problem and has edited the complete manuscript. Abbas Kasaeipoor and Taher Armaghani conducted the numerical solution and post processed the results. Muneer Ismael wrote the sections and analyzed the physical behaviors of some results of this manuscript. All authors contributed in directing the conclusions and assisted in revising the manuscript according to the Reviewers comments. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

B0Magnetic field strength, TGreek symbols
C p Specific heat, J·kg−1·K−1 α Thermal diffusivity, m2·s−1
g Gravitational acceleration, m·s2 β Thermal expansion coefficient, K−1
H Length of heat source, mεperformance criterion ( S m / N u m )
H a Hartmann number, B 0 L σ f / ρ f ν f ϕ Solid volume fraction
kThermal conductivity, W·m−1·K−1 σ Effective electrical conductivity, μ·S/cm
L Length of cavity, m κ b Boltzmann constant, J·K−1
N u Local Nusselt number θ Dimensionless temperature, ( T T c ) / ( T h T c )
NumAverage Nusselt number of heat source μ Dynamic viscosity, N·S·m−2
p Fluid pressure, Pa ν Kinematic viscosity, m2·s−1
P Dimensionless pressure, p L / ρ n f α f 2 ρ Density, kg·m3
Pr Prandtl number, νf/αfSubscripts
T Temperature, K c Cold
TcCold wall temperature, K f Pure fluid
T h Hot wall temperature, K h hot wall
u , v Velocity components in x, y directions, m·s−1 m Average
U , V Dimensionless Velocity components, ( u , v ) L / α f n f Nanofluid
x , y Cartesian coordinates, m s Nanoparticle
X,YDimensionless coordinates, (x,y)/L

References

  1. Bejan, A. A study of entropy generation in fundamental convective heat transfer. J. Heat Transf. 1979, 101, 718–725. [Google Scholar] [CrossRef]
  2. Bejan, A. Second-law analysis in heat transfer and thermal design. Adv. Heat Transf. 1982, 15, 1–58. [Google Scholar] [CrossRef]
  3. Bejan, A. Entropy Generation Minimization; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  4. Mahmud, S.; Fraser, R.A. Second law analysis of heat transfer and fluid flow inside a cylindrical annular space. Exergy Int. J. 2002, 2, 322–329. [Google Scholar] [CrossRef]
  5. Mahmud, S.; Islam, A.K.M.S. Laminar free convection and entropy generation inside an inclined wavy enclosure. Int. J. Therm. Sci. 2003, 42, 1003–1012. [Google Scholar] [CrossRef]
  6. Bouabid, M.; Magherbi, M.; Hidouri, N.; Ben Brahim, A. Entropy generation at natural convection in an inclined rectangular cavity. Entropy 2011, 13, 1020–1033. [Google Scholar] [CrossRef]
  7. Cheng, X.; Liang, X. Discussion on the applicability of entropy generation minimization to the analyses and optimizations of thermodynamic processes. Energy Convers. Manag. 2013, 73, 121–127. [Google Scholar] [CrossRef]
  8. Baytas, A.C. Entropy generation for natural convection in an inclined porous cavity. Int. J. Heat Mass Transf. 2000, 43, 2089–2099. [Google Scholar] [CrossRef]
  9. Al-Hadhrami, A.K.; Elliott, L.; Ingham, D.B. A New model for viscous dissipation in porous media across a range of permeability values. Transp. Porous Media 2003, 53, 117–122. [Google Scholar] [CrossRef]
  10. Mahmud, S.; Fraser, R.A. Magneto-hydrodynamic free convection and entropy generation in a square porous cavity. Int. J. Heat Mass Transf. 2004, 47, 3245–3256. [Google Scholar] [CrossRef]
  11. Hooman, K.; Ejlali, A. Entropy generation for forced convection in a porous saturated circular tube with uniform wall temperature. Int. Commun. Heat Mass Transf. 2007, 34, 408–419. [Google Scholar] [CrossRef]
  12. Zahmatkesh, I. On the importance of thermal boundary conditions in heat transfer and entropy generation for natural convection inside a porous enclosure. Int. J. Therm. Sci. 2008, 47, 339–346. [Google Scholar] [CrossRef]
  13. Varol, Y.; Oztop, H.F.; Pop, I. Entropy generation due to natural convection in non-uniformly heated porous isosceles triangular enclosures at different positions. Int. J. Heat Mass Transf. 2009, 52, 1193–1205. [Google Scholar] [CrossRef]
  14. Basak, T.; Kaluri, R.S.; Balakrishnan, A.R. Entropy generation during natural convection in a porous cavity: Effect of thermal boundary conditions. Numer. Heat Transf. 2012, 62, 336–364. [Google Scholar] [CrossRef]
  15. Choi, S.U.S.; Eastman, J.A. Enhancing Thermal Conductivity of Fluid with Nanoparticles. In Proceedings of 1995 ASME International Mechanical Engineering Congress and Exhibition, San Francisco, CA, USA, 12–17 November 1995.
  16. Aminossadati, S.M.; Ghasemi, B. The effects of orientation of an inclined enclosure on laminar natural convection. Int. J. Heat Technol. 2005, 23, 43–49. [Google Scholar]
  17. Ghasemi, B.; Aminossadati, S.M. Brownian motion of nanoparticles in a triangular enclosure with natural convection. Int. J. Therm. Sci. 2010, 49, 931–940. [Google Scholar] [CrossRef]
  18. Abu-Nada, E.; Chamkha, A.J. Effect of nanofluid variable properties on natural convection in enclosures filled with a CuO–EG–water nanofluid. Int. J. Therm. Sci. 2010, 49, 2339–2352. [Google Scholar] [CrossRef]
  19. Nemati, H.; Frhadi, M.; Sedighi, K.; Ashorynejad, H.R.; Fattahi, E. Magnetic field effects on natural convection flow of nanofluid in a rectangular cavity using the Lattice Boltzmann model. Sci. Iran. 2012, 19, 303–310. [Google Scholar] [CrossRef]
  20. Mahmoudi, A.H.; Pop, I.; Shahi, M. Effect of magnetic field on natural convection in a triangular enclosure filled with nanofluid. Int. J. Therm. Sci. 2012, 59, 126–140. [Google Scholar] [CrossRef]
  21. Habibi, M.M.; Pop, I. Natural convection flow and heat transfer in an eccentric annulus filled by Copper nanofluid. Int. J. Heat Mass Transf. 2013, 61, 353–364. [Google Scholar] [CrossRef]
  22. Sheikhzadeh, G.A.; Dastmalchi, M.; Khorasanizadeh, H. Effects of nanoparticles transport mechanisms on Al2O3-water nanofluid natural convection in a square enclosure. Int. J. Therm. Sci. 2013, 66, 51–62. [Google Scholar] [CrossRef]
  23. Sheikholeslami, M.; Gorji-Bandpy, M.; Seyyedi, S.M.; Ganji, D.D.; Rokni, H.B.; Soleimani, S. Application of LBM in simulation of natural convection in a nanofluid filled square cavity with curve boundaries. Powder Technol. 2013, 247, 87–94. [Google Scholar] [CrossRef]
  24. Sheikholeslami, M.; Gorji-Bandpy, M.; Ellahi, R.; Zeeshan, A. Simulation of MHD CuO-water nanofluid flow and convective heat transfer considering Lorentz forces. J. Magn. Magn. Mater. 2014, 369, 69–80. [Google Scholar] [CrossRef]
  25. Kashani, S.; Ranjbar, A.A.; Mastiani, M.; Mirzaei, H. Entropy generation and natural convection of nanoparticle-water mixture (nanofluid) near water density inversion in an enclosure with various patterns of vertical wavy walls. Appl. Math. Comput. 2014, 226, 180–193. [Google Scholar] [CrossRef]
  26. Cho, C.-C.; Yau, H.-T.; Chiu, C.-H.; Chiu, K.-C. Numerical investigation into natural convection and entropy generation in a nanofluid-filled U-shaped cavity. Entropy 2015, 17, 5980–5994. [Google Scholar] [CrossRef]
  27. Ting, T.W.; Hung, Y.M.; Guo, N. Entropy generation of viscous dissipative nanofluid flow in thermal non-equilibrium porous media embedded in microchannels. Int. J. Heat Mass Transf. 2015, 81, 862–877. [Google Scholar] [CrossRef]
  28. Mahmoudi, A.; Mejri, I.; Abbassi, M.A.; Omri, A. Analysis of the entropy generation in a nanofluid-filled cavity in the presence of magnetic field and uniform heat generation/absorption. J. Mol. Liq. 2014, 198, 63–77. [Google Scholar] [CrossRef]
  29. Aghaei, A.; Khorasanizadeh, H.; Sheikhzadeh, G.; Abbaszadeh, M. Numerical study of magnetic field on mixed convection and entropy generation of nanofluid in a trapezoidal enclosure. J. Magn. Magn. Mater. 2016, 403, 133–145. [Google Scholar] [CrossRef]
  30. Kefayati, G.H.R. Simulation of heat transfer and entropy generation of MHD natural convection of non-Newtonian nanofluid in an enclosure. Int. J. Heat Mass Transf. 2016, 92, 1066–1089. [Google Scholar] [CrossRef]
  31. Ellahi, R.; Hassan, M.; Zeeshan, A. Shape effects of nanosize particles in Cu-H2O nanofluid on entropy generation. Int. J. Heat Mass Transf. 2015, 81, 449–456. [Google Scholar] [CrossRef]
  32. Ismael, M.A.; Armaghani, T.; Chamkha, A.J. Conjugate heat transfer and entropy generation in a cavity filled with a nanofluid-saturated porous media and heated by a triangular solid. J. Taiwan Inst. Chem. Eng. 2015, in press. [Google Scholar] [CrossRef]
  33. Biserni, C.; Rocha, L.A.O.; Stanescu, G.; Lorenzini, E. Constructal H-shaped cavities according to Bejan’s theory. Int. J. Heat Mass Transf. 2007, 50, 2132–2138. [Google Scholar] [CrossRef]
  34. Mahmoodi, M. Numerical simulation of free convection of a nanofluid in L-shaped cavities. Int. J. Therm. Sci. 2011, 50, 1731–1740. [Google Scholar] [CrossRef]
  35. Mahmoodi, M.; Hashemi, S.M. Numerical study of natural convection of a nanofluid in C-shaped enclosures. Int. J. Therm. Sci. 2012, 55, 76–89. [Google Scholar] [CrossRef]
  36. Cho, C.-C.; Yau, H.-T.; Chen, C.-K. Enhancement of natural convection heat transfer in a U-shaped cavity filled with Al2O3-water nanofluid. Therm. Sci. 2012, 16, 1317–1323. [Google Scholar] [CrossRef]
  37. Mansour, M.A.; Bakeir, M.A.; Chamkha, A. Natural convection inside a C-shaped nanofluid-filled enclosure with localized heat sources. Int. J. Numer. Method Heat Fluid Flow 2014, 24, 1954–1978. [Google Scholar] [CrossRef]
  38. Mojumder, S.; Saha, S.; Saha, S.; Mamun, M.A.H. Effect of magnetic field on natural convection in a C-shaped cavity filled with ferrofluid. Procedia Eng. 2015, 105, 96–104. [Google Scholar] [CrossRef]
  39. Kasaeipoor, A.; Ghasemi, B.; Aminossadati, S.M. Convection of Cu-water nanofluid in a vented T-shaped cavity in the presence of magnetic field. Int. J. Therm. Sci. 2015, 94, 50–60. [Google Scholar] [CrossRef]
  40. Al-Zamily, A.M.J. Effect of magnetic field on natural convection in a nanofluid-filled semi-circular enclosure with heat flux source. Comput. Fluids 2014, 103, 71–85. [Google Scholar] [CrossRef]
  41. Sheikholeslami, M.; Ganji, D.D.; Gorji-Bandpy, M.; Soleimani, S. Magnetic field effect on nanofluid flow and heat transfer using KKL model. J. Taiwan Inst. Chem. Eng. 2014, 45, 795–807. [Google Scholar] [CrossRef]
  42. Öğüt, E.B. Natural convection of water-based nanofluids in an inclined enclosure with a heat source. Int. J. Therm. Sci. 2009, 48, 2063–2073. [Google Scholar] [CrossRef]
  43. Masoumi, N.; Sohrabi, N.; Behzadmehr, A. A new model for calculating the effective viscosity of nanofluids. J. Phys. D 2009, 42, 055501. [Google Scholar] [CrossRef]
  44. Chon, C.H.; Kihm, K.D.; Lee, S.P.; Choi, S.U.S. Empirical correlation finding the role of temperature and particle size for nanofluid (Al2O3) thermal conductivity enhancement. Appl. Phys. Lett. 2005, 87, 153107. [Google Scholar] [CrossRef]
  45. Mintsa, H.A.; Roy, G.; Nguyen, C.T.; Doucet, D. New temperature dependent thermal conductivity data for water-based nanofluids. Int. J. Therm. Sci. 2009, 48, 363–371. [Google Scholar] [CrossRef]
  46. Popa, C.V.; Kasaeian, A.B.; Nasiri, S.; Korichi, A.; Polidori, G. Natural convection heat and mass transfer modeling for Cu/water and CuO/water nanofluids. Adv. Mech. Eng. 2013. [Google Scholar] [CrossRef]
  47. Maxwell, J.C. A Treatise on Electricity and Magnetism, 2nd ed.; Oxford University Press: Cambridge, UK, 1904. [Google Scholar]
  48. Patankar, S. Numerical Heat Transfer and Fluid Flow; CRC Press: Boca Raton, FL, USA, 1980. [Google Scholar]
  49. Pirmohammadi, M.; Ghassemi, M. Effect of magnetic field on convection heat transfer inside a tilted square enclosure. Int. Commun. Heat Mass Transf. 2009, 36, 776–780. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Chamkha, A.; Ismael, M.; Kasaeipoor, A.; Armaghani, T. Entropy Generation and Natural Convection of CuO-Water Nanofluid in C-Shaped Cavity under Magnetic Field. Entropy 2016, 18, 50. https://doi.org/10.3390/e18020050

AMA Style

Chamkha A, Ismael M, Kasaeipoor A, Armaghani T. Entropy Generation and Natural Convection of CuO-Water Nanofluid in C-Shaped Cavity under Magnetic Field. Entropy. 2016; 18(2):50. https://doi.org/10.3390/e18020050

Chicago/Turabian Style

Chamkha, Ali, Muneer Ismael, Abbas Kasaeipoor, and Taher Armaghani. 2016. "Entropy Generation and Natural Convection of CuO-Water Nanofluid in C-Shaped Cavity under Magnetic Field" Entropy 18, no. 2: 50. https://doi.org/10.3390/e18020050

APA Style

Chamkha, A., Ismael, M., Kasaeipoor, A., & Armaghani, T. (2016). Entropy Generation and Natural Convection of CuO-Water Nanofluid in C-Shaped Cavity under Magnetic Field. Entropy, 18(2), 50. https://doi.org/10.3390/e18020050

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