This article focuses on the exact solution regarding convective heat transfer of a magnetohydrodynamic (MHD) Jeffrey fluid over a stretching sheet. The effects of joule and viscous dissipation, internal heat source/sink and thermal radiation on the heat transfer characteristics are taken in account in the presence of a transverse magnetic field for two types of boundary heating process namely prescribed power law surface temperature (PST) and prescribed heat flux (PHF). Similarity transformations are used to reduce the governing non-linear momentum and thermal boundary layer equations into a set of ordinary differential equations. The exact solutions of the reduced ordinary differential equations are developed in the form of confluent hypergeometric function. The influence of the pertinent parameters on the temperature profile is examined. In addition the results for the wall temperature gradient are also discussed in detail.

It is well known that investigation regarding boundary layer flows over continuous stretching surfaces finds numerous important applications occurring in many engineering process. Drawing of plastic films and wires, manufacture of foods, crystal growing, liquid films in condensation process etc. are a few examples of practical importance. Further it has several practical applications in the field of metallurgy and chemical engineering such as extrusion process, heat materials traveling between a feed roll and a wind up roll. Sakiadis1,2 performed the pioneering works by studying the flow induced by a moving surface. The closed form exponential solution for the flow by a linear stretching is examined by Crane.3 Later on Vajravelu and Roper4 analyzed the flow and heat transfer over a stretching sheet with prescribed surface temperature. Many researchers investigated the stretching sheet problems by considering different effects (some studies can be cited5–9).

A considerable interest has been revealed by different researchers in studying the flow and heat transfer of an electrically conducting and heat generating fluid due to its important applications in the process of purification of molten metals from non-metallic inclusions. In polymer technology and specific metallurgical operations the rate of cooling of continuous stretched strips is controlled by applying the principles of MHD techniques. This application is well marked in case of drawing, annealing, and thinning of copper wires. Chen10 examined the combined effects of viscous dissipation joule heating for the momentum and thermal transport of MHD Newtonian fluid. Liu11 studied the flow and heat transfer in an electrically conducting second grade fluid with transverse magnetic field with power law surface temperature. Parida et al.12 have examined the MHD heat and mass transfer in a rotating system with periodic suction. Anjli devi and Ganja13 have studied the MHD flow over a stretching porous surface with viscous dissipation.

At high operating temperatures in engineering process radiation may play a dynamic role under many non-isothermal situations. In polymer processing industry if the entire system involving the polymer extrusion process is placed in a thermally controlled environment then thermal radiation might be great important in controlling heat transfer process. Hayat et al.14 studied the radiative flow of Jeffrey fluid in a porous medium with power law heat flux and heat source. Vyas and Rai15 examined Radiative flow with variable thermal conductivity over a non-isothermal stretching sheet in a porous medium. Singh16 studied the magneto convection flow of viscoelastic fluid over stretching sheet with heat source and radiation effects. Many other researchers have studied the flow and heat transfer of fluids by taking the radiation heat source/sink effects.17–21 

Motivated by above studies we intend to investigate the flow and heat transfer of an MHD Jeffrey fluid over a stretching sheet subject to power law temperature in the heat source/ sink with radiation. Also we are considering two general cases of non-isothermal boundary conditions, namely (1) the sheet with prescribed power law surface temperature (PST case) (2) and the sheet with prescribed power law heat flux (PHF case). Using appropriate similarity transformations the highly non-linear partial differential equations are reduced to a set of non-linear ordinary differential equations. We derived closed form analytical solution for non-dimensional velocity and temperature distribution in the form of confluent hypergeometric function (Kummer’s function).22 The effects of non-dimensional parameters, magnetic number, Eckert number, radiation parameter, heat source/sink parameter and Deborah number on the temperature profile are explored graphically and discussed in details.

The constitutive equations for a Jeffrey fluid are given by23 

τ = p I + S ,
(1)

with S as the extra stress tensor defined by

S = μ 1 + λ R 1 + λ 1 R 1 t + V . R 1 ,
(2)

where τ is the Cauchy stress tensor, p the pressure, μ the dynamic viscosity, λ and λ1 are the material parameters of the Jeffrey fluid and R1 the Rivlin-Ericksen tensor defined by

