OJFDOpen Journal of Fluid Dynamics2165-3852Scientific Research Publishing10.4236/ojfd.2013.32008OJFD-32661ArticlesPhysics&Mathematics Unsteady Incompressible Viscoelastic Flow of a Generalised Maxwell Fluid between Two Rotating Infinite Parallel Coaxial Circular Disks himanBose1*UmaBasu1*Department Applied Mathematics, University of Calcutta, Kolkata, India* E-mail:dhimanbose09@gmail.com(HB);ubappmath@caluniv.ac.in(UB);0506201303025763March 23, 2013April 24, 2013 May 1, 2013© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

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.

Generalised Maxwell Fluid; Laplace Transform; Finite Fourier Sine Transform; Mittag-Leffler Function; Fractional Derivative
1. 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 ﬁrst formulated by von Kármán . He considered the flow of a viscous incompressible fluid under the influence of a rotating disk. Later Cochran  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  studied the steady flow between two porous rotating discs in the presence of transverse magnetic field. Hossain and Wilson  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  have studied unsteady flows of a visco-elastic fluid with the fractional Maxwell model between two parallel plates. Maji, Ghara, Jana and Das  have considered unsteady MHD flow between two eccentric rotating disks. Liu, Zheng, Zhang and Zong  discussed the oscillating flows and heat transfer of a Generalised Oldroyed-B fluid in the presence of magnetic field. Kempegowda and Balagondar  have worked out exact solutions of nonNewtonian 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.

2. Generalised Maxwell Model and Basic Equation

The constitutive equation of a Generalised incompressible Maxwell fluid can be written as,

where, is the shear stress, is a relaxation parameter, G is the shear modulus, α and β are fractional parameters such that and is the shear strain. and are Caputo operators given by

For the Equation (1) gives Ordinary Maxwell fluid model and for , , a Classical Newtonian fluid model is recovered.

The Equation (1) can be rewritten as

where is the shear rate.

The equation of motion in the absence of the body force can be written as

where is the density of the fluid, is the fluid velocity, is the material derivative, is the stress tensor.

The equation of continuity is given by

3. 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 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 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 at  at “s” is some constant.

The initial condition is given by Now let us introduce the dimensionless variables Then the governing Equation (8) in non-dimensional variables is given by (for simplicity the dimensionless mark “'” will be neglected hereinafter).

where, The boundary conditions in non-dimensional variable becomes and Let us consider the transformation given by

Then in terms of new variable the governing equation becomes

Subject to the boundary conditions and and initial condition Taking Laplace transformation and using initial condition we get from Equation (11)

where, is the Laplace transformation of defined by where “p” is Laplace transform parameter.

Taking finite Fourier sine transformation we get from the Equation (12)

where, is the finite Fourier sine transformation of defined by where Taking Laplace transformation of the boundary conditions we get, and 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 (16)

Now we have an important Laplace transformation of the nth order derivative of Mittag-Leffler function given by

Taking inverse Laplace transformation we get from Equation (16)

Taking inverse finite Fourier sine transformation we get from Equation (19) Changing the variable to by the transformation we get the expression for the velocity field as follows,

4. Limiting Cases

Case-I If then the equation of motion is

Subject to the boundary conditions and 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 (22)

Case-II If , the equation of motion is given by

subject to the boundary condition and 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

5. 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-

ferent values of material parameter “η” in Figure 6. It can be noticed that the velocity decreases with the increasing values of the parameter “η” and the point of maximum velocity of the curve shifts towards the lower disks. The velocities are maximum near the lower disk for the cases. In Figure 7, as “z” increases from 0 to 0.5 the fluid velocity decreases negatively whereas the velocity increases positively as “z” increases from 0.6 to 1.0. The velocity is depicted against the distance from the lower disk for three different cases namely Case-I Ordinary Maxwell Fluid, Case-II Classical Newtonian Fluid and Case III Generalised Maxwell fluid in Figure 8. In CaseI and in Case-II and in Case-III . In the case for Ordinary Maxwell Fluid, there is a point of local maximum (near

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.

REFERENCESReferencesT. Von Kármán, “über Laminare und Turbulente Reibung,” Zeitschrift für Angewandte Mathematik und Mechanik, Vol. 1, No. 4, 1921, pp. 233-255. doi:10.1002/zamm.19210010401W. G. Cochran, “The Flow Due to a Rotating Disc,” Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 30, No. 3, 1934, pp. 365-375. doi:10.1017/S0305004100012561M. A. Hossain and A. F. M. A Rahman, “The Steady Flow between Two Porous Rotating Discs in Presence of Transverse Magnetic Field,” Indian Journal of Pure and Applied Mathmatics, Vol. 15, No. 2, 1984, pp. 187-194.M. A. Hossain, A. Hossain and M. Wilson, “Unsteady Flow of Viscous Incompressible Fluid with Temperature-Dependent Viscosity Due to a Rotating Disc in Presence of Transverse Magnetic Field and Heat Transfer,” International Journal of Thermal Sciences, Vol. 40, No. 1, 2001, pp. 11-20. doi:10.1016/S1290-0729(00)01183-2T. Wenchang, P. Wenxiao and X. Mingyu, “The Unsteady Flows of a Visco-Elastic Fluid with the Fractional Maxwell Model between Two Parallel Plates,” International Journal of Non-Linear Mechanics, Vol. 38, No. 5, 2003, pp. 645-650. doi:10.1016/S0020-7462(01)00121-4S. L.Maji, N. Ghara, R. N. Jana and S. Das, “Unsteady MHD Flow between Two Eccentric Rotating Disks,” Journal of Physical Science, Vol. 13, 2009, pp. 87-96.Y. Liu, L. Zheng, X. Zhang and F. Zong, “The Oscillating Flows and Heat Transfer of a Generalized Oldroyed-B Fluid in Magnetic Field,” International Journal of Applied Mathematics, Vol. 40, 2010.M. Kempegowda and P. M. Balagondar, “The Exact Solutions of Non-Newtonian Fluid Flow between Two Moving Parallel Disks and Stability Analysis,” Applied Mathrmatical Sciences, Vol. 6, No. 37, 2012, pp. 1827-1835.