Entropy Generation in Viscoelastic Fluid Over a Stretching Surface

Due to the increasing importance in processing industries and elsewhere when materials whose flow behavior cannot be characterized by Newtonian relationships, a new stage in the evolution of fluid dynamics theory is in progress. An intensive effort, both theoretical and experimental, has been devoted to problems of non-Newtonian fluids. The study of MHD flow of viscoelastic fluids over a continuously moving surface has wide range of applications in technological and manufacturing processes in industries. This concerns the production of synthetic sheets, aerodynamic extrusion of plastic sheets, cooling of metallic plates, etc. (Crane, 1970) considered the laminar boundary layer flow of a Newtonian fluid caused by a flat elastic sheet whose velocity varies linearly with the distance from the fixed point of the sheet. (Chang, 1989; Rajagopal et al., 1984) presented an analysis on flow of viscoelastic fluid over stretching sheet. Heat transfer cases of these studies have been considered by (Dandapat & Gupta, 1989, Vajravelu & Rollins, 1991), while flow of viscoelastic fluid over a stretching surface under the influence of uniform magnetic field has been investigated by (Andersson, 1992). Thereafter, a series of studies on heat transfer effects on viscoelastic fluid have been made by many authors under different physical situations including (Abel et al., 2002, Bhattacharya et al., 1998, Datti et al., 2004, Idrees & Abel, 1996, Lawrence & Rao, 1992, Prasad et al., 2000, 2002). (Khan & Sanjayanand, 2005) have derived similarity solution of viscoelastic boundary layer flow and heat transfer over an exponential stretching surface. (Cortell, 2006) have studied flow and heat transfer of a viscoelastic fluid over stretching surface considering both constant sheet temperature and prescribed sheet temperature. (Abel et al., 2007) carried out a study of viscoelastic boundary layer flow and heat transfer over a stretching surface in the presence of non-uniform heat source and viscous dissipation considering prescribed surface temperature and prescribed surface heat flux. (Khan, 2006) studied the case of the boundary layer problem on heat transfer in a viscoelastic boundary layer fluid flow over a non-isothermal porous sheet, taking into account the effect a continuous suction/blowing of the fluid, through the porous boundary. The effects of a transverse magnetic field and electric field on momentum and heat transfer characteristics in viscoelastic fluid over a stretching sheet taking into account viscous dissipation and ohmic dissipation is presented by (Abel et al., 2008). (Hsiao, 2007) studied


