Numerical Solution for Variable Accelerated Flow Subject to Slip Effect

In this paper, we examine the unsteady magneto hydrodynamic (MHD) flow generated by a disc that is making non-coaxial rotations with a third grade fluid at infinity and moving with a variable acceleration. The fluid is assumed to satisfy slip boundary condition on the disc. The governing equations are three dimensional and highly non-linear in nature. The assumed slip boundary condition is non-linear as well. The governing equations are transformed to a nonlinear boundary value problem which is solved numerically. Comparison of this generalized problem with uniformly accelerated disk satisfying no slip condition is made. Variations of the characterizing dimensionless parameters such as slip parameter λ, acceleration parameter c, unsteady parameter τ, third grade parameter β, suction parameter S, and magnetic parameter N on the flow field are discussed and analyzed graphically.


Introduction
Considerable interest has been shown in the literature to investigate the steady/unsteady flows caused by eccentric rotations of disk and viscous/rheological fluid at infinity incorporating various parameters of physical interest.Erdogan [1] analyzed the transient flow resulting from eccentric rotations of disc and fluid at infinity.Some MHD boundary layer flows are mentioned in references [2]- [7] that helped to improve our understanding of astrophysical, geophysical and engineering problems.Hayat et al. [8] extended the work of Ref. [1] for porous disk and applied magnetic field.The transient flow caused by eccentric rotations of a porous disk executing oscillations and a fluid at infinity was ex-amined by Erdogan [9].Non-Newtonian fluid flows [10]- [16] are also very important in real world applications.Siddiqui et al. [17] and Hayat et al. [18] generalized the flow analysis of ref. [8] and [9] for non-Newtonian second and third grade MHD fluids respectively.Asghar et al. [19] [20] discussed the flow generated by eccentric rotations of constantly accelerated disk and non-Newtonian fluid at infinity using no slip boundary condition.The constant acceleration corresponds to the linear velocity; however the nonlinear velocity; that is, variable acceleration will be both interesting and realistic.
Although no slip boundary condition is well established in fluid mechanics; there are situations where this condition is not adequate in that the flow may exhibits slip effect [21] [22] [23] [24] [25].It has been experimentally verified that the slip effects become important in nanochannels and micro channels (see [26]).Concerning non-Newtonian fluids, no slip condition is not adequate while considering polymer melts.What we have chosen here is to study the unsteady flow generated by two independent mechanisms of non-coaxial rotation and the non-linear velocity (variable acceleration) of the disc with the slip condition.
Crank-Nicolson scheme is used for the numerical results of the nonlinear governing differential equation and non-linear boundary conditions.The influence of parameters of importance is analyzed and interpreted physically putting emphasis on acceleration parameter and slip parameter.

Problem Statement
Consider the unsteady flow of an incompressible third grade fluid filling the semi-infinite space z > 0 above a disk at z = 0.The disk is taken to be porous and non-conducting.Both cases of suction and injection/blowing are to be analyzed.
The flow is generated by eccentric rotations of variable accelerated disk and third grade fluid at infinity.The disk and the fluid at infinity initially rotate about z'-axis with angular velocity Ω (see Figure 1).At t = 0, the disk suddenly starts rotating about z-axis with same angular velocity Ω and moving linearly with variable acceleration along x-axis.However, the fluid at t = 0 continues to rotate about the z'-axis.The axes of rotation of the disk and fluid at infinity are taken in the plane x = 0.The distance between the axes of rotations is denoted by l.The fluid is electrically conducting and a constant magnetic field B 0 is applied in the transverse direction.We consider a slip between the velocity of the fluid at the surface of the disk and the velocity of the disk.The relative velocity between the fluid and the wall is assumed to be proportional to the shear rate at the wall.
The Cauchy stress tensor T (Fosdick and Rajagopal [27]) for a thermodynamically compatible third grade fluid is given by ( ) in which p is the hydrostatic pressure, µ is the dynamic viscosity of the fluid.The equations governing the present flow are where ρ is the density and σ is the electrical conductivity of the fluid.
The velocity field , , u v w = V for our problem is given by satisfying the equation of continuity (2) identically.Here W 0 > 0 denotes the suction and W 0 < 0 the injection velocity.
The initial and boundary conditions of the problem are given by , as , for all , , for 0, and at 0, where λ 1 is a slip parameter, c 0 is a constant having dimension L/T 3 and the shear stresses are Open Journal of Fluid Dynamics Substituting (4) into (3), gives ( ) where the modified pressure P is expressed as ( ) Equation ( 12) shows that P is not a function of z.We can eliminate P from Equations ( 10)-( 11) by differentiating and then integrating with respect to z, and combining the resulting equations to obtain the following equation ( ) where Equation ( 15) together with the Conditions ( 16) and ( 17), in dimensionless variables, can be expressed as ( ) where

