Research Article  Open Access
Tapas Ray Mahapatra, Sabyasachi Mondal, Dulal Pal, "Heat Transfer due to Magnetohydrodynamic StagnationPoint Flow of a PowerLaw Fluid towards a Stretching Surface in the Presence of Thermal Radiation and Suction/Injection", International Scholarly Research Notices, vol. 2012, Article ID 465864, 9 pages, 2012. https://doi.org/10.5402/2012/465864
Heat Transfer due to Magnetohydrodynamic StagnationPoint Flow of a PowerLaw Fluid towards a Stretching Surface in the Presence of Thermal Radiation and Suction/Injection
Abstract
An analysis is made on the study of twodimensional MHD (magnetohydrodynamic) boundarylayer stagnationpoint flow of an electrically conducting powerlaw fluid over a stretching surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnationpoint in the presence of thermal radiation and suction/injection. The paper examines heat transfer in the stagnationpoint flow of a powerlaw fluid except when the ratio of the free stream velocity and stretching velocity is equal to unity. The governing partial differential equations along with the boundary conditions are first brought into a dimensionless form and then the equations are solved by RungeKutta fourthorder scheme with shooting techniques. It is found that the temperature at a point decreases/increases with increase in the magnetic field when free stream velocity is greater/less than the stretching velocity. It is further observed that for a given value of the magnetic parameter , the dimensionless rate of heat transfer at the surface and decreases/increases with increase in the powerlaw index . Further, the temperature at a point in the fluid decreases with increase in the radiation parameter when free stream velocity is greater/less than the stretching velocity.
1. Introduction
Recently, a considerable amount of interest has been evinced in the study of MHD (magnetohydrodynamic) stagnationpoint flow and heat transfer of a powerlaw fluid towards a stretching surface in view of its various industrial applications such as extrusion of polymer sheets, continuous stretching, manufacturing plastic films, and artificial fibres. In a meltspinning process, the extrudate from the die is usually drawn and simultaneously stretched into a sheet which is then solidified through quenching or gradual cooling by direct contact with water. The qualities of the final product depend on the rate of heat transfer at the stretching surface. It is, therefore, of great importance to know the flow behavior over a stretching surface which determines the rate of cooling.
The results of these studies are of great technical importance; for example, in the prediction of skinfriction, as well as heat/mass transfer near stagnation regions of bodies in high speed flows. Sapunkov [1] investigated the twodimensional orthogonal stagnationpoint flow of an incompressible electrically conducting powerlaw fluid towards a rigid surface in the presence of a uniform transverse magnetic field. Crane [2] studied the steady flow of an incompressible viscous fluid over a deformable sheet, the flow being caused solely by the stretching of the sheet in its own plane with a velocity varying linearly with the distance from a fixed point. Pavlov [3] gave an exact similarity solution to the MHD boundary layer equations for the steady twodimensional flow of an electrically conducting incompressible fluid due to the stretching of a plane elastic surface in the presence of a uniform transverse magnetic field. A comprehensible discussion on the potential engineering application of nonNewtonian powerlaw electrically conducting fluids permeated by magnetic fields was presented by Martinson and Pavlov [4]. Heat transfer in the flow over a stretching surface maintained at constant as well as variable temperature was analyzed by P. S. Gupta and A. S. Gupta [5] and Carragher and Crane [6]. Massoudi and Ramezan [7] used a perturbation technique to study the stagnation point flow and heat transfer of a second grade nonNewtonian fluid. Their analysis is valid only for small values of the parameter that determines the behavior of the nonNewtonian fluid. Later, Massoudi and Ramezan [8] extended the problem to nonisothermal surfaces. MHD flow of a powerlaw fluid over a stretching surface was examined by Andersson et al. [9]. Garg [10] improved the solution obtained by Massoudi and Ramezan [7] by numerically computing the flow characteristics for any value of the nonNewtonian parameter using a pseudosimilarity solution. Mahapatra and Gupta [11] analyzed the steady twodimensional orthogonal stagnationpoint flow of an incompressible viscous electrically conducting fluid towards a stretching surface. Steady twodimensional orthogonal stagnationpoint flow of an incompressible viscous fluid towards a stretching surface was studied by Mahapatra and Gupta [12] who also investigated heat transfer in this flow. Many authors presented some mathematical results on flow and heat transfer and a good amount of references can be found in the papers by Liao and Pop [13] and Nazar et al. [14]. Prasad et al. [15] studied MHD powerlaw fluid flow and heat transfer over a nonisothermal stretching sheet. Recently, Mahapatra et al. [16] studied the steady twodimensional orthogonal stagnationpoint flow of an electrically conducting powerlaw fluid towards a stretching surface when the surface is stretched in its own plane with velocity proportional to the distance from the stagnationpoint.
A new dimension is added to the study of flow and heat transfer over a stretching sheet by consisting the effect of thermal radiation. It is well known that the effect of thermal radiation is important in space technology and high temperature processes. Thermal radiation also plays an important role in controlling heat transfer process in polymer processing industry. The effect of radiation on heat transfer problems has been studied by Hossain and Takhar [17]. Recently, Pal [18] has investigated heat and mass transfer in stagnationpoint flow towards a stretching surface in the presence of buoyancy force and thermal radiation. The objective of this paper is to study the heat transfer in steady twodimensional orthogonal stagnationpoint flow of an electrically conducting powerlaw fluid towards a stretching surface in the presence of a magnetic field and thermal radiation. Temperature of the wall and temperature of free stream fluid are kept constant in such a way that the former is higher than the later.
2. Flow Analysis
Consider twodimensional steady stagnationpoint flow of an electrically conducting powerlaw fluid towards a flat surface coinciding with the plane in the presence of a uniform transverse magnetic field , the flow being confined to the region . Two equal and opposite forces are applied on the stretching surface along the axis so that the surface is stretched keeping the origin fixed as shown in Figure 1. The MHD equations for steady twodimensional stagnationpoint flow in the boundary layer towards the stretching surface are as follows: where the induced magnetic field is neglected (which is justified for MHD flow at small magnetic Reynolds number by Shercliff [19]). It is also assumed that the external electric field is zero and the electric field due to polarization of charges is negligible. Here and are the velocity components along the and directions, respectively. Further , , , and denote density of the fluid, electrical conductivity, applied magnetic field, and shear stress, respectively. In (2), stands for the stagnationpoint velocity in the inviscid free stream. Here the stress tensor is defined by (Mahapata and Gupta [12]) where denotes the rate of strain tensor, is the consistency coefficient and is the powerlaw index. Fluids obeying constitutive equation (3) are called powerlaw fluid. It is called pseudoplastic powerlaw fluid when and dilatant powerlaw fluid when since the apparent viscosity decreases (shearthinning) or increases with the increase in shear rate (shearthickening) accordingly as or . It is called Newtonian fluid when . Here we denote and . In the present problem, we have when and when , which gives the shear stress as When , the momentum equation (2) becomes and when , (2) becomes Here is the externally imposed magnetic field in the direction. The appropriate boundary conditions are where and are positive constants and , where is constant.
We introduce the variables as where stream function is defined by the following relations Using (9), (6) and (7) now become, when as and when , we have From (8), we get the dimensionless boundary conditions as where the prime denotes differentiation with respect to .
Here that is, injection, suction and absence of suction/injection parameter and (=) is the magnetic parameter. The skin friction coefficient at the wall is given by where is local Reynolds number based on the sheet velocity . The equations (11) and (12) are solved numerically as described in Mahapatra et al. [16] for different values of powerlaw index in the range and for different values of and .
3. Heat Transfer
Let us consider the heat transfer equation in the stagnationpoint flow of electrically conducting powerlaw fluid towards a stretching surface. Using boundary layer approximations and neglecting viscous and Ohmic dissipation, the transfer of heat in the steady twodimensional flow of a powerlaw fluid can be expressed in the form of the energy equation in the presence of thermal radiation is given by where , , , and denote temperature, thermal diffusivity, the radiation heat flux, and specific heat at constant pressure of the fluid. In addition, the radiation heat flux is considered according to Rosseland approximation such that where and are the StefanBoltzmann constant and the mean absorption coefficient, respectively. As done by Raptis [20], the fluidphase temperature differences within the flow are assumed to be sufficiently small so that may be expressed as a linear function of temperature. This is done by expanding in a Taylor series about the free stream temperature and neglecting higher order terms to yield, Using (16) and (17), we obtain from (15) The appropriate boundary conditions are where and are constants with . Now the energy equation (18), with the above boundary condition, becomes where subject to Here and are the radiation parameter and the local Prandtl number, respectively. The Equations (20) and (22) are solved numerically using RungeKuttta fourthorder method with shooting technique.
4. Discussion of the Results
The nonlinear boundary value problem (11), (12), and (20) cannot be solved in closed form, so these equations are solved numerically using fourthorder RungeKutta method with shooting technique. Extensive calculations have been performed to obtain the flow and temperature profiles for various values of physical parameters such as , , , , and so forth as taken by Mahapatra and Gupta [11], Mahapatra et al. [16]. We have compared numerical result for with those obtained by Mahapatra et al. [16] in Table 1 for various values of , , and in absence of suction/injection. Table 2 shows the computed values of for various values of , , and when and . From this table, it is seen that is negative when and is positive when . It is also seen that decreases with increase in for fixed for both and . It is also observed that increases with increase in for fixed for both and . The computed values of for various values of , , and with , when are presented in Table 3. From this table, it is seen that increases with increasing when for fixed values of , , and for . The same trend is observed by increasing the value of when for and opposite trend is observed when for . It is also seen that increases with increasing when for fixed values of , and when and opposite trend is seen for when .

