MHD Slip Flow past a Shrinking Sheet

Abstract

An analysis is made for the steady two-dimensional laminar boundary layer flow of a viscous, incompressible, electrically conducting fluid near a stagnation point past a shrinking sheet with slip in the presence of a magnetic field. The governing boundary layer equations are transformed to ordinary differential equations by taking suitable similarity variables and solved numerically by Shooting method. The effects of the various parameters such as velocity ratio parameter, slip parameter, Prandtl number, Eckert number and magnetic parameter for velocity and temperature distributions have been discussed in detail through graphical representation.

Share and Cite:

Chaudhary, S. and Kumar, P. (2013) MHD Slip Flow past a Shrinking Sheet. Applied Mathematics, 4, 574-581. doi: 10.4236/am.2013.43081.

1. Introduction

The viscous flow and heat transfer in the boundary layer region due to a stretching sheet has wide theoretical and technical applications in manufacturing process and in industries such as extraction of polymer sheet, paper production, wire drawing, glass-fiber production, the cooling and drying of paper and textiles. The study of heat transfer and flow field is necessary for determining the quality of the final products of such process. Crane [1] was first who considered steady boundary layer flow of a viscous incompressible fluid over a linearly stretching plate and gave an exact similarity solution in closed analytical form. The effects of heat and mass transfer and magnetic field under various physical conditions have been investigated by several authors such as Chen and Char [2], Chiam [3], Andersson [4], Ariel et al. [5], Jat and Chaudhary [6,7], Wang [8], Fang et al. [9], Nadeem et al. [10], Bhattacharyya and Layek [11] and recently Bhattacharyya et al. [12] studied the slip effects on boundary layer stagnation-point flow and heat transfer towards a shrinking sheet.

Based on the above mentioned investigations and applications, this paper is concerned with a steady, two dimensional stagnation flow of an electrically conducting fluid past a shrinking sheet in the presence of a magnetic field. The results of velocity and temperature distributions for different parameters such as the velocity ratio parameter, the slip parameter, the Prandtl number, the Eckert number and the magnetic parameter were obtained.

2. Formulation of the Problem

Consider the steady two-dimensional laminar flow of a viscous incompressible electrically conducting fluid towards a linearly shrinking sheet such that the sheet is shrinked in its own plane with velocity proportional to the distance from the stagnation point in the presence of an externally applied normal magnetic field of constant strength. The shrinking surface has a linear velocity and uniform temperature, while the velocity of the flow external to the boundary layer is and temperature as shown in Figure 1. Therefore, under the usual boundary layer approximations, the governing equations of motion are:

(1)

(2)

(3)

where is the coefficient of kinematic viscosity, the electrical conductivity, the magnetic permeability, the density, the specific heat at constant pressure, the thermal conductivity and the coefficient of viscosity. The other symbols have their usual

Figure 1. Physical model of two-dimensional stagnation point flow past a shrinking sheet.

meanings.

The boundary conditions are:

(4)

where is a proportionality constant of the velocity of shrinking sheet, is a slip length and is a constant proportional to the free stream velocity for away from the sheet.

3. Analysis of the Velocity and the Thermal Boundary Layers

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 Magnetic parameter,

is the Prandtl number and

is the Eckert number.

The corresponding boundary conditions are:

(11)

where is the velocity ratio parameter and

is the slip parameter.

For numerical solution of the Equations (9) and (10), we apply 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)

4. Results and Discussion

The Equation (14) is that obtained by Bhattacharyya et al. [12] for the non-magnetic case and the remaining equations are ordinary linear differential equations and have been solved numerically by Newton’s shooting method with fourth-order Runge-Kutta integration scheme for various values of the parameters. The velocity and temperature distributions for various values of parameters are shown in Figures 2-6 respectively.

Figure 2. Velocity distribution against for various values of and with.

Figure 3. Velocity distribution against for various values of and with.

Figure 4. Velocity distribution against for various values of and with.

Figure 5. Temperature distribution against for various values of, and with and.

Figure 6. Temperature distribution against for various values of and with, and.

The velocity profiles for different values of the velocity ratio parameter, the slip parameter

and the magnetic parameter are shown in Figures 2-4. It is observed that the velocity boundary layer thickness increases with the increasing values of the velocity ratio parameter and the slip parameter, whereas it decreases as the magnetic parameter increases for a fixed.

