**Journal of Applied Mathematics and Physics**

Vol.03 No.07(2015), Article ID:58436,9 pages

10.4236/jamp.2015.37112

Unsteady Magnetohydrodynamic Boundary Layer Flow near the Stagnation Point towards a Shrinking Surface

Santosh Chaudhary, Pradeep Kumar

Department of Mathematics, Malaviya National Institute of Technology, Jaipur, India

Email: d11.santosh@yahoo.com, pradeep17matrix@gmail.com

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 23 June 2015; accepted 27 July 2015; published 30 July 2015

ABSTRACT

The unsteady two-dimensional, laminar flow of a viscous, incompressible, electrically conducting fluid towards a shrinking surface in the presence of a uniform transverse magnetic field is studied. Taking suitable similarity variables, the governing boundary layer equations are transformed to ordinary differential equations and solved numerically by a perturbation technique for a small magnetic parameter. The effects of various parameters such as unsteadiness parameter, velocity parameter, magnetic parameter, Prandtl number and Eckert number for velocity and temperature distributions along with local Skin friction coefficient and local Nusselt number have been discussed in detail through numerical and graphical representations.

**Keywords:**

Unsteady Flow, Magnetohydrodynamic, Boundary Layer, Stagnation Point, Shrinking Surface

1. Introduction

Stagnation flow of an incompressible viscous fluid over a shrinking sheet has many important practical applications in engineering and industrial processes, such as the extrusion of a polymer in a melt-spinning process, continuous casting of metals, the aerodynamic extrusion of plastic sheets, the cooling of metallic sheets or electronic chips and many others. In all these cases, a study of fluid flow and heat transfer has important significance because the quality of the final product depends on the rate of cooling and the process of stretching.

In recent years, the boundary layer flow due to a shrinking sheet has attraction of many researchers because of its useful applications. A very interesting example in which the shrinking sheet situation occurs is of a rising shrinking balloon. Shrinking film is also a common application of shrinking sheet problems in engineering and industries. Shrinking film is very useful in packaging of bulk products because it can be unwrapped easily with adequate heat.

From the stretching case, the flow of shrinking sheet is different and the fluid is attracted towards a slot. Physically, the generated velocity at shrinking sheet has an unsteady flow due to the application of inadequate suction and is not confined within the boundary layer.

In view of all these applications, Sakiadis [1] initiated the study of boundary layer flow over a continuous solid surface moving with constant speed. Later Crane [2] considered the problem of the flow over a linearly stretching sheet in an ambient fluid and gave a similarity solution in closed analytical form for the steady two-dimensional problem. Gupta and Gupta [3] , and Vleggaar [4] have investigated the solution of stretching flow problems at the constant surface temperature while Soundalgekar and Ramana [5] and Grubka and Bobba [6] have analysed the solution of stretching flow problems with a variable surface temperature. Many researchers such as Magyari and Keller [7] , Elbashbeshy and Bazid [8] , Jat and Chaudhary [9] -[11] , Bachok et al. [12] and Zheng et al. [13] have analyzed the stretching sheet problems with different aspects of fluid, such as the heat transfer, the permeability of the surface and the unsteadiness flow. Mahapatra and Nandy [14] [15] studied the unsteady stagnation-point flow and heat transfer over an unsteady shrinking sheet. Recently Aly and Vajravelu [16] and Chaudhary and Kumar [17] discussed numerical solutions of boundary layer flow problems over different surfaces in a porous medium. More recently Nandy et al. [18] and Rosca and Pop [19] investigated the unsteady boundary layer flow over a permeable stretching or shrinking surface.

Realizing the increasing technical applications of the magnetohydrodynamic effects, the aim of the present work is concerned with a steady, two-dimensional unsteady stagnation flow of an electrically conducting fluid over a shrinking surface in the presence of a uniform transverse magnetic field.

2. Formulation of the Problem

