Unsteady Incompressible Viscoelastic Flow of a Generalised Maxwell Fluid between Two Rotating Infinite Parallel Coaxial Circular Disks

The unsteady incompressible viscous flow of a Generalised Maxwell fluid between two coaxial rotating infinite parallel circular disks is studied by using the method of integral transforms. The motion of the fluid is created by the rotation of the upper and lower circular disks with different angular velocities. A fractional calculus approach is utilized to determine the velocity profile in series form in terms of Mittag-Leffler function. The influence of the fractional as well as the material parameters on the velocity field is illustrated graphically.


Introduction
The study of fluid flow between two parallel disks is of practical importance in many fields such as machine storage devices, computer devices, crystal growth processes, turbine engines, radial diffusers, lubrication, viscometry etc.The rotating disc problem was first formulated by von Kármán [1].He considered the flow of a viscous incompressible fluid under the influence of a rotating disk.Later Cochran [2] obtained asymptotic solutions to the steady hydro-dynamic problem formulated by von Karman.It is found that the disc acts are like a centrifugal fan, the fluid near the disc being thrown radially outwards.Hossain and Rahman [3] studied the steady flow between two porous rotating discs in the presence of transverse magnetic field.Hossain and Wilson [4] investigated unsteady flow of viscous incompressible fluid with temperature-dependent viscosity due to a rotating disc in the presence of transverse magnetic field and heat transfer.Wenchang, Wenxiao and Mingyu [5] have studied unsteady flows of a visco-elastic fluid with the fractional Maxwell model between two parallel plates.Maji, Ghara, Jana and Das [6] have considered unsteady MHD flow between two eccentric rotating disks.Liu, Zheng, Zhang and Zong [7] discussed the oscillating flows and heat transfer of a Generalised Oldroyed-B fluid in the presence of magnetic field.Kempegowda and Balagon-dar [8] have worked out exact solutions of non-Newtonian fluid flow between two moving parallel disks with stability analysis.
In the present paper we have considered unsteady incompressible visco-elastic flow of a generalized Maxwell fluid between two rotating infinite coaxial circular disks.In the aforesaid problems, time derivative of integer order has been considered in the Navier-Stokes equation but in the present problem we have considered the constitutive equation for Maxwell fluid with fractional order time derivative instead of integer order time derivative.In the constitutive equation the time derivative of integer order is replaced by the Caputo fractional calculus operator.We have obtained the analytical solution to the velocity field in series involving Mittag-Leffler function and illustrated graphically the dependence of the velocity field on the fractional and material parameters.

Generalised Maxwell Model and Basic Equation
The constitutive equation of a Generalised incompressible Maxwell fluid can be written as, where, τ is the shear stress, The equation of motion in the absence of the body force can be written as where ρ is the density of the fluid, V is the fluid velocity, D Dt is the material derivative, τ is the stress tensor.
The equation of continuity is given by 0