R 1 = V + V t .
(3)

We consider the steady two dimensional boundary layer flow of an electrically conducting Jeffrey fluid past a stretching sheet in the presence of heat source and chemical reaction. The flow is generated due to linear stretching of the sheet caused by simultaneous application of equal and opposite forces along the x -axis while keeping the origin fixed (Fig. 1). A uniform magnetic field of strength Bo is imposed to the flow perpendicular to the sheet.

FIG. 1.

Schematic diagram of flow towards a stretching sheet.

FIG. 1.

Schematic diagram of flow towards a stretching sheet.

Close modal

The boundary layer equations which govern the flow and heat transfer of a Jeffrey fluid over a stretching surface are

u x + v y = 0 ,
(4)
u u x + v u y = υ 1 + λ 2 u y 2 + λ 1 u 3 u x y 2 + v 3 u y 3 u x 2 u y 2 + u y 2 u x y σ B o 2 ρ u ,
(5)
ρ c p u T x + v T y = k 2 T y 2 + μ 1 + λ u y 2 + λ 1 u y y u u x + v u y + q T T + σ B o 2 u 2 q r y ,
(6)

where u and v are the velocity components in the x and y direction respectively, T the fluid temperature, υ the kinematic viscosity, ρ the fluid density, cp the specific heat λ indicates the ratio of relaxation and retardation times and λ1 the relaxation time, qr the radiative heat flux and q is the rate of volumetric heat generation/absorption. The last three terms in Eq. (6) are the work done due to deformation, internal heat generation/absorption and the Joule heating. By using Rosseland approximation for thermal radiation24 the radiative heat flux can be written as

q r = 4 σ 1 3 α R T 4 y ,
(7)

where σ1 and αr are the Stephan-Boltzmann constant and the mean absorption coefficient. We assume that the difference in the temperature within the flow is such that T4 can be expressed as a linear combination of the temperature. This can be accomplished by expanding T4 in the Taylor series about T and neglecting the higher order terms as

T 4 4 T 3 T 3 T 4 .
(8)

Using Eq. (7) in Eq. (6), the energy equation becomes

u T x + v T y = k ρ c p 2 T y 2 + μ ρ c p ( 1 + λ ) u y 2 + λ 1 u y y u u x + v u y + q ρ c p T T + σ B o 2 ρ c p u 2 + 16 σ 1 T 3 3 ρ c p α R 2 T y 2 .
(9)

In this paper we intend to investigate the thermal transport phenomenon for two non-isothermal conditions, namely prescribed surface temperature (PST) and prescribed surface heat flux (PHF). Correspondingly, we consider the following boundary conditions for the flow and temperature

u = u w = c x , v = 0 , at  y = 0 ,
(10)
u = 0 , u y 0 as  y
(11)
T = T w ( x ) = T + A 1 x l 2 PST case q w ( x ) = k T y = A 2 x l 2 PHF case at  y = 0
(12)
T T as  y ,
(13)

where A1 and A2 are the constants depending upon the properties of the fluid, l the characteristic length, Tw and T are the temperatures at the wall and far away from the wall respectively.

To facilitate the analysis we introduce the following similarity transformation and dimensionless variable

η = c υ y , u = cxf ( η ) , v = c υ f ( η ) .
(14)

The governing equation of motion (4) reduces to

f + 1 + λ ff f 2 + β f 2 ff i v Mf = 0 ,
(15)

where prime denotes differentiation with respect to η. Furthermore β = λ1c is the Deborah number and M = σ B o 2 ρ c the magnetic parameter. The corresponding boundary conditions are then

f ( η ) = 1 , f ( η ) = 0  at  η = 0 ,
(16)
f ( η ) 0  as  η .
(17)

Lawrence and Rao25 have given a general method and obtained all non-unique solution of the reduced Eq. (14), among all these solutions of the form, we consider the realistic solution of the form

f ( η ) = 1 exp α η α ,
(18)

where

α = 1 + λ + M 1 + β .
(19)

In PST case we define the following dimensionless temperature variable as

T ( η ) = T + Δ T θ ( η ) ,
(20)

