Journal of the Egyptian Mathematical Society (2014) xxx, xxx–xxx Egyptian Mathematical Society Journal of the Egyptian Mathematical Society www.etms-eg.org www.elsevier.com/locate/joems ORIGINAL ARTICLE Steady mixed convection stagnation point flow of MHD Oldroyd-B fluid over a stretching sheet M. Sajid a, B. Ahmed b, Z. Abbas c,* a Theoretical Physics Division, PINSTECH, P.O. Nilore, Islamabad 44000, Pakistan b Department of Basic Sciences, Riphah International University, Islamabad 44000, Pakistan c Department of Mathematics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan Received 3 June 2013; revised 30 April 2014; accepted 29 May 2014 KEYWORDS Abstract This study deals with the steady mixed convection stagnation point flow of an incom- Mixed convection; pressible Oldroyd-B fluid over the stretching sheet in the presence of a constant applied magnetic Stagnation point flow; field. It is assumed that the surface temperature varies linearly with the distance from the stagnation Oldroyd-B fluid; point. A coupled system of non-linear differential equations is developed by employing the similar- Magnetic field; ity transformations. To analyze the behavior of the velocity, temperature, skin friction coefficient Finite difference method and rate of heat transfer through the wall, a numerical solution is developed using finite difference scheme. The obtained results are used to discuss the influence of pertinent parameters of interest. 2010 MATHEMATICS SUBJECT CLASSIFICATION: 65N06; 76D05; 76R99 ª 2014 Production and hosting by Elsevier B.V. on behalf of Egyptian Mathematical Society. 1. Introduction process, manufacturing of sheets, coating and paper produc- tion. The initial work was investigated by Sakiadis [1] for The analysis regarding the fluid flow and heat transfer over a two-dimensional boundary layer flow when the plate is moving stretching sheet is a subject of interest for many researchers with constant velocity. The Sakiadis’s problem for heat transfer working in the area of the two-dimensional boundary layer analysis was studied by Erickson et al. [2]. In the above mention flows. The reasons lie in the fact that such kinds of investiga- studies the velocity of the sheet is assumed to be constant. This tions have applications in the manufacturing industry. assumption of constant velocity is adequate when we are inter- For example, in many manufacturing processes such as glass ested in the analysis of continuous extrusion of polymer sheets. fiber production, hot rolling, continuous casting, extrusion Due to the flexibility of polymer materials, a stretching may occur. Crane [3] found a closed form solution for boundary layer * Corresponding author. Tel.: +92 62 9255480. flow by imposing the condition of stretching wall. He assumed a E-mail address: [email protected] (Z. Abbas). linear variation of stretching velocity with respect to the distance Peer review under responsibility of Egyptian Mathematical Society. from the origin. The Crane’s problem for viscous fluid was extended in different directions in the literatures [4–9],etc. In above mentioned Refs. [3–9], the fluid flow is induced due to the motion of the surface. But in many real situations Production and hosting by Elsevier the thermal buoyancy is also play part for the occurrence of http://dx.doi.org/10.1016/j.joems.2014.05.013 1110-256X ª 2014 Production and hosting by Elsevier B.V. on behalf of Egyptian Mathematical Society. Please cite this article in press as: M. Sajid et al., Steady mixed convection stagnation point flow of MHD Oldroyd-B fluid over a stretching sheet, Journal of the Egyptian Mathematical Society (2014), http://dx.doi.org/10.1016/j.joems.2014.05.013 2 M. Sajid et al. @u @v fluid flow. Therefore, the flow behavior is analyzed under the þ ¼ 0; ð1Þ influence of both the mechanism; the motion of solid surface @x @y and buoyancy. These buoyancy forces occurs due to the heat- @u @u @2u @2u @2u @2u ing and cooling of stretching sheets resulting in producing u þ v ¼K u2 þ 2uv þv2 þa2x þm @x @y 1 @x2 @x@y @y2 @y2 changes to both flow and temperature fields as discussed by Chen and Strobel [10]. The literature survey indicates that @3v @3u @u @2u @u @2v @u þmK u þ v À À þ rB2 ax À u ÀK v work in this direction is carried out by many workers in the 2 2 @y3 @x @y2 @y @y2 0 1 @y @x@ y field [11–14] and reference therein. On the other hand, the con- @T @T @u þgb ðT ÀT ÞþK þ v À ðTÀ T Þ ; ð2Þ stitutive relationships for rate type fluids are implicit and elim- 1 1 @x @y @x 1 ination of stress components from the equation of motion is not straightforward. This fact is the major cause for the lack @T @T k @2T l @u 2 u v ; 3 of literature on the two-dimensional flow of rate type fluids. þ ¼ 2 þ ð Þ @x @y qcp @y qcp @y The simplest class of rate type fluids is the Maxwell fluid and one can find a number of articles regarding the stretching flow where u and v are the velocity components in the x and y direc- of Maxwell fluid [15–20] and reference there in. However, for tions respectively, m is the kinematics viscosity of fluid, p is the an Oldroyd-B fluid there are only two studies regarding the pressure, q is the density of the fluid, K 1 is the relaxation time, flow over a stretching sheet [21,22]. K2 is the retardation time, g is the gravitational acceleration, b On the other hand, an extensive literature can be found on is the thermal expansion coefficient, cp is the specific heat, k is the mixed convection flow of Non-Newtonian fluids over the thermal diffusivity and T is the temperature. continuously moving surfaces. Mixed convection flow of a The appropriate boundary conditions applicable to the non-Newtonian micropolar fluid was discussed by Takhar present flow problem are: et al. [23]. Mushtaq et al. [24] considered the steady mixed con- u ¼ uwðxÞ¼cx; v ¼ 0; T ¼ TwðxÞ¼T1 þ bx at y ¼ 0; vection flow of a second grade fluid by considering the case of @u variable surface temperature. The steady mixed Convection u ¼ ueðxÞ¼ax; ! 0; T ¼ T1 as y !1: ð4Þ flow of a micropolar fluid near the stagnation point on a @y vertical surface is studied by Lok et al. [25]. In 2008 , Hayat in which both a and c have the dimension of (time)À1 and b is a et al. [26] investigated the series solution for the mixed convec- positive constant. Here, the forth condition in Eq. (4) is the tion flow of a micropolar fluid over a non-linear stretching augmented condition discussed by Grag and Rajagopal [29] . sheet using homotopy analysis method. Recently, Hsiao [27] Introducing the standard similarity transformations for a studied the heat, mass transfer and mixed convection for stretching flow MHD flow of a viscoelastic fluid past a continuously moving rffiffiffi pffiffiffiffiffi T À T c surface with Ohmic dissipation numerically. Very recently, u ¼ cxf0ðgÞ; v ¼ cmfðgÞ; hðgÞ¼ 1 ; g ¼ y: the mixed convection in the stagnation-point flow of a Maxwell Tw À T1 m fluid toward a vertical stretching sheet has been investigated by ð5Þ Abbas et al. [28]. They have discussed the results both analyti- Governing flow problem takes the form cally using homotopy analysis method (HAM) and numerically 000 02 00 a2 2 a 0 00 0 using finite difference method. To best of our knowledge, no f À f þ ff þ 2 þ M ð À f À k1ff ÞþArðÞh À k1fh c c ÀÁ ð6Þ such attention has been given for the mixed convection flow 0 00 2 000 002 iv þk1ðÞþ2ff f À f f k2 f À ff ¼ 0; of an Oldroyd-B fluid. Our aim was to discuss the mixed con- 00 0 0 002 vection in the stagnation-point flow of an Oldroyd-B fluid over h þ Prðfh À hf ÞþPrEcf ¼ 0; ð7Þ a stretching sheet. As a first step the boundary layer equations f ¼ 0; f0 ¼ 1; h ¼ 1atg ¼ 0; under these assumptions have been developed and then a a f0 ¼ ; f00 ¼ 0; h ¼ 0asg !1: ð8Þ numerical solution by employing a finite difference method is c presented for the transformed nonlinear coupled ordinary Here prime denotes derivative with respect to g and k1 = K1c differential equations. The effects of various involving physical and k2 = K2c are the dimensionless material parameters of the parameters on the flow and temperature distributions are 2 2 fluid and Ar = Gr/Re is the Archimedes number, Gr = gbb/m discussed through graphs and tables. is the Grashof number, Re = c/m is the Reynolds number and 2 2 M ¼ rB0=cq is the Hartmann number. 2. Formulation of the problem 3. Results and discussion Consider an incompressible, two-dimensional Oldroyd-B fluid in the region of a stagnation point over a semi-infinite, imper- The numerical procedure explained in detail by Sajid et al. [22] meable stretching sheet at y = 0. The assumed stretching is adopted for the solution of the transformed Eqs. (6)–(8). The velocity and wall temperature of the sheet is uw(x)=cx, with effects of the involving parameters for example, the dimension- c > o and Tw(x)=bx, with b > 0, respectively. Far away less relaxation/retardation times (k1,k2), the Hartmann number from the plate the external flow velocity is ue(x)=ax, with M, the ratio of the external flow rate to the stretching rate a/c , a > 0 and we assume that body attains a uniform temperature the Archimedes number Ar , the Prandtl number Pr and the 0 T1 outside the boundary layer.