ISRN Chemical Engineering

Volume 2012 (2012), Article ID 941587, 9 pages

doi:10.5402/2012/941587

Research Article

Parametric Study of Entropy Generation in a Fluid with Internal Heat Generation between Two Rotating Cylinders Subjected to Convective Cooling at the Surface

M. El Haj Assad¹ and Hakan F. Oztop²

¹Department of Energy Technology, School of Science and Technology, Aalto University, P.O. Box 11400, 00076 Aalto, Finland

²Department of Mechanical Engineering, Faculty of Technology, Firat University, 23119 Elazig, Turkey

Received 18 October 2012; Accepted 5 November 2012

Academic Editors: C.-T. Hsieh and P. L. Maffettone

Copyright © 2012 M. El Haj Assad and Hakan F. Oztop. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Entropy generation analysis of a steady fluid flow with internal heat generation between two rotating cylinders was presented. The surface of the inner cylinder was kept at constant temperature while the surface of the outer cylinder was exposed to convection cooling. Analytical expressions for velocity and temperature distributions within the fluid were obtained. The effects of velocity ratio, Biot number, Brinkman number and other dimensionless parameters on the temperature distribution and local and total entropy generation rates were investigated and the results were presented graphically.

1. Introduction

Heat transfer analysis in annuli with rotating inner and outer cylinders is an important research field and has a wide range of many engineering applications such as swirl nozzles, commercial viscometers, journal bearings, chemical mixing devices, and electric motors.

Entropy generation and its minimization have been considered as an effective tool to improve the performance of any heat transfer process. Entropy generation minimization of diabatic distillation column with trays has been investigated using a new approach [1] in which the exchanged heat has been considered as a control variable instead of temperature. A numerical study of mixed convection and entropy generation in Poiseulle-Benard channel in different angles has been presented [2]. A numerical study of entropy generation using 3D in the case of liquid metal laminar natural convection in a differentially heated cubic cavity and in the presence of an external magnetic field orthogonal to the isothermal walls has been carried out [3]. A numerical analysis has been presented to study laminar-forced convection and entropy generation in a helical coil with constant wall heat flux [4] where optimal Reynolds numbers were found to be related to the wall heat flux. Entropy generation analysis has been applied to a natural convection cooling of heat source mounted inside a cavity [5]. Second law analysis has been used to obtain analytically the entropy generation for a laminar flow in semicylindrical ducts subjected to constant wall heat flux [6]. Entropy generation due to convection and conduction mode of heat transfer (conjugate natural convection) in enclosures bounded by vertical solid walls with different thicknesses has been studied [7]. Entropy generation has been obtained analytically by using second law analysis for a laminar flow in rectangular ducts with semicircular ends [8]. In the analysis [8], the effects of constant wall heat flux and constant wall temperature as boundary conditions have been studied.

Entropy generation analysis in an annuli flow with inner cylinder at rest, and rotating outer cylinder has been presented [9] where the entropy generation in the system due to the heat transfer has been considered and the entropy generation due to fluid friction has been neglected. An experimental study of heating fluid between two concentrating cylinders with cavities with rotating inner cylinder and outer cylinder at rest has been conducted [10]. Entropy generation for nonlinear viscoelastic fluid in concentric rotating cylinders subjected to different boundary conditions has been investigated [11] in which the results show that entropy generation number increases with increasing Brinkman number. The problem of natural convection of gas in an annulus between two horizontal coaxial cylinders at rest with uniform internal heat generation and different temperatures at the boundaries has been introduced and investigated numerically [12]. The effect of variable viscosity of fluid flow on entropy generation between two concentric cylinders at rest with a convective cooling at the outer cylinder has been investigated [13]. Entropy generation analysis inside concentric annuli with relative rotation has been presented [14] without investigating the effect of internal heat generation within the fluid and without considering the exact fluid temperature in calculating the entropy generation rate. More research studies related to entropy generation can be found in the literature [15–18].

As can be seen from the above research works related to annular flow between two rotating cylinders, the effect of internal heat generation on local and total entropy generation rates has not yet been investigated in the literature and this study is considered as a first approach to consider that effect.

2. Mathematical Formulation