with ΔT = TwT. The dimensionless energy equation can be written as

1 + N θ + Pr f θ + Pr Q 2 f θ = Pr E c M f 2 + 1 1 + λ f 2 + β f f f f f i v ,
(21)

subject to the boundary conditions

θ ( 0 ) = 1 , θ ( ) = 0 ,
(22)

here the dimensionless parameters Pr, N, Q and Ec are the Prandtl number, radiation parameter, internal heat generation/absorption parameter and Eckert number, respectively and are defined as

Pr = ρ c p ν k , N = 16 σ 1 T 3 3 α R k , Q = q c ρ c p , E c = c 2 l 2 A 1 c p .
(23)

By introducing the new variable

ξ = Pr * e α η ,
(24)

where Pr * = Pr α 2 1 + N is the modified Prandlt number. Making use of Eqs. (18) and (24) in Eqs. (21) and (22) the transform problem is

ξ d 2 θ d ξ 2 + 1 Pr * ξ d θ d ξ + 2 + Pr * Q ξ θ = E c α 2 Pr * M α 2 + 1 1 + λ 1 + β ,
(25)
θ ξ = Pr * = 1 , θ ( ξ = 0 ) = 0 .
(26)

Eqs. (25) and (26) constitute a non- homogeneous boundary value problem. Denoting the solution of homogenous part of Eq. (25) by θc and non-homogenous part by θp, so we can write

θ ( ξ ) = θ c ( ξ ) + θ p ( ξ ) .
(27)

A closed form particular solution of Eq. (24) can be written as

θ p ( ξ ) = E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q ξ Pr * 2 ,
(28)

Making use of the boundary condition (26) the homogenous part of the temperature can be written in the form of confluent hypergeometric function (Kummer’s function) as

θ c ( ξ ) = C 1 M χ 2 , Z + 1 , ξ .
(29)

Using the boundary conditions (26) the solution of the Eq. (25) is determined to be

θ ( ξ ) = 1 + E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q ξ Pr * χ M χ 2 , Z + 1 , ξ M χ 2 , Z + 1 , Pr * E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q ξ Pr * 2 ,
(30)

where

χ = Pr * + Pr * 2 4 Pr * Q 2 , Z = Pr * 2 4 Pr * Q ,
(31)

and M(a,   b,   ξ) is the Kummer’s function (Abramowitz and Stegun22) and it is defined as

M ( a , b , ξ ) = 1 + n = 1 a n ξ n b n n ! ,
(32)
( a ) n = a ( a + 1 ) ( a + 2 ) ( a + n 1 ) ,
( b ) n = b ( b + 1 ) ( b + 2 ) ( b + n 1 ) .

The temperature profile in term of η can be written by invoking Eq. (24) in Eq. (30),

θ ( η ) = 1 + E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q e χ α η M χ 2 , Z + 1 , Pr * e α η M χ 2 , Z + 1 , Pr * E Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q e 2 α η .
(33)

Note that this solution also satisfies the boundary conditions given in Eq. (22)

The local Nusselt number is defined as

N u x = h x k = q w T w T x k .
(34)

In terms of dimensionless parameters it reduces to

N u x Re x 1 2 = θ ( 0 ) .
(35)

The dimensionless wall temperature gradient can be written as

θ ( 0 ) = α 1 + E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q Pr * χ 2 Z + 1 M χ 1 , Z + 2 , Pr * M χ 2 , Z + 1 , Pr * χ + 2 α E c Pr 1 + M α 2 1 + N 4 2 Pr * + Pr * Q .
(36)

In PHF case we define the following dimensionless temperature variable as

Φ = T T q w x / k Re x 1 2 .
(37)

The transformed energy equation along with the boundary conditions can be written as

1 + N Φ + Pr f Φ + Pr Q 2 f Φ = Pr E c M f 2 + 1 1 + λ f 2 + β f f f f f ,
(38)
Φ ( 0 ) = 1 , Φ ( ) = 0 ,
(39)

where the Eckert number in PHF case is defined as E c = k c 2 l 2 c / υ / A 2 c p while the other parameters are same as defined in PST case. Making use of Eqs. (18) and (24) in Eqs. (38) and (39) we obtained the following

