Appl. Math. Mech. -Engl. Ed. DOI 10.1007/s10483-014-1836-9 c Shanghai University and Springer-Verlag Berlin Heidelberg 2014

Applied Mathematics and Mechanics (English Edition)

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet∗ N. S. AKBAR1 ,

Z. H. KHAN2,3 ,

R. U. HAQ4 ,

S. NADEEM4

(1. Department of Basic Science & Humanities, College of Electrical & Mechanical Engineering (CEME), National University of Sciences and Technology, Islamabad 46000, Pakistan; 2. School of Mathematical Sciences, Peking University, Beijing 100871, P. R. China; 3. Department of Mathematics, University of Malakand, Khyber Pakhtunktwo 18800, Pakistan; 4. Department of Mathematics, Quaid-i-Azam University, Islamabad 44000, Pakistan)

Abstract The present article investigates the dual nature of the solution of the magnetohydrodynamic (MHD) stagnation-point flow of a Prandtl fluid model towards a shrinking surface. The self-similar nonlinear ordinary differential equations are solved numerically by the shooting method. It is found that the dual solutions of the flow exist for certain values of the velocity ratio parameter. The special case of the first branch solutions (the classical Newtonian fluid model) is compared with the present numerical results of stretching flow. The results are found to be in good agreement. It is also shown that the boundary layer thickness for the second solution is thicker than that for the first solution. Key words stagnation-point flow, shrinking sheet, Prandtl fluid, magnetohydrodynamic (MHD), dual solution, shooting method Chinese Library Classification O361.3 2010 Mathematics Subject Classification

1

76E25, 76W05

Introduction

A number of studies have been reported in the literature focusing on the stagnation-point flow towards a stretching sheet because of its industrial and engineering applications such as extrusion, paper production, insulating materials, glass drawing, and continuous casting. Hiemenz[1] initiated the two-dimensional stagnation-point flow. He discovered that the NavierStokes equations governing the flow can be transformed into an ordinary differential equation (ODE) of third-order by using the similarity transformation. The boundary layer flow over a stretching surface was studied by Sakiadis[2] . He modeled the laminar boundary-layer behavior on a moving continuous flat surface and presented the numerical solutions for the boundarylayer equations. Later on, this idea was extended by Crane[3] for both linear and exponentially stretching sheets. Chiam[4] discussed the steady two-dimensional stagnation-point flow of a viscous fluid towards a stretching sheet. He made the analysis with the assumption that the sheet was stretched in its own plane with a velocity proportional to the distance from the stagnation point. Wang[5] studied the free convection on a vertical stretching surface. Nazar et al.[6] ∗ Received Jul. 23, 2013 / Revised Dec. 27, 2013 Corresponding author Z. H. KHAN, Ph. D., E-mail: [email protected]

2

N. S. AKBAR, Z. H. KHAN, R. U. HAQ, and S. NADEEM

analyzed the unsteady two-dimensional stagnation-point flow of an incompressible viscous fluid over a flat deformable sheet. The magnetohydrodynamic (MHD) stagnation-point flow towards a stretching vertical sheet was discussed by Ishak et al.[7] . The two-dimensional stagnationpoint flow of a viscoelastic fluid was studied by Sadeghy et al.[8] assuming that the fluid obeyed the upper-convected Maxwell (UCM) model. The boundary-layer hypothesis was used to simplify the equations of motion which were further reduced to a single nonlinear third-order ODE with the idea of the stream function coupled with the technique of the similarity solution. The resulting equations were solved using the Chebyshev pseudo-spectral collocation-point method. Attia[9] made an analysis of the steady hydromagnetic laminar three-dimensional stagnationpoint flow of an incompressible viscous fluid impinging on a permeable stretching surface with heat generation or absorption. MHD flows play an important role in the motion of fluids. Normally, a uniform magnetic field is applied normal to the plate which is maintained at a constant temperature. The steady MHD mixed convection flow of a viscoelastic fluid in the neighborhood of two-dimensional stagnationpoints with a magnetic field has been investigated by Kumari and Nath[10] considering the UCM model. The boundary layer theory could simplify the equations of motion, the induced magnetic field, and the energy into three coupled nonlinear ODEs. These equations were finally solved by using the finite difference method. The results indicated that the increase in the elasticity number causes the reduction in the surface velocity gradient and the surface heat transfer. The dual solutions of the boundary layer flow over moving surfaces are of practical importance in engineering analysis. It gives the possibility to determine the most realistic and physically meaningful solutions. Extensive studies have been carried out to investigate the multiple solutions of the boundary layer flows driven by moving surfaces. Kemp and Acrivos[11] found the dual solutions for a moving-wall boundary layer reverse flow. Riley and Weidman[12] presented multiple solutions for the Falker-Skan equation with a stretching boundary. Ingham[13] presented the dual solutions of a steady mixed convection boundary layer assisting flow over a moving vertical flat plate. Ridha and Curie[14] found the dual solutions for both assisting and opposing flows. Recently, Mahapatra et al.[15] showed that the upper branch solution was always stable whereas the lower branch was unstable. Moreover, the stable branch was physically meaningful. Makinde et al.[16] showed that the dual solutions existed for the stagnation-point flow of an electrically conducting nanofluid towards a vertically stretching sheet. In the present article, we discuss the stagnation-point flow of a Prandtl fluid towards a shrinking sheet with a magnetic field. To the best of the authors’ knowledge, no such investigation has been carried out yet. The main objective of the article is to discuss the dual solutions for both stagnation points and shrinking flows.

