Research Article  Open Access
Effects of Magnetic Field and Thermal Radiation on Stagnation Flow and Heat Transfer of a PowerLaw Fluid over a Shrinking Sheet
Abstract
An analysis is made on the steady twodimensional boundary layer magnetohydrodynamic (MHD) stagnationpoint flow and radiative heat transfer of an electrically conducting powerlaw fluid over a shrinking sheet which is shrunk in its own plane with a velocity proportional to the distance from a fixed point. The similarity transformations are used to transform the boundary layer equations into a system of nonlinear ordinary differential equations which are then solved numerically using shooting technique. It is found that multiple solutions exist for a certain range of the ratio of the shrinking velocity to the free stream velocity (i.e., α) which again depends on the magnetic parameter (M) and the powerlaw index parameter (n). The results pertaining to the present study indicate that as the strength of the magnetic parameter increases, the range of α where similarity solutions exist gradually increases. It is also observed that the temperature at a point decreases with increase in M for the first solution branch, whereas it increases with increase in M for the second solution branch. The reported results are in good agreement with the available published work in the literature.
1. Introduction
From a technological point of view, nonlinear fluid rheology is of special interest and has practical applications. Hence the study of nonNewtonian fluid is very much important. A large number of industrial fluids such as molten plastics, artificial fibres, food stuff, polymers, and slurries are nonNewtonian in their flow characteristics. Many of the inelastic nonNewtonian fluids encountered in chemical engineering processes and biochemical industries are known to follow the empirical Ostwaldde Waele model (see Metzner [1]) or socalled powerlaw model. This model is described by a simple nonlinear equation of state for inelastic fluids which includes linear Newtonian fluids as a special case.
Also, the study of magnetohydrodynamic (MHD) flow of an electrically conducting fluid is of considerable interest in modern metallurgical and metalworking processes. This type of flow has attracted the interest of many researchers due to its application in many engineering problems such as MHD generators, plasma studies, nuclear reactors, and geothermal energy extractions. Also by the application of magnetic field, hydromagnetic techniques are used for the purification of molten metals from nonmetallic inclusions. Therefore the type of problem that we are dealing with is very useful to polymer technology and metallurgy.
Stagnationpoint flow is a topic of significance in fluid mechanics as it appears in virtually all flow fields of science and engineering. In some cases flow is stagnated by a solid wall, while in others a free stagnation point or a line exists in the interior of the fluid domain. Hiemenz [2] was the first to discover that the stagnationpoint flow can be analyzed by the NavierStokes equations for the twodimensional case. The axisymmetric case was studied by Homann [3]. Thereafter numerous investigations were made on the stretching sheet problem with linear stretching in different directions [4–16].
Recently, the boundary layer flow due to a shrinking sheet has attracted considerable interest. One of the common applications of shrinking sheet problems in industries and engineering is shrinking film. In packaging of bulk products, shrinking film is very useful as it can be unwrapped easily with adequate heat. Shrinking problem can also be applied to study the capillary effects in small pores, the shrinkingswell behaviour, and the hydraulic properties of agricultural clay soils. The associated changes in hydraulic and mechanical properties of such soils will seriously hamper predictions of the flow and transport processes which are essential for agricultural development and environmental strategies.
For this flow configuration, the fluid is attracted towards a slot, and the flow is quite different from the stretching case. From a physical point of view, vorticity generated at the shrinking sheet is not confined within a boundary layer, and a steady flow is not possible unless either a stagnationpoint flow towards the sheet is applied, or adequate suction is applied at the sheet. This type of shrinking flow is essentially a backward flow as discussed by Goldstein [17]. For a backward flow configuration, the fluid losses memory of the perturbation introduced by the slot. As a result, the flow induced by the shrinking sheet shows quite distinct physical phenomena from the stretching case.
Miklavcic and Wang [18] investigated both twodimensional and axisymmetric viscous flow induced by a shrinking sheet in the presence of uniform suction. The above shrinking sheet problem with a powerlaw surface velocity was analyzed by Fang [19]. The unsteady viscous flow over a continuously shrinking surface with mass suction was also investigated by Fang et al. [20]. Recently, Fang and Zhang [21] gave an exact solution of MHD boundary layer equations in closed analytical form for the flow over a shrinking sheet in the presence of suction at the surface. Hayat et al. [22] analyzed the MHD flow of a second grade fluid over a shrinking sheet. Steady twodimensional and axisymmetric stagnationpoint flow with heat transfer on a shrinking sheet were investigated by Wang [23]. Mahapatra et al. [24] investigated the flow and heat transfer of the steady twodimensional MHD stagnationpoint flow of an electrically conducting incompressible viscous fluid over a shrinking sheet, the flow being permeated by a uniform transverse magnetic field. Note that with an added stagnationpoint flow to contain the vorticity, similarity solution is possible even in the absence of suction at the surface. Recently, Fang et al. [25] obtained a closed form of analytical solution for the flows of nonNewtonian powerlaw fluids over a shrinking sheet with mass transfer.
Owing to wide range of technological applications, flow of nonNewtonian powerlaw fluid over a stretching/shrinking sheet has attracted many attention. To the best of author’s knowledge, described the frequent mentioned literature, there is no investigation on the MHD stagnation point flow and radiative heat transfer of powerlaw fluid over a shrinking sheet. In this paper, we investigate the steady twodimensional MHD stagnationpoint flow and heat transfer of an electrically conducting powerlaw fluid over a shrinking sheet, the flow being permeated by a uniform transverse magnetic field. Representative results for wall shear stress and velocity profiles are present for some values of the governing parameters. The temperature distribution in the flow is determined when the surface is held at a constant temperature, and the surface heat flux is also calculated. It is hoped that the obtained results will not only present useful information for applications but also serves as a complement to the previous studies.
2. Flow Analysis
Consider the steady twodimensional stagnationpoint flow of an electrically conducting incompressible powerlaw fluid towards a surface which is shrunk in its own plane with a velocity proportional to the distance from the stagnation point. Here axis runs along the shrinking surface in the direction opposite to sheet motion. A uniform magnetic field of strength is applied in a direction normal to the surface that is parallel to axis. We consider the twodimensional irrotational stagnationpoint flow at infinity given by , , where (>0) is the strength of the stagnationpoint flow. On the sheet, the velocities are , , where (<0) is the shrinking rate (stretching rate if ). Here shrinking of the sheet is along the negative direction of axis. Figure 1 describes the physical model and the coordinate system, where the and axes are measured along the surface of the sheet and normal to it, respectively.
The rheological equation of state for an incompressible powerlaw fluid (Wilkinson [26]) is defined as where is the deviatoric part of the stress tensor . Here is the pressure, and is the rateofstrain tensor given by where is the velocity components. Here >0is called the consistency coefficient of the fluid, and (>0) is the powerlaw index. Fluids obeying constitutive equation (1) are called powerlaw fluids, and for such fluids the linear relationship between the stress tensor and the rateofstrain tensor does not hold. We note that for , the fluid is called pseudoplastic powerlaw fluid; for , it is called dilatant powerlaw fluid, and when , the fluid is the Newtonian fluid.
Using boundary layer approximation, the equations for the steady MHD twodimensional flow are in usual notation where is the density and is the pressure of the fluid. The last term in (3) stands for the component of the Lorentz force per unit mass where is the electric current density and is the imposed uniform magnetic field (acting along direction). In writing (3), we have neglected the induced magnetic field in comparison with the imposed magnetic field since the magnetic Reynolds number for the flow is assumed to be very small. Such an assumption is justified for flow of electrically conducting fluids such as liquid metals, for example, mercury and liquid Sodium (Shercliff [27]).
The equation of continuity is Clearly the electric current in this twodimensional flow acts parallel to the axis which is normal to plane. So by Ohm’s law where is the electrical conductivity (assumed constant) of the fluid and is the electric field along axis. It is assumed that the electric current lines in this twodimensional configuration are closed in a manner consistent with the experiment. Now for the flow in the steady state, Maxwell’s equation gives where is the electric field which acts parallel to the direction. This gives and , so that is a function of only.
Since the induced magnetic field is neglected in view of our assumption , electric current in the flow is determined from Ohm’s law and not from , being the magnetic permeability. But the consequence of this equation must be satisfied (Shercliff [27]). Since the physical variables are functions of and , then gives from (6) that = constant. Then using (6), we find from (3) that Now since from (4), the pressure is independent of , and the pressure gradient can now be obtained from (8) in the free stream as where is the free stream velocity and use is made of the fact = constant. Hence eliminating from (8) and (9), we get Note that the electric field does not directly affect the boundary layer equation, although it has an effect on the relation between the free stream velocity and the pressure distribution given by (9). Thus the governing equations for the velocity components and are (5) and (10). In the present problem we have . Hence from (1) and (3), we have the shear stress as Thus the momentum equation (10) reduces to the form The boundary conditions are where is for stretching and is for shrinking. We introduce the similarity transformations where the stream function is defined by the following relations:
With the values of and obtained from (14)(15), we find that the equation of continuity (5) is identically satisfied. Substituting (14)(15) into (12), we obtain the following nonlinear ordinary differential equation: where ) is the dimensionless magnetic parameter and a prime denotes differentiation with respect to the independent variable . From (13), we get the dimensionless boundary conditions as Equation (16) along with the boundary conditions (17) holds for both (i.e., for stretching sheet) and (i.e., for shrinking sheet). The skin friction coefficient at the wall is given by where is the local Reynolds number based on the sheet velocity .
3. Heat Transfer
To determine the temperature distribution in the above flow field, we solve the following energy equation using boundary layer approximation, neglecting viscous and ohmic dissipation, where , , and denote the thermal diffusivity, the specific heat at constant pressure, and the temperature of the fluid, respectively, and is the radiative heat flux. Using Rosseland’s approximation for radiation (Brewster [28]), we can write where is the StefanBoltzmann constant and is the mean absorption coefficient. Assuming the temperature difference within the flow is such that may be expanded in a Taylor series about and neglecting higher order terms, we get . Hence (19) becomes The appropriate boundary conditions are where and are constants with . We introduce the dimensionless temperature as
Now the energy equation (19), with the above boundary condition, becomes where is the generalized Prandtl number and = is the thermal radiation parameter. The boundary conditions for are obtained from (22) and (23) as It is interesting to note that when (i.e., ), then (16) together with the boundary conditions (17) admits of an exact analytical solution . Hence in this case, the energy equation (24) reduces to the form where is a positive constant. Solution of (26) can be written as where and are two constants of integration and denotes the error function. Using the boundary conditions (25), we have the solution as
4. Numerical Solutions
The governing momentum equation (16) along with the boundary conditions (17) and the thermal equation (24) together with the boundary conditions (25) are solved numerically. To do this, we first transform the nonlinear differential equation (16) to a system of three firstorder differential equations, which are solved by means of a standard fourthorder RungeKutta scheme coupled with a conventional shooting procedure. The step size and the convergence criteria are used in this program. In practice, the value must be replaced by an approximation , where the value of is chosen; in such a way, all profiles satisfy the far field boundary conditions asymptotically. We ran our bulk computations with the value , which was sufficient to achieve the far field boundary conditions asymptotically for all values of the parameters considered. Then a Newton iteration procedure is employed to assure quadratic convergence of the iterations to satisfy the boundary condition .
Once the numerical values of ) are known, we have solved the thermal equation (24) together with the boundary condition (25) using the same technique for obtaining the numerical values of . An explanation of this algorithm is given in Gladwell and Sayers [29].
5. Results and Discussion
In order to validate the method used in this study and to judge the accuracy of the present analysis, comparison with the results of Wang [23] for the skin friction coefficient for Newtonian fluid (taking ) in the absence of magnetic field () for stagnation point flow towards a stretching/shrinking sheet is made and found excellent agreement, as shown in Tables 1 and 2. Therefore we are confident that the present results are accurate.