ξ d 2 Φ d ξ 2 + 1 Pr * ξ d Φ d ξ + 2 + Pr * Q ξ Φ = E α 2 Pr * M α 2 + 1 1 + λ 1 + β ,
(40)
Φ ξ = Pr * = 1 α Pr * , Φ ( ξ = 0 ) = 0 .
(41)

The analytical solution of Eq. (40) subject to the boundary condition in Eq. (41) by following the same procedure as describe in PST case, is given by

Φ ( ξ ) = 2 E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q + 1 / α ξ Pr * χ M χ 2 , Z + 1 , ξ χ M χ 2 , Z + 1 , Pr * Pr * χ 2 Z + 1 M χ 1 , Z + 2 , Pr * E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q ξ Pr * 2 .
(42)

The temperature profile in term of η can be written by substituting Eq. (24) in Eq. (42)

Φ ( η ) = 2 E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q + 1 / α e χ α η M χ 2 , Z + 1 , Pr * e α η χ M χ 2 , Z + 1 , Pr * Pr * χ 2 Z + 1 M χ 1 , Z + 2 , Pr * E c Pr * M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q e 2 α η .
(43)

The expression for the dimensionless wall temperature is obtained as

Φ ( 0 ) = 2 E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q + 1 / α M χ 2 , Z + 1 , Pr * χ M χ 2 , Z + 1 , Pr * Pr * χ 2 Z + 1 M χ 1 , Z + 2 , Pr * E c Pr M α 2 + 1 1 + λ 1 + β 1 + N 4 2 Pr * + Pr * Q .
(44)

For PHF case the local Nusselt number can be expressed as

N u x Re x 1 2 = 1 / Φ ( 0 ) .
(45)

In this article an MHD boundary layer problem for momentum and heat transfer in Jeffrey fluid flow over a non-isothermal stretching sheet in the presence of dissipative energy, thermal radiation and internal heat source/ sink is investigated. The highly non-linear governing partial differential equations are reduced into a set of non-linear ordinary differential equations by applying suitable similarity transformations and their analytic solutions are obtained in the form of confluent hypergeometric function.

This section highlights the influence of pertinent parameters on the temperature profiles. Figures 2(a) and 2(b) describe the effect of ratio of relaxation and retardation times parameter λ on temperature profiles θ(η) and Φ(η) in both PST and PHF cases, respectively for the selected values of Pr, M, Q, N, Ec and β. It is quite clear that an increase in the value of λ results in an increase in the temperature profiles θ(η) and Φ(η) for both PST and PHF cases. Thus broadened the thermal boundary layer thickness such that the wall temperature gradient decreases in PST case and the surface temperature increases in PHF case. In order to see the effects of Eckert number Ec figures 3(a) and 3(b) are displayed for both PST and PHF cases on the temperature distribution θ(η) and Φ(η). It is evident from these figures that the thermal boundary is enhanced by increasing the values of Eckert number Ec so the dissipative energy becomes more important with an enlargement of temperature profile. Also the energy dissipation represents an increase in temperature gradient in PST case and the wall temperature in PHF case. Figures 4(a) and 4(b) are plotted to see the effects of radiation parameter N on the temperature profiles θ(η) and Φ(η) for both PST and PHF cases. It is obvious that by increasing the radiation parameter N the temperature profiles increases in both cases. Thus the thermal boundary layer thickness enhanced in the presence of thermal radiation such that the wall temperature decreases in PST case and the surface temperature increases in PHF case. This result indicates that thermal radiation is to reduce the heat transfer rate. The influence of heat generation/absorption parameter Q on the temperature profiles θ(η) and Φ(η) is displayed in figures 5(a) and 5(b) for both PST and PHF cases respectively. It is worth noting that for the case of Q > 0 the energy may generates in the boundary layer, in consequences, the temperature profiles θ(η) and Φ(η) of the fluid temperature will increase. In the other hand, for Q < 0 (heat sink) energy absorbs in the boundary layer which causes the temperature profile θ(η) and Φ(η) to decreases in both PST and PHF cases. The effects of the magnetic number M on temperature profile for both PST and PHF cases are plotted in figures 6(a) and 6(b). An increase in the magnetic number M enhances the boundary layer thickness for both cases. Which occurs due to the fact that a resistive Lorentz force is produced by applying the transverse magnetic field to the electrically conducting fluid. This force has a tendency to slow down the motion of the fluid in the boundary layer and increases the temperature distribution. Figures 7(a) and 7(b) demonstrate the effects of Deborah number β on the temperature distribution θ(η) and Φ(η) respectively. For both PST and PHF cases, keeping different values of other parameters fixed. As we increase the value of Deborah number β the temperature distribution θ(η) and Φ(η) decreases in both PST and PHF cases. Since the Deborah number β depend on the retardation time λ1, physically a large retardation time of any material makes it less viscous, which may results an increase in its motion, which consequently reduces the thermal boundary layer thickness and temperature distribution.