2

Mathematical formulation

We consider a two-dimensional stagnation-point flow of an incompressible Prandtl fluid over a wall coinciding with the plane y = 0. The flow is confined to the plane y > 0, and the sheet shrinks linearly along the x-axis. Moreover, we consider the effects of stagnation away from the x-axis, and the uniform magnetic field is applied normal to the fluid flow (see Fig. 1). The extra stress tensor for the Prandtl fluid is defined as[17]   2  2  1  2 Aarcsin C1 ∂∂uy¯¯ + ∂∂¯xv¯ ∂u ¯ , τ=  ∂ u¯ 2  ∂¯v 2  1 ∂ y¯ 2 + ∂y ¯ ∂x ¯

(1)

where A and C are material constants of the Prandtl fluid model, respectively. The flow equations for the Prandtl fluid model after applying the boundary layer approximations can be

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet

Fig. 1

3

Flow configuration and coordinate system

defined as follows: ∂u ∂v + = 0, ∂x ∂y u

(2)

∂u ∂u ∂ue νA ∂ 2 u νA +v = ue + + ∂x ∂y ∂x C ∂y 2 2C 3



∂u ∂y

2

∂ 2 u σB02 (u − ue ), − ∂y 2 ρ

(3)

where u and v are the velocity components along the x- and y-axes, respectively. ν is the kinematic viscosity, B0 is the magnetic field, σ is the permeability of the fluid, and ρ is the density of the fluid. The flow velocity outside the boundary layer (inviscid fluid) is ue (x) = ax, where a > 0. From Eq. (3), we can see that most of the liquid metals possess low conductivity which may lead to small electrical currents generated by the fluid flow in the presence of a magnetic field. In other words, the magnetic Reynolds number is much smaller than one[18–19] . The corresponding boundary conditions are  u = uw (x) = bx, v = vw at y = 0, (4) u → ax as y → ∞, where a > 0 is the constant, and uw is the velocity at wall. Introduce the following similarity transformations:  √ a y, u = axf  (η), v = − aνf (η). η= ν

(5)

By transformations (5), Eqs. (2)−(4) take the forms αf  − (f  )2 + f f  + βf  (f  )2 + M 2 (1 − f  ) + 1 = 0,  f = s, f  = λ = b/a at η = 0, f  → 1 at η → ∞. In the above equations, α = b a

A C

is the Prandtl parameter, β =

a3 x2 A 2Cν is the σB 2 M 2 = ρa0 is

(6) (7) elastic parameter,

λ = is the stretching (λ > 0) or shrinking (λ < 0) parameter, the magnetic field parameter, and s is the mass transfer parameter with s > 0 for suction and s < 0 for injection. After using the boundary layer approximations, the wall shear stress τw can be given by  3 ∂u A A ∂u + . (8) τw = C ∂y 6C 3 ∂y

4

N. S. AKBAR, Z. H. KHAN, R. U. HAQ, and S. NADEEM