In this section, We developed a mathematical model for the flow and thermal characteristics of a fluid between two rotating cylinders which rotate in the same direction. The schematic diagram of the physical problem is illustrated in Figure 1. It is assumed that the fluid thermal properties are temperature independent. The flow considered in the present is one-dimensional, steady, laminar and purely tangential. The radius of the inner cylinder is and of the outer cylinder is . It is assumed that the outer cylinder wall thickness is negligible and has the same thermal conductivity as that of the fluid. The inner and the outer cylinders are rotating at angular velocities and , respectively. The inner fluid surface is kept at temperature while the outer surface is subjected to convection by a coolant at temperature .

Figure 1: Schematic diagram of concentric rotating annuli.

2.1. Velocity Distribution

Under the above flow conditions, the momentum equation can be written as

where is the tangential velocity. The boundary conditions are

Consider the following dimensionless parameters:

The momentum equation of (1) can be written in a dimensionless form as

Consider the following dimensionless boundary conditions:

The solution of (5) subjected to the boundary conditions in (6) is obtained as

2.2. Temperature Distribution

The energy equation of the fluid with internal heat generation between the two rotating cylinders can be expressed as

where is the internal heat generation, is the thermal conductivity, is the dynamic viscosity, and is the viscous dissipation. The viscous dissipation is

The thermal boundary conditions are

where is the convection heat transfer coefficient. In order to express (8) in a dimensionless form, the following dimensionless parameters are introduced as

where is the Biot number, is the Brinkman number, is the Eckert number, is the Prandtl number, is the dimensionless internal heat generation, and is the specific thermal capacity. Using the parameters of (11), (8) can be rewritten as

with the dimensionless thermal boundary conditions as

The solution of (12) can be obtained as

where

3. Entropy Generation Analysis

The local entropy generation rate due to fluid temperature gradient and fluid friction (viscous dissipation) is given as [19]

The first term of (18) corresponds to the heat transfer irreversibility and the second term corresponds to the irreversibility due to viscous dissipation. Equation (18) considers the exact local temperature of the fluid obtained by (14). Using the dimensionless parameters defined in (4) and (11), the entropy generation number can be obtained from (2) as

where is the temperature ratio of the coolant fluid to the inner cylinder surface temperature. The total entropy generation is expressed as

where is the cross-sectional area of the gab between the cylinders and is the length of the rotating cylinders. Substituting (19) into (20), the total entropy generation rate can be obtained as

Equation (21) can be written in a dimensionless form as

where is the dimensionless total entropy generation rate.

4. Results and Discussion

Entropy generation of a steady fluid flow with internal heat generation between two rotating cylinders is presented for different parameters. The tangential velocity profiles along the radial direction for the radius ratio () and are presented in Figures 2(a) and 2(b), respectively, at different velocity ratio () values. The figures show that the velocity decreases as increases for and but the velocity has a minimum value for . This minimum velocity can be obtained from the velocity expression given in (7) at . It should be remembered that maximum value of is equal to and the value of should always be larger than 1 to insure that such a minimum in velocity exists.

(a)

(b)

Figure 2: Velocity distribution ((a) ; (b) ).

The dimensionless temperature distributions of the rotating fluid as function of are shown in Figures 3–5 for different set of parameters such as Brinkman number (), Biot number (), and dimensionless internal heat generation () at different velocity ratio values. Figure 3 shows the effect of Brinkman number on the temperature variation for , , and . The figure shows that as Brinkman number increases the temperature increases. The difference between the temperature profiles widens as Brinkman increases, and this difference is very clear in Figure 3(c) for . In all Figures 3(a)–3(c), the temperature has a maximum value at a certain location. The temperature of the outer rotating fluid surface increases as decreases for the same and values but this temperature increases as Brinkman increases for the same and values.

(a)

(b)

(c)

Figure 3: Temperature distribution ((a) ; (b) ; (c) ).

(a)

(b)

Figure 4: Temperature distribution ((a) ; (b) ).

(a)

(b)

Figure 5: Temperature distribution ((a) ; (b) ).

The effect of Biot number on the fluid temperature distribution as function of is given in Figure 4 for , , and . It is shown in Figure 4 that an increase in Biot number results in a decrease in fluid temperature at the outer surface due to better convection cooling. In Figure 4, the fluid temperature has a maximum value for all and values considered in this study. This maximum temperature is closer to the outer surface at lower Biot number due to the lower heat convection to the surroundings. When Biot number increases in more heat will be transferred to the surroundings resulting in a higher fluid temperature gradient at the outer fluid surface as can be seen in Figures 4(a) and 4(b). The maximum temperature is shifted closer to the inner surface as Biot number increases.

