Unsteady Incompressible Flow of a Generalized Oldroyd-B Fluid between Two Oscillating Infinite Parallel Plates in Presence of a Transverse Magnetic Field

In this paper an attempt has been made to study the unsteady incompressible flow of a generalized Oldroyd-B fluid between two oscillating parallel plates in presence of a transverse magnetic field. An exact solution for the velocity field has been obtained by means of Laplace and finite Fourier sine transformations in series form in terms of Mittage-Leffler function. The dependence of the velocity field on fractional as well as material parameters has been illustrated graphically. The velocity fields for the classical Newtonian, generalized Maxwell, generalized second grade and ordinary Oldroyd-B fluids are recovered as limiting cases of the flow considered for the generalized Oldroyd-B fluid.


Introduction
The magneto hydrodynamic flow problem between two parallel plates has shown immense attention during the last several decades.The study has significant applications in the field of hydrodynamical machines and apparatus, magnetic storage devices, computer storage devices, lubrication, crystal growth processes, radial diffusers, MHD pumps, MHD power generators, purification of crude oil, petroleum industries etc. Bandelli et al. [1] discussed start-up flows of second grade fluids in domains with one finite dimension.Fetecau et al. [2] investigated exact solutions for the flow of a generalized Oldroyd-B fluid induced by a constantly accelerating plate between two side walls perpendicular to the plate.Hayat et al. [3] made homotopy analysis of MHD boundary layer flow of an upper-convected Maxwell fluid.Jamil and Khan [4] studied slip effects on fractional viscoelastic fluids.Shen et al. [5] studied the Rayleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative model.Vieru et al. [6] discussed the flow of a generalized Oldroyd-B fluid due to a constantly accelerating plate.Wenchang et al. [7] investigated unsteady flows of a viscoelastic fluid with the fractional Maxwell model between two parallel plates.Vieru et al. [8] studied the unsteady flow of a generalized Oldroyd-B fluid due to an infinite plate subject to a time-dependent shear stress.
In the present paper we consider the flow of a generalized Oldroyd-B fluid between two oscillating infinite parallel plates in presence of transverse magnetic field.We have formulated the expression for the velocity field for the said flow in terms of Mittage-Leffler function.In the constitutive equation of the fluid model, the time derivative of integral order has been replaced by Riemann-Liouville fractional calculus operator.The exact solution for the velocity field is obtained by using the method of integral transformations and the dependence of the said field on the material as well as fractional calculus parameters is illustrated graphically.

Mathematical Formulation and Basic Equation
Let us consider an incompressible generalized Oldroyd-B fluid bounded by two infinite parallel plates as shown in Figure 1.The plates are initially at rest and at 0 t the plates start to oscillate in its plane with the velocity ( ) where V is the fluid velocity.Due to the shear, the fluid is moved gradually.We have taken Cartesian coordinate system.x-and y-coordinates are taken along and perpendicular to the parallel plates respectively.Accordingly, the initial condition is given by ( ) ,0 0 u y = , 0 1 y ≤ ≤ and the boundary conditions are given by ( ) ( ) , cos u d t V t ω = .We take the velocity and stress of the form where ( ) , u y t is the velocity component in the x-direction.The constitutive relationship for the fluid associated with the present problem is given by, In the relation (2), t D α and t D β are Caputo operators defined by We consider a generalized Oldroyd-B fluid between two infinite parallel disks in presence of an imposed magnetic field 0 B that acts in the direction of the positive y-axis.Then in the presence of the body force where "σ " is constant and " ρ " is the density of the fluid.
Eliminating xy S between the Equations ( 2) and ( 4) we have the governing equation we get the governing equation in non-dimensional quantities as Taking finite Fourier sine transformation we get from Equation ( 6) is the finite Fourier sine transformation of ( ) , u y t .
Using the boundary conditions (8) the Equation (9) can be rewritten as Taking Laplace transformation and using ( ) Now in order to avoid the lengthy procedure of residues and contour integrals, we rewrite the Equation (11) into series form as ) Now we have an important Laplace transformation of the nth order derivative of Mittage-Leffler function given by ( ) ( ) ( ) where Taking Laplace Inverse transformation we get from the Equation (12) Taking inverse finite Fourier sine transformation we get the velocity profile from the Equation (18) as ( )