The coefficient of the skin friction is defined by Cf =

τw . ρu2w

(9)

In the dimensionless form, the skin friction is defined as 1

Rex2 Cf = (αf  (η) + β(f  (η))3 )η=0 .

3

(10)

Numerical method for solutions

Numerical solutions to the governing ordinary differential equation (6) subject to the boundary conditions (7) are obtained by the shooting method. First, we convert the boundary value problem (BVP) into the initial value problem (IVP) and assume a suitable finite value for the far field boundary condition, i.e., η → ∞, saying η∞ . Then, we set the following first-order system: f  = p,

f  = p = q,

f  = q  =

p2 − f q − M 2 (1 − p) − 1 α + βq 2

(11)

with the boundary conditions f (0) = s,

f  (0) = p (0) = λ.

(12)

To solve (11) with (12) as an IVP, the values for q(0), i.e., f  (0), are needed. However, no such a value is given prior to the computation. The initial guess values of f  (0) are chosen, and the fourth-order Runge-Kutta method is used to obtain a solution. We compare the calculated values of f  (η) at the far field boundary condition η∞ (= 20) with the given boundary condition f  (η) → ∞, and the values of f  (0) are adjusted by the Secant method for better approximation. The step-size is taken as Δη = 0.01, and the accuracy to the fifth decimal place is regarded as the criterion of convergence. It is important to note that the dual solutions are obtained by setting two different initial guesses for the values of f  (0), where both profiles (first and second solutions) satisfy the far field boundary condition (7) asymptotically but with different shapes. To validate the accuracy of the proposed numerical scheme, the obtained results corresponding to the skin-friction coefficient f  (0) (for the Newtonian fluid case, i.e., β = 0) are compared with those available in the literatures[20–22] in Tables 1 and 2 for the stretching case (λ > 0) and the shrinking case (λ < 0). Table 1

Comparison of values of coefficient of skin-friction (with M = 0) for stretching sheet with different values of λ

Present results λ

α = 1, β = 0

Results of Mahapatra and Nandy[22]

Results of Wang[21]

Results of Lok et al.[20]

0.0

1.232 6

1.232 6

1.232 6

−

0.1

1.146 6

1.146 6

1.146 6

−

0.2

1.051 1

1.051 1

1.051 1

−

0.5

0.713 3

0.713 3

0.713 3

0.713 3

1.0

0.000 0

0.000 0

0.000 0

−

2.0

−1.887 3

−1.887 3

−1.887 3

−1.887 3

5.0

−10.264 7

−10.264 7

−10.264 8

−10.264 8

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet Table 2

−0.25 −0.50 −0.75 −1.0 −1.10 −1.15 −1.20 −1.246

4

Comparison of values of coefficient of skin-friction (with M = 0) for shrinking sheet with different values of λ Present results for α=1, β=0 First solution Second solution

λ

5

1.402 1.495 1.489 1.328 1.186 1.082 0.932 0.554

2 6 3 8 7 2 5 3

− − − 0.000 0.049 0.116 0.233 0.554

Results of Ref. [22] First solution Second solution 1.402 1.495 1.489 1.328 1.186 1.082 0.932 0.584

0 2 7 6 2

2 7 3 8 7 2 4 4

− − − 0.000 0.049 0.116 0.233 0.554

0 2 7 6 2

Results and discussion

In the present section, we discuss the effects of the physical parameters such as the suction/injection parameter s, the Prandtl parameter α, the elastic parameter β, and the magnetic field parameter M on both the velocity profile and the skin friction coefficient. The effects of the various fluid flow parameters on the skin friction coefficient are presented in Figs. 2−4. From Fig. 2, it can be observed that for fixed suction and MHD parameters, the local skin friction increases as the elastic parameter (β) decreases. Meanwhile, the existence of the dual solution for the shrinking sheet case is found when λc < λ < 0, where the critical values of λ are found to be λc ≈ −4.69, −4.78, and −5.47 for β = 0.3, 0.2, and 0.1, respectively. Moreover, as the elastic parameter increases, the curve for the local skin friction first decreases gradually, and then decreases to zero at λc = −4. It is also observed that for decreasing values of the elastic parameter β, both the branches of the skin friction coefficient exhibit the same increasing tendency.

Fig. 2