The maximum temperature is an important factor because the contribution of thermal effect on entropy generation will be equal to zero at that maximum value. Figures 5(a) and 5(b) show the temperature distribution for the case with no internal heat generation () and with internal heat generation (), respectively, for , , and . It is shown in Figure 5(b) that the internal heat generation results in higher fluid temperature. The maximum temperature value is reached at higher value (closer to the outer surface) as increases.

The variations of local entropy generation number with are shown in Figures 6–8 for the same set of parameters considered in Figures 3–5. Figures 6(a) and 6(c) correspond to Figures 3(a) and 3(c), Figures 7(a) and 7(b) correspond to Figures 4(a) and 4(b), and Figures 8(a) and 8(b) correspond to Figures 5(a) and 5(b). Figure 9 shows the effect of the temperature ratio on the local entropy generation number for , , , and . Figures 6–9 show that the local entropy generation number decreases to reach a minimum value at a certain location and then increases very slowly as the outer surface is approached. It can be seen from (19) that even if we have maximum temperature and minimum velocity at the same location, this will not result in zero entropy generation, that is, .

(a)

(b)

(c)

Figure 6: Local entropy generation number distribution ((a) ; (b) ; (c) ).

(a)

(b)

Figure 7: Local entropy generation number distribution ((a) ; (b) ).

(a)

(b)

Figure 8: Local entropy generation number distribution ((a) ; (b) ).

(a)

(b)

Figure 9: Local entropy generation number distribution ((a) ; (b) ).

Figures 6–9 also show that the effect of on local entropy generation number is higher at the inner fluid surface and this effect vanishes as the outer surface is reached. The total entropy generation variations with are shown in Figures 10–13 for different set of parameters. Figure 10 corresponds to Figure 6, Figure 11 corresponds to Figure 7 except one more case for is added, and Figure 12 corresponds to Figure 8 except one more case for is added. Figure 13 corresponds to Figure 9 except one more case for is added. It is shown in these figures that increasing Brinkman number results in an increase in the total entropy generation, increasing Biot number results in a decrease in the total entropy generation, decreasing internal heat generation yields in lower total entropy generation, and decreasing temperature ratio value results in an increase in the total entropy generation. Figures 10–13 show that there exists a minimum total entropy generation for all the set of parameters considered in this study and this minimum occurs at .

Figure 10: Effect of on dimensionless total entropy generation rate.

Figure 11: Effect of on dimensionless total entropy generation rate.

Figure 12: Effect of on dimensionless total entropy generation rate.

Figure 13: Effect of on dimensionless total entropy generation rate.

5. Conclusions

A mathematical model of a fluid with internal heat generation between two rotating cylinders has been developed for velocity, temperature, local entropy generation rate, and total entropy generation rate. It has been proved that the local and the total entropy generation rates have minimum values for all the parameters considered in this work. It is concluded that decreasing Brinkman number, increasing Biot number, decreasing the internal heat generation, and increasing the temperature ratio result in lower total entropy generation rate. It is believed that the results of this study will be of great importance in engineering applications.

Nomenclature

: | Cross-sectional area ()
: | Biot number
: | Brinkman number
: | Specific heat capacity ()
: | Constant  (15)
: | Constant  (16)
: | Constant  (17)
: | Eckert number
: | Cylinder length (m)
: | Entropy generation number
: | Dimensionless total entropy generation rate
: | Internal heat generation ()
: | Dimensionless internal heat generation
: | Prandtl number
: | Radial distance (m)
: | Dimensionless radial distance
: | Total entropy generation rate ()
: | Local entropy generation rate ()
: | Temperature ratio
: | Temperature (K)
: | Tangential velocity ()
: | Dimensionless velocity
: | Velocity ratio
: | Radius ratio.

Greek Letters

: | Convective heat transfer coefficient ()
: | Viscous dissipation function  (9)
: | Thermal conductivity ()
: | Dynamic viscosity ()
: | Dimensionless temperature
: | Angular velocity ().

Subscripts

: | Inner cylinder surface
: | Outer cylinder surface
: | Coolant at the outer cylinder surface.

References

