Open access peer-reviewed chapter

Entropy Generation in Viscoelastic Fluid Over a Stretching Surface

By Saouli Salah and Aïboud Soraya

Submitted: November 16th 2010Reviewed: April 15th 2011Published: October 10th 2011

DOI: 10.5772/21859

Downloaded: 2226

1. 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, 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, 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 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, 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.

2. 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 u,vin these directions.

Under the usual boundary layer approximations, the flow is governed by the following equations:

ux+vy=0E1
uux+vuy=ν2uy2k0(u3uxy2+v3uy3uy2uxy+uy2uy2)σB02ρuE2

The constant k0=α1ρis the viscoelastic parameter.

The boundary conditions are given by

y=0,u=uP=λx,v=0E3
y=,u=0,uy=0E4

The heat transfer governing boundary layer equation with temperature-dependent heat generation (absorption) is

ρCP(uTx+vTy)=k2Ty2+Q(TT)E5

The relevant boundary conditions are

y=0,T=TP=A(xl)2+TE6
y=,T=TE7

3. Analytical solution

The equation of continuity is satisfied if we choose a dimensionless stream function Ψ(x,y)such that

u=Ψy,v=ΨxE8

Introducing the similarity transformations

η=yλν,Ψ(x,y)=xνλf(η)E9

Momentum equation (2) becomes

f2(η)f(η)f(η)=f(η)k0λν[2f(η)f(η)f(η)fIV(η)f2(η)]Mnf(η)E10

where

Mn=σB02ρλE11

Now let us seek a solution of Eq. (8) in the form

f(η)=eαη,α0E12

which is satisfied by the following boundary conditions:

η=0,f(0)=0,f(0)=1E13
η=,f()=0,f()=0E14

On substituting (9) into (8) and using boundary conditions (12) and (13) the velocity components take the form

u=λxf(η)E15
v=νλf(η)E16

Where k1=k0λνis the viscoelastic parameter, and

α=1+Mn1k1E17

Defining the dimensionless temperature

Θ(η)=TTTPTE18

and using (11), (14), (15), Eq. (17) and the boundary conditions (6) and (7) can be written as

Θ(η)+Prα(1eαη)Θ(η)(2Preαηβ)Θ(η)=0E19
η=0,Θ(0)=1E20
η=,Θ()=0E21

Where and β=Qνλkare respectively the Prandtl number and the heat/sink parameter.

Introducing the variable

ξ=Prα2eαηE22

And inserting (21) in (18) we obtain

ξΘ(ξ)+Prα(1Prα2+ξ)Θ(ξ)(2βα2ξ)Θ(ξ)=0E23

And (19) and (20) transform to

ξ=Prα2,Θ(Prα2)=1E24
ξ=0,Θ(0)=0E25

The solution of Eq. (22) satisfying (23) and (24) is given by

Θ(ξ)=(α2Prξ)a+bM(a+b2,2b+1,ξ)M(a+b2,2b+1,Prα2)E26

The solution of (25) in terms of ηis written as

Θ(ξ)=eα(a+b)ηM(a+b2,2b+1,Prα2eαη)M(a+b2,2b+1,Prα2)E27

wherea=Pr2α2, b=Pr24α2β2α2and M(a+b2,2b+1,Prα2eαη)is the Kummer’s function.

4. Second law analysis

According to (Woods, 1975), the local volumetric rate of entropy generation in the presence of a magnetic field is given by

SG=kT2[(Tx)2+(Ty)2]+μT(uy)2+σB02Tu2E28

Eq. (27) clearly shows contributions of three sources of entropy generation. The first term on the right-hand side of Eq. (27) 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. It is appropriate to define dimensionless number for entropy generation rateNS. This number is defined by dividing the local volumetric entropy generation rate SGto a characteristic entropy generation rateSG0. For prescribed boundary condition, the characteristic entropy generation rate is

SG0=k(ΔT)2l2T2E29

therefore, the entropy generation number is

NS=SGSG0E30

using Eq. (11), (26) and (27), the entropy generation number is given by

NS=4X2Θ2(η)+RelΘ2(η)+RelBrΩf2(η)+BrHa2Ωf2(η)E31

where Reland Brare respectively the Reynolds number and the Brinkman number. ΩandHa, are respectively the dimensionless temperature difference and the Hartman number. These number are given by the following relationships

Rel=ullν,Br=μuP2kΔT,Ω=ΔTT,Ha=B0lσμE32

5. 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. Figs. 1 and 2 show the variations of the longitudinal velocity f(η)and the transverse velocity f(η)as function of ηfor several values of magnetic parameterMn. It can be observed that f(η)decreases with ηand f(η)increases with ηasymptotically for Mnkeeping constant. For a fixed positionη, both f(η)and f(η)decreases withMn, thus the presence of the magnetic field decreases the momentum boundary layer thickness and increase the power needed to stretch the sheet.