Introduction
Due to the increasing importance in processing industries and elsewhere when materials whose flow behavior cannot be characterized by Newtonian relationships, a new stage in the evolution of fluid dynamics theory is in progress.An intensive effort, both theoretical and experimental, has been devoted to problems of non-Newtonian fluids.The study of MHD flow of viscoelastic fluids over a continuously moving surface has wide range of applications in technological and manufacturing processes in industries.This concerns the production of synthetic sheets, aerodynamic extrusion of plastic sheets, cooling of metallic plates, etc. (Crane, 1970) considered the laminar boundary layer flow of a Newtonian fluid caused by a flat elastic sheet whose velocity varies linearly with the distance from the fixed point of the sheet.(Chang, 1989;Rajagopal et al., 1984) presented an analysis on flow of viscoelastic fluid over stretching sheet.Heat transfer cases of these studies have been considered by (Dandapat & Gupta, 1989, Vajravelu & Rollins, 1991), while flow of viscoelastic fluid over a stretching surface under the influence of uniform magnetic field has been investigated by (Andersson, 1992).Thereafter, a series of studies on heat transfer effects on viscoelastic fluid have been made by many authors under different physical situations including (Abel et al., 2002, Bhattacharya et al., 1998, Datti et al., 2004, Idrees & Abel, 1996, Lawrence & Rao, 1992, Prasad et al., 2000, 2002).(Khan & Sanjayanand, 2005) have derived similarity solution of viscoelastic boundary layer flow and heat transfer over an exponential stretching surface.(Cortell, 2006) have studied flow and heat transfer of a viscoelastic fluid over stretching surface considering both constant sheet temperature and prescribed sheet temperature.(Abel et al., 2007) carried out a study of viscoelastic boundary layer flow and heat transfer over a stretching surface in the presence of non-uniform heat source and viscous dissipation considering prescribed surface temperature and prescribed surface heat flux.(Khan, 2006) studied the case of the boundary layer problem on heat transfer in a viscoelastic boundary layer fluid flow over a non-isothermal porous sheet, taking into account the effect a continuous suction/blowing of the fluid, through the porous boundary.The effects of a transverse magnetic field and electric field on momentum and heat transfer characteristics in viscoelastic fluid over a stretching sheet taking into account viscous dissipation and ohmic dissipation is presented by (Abel et al., 2008).(Hsiao, 2007) studied the conjugate heat transfer of mixed convection in the presence of radiative and viscous dissipation in viscoelastic fluid past a stretching sheet.The case of unsteady magnetohydrodynamic was carried out by (Abbas et al., 2008).Using Kummer's funcions, (Singh, 2008) carried out the study of heat source and radiation effects on magnetohydrodynamics flow of a viscoelastic fluid past a stretching sheet with prescribed power law surface heat flux.The effects of non-uniform heat source, viscous dissipation and thermal radiation on the flow and heat transfer in a viscoelastic fluid over a stretching surface was considered in (Prasad et al., 2010).The case of the heat transfer in magnetohydrodynamics flow of viscoelastic fluids over stretching sheet in the case of variable thermal conductivity and in the presence of non-uniform heat source and radiation is reported in (Abel & Mahesha, 2008).Using the homotopy analysis, (Hayat et al., 2008) looked at the hydrodynamic of three dimensional flow of viscoelastic fluid over a stretching surface.The investigation of biomagnetic flow of a non-Newtonian viscoelastic fluid over a stretching sheet under the influence of an applied magnetic field is done by (Misra & Shit, 2009).(Subhas et al., 2009) analysed the momentum and heat transfer characteristics in a hydromagnetic flow of viscoelastic liquid over a stretching sheet with non-uniform heat source.(Nandeppanavar et al., 2010) analysed the flow and heat transfer characteristics in a viscoelastic fluid flow in porous medium over a stretching surface with surface prescribed temperature and surface prescribed heat flux and including the effects of viscous dissipation.(Chen, 2010) studied the magneto-hydrodynamic flow and heat transfer characteristics viscoelastic fluid past a stretching surface, taking into account the effects of Joule and viscous dissipation, internal heat generation/absorption, work done due to deformation and thermal radiation.(Nandeppanavar et al., 2011) considered the heat transfer in viscoelastic boundary layer flow over a stretching sheet with thermal radiation and non-uniform heat source/sink in the presence of a magnetic field Although the forgoing research works have covered a wide range of problems involving the flow and heat transfer of viscoelastic fluid over stretching surface they have been restricted, from thermodynamic point of view, to only the first law analysis.The contemporary trend in the field of heat transfer and thermal design is the second law of thermodynamics analysis and its related concept of entropy generation minimization.Entropy generation is closely associated with thermodynamic irreversibility, which is encountered in all heat transfer processes.Different sources are responsible for generation of entropy such as heat transfer and viscous dissipation (Bejan, 1979(Bejan, , 1982)).The analysis of entropy generation rate in a circular duct with imposed heat flux at the wall and its extension to determine the optimum Reynolds number as function of the Prandtl number and the duty parameter were presented by (Bejan, 1979(Bejan, , 1996)).(Sahin, 1998) introduced the second law analysis to a viscous fluid in circular duct with isothermal boundary conditions.In another paper, (Sahin, 1999) presented the effect of variable viscosity on entropy generation rate for heated circular duct.A comparative study of entropy generation rate inside duct of different shapes and the determination of the optimum duct shape subjected to isothermal boundary condition were done by (Sahin, 1998).(Narusawa, 1998) gave an analytical and numerical analysis of the second law for flow and heat transfer inside a rectangular duct.In a more recent paper, (Mahmud & Fraser, 2002a, 2002b, 2003) applied the second law analysis to fundamental convective heat transfer problems and to non-Newtonian fluid flow through channel made of two parallel plates.The study of entropy generation in a falling liquid film along an inclined heated plate was carried out by (Saouli & Aïboud-Saouli, 2004).As far as the effect of a magnetic field on the entropy generation is Entropy Generation in Viscoelastic Fluid Over a Stretching Surface 201 concerned, (Mahmud et al., 2003) studied the case of mixed convection in a channel.The effects of magnetic field and viscous dissipation on entropy generation in a falling film and channel were studied by (Aïboud-saouli et al., 2006(Aïboud-saouli et al., , 2007)).The application of the second law analysis of thermodynamics to viscoelastic magnetohydrodynamic flow over a stretching surface was carried out by (Aïboud & Saouli 2010a, 2010b).The objective of this paper is to study the entropy generation in viscoelastic fluid over a stretching sheet with prescribed surface temperature in the presence of uniform transverse magnetic field.

Formulation of the problem
In two-dimensional Cartesian coordinate system   , x y we consider magneto-convection, steady, laminar, electrically conduction, boundary layer flow of a viscoelastic fluid caused by a stretching surface in the presence of a uniform transverse magnetic field and a heat source.The x -axis is taken in the direction of the main flow along the plate and the y -axis is normal to the plate with velocity components , uv in these directions.Under the usual boundary layer approximations, the flow is governed by the following equations: The constant is the viscoelastic parameter.
The boundary conditions are given by 0, , 0 The heat transfer governing boundary layer equation with temperature-dependent heat generation (absorption) is The relevant boundary conditions are