FIG. 2.

Influence of λ on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 2.

Influence of λ on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal
FIG. 3.

Influence of the Eckert number Ec on the temperatures profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 3.

Influence of the Eckert number Ec on the temperatures profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal
FIG. 4.

Influence of N on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 4.

Influence of N on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal
FIG. 5.

Influence of Q on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 5.

Influence of Q on the temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal
FIG. 6.

Influence of M on temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 6.

Influence of M on temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal
FIG. 7.

Influence of β on temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

FIG. 7.

Influence of β on temperature profiles: (a) θ(η) for PST case and (b) Φ(η) for PHF case.

Close modal

Table I displays the effect of various parameters on dimensionless temperature at the wall. From this table we can observe that an increasing values of the magnetic parameter M, heat generation/absorption parameter Q, Eckert number Ec, relaxation times parameter λ, and radiation number N is to increase the wall temperature gradient θ′(0) in PST case and the wall temperature Φ(0) in PHF case. Also we observe that the effect of Deborah number β is to increase the wall temperature gradient θ′(0) in PST case and the wall temperature Φ(0) in PHF case.

TABLE I.

Wall temperature gradient θ′(0) for PST case and wall temperature Φ(0) for PHF case.

M Q E β λ N Pr θ′(0)PST case Φ(0)PHF case
0.3  0.1  0.1  0.1  0.2  0.2  1.04481  0.95875 
0.5              1.00949  0.99104 
0.7              0.97464  1.02458 
0.3  -0.1            1.15458  0.87057 
  0.0            1.10430  0.90889 
  0.1            1.04481  0.95875 
0.3  0.1  0.0          1.08648  0.92040 
    0.2          1.00314  0.99711 
    0.4          0.91978  1.07383 
0.3  0.1  0.1  0.0        1.02542  0.97616 
      0.3        1.07510  0.93288 
      0.6        1.10750  0.90686 
0.3  0.1  0.1  0.1  0.0      1.06961  0.93754 
        0.4      1.01989  0.98121 
        0.8      0.96770  1.03219 
0.3  0.1  0.1  0.1  0.2  0.0    1.18053  0.85308 
          0.3    0.98829  1.01140 
          0.6    0.84836  1.17204 
0.3  0.1  0.1  0.1  0.2  0.2  0.7  0.80269  1.23660 
            1.0  1.04481  0.95875 
            1.2  1.18053  0.853084 
M Q E β λ N Pr θ′(0)PST case Φ(0)PHF case
0.3  0.1  0.1  0.1  0.2  0.2  1.04481  0.95875 
0.5              1.00949  0.99104 
0.7              0.97464  1.02458 
0.3  -0.1            1.15458  0.87057 
  0.0            1.10430  0.90889 
  0.1            1.04481  0.95875 
0.3  0.1  0.0          1.08648  0.92040 
    0.2          1.00314  0.99711 
    0.4          0.91978  1.07383 
0.3  0.1  0.1  0.0        1.02542  0.97616 
      0.3        1.07510  0.93288 
      0.6        1.10750  0.90686 
0.3  0.1  0.1  0.1  0.0      1.06961  0.93754 
        0.4      1.01989  0.98121 
        0.8      0.96770  1.03219 