M. D. Spasojević, M. R. Janković, and D. D. Djaković, “A new approach to entropy production minimization in diabatic distillation column with trays,” Thermal Science, vol. 14, no. 2, pp. 317–328, 2010.

M. Nourollahi, M. Farhadi, and K. Sedighi, “Numerical study of mixed convection and entropy generation in the poiseulle-benard channel in different angles,” Thermal Science, vol. 14, no. 2, pp. 329–340, 2010.

L. Kolsi, A. Abidi, M. N. Borjini, and H. B. Aïssia, “The effect of an external magnetic field on the entropy generation in three-dimensional natural convection,” Thermal Science, vol. 14, no. 2, pp. 341–352, 2010.

T. H. Ko, “Numerical investigation of laminar forced convection and entropy generation in a helical coil with constant wall heat flux,” Numerical Heat Transfer A, vol. 49, no. 3, pp. 257–278, 2006.

A. M. Houshang, M. Shahi, and F. Talebi, “Entropy generation due to natural convection in a partially open cavity with a thin heat source subjected to a nanofluid,” Numerical Heat Transfer A, vol. 61, no. 4, pp. 283–305, 2012.

H. F. Oztop, “Effective parameters on second law analysis for semicircular ducts in laminar flow and constant wall heat flux,” International Communications in Heat and Mass Transfer, vol. 32, no. 1-2, pp. 266–274, 2005.

Y. Varol, H. F. Oztop, and A. Koca, “Entropy generation due to conjugate natural convection in enclosures bounded by vertical solid walls with different thicknesses,” International Communications in Heat and Mass Transfer, vol. 35, no. 5, pp. 648–656, 2008.

H. F. Öztop, I. Daǧtekin, and A. Z. Şahin, “Second law analysis of fully developed laminar flow for rectangular ducts with semicircular ends,” International Communications in Heat and Mass Transfer, vol. 36, no. 7, pp. 725–730, 2009.

B. S. Yilbas, “Entropy analysis of concentric annuli with rotating outer cylinder,” Exergy an International Journal, vol. 1, pp. 60–66, 2001.

H. Hamakawa, H. Mori, M. Iino, M. Hori, M. Yamasaki, and T. Setoguchi, “Experimental study of heating fluid between two concentric cylinders with cavities,” Journal of Thermal Science, vol. 17, no. 2, pp. 175–180, 2008.

M. Mirzazadeh, A. Shafaei, and F. Rashidi, “Entropy analysis for non-linear viscoelastic fluid in concentric rotating cylinders,” International Journal of Thermal Sciences, vol. 47, no. 12, pp. 1701–1711, 2008.

N. A. Roschina, A. V. Uvarov, and A. I. Osipov, “Natural convection in an annulus between coaxial horizontal cylinders with internal heat generation,” International Journal of Heat and Mass Transfer, vol. 48, no. 21-22, pp. 4518–4525, 2005.

M. S. Tshehla and O. D. Makinde, “Analysis of entropy generation in a variable viscosity fluid flow between two concentric pipes with convective cooling at the surface,” International Journal of Physical Sciences, vol. 6, no. 25, pp. 6053–6060, 2011.

S. Mahmud and R. A. Fraser, “Analysis of entropy generation inside concentric cylindrical annuli with relative rotation,” International Journal of Thermal Sciences, vol. 42, no. 5, pp. 513–521, 2003.

M. E. H. Assad, “Effect of maximum and minimum heat capacity rate on entropy generation in a heat exchanger,” International Journal of Energy Research, vol. 34, no. 14, pp. 1302–1308, 2010.

O. D. Makinde and A. Aziz, “Second law analysis for a variable viscosity plane Poiseuille flow with asymmetric convective cooling,” Computers and Mathematics with Applications, vol. 60, no. 11, pp. 3012–3019, 2010.

E. M. A. Mokheimer and A. H. Al-Salim, “Entropy generation due to mixed convection between vertical parallel plates under isothermal boundary conditions,” International Journal of Exergy, vol. 6, no. 5, pp. 671–697, 2009.

M. El Haj Assad, “Entropy generation analysis in a slab with non-uniform heat generation subjected to convection cooling,” International Journal of Exergy, vol. 9, no. 3, pp. 355–369, 2011.

A. Bejan, Entropy Generation Minimization, CRC Press, Boca Raton, Fla, USA, 1996.