The temperature profiles for different values of the velocity ratio parameter, the slip parameterthe Prandtl number, the Eckert number and the magnetic parameter are plotted in Figures 5 and 6. It is observed that for the slip parameter the thermal boundary layer thickness decreases with the increasing values of the velocity ratio parameter, the Prandtl number and the Eckert number and the reverse phenomenon is observed for the magnetic parameter.

5. Conclusions

In this paper, the stagnation flow for two-dimensional electrically conducting fluid past a shrinking sheet with slip boundaries in the presence of a magnetic field is studied. The similarity equations are derived and solved numerically. It is found that the velocity boundary layer thickness increases with the increasing values of the velocity ratio parameter, the slip parameter. Further we observed that the magnetic parameter decreases with the increasing value of the velocity boundary layer thickness but the reverse phenomenon occurs for thermal boundary layer thickness. Also it observed that the thermal boundary layer thickness decreases with increasing values of the velocity ratio parameter, the Prandtl number and the Eckert number.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] L. J. Crane, “Flow Past a Stretching Plate,” Zeitschrift für angewandte Mathematik und Physik, Vol. 21, No. 4, 1970, pp. 645-647. doi:10.1007/BF01587695 [2] C. K. Chen and M. I. Char, “Heat Transfer of a Continuous, Stretching Surface with Suction or Blowing,” Journal of Mathematical Analysis and Applications, Vol. 135, No. 2, 1988, pp. 568-580. [3] T. C. Chiam, “Stagnation-Point Flow towards a Stretching Plate,” Journal of the Physical Society of Japan, Vol. 63, No. 6, 1994, pp. 2443-2444. doi:10.1143/JPSJ.63.2443 [4] H. I. Andersson, “Slip Flow Past a Stretching Surface,” Acta Mechanica, Vol. 158, No. 1-2, 2002, pp. 121-125. doi:10.1007/BF01463174 [5] P. D. Ariel, S. Hayat and S. Asghar, “The Flow of an Elastico-Viscous Fluid Past a Stretching Sheet with Partial Slip,” Acta Mechanica, Vol. 187, No. 1, 2006, pp. 29-35. doi:10.1007/s00707-006-0370-3 [6] R. N. Jat and S. Chaudhary, “MHD Stagnation Flows with Slip,” Il Nuovo Cimento, Vol. 122B, No. 8, 2007, pp. 823-831. [7] R. N. Jat and S. Chaudhary, “Magnetohydrodynamic Boundary Layer Flow near the Stagnation Point of a Stretching Sheet,” Il Nuovo Cimento, Vol. 123B, No. 5, 2008, pp. 555-566. [8] C. Y. Wang, “Stagnation Flow towards a Shrinking Sheet,” International Journal of Non-Linear Mechanics, Vol. 43, No. 5, 2008, pp. 377-382. doi:10.1016/j.ijnonlinmec.2007.12.021 [9] T. Fang, S. Yao, A. Zhang and A. Aziz, “Viscous Flow over a Shrinking Sheet with a Second Order Slip Flow Model,” Communications in Nonlinear Science and Numerical Simulation, Vol. 15, No. 7, 2010, pp. 1831-1842. doi:10.1016/j.cnsns.2009.07.017 [10] S. Nadeem, A. Hussain and M. Khan, “HAM Solutions for Boundary Layer Flow in the Region of the Stagnation Point towards a Stretching Sheet,” Communications in Nonlinear Science and Numerical Simulation, Vol. 15, No. 3, 2010, pp. 475-481. doi:10.1016/j.cnsns.2009.04.037 [11] K. Bhattacharyya and G. C. Layek, “Effects of Suction/ Blowing on Steady Boundary Layer Stagnation-Point Flow and Heat Transfer towards a Shrinking Sheet with Thermal Radiation,” International Journal of Heat and Mass Transfer, Vol. 54, No. 1-3, 2011, pp. 302-307. doi:10.1016/j.ijheatmasstransfer.2010.09.043 [12] K. Bhattacharyya, S. Mukhopadhyay and G. C. Layek, “Slip Effects on Boundary Layer Stagnation-Point Flow and Heat Transfer towards a Shrinking Sheet,” International Journal of Heat and Mass Transfer, Vol. 54, No. 1, 2011, pp. 308-313. doi:10.1016/j.ijheatmasstransfer.2010.09.041