Consider an unsteady two-dimensional steady flow of a viscous incompressible electrically conducting fluid near a stagnation point over a continuously shrinking surface placed in the plane of a Cartesian coordinate system in the presence of time dependent free stream. The x-axis is taken along the shrinking surface in the direction of motion and y-axis is perpendicular to it. A uniform magnetic field of constant strength is assumed to be applied normal to the shrinking surface (Figure 1). The surface is assumed to be

Figure 1. Coordinate system for the shrinking surface.

highly elastic and is shrinking in the x-direction with a velocity is and surface temperature while the velocity of the flow, external to the boundary layer is and temperature. Therefore, under the usual boundary layer and Boussinesq approximations, the systems of boundary layer equations are given by

(1)

(2)

(3)

where is the coefficient of kinematic viscosity, the electrical conductivity, the magnetic permea-

bility, the density, Pr the Prandtl number, the coefficient of viscosity, the specific heat at constant pressure and the thermal conductivity. The other symbols have their usual meanings.

The boundary conditions are

(4)

where c is a constant, is the shrinking rate and a is the strength of the stagnation point flow.

3. Analysis

The continuity Equation (1) is identically satisfied by stream function, defined as

(5)

For the solution of the momentum and the energy Equations (2) and (3), the following dimensionless variables are defined:

(6)

(7)

(8)

Equations (5) to (8), transform Equations (2) and (3) into

(9)

(10)