Figure 2 shows the variation of component of velocity for different values of and for fixed values of and and . It is clear from this figure that when , the velocity at a point in the neighbourhood of stretching plate decreases with the increase in the suction parameter , whereas when the opposite trend is observed. In both cases velocity at infinity coincides with the values of for all values of which is also seen from this figure. Figure 3 shows that the temperature at a given point decreases with increase in for both the cases when and for fixed values of , , , and . Further, this figure shows that the temperature at a given point decreases with increasing for fixed values of , , , and .
Figure 4 shows the variation of component velocity profile for different values of and for fixed values of . It is clear from this figure that when , the velocity at a point in a neighbourhood of the stretching plate increases with increase in . When the opposite trend is observed. In both cases velocity at infinity coincides with the values of for all which is also seen from this figure. It can also be seen that when , the flow has a boundary layer structure and the thickness of this boundary layer structure increases with increase in . It is interesting to note that when , the flow has an inverted boundary layer structure. This is due to the fact when , the stretching velocity of the surface () exceeds the velocity of the external stream ().
Figures 5, 6, and 7 demonstrate that the temperature at a given point decreases with increasing for both and , this is due to the fact that the thermal boundary layer thickness decreases with increase in the value of Prandtl number for fixed values of , , and when . Figures 8, 9, and 10 show the variation of temperature profiles for various values of and in the thermal boundary layer for fixed values of , , and when , respectively. From this figures, it is observed that the temperature profiles decreases with increasing for both and . Further, these figures show that the temperature at a given point decreases with increasing for fixed values of , , , and . Figures 11, 12, and 13 reveal that the temperature at a given point decreases with increasing radiation parameter for both and for fixed values of , , and when . Due to the fact that the divergence of the radiation heat flux decreases with increases in the Rosseland radiative absorption coefficient increases (see (16)) which in turn decreases the rate of radiative heat transfer to the fluid which causes the fluid temperature to decrease. Figures 11â€“13 also depict that the temperature at a given point decreases with increasing for the fixed values of , , , and . Figures 11â€“13 also depict that the temperature at a given point decreases with increasing for the fixed values of , , , and .
Figure 14 shows that the temperature at a given point increases with increase in when (=0.2) for fixed values of , , and . This is due to the fact that as the strength of magnetic field increases the Lorentz force which opposes the flow in the boundary layer that is responsible for increase of the temperature in the thermal boundary layer. Further, it is observed from this figure that the temperature at a given point decreases with increase in when for fixed values of , , and . This paradoxical result can be explained by the fact that when , the velocity of the stretching sheet exceeds the velocity of the inviscid stream and an inverted boundary layer is formed near the surface. Thus one would expect that the velocity at a point in this boundary layer increases with increase in and thereby decreases the temperature in the thermal boundary layer.
5. Conclusions
The present paper deals with the analysis of the heat transfer in the steady twodimensional MHD stagnationpoint flow of an electrically conducting powerlaw fluid over a stretching surface. The sheet is stretched in its plane with a velocity , being the distance along the plate from the stagnation point. The following conclusions can be drawn from the numerical results obtained. (i)When , the velocity at a point in a neighbourhood of the stretching plate increases with increase in the powerlaw index but velocity at a point decreases with the increase in powerlaw index, beyond certain distance from stretching sheet, whereas opposite trend is observed when