Mathematical Formulation
Let an incompressible viscous Generalised Maxwell fluid be bounded by two coaxial infinite parallel circular disks at a distance "d" apart and the fluid as well as the disks are initially at rest as shown in Figure 1.Suddenly the lower and the upper disks begin to rotate with velocities Ω and sΩ respectively and as a consequence the fluid is set in motion.Here we take the cylindrical coordinate system ( ) , where r, θ and z-coordinates are taken in the radial, cross radial and the direction joining the centers of the circular disks.We take the velocity profile of the form ( ) where e θ is the unit vector in the θ-direction.For such flows the constraint of incompressibility is automatically satisfied.For this problem the constitutive relationship becomes The momentum equation is Eliminating zθ τ between the Equations ( 6) and (7) we get the basic equation as The Equation ( 8) is the governing equation of the flow of a Generalised Maxwell fluid between two rotating infinite parallel circular disks considered in the present problem.
The boundary conditions are given by 0, , 0 The initial condition is given by ( ) Then the governing Equation ( 8) in non-dimensional variables is given by (for simplicity the dimensionless mark "'" will be neglected hereinafter). where, Let us consider the transformation given by Then in terms of new variable the governing equation becomes Subject to the boundary conditions ( ) ( ) Taking Laplace transformation and using initial condition we get from Equation (11) where, ( ) where "p" is Laplace transform parameter.
Taking finite Fourier sine transformation we get from the Equation (12) V n p * is the finite Fourier sine transformation of ( ) where 1, 2, 3, n =  Taking Laplace transformation of the boundary conditions we get, ( ) ( ) Using the above conditions we get from Equation (13) The Equation ( 14) can be written as In order to avoid the lengthy procedure of residues and contour integrals, we rewrite the Equation (15) into series form given by ( ) ) Taking inverse Laplace transformation we get from Equation ( 16) Taking inverse finite Fourier sine transformation we get from Equation (

Limiting Cases
Subject to the boundary conditions ( ) and initial condition ( ) The Equation ( 21) is the damped wave equation and it represents the governing equation of an Ordinary Maxwell fluid.
Then we get the velocity profile from the Equation ( 20) as subject to the boundary condition ( ) and initial condition ( ) The Equation ( 23) is diffusion equation and it represents the governing equation of a Classical Newtonian Fluid.
Then we get the velocity profile from the Equation (20) as ( )

Conclusions and Numerical Results
In the present paper we have found out the analytical solution to the velocity field by integral transform in series form in terms of Mittage-Leffler function for the unsteady incompressible flow of a Generalised Maxwell fluid between two rotating infinite parallel coaxial circular disks.We have got the solutions to the velocity fields for ordinary Maxwell fluid and Classical Newtonian fluid as the limiting cases of the solution of Generalised Maxwell fluid.In the constitutive equation for the Maxwell fluid the time derivative of integer order is replaced by Riemann-Liouville operator.The dependence of the velocity field on the fractional as well as material parameters has been illustrated graphically.
In Figure 2 the velocity is depicted against the distance from the lower disk along the direction of the common axis of rotation for different values of the fractional calculus parameter α.From the Figure it is evident that as the value of the parameter α increases the fluid velocity in the θ-direction increases.It can be noticed that the point of maximum velocity of the velocity curve gradually shifts towards the lower disk as α increases.In Figure 3 the velocity is plotted against the distance from the lower disk for different values of fractional calculus parameter β.As the value of β increases, the fluid velocity increases near the lower disk whereas the velocity decreases near the upper disk.The velocity is depicted against the distance from the lower disk for different values of parameters "α" and "s" in Figure 4.As "α" and "s" increase simultaneously the fluid velocity increases and is maximum near the midpoint region between the two parallel disks.It is evident from the figure that the points of maximum velocity of the velocity curve gradually shift towards the lower disk as in the case in Figure 1.
In Figure 5 the velocity is plotted against the distance from the lower disk for different values of the material parameter "ζ".It can be observed that as "ζ" increases the fluid velocity decreases near the lower disk whereas the velocity increases near the upper disk.The velocity is depicted against the distance from the lower disk for dif-   the lower disk) at which the velocity gradient is zero.The velocity curves for Classical Newtonian and Generalised Maxwell Fluids are almost parallel to the horizontal axis compared to the velocity curve for the Ordinary Maxwell Fluid.It can be seen that the velocity curve for the Generalised Maxwell Fluid has a point of local maximum.

ε
parameter, G is the shear modulus, α and β are fractional parameters such that 0 1 α β ≤ ≤ ≤ and ε is the shear strain.t D α and t D β are Caputo operators given by is the shear rate.

Figure 1 .
Figure 1.Geometry of the problem.
Now we have an important Laplace transformation of the nth order derivative of Mittag-Leffler function

Figure 2 .Figure 3 .
Figure 2. The velocity profile is depicted against the distance from the lower disk along the direction of the common axis of rotation for different values of fractional parameter α. β = 1.0, ζ = 0.03, η = 0.02, s = 0.2, t = 0.2.

Figure 4 .Figure 5 . 8 .
Figure 4.The velocity profile is depicted against the distance from the lower disk along the direction of the common axis rotation for different values of parameters α and s. β = 1.0, ζ = 0.03, η = 0.02, t = 0.2.

Figure 6 .Figure 7 .Figure 8 .
Figure 6.The velocity profile is depicted against the distance from the lower disk along the direction of the common axis of rotation for different values of η. α = 0.2, β = 0.8, ζ = 0.03, s = 0.2, t = 0.2.