A computational study has been performed to make a computational analysis of natural convection and entropy generation in a sharp edged finned cavity. Three dimensional analysis has been done by solving governing equations with a written computational code in Fortran. The study is performed for fin inclination angles from −60° to 60°, Ra = 10^{5}, Pr = 0.7, R_{c} (conductivities ratio) changes from 0.01 to 100 and irreversibility coefficient is taken as . It is observed that higher values of thermal conductivity ratio (Rc ⩾ 1) do not affect entropy generation due to heat transfer and Bejan number.
Entropy generation; Inclined fin; Natural convection; Three dimensional
Be Bejan number
Cp specific heat at constant pressure (J/kg K)
g gravitational acceleration (m/s^{2})
k thermal conductivity (W/m K)
l cavity width
n unit vector normal to the wall
N_{s} dimensionless local generated entropy
Nu local Nusselt number
Pr Prandtl number
Ra Rayleigh number
Rc thermal conductivity ratio
 generated entropy (kJ/kg K)
t dimensionless time
T dimensionless temperature
 cold temperature (K)
 hot temperature (K)
T_{o} bulk temperature (K)
 dimensionless velocity vector
x, y, z dimensionless Cartesian coordinates
 thermal diffusivity (m^{2}/s)
 thermal expansion coefficient (1/K)
 irreversibility coefficient
 fin inclination
 density (kg/m^{3})
 dynamic viscosity (kg/m s)
 kinematic viscosity (m^{2}/s)
 dimensionless vector potential ()
 dimensionless vorticity
 dimensionless temperature difference
av average
x, y, z Cartesian coordinates
fr friction
f fluid
av average
nf nanofluid
s solid
th thermal
tot total
′ dimensional variable
Natural convection heat transfer occurs in many field of engineering such as cooling of electronic equipment, solar collectors, heating of buildings, heat exchangers and many other applications. Control of energy is extremely important to save energy consumption. Thus, entropy generation analysis is highly important to answer the question of why or where energy is consumed in the system?
Using of passive systems such as baffle, fin or obstacle is a simple way to control heat and fluid flow due to natural convection in cavities. Varol et al. [1] studied the laminar natural convection heat transfer in an inclined fin located cavity by solving two dimensional governing equations. Also, they performed an experiment to compare their results and found the inclination angle of the fin is the important parameter to control heat and fluid flow. As a similar work, Ozgen et al. [2] tested the effects of Prandtl number and baffle location on natural convection. Li and Byon [3] studied the orientation effect on the thermal performance of radial heat sinks with a circular base, concentric ring, and rectangular fins both experimentally and numerically. They observed that the radial heat sink in orientation 0 is found to have slightly better thermal performance than that of 180. Saravanan and Sivaraj [4] performed a numerical work on combined natural convection and thermal radiation in a square cavity with a nonuniformly heated plate. Effects of heat generating baffles on natural convection in a square cavity were studied by Saravanan and Vidhya Kumar [5]. In their study, the baffles were kept centrally separated by a distance d in the cavity and assumed to generate heat uniformly. They observed that flow inhibition mechanism was dominating only when the baffles are away from the cavity walls. Other related works on twodimensional analysis of natural convection in inclined fin attached cavities are studied in Refs. [6], [7], [8] and [1].
Importance of analyses of entropy generation from the energy efficiency point of view is presented in reviewable work of Öztop and AlSalem [9] and Sciacovelli et al. [10]. They mainly focused on entropy generation in natural and mixed convection heat transfer. Shuja et al. [11] solved a problem on natural convection in a square cavity with a heat generating body by considering entropy generation. Khorasanizadeh et al. [12] performed a numerical study of natural convection and entropy generation of Cuwater nanofluid within an enclosure with a conductive baffle embedded on bottom hot wall. They solved the governing equations numerically with Finite Volume Method by using the SIMPLER algorithm. They showed that assuming a constant irreversibility factor, χ, with change of Ra and ϕ is not correct.
Entropy generation due to natural convection in an enclosure heated locally from below with two isoflux sources was investigated by Mukhopadhyay [13]. They studied unequal heater length and heater strength, the effects of heater length and strength ratios. It is found that minimum entropy generation rate was achieved for the same condition at which the minimum peak heater temperature was obtained. Famouri and Hooman [14] studied the entropy generation in twodimensional domain to investigate the entropy generation for natural convection by heated partitions in a cavity.
Dalir [24] studied numerically the entropy generated by laminar forced convection flow and heat transfer of an incompressible Jeffrey nonNewtonian fluid over a linearly stretching, impermeable and isothermal sheet. Results indicate that the generated entropy increases with the increase in Deborah number and decreases by the increase in ratio of relaxation to retardation times. Das et al. [25] investigated the entropy generation in a MHD pseudoplastic nanofluid flow through a porous channel with convective heating. Three different types of nanoparticles were considered and several parameters were varied to find the minimum of entropy generation.
Studies on three dimensional domain of natural convection are very limited. In this context, Wakashima and Saitoh [15] and Fusegi et al. [16] made studies to analyze the natural convection in three dimensional domains.
The main aim of this work was to study the effects of presence of inclined fin on natural convection heat transfer fluid flow and heat transfer in a three dimensional domain of cavity. Entropy generation is calculated by using obtained data of velocity and temperature. Entropy generation due to heat transfer, fluid flow and total entropy generation are presented via contours and variation according to thermal conductivity ratio.
The considered physical model is presented in Fig. 1(a) (on the left) and Fig. 1(b) (on the right). As seen from the figure, the cavity is heated from the bottom wall and cooled from the top while vertical walls are adiabatic. A sharp edged fin is located to left vertical wall to control heat and fluid flow. This fin is conductive and its conductivity changed from R_{c} = 0.01 to 100. Inclination angle of the fin changes from −60° to 60° and thickness and length of the fin are taken as constant.

