Thermal Radiation Effects on MHD Boundary Layer Flow over an Exponentially Stretching Surface ()
1. Introduction
The study of boundary layer flow and its applications are vital for advancement in the field of technology and engineering. The computation and computer coordinated applications of flow over a stretching surface are playing a pivotal role in different realm of industrial products of aerodynamics, polymers and metallurgy, such as liquid films in condensation process, artificial fibers, glass fiber, metal spinning, the cooling process of metallic plate in a cooling bath and glass, wire drawing, paper production, aerodynamic extrusion of plastic sheets, crystal growing, cable coating and many others, to get end product of desired quality and parameters. Sakiadis [1] probably was the first who investigated boundary layer flow on a moving continuous solid surface. Crane [2] extended this concept to a linearly stretching plate whose velocity is linearly proportional to the distance from the slit and produced an exact analytical solution for the steady two-dimensional flow problems. Gupta and Gupta [3] , Carragher and Crane [4] , Grubka and Bobba [5] , Chen and Char [6] , Ali [7] , Andersson [8] , Ariel et al. [9] , Ishak et al. [10] , Jat and Chaudhary [11] [12] , Wang [13] and Nadeem et al. [14] analyzed the effects of heat transfer on a stretching surface taking into account different aspects of the problem.
Boundary layer flow and heat transfer over an exponentially stretching surface have wider applications in technology such as in case of annealing and thinning of copper wires. Magyari and Keller [15] obtained analytical and numerical solutions for boundary layer flow over an exponentially stretching continuous surface with an exponential temperature distribution. Many other problems on exponentially stretching surface under different physical situations were observed by Elbashbeshy [16] , Partha et al. [17] , Khan [18] , Sanjayanand and Khan [19] and El-Aziz [20] .
At higher operating temperature, the effects of thermal radiation and heat transfer play a pivotal role on the fluid flow problem of boundary layer. The application of controlled heat transfer in polymer industries is very important to get final product of desired parameters. The modern system of electric power generation, plasma, space vehicles, astrophysical flows and cooling of nuclear reactors are governed by applications of thermal radiation and heat transfer of fluid flow. Elbashbeshy [21] determined the effect of radiation on flow of an incompressible fluid along a heated horizontal stretching sheet. Sajid and Hayat [22] extended this concept by investigating the influence of thermal radiation on the boundary layer flow over an exponentially stretching sheet and solved the problem analytically. Recently, Bidin and Nazar [23] , Jat and Chaudhary [24] , Nadeem et al. [25] and Mukhopadhyay and Gorla [26] investigated various aspects of such problem either analytically or numerically.
With reference to above significant studies and in view of importance of MHD applications in various field of technologies, the objective of present paper is to investigate the effect of thermal radiation on an electrically conducting two-dimensional boundary layer incompressible viscous fluid flow over an exponentially stretching surface in the presence of uniform magnetic field by using Rosseland approximation. Numerical results of the momentum and energy equations are computed by using shooting method. The promising results of velocity and temperature distributions, local skin-friction coefficient and surface heat transfer are discussed for various physical parameters and simplified their effects for different conditions.
2. Problem Formulation
Consider the steady two-dimensional laminar boundary layer flow
of a viscous incompressible electrically conducting radiative fluid over continuous exponentially stretching surface in the presence of an externally applied normal magnetic field of constant strength
. The x-axis is taken along the stretching surface in the direction of motion and y-axis is taken perpendicular to it. The stretching surface has a uniform temperature
and a linear velocity
while temperature of flow external to the boundary layer is
. The system of governing boundary layer equations (which model Figure 1) are given by:
![]()
Figure 1. Sketch of the physical problem.
(1)
(2)
(3)
where T0 is the reference temperature, L is the reference length, U0 is the reference velocity,
is the coef-
ficient of kinematic viscosity,
is the coefficient of viscosity,
is the fluid density,
is the electrical conductivity, Cp is the specific heat at constant pressure, T is the temperature,
is the thermal conductivity and qr is the radiative heat flux. The other symbols have their usual meanings.
The boundary conditions are:
(4)
By using Rosseland approximation of the radiation for an optically thick boundary layer, the radiative heat flux qr is expressed (Bidin and Nazar [23] ) as:
(5)
where
is the Stefan-Boltzmann constant and
is the mean absorption coefficient. The above radiative heat flux qr is effective at a point away from boundary layer surface in an intensive absorption flow. Considering that the temperature variation within the flow is very small, the T4 may be expressed as a linear function of temperature T. Expanding T4 by Taylor’s series about temperature
and neglecting higher-order terms, hence
(6)
Using Equation (5) and (6), equation (3) is reduced to:
(7)
3. Similarity Analysis
The continuity Equation (1) is identically satisfied if we defined stream function
as:
(8)
For the solution of momentum and energy Equations (2) and (7), introducing the following dimensionless variables:
(9)
(10)
(11)
Using Equations (8) to (11), Equations (2) and (7) are reduced to:
(12)
(13)
The boundary conditions are:
(14)
where prime (') denote differentiation with respect to η,
is the Magnetic parameter,
is the Radiation parameter,
is the Prandtl number and
is the Eckert number.
4. Numerical Solution of the Problem
For numerical solution of the Equations (12) and (13), we use the following power series in terms of small magnetic parameter M as:
(15)
(16)
Substituting the values of
and
from Equations (15) and (16) and its derivatives in Equations (12) and (13), and then equating the coefficients of like powers of M, we get the following set of equations:
(17)
(18)
(19)
(20)
(21)
(22)
The corresponding boundary conditions are:
(23)
The Equation (17) is same as that obtained by Bidin and Nazar [23] for non-magnetic case and the remaining equations from (18) to (22) are ordinary linear differential equations and have been solved numerically by Shooting method with boundary condition (23). The velocity and temperature distributions for various values of parameters are shown in Figures 2-6 respectively.
5. Local Skin Friction Coefficient and Local Nusselt Number
The important physical quantities are the local skin-friction coefficient Cf and the local Nusselt number Nu, which are defined as:
![]()
Figure 2. Velocity distribution against η for various values of M.
![]()
Figure 3. Temperature distribution against η for various values of M with K = 0.5, Pr = 1 and Ec = 0.0.
![]()
Figure 4. Temperature distribution against η for various values of K with M = 0.04, Pr = 1 and Ec = 0.0.
![]()
Figure 5. Temperature distribution against η for various values of Pr with M = 0.04, K = 0.5 and Ec = 0.0.
![]()
Figure 6. Temperature distribution against η for various values of Ec with M = 0.04, K = 0.5 and Pr = 1.
(24)
and
(25)
In the present case which can be expressed in dimensionless form as:
(26)
and
(27)
where
is the surface shear stress and
is the local Reynolds number. The numerical
values of
and
are proportional to the local skin-friction coefficient Cf and local Nusselt number Nu at the surface respectively and these are presented by Table 1 for various values of the physical parameters.
6. Results and Discussion
Figure 2 shows variation of velocity distribution
against
for various values of the magnetic parameter M. This figure shows that the fluid velocity decreases with increasing value of the magnetic parameter M, due to the effect of Lorentz force produced by transverse magnetic field causes deceleration of fluid velocity.
Figures 3-6 show the temperature distributions
against
for various values of the magnetic parameter M, the radiation parameter K, the Prandtl number Pr and the Eckert number Ec. It is observed from these figures that the temperature distribution
increases with increasing value of any parameter, such as the magnetic parameter M, the radiation parameter K and the Eckert number Ec. However, it decreases with increasing value of the Prandtl number Pr. An increasing Prandtl number Pr, causes decrease in thermal boundary layer of fluid flow.
The values of the local skin-friction coefficient Cf and the local Nusselt number Nu in terms of
and
respectively, are presented in the Table 1, for various values of the magnetic parameter M, the radiation parameter K and the Prandtl number Pr, with the Eckert number Ec = 0.0. It is significant that the local skin-friction coefficient Cf and the local Nusselt number Nu decreases with increasing value of the magnetic parameter M. Moreover, the local Nusselt number Nu decreases with increasing value of the radiation parameter K, whereas the reverse phenomena occurs for the Prandtl number Pr. Further, Table 1 shows that all values of
and
are negative, corresponding to various values of physical parameters. A negative sign of
implies the exertion of drag force on the surface and a negative sign of
implies heat transfer from the surface.
![]()
Table 1. Variation of surface shear stress f'' (0) with M and surface heat transfer rate θ' (0) with M, K, Pr and Ec = 0.0.
7. Conclusion
The characteristic relationships among various parameters influencing viscous incompressible electrically conducting fluid over an exponentially stretching surface in the presence of a uniform magnetic field with thermal radiation have been analyzed and illustrated graphically. The similarity equations are determined and solved numerically by shooting method. It is observed that thickness of the velocity boundary layer, the local skin-friction coefficient and the local Nusselt number decreases with increasing value of the magnetic parameter. However, thickness of the thermal boundary layer increases with increasing value of the magnetic parameter. Further, it is observed that thickness of thermal boundary layer increases with increasing value of the radiation parameter or the Eckert number, whereas, reverse phenomenon observed for the Prandtl number. Moreover, the local Nusselt number decreases with increasing value of the radiation parameter, while reverse behaviour observed for the Prandtl number.