0.3  0.1  0.1  0.1  0.2  0.0    1.18053  0.85308 
          0.3    0.98829  1.01140 
          0.6    0.84836  1.17204 
0.3  0.1  0.1  0.1  0.2  0.2  0.7  0.80269  1.23660 
            1.0  1.04481  0.95875 
            1.2  1.18053  0.853084 

The influence of Nusselt number in terms of −θ′(0) in PST case and 1/Φ(0) in PHF case with different values observed through Table I. A remarkable decrease is noticed in the value of Nusselt number with an increase in the value of magnetic parameter M for both PST and PHF cases. The effect of heat source/ sink parameter Q is noticed. The heat transfer rate increases as Q < 0 and decreases as Q > 0. An increasing value of Eckert number Ec results a decrease in heat transfer rate for both cases. This is due to the fact that large value of Eckert number Ec indicates more energy dissipating in the boundary layer, consequently the fluid temperature increases and the heat transfer rate decreases. A decreasing effect in the value of Nusselt number is observed by increasing the value of radiation parameter N. Furthermore, by increasing the value of Prandtl number Pr the Nusselt number increases. As the fluid with larg-e Prandtl number Pr has higher heat capacity and hence increases the heat transfer rate. Moreover, a decrease in the heat transfer rate is seen from the Table I with an increase in the values of ratio of relaxation and retardation times parameter λ, while a quite opposite trend is observed in case of Deborah number β.

An analytic investigation has been carried out on an MHD Jeffrey fluid past over a stretching sheet with viscous dissipation as well as energy dissipation due to joule heating, internal heat generation/absorption and thermal radiation. By invoking suitable transformation highly nonlinear partial differential equations were reduced into non-linear ordinary differential equation. The closed form solution of the reduced heat equation is obtained in the form of Kummer’s function for two types of thermal boundary condition that is prescribed surface temperature (PST) and prescribed heat flux (PHF). The effects of various physical parameters such as magnetic field, thermal radiation, heat source/sink, Eckert number, relaxation time parameter on thermal behavior are examined and discussed graphically in details.