Figure 1. (a) Three dimensional physical model, and (b) twodimensional view.

The 3D vorticityvector potential formalism is used as numerical method, in order to eliminate the pressure term, which is delicate to treat. and are respectively defined by the two following relations:

(1) 
After nondimensionalization the system of equations controlling the phenomenon is as follows:

(2) 

(3) 

(4) 

The energy equation (conduction) needs to be solved in the solid portion of the domain, and the fin conductivity k_{s} is assumed constant. At the solid–fluid interface the temperature and heat flux must be continuous. The latter requirement is mathematically expressed as

(5) 
where is the thermal conductivity ratio between the material of the fin and the medium that fills the cavity. There is no flow in the fraction of the volume occupied by the blinds (solid). The control volume finite difference method is used to discretize Eqs. (1), (2), (3) and (4). The centraldifference scheme for treating convective terms and the fully implicit procedure to discretize the temporal derivatives are retained. The grid is uniform in all directions with additional nodes on boundaries. The successive relaxation iterating scheme is used to solve the resulting nonlinear algebraic equations. The time step 10^{−4} and spatial mesh 81 × 81 × 81 are retained to carry out all numerical tests. The solution is considered acceptable when the following convergence criterion is satisfied for each step of time as

(6) 
Boundary conditions for considered model are given as follows:
Temperature boundary conditions

on all other walls (adiabatic).
Vorticity

Vector potential

Velocity

Theory on entropy generation is given in the literature [17], [18], [19], [20], [21], [22] and [23]. The generated entropy is written in the following form:

The first term represents the generated entropy due to temperature gradient and the second due to the friction effects.

The dissipation function is written in incompressible flow as follows:

(5) 
from where the generated entropy is written as

(6) 
After nondimensionalization, the generated entropy number (dimensionless local entropy generated) is written in the following way as

(7) 
from where

(8) 
with is the irreversibility coefficient.
The first term of N_{s} represents the local irreversibility due to the temperatures gradients, and it is noted N_{S}_{}_{th}. The second term represents the contribution of the viscous effects in the irreversibility and it is noted N_{S}_{}_{fric}. N_{s} give a good idea on the profile and the distribution of the generated local dimensionless entropy. The total dimensionless generated entropy is written as follows:

(9) 
Bejan number (Be) is the ratio of heat transfer irreversibility to the total irreversibility due to heat transfer and fluid friction:

(10) 
The local and average Nusselt numbers at the cold wall are given by the following:

(11) 
Validation of the code is performed against data in published studies. Thus, results are compared with studies by Wakashima and Saitoh [15] and Fusegi et al. [16] for differentially heated cubic cavity and Pr = 0.71 and listed in Table 1. The table shows that the obtained results are acceptable when compared with results in the literature. Also, obtained results are compared by 2D solutions and compared with the literature as seen in Fig. 2. As seen from the figures there is a good agreement between solutions.
Ra  Authors  (center)  (center)  (y)  (x)  

10^{4}  Present work  0.05528  1.1063  0.199(0.826)  0.221(0.112)  2.062 
Wakashima and Saitoh [15]  0.05492  1.1018  0.198(0.825)  0.222(0.117)  2.062  
Fusegi et al. [16]  –  –  0.201(0.817)  0.225(0.117)  2.1  
10^{5}  Present work  0.034  0.262  0.143(0.847)  0.245(0.064)  4.378 
Wakashima and Saitoh [15]  0.03403  0.2573  0.147(0.85)  0.246(0.068)  4.366  
Fusegi et al. [16]  –  –  0.147(0.855)  0.247(0.065)  4.361  
10^{6}  Present work  0.01972  0.1284  0.0832(0.847)  0.254(0.032)  8.618 
Wakashima and Saitoh [15]  0.01976  0.1366  0.0811(0.86)  0.2583(0.032)  8.6097  
Fusegi et al. [16]  –  –  0.0841(0.856)  0.259(0.033)  8.77 