where a prime (') denotes differentiation with respect to, is the unsteadiness parameter, is the velocity parameter, is the Magnetic parameter and is the Eckert number.

The corresponding boundary conditions are

(11)

For numerical solution of the Equations (9) and (10), through a perturbation technique, by assuming the following power series in a small magnetic parameter as

(12)

(13)

Substituting Equations (12) and (13) and its derivatives in Equations (9) and (10) and then equating the coefficients of like powers of, we get the following set of equations

(14)

(15)

(16)

(17)

(18)

(19)

with the boundary conditions

(20)

The Equation (14) is obtained by Mahapatra and Nandy [14] for the non-magnetic case and the remaining equations are ordinary linear differential equations and have been solved numerically by Runge-Kutta method of fourth order. The velocity and temperature distributions for various values of parameters are shown in Figure 2 to Figure 6 respectively.

4. Skin Friction and Surface Heat Transfer

The physical quantities of interest, the local skin friction coefficient and the local Nusselt number Nu i.e. surface heat transfer are given by:

(21)

Figure 2. Velocity distribution against η for various values of β and Re_{m} with.

Figure 3. Temperature distribution against η for various values of Re_{m} and Pr with, and.

Figure 4. Temperature distribution against η for various values of β and Re_{m} with, and.

Figure 5. Temperature distribution against η for various values of and Re_{m} with, and.

Figure 6. Temperature distribution against η for various values of Re_{m} and Ec with, and.

and

(22)

which, in the present case can be expressed in the following forms

(23)

and

(24)

where is the local Reynolds number.

Numerical values of the functions and, which are proportional to local skin friction and local heat transfer rate at the surface respectively for various values of the parameter are presented in Table 1 and Table 2.

Table 1. Numerical values of for various values of the parameters β, α & Re_{m}.

Table 2. Numerical values of for various values of the parameters β, α, Re_{m}, Pr & Ec.

5. Results and Discussion

Figure 2 shows the variation of velocity distribution against η for various values of the unsteadiness parameter β, the velocity parameter α and the magnetic parameter Re_{m}. It may be observed that, for the fixed value of the velocity parameter α velocity distribution increases with the decreasing value of the unsteadiness parameter β, and opposite phenomenon occur for the magnetic parameter Re_{m}, for a fixed η.

Figure 3 to Figure 6 show the variation of temperature distribution against η for the various values of the parameters such as the unsteadiness parameter β, the velocity parameter α, the magnetic parameter Re_{m}, the Prandtl number Pr and the Eckert number Ec. From these figures it may be observed that the temperature distribution decreases with increasing values of the unsteadiness parameter β, the velocity parameter α, the magnetic parameter Re_{m}, the Prandtl number Pr and the Eckert number Ec.

In Table 1, the numerical values of the function for various values of the unsteadiness parameter β, the velocity parameter α and the magnetic parameter Re_{m} are given. It may be observed from the table that the boundary values for the non-magnetic flow are the same as those obtained by Mahapatra and Nandy [14] . The value of the function decreases with the increasing values of the unsteadiness parameter β and the magnetic parameter Re_{m} respectively taking other parameters constant and reverse phenomenon occurs for the velocity parameter α.

In Table 2, the numerical values of the function for the different values of the unsteadiness parameter β, the velocity parameter α, the magnetic parameter Re_{m}, the Prandtl number Pr and the Eckert number Ec are given. It may be observed from the table that the boundary values for the non-magnetic flow are same as those obtained by Mahapatra and Nandy [14] . The value of the function increases with the increasing value of the unsteadiness parameter β, considering other parameters constant and same phenomenon occurs for the velocity parameter α, the magnetic parameter Re_{m}, the Prandtl number and the Eckert number Ec. It is further observed that the function decreases with an increasing value of the Prandtl number for fixed other parameters.

6. Conclusions

The present work extends the two-dimensional unsteady stagnation flow of an electrically conducting fluid, over shrinking surface in the presence of magnetic field. Under some special conditions, the problem will reduce the results obtained by previous researchers. The effects of different parameters such as the unsteadiness parameter, the velocity parameter, the magnetic parameter, the Prandtl number and the Eckert number are studied in detail. The velocity as well as thermal boundary layer thickness decreases with the increasing values of the unsteadiness parameter, the velocity parameter, the magnetic parameter, the Prandtl number and the Eckert number whereas in the velocity reverse phenomenon occurs for the magnetic parameter. From the results it can be concluded that skin friction and Nusselt number vary according to the velocity and thermal boundary layers thickness respectively with different parameters.

Cite this paper

SantoshChaudhary,PradeepKumar, (2015) Unsteady Magnetohydrodynamic Boundary Layer Flow near the Stagnation Point towards a Shrinking Surface. *Journal of Applied Mathematics and Physics*,**03**,921-930. doi: 10.4236/jamp.2015.37112

References

- 1. Sakiadis, B.C. (1961) Boundary-Layer Behaviour on Continuous Solid Surfaces: I. Boundary-Layer Equations for Two Dimensional and Axisymmetric Flow. Journal of American Institute of Chemical Engineers (AIChE), 7, 26-28. http://dx.doi.org/10.1002/aic.690070108
- 2. Crane, L.J. (1970) Flow past a Stretching Plate. Zeitschrift für Angewandte Mathematik und Physik, 21, 645-647. http://dx.doi.org/10.1007/BF01587695
- 3. Gupta, P.S. and Gupta, A.S. (1977) Heat and Mass Transfer on a Stretching Sheet with Suction or Blowing. The Canadian Journal of Chemical Engineering, 55, 744-746.

http://dx.doi.org/10.1002/cjce.5450550619 - 4. Vleggaar, J. (1977) Laminar Boundary-Layer Behaviour on Continuous, Accelerating Surfaces. Chemical Engineering Science, 32, 1517-1525. http://dx.doi.org/10.1016/0009-2509(77)80249-2
- 5. Soundalgekar, V.M. and Ramana Murty, T.V. (1980) Heat Transfer in Flow past a Continuous Moving Plate with Variable Temperature. Warme-und Stoffubertragung, 14, 91-93.

http://dx.doi.org/10.1007/BF01806474 - 6. Grubka, L.J. and Bobba, K.M. (1985) Heat Transfer Characteristics of a Continuous, Stretching Surface with Variable Temperature. Journal of Heat Transfer, 107, 248-250.

http://dx.doi.org/10.1115/1.3247387 - 7. Magyari, E. and Keller, B. (1999) Heat and Mass Transfer in the Boundary Layers on an Exponentially Stretching Continuous Surface. Journal of Physics D: Applied Physics, 32, 577-585.

http://dx.doi.org/10.1088/0022-3727/32/5/012 - 8. Elbashbeshy, E.M.A. and Bazid, M.A.A. (2003) Heat Transfer over an Unsteady Stretching Surface with Internal Heat Generation. Applied Mathematics and Computation, 138, 239-245.

http://dx.doi.org/10.1016/S0096-3003(02)00106-6 - 9. Jat, R.N. and Chaudhary, S. (2008) Magnetohydrodynamic Boundary Layer Flow near the Stagnation Point of a Stretching Sheet. Il Nuovo Cimento, 123B, 555-566.
- 10. Jat, R.N. and Chaudhary, S. (2009) Unsteady Magnetohydrodynamic Boundary Layer Flow over a Stretching Surface with Viscous Dissipation and Joule Heating. Il Nuovo Cimento, 124B, 53-59.
- 11. Jat, R.N. and Chaudhary, S. (2010) Radiation Effects on the MHD Flow near the Stagnation Point of a Stretching Sheet. Zeitschrift für Angewandte Mathematik und Physik, 61, 1151-1154.

http://dx.doi.org/10.1007/s00033-010-0072-5 - 12. Bachok, N., Ishak, A. and Nazar, R. (2011) Flow and Heat Transfer over an Unsteady Stretching Sheet in a Micropolar Fluid. Meccanica, 46, 935-942. http://dx.doi.org/10.1007/s11012-010-9353-5
- 13. Zheng, L., Wang, L. and Zhang, X. (2011) Analytic Solutions of Unsteady Boundary Flow and Heat Transfer on a Permeable Stretching Sheet with Non-Uniform Heat Source/Sink. Communications in Nonlinear Science and Numerical Simulation, 16, 731-740.

http://dx.doi.org/10.1016/j.cnsns.2010.05.022 - 14. Mahapatra, T.R. and Nandy, S.K. (2011) Unsteady Stagnation-Point Flow and Heat Transfer over an Unsteady Shrinking Sheet. International Journal of Applied Math and Mechanics, 7, 11-26.
- 15. Mahapatra, T.R. and Nandy, S.K. (2013) Slip Effects on Unsteady Stagnation-Point Flow and Heat Transfer over a Shrinking Sheet. Meccanica, 48, 1599-1606.

http://dx.doi.org/10.1007/s11012-012-9688-1 - 16. Aly, E.H. and Vajravelu, K. (2014) Exact and Numerical Solutions of MHD Nano Boundary Layer Flows over Stretching Surfaces in a Porous Medium. Applied Mathematics and Computation, 232, 191-204.

http://dx.doi.org/10.1016/j.amc.2013.12.147 - 17. Chaudhary, S. and Kumar, P. (2014) MHD Forced Convection Boundary Layer Flow with a Flat Plate and Porous Substrate. Meccanica, 49, 69-77. http://dx.doi.org/10.1007/s11012-013-9773-0
- 18. Nandy, S.K., Sidui, S. and Mahapatra, T.R. (2014) Unsteady MHD Boundary Layer Flow and Heat Transfer of Nanofluid over a Permeable Shrinking Sheet in the Presence of Thermal Radiation. Alexandria Engineering Journal, 53, 929-937. http://dx.doi.org/10.1016/j.aej.2014.09.001
- 19. Rosca, N.C. and Pop, I. (2015) Unsteady Boundary Layer Flow over a Permeable Curved Stretching/Shrinking Surface. European Journal of Mechanics—B/Fluids, 51, 61-67.

http://dx.doi.org/10.1016/j.euromechflu.2015.01.001