Entropy generation in hydromagnetic and thermal boundary layer flow due to radial stretching sheet with Newtonian heating

Document Type : Full Lenght Research Article

Authors

1 University of Gour Banga, Malda 732 103, WB, India

2 Department of Applied Mathematics, Vidyasagar University, Midnapore 721 1102, India

3 Faculty of Military Science, Stellenbosch University, Private Bag X2, Saldanha 7395, South Africa

Abstract

The entropy generation during hydromagnetic boundary layer flow of a viscous incompressible electrically conducting fluid due to radial stretching sheet with Newtonian heating in the presence of a transverse magnetic field and the thermal radiation has been analyzed. The governing equations are then solved numerically by using the fourth order Runge-Kutta method with shooting technique. The effects of the pertinent parameters on the fluid velocity, temperature, entropy generation number, Bejan number, as well as the shear stress at the surface of the sheet are discussed graphically and quantitatively. It is examined that because the presence of a magnetic field, the entropy generation in a thermal system can be controlled and reduced.

Keywords

Main Subjects


 

 

Entropy generation in hydromagnetic and thermal boundary layer flow due to radial stretching sheet with Newtonian heating

 

 

Sanatan Das1* , Rabindra Nath Jana and  Oluwole Daniel Makinde3

 

I*Department of Mathematics, University of Gour Banga, Malda 732 103, India

2Department of Applied Mathematics, Vidyasagar University, Midnapore 721 102, India

3Faculty of Military Science, Stellenbosch University, Private Bag X2, Saldanha 7395, South Africa

 

Corresponding author e-mail: tutusanasd@yahoo.co.in

 

 

 


Journal of Heat and Mass Transfer Research

 

Journal homepage: http://jhmtr.journals.semnan.ac.ir

 

A B S T R A C T

The entropy generation during hydromagnetic boundary layer flow of a viscous incompressible electrically conducting fluid due to radial stretching sheet with Newtonian heating in the presence of a transverse magnetic field and the thermal radiation has been analyzed. The governing equations are then solved numerically by using the fourth order Runge-Kutta method with shooting technique. The effects of the pertinent parameters on the fluid velocity, temperature, entropy generation number, Bejan number, as well as the shear stress at the surface of the sheet are discussed graphically and quantitatively. It is examined that because the presence of a magnetic field, the entropy generation in a thermal system can be controlled and reduced.

 

 

© 2016 Published by Semnan University Press. All rights reserved.

 

 

PAPER INFO

History:

Received 3 January 2014

Received in revised form 4 June 2014

Accepted 29 September 2015

 

Keywords:

Hydromagnetic flow,

Newtonian heating, Entropy generation, radial stretching sheet

 

Journal of Heat and Mass Transfer Research 2(2015) 51-61

 

 

 

 

 

 


 

1. Introduction

 The boundary layer flow of a continuously stretching sheet has attracted considerable attention in recent years due to its numerous applications in industry. It occurs frequently in manufacturing involving the hot metal rolling, wire drawing, glass-fiber production, paper production, drawing of plastic films, and metal spinning, as well as metal and polymer extrusion processes. The boundary layer flow and the heat transfer over a continuous stretching surface have been discussed by a significant number of researchers [1-7] due to its wide applications in industrial and engineering processes. However, less attention has been paid to the boundary layer flow over a radial stretching sheet. Hayat et al. [8] have studied an axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet. Sahoo and Sharma [9] have investigated the MHD flow and heat transfer from a continuous surface in a uniform free stream of a non-Newtonian fluid. Sajid et al. [10] have presented an unsteady axisymmetric flow and heat transfer over a radial stretching sheet. Ariel [11] has presented the extended homotopy perturbation method and computation of flow past a stretching sheet. Sahoo [12] has examined the effects of the partial slip on axisymmetric flow of an electrically conducting viscoelastic fluid past a stretching sheet. Effects of slip, viscous dissipation and Joule heating on the MHD flow and heat transfer of a second grade fluid past a radial stretching sheet has been investigated by Sahoo [13].