Variations of skin friction with λ for different values of β

Fig. 3

Variations of skin friction with λ for different values of M

It is clear from Figs. 2−4 that dual solutions exist beyond the critical value (turning point) λc . Large imposition of suction is required so that dual solutions are possible for the flow with large magnetic parameters. In reality, between these two solutions, only one solution is stable. The first solution is assumed to be physically stable because its solution is the continuation of the case of injection (s < 0). The second has negative values of the skin friction coefficient. This solution shows the occurrence of flow separation and reverse, which cause the difficulties in the numerical computation. Merkin[23] , Weidman et al.[24] , Paullet and Weidman[25] , and Harris et al.[26] have presented the mathematical proof of the conjecture of dual numerical

6

N. S. AKBAR, Z. H. KHAN, R. U. HAQ, and S. NADEEM

solutions. They performed stability analyses and revealed that the solutions along the upper branch (first solution) are linearly stable while those on the lower branch (second solutions) are linearly unstable. The effects of the magnetic field parameter M and the mass transfer parameter s are presented in Figs. 3 and 4, respectively, while the rest parameters are fixed. It is observed from Fig. 3 that in the present region for the values of M , we can have dual solutions of the skin friction coefficient against λ. Figure 3 shows that in both cases, lower and upper branches give the same increasing behavior for higher values of M . It is observed from Fig. 4 that both curves have the same increasing behavior for higher values of s. Finally, we show that for the stretching/shrinking parameter λ, dual solutions exist.

Fig. 4

Variations of skin friction with λ for different values of s

Fig. 5

Velocity distribution for different values of α

The variations of the velocity profile f  (η) with various fluid flow parameters are plotted in Figs. 5−7. The effects of the Prandtl fluid parameter α and the elastic parameter β are plotted in Figs. 5−6. The dual velocity profiles show that the velocity decreases with the increases in α and β while conversely increases for the second solution. It is also noted that the boundary layer is thicker for the second solution in comparison with that for the first solution. The effects of the suction/injection parameter s on the dual velocity profiles are shown in Fig. 7. As the mass transfer parameter s increases in the magnitude, both the profiles increase. For large values of the far field boundary η∞ , the second solution velocity profiles decrease with the increase in s.

Fig. 6

Velocity distribution for different values of β

Fig. 7

Velocity distribution for different values of s

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet

5

7

Conclusions

We theoretically study the boundary layer stagnation-point flow of the Prandtl fluid model towards a shrinking sheet in the presence of a magnetic field. We discuss the effects of the suction/injection parameter s, the Prandtl parameter α, the elastic parameter β, and the magnetic parameter M on the fluid flow. The key conclusions of this study are as follows. (i) We show that dual solutions exist for the proposed Prandtl fluid flow model. (ii) The skin friction increases with the increases in the suction/injection parameter s and the magnetic parameter M , whereas with the increase in the elastic parameter β, the skin friction decreases. (iii) It has been observed that the boundary layer thickness for the second solution is thicker than that for the first solution. (iv) For the increasing suction/injection parameter s, the boundary layer thickness for both solutions decreases.