Figure 2. Comparison of results of present code with those of Varol and Ozgen (2D) [7].

A computational study is performed in this work is to make three dimensional analysis of natural convection and entropy generation. A computer code was written by using Fortran platform. All results are presented for Ra = 10^{5}, Pr = 0.07, and thermal conductivity ratio changes from 0.01 to 100.
Particle trajectories are presented in Fig. 3 for different inclination angles of the fin and Rc = 1.0. This figure gives an idea on particle transport way inside the cavity. As seen from the figures, particle trajectories are continued for θ = 0°. These figures show that twodimensional solutions are not enough for solution of cubic cavity. Complexity of trajectories is dependent on inclination angle of the fin and it becomes messy behind the fin due to narrow space. Isosurfaces of temperature for Rc = 1.0 at different inclination angles of the fin are given in Fig. 4. This figure gives an idea on three dimensional distribution of temperature. The distribution of temperature is almost parallel to horizontal walls near the walls. However, they exhibit a parallel distribution to each other near the top side of the fin for θ = 60° and 30°. In case of θ = −30° and θ = −60°, temperature distribution behaves almost as empty cavity due to values of thermal conductivity value. As shown in Fig. 5 the fin is the main factor to manager of the flow and heat transfer inside the cavity. This is clearer from the velocity projections in Fig. 5 at different inclination angles and different thermal conductivity ratios. In case of θ = 60°, the flow is captured behind the fin and dimension of its main eye increases with increase in thermal conductivity due to increase in flow strength. The main flow does not change with thermal conductivity ratio. Similar distribution is observed for θ = 30° and velocity distribution is completely symmetric for θ = 0°. The fin separates the flow two ways in case of θ = −30°. A number of cells are affected from the thermal conductivity and almost a triangular crosssectional cavity is formed behind the fin for θ = −60° and the flow behaves as almost empty cavity at the remaining part and the main flow is not affected by changing the thermal conductivity ratio.

Figure 3. Particle trajectories for R_{c} = 1 at different inclination angles of the fin.


Figure 4. Isosurfaces of temperature for R_{c} = 1 at different inclination angles of the fin.


Figure 5. Velocity projection in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.

Fig. 6 illustrates the isotherms for different inclination angles of the fin and thermal conductivity ratio. In case of Rc = 1.0, the fin makes minor effects on temperature distribution inside the cavity. For θ = 60°, isotherms are cumulated on fin at Rc = 0.01 and the fluid behind the fin becomes stagnant due to low conductivity of the fin. This situation is valid for the case of θ = 30°. For the case of high conductivity, inclination angle of the fin becomes an important parameter. For θ = 0°, symmetric distribution is observed for temperature for all values of thermal conductivity due to bottom heating of the cavity.

Figure 6. Isotherms in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.

Fig. 7 is plotted to show the local entropy generation due to heat transfer for different inclination angles of the fin and different thermal conductivity ratios. As seen from the figure, inclination angle of the fin plays an important role on distribution of local entropy generation. In case of θ = −60° and Rc = 0.01, local entropy generation contours are accumulated inside the fin due to lower thermal conductivity ratio. For Rc = 1.0, the location of the fin becomes insignificant on entropy generation. However, edge of the fin plays an important role on entropy generation for Rc = 100. Almost symmetric distribution is observed for θ = 0°, and plume like distribution is observed for all values of thermal conductivity ratios due to separation effect of the fin to cavity. For the positive values of inclination angles, the flow inside the cavity is deferred to right top corner and energy is lost at that part mostly.

Figure 7. Local entropy generation due to heat transfer in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.

Local entropy generation due to friction is illustrated in Fig. 8. Effects of inclination angle of the fin and thermal conductivity ratio are evaluated in this figure. It can be seen that thermal conductivity ratio is not effective on entropy distribution and contours are cluster middle of the right vertical wall and bottom of the inclined fin for θ = 60°. Bottom edge of the fin affects entropy generation due to moving flow from the bottom and flow separation at that point for θ = 30°. Similar situation is observed for θ = 0° but top and bottom edges of the fin become effective. When the fin turns the negative axis only top edge becomes effective on entropy generation. For θ = −60°, thermal conductivity ratio plays an important role and cumulated entropy contour moves from bottom to edge of the fin while thermal conductivity increases.