Numerical Solution
We obtain a numerical solution of the problem given by the governing Equation ( 18) together with Conditions ( 19)- (21).The discretized form of Equation ( 18) is presented using a finite difference scheme in ref. [19].However, the boundary Condition ( 19) is discretized as where ) ( ) .
To evaluate 0, 1 j F + , we firstly take in the system of algebraic equations and the solution of the system is sought which results in known values of , 1 ; 1, 2, 3, , 1 . Secondly, we update 0, 1 j F + by using iterative method as follows: where and this iterative procedure is continued until and by using iterative method as described above with 0 0,0 0 F = as initial guess. For where , , , , , , , , , , , ( ) .
In matrix form, the above set of M-1equationscan be finally written as .

Discussion
Recalling that we have introduced the new features of the slip boundary condition and the variable acceleration of the disc in non-coaxial rotation; we will lay emphasis to investigate the variation of slip parameter λ and the acceleration parameter c on the velocity profiles.However, for completeness we will show the variations of other parameters as well to look for the joint effects of already discussed and the new phenomena.Firstly, in order to ensure reliability of our numerical scheme; the numerical and the analytical results [20] are compared in Table 1, showing a good agreement.
Table 1.τ = 1, α = 0.1, β = 0, c = 1, S = 0.5, N = 0.5.with c.Furthermore, it is evident from Figure 4(b) that the velocity profile g/Ωl increases with increasing c close to the disk while it decreases away from the disk.The momentum boundary layer thickness for the velocity profile g/Ωl, decreases with increasing c.

Time Evolution of the Velocity Field
The evolution of the velocity profiles f/Ωl and g/Ωl for different values of time .This leads us to believe that; as the time evolves, the velocity profile f/Ωl increases close to the disk and decreases nominally farther from the disk.We note from Figure 5(a) that thickness of the boundary layer increases as the time τ progresses.It is evident from Figure 5(b) that the velocity profile g/Ωl and thickness of the boundary layer are increasing with time τ .otherwise the behavior is almost the same as has been seen in Figure 6.

Variation of Other Physical Parameters
The effects of injection/suction can be seen from .This figure also shows that increasing porosity   parameter S the velocity profile f/Ωl increases whereas the velocity profile g/Ωl decreases.However the boundary layer thickness reduces for f/Ωl and g/Ωl.The Open Journal of Fluid Dynamics effects of magnetic field parameter N are identical to that of suction (see Figure 9).
To see the effects of rotation Ω; we draw the Figure 10            The effects of these quantities on g/Ωl are exactly the opposite.The effects of suction and injuction/blowing on the boundary layer thickness remain the same irrespective of third grade fluid and the slip condition.That is, the suction causes thinning and blowing causes thickening of the boundary layer thickness.We remember that the magnetic field is applied parallel to the axis of rotation.This causes a decrease in the thickness of boundary layer for f/Ωl and an increase for the velocity profile g/Ωl.The comparison of the velocity profile for variable acceleration with constant acceleration of the disc discussed earlier [2] reveals important observation.For time τ < 1, the velocity profiles for constant accelerated flow are greater than that of variable accelerated flow.However, for the time τ > 1, the velocity profiles for variable accelerated flow are much larger than constant accelerated flow.

Figure 2 and
Figure 2 and Figure 3 show the effects of slip parameter λ on the flow pattern when 0.2 α = , 1.0 β = , 0.25 c = , 0 S N = = for 0.5 τ = and 1.5 τ = .We note that increasing in λ enhances the velocity profile f/Ωl; whereas the velocity profile f/Ωl increases close to the disk and decreases away from the disk.It is further observed that the slip reduces the boundary layer thickness for the two velocity profiles.Thus the slip condition helps to control the boundary layer for rotating disc.

Figure 4 (Figure 2 .
Figure 4(a) reveals that the velocity profile f/Ωl increases with small increase in acceleration parameter c keeping the values of the parameters 0.2 α = ,
parameter τ in the case of rigid plate and non-conducting third grade fluid β , are shown in Figure5.Here we take 0

Figure 8 .
The fluid velocities f/Ωl and g/Ωl are found to increase and decrease respectively at the surface of the disk; for slip parameter 0
Figure 10 and Figure 11.Remembering that this problem is essentially a three dimensional; three dimensional velocity profiles are plotted in Figures 12-16.These streamlines are drawn by taking into consideration various physical parameters mentioned in Figure captions.
In the present work, the unsteady three-dimensional MHD flow caused by a porous disk rotating non-coaxially with fluid at infinity is investigated numerically.Both the disk and the fluid at infinity are rotating with the same angular velocity.Particular emphasis has been given to non-Newtonian third grade fluid, slip boundary condition and variable acceleration of the disk.We observe that H. A. Ashi DOI: 10.4236/ojfd.2017.74033497 Open Journal of Fluid Dynamics