References [1] Hiemenz, K. Die Grenzschicht an einem in den gleichf¨ ormigen Fl¨ uessigkeitsstrom eingetauchten geraden Kreiszylinder. Dinglers Polytechnisches Journal, 321, 321–410 (1911) [2] Sakiadis, B. C. Boundary layer behavior on continuous solid flat surfaces. Journal of AICHE, 7, 26–28 (1961) [3] Crane, L. Flow past a stretching plate. Zeitschrift f¨ ur Angewandte Mathematik und Physik (ZAMP ), 21, 645–647 (1970) [4] Chiam, T. C. Stagnation-point flow towards a stretching plate. Journal of the Physical Society of Japan, 63, 2443–2444 (1994) [5] Wang, C. Y. Free convection on a vertical stretching surface. Journal of Applied Mathematics and Mechanics (ZAMM ), 69, 418–420 (1989) [6] Nazar, R., Amin, N., Filip, D., and Pop, I. Unsteady boundary layer flow in the region of the stagnation-point on a stretching sheet. International Journal of Engineering Science, 42, 1241– 1253 (2004) [7] Ishak, A., Nazar, R., and Pop, I. Magnetohydrodynamic stagnation point flow towards a stretching vertical sheet. Magnetohydrodynamic, 42, 17–30 (2006) [8] Sadeghy, K., Hajibeygib, H., and Taghavia, S. M. Stagnation-point flow of upper-convected Maxwell fluids. International Journal of Non-Linear Mechanics, 41, 1242–1247 (2006) [9] Attia, H. A. Axisymmetric stagnation point flow towards a stretching surface in the presence of a uniform magnetic field with heat generation. Tamkang Journal of Science and Engineering, 1, 11–16 (2007) [10] Kumari, M. and Nath, G. Steady mixed convection stagnation-point flow of upper convected Maxwell fluids with magnetic field. International Journal of Non-Linear Mechanics, 44, 1048– 1055 (2009) [11] Klemp, J. B. and Acrivos, A. A moving-wall boundary layer with reverse flow. Journal of Fluid Mechanics, 76, 363–381 (1976) [12] Riley, N. and Weidman, P. D. Multiple solutions of the Falkner-Skan equation for flow past a stretching boundary. SIAM Journal on Applied Mathematics, 49, 1350–1358 (1989) [13] Ingham, D. B. Singular and non-unique solutions of the boundary-layer equations for the flow due to free convection near a continuously moving vertical plate. Zeitschrift f¨ ur Angewandte Mathematik und Physik (ZAMP ), 37, 559–572 (1986) [14] Ridha, A. and Curie, M. Aiding flows non-unique similarity solutions of mixed-convection boundary-layer equations. Zeitschrift f¨ ur Angewandte Mathematik und Physik (ZAMP ), 47, 341– 352 (1996) [15] Mahapatra, T. R., Nandy, S. K., Vajravelu, K., and van Gorder, R. A. Stability analysis of the dual solutions for stagnation-point flow over a non-linearly stretching surface. Meccanica, 47, 1623–1632 (2012)

8

N. S. AKBAR, Z. H. KHAN, R. U. HAQ, and S. NADEEM

[16] Makinde, O. D., Khan, W. A., and Khan, Z. H. Buoyancy effects on MHD stagnation point flow and heat transfer of a nanofluid past a convectively heated stretching/shrinking sheet. International Journal of Heat and Mass Transfer, 62, 526–533 (2013) [17] Akbar, N. S., Nadeem, S., and Lee, C. Peristaltic flow of a Prandtl fluid model in an asymmetric channel. International Journal of Physical Sciences, 7, 687–695 (2012) [18] Haverkort, J. W. and Peeters, T. W. J. Magnetohydrodynamics of insulating spheres. Magnetohydrodynamics, 45, 111–126 (2009) [19] Haverkort, J. W. and Peeters, T. W. J. Magnetohydrodynamic effects on insulating bubbles and inclusions in the continuous casting of steel. Metallurgical and Materials Transactions B, 41, 1240–1246 (2010) [20] Lok, Y. Y., Amin, N., and Pop, I. Non-orthogonal stagnation point towards a stretching sheet. International Journal of Non-Linear Mechanics, 41, 622–627 (2006) [21] Wang, C. Y. Stagnation flow towards a shrinking sheet. International Journal of Non-Linear Mechanics, 43(5), 377–382 (2008) [22] Mahapatra, T. R. and Nandy, S. K. Stability of dual solutions in stagnation-point flow and heat transfer over a porous shrinking sheet with thermal radiation. Meccanica, 48, 23–32 (2013) [23] Merkin, J. H. On dual solutions occurring in mixed convection in a porous medium. Journal of Engineering Mathematics, 20, 171–179 (1985) [24] Weidman, P. D., Kubitschek, D. G., and Davis, A. M. J. The effect of transpiration on selfsimilar boundary layer flow over moving surfaces. International Journal of Engineering Science, 44, 730–737 (2006) [25] Paullet, J. and Weidman, P. Analysis of stagnation point flow towards a stretching sheet. International Journal of Non-Linear Mechanics, 42(9), 1084–1091 (2007) [26] Harris, S. D., Ingham, D. B., and Pop, I. Mixed convection boundary-layer flow near the stagnation point on a vertical surface in a porous medium: Brinkman model with slip. Transport in Porous Media, 77, 267–285 (2009)