The second law of thermodynamics is more reliable than the first law of thermodynamics  due to the limitation of the efficiency of the first law in heat transfer engineering systems (Oztop and Al-Salem [14]). In order to access the best design of the thermal systems, one can employ the second law of thermodynamics by minimizing the irreversibility. Entropy generation is a criterion of the destruction of the available system work. The evaluation of the entropy generation is carried out to improve system performance. Heat transfer, mass transfer, viscous dissipation, etc. can be used as sources of entropy generation. Entropy generation can be used as a quantitative measure of irreversibilities that are associated with a process, because of this fact that the greater the entropy generation indicates the greater the extent of irreversibilities. In many engineering and industrial processes, entropy production, destroys the available energy in the system. It is therefore imperative to determine the rate of entropy generation in a system, in order to optimize energy in the system for efficient operation in the system. According to the second law of thermodynamics, all the flows and heat transfer processes undergo changes that are irreversible. These irreversible changes are mostly caused by the energy losses during the processes. Although measures can be taken to reduce these irreversible effects, it is impossible to recover all of the lost energy. This process causes the entropy of the system to increase. Due to this, the entropy generation rate is used as a standard metric to study the irreversibile effects. This method was proposed by Bejan [15,16].

In recent years, many papers have been published on the applications and entropy generation rates of the second law of thermodynamics. Odat et al. [17] have explored the entropy generation effects in the laminar flow past a flat plate under the influence of the magnetic field and found that the entropy generation rate increases with the magnetic field intensity. This study has revealed that the magnetic field is one of the causes that are responsible for the entropy production in the systems. Saouli and Aiboud-Saouli [18] investigated the second law analysis of laminar falling liquid film along an inclined heated plate. Esfahani and Jafarian [19] have presented the entropy generation analysis of a flat plate boundary layer with different solution techniques. An irreversibility analysis for gravity driven non-Newtonian liquid film along an inclined isothermal plate has been presented by Makinde [20]. Arikoglu et al. [21] have examined the effect of slip on entropy generation in a single rotating disk in MHD flow. Aiboud and Saouli [22] have illustrated the application of the second law analysis of thermodynamics to viscoelastic magneto-hydrodynamic flow over a stretching surface analytically by using Kummer's functions. Makinde [23] has conducted a thermodynamic analysis on a gravity-driven liquid film along an inclined heated plate. He assumed viscosity to be a variable quantity and considered the convective cooling effect. Makinde [24] has studied a variable viscosity boundary layer flow over a flat plate under the effects of thermal radiation and Newtonian heating, as well as explored entropy generation effects in this flow. Makinde [25] has examined the entropy generation on magnetohydrodynamic flow and heat transfer over a flat plate with a convective boundary condition. The effect of viscous dissipation and thermal radiation on entropy generation in Blasius flow has been displayed numerically by Butt et al. [26]. Rashidi et al. [27] have analyzed the entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. The results showed that as the thermal radiation parameter increases, the generated entropy decreases. Butt and Ali [28] have illustrated the effects of magnetic field on entropy generation in the flow and heat transfer due to a radially stretching surface. Butt and Ali [29,30] have carried out the entropy analysis of flow and heat transfer caused by a moving surface. Recently, Butt and Ali [31] have presented the entropy analysis of magnetohydrodynamic flow and heat transfer over a convectively heated radially stretching surface.

The aim of the present paper is to explore the entropy generation in MHD boundary layer flow of a viscous incompressible electrically conducting fluid due to a radial stretching sheet in the presence of a transverse magnetic field by taking the convective boundary condition into account. The partial differential equations governing the flow are reduced to nonlinear ordinary differential equations, which are solved numerically by shooting technique by using fourth order Runge-Kutta method. A systematic study of the effects of the various pertinent parameters on the flow and heat transfer characteristic is carried out with the help of graphs and table.

2. Mathematical formulation

