Applied Mathematics
Vol.06 No.01(2015), Article ID:53065,9 pages
10.4236/am.2015.61011

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

Dhiman Bose, Uma Basu

Department of Applied Mathematics, University of Calcutta, Calcutta, India

Email: dhimanbose09@gmail.com

Received 27 November 2014; revised 18 December 2014; accepted 30 December 2014

ABSTRACT

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.

Keywords:

Oldroyd-B Fluid, Exact Solution, Mittage-Lefller Function, Fractional Derivative, Transverse Magnetic Field

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

2. 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 the plates start to oscillate in its plane with the velocity and where 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, and the boundary conditions are given by,.

We take the velocity and stress of the form

(1)

where is the velocity component in the x-direction.

The constitutive relationship for the fluid associated with the present problem is given by,

(2)

In the relation (2), and are Caputo operators defined by

(3)

Figure 1. Geometry of the problem.

According to our problem,

We consider a generalized Oldroyd-B fluid between two infinite parallel disks in presence of an imposed magnetic field that acts in the direction of the positive y-axis. Then in the presence of the body force, the momentum equation is given by

(4)

where “” is constant and “” is the density of the fluid.

Eliminating between the Equations (2) and (4) we have the governing equation

(5)

where is the kinematic viscosity and

Introducing the non-dimensional quantities,

we get the governing equation in non-dimensional quantities as

(6)

(Omitting the dimensionless mark “*”)

(7)

and the boundary conditions, for

(8)

Taking finite Fourier sine transformation we get from Equation (6)

(9)

where is the finite Fourier sine transformation of.

Using the boundary conditions (8) the Equation (9) can be rewritten as

(10)

Taking Laplace transformation and using we get from the above equation

(11)

Now in order to avoid the lengthy procedure of residues and contour integrals, we rewrite the Equation (11) into series form as

(12)

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

(13)

where

(14)

Taking Laplace Inverse transformation we get from the Equation (12)

(15)

Taking inverse finite Fourier sine transformation we get the velocity profile from the Equation (18) as

(16)

3. Limiting Cases

Case-I If, then the equation of motion is given by

(17)

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

(18)

Case-II If, then the equation is given by

(19)

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

(20)

Case-III If, then the equation of motion is given by

(21)

subject to the initial and boundary conditions given by the Equations (7) and (8) respectively.

The Equation (21) is the governing equation for a generalized second grade fluid and the velocity field is given by

(22)

Case-IV If, then the equation of motion is given by

(23)

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

(27)

4. 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 2 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 minimum 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 2. The velocity is depicted against the distance from the lower plate for different values of fractional calculus parameter in Figure 3. 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 of the lower plate in Figure 5. It is evident from the figure that the fluid velocity decreases for higher values of the parameter. The initial points of the velocity curves near the lower plate are different in domain of spatial

Figure 2. The velocity profile is depicted against the distance from the lower plate for different values of the fractional calculus parameter., , , , , , , _______, ________, _______.

Figure 3. The velocity profile is depicted against the distance from the lower plate for different values of the fractional calculus parameter., , , , , , , ________, _______, _______.

Figure 4. The velocity is depicted against the distance from the lower plate for different values of parameter., , , , , , , _______, _______, ________.

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

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

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 and, , , , , , , _______, ________, _______.

variable Y for different values of frequency of oscillation of the lower plate. There is negative velocity in Figure 5 near the lower plate for values for higher frequency of oscillation of the lower plate. In Figure 6 as the frequency of the oscillation 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 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 Figure 6 for higher values of the frequency of oscillation of the upper plate. In Figure 7 the 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.. The fluid velocity increases with increasing values of and the velocity curves are oscillatory in nature. It can be noticed that for equal frequency of oscillations of the two plates, the heights of the initial and terminal points on the velocity curve in the domain of spatial variable are equal.

References

1. Bandeli, R. and Rajagopal, K.R. (1995) Start-Up Flows of Second Grade Fluids in Domains with One Finite Dimension. International Journal of Non-Linear Mechanics, 30, 817-839. http://dx.doi.org/10.1016/0020-7462(95)00035-6
2. Fetecau, C., Fetecau, C., Karman, M. and Vieru, D. (2009) 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. Journal of Non- Newtonian Fluid Mechanics, 156, 189-201. http://dx.doi.org/10.1016/j.jnnfm.2008.06.005
3. Hayat, T. and Sajid, M. (2007) Homotopy Analysis of MHD Boundary Layer Flow of an Upper-Convected Maxwell Fluid. International Journal of Engineering Science, 45, 393-401. http://dx.doi.org/10.1016/j.ijengsci.2007.04.009
4. Jamil, M. and Khan, N.A. (2011) Slip Effects on Fractional Viscoelastic Fluids. International Journal of Differential Equations, 2011, Article ID 193813. http://dx.doi.org/10.1155/2011/193813
5. Shen, F., Tan, W., Zhao, Y. and Masuoka, T. (2006) The Rayleigh-Stokes Problem for a Heated Generalised Second Grade Fluid with Fractional Derivative Model. Nonlinear Analysis: Real World Application, 7, 1072-1080. http://dx.doi.org/10.1016/j.nonrwa.2005.09.007
6. Vieru, D., Fetecau, C. and Fetecau, C. (2008) Flow of a Generalized Oldroyd-B Fluid Due to a Constantly Accelerating Plate. Applied Mathematics and Computation, 201, 834-842. http://dx.doi.org/10.1016/j.amc.2007.12.045
7. Wenchang, T., Wenxiao, P. and Mingyu, X. (2003) A Note on Unsteady Flows of a Viscoelastic Fluid with the Fractional Maxwell Model between Two Parallel Plates. International Journal of Non-Linear Mechanics, 38, 645-650. http://dx.doi.org/10.1016/S0020-7462(01)00121-4
8. Vieru, D., Fetecau, C. and Fetecau, C. (2010) Unsteady Flow of a Generalized Oldroyd-B Fluid Due to an Infinite Plate Subject to a Time-Dependent Shear Stress. Canadian Journal of Physics, 88, 675-687. http://dx.doi.org/10.1139/P10-055