The effects of the viscoelastic parameter k1on the longitudinal velocity f(η)and the transverse velocity f(η)are illustrated on figs. 3 and 4. As it can be seen, for a fixed value ofη, both f(η)and f(η)decrease as viscoelastic parameter rises. This can be explained by the fact that, as the viscoelastic parameter increases, the hydrodynamic boundary layer adheres strongly to the surface, which in turn retards the flow in the longitudinal and the transverse directions.

Figure 1.

Effect of the magnetic parameter on the longitudinal velocity.

Figure 2.

Effect of the magnetic parameter on the transverse velocity.

Figure 3.

Effect of the viscoelastic parameter on the longitudinal velocity.

Figure 4.

Effect of the viscoelastic parameter on the transverse velocity.

Fig. 5 depicts the temperature profiles Θ(η)as function of ηfor different values of the Prandtl numberPr. 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.

Figure 5.

Effect of the Prandtl number on the temperature.

The temperature profiles Θ(η)as function of ηfor different values of the magnetic Mnare plotted in fig. 6. An increase in the magnetic parameter Mnresults 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 Θ(η)augments with the heat source/ sink parameterβ. This is due to the fact that the increase of the heat source/sink parameter means an increase of the heat generated inside the boundary layer leading to higher temperature profile.

The influence of the magnetic parameter Mnon the entropy generation number NSis shown on fig. 8. The entropy generation number NSdecreases with ηfor Mnkeeping constant. For

Figure 6.

Effect of the magnetic parameter on the temperature.

Figure 7.

Effect of the heat source/sink parameter on the temperature.

Figure 8.

Effect of the magnetic parameter on the entropy generation number.

fixed value ofη, the entropy generation number increases with the magnetic parameter, because the presence of the magnetic field creates more entropy in the fluid. Moreover, the stretching surface acts as a strong source of irreversibility.

Figure 9.

Effect of the Prandtl number on the entropy generation number.

Fig. 9 illustrates the effect of the Prandtl number Pron the entropy generation numberNS. 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.

Figure 10.

Effect of the Reynolds number on the entropy generation number.

The influence of the Reynolds number Relon 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.

Figure 11.

Effect of the dimensionless group on the entropy generation number.

The effect of the dimensionless group parameter BrΩ1on the entropy generation number NSis depicted in fig. 11. The dimensionless group determines the relative importance of 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.

Figure 12.

Effect of the Hartman number on the entropy generation number.

The effect of the Hartman number Haon the entropy generation number NSis 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.

6. 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 Prandtl 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:

  1. The velocities depend strongly on the magnetic and the viscoelastic parameters.

  2. The temperature varies significantly with the Prandlt number, the magnetic parameter and the heat source/sink parameter.

  3. 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.

  4. The surface acts as a strong source of irreversibility.

7. Nomenclature

Aconstant, KB0uniform magnetic field strength, Wb.m-2BrBrinkman number,Br=μu02kΔTCPspecific heat of the fluid, J.kg-1.K-1fdimensionless functionHaHartman numberHa=B0lσμkthermal conductivity of the fluid, W.m-1.K-1k1viscoelastic parameter, k1=k0λνk0viscoelastic parameter, m2lcharacteristic length, mMKummer’s functionMnmagnetic parameter, Mn=σB02ρλNSentropy generation number, NS=SGSG0PrPrandlt number, Pr=μCPkQrate of internal heat generation or absorption, W.m-3.K-1 RelReynolds number based on the characteristic length,Rel=ullνSGlocal volumetric rate of entropy generation, W.m-3.K-1 SG0characteristic volumetric rate of entropy generation, W.m-3.K-1Ttemperature, Kuaxial velocity, m.s-1ulplate velocity based on the characteristic length, m.s-1uPplate velocity, m.s-1vtransverse velocity, m.s-1xaxial distance, mXdimensionless axial distance,X=xlytransverse distance, mαpositive constantβheat source/sink parameter,β=Qνλkλproportional constant, s-1ηdimensionless variable, η=yλνξdimensionless variable, ξ==Prα2eαημdynamic viscosity of the fluid, kg.m-1.s-1νkinematic viscosity of the fluid, m2.s-1ΔTtemperature difference, ΔT=TpTΩdimensionless temperature difference,Ω=ΔTTΘdimensionless temperature, Θ=TTTpTρdensity of the fluid, kg.m-3σelectric conductivity, Ω-1.m-1subscriptsP platefar from the sheet

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Saouli Salah and Aïboud Soraya (October 10th 2011). Entropy Generation in Viscoelastic Fluid Over a Stretching Surface, Thermodynamics - Systems in Equilibrium and Non-Equilibrium, Juan Carlos Moreno-Piraján, IntechOpen, DOI: 10.5772/21859. Available from:

chapter statistics

2226total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

From Particle Mechanics to Pixel Dynamics: Utilizing Stochastic Resonance Principle for Biomedical Image Enhancement

By V.P.Subramanyam Rallabandi and Prasun Kumar Roy

Related Book

First chapter

Thermodynamics of Ligand-Protein Interactions: Implications for Molecular Design

By Agnieszka K. Bronowska

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us