Page 1
Study of the effects of rotation and magnetic field on thestagnation flow of a second grade fluid towards a
stretching/shrinking sheet using HPM and BPES protocols
Abstract
In the present investigation we have presented the effects of rotation and
magnetic field on the stagnation point flow of a second grade fluid towards a
porous stretching /shrinking sheet. The governing equations of second grade
fluid are simplified by using boundary layer approach and similarity
transformation. The reduced nonlinear coupled differential equations are
solved using the Boubaker Polynomials Expansion Scheme BPES and the
Homotopy Analysis Method HAM. The results of various physical
parameters are discussed through graphs. The convergence of the solution is
discussed by h-curves, homotopy pade approximation and residual error.
Keywords: Stagnation flow; stretching/shrinking sheet; rotating frame; HAM solutions; BPES solutions.
1
Page 2
1. Introduction
In a number of applications, rotation plays a significant role such as in cosmical fluid
dynamics. Similarly a great deal of metrology depends upon the dynamics of the
revolving fluid. The large scale and moderate motion of the atmosphere are greatly
affected by the vorticity of the earth’s rotation. In the case of infinite fluid rotating as a
rigid body about an axis, the amount of energy possessed by the fluid is infinite and it is
of great interest to know how small disturbances propagate in such a fluid. Greenspan
and Howard [1] have initiated the study of the dynamics of spin up of an incompressible,
homogenous Newtonian rotating fluid. They have presented a detailed mathematical and
physical analysis of the transient process by which the fluid adjusts to a small change in
the rotation rate of its boundary. It has been demonstrated that the Ekman boundary layer
is established on the horizontal boundary surfaces and it is primarily responsible for the
adjustment process. Furthermore, the Ekman layer produces a secondary interior
circulation throughout the fluid which transports angular momentum. After initiation of
Greenspan and Howard [1], a large number of studies have been presented by various
researchers keeping different fluids and its geometries for rotating flows [2-6].
The analysis of the effects of rotation and magnetic field in the fluid flows has been an
active area of research because of its geophysical and technological importance. It is well
known that a number of astronomical bodies (e.g. the Sun, Earth, Jupiter, magnetic stars,
Pulsars) possess fluid interiors and magnetic fields. Changes in rotation rate of such are
discussed elsewhere [7, 8].
Boundary layer behavior over a stretching surface is another area of research which has
attracted the attention of many researchers due to its important applications in
engineering. Such applications includes aerodynamics extrusion of plastic sheets, the
boundary layer along a liquid film condensation process, the cooling process of metallic
plate in a cooling bath and in the glass and polymer industries. After the initiation of
2
Page 3
Sakiadis [9] several authors consider the boundary layer over a stretching surface [10-
15].
In almost all the fields of engineering and scientific interest, stagnation point flow is an
important phenomenon. These flows may be characterized inviscid or viscous, steady or
unsteady, two dimensional or three dimensional, symmetric or asymmetric, normal or
oblique, homogenous or two fluids and forward or reverse. Some important studies to the
topic are described in refs. [16-22].
Motivated from the above analysis, the aim of this paper is to investigate the effects of
rotation and magnetic field on the stagnation point flow of a second grade fluid towards a
similarity transformation and then the reduced problem have been solved using the
Boubaker Polynomials Expansion Scheme BPES and the Homotopy Analysis Method
HAM. Homotopy analysis method is the successful tool for solving nonlinear ordinary
and partial differential equation and some recent related developments on the topic are
defined in [23-27]. The Boubaker Polynomials Expansion Scheme BPES is also aprotocol
which has been significantly tested in severel applied physics fields [28-38]. The solution
for both stretching and shrinking phenomena are presented. The results of the necessary
physical parameters are discussed through graphs. The convergence of the solution have
been discussed through h-curves, homotopy pade approximation and residual error.
2. Mathematical formulation
We consider an incompressible, stagnation point flow of a second grade fluid
towards a porous stretching/shrinking sheet in a rotating frame in the presence of the
magnetic field. We are choosing the cartesian coordinate system in which (u; v;w) are the
velocity components along x; y and z-axis with Ω being the angular velocity of the
rotating fluid in the z-direction. The potential stagnation flow at infinity is taken by u =
ax, w = -az, v=0 (in which a is the strength of the stagnation flow) where as the velocities
on the shrinking(stretching) rate are considered u = bx, w =-W, v=0, here b > 0 is the
stretching rate (shrinking if b < 0). A constant magnetic field B0 is imposed along z-axis
vertical to the sheet. The governing boundary layer equations for the flow problem are:
3
Page 4
We introduce the similarity variables and non-dimensional variables as follows:
Making use of Eq. (4); the incompressibility condition is automatically satisfied and the momentum equations take the form:
where , , and prime denotes the derivative with respect to .
The boundary conditions for the problem under consideration are:
4
Page 5
where
3. HAM Solution
To seek HAM solution of Eqs. (5) and (6) we select:
as the initial approximations of f, h and:
as auxiliary linear operators which satisfy:
where Ci (i = 1..5) are arbitrary constants. If p [0, 1] is an embedding parameter and
(i = 1.. 2) are non-zero auxiliary parameters then the zeroth order and mth order
deformation problems are:
Zeroth order deformation problems:
5
Page 6
where:
mth order deformation problems:
6
Page 7
where :
The symbolic software MATHEMATICA isused to get the solutions of Eqs. (21) and
(22) up to first few order of approximations. It is found that the solution for f, h are given
by:
in which the coefficients and of and can be found by using given
boundary conditions and by the initial guess approximations in Eqs. (8) and (9) and
numerical data is presented through graphs.
4. BPES Solution
To seek Boubaker polynomials expansion scheme BPES [28-38] solution of Eqs. (5) and
(6) we set:
(29)
7
Page 8
and:
(30)
where are the 4k-order Boubaker polynomials, are minimal positive roots,
is a prefixed integer, and are unknown pondering real coefficients.
The main advantage of these formulations is the fact of verifying the four boundary
conditions in Eq. 1, in advance to problem resolution.
Due to the properties of the Boubaker polynomials [29-32] , and since are roots,
the following conditions stand, for i.e. :
(31)
By introducing expressions (29) and (30) in Eq. (5-8), and by majoring and integrating
along the interval , and are confined, through the coefficients and
, to be weak solutions of the system (32):
(32)
which gives:
(33)
with
8
Page 9
The set of solutions and is the one which minimizes the global
Minimum Square function :
(34)
The correspondent solutions are represented along with the HPM solution (27-28).
5. Convergence of the analytic solution
The analytic series solutions given in Eqs. (27) and (28) of the considered problems are
found by homotopy analysis method. These expressions contain non-zero auxiliary
parameters (i = 1.. 2) which can adjust and control the convergence of the solutions.
To ensure the convergence of the solutions in the admissible range of the values of the
auxiliary parameters one can draw h curves. For the present cases, 20th order -curves
are plotted in Figs. (1) and (2). It is seen from these figures that the ranges for the values
of are -0.7 ≤ ≤ -0.2 and -0.95 ≤ ≤ -0.25. The residual error of the problem is
calculated and shown in Fig (3).
6. Results and discussion
This section deals with the influence of stretching/shrinking parameter
α, magnetic field parameter M, rotation parameter and suction parameter s on and h.
The velocities and h against for various values of α are presented in Figs (4) and (5).
It is observed that magnitude of and h increases with the increase in α. However, the
layer thickness decreases. The variation of M on and h are shown in Figs (6) and (7). It
is observed that the magnitude of increases with the increase in M whereas h decreases
with the increase in M however boundary layer thickness decreases for both the velocity
cases. The magnitude of increases while the magnitude of h decreases with the
9
Page 10
increases in (see Figs (8) and (9)). Figures (10) and (11) show that the suction causes
boundary layer reduction.
It is observed from Figs (12) and (13) that increases with the increase in β whereas h
increases with increase in β up to β = 0.09 then it gives the reverse behavior.
References
[1] H. P. Greenspan and L. N. Howard, On a time dependent motion of a rotating
fluid, J. Fluid Mech. 17; 385(1963):
[2] B. S. Mazumder, An exact solution of oscillatory couette flow in a rotating
system,ASME J. Appl. Mech. 58; 1104(1991):
[3] R. Ganapathy, A note on oscillatory couette áow in a rotating system, ASME J.
Appl. Mech. 61, 208 (1994)
[4] P. Puri, Rotating áow of an elastico-viscous áuid on an oscillating plate,
Z.Angew. Math. Mech. (ZAMM). 54; 743(1974).
[5] P. Puri and P. K. Kulshrestha, Rotating áow of non-Newtonion áuids, Appl. Anal.
4; 131(1974).
[6] K. Vajravelu and L. Debnath, A study of non-linear convective áows in rotating
wavy channels, Z.Angew. Math. Mech. (ZAMM). 64; 303(1984).
[7] S. Chandrasekhar, The instability of a layer heated below and subjected to a
coriolis force, Proc. Roy. Soc. Lond. A 217; 306(1953).
[8] S. Chandrasekhar, Hydrodynamic and hydromagnetic stability, Oxford University
Press, 1961.
[9] B. C. Sakiadis, Boundary layer behaviour on continuous solid surfaces, AIChE J 7
(1961) 26-28.
10
Page 11
[10] H. Xu and S. J. Liao, Dual solutions of boundary layer flow over an upstream
moving plate, Comm. Nonlinear Sci. Num. Sim. 13 (2008) 350-358.
[11] S. J. Liao, An analytic solution of unsteady boundary layer flows caused by an
impulsive stretching plate, Comm. Nonlinear Sci. Num. Sim. 11 (2006) 326-339.
[12] S. J. Liao, A new branch of solutions of boundary layer áows over an
impermeable stretching plate, Int. J. Heat Mass Transfer, 48 (2005) 2529- 2539.
[13] M. S. Abel, Mahantesh, M. Nandeppanavar, Heat transfer in MHD viscoelastic
boundary layer áow over a stretching sheet with non-uniform heat source/sink,
Comm. Non-linear Sci. Num. Sim. 14 (2009) 2120- 2131.
[14] A. Ishak, R. Nazar, I. Pop, MHD boundary-layer flow of a micropolar
áuid past a wedge with constant wall heat flux, Comm. Nonlinear Sci.
Num. Sim. 14 (2009) 109-118.
[15] S. Bose, S. Chakraborty, A boundary layer analysis of electro-magneto-
hydrodynamic forced convective transport over a melting slab, Int. J. Heat Mass
Transfer 51 (2008) 5465-5474.
[16] A. Ishak, R. Nazar, I. Pop, Dual solutions in mixed convection boundary layer
flow of micropolar fluids, Comm. Nonlinear Sci. Num. Sim. 14 (2009) 1324-
1333.
[17] C. Y. Wang, Stagnation flow towards a shrinking sheet, Int. J. Nonlinear. Mech.
43 (2008) 377-382:
[18] C. Y. Wang, Stagnation slip flow and heat transfer on a moving plate, Che. Eng.
Sci. 61 (2006) 7668-7672:
11
Page 12
[19] C. Y. Wang, Off-centered stagnation flow towards a rotating disc, Int. J. Eng.
Sci. 46 (2008) 391-396.
[20] A. Ishak, Roslinda Nazar, Norihan M. Arifin, Ioan Pop, Dual solutions in mixed
convection flow near a stagnation point on a vertical porous plate, Int. J. Thermal
Sci., 47(2008) 417-422:
[21] A. Ishaka, Roslinda Nazar, Ioan Pop, Stagnation áow of a micropolar áuid
towards a vertical permeable surface,Int.Comm. Heat Mass Transfer 35 (2008)
276-281.
[22] C. Y. Wang, Similarity stagnation point solutions of the Navier-Stokes
equations review and extension, Eur. J.Mech. B/Fluids 27 (2008) 678-683.
[23] S. Abbasbandy, Approximate solution for the nonlinear model of di§usion and
reaction in porous catalysts by means of the homotopy analysis method,
Chem. Eng. J. 136 (2008)144-150.
[24] R. Ellahi, T. Hayat, T. Javed, S. Asghar, On the analytic solution of
nonlinear flow problem involving Oldroyd 8-constant fluid, Math. Comp.
Modelling 48 (2008) 1191-1200.
[25] S. Abbasbandy, Soliton solutions for the Fitzhugh-Nagumo equation with the
homotopy analysis method, Appl. Math. Modelling 32 (2008) 2706-2714.
[26] S. Abbasbandy, Y. Tan, S.J. Liao, Newton-homotopy analysis method for
nonlinear equations, Appl. Math. Comp. 188 (2007)1794-1800.
[27] S. Abbasbandy, The application of homotopy analysis method to nonlinear
equations arising in heat transfer, Phy. Lett. A 360 (2006) 109-113.
12
Page 13
[28] Oyodum, O.D., Awojoyogbe, O.B., Dada, M., Magnuson, J., 2009. On the earliest
definition of the Boubaker polynomials. Eur. Phys. J. Appl. Phys. 46, 21201–
21203.
[29] Awojoyogbe, O.B., Boubaker, K., 2008. A solution to Bloch NMR flow equations
for the analysis of homodynamic functions of blood flow system using m-
Boubaker polynomials. Curr. Appl. Phys. 9, 278–283.
[30] Ghanouchi, J., Labiadh, H., Boubaker, K., 2008. An Attempt to solve the heat
transfer equation in a model of pyrolysis spray using 4q-order Boubaker
polynomials. Int. J. Heat Technol. 26, 49–53.
[31] Slama, S., Bouhafs, M., Ben Mahmoud, K.B., Boubaker, A., 2008. Polynomials
solution to heat equation for monitoring A3 point evolution during resistance spot
welding. Int. J. Heat Technol. 26, 141–146.
[32] Slama, S., Bessrour, J., Bouhafs, M., BenMahmoud, K. B., 2009a .Numerical
distribution of temperature as a guide to investigation of melting point maximal
front spatial evolution during resistance spot welding using Boubaker
polynomials.Numer.HeatTransferPartA55,401–408.
[33] Fridjine, S., Amlouk, M., 2009. A new parameter: an ABACUS for optimizing
functional materials using the Boubaker polynomials expansion scheme. Mod.
Phys. Lett. B 23, 2179–2182.
[34] Tabatabaei, S., Zhao, T., Awojoyogbe, O., Moses, F., 2009. Cut-off cooling
velocity profiling inside a keyhole model using the Boubaker polynomials
expansion scheme. Heat Mass Transfer 45, 1247–1251.
13
Page 14
[35] Zhao, T.G., Wang, Y.X., Ben Mahmoud, K.B., 2008. Limit and uniqueness of the
Boubaker–Zhao polynomials imaginary root sequence. Int. J. Math. Comput. 1,
13–16.
[36] Belhadj, A., Onyango, O., Rozibaeva, N., 2009. Boubaker polynomials expansion
scheme-related heat transfer investigation inside keyhole model. J. Thermo- phys.
Heat Transfer 23, 639–640.
[37] Guezmir, N., Ben Nasrallah, T., Boubaker, K., Amlouk, M., Belgacem, S., 2009.
Optical modeling of compound CuInS2 using relative dielectric function approach
and Boubaker polynomials expansion scheme BPES. J. Alloys Compd. 481, 543–
548.
[38] Dubey, B., Zhao, T.G., Jonsson, M., Rahmanov, H. 2010. A solution to the
accelerated-predator-satiety Lotka–Volterra predator–prey problem using
Boubaker polynomial expansion scheme. J. Theor. Biology 264, 154-160.
14