Consider a steady two-dimensional boundary layer flow due to the stretching of the sheet along the radial direction with the velocity , where  is a constant, as shown in Fig. 1. The sheet is located in the plane  and the fluid is confined to the region . The flow is due to the stretching of the sheet. A uniform magnetic field of strength  is applied perpendicularly to the sheet, i.e. in the -direction. The lower surface of the sheet is heated by convection from a hot fluid at temperature  which provides a heat transfer coefficient  while the temperature of the ambient cold fluid is . The cold fluid on the upper side of the sheet is assumed to be an electrically conducting Newtonian fluid with constant fluid property. There is no external electric field. Also, it is assumed that the magnetic Reynolds number is small enough so that the induced magnetic field can be neglected.

An order-of-magnitude analysis of the momentum equation (normal to the sheet) by using the usual boundary layer approximations shows that the pressure gradient is constant. Thus neglecting the pressure gradient, the continuity, momentum and energy equations in a viscous MHD incompressible boundary layer flow can be written, respectively, as follows

 

(1)

 

(2)

 

(3)

 

Where  and  are the velocity components along the  and -directions, respectively,  the temperature of the fluid,  the electric conductivity of the fluid,  the kinematic viscosity,  the fluid density,  the thermal conductivity,  the specific heat at constant pressure and  the radiative heat flux.

The appropriate boundary conditions are

 

 

 

(4)

Where  is the thermal conductivity.

 

For an optically thick fluid, in addition to emission there is also self absorbed and usually the absorption coefficient is large and dependent on the wavelength so that we can apply Rosseland approximation for radiative flux. The Rosseland approximation [32] applies to optically thick medium and taking into account the model proposed by Magyari and Pantokratoras [33], the net radiation heat flux  is given by

 

 

(5)

 

Where  is the Stefan-Boltzmann constant and  the Rosseland mean absorption coefficient. The Rosseland mean absorption coefficients for the optically thick media are expressed as function of the thermodynamic properties of the media. It is assumed that the temperature difference between the fluid temperature and the free stream temperature  is small, so that the term  may be expressed as a linear function of temperature. This is done by expanding  in a Taylor series about a free stream temperature  as follows:

 

(6)

Neglecting higher-order terms in the equation (6) beyond the first order in we get

 

(7)

 

The use of the equation (5) and (7), the equation (3) becomes

 

 

 

(8)

The following similarity variables are introduced (Butt and Ali  [28]):

 

 

 

(9)

 

 

 

Fig.1: Geometry of the problem

 

where  is the independent similarity variable,  the dimensionless stream function and  the dimensionless temperature.

The use of equation (9), equations (2) and (8) reduce to

 

(10)

 

(11)

Where  is the magnetic parameter,  the radiation parameter, and  the Prandtl number which measures the ratio of momentum diffusivity to the thermal diffusivity.

The thermal radiation is quite significant and the quality of the final product can be controlled by the control of cooling rate via the radiation parameter. In polymer industry, the thermal radiation effect may play an important role in the control of heat transfer process if the process is directed in a thermally controlled environment.

The corresponding boundary conditions are

 

 

(12)

Where  is the surface convection parameter or so-called Biot number. When , the convective boundary condition reduces to a uniform surface temperature boundary condition. For uniformly heated ( ) radial stretching sheet and in the absence of thermal radiation ( ) the present problem reduces to the Butt and Ali [28].

3. Numerical method for solution

 

The governing non-linear ordinary differential equations (10) and (11) cannot be solved analytically. This set of equations with the boundary conditions (12) must be solved numerically by applying most efficient fourth-order Runge-Kutta integration scheme with shooting algorithm. Equations (10) and (11) and boundary conditions (12) are reduced to a set of simultaneous first order differential equations by setting , , ,  and  as follows

 

 

 

 

(13)

With the boundary conditions

 

(14)

