A note on convective heat transfer of an MHD Jeffrey fluid over a stretching sheet

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. Sakiadis^1,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.³ Later on Vajravelu and Roper⁴ 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 cited^5–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. Chen¹⁰ examined the combined effects of viscous dissipation joule heating for the momentum and thermal transport of MHD Newtonian fluid. Liu¹¹ 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.¹² have examined the MHD heat and mass transfer in a rotating system with periodic suction. Anjli devi and Ganja¹³ 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.¹⁴ studied the radiative flow of Jeffrey fluid in a porous medium with power law heat flux and heat source. Vyas and Rai¹⁵ examined Radiative flow with variable thermal conductivity over a non-isothermal stretching sheet in a porous medium. Singh¹⁶ 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).²² 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 by²³ 

with S as the extra stress tensor defined by

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

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 B_o is imposed to the flow perpendicular to the sheet.

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

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, c_p the specific heat λ indicates the ratio of relaxation and retardation times and λ₁ the relaxation time, q_r 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 radiation²⁴ the radiative heat flux can be written as

where σ₁ 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 T⁴ can be expressed as a linear combination of the temperature. This can be accomplished by expanding T⁴ in the Taylor series about T_∞ and neglecting the higher order terms as

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

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

where A₁ and A₂ are the constants depending upon the properties of the fluid, l the characteristic length, T_w 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

The governing equation of motion (4) reduces to

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

Lawrence and Rao²⁵ 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

where

In PST case we define the following dimensionless temperature variable as

with ΔT = T_w − T_∞. The dimensionless energy equation can be written as

subject to the boundary conditions

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

By introducing the new variable

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

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

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

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

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

where

and M(a_∘,   b_∘,   ξ) is the Kummer’s function (Abramowitz and Stegun²²) and it is defined as

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

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

The local Nusselt number is defined as

In terms of dimensionless parameters it reduces to

The dimensionless wall temperature gradient can be written as

In PHF case we define the following dimensionless temperature variable as

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

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

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

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

The expression for the dimensionless wall temperature is obtained as

For PHF case the local Nusselt number can be expressed as

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, E_c 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 E_c 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 E_c 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 λ₁, 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.

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

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 E_c results a decrease in heat transfer rate for both cases. This is due to the fact that large value of Eckert number E_c 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.