Figure 8. Local entropy generation due to friction in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.

Local total entropy generation is shown in Fig. 9 in z = 0.5 plan for different inclination angles of the fin and different thermal conductivity ratios. The figure shows that edge of the fin becomes effective on entropy generation for all parameters due to flow separation at that point and impinged flow onto the fin. Also, entropy generation changes place to place inside the cavity based on flow distribution and position of the fin. As obtained from the flow and temperature distribution, local entropy distribution shows symmetrical distribution for θ = 0°. The inclination angle of the fin can be a control parameter for energy saving. For example, entropy generation values become stronger near the edge and top side of the cavity for θ = 60 and Rc = 100. Edge effect of the fin becomes minimum for Rc = 1.0 for all values of inclination.

Figure 9. Local total entropy generation in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.

Fig. 10 presents the distribution of local Bejan number which is calculated via Eq. (10). This distribution is given in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios. As indicated earlier that the presence of the fin is a power source for entropy generation, Fig. 11 illustrates the variation of average Nusselt number with the considered parameters of fin angle and thermal conductivity ratio. The average heat transfer shows almost same results for −30° and 30° due to the occurrence of the same effects. In case of Rc = 1.0, a peak is observed. However, higher heat transfer is formed for other inclination angle of the fin and maximum heat transfer is formed for −60 degree at Rc = 100.

Figure 10. Local Bejan number in z = 0.5 plan at different inclination angles of the fin and different thermal conductivity ratios.


Figure 11. Variation of average Nusselt number versus R_{c} at different inclination angles.

Fig. 12 displays the variation of entropy generation due to heat transfer for different parameters. S_{th} value is decreased until Rc = 1.0 and it goes constant due to higher value of thermal conductivity ratio. The minimum value is formed for −30° due to motionless flow under the fin. Entropy generation due to friction for different values of thermal conductivity and inclination angle is presented in Fig. 13. The figures show that position of the fin becomes very effective on entropy generation due to friction. Values are decreased with increase in thermal conductivity almost linearly except −60° of the fin angle.

Figure 12. Entropy generation due to heat transfer versus R_{c} at different inclination angles.


Figure 13. Entropy generation due to friction versus R_{c} at different inclination angles.

Variation of total entropy generation with different inclination angles of the fin and thermal conductivity ratio is presented in Fig. 14. The figure gives that these two parameters are good control parameter or optimization parameter for energy consumption. In other words, maximization or minimization is possible inside the cavity by using fin. As seen from the figure, total entropy generation decreases up to Rc = 1.0 and it becomes constant after these value. It is noted that the fin is more conductive for higher value of Rc. The highest entropy generation occurs for −30 and the lowest value is observed for 0°. For all values there is a minimum point around Rc = 1.0.

Figure 14. Total entropy generation versus R_{c} at different inclination angles.

Fig. 15 illustrates the variation of Bejan number versus Rc for different inclination angles of the conductive fin. As well known from the literature that Bejan number gives the ratio of heat transfer and fluid friction, thus, as seen from the figure, Bejan number decreases with thermal conductivity ratio till Rc = 1.0 and it becomes almost constant for all values of thermal conductivity ratio.

Figure 15. Variation of Bejan number versus R_{c} at different inclination angle.

A three dimensional computational study has been performed in this work in an inclined fin located cavity to obtain heat transfer, temperature distribution, fluid flow and entropy generation. In this work, parameters are taken as Ra = 10^{5}, Pr = 0.7, R_{c} (conductivities ratio) changes from 0.01 to 100 and irreversibility coefficient is taken as . Firstly, three dimensional solution in a cubic cavity brings some advantages according to two dimensional solution to see the effects of heat and fluid flow in different directions. The main effective parameter on heat and fluid flow is the inclination angle of the fin. There is symmetrical distribution on temperature distribution and flow field in a cavity when θ = 0°. Edge of the fin is extremely effective on entropy generation and Bejan number for all studied parameters. Entropy generation distribution is also obtained as symmetric for all values of thermal conductivity ratio. But higher values of thermal conductivity ratio (Rc ⩾ 1) are not effective on entropy generation. For Rc = 1, there is a maximum or minimum value on heat transfer according to inclination angle of the fin. The highest heat transfer is seen for the −60° of the fin angle and Rc = 100. It is an interesting result that heat transfer gives same results for −30° and 30° of the fin angle. Both entropy generation and heat transfer become higher for positive values of inclination angle. Entropy generation due to heat transfer becomes dominant onto friction one.
Published on 12/04/17
Licence: Other
Are you one of the authors of this document?