Limiting Cases
then the equation of motion is given by ( subject to the initial and boundary conditions given by the Equations ( 7) and ( 8) respectively.The Equation (17) represents the governing equation of a classical Newtonian fluid and the corresponding velocity field is given by .
then the equation is given by 1 , subject to the initial and boundary conditions given by the Equations ( 7) and ( 8).The Equation ( 19) represents the governing equation of a generalized Maxwell fluid and the corresponding velocity field is given by then the equation of motion is given by subject to the initial and boundary conditions given by the Equations ( 7) and ( 8) respectively.The Equation ( 23) represents the governing equation of an ordinary Oldroyd-B fluid and the corresponding velocity field is given by

Conclusions and Numerical Results
In this paper we have presented the flow of a generalized Oldroyd-B fluid between two oscillating infinite parallel plates.The velocity field has been determined by means of Laplace and finite Fourier sine transformations in series form in terms of Mittage-Leffler function.The dependence of the velocity field on the fractional calculus parameters and material parameters has been illustrated graphically.The solutions for the four limiting cases have been discussed from the solution of the flow problems of a generalized Oldroyd-B fluid. In

Figure 1 .
Figure 1.Geometry of the problem.

Figure 2 2 . 3 .
the velocity is depicted against the distance from the lower plate for different values of the fractional calculus parameter α .As α increases, the fluid velocity increases and there are points of local mini- mum and local maximum in the velocity curves which are oscillatory in nature.Negative velocity can be observed near the upper plate for values of α near zero in Figure The velocity is depicted against the distance from the lower plate for different values of fractional calculus parameter β in Figure As β increases, the fluid velocity decreases, which is opposite to the case in Figure 2 and the points of local minimum and maximum can be observed in the velocity curves.The flow patterns are oscillatory in nature.In Figure 4 the velocity is plotted against the distance from the lower plate for different values of the parameter M. As the parameter M takes increasing values, the fluid velocity increases and the velocity curves are oscillatory in nature.The velocity profile is plotted against the distance from the lower plate for different values of the frequency of oscillation 1 ω of the lower plate in Figure 5.It is evident from the figure that the fluid velocity decreases for higher values of the parameter 1 ω .The initial points of the velocity curves near the lower plate are different in domain of spatial

Figure 4 .
Figure 4.The velocity is depicted against the distance from the lower plate for different values of parameter M . 1 1.2 ω = ,

Figure 5 .
Figure 5.The velocity is depicted against the distance from the lower plate for different values of the parameter 1 ω .

Figure 6 .
Figure 6.The velocity profile is depicted against the distance from the lower plate for different values of the parameter 2 ω .

Figure 7 .
Figure 7.The velocity profile is depicted against the distance from the lower plate for different values of the fractional calculus parameter α and for equal values of 1 ω and 2 ω , 1

Y for different values of frequency of oscillation 1 ωFigure 5 2 ω 2 ω
Figure5near the lower plate for values for higher frequency of oscillation of the lower plate.In Figure6as the frequency of the oscillation 2 ω of the upper plate changes, the terminal points of the velocity curves near the upper plate differ.The fluid velocity decreases with the increase of the frequency of oscillation 2 ω of the upper plate near that one.It can be noticed that there are points of local minimum and maximum for velocity curves for all the three cases.Negative velocity can be observed near the upper plate in Figure6for higher values of the frequency of oscillation of the upper plate.In Figure7the velocity profile is depicted against the distance from the lower plate for different values of the parameter α in which the frequencies of oscillations of the plates are equal i.e. 1 , ,