1.
B. C.
Sakiadis
, “
Boundary layer behavior on continuous solid surfaces: I. Boundary layer equations for two-dimensional and axisymmetric flow
,”
American Institute of Chemical Engineering Journal
7
,
26
-
28
(
1961
).
2.
B. C.
Sakiadis
, “
Boundary layer behavior on continuous solid surfaces: II. Boundary layer on a continuous flat surface
,”
American Institute of Chemical Engineering Journal
7
,
221
-
225
(
1961
).
3.
L. J.
Crane
, “
Flow past a stretching plate
,”
Zeitschrift für angewandte Mathematik und Physik
21
,
645
-
655
(
1970
).
4.
K.
Vajravelu
and
T.
Roper
, “
Flow and heat transfer in a second grade fluid over a stretching sheet
,”
International Journal of Non-linear Mechanics
34
,
1031
-
1036
(
1999
).
5.
C. Y.
Wang
, “
The three-dimensional flow due to a stretching flat surface
,”
Physics of Fluids
27
,
1915
-
1917
(
1984
).
6.
A.
Shahzad
,
R.
Ali
, and
M.
Khan
, “
On the exact solution for axissymetric flow and heat transfer over a nonlinear radially stretching sheet
,”
Chinese Physics Letters
29
,
084705
(
2012
).
7.
A.
Shahzad
and
R.
Ali
, “
MHD flow of a non-Newtonian Power law fluid over a vertical stretching sheet with the convective boundary condition
,”
Walailak Journal of Science and Technology (WJST)
10
(
1
),
43
-
56
(
2012
).
8.
T.
Fang
,
J.
Zhang
, and
S.
Yao
, “
Slip MHD viscous flow over a stretching sheet-An exact solution
,”
Communication in Nonlinear Science and Numerical Simulation
14
,
3731
-
3737
(
2009
).
9.
R.
Cortell
, “
Effects of viscous dissipation and work done by deformation on the MHD flow and heat transfer of a viscoelastic fluid over a stretching sheet
,”
Physics Letters A
357
,
298
-
305
(
2006
).
10.
C.-H.
Chen
, “
Combined effects of Joule heating and viscous dissipation on magnetohydrodynamic flow past a permeable, stretching surface with free convection and radiative heat transfer
,”
ASME Journal of Heat Transfer
132
,
article no: 064503
(
2010
).
11.
I.-C.
Liu
, “
Flow and heat transfer of an electrically conducting fluid of second grade over a stretching sheet subject to a transverse magnetic field
,”
International Journal of Heat and Mass Transfer
47
,
4427
-
4437
(
2004
).
12.
S. K.
Parida
,
M.
Acharya
,
G. C.
Dash
, and
S.
Panda
, “
MHD heat and mass transfer in a rotating system with periodic suction
,”
Arabian Journal of Science and Engineering
36
,
1139
-
1151
(
2011
).
13.
S. P.
Anjali
and
B.
Ganga
, “
Viscous dissipation effects on nonlinear MHD flow in a porous medium over a stretching porous surface
,”
International Journal of Applied Mathematics and Mechanics
5
,
45
-
49
(
2009
).
14.
T.
Hayat
,
S. A.
Shehzad
,
M.
Qasim
, and
S.
Obaidat
, “
Radiative flow of Jeffrey fluid in a porous medium with power law heat flux and heat source
,”
Nuclear Engineering and Design
243
,
15
-
19
(
2012
).
15.
P.
Vyas
and
A.
Rai
, “
Radiative flow with variable thermal conductivity over a non-isothermal stretching sheet in a porous medium
,”
International Journal of Contemporary Mathematical Sciences
5
,
2685
-
2698
(
2010
).
16.
A. K.
Singh
, “
Heat source and radiation effect on magneto-convection flow of a viscoelastic fluid past a stretching sheet: Analysis with Kummer’s function
,”
International Communication in Heat and Mass Transfer
35
,
637
-
642
(
2008
).
17.
S. K.
Khan
, “
Heat transfer in a viscoelastic fluid flow over a stretching surface with heat source/sink, suction/blowing and radiation
,”
International Journal of Heat and Mass Transfer
49
,
628
-
639
(
2006
).
18.
R. C.
Bataller
, “
Viscoelastic fluid flow and heat transfer over a stretching sheet under the effects of a non-uniform heat source, viscous dissipation and thermal radiation
,”
International Journal of Heat and Mass Transfer
50
,
3152
-
3162
(
2007
).
19.
S.
Nadeem
,
S.
Zaheer
, and
T.
Fang
, “
Effects of thermal radiation on the boundary layer flow of a Jeffrey fluid over an exponentially stretching surface
,”
Numerical. Algorithms Journal
57
,
187
-
205
(
2011
).
20.
M. M.
Nandeppanavar
,
K.
Vajravelu
, and
M. S.
Abel
, “
Heat transfer in MHD viscoelastic boundary layer flow over a stretching sheet with thermal radiation and non-uniform heat source/sink
,”
Communication in NonlinearScienceand NumericalSimulation
16
,
3578
-
3590
(
2011
).
21.
R.
Cortell
, “
Combined effect of viscous dissipation and thermal radiation on fluid flows over a non-linearly stretched permeable wall
,”
Meccanica
47
,
769
-
781
(
2012
).
22.
M.
Ambromowitz
and
F.
Stegun
,
Handbook of Mathematical Functions
(
Dover
,
Newyork
,
1965
), pp.
556
-
563
.
23.
S.
Nadeem
and
N. S.
Akbar
, “
Peristaltic flow of a Jeffrey fluid with variable viscosity in an asymmetric channel
,”
Zeitschrift für Naturforschung A
64
,
713
-
722
(
2009
).
24.
M. A.
Seddeek
and
M. S.
Abdelmeguid
, “
Effects of radiation and thermal diffusivity on the heat transfer over a stretching with variable heat flux
,”
Physics Letters A
348
,
172
-
179
(
2006
).
25.
P. S.
Lawrence
and
B. N.
Rao
, “
The non-uniqueness of MHD flow of a visco-elastic fluid past a stretching sheet
,”
Acta Mechanica
112
,
223
-
228
(
1995
).