Where  and  are determined such that  and . The essence of this method is to reduce the boundary value problem to an initial value problem and then to use the shooting numerical technique to guess  and  until the boundary conditions  and  are satisfied. It is important to note that the infinity ( ) in the above equations represents the boundary layer thickness. We compare the calculated values for  and  at  with the given boundary condition  and  and adjust the estimated values,  and , to give a better approximation for the solution. The accuracy of the assumed missing initial condition is checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. The step size  is used to obtain the numerical solution with seven-decimal place accuracy as the criterion of convergence. The numerical computations are done by MATLAB built-in routine. The method is adequately explained in literature and it has second order convergence, unconditionally stable. It gives accurate result for boundary layer equations. In the present study, a uniform grid is used which is concentrated toward the wall.

4. Results and discussion

 In order to gain a clear physical insight of the problem, we have discussed the effects of different values of magnetic parameter , radiation parameter , Prandtl number  and the Biot number on the velocity, temperature and shear stress at the surface of the sheet. The default values of the other parameters are mentioned in the description of the respected figures. Fig.2 shows that the fluid velocity  decreases with an increase in magnetic parameter . A drag-like Lorentz force is created by the application of the transverse magnetic field of the electrically conducting fluid. This force has the tendency to slow down the fluid flow. Fig.3 shows that the fluid temperature  increases by increasing the magnetic parameter  in the boundary layer region and the thermal boundary layer thickness increases as the magnetic field becomes stronger. It is important to notice that large resistances on the fluid particles, which cause heat to be generated in the fluid, as the transverse applied magnetic field increases. It  can be observed from Fig.4 that the fluid temperature increases by increasing values of radiation parameter . The increase in radiation parameter means the release of heat energy from the flow region and so the fluid temperature decreases. A decrease in the values of  for given  and  means a decrease in the Rosseland radiation absorptivity . Since the divergence of the radiative heat flux  increases,  decreases which in turn causes to increase the rate of radiative heat transfer to the fluid and hence the fluid temperature increases. In the presence of thermal radiation, the thermal boundary layer always found to be increased. This means that the thermal boundary layer increases and there is more uniform temperature distribution across the boundary layer.

The effect of Prandtl number Pr on the heat transfer process is shown in the Fig.5. This figure shows that an increase in Prandtl number number Pr results in a decrease in the temperature distribution because, the thermal boundary layer thickness decreases  by increasing  the Prandtl number Pr. In  summary, an increase in the Prandtl number means slow rate of the thermal diffusion. The graph also shows that the wall temperature decreases, as the values of Prandtl number Pr increase. It is also found that the thermal boundary layer thickness reduces as Pr increases. Fig.6 demonstrates the effects of the Biot number  on fluid temperature. It  can be seen that the temperature profiles within the boundary layer increase  by increasing  in Biot number, leading to an increase in the thermal boundary layer thickness. Biot number is the ratio of the hot fluid side convection resistance to the cold fluid side convection resistance of a surface. For fixed cold fluid properties, Biot number  is directly proportional to the heat transfer coefficient   associated with the hot fluid. The thermal resistance on the hot fluid side is inversely proportional to . Thus, as  increases, the hot fluid side convection resistance decreases and consequently, the surface temperature increases. It is also noticed that for large values of , i.e. , the temperature profile attains its maximum value 1; thus the convective boundary condition becomes the prescribed surface temperature case. The thermal boundary layer thickness increases as  increases.

 

 

Fig.2: Velocity profiles for different

 

 

 

Fig.3: Temperature profiles for different  when ,    and

 

 

 

Fig.4: Temperature profiles for different  when ,    and

 

 

Fig.5: Temperature profiles for different  when ,    and

 

 

Fig.6: Temperature profiles for different  when ,     and

At the engineering point of view, the shear stress at the surface of the sheet is an important characteristic in the heat transfer studies, since it is directly related to the heat transfer coefficients. The increased shear stress is generally a disadvantage in the technical applications, while the increased heat transfer can be exploited in some applications such as heat exchangers, but should be avoided in other applications such as a gas turbine, for instance. The numerical values of the surface temperature  and the shear stress  at the surface of the sheet  are entered in the Table 1 for several values of , , Pr and . It  can be seen from the Table 1 that the surface temperature  increases by increasing values of either  or  or  while it decreases  by increasing  in . This is  because of the fact that the Biot number  is directly proportional to the heat transfer coefficient associated with hot fluid. The thermal resistance on the hot fluid side is inversely proportional to . Thus, as  increases, the hot fluid side convection resistance decrease; as a result, the surface temperature  increases. The shear stress  increases by increasing values of .

4.1. Entropy generation

Non-equilibrium conditions happen due to the exchange of the energy and momentum within the fluid and at solid boundaries that result in entropy generation. According to Woods [34], the local volumetric rate of entropy generation for a viscous incompressible conducting fluid in the presence of a magnetic field is given by

 

 

 

(15)

 

The first term in equation (15) is the irreversibility due to the heat transfer, the second term is the entropy generation due to viscous dissipation and the third term is a local entropy generation due to the effect of magnetic field (Joule

The dimensionless entropy generation number may be defined by the following relationship:

 

 

(16)

The use of (9), the entropy generation number in dimensionless form is

 

 

 

(17)

Where  is the Reynolds number,   is the Brinkmann number and  the non-dimensional temperature difference.

The entropy generation number  can be written as a summation of the entropy generation due to heat transfer denoted by  and the entropy generation due to fluid friction with magnetic field denoted by  given as

 

 

(18)

In order to obtain an idea of whether entropy generation due to the heat transfer dominates over entropy generation due to the fluid friction and the magnetic field, or vice versa, the Bejan number  is defined to be the ratio of entropy generation due to the heat transfer to the entropy generation number [35].

 

 

 

(19)

Where  is the irreversibility ratio. Heat transfer dominates for  and fluid friction with magnetic effects dominates when . The contribution of both heat transfer and fluid friction to entropy generation is equal when . The Bejan number  takes the values between 0 and 1 (see Cimpean et al. [36]). The value of  is the limit at which the heat transfer irreversibility dominates,  is the opposite limit at which the irreversibility is dominated by the combined effects of fluid friction and magnetic field and  is the case in which the heat transfer and fluid friction with magnetic field entropy production rates are equal. Further, the behavior of the Bejan number  is studied for the optimum values of the parameters at which the entropy generation takes its minimum.

 

The influences of the different governing parameters on entropy generation and Bejan number are presented in Figs. 7-14. It  can be seen from Fig.7 that the entropy generation number  increases  by increasing  in magnetic parameter . An increase in the magnetic field intensity causes an increase in the entropy generation. It reveals that the magnetic field is a source of entropy generation in addition to the fluid friction and the heat transfer. Also, it  can be seen that the entropy generation effects are prominent on the surface of the stretching sheet and in the region close to it. However, in the free stream region the entropy effects are negligible. This implies that in order to control the entropy which is generated in boundary layer flow, the value of the magnetic parameter should be reduced, which is an issue of interest in nuclear-MHD propulsion. Fig.8 shows that the entropy generation number  increases by increasing Reynolds number  due to higher heat transfer rates at the sheet surface. An increase in Reynolds number , the entropy effects due to heat transfer entropy effects become prominent and fluid friction and magnetic field are lessen near the stretching sheet surface. However, as the distance increases from the sheet surface, these effects are negligible. Fig.9 illustrates the effects of the Biot number  on the entropy generation. Increasing  enhances the entropy generation. It means that the convective surface boundary condition acts as a strong source of irreversibility. Therefore, the entropy can be minimized by reducing the convection through the boundary. The closeness of the curves in Fig.9 can be attributed that the entropy effects are dominated by the entropy effects due to fluid friction and magnetic field. It can be observed from Fig.10 that entropy generation number increases by increasing the group parameter  due to the viscous heating effects. An increase in the values of the group parameter  due to the combined effects of viscous heating and temperature difference yields a higher entropy generation number.

Fig.11 shows that the Bejan number  increases by increasing the values of magnetic parameter . For large values of , the entropy effects due to fluid friction and magnetic field are fully dominated by heat transfer entropy effects near the sheet surface. Fig.12 illustrates  the effects of the radiation parameter  on the Bejan number . By increasing values of , the entropy effects due to heat transfer become strong and hence Bejan number  increases. This is due to the fact that the effect of the radiation parameter  is to enhance the temperature significantly in the flow region. The increase in radiation parameter means the generation of heat energy in the flow region. An increase in Biot number  leads to increase the Bejan number  as shown in Fig.13. Also, an increase in the values of the Biot number results in an increase in the dominant effect of heat transfer irreversibility at the sheet surface. This means that the sheet surface acts as a strong source of irreversibility. Fig.14 reveals that the Bejan number  decreases  by increasing  group parameter  . This is quite true as higher values of , which increase the magnitude of fluid friction with magnetic field irreversibility  but has no effect on the heat transfer irreversibility , increases the values of  leading to lower Bejan number. The group parameter is an important dimensionless number for irreversibility analysis. It determines the relative importance of viscous effects to that of temperature gradient entropy generation. The graphs of the Bejan number are useful to obtain an idea on whether heat transfer irreversibility dominates fluid friction irreversibility or vice versa.

 

5. Conclusion

Entropy generation analysis in hydromagnetic boundary layer flow of a viscous incompressible electrically conducting fluid due to a radially stretching sheet with Newtonian heating in the presence of a transverse magnetic field has been carried out. The velocity and temperature profiles are obtained numerically and used to compute the entropy generation number. The effects of the pertinent parameters on velocity and temperature profiles are presented graphically.

The influences of the same parameters and the dimensionless group parameter on the entropy generation rate and Bejan number are also discussed.

From the results the following conclusions could be drawn:

  • · The magnetic field retards the fluid velocity while it causes to increase the fluid temperature.
  • · The fluid temperature decreases by increasing the values of the radiation parameter, leading to a decrease in the thermal boundary layer thickness.
  • · It is also found that the thermal boundary layer thickness reduces as Prandtl number increases.
  • · The fluid temperature increases by increasing values of Biot number, leading to an increase in the thermal boundary layer thickness.
  • · The surface of the sheet acts as strong source of entropy and the heat transfer irreversibility.

 

 

 

Table 1. The shear stress  and the temperature  at the sheet .

 

 

   

 

 

 0.5

 

0.72

0.1

1.57758

 0.08345

1

 

0.72

0.1

1.70391

 0.08358

2

 

0.72

0.1

1.94118

 0.08384

 

 

1

0.1

 

 

 

 0.08018

1

0.5

 

 

 

 0.08283

1

1

 

 

 

 0.08474

 

 

 

 

0.72

 

 

 0.08358

 

 

2

 

 

 0.07279

 

 

 

 

 0.06597

 

 

 

 

 

0.1

 

 0.08358

 

 

 

1

 

 0.47701

 

 

 

5

 

 0.82016

                           


 

Fig.7:  for different  when ,  and

 

Fig.8:  for different  when ,  and

 

Fig.9:  for different  when ,  and

 

Fig.10:  for different  when ,  and

 

Fig.11: Bejan number  for different  when ,  and

 

Fig.12: Bejan number  for different  when ,  and

 

Fig.13: Bejan number  for different  when ,  and

 

Fig.14: Bejan number  for different  when  and

  • · The optimum design and efficient performance of a flow system or a thermally designed system can be improved by choosing the appropriate values of the physical parameters. This will be enabled to reduce the effects of entropy generated within the system.

 

Acknowledgement

The authors are thankful to the honourable reviewers for constructive suggestions which helped a lot in improving the manuscript.

 

Nomenclature

 

   

 a constant

 

 applied magnetic field

 

 Bejan number

 

 Biot number

 

 Brinkmann number

 

 specific heat at constant pressure

 

 volumetric rate of entropy generation

 

 non-dimensional steam function

 

 first order derivative with respect to

 

 second order derivative with respect to

 

 third order derivative with respect to

 

 heat transfer coefficient

 

 thermal conductivity

 

 Rosseland mean absorption coefficient

 

 magnetic parameter

 

 entropy generation due to heat transfer

 

 entropy generation due to fluid friction and magnetic field

 

 entropy generation number

 

 Prandtl number

 

 radiative heat flux

 

 radiation parameter

 

 Reynolds number

 

 fluid temperature

 

 hot fluid temperature

 

 free stream temperature

,  

 velocity components in  and -directions

 

 Stretching velocity of the sheet in radial direction

 

 Greek symbols

 

 similarity variable

 

 kinematic viscosity

 

 irreversibility distribution ratio

 

 the fluid density

 

 Stefan-Boltzman constant

 

 non-dimensional temperature difference

 

 non-dimensional temperature

 

 first order derivative with respect to

 

 second order derivative with respect to

References

  [1] Carragher, P., Crane, L.J.(1982). Heat transfer on a continuous stretching sheet. Z. Angew Math. Mech. 62, 564-565.

  [2] Wang, C.Y. (1990).  Liquid film on an unsteady stretching sheet. Q. Appl. Math. 48, 601-610.

  [3] Ariel, P. D. (2001). Axisymmetric flow of a second grade fluid past a stretching sheet. Int. J. Eng. Sci. 39(5), 529-553.

  [4] Chiam, T.C. (1994). Stagnation point flow towards a stretching plate. J. Phys. Soc. Jpn. 63, 2443-2444.

  [5] Liu, I. C. (2004). Flow and heat transfer of an electrically conducting fluid of second grade over a stretching sheet subject to a transverse magnetic field. Int. J. Heat Mass Transfer 47(19–20), 4427-4437.

  [6] Wang, C.Y . (2007). Natural convection on a vertical radially stretching sheet. Math. Anal. App. 332, 877- 883.

  [7] Ariel, P.D. (2007). Axisymmetric flow due to a stretching sheet with partial slip. Comput. Math. Appl.  54, 1169-1183.

  [8] Hayat, T., Sajid, M. (2007). Analytic solution for axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet. Int. J. Heat Mass Transf  50, 75-84.

  [9] Sahoo, B., Sharma, H. G. (2007).  MHD flow and heat transfer from a continuous surface in a uniform free stream of a non-Newtonian fluid. Appl. Math. Mech.-Eng. Ed. 28(11),1467-1477.

 [10] Sajid, M., Ahmad, I., Hayat, T., Ayub, M. (2008). Series solution for unsteady axisymmetric flow and heat transfer over a radially stretching sheet. Commun. Nonlinear Sci. Numer. Simul. 13(10), 2193-2202.

 [11] Ariel, P.D. (2009). Extended homotopy perturbation method and computation of flow past a stretching sheet. Comput. Math. Appl. 58, 2402-2409.

 [12] Sahoo, B. (2009). Effects of partial slip on axisymmetric flow of an electrically conducting viscoelastic fluid past a stretching sheet. Cent. Eur. J. Phys., DOI; 10.2478/s11534-009-0105-x.

 [13] Sahoo, B. (2010).  Effects of slip, viscous dissipation and joule heating on the MHD flow and heat transfer of a second grade fluid past a radially stretching sheet. Appl. Math. Mech.- Eng. Ed. 31, 159-173.

 [14] Oztop, H. F., Al-Salem, K. A. (2012).  review on entropy generation in natural and mixed convection heat transfer for energy systems, Renewable Sustainable Energy Rev. 16 (1), 911-920.

 [15] Bejan, A. (1980). Second law analysis in heat transfer. Energy Int. J. 5, 721-732.

 [16] Bejan, A . (1994). Entropy Generation Through Heat and Fluid Flow, Wiley, Canada.

 [17] Odat, M. Q. A., Damseh, R. A., Nimr, M. A. A. (2004). Effect of magnetic field on entropy generation due to laminar forced convection past a horizontal flat plate. Entropy 4, 293-303.

 [18] Saouli, S., Aiboud-Saouli, S. (2004). Second law analysis of laminar falling liquid film along an inclined heated plate. Int. Comm. Heat Mass Transfer 31(6), 879-886.

 [19] Esfahani, J. A., Jafarian, M. M. (2005).  Entropy generation analysis of a flat plate boundary layer with various solution methods. Scientia Iranica 12, 233-240.

 [20] Makinde, O. D. (2006). Irreversibility analysis for gravitry driven non-Newtonian liquid film along an inclined isothermal plate. Physica Scripta 74, 642-645.

 [21] Arikoglu, A., Ozkol. I., Komurgoz, G. (2008).Effect of slip on entropy generation in a single rotating disk in MHD flow. Appl. Energy 85(12), 1225-1236.

 [22] Adboud, S., Saouli, S. (2010).Entropy analysis for viscoelastic magneto hydrodynamic flow over a stretching surface. Int. J. Non Linear Mech. 45 (5), 482-489.

 [23] Makinde, O. D. (2010).Thermodynamic second law analysis for a gravity driven variable viscosity liquid film along an inclined heated plate with convective cooling. J. Mech. Sci. Tech. 24 (4), 899-908.

 [24] Makinde, O. D. (2011).Second law analysis for variable viscosity hydromagnetic boundary layer flow with thermal radiation and Newtonian heating, Entropy 13, 1446-1464.

 [25] Makinde, O. D. (2012).Entropy analysis for MHD boundary layer flow and heat transfer over a flat plate with a convective surface boundary condition. Int. J. Exergy 10 (2), 142 - 154.

 [26] Butt, A.S., Munawar, S., Ali, A., Mehmood, A. (2012).Entropy generation in the Blasius flow under thermal radiation. Phys. Scr. 85 (3), 035008.

 [27] Rashidi, S., Abelman, N., Freidooni, M. (2013). Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid, Int. J. Heat and Mass Transfer 62, 515-525.

 [28] Butt, A.S., Ali, A. (2013). Effects of magnetic field on entropy generation in flow and heat transfer due to a radially stretching surface. Chin. Phys. Lett. 30(2), 024701.

 [29] Butt, A.S., Ali, A.( 2014). Entropy analysis of flow and heat transfer caused by a moving plate with thermal radiation. J. Mech. Sci. Tech. 28 (1), 343-348.

 [30] Butt, A.S., Ali, A . A computational study of entropy generation in magnetohydrodynamic flow and heat transfer over an unsteady stretching permeable sheet. Eur. Phys. J. Plus 129(3), 1-13.

 [31] Butt, A.S., Ali, A.( 2014). Entropy analysis of magnetohydrodynamic flow and heat transfer over a convectively heated radially stretching surface. J. Taiwan Institute of Chem. Eng. 45, 1197-1203.

 [32] Rosseland, S. (1931). Astrophysik und atom-theoretische Grundlagen, Springer-Verlag, Berlin.

 [33] Magyari, E., Pantokratoras, A. (2011). Note on the effect of thermal radiation in the linearized Rosseland approximation on the heat transfer characteristics of various boundary layer flows. Int. Commun. Heat and Mass Trans. 38, 554-556.

 [34] Woods, L. C. (1975). Thermodynamics of Fluid Systems. Oxford University Press, Oxford, UK.

 [35] Paoletti, S, Rispoli, F, Sciubba, E.( 1989). Calculation of exergetic losses in compact heat exchanger passages. ASME AES 10, 21-29.

 [36] Cimpean, D., Lungu, N., Pop, I. (2008). A problem of entropy generation in a channel filled with a porous medium. Creative Math. Inf. 17, 357-362.