Analytical solution
The equation of continuity is satisfied if we choose a dimensionless stream function   Introducing the similarity transformations where Now let us seek a solution of Eq. ( 6) in the form which is satisfied by the following boundary conditions: On substituting (7) into (6) and using boundary conditions (10a) and (10b) the velocity components take the form Where is the viscoelastic parameter, and Defining the dimensionless temperature Entropy Generation in Viscoelastic Fluid Over a Stretching Surface 203 and using ( 9), ( 11), (12), Eq. ( 14) and the boundary conditions (5a) and (5b) can be written as       Pr 12 P r 0 And ( 16a) and (16b) transform to The solution of Eq. ( 18) satisfying (19a) and (19b) is given by The solution of (20) in terms of  is written as

Second law analysis
According to (Woods, 1975), the local volumetric rate of entropy generation in the presence of a magnetic field is given by Eq. ( 22) clearly shows contributions of three sources of entropy generation.The first term on the right-hand side of Eq. ( 22) is the entropy generation due to heat transfer across a finite temperature difference; the second term is the local entropy generation due to viscous dissipation, whereas the third term is the local entropy generation due to the effect of the magnetic field.
using Eq. ( 9), ( 21) and ( 22), the entropy generation number is given by where Re l and Br are respectively the Reynolds number and the Brinkman number. and Ha , are respectively the dimensionless temperature difference and the Hartman number.These number are given by the following relationships

Results and discussion
The flow and heat transfer in a viscoelastic fluid under the influence of a transverse uniform magnetic field has been solved analytically using Kummer's functions and analytic expressions of the velocity and temperature have been used to compute the entropy generation.The entropy generation number is higher for higher Prandtl number near the surface, but as  increases, the entropy generation number shows different variation.This is due to the fact that according to fig. 6, the temperature profiles decrease sharply with the increase of the Prandlt number.viscous effect.For a given  , the entropy generation number is higher for higher dimensionless group.This is due to the fact that for higher dimensionless group, the entropy generation numbers due to the fluid friction increase.

Conclusion
The velocity and temperature profiles are obtained analytically and used to compute the entropy generation number in viscoelastic magnetohydrodynamic flow over a stretching surface The effects of the magnetic parameter and the viscoelastic parameter on the longitudinal and transverse velocities are discussed.The influences of the number, the magnetic parameter and the heat source/sink parameter on the temperature profiles are presented.
As far as the entropy generation number is concerned, its dependence on the magnetic parameter, the Prandlt, the Reynolds, the Hartmann numbers and the dimensionless group are illustrated and analyzed.
From the results the following conclusions could be drawn: a.The velocities depend strongly on the magnetic and the viscoelastic parameters.b.The temperature varies significantly with the Prandlt number, the magnetic parameter and the heat source/sink parameter.c.The entropy generation increases with the increase of the Prandlt, the Reynolds, the Hartmann numbers and also with the magnetic parameter and the dimensionless group.d.The surface acts as a strong source of irreversibility.

Figs
Fig. 1.Effect of the magnetic parameter on the longitudinal velocity.

Fig. 5
Fig. 5 depicts the temperature profiles    as function of  for different values of the Prandtl number Pr .As it can be noticed,     decreases with  whatever is the value of the Prandtl number, for a fixed value of  , the temperature     decreases with an increase in Prandtl number which means that the thermal boundary layer is thinner for large Prandtl number.

Fig. 5 .
Fig. 5. Effect of the Prandtl number on the temperature.The temperature profiles     as function of  for different values of the magnetic Mn are plotted in fig.6.An increase in the magnetic parameter Mn results in an increase of the temperature; this is due to the fact that the thermal boundary layer increases with the magnetic parameter.Fig. 7 represents graphs of temperature profiles     as function of  for various values of the heat source/sink parameter  .For fixed value of  , the temperature   

Fig. 8 .
Fig. 8. Effect of the magnetic parameter on the entropy generation number.

Fig. 10 .Fig. 11 .
Fig. 10.Effect of the Reynolds number on the entropy generation number.The influence of the Reynolds number Re l on the entropy generation number is plotted on fig.10.For a given value of  , the entropy generation number increases as the Reynolds number increases.The augmentation of the Reynolds number increases the contribution of the entropy generation number due to fluid friction and heat transfer in the boundary layer.

Fig. 12 .
Fig. 12.Effect of the Hartman number on the entropy generation number.The effect of the Hartman number Ha on the entropy generation number S N is plotted in fig.12.For a given  , as the Hartman number increases, the entropy generation number increases.The entropy generation number is proportional to the square of Hartman number which proportional to the magnetic field.The presence of the magnetic field creates additional entropy.
It is appropriate to define dimensionless number for entropy generation rate S N .This number is defined by dividing the local volumetric entropy generation rate