From a physical point of view, vorticity generated at the shrinking sheet is not confined within a boundary layer, and a steady solution is not possible unless sufficient stagnation flow is added to the free stream. Even with an added stagnation flow, a steady solution is possible only when the ratio of shrinking velocity and free stream velocity (i.e., ) becomes less than a certain numerical value which again depends on the magnetic parameter () and the powerlaw index parameter (). Note that (16) is nonlinear. It may be possible to have multiple or bifurcated solutions. Our numerical results reveal that for Newtonian fluid (i.e., ), the solution of (16) satisfying boundary conditions (17) is unique for . But dual solutions exist for , and no similarity solution is found for . These results agree well with those of (Wang [23]) who studied steady stagnationpoint flow over a shrinking sheet. The ranges of for unique solution, dual solutions, and nonexistence of the solution for different values of and are present in concise form in Table 3. From Table 3 it is seen that as the strength of the magnetic parameter increases, the range of where similarity solutions exist gradually increases. The novel result which emerges from the above analysis is that the range of where dual solutions exist progressively decreases with increase in . It is also observed that as the powerlaw index increases, the range of where two solutions exist decreases.

Figure 2 shows the variation of and the horizontal component of velocity, with for several values of the magnetic parameter with fixed values of and for the first and second solution branches. In the following discussion, we identify the first solution branch as the solution with higher values of for a given and second solution branch as the lesser values of for the same (in the dual solution range). It is observed that for the first solution branch, the horizontal velocity component at a point increases with increase in , but for the second solution branch it decreases with increase in . Figure 3 shows the variation of with for different values of the powerlaw index for fixed value of the shrinking rate parameter and for both the solution branches. The figure shows that decreases with increase in for the first solution branch, but for the second solution branch it increases with increase in .
Figure 4 displays the variation of with for different values of for both the solution branches with fixed values of , , and . It is observed that temperature at a point decreases with increase in for the first solution branch and temperature at a point increases with increase in for the second solution branch. From a physical point of view, this follows from the fact that as the magnetic parameter increases, the boundary layer becomes thicker for the first branch, while it becomes thinner for the second branch. Figure 5 displays the variation of with for several values of α (<0) with fixed values of , , and for the first and the second solution branches. The figure reveals that temperature at a point increases with increase in for the first branch and decreases with increase in for the second branch. From a physical point of view, this follows from the fact that the boundary layer becomes thicker as increases for the first branch and it becomes thinner for the second branch as increases.
Figure 6 shows the variation of with for several values of Pr with fixed values of , , and for both the solution branches. The figure indicates that as Pr increases, temperature at a point decreases for the first solution branch. For the second branch temperature at a point increases with increase in Pr in a certain neighbourhood of the shrinking sheet but beyond this neighbourhood temperature at a point decreases with increase in Pr. Figure 7 shows the variation of temperature profiles for various values of the radiation parameter for fixed values of the other parameters. The figure reveals that an increase in radiation parameter leads to an increase in the temperature at a point for both the solution branches except in a small region (in case of second solution).
Figure 8 displays the variation of the dimensionless surface heat flux with α (<0) for different values of the magnetic parameter . The figure reveals that for fixed values of and Pr, the surface heat flux increases with increase in for the first solution branch. Interestingly for the second solution branch for the same values of and Pr, the surface heat flux is almost zero for different values of . Figure 9 shows the variation of the dimensionless surface heat flux with for different values of Pr with fixed and . It is seen that up to certain value of , increases with increase in Pr, and beyond this , decreases with increase in Pr for the first solution branch. Again for the second solution branch, decreases as Pr increases. It is to be noted that for higher values of Prandtl number, is almost zero for the second solution branch.
6. Conclusion
The present paper deals with the analysis of the steady twodimensional MHD stagnationpoint flow and heat transfer of a powerlaw fluid over a continuously shrinking surface in the presence of thermal radiation. In comparison to the stretching sheet, the shrinking sheet has some interesting characteristics. Greatly different solution behaviour with multiple solution branches (for certain value of the shrinking rate parameter) is found compared with the corresponding stretching sheet problem. The following conclusions can be drawn from the numerical results obtained. (i)As the strength of the magnetic parameter increases, the range of where similarity solutions exist gradually increases. (ii)Dual solutions exist for certain parameter domain that depends on , , and . (iii)Velocity at a point increases with increase in for the first solution, and it decreases with increase in for the second solution. (iv)Temperature at a point decreases with increase in for the first solution and increases with increase in for the second solution. (v)An increase in radiation parameter leads to an increase in the temperature at a point.
Acknowledgment
The author appreciates the comments of the reviewer, which have led to improvement of the paper.
References
 A. B. Metzner, “Heat transfer in nonNewtonian fluids,” Advances in Heat Transfer, vol. 2, pp. 357–397, 1965. View at: Publisher Site  Google Scholar
 K. Hiemenz, “Die Grenzschicht in Einem in Dem Gleichformingen Flussigkeitsstrom Eingetauchten Gerade Kreiszlinder,” Dingler's Polytechnic Journal, vol. 326, pp. 321–410, 1911. View at: Google Scholar
 F. Homann, “Der Einfluss Grosser Zahighkeit Bei Der Stromung um den Zylinder und um die Kugel,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 16, pp. 153–164, 1936. View at: Google Scholar
 L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik, vol. 21, pp. 645–647, 1970. View at: Google Scholar
 K. B. Pavlov, “Magnetohydrodynamic flow of an incompressible viscous fluid caused by the deformation of a plane surface,” Magnitnaya Gidrodinamika, vol. 4, pp. 146–147, 1974. View at: Google Scholar
 H. I. Andersson, “MHD flow of a viscoelastic fluid past a stretching surface,” Acta Mechanica, vol. 95, no. 1–4, pp. 227–230, 1992. View at: Publisher Site  Google Scholar
 H. I. Andersson and B. S. Dandapat, “Flow of a powerlaw fluid over a stretching sheet,” Stability and Applied Analysis of Continuous Media, vol. 1, pp. 339–347, 1991. View at: Google Scholar
 H. I. Andersson, K. H. Bech, and B. S. Dandapat, “Magnetohydrodynamic flow of a powerlaw fluid over a stretching sheet,” International Journal of NonLinear Mechanics, vol. 27, no. 6, pp. 929–936, 1992. View at: Google Scholar
 T. C. Chiam, “Stagnationpoint flow towards a stretching plate,” Journal of the Physical Society of Japan, vol. 63, no. 6, pp. 2443–2444, 1994. View at: Google Scholar
 K. Vajravelu, “Flow and heat transfer in a saturated porous medium over a porous medium over a stretching,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 74, pp. 605–614, 1994. View at: Google Scholar
 K. Vajravelu and J. Nayfeh, “Convective heat transfer at a stretching sheet,” Acta Mechanica, vol. 96, no. 1–4, pp. 47–54, 1993. View at: Publisher Site  Google Scholar
 T. R. Mahapatra and A. S. Gupta, “Heat transfer in stagnationpoint flow towards a stretching sheet,” Heat and Mass Transfer, vol. 38, no. 6, pp. 517–521, 2002. View at: Publisher Site  Google Scholar
 S. Abel, K. V. Prasad, and A. Mahaboob, “Buoyancy force and thermal radiation effects in MHD boundary layer viscoelastic fluid flow over continuously moving stretching surface,” International Journal of Thermal Sciences, vol. 44, no. 5, pp. 465–476, 2005. View at: Publisher Site  Google Scholar
 T. R. Mahapatra, S. K. Nandy, and A. S. Gupta, “Magnetohydrodynamic stagnationpoint flow of a powerlaw fluid towards a stretching surface,” International Journal of NonLinear Mechanics, vol. 44, no. 2, pp. 123–128, 2009. View at: Publisher Site  Google Scholar
 K. V. Prasad, D. Pal, and P. S. Datti, “MHD powerlaw fluid flow and heat transfer over a nonisothermal stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2178–2189, 2009. View at: Publisher Site  Google Scholar
 K. V. Prasad and K. Vajravelu, “Heat transfer in the MHD flow of a power law fluid over a nonisothermal stretching sheet,” International Journal of Heat and Mass Transfer, vol. 52, no. 2122, pp. 4956–4965, 2009. View at: Publisher Site  Google Scholar
 S. Goldstein, “On backward boundary layers and flow in converging passages,” Journal of Fluid Mechanics, vol. 21, no. 1, pp. 33–45, 1965. View at: Publisher Site  Google Scholar
 M. Miklavcic and C. Y. Wang, “Viscous flow due to a shrinking sheet,” Quarterly of Applied Mathematics, vol. 64, no. 2, pp. 283–290, 2006. View at: Google Scholar
 T. Fang, “Boundary layer flow over a shrinking sheet with powerlaw velocity,” International Journal of Heat and Mass Transfer, vol. 51, no. 2526, pp. 5838–5843, 2008. View at: Publisher Site  Google Scholar
 T. Fang, J. Zhang, and S. Yao, “Viscous flow over an unsteady shrinking sheet with mass transfer,” Chinese Physics Letters, vol. 26, no. 1, Article ID 014703, 2009. View at: Publisher Site  Google Scholar
 T. Fang and J. Zhang, “Closedform exact solutions of MHD viscous flow over a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2853–2857, 2009. View at: Publisher Site  Google Scholar
 T. Hayat, Z. Abbas, and M. Sajid, “On the analytic solution of magnetohydrodynamic flow of a second grade fluid over a shrinking sheet,” Journal of Applied Mechanics, vol. 74, no. 6, pp. 1165–1171, 2007. View at: Publisher Site  Google Scholar
 C. Y. Wang, “Stagnation flow towards a shrinking sheet,” International Journal of NonLinear Mechanics, vol. 43, no. 5, pp. 377–382, 2008. View at: Publisher Site  Google Scholar
 T. R. Mahapatra, S. K. Nandy, and A. S. Gupta, “Momentum and heat transfer in MHD stagnationpoint flow over a shrinking sheet,” Journal of Applied Mechanics, vol. 78, no. 2, Article ID 021015, 2011. View at: Publisher Site  Google Scholar
 T. Fang, H. Tao, and Y. Zhong, “NonNewtonian powerlaw flud flow over a shrinking sheet,” Chinese Physics Letters, vol. 29, no. 11, Article ID 114703, 2012. View at: Publisher Site  Google Scholar
 W. L. Wilkinson, NonNewtonian Fluids, Pergamon Press, London, UK, 1960.
 J. A. Shercliff, A Textbook of Magnetohydrodynamics, Pergamon Press, Oxford, UK, 1965.
 M. Q. Brewster, Thermal Radiative Transfer Properties, John Wiley & Sons, New York, NY, USA, 1972.
 I. Gladwell and D. K. Sayers, Computational Techniques for Ordinary Differential Equations, Academic Press, London, UK, 1980.
Copyright
Copyright © 2013 Samir Kumar Nandy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.