Hindawi Publishing Corporation International Journal of Differential Equations Volume 2011, Article ID 193813,19 pages doi:10.1155/2011/193813

Research Article

Slip Effects on Fractional Viscoelastic Fluids

Muhammad Jamil1,2 and Najeeb Alam Khan3

1 Abdus Salam School of Mathematical Sciences, GC University, Lahore 54600, Pakistan

2 Department of Mathematics, NED University of Engineering and Technology, Karachi 75270, Pakistan

3 Department of Mathematics, University of Karachi, Karachi 75270, Pakistan

Correspondence should be addressed to Muhammad Jamil, jqrza26@yahoo.com Received 23 May 2011; Accepted 7 September 2011 Academic Editor: Wen Chen

Copyright © 2011 M. Jamil and N. A. Khan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Unsteady flow of an incompressible Maxwell fluid with fractional derivative induced by a sudden moved plate has been studied, where the no-slip assumption between the wall and the fluid is no longer valid. The solutions obtained for the velocity field and shear stress, written in terms of Wright generalized hypergeometric functions p'Vq, by using discrete Laplace transform of the sequential fractional derivatives, satisfy all imposed initial and boundary conditions. The no-slip contributions, that appeared in the general solutions, as expected, tend to zero when slip parameter is 0 ^ 0. Furthermore, the solutions for ordinary Maxwell and Newtonian fluids, performing the same motion, are obtained as special cases of general solutions. The solutions for fractional and ordinary Maxwell fluid for no-slip condition also obtained as limiting cases, and they are equivalent to the previously known results. Finally, the influence of the material, slip, and the fractional parameters on the fluid motion as well as a comparison among fractional Maxwell, ordinary Maxwell, and Newtonian fluids is also discussed by graphical illustrations.

1. Introduction

There are many fluids in industry and technology whose behavior cannot be explained by the classical linearly viscous Newtonian model. The departure from the Newtonian behavior manifests itself in a variety of ways: non-Newtonian viscosity (shear thinning or shear thickening), stress relaxation, nonlinear creeping, development of normal stress differences, and yield stress [1]. The Navier-Stokes equations are inadequate to predicted the behavior of such type of fluids; therefore, many constitutive relations of non-Newtonian fluids are proposed [2]. These constitutive relations give rise to the differential equations, which, in general, are more complicated and higher order than the Navier-Stokes equations. Therefore, it is difficult to obtain exact analytical solutions for non-Newtonian fluids [3]. Modeling of the equation governing the behaviors of non-Newtonian fluids in different

circumstance is important from many points of view. For examples, plastics and polymers are extensively handeled by the chemical industry, whereas biological and biomedical devices like hemodialyser make use of the rheological behavior [4]. In general, the analysis of the behavior of the fluid motion of non-Newtonian fluids tends to be much more complicated and subtle in comparison with that of the Newtonian fluids [5].

The fractional calculus, almost as old as the standard differential and integral one, is increasingly seen as an efficient tool and subtle frame work within which useful generalization is quite long and arguments almost yearly. It includes fractal media, fractional wave diffusion, fractional Hamiltonian dynamics, and biopolymer dynamics as well as many other topics in physics. Fractional calculus is useful in the field of biorheology and bioengineering, in part, because many tissue-like materials (polymers, gels, emulsions, composites, and suspensions) exhibit power-law responses to an applied stress or strain [6, 7]. An example of such power-law behavior in elastic tissue was observed recently for viscoelastic measurements of the aorta, both in vivo and in vitro [8, 9], and the analysis of these data was most conveniently performed using fractional order viscoelastic models. The starting point of the fractional derivative model of non-Newtonian model is usually a classical differential equation which is modified by replacing the time derivative of an integer order by the so-called Riemann-Liouville/Caputo fractional calculus operators. This generalization allows one to define precisely noninteger order integrals or derivatives. In general, fractional model of viscoelastic fluids is derived from well-known ordinary model by replacing the ordinary time derivatives, to fractional order time derivatives and this plays an important role to study the valuable tool of viscoelastic properties. We include here some investigation [10-18] in which the fractional calculus approach has been adopted for the flows of non-Newtonian fluids. Furthermore, the one-dimensional fractional derivative Maxwell model has been found very useful in modeling the linear viscoelastic response of some polymers in the glass transition and the glass state [19]. In other cases, it has been shown that the governing equations employing fractional derivatives are also linked to molecular theories [20]. The use of fractional derivatives within the context of viscoelasticity was firstly proposed by Germant [21]. Later, Bagley and Torvik [22] demonstrated that the theory of viscoelasticity of coiling polymers predicts constitutive relations with fractional derivatives, and Makris et al. [23] achieved a very good fit of the experimental data when the fractional derivative Maxwell model has been used instead of the Maxwell model for the silicon gel fluid. Furthermore, it is worth pointing out that Palade et al. [24] developed a fully objective constitutive equation for an incompressible fluidreducible to the linear fractional derivative Maxwell model under small deformations hypothesis.

A general view of the literature shows that the slip effects on the flows of non-Newtonian fluids has been given not much attention. Especially, polymer melts exhibit a macroscopic wall slip. The fluids exhibiting a boundary slip are important in technological applications, for example, the polishing of artificial heart valves, rarefied fluid problems, and flow on multiple interfaces. In the study of fluid-solid surface interactions, the concept of slip of a fluid at a solid wall serves to describe macroscopic effects of certain molecular phenomena. When the molecular mean free path length of the fluid is comparable to the distance between the plates as in nanochannels or microchannels, the fluid exhibits non-continuum effects such as slipflow, as demonstrated experimentally by Derek et al. [25]. Experimental observations show that [26-28] non-Newtonian fluids, such as polymer melts, often exhibit macroscopic wall slip, which, in general, is described by a nonlinear and nonmonotone relation between the wall slip velocity and the traction. A more realistic class of slip flows are those in which the magnitude of the shear stress reaches some critical value,

here called the slip yield stress, before slip occurs. In fact, some experiments show that the onset slip and slip velocity may also depend on the normal stress at the boundary [26, 29]. Much of the research involving slip presumes that the slip velocity depends on the shear stress. The slip condition is an important factor in sharskin, spurt, and hysteresis effects, but the existing theory for non-Newtonian fluids with wall slippage is scant. We mention here some recent attempt regarding exact analytical solutions of non-Newtonian fluids with slip effects [30-36].

The objective of this paper is twofold. Firstly, is to give few more exact analytical solutions for viscoelastic fluids with fractional derivative approach, which is more natural and appropriate tool to describe the complex behavior of such fluids. Secondly, is to study the slip effects on viscoelastic fluid flows, which is important due to their practical applications. More precisely, our aim is to find the velocity field and the shear stress corresponding to the motion of a Maxwell fluid due to a sudden moved plate, where no-slip assumption is no longer valid. However, for completeness, we will determine exact solutions for a larger class of such fluids. Consequently, motivated by the above remarks, we solve our problem for Maxwell fluids with fractional derivatives. The general solutions are obtained using the discrete Laplace transforms. They are presented in series form in terms of the Wright generalized hypergeometric functions p¥q and presented as sum of the slip contribution and the corresponding no-slip contributions. The similar solutions for ordinary Maxwell fluids can easily be obtained as limiting cases of general solutions. The Newtonian solutions are also obtained as special cases of fractional and ordinary Maxwell fluids. Furthermore, the solutions for fractional and ordinary Maxwell fluid for no-slip condition also obtained as a special cases, and they are similar with previously known results in the literature. Finally, the influence of the material, slip and fractional parameters on the motion of fractional and ordinary Maxwell fluids is underlined by graphical illustrations. The difference among fractional Maxwell, ordinary Maxwell, and Newtonian fluid models is also highlighted.

2. The Differential Equations Governing the Flow

The equations governing the flow of an incompressible fluid include the continuity equation and the momentum equation. In the absence of body forces, they are

V-V = 0, V-T = p^ + p(V -V)V, (2.1)

where p is the fluid density, V is the velocity field, t is the time, and V represents the gradient operator. The Cauchy stress T in an incompressible Maxwell fluid is given by [10,11,14-17]

T = -pI + S, S + X(js - LS - SLT) = [iA, (2.2)

where -pI denotes the indeterminate spherical stress due to the constraint of incompress-ibility, S is the extrastress tensor, L is the velocity gradient, A = L + LT is the first Rivlin Ericsen tensor, i is the dynamic viscosity of the fluid, 1 is relaxation time, the superscript T indicates the transpose operation, and the superposed dot indicates the material time derivative. The model characterized by the constitutive equations (2.2) contains as special

case the Newtonian fluid model for 1 ^ 0. For the problem under consideration, we assume a velocity field V and an extrastress tensor S of the form

V = V(y,t) = u(y,t)i, S = S(y,t), (2.3)

where i is the unit vector along the x-coordinate direction. For these flows, the constraint of incompressibility is automatically satisfied. If the fluid is at rest up to the moment t = 0, then

V(y, 0) = 0, S(y, 0) = 0, (2.4)

and (2.1)-(2.3) yield the meaningful equation

( -4)^=-p(1 -4)dx-(1 -4>(y-t)=

where t(y, t) = Sxy(y, t) is the nonzero shear stress and v = »/p is the kinematic viscosity of the fluid.

The governing equations corresponding to an incompressible Maxwell fluid with fractional derivatives, performing the same motion in the absence of a pressure gradient in the flow direction, are (cf. [4,15,17])

(1 - x-d;) = v^, (1 - x-d;)t (y,0 = , (2.6)

where a is the fractional parameter and the fractional differential operator so-called Caputo fractional operator Da is defined by [37,38]

D?f (t) =

f '(t ) r(1 - a)} 0 (t - t ) df (t) dt '

■dT, 0 < a < I,

a = 1,

and r(^) is the Gamma function. In the following, the system of fractional partial differential equations (2.6), with appropriate initial and boundary conditions, will be solved by means of Fourier sine and Laplace transforms. In order to avoid lengthy calculations of residues and contour integrals, the discrete inverse Laplace transform method will be used [10-18].

3. Statement of the Problem

Consider an incompressible Maxwell fluid with fractional derivatives occupying the space lying over an infinitely extended plate which is situated in the (x, z) plane and perpendicular to the y-axis. Initially, the fluid is at rest, and at the moment t = 0-, the plate is impulsively brought to the constant velocity U in its plane. Here, we assume the existence of slip boundary between the velocity of the fluid at the wall u(0,t) and the speed of the wall, and

the relative velocity between u(0,t) and the wall is assumed to be proportional to the shear rate at the wall. Due to the shear, the fluid above the plate is gradually moved. Its velocity is of the form (2.3) i while the governing equations are given by (2.6). The appropriate initial and boundary conditions are [39]

u{y,0) = dud0 = 0; t{y,0) = 0, y> 0, (3.1)

u(0,t) = UH(t) + 9H(t)

öu(y, t)

where H(t) is the Heaviside function and 0 is the slip strength or slip coefficient. If 0 = 0, then the general assumed no-slip boundary condition is obtained. If 0 is finite, fluid slip occurs at the wall, but its effect depends upon the length scale of the flow. Furthermore, the natural conditions

öu(y, t) dy

■ 0 as y —> œ, t > 0

have to be also satisfied. They are consequences of the fact that the fluid is at rest at infinity, and there is no shear in the free stream.

4. Solution of the Problem

4.1. Calculation of the Velocity Field

Applying the Laplace transform to (2.6)i, using the Laplace transform formula for sequential fractional derivatives [37, 38], and taking into account the initial conditions (3.1) 1/2, we find that

öü(y,q) q( 1 + Xaqa) dy2 v

u{yq) = a

subject to the boundary conditions ïï((M) = U +

, u(y,q),

ö^y,q) dy

0 as y —> œ,

where u(y,q) is the image function of u(y,t) and q is a transform parameter. Solving (4.1) and (4.2), we get

q{ 1 + d[q( 1 + Xaqa)/v\1/2}

q(1 +Xaqa)

In order to obtain u(y,t) = L-1{u(y,q)} and to avoid the lengthy and burdensome calculations of residues and contours integrals, we apply the discrete inverse Laplace transform method [10-18]. However, for a suitable presentation of the velocity field, we firstly rewrite (4.3) in series form

U,, TVf aJ!a \ V r(m - (k/2))(-1-a)m 1

l(y,q) u + d]j v J E; mir(-fc/2) q-(a+1)(k/2)+am+1

# ( v k+m

oo oo m

^ia -in ir(n - ((k+m)/2))(-ra)n 1

k=0 m=i m! \ \v J n=0 n!r(-(k - m)/2) q-(a-i)(k-m/2)-an-i' k (4.4)

where we use the fact that

k r(a - 1) = r(k - a) (4 5) ( 1) r(a - k -1) r(-a) ' (4.5)

Inverting (4.4) by means of discrete inverse Laplace transform, we find that u(y,i) = UHO-UH(f)i(-ejv)V<a-1»k/2>gm^^iZn - 1)

/ I-\ k-m

v v ...m f ta \

W ^ \ t-(a-1)((k-m)/2)

m! \ V v

UH(t)£ekZm (-VV ) t

k=0 m=1

1_r(n - ((k + m)/2))(-ta/Xa)n_

x n=0n!r(-(k + m)/2)r(-((k + m)/2)(a + 1) + an + 1)'

In term of Wright generalized hypergeometric function [40], we rewrite the above equation as a simpler form

u(y,t) = UH(t) + UH(t)^[ -0V— ) t-(a+1)k/11^2

k=\ V V /

(-k/2,1)

(-k/2,0)'(-(k/2)(a+1)+1,a)

/ i-\ k+m

o o -..mi \a \

-UH(t)Zo^HV J t-(a+1)(k+m/2) lT2

k=0 m=1 m' \ V v J

(-(k+m)/2,1)

(-(k+m)/2'0),(-((k+m)/2)(a+1) + 1 ,a)

where the Wright generalized hypergeometric function p¥q is defined as

(ai,Ai),...(ap,Ap ) (bi,Bi ),...(b„B,)

(z)"nP=ir(a; + Ajn)

n!H J=i^ bj + Bjn)

In order to justify the initial conditions (3.1)1/2, we use the initial value theorem of Laplace transform [41]

u(y, o) = Jim u(:y,t) = lhmJ qu(y>q)\ = 0,

dtu(y, 0) = limdtu(y, t) = lim \q2u(y,q) - qu(y, 0) = 0.

t^0 qL J

Furthermore, to justify the boundary condition (3.2), we have

du{y,t)

= UH(t)fi ( -0

-(a+i)((k+i)/2)

= umz -M - t-

-(a+i)(k/2)

(—(k+i)/2,i)

(-(k+i)/2,0),(-((k+i)/2)(a+i) + i,a)

(—(k+i)/2,i)

(-(k+i) /2,0),(-((k+i) /2) (a+i)+i,a)

i(0,t) = UH(t) + UH(t)% ( -0

-(a+i)(k/2)

(-(k/2),1)

(-k/2,0),(-(k/2)(a+i) + i,a)

(4.10)

It is easy to see that the exact solution (4.7) satisfies the boundary condition (3.2).

4.2. Calculation of the Shear Stress

Applying the Laplace transform to (2.6)2 and using the initial condition (3.i)3, we find that

%y'q) =

P du(y,q) i + Xaqa dy '

(4.11)

where t(y,q) is the Laplace transform of t(y,t). Using (4.3) in (4.11), we find that

T{y,q) = -

pU[q(i + Xaqa)]

v{ i + 0[q( i + Xaqa)/v\1/2 }

q(i+Xaqa)

(4.i2)

in order to obtain t(y, t) under the suitable form, we write (4.12) in series form

^ = - »U ^ g -V ^

v k=0 m=0m!

(4.13)

r(n - (k - m - 1)/2)(-1-a)n

n!r(-(k - m - 1)/2) q-(a-1)((k-m-1)/2)+

Inverting (4.13) by means of the discrete inverse Laplace transform, we find the shear stress t(y, t) under simple form

v V y.m

-(y,t) = -pUH(t)£ ek^y-

k=0 m=0 '"■■

-(a-1)((m-k-1)/2)-1

r(n - (k - m - 1)/2)(-ta/1a)n

n!r(-(k - m - 1)/2)r(-(a - 1)(k - m - 1)/2 - an)'

(4.14)

or equivalently

v v ym -(yj) = -puh (t)£ ekYjm

k=0 m=0

-(a-1)((m-k-1)/2)-1

tL ~1a

(-(k-m-1)/2,1)

(-(k-m-1)/2,0),(-((k-m-1)/2)(a-1),a)

(4.15)

5. The Special Cases

5.1. Ordinary Maxwell Fluid with Slip Effects

Making a ^ 1 into (4.7) and (4.15), we obtain the velocity field

i(y,t) = UH(t) - UH(t)g y-6)!1 y t-k 1^2

t (-k/2,1)

-1 (-k/2,0),(-k-1,1)

UH(t)£ ekZ

k=0 m=1

-(k-m)

(-(k-m)/2,1)

(-(k-m)/2,0),(-(k-m)-1,1)

and the associated shear stress

œ œ ^ m

ïkSrV

-(y,t) = -pUH(t)£e^m( -yx ) r(k+m) ^

k=0 m=0\ V v

(-(k+m-1)/2,1) (-(k+m-1)/2,0),(-(k+m) + 1,1)

corresponding to an ordinary Maxwell fluid performing the same motion.

5.2. Fractional Maxwell Fluid without Slip Effects

Making 0 ^ 0 into (4.7) and (4.15), we obtain the solutions for velocity field

i{yt = uH(t)+uH(t)£y— ( -yv ) t

-(a+1)(m/2)

(-m/2,1)

(-m/2,0),(-(m/2)(a+1) + 1,a)

and the associated shear stress

t(y,t)=-puH(t)zm) t-(«+1)((m-1)/2)-11^2

m=0 \ " /

(-(m-1)/2,1)

(-(m-1) /2,0),(-((m-1) /2) (a+1),a)

they are equivalent to the known solutions obtained in [42, 43] for Sokes' first problem of fractional Maxwell fluid.

5.3. Ordinary Maxwell Fluid without Slip Effects

Making a ^ i into (5.3) and (5.4), we recover the solutions for velocity field shear stress for Stokes' first problem of ordinary Maxwell fluid.

5.4. Newtonian Fluid with Slip Effects

Finally, making X ^ 0 into (4.3) and (4.i2), the solutions for a Newtonian fluid with slip effects are obtained

„(y,0 = UW./2,.( --V) + Ug( — ) W./2,-(k/2,+,v ^

in which

Wa,b (z) = X

=0n!r(an + b) '

0 = 0.5 eœ t = 2s t=3s

t = 4s

0 = 0.5 eee t = 2s t = 3s

t = 4s

Figure 1: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, v = 0.379, » = 33,1 = 1.5, a = 0.2, 6 = 0.5, and different values of t.

is the Wright function [40]. Using the definition of Wright function and the series expression of error function, we can easily prove that

W-1/2,1 (-z) = erfc(2)' W-1/2,1/2(-z) = exp^-•

Substituting (5.7) into (5.5), we can reduce to

(y,t) = UN (y,t) + UjT (-"L) W-1/2,-(k/2)+^-"=)'

-JL- \

uU œ / 0 \K

■ty,t) = Ny'O - E(""^j W-1/2,-(k-1)/2

un (y, t) = U erfc tn (y,t) = -

uu „

"nvt6*Pl 4vt

are classical solutions for Stokes' first problem of Newtonian fluid [42,44].

6. Numerical Results and Conclusions

In this paper, the unsteady flow of fractional Maxwell fluid over an infinite plate, where the no-slip assumption between the wall and the fluid is no longer valid, is studied by means

eee t = 2s t=3s

t = 4s

eee t = 2s

t = 4s

Figure 2: Profiles of the velocity field u(y, t) and the shear stress t (y, t) for fractional Maxwell fluid given by (5.3) and (5.4), for U = 1, y = 1, v = 0.379, p = 33,1 = 1.5, a = 0.2, 0 = 0.0, and different values of t.

e=o e = o.5

0.25 3.19 y=1

e=o e = o.5

6.13 t

9.o6 12

Figure 3: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, v = 0.379, u = 33,1 = 1.5, a = 0.2, and different values of 0.

of the discrete Laplace transforms. The motion of the fluid is due to the plate that at time t = 0+ is suddenly moved with a constant velocity U in its plane. Closed-form solutions are obtained for the velocity u(y,t) and the shear stress t (y,t) in series form in terms of the Wright generalized hypergeometric functions. These solutions, presented as a sum of the slip contribution and the corresponding no-slip contributions, satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Maxwell fluids are also

eee 0 = 0 *** 0 = 0.5

-1.9 -

3.19 6.13 t

eee 0=0 *** 0 = 0.5

9.06 12

Figure 4: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 2, v = 0.379, » = 33,1 = 1.5, a = 0.2, and different values of 6.

0 = 0.5 ooo A = 1.5 *** A=4

e-o-9- A = 12

0 = 0.5 000 A = 1.5 A=4

A = 12

Figure 5: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, v = 0.379, » = 33, a = 0.2, 6 = 0.5, t = 5 s, and different values of 1.

obtained from general solutions for a ^ 1. In the special case when 6 ^ 0, the general solution reduces to previously known results for Stokes' first problem.

In order to reveal some relevant physical aspects of the obtained results, the diagrams of the velocity field u(y,t) and the shear stress t(y,t) have been drawn against y and t for different values of t, material constants v, 1, slip parameter 6, and fractional parameter a. From all figures, it is clear that increasing the slip parameter at the wall the velocity decreases at the wall. Figures 1 and 2 are prepared to show the effect of time on velocity and shear

OOO A = 1.5 ooo A = 12

Figure 6: Profiles of the velocity field u(y, t) and the by (5.3) and (5.4), for U = 1, y = 1, v = 0.379, p = 33,

ooo A = 1.5 ooo A = 12

*** A=4

shear stress t (y, t) for fractional Maxwell fluid given a = 0.2, 0 = 0.0, t = 5 s, and different values of 1.

0 = 0.5 eee a = 0.2 *** a = 0.4

«-a-* a = 0.6

-8.7 -

0 = 0.5 eee a = 0.2 *** a = 0.4

«-e-* a = 0.6

Figure 7: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, v = 0.379, p = 33,1 = 1.5, 0 = 0.5, t = 5 s, and different values of a.

stress profiles with and without slip effects. It is clear that velocity and shear stress (in absolute value) are smaller when slip parameter is nonzero. It is also noted that velocity on the whole flow domain while the shear stress (in absolute value) on large part of flow domain are increasing functions of time t. Figures 3 and 4 are sketched to see the influence of slip effects on fluid motion for two different values of y. It is noted that velocity and shear stress (in absolute value), as expected are decrease when slip parameter 0 and y increase.

eee a = 0.2 a = 0.4

a = 0.6

a = 0.2 a = 0.4

a = 0.6

Figure 8: Profiles of the velocity field u(y,t) and the shear stress t(y,t) for fractional Maxwell fluid given by (5.3) and (5.4), for U = 1, y = 1, v = 0.379,» = 33,1 = 1.5, 6 = 0.0, t = 5 s, and different values of a.

0 = 0.5 eee v = 0.2 v = 0.3

v = 0.4

0 = 0.5

v = 0.2 v = 0.3

v = 0.4

Figure 9: Profiles of the velocity field u(y,t) and the shear stress t(y,t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, p = 87.071,1 = 1.5, a = 0.2, 6 = 0.5, t = 5 s, and different values of v.

The influence of relaxation time and kinematic viscosity v on fluid motion are presented in Figures 5-8. As expected, the two material parameter have opposite effects on fluid motion. For instance, the velocity and shear stress decreases with respect to 1. More important for us is to see the effects of fractional and slip parameter on fluid motion. It is observed that velocity and shear stress either slip effects present or not are increasing functions of fractional parameter, as shown in Figures 9 and 10. The effect of slip parameter is clear from Figure 11.

eee v = 0.2

v = 0.3

v = 0.4

v = 0.2

v = 0.3

v = 0.4

Figure 10: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (5.3) and (5.4), for U = 1, y = 1, p = 87.071,1 = 1.5, a = 0.2, 0 = 0.0, t = 5 s, and different values of v.

eee 0 = 0 0 = 0.5

eee 0=0 0 = 0.5

o-o-e- 0 = 1

Figure 11: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell fluid given by (4.7) and (4.15), for U = 1, y = 1, v = 0.379, p = 33,1 = 1.5, a = 0.2, t = 5 s, and different values of 0.

Finally, for comparison, the velocity field and the shear stress corresponding to the three models (fractional Maxwell, ordinary Maxwell, and Newtonian) are together depicted in Figures 12-14 for three different values of slip parameter and the same values of t and of the material constants. It is clearly seen from Figures 12 and 13 that the ordinary Maxwell fluid swiftest and the fractional Maxwell fluid is the slowest near the moving plate for slip parameters 0 = 0.0 and 0 = 0.2. However, the monotonicity is change on large part of the flow domain. The shear stress corresponding to ordinary Maxwell fluid is highest near the

Fractional Maxwell eee Ordinary Maxwell Newtonian

*** Fractional Maxwell Ordinary Maxwell Newtonian

Figure 12: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell, ordinary Maxwell and Newtonian fluids, for U = 1, y = 1, v = 0.379,» = 33,1 = 1.5, a = 0.2, t = 2 s, and 6 = 0.0.

*** Fractional Maxwell Ordinary Maxwell Newtonian

0 = 0.2

*** Fractional Maxwell Ordinary Maxwell Newtonian

Figure 13: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell, ordinary Maxwell and Newtonian fluids, for U = 1, y = 1, v = 0.379,» = 33,1 = 1.5, a = 0.2, t = 2 s, and 6 = 0.2.

moving plate. For higher values of slip parameter the fractional Maxwell fluid is swiftest and the ordinary Maxwell is slowest and shear stress corresponding to fractional Maxwell fluid is largest on the whole flow domain as it is clear from Figure 14. It is important to note the difference between fractional and ordinary Maxwell fluid that, when slip effect is not present, the ordinary Maxwell fluid have oscillating behavior near the moving plate as

0.125 -1-1-1--0

0.042 -

±±± Fractional Maxwell eee Ordinary Maxwell ^^ Newtonian

Fractional Maxwell Ordinary Maxwell Newtonian

Figure 14: Profiles of the velocity field u(y, t) and the shear stress t(y, t) for fractional Maxwell, ordinary Maxwell and Newtonian fluids, for U = 1, y = 1, v = 0.379, f = 33, X = 1.5, a = 0.2, t = 2 s, and 0 = 5.0.

shown in Figure 12, which is the natural one, ordinary Maxwell fluid being the viscoelastic fluid. However the fractional Maxwell fluid have no oscillation. The units of the material constants in all figures are SI units.

Acknowledgments

The authors would like to express their sincere gratitude to the editor and referees for their careful assessment and fruitful remarks and suggestions regarding the initial version of the paper. M. Jamil is highly thankful and grateful to the Abdus Salam School of Mathematical Sciences, GC University, Lahore, Pakistan, Department of Mathematics, NED University of Engineering & Technology, Karachi 75270, Pakistan, and also Higher Education Commission of Pakistan for generous support and facilitating this research work. The author N. A. Khan is highly thankful and grateful to the Dean of Faculty of Sciences, University of Karachi, Karachi 75270, Pakistan for generous support and facilitating this research work.

References

[1] C. Fetecau, M. Khan, C. Fetecau, and H. T. Qi, "Exact solutions for the flow of a generalized Oldroyd-B fluid induced by a suddenly moved plate between two side walls perpendicular to the plate," Proceeding of the Romanian academy, Series A, vol. 11, pp. 3-10, 2010.

[2] J. Zierep and C. Fetecau, "Energetic balance for the Rayleigh-Stokes problem of a Maxwell fluid," International Journal of Engineering Science, vol. 45, no. 2-8, pp. 617-627, 2007.

[3] T. Hayat, A. H. Kara, and E. Momoniat, "Travelling wave solutions to Stokes' problem for a fourth grade fluid," Applied Mathematical Modelling, vol. 33, no. 3, pp. 1613-1619, 2009.

[4] L. Zheng, F. Zhao, and X. Zhang, "Exact solutions for generalized Maxwell fluid flow due to oscillatory and constantly accelerating plate," Nonlinear Analysis: Real World Applications, vol. 11, no. 5, pp. 3744-3751, 2010.

[5] T. Hayat, S. Nadeem, and S. Asghar, "Periodic unidirectional flows of a viscoelastic fluid with the fractional Maxwell model," Applied Mathematics and Computation, vol. 151, no. 1, pp. 153-161, 2004.

R. S. Lakes, Viscoelastic Solids, CRC Press, Boca Raton, Fla, USA, 1999. R. G. Larson, The Structure and Rheology of Complex Fluids, Oxford University Press, NewYork, NY, USA, 1999.

D. Craiem, F. J. Rojo, J. M. Atienza, R. L. Armentano, and G. V. Guinea, "Fractional-order viscoelasticity applied to describe uniaxial stress relaxation of human arteries," Physics in Medicine and Biology, vol. 53, no. 17, pp. 4543-4554, 2008.

D. Craiem and R. L. Armentano, "A fractional derivative model to describe arterial viscoelasticity," Biorheology, vol. 44, no. 4, pp. 251-263, 2007.

M. Khan, S. Hyder Ali, C. Fetecau, and H. Qi, "Decay of potential vortex for a viscoelastic fluid with fractional Maxwell model," Applied Mathematical Modelling, vol. 33, no. 5, pp. 2526-2533, 2009. A. Mahmood, S. Parveen, A. Ara, and N. A. Khan, "Exact analytic solutions for the unsteady flow of a non-Newtonian fluid between two cylinders with fractional derivative model," Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 8, pp. 3309-3319, 2009.

S. Wang and M. Xu, "Axial Couette flow of two kinds of fractional viscoelastic fluids in an annulus," Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 1087-1096, 2009.

C. Fetecau, M. Jamil, C. Fetecau, and D. Vieru, "The Rayleigh-Stokes problem for an edge in a generalized Oldroyd-B fluid," Zeitschrift fur Angewandte Mathematik und Physik, vol. 60, no. 5, pp. 921-933, 2009.

C. Fetecau, A. Mahmood, and M. Jamil, "Exact solutions for the flow of a viscoelastic fluid induced by a circular cylinder subject to a time dependent shear stress," Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 12, pp. 3931-3938, 2010.

L. C. Zheng, K. N. Wang, and Y. T. Gao, "Unsteady flow and heat transfer of a generalized Maxwell fluid due to a hyperbolic sine accelerating plate," Computers and Mathematics with Applications, vol. 61, no. 8, pp. 2209-2212, 2011.

C. Fetecau, M. Jamil, and A. Mahmood, "Flow of fractional Maxwell fluid betweencoaxial cylinders," Archive of Applied Mechanics, vol. 81, pp. 1153-1163, 2011.

M. Jamil, A. Rauf, A. A. Zafar, and N. A. Khan, "New exact analytical solutions for Stokes firstproblem of Maxwell fluid with fractional derivative approach," Computers & Mathematics with Applications, vol. 62, pp. 1013-1023, 2011.

M. Jamil, N. A. Khan, and A. A. Zafar, "Translational flows of an Oldroyd-B fluid with fractional derivatives," Computers and Mathematics with Applications, vol. 62, no. 3, pp. 1540-1553, 2011. A. Heibig and L. I. Palade, "On the rest state stability of an objective fractional derivative viscoelastic fluid model," Journal of Mathematical Physics, vol. 49, no. 4, Article ID 043101, 2008. C. Friedrich, "Relaxation and retardation functions of the Maxwell model with fractional derivatives," Rheologica Acta, vol. 30, no. 2, pp. 151-158,1991.

A. Germant, "On fractional differentials," Philosophical Magazine, vol. 25, pp. 540-549,1938.

R. L. Bagley and P. J. Torvik, "A theoretical basis for the applications of fractional calculus to

viscoelasticity," Journal of Rheology, vol. 27, no. 3, pp. 201-210,1983.

N. Makris, G. F. Dargush, and M. C. Constantinou, "Dynamic analysis of generalized viscoelastic fluids," Journal of Engineering Mechanics, vol. 119, no. 8, pp. 1663-1679,1993.

L. I. Palade, P. Attane, R. R. Huilgol, and B. Mena, "Anomalous stability behavior of a properly invariant constitutive equation which generalises fractional derivative models," International Journal of Engineering Science, vol. 37, no. 3, pp. 315-329,1999.

C. Derek, D. C. Tretheway, and C. D. Meinhart, "A generating mechanism for apparent fluid slip in hydrophobic microchannels," Physics of Fluids, vol. 16, no. 5, article 1509, 7 pages, 2004. G. V. Vinogradov and L. I. Ivanova, "Wall slippage and elastic turbulence of polymers in the rubbery state," Rheologica Acta, vol. 7, no. 3, pp. 243-254,1968.

S. Luk, R. Mutharasan, and D. Apelian, "Experimental observations of wall slip: tube and packed bed flow," Industrial and Engineering Chemistry Research, vol. 26, no. 8, pp. 1609-1616,1987. K. B. Migler, H. Hervet, and L. Leger, "Slip transition of a polymer melt under shear stress," Physical Review Letters, vol. 70, no. 3, pp. 287-290,1993.

S. G. Hatzikiriakos and J. M. Dealy, "Wall slip of molten high density polyethylenes. II. Capillary rheometer studies," Journal of Rheology, vol. 36, pp. 703-741,1992.

A. R. A. Khaled and K. Vafai, "The effect of the slip condition on Stokes and Couette flows due to an oscillating wall: exact solutions," International Journal of Non-Linear Mechanics, vol. 39, no. 5, pp. 795-809, 2004.

[31] T. Hayat, M. Khan, and M. Ayub, "The effect of the slip condition on flows of an Oldroyd 6-constant fluid," Journal of Computational and Applied Mathematics, vol. 202, no. 2, pp. 402-413, 2007.

[32] M. Khan, "Partial slip effects on the oscillatory flows of a fractional Jeffrey fluid in a porous medium," Journal of Porous Media, vol. 10, no. 5, pp. 473-487, 2007.

[33] R. Ellahi, T. Hayat, and F. M. Mahomed, "Generalized couette flow of a third-grade fluid with slip: the exact solutions," Zeitschrift fur Naturforschung—Section A, vol. 65, no. 12, pp. 1071-1076, 2010.

[34] T. Hayat, S. Najam, M. Sajid, M. Ayub, and S. Mesloub, "On exact solutions for oscillatory flows in a generalized Burgers fluid with slip condition," Zeitschrift fur Naturforschung—Section A, vol. 65, no. 5, pp. 381-391, 2010.

[35] T. Hayat, S. Zaib, C. Fetecau, and C. Fetecau, "Flows in a fractional generalized Burgers' fluid," Journal of Porous Media, vol. 13, no. 8, pp. 725-739, 2010.

[36] L. Zheng, Y. Liu, and X. Zhang, "Slip effects on MHD flow of a generalized Oldroyd-B fluid withfractional derivative," Nonlinear Analysis: Real World Applications, vol. 13, no. 2, pp. 513-523,2012.

[37] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.

[38] F. Mainardi, Frcational Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathemtical Models, Imperial College Press, London, UK, 2010.

[39] I. C. Christov, "Stokes' first problem for some non-Newtonian fluids: results and mistakes," Mechanics Research Communications, vol. 37, no. 8, pp. 717-723, 2010.

[40] A. M. Mathai, R. K. Saxena, and H. J. Haubold, "Solutions of certain fractional kinetic equations and a fractional diffusion equation," Journal of Mathematical Physics, vol. 51, no. 10, Article ID 103506, 2010.

[41] L. Debnath and D. Bhatta, Integral Transforms and Their Applications, Chapman & Hall/CRC, 2nd edition, 2007.

[42] W. Tan and M. Xu, "Plane surface suddenly set in motion in a viscoelastic fluid with fractional Maxwell model," Acta Mechanica Sinica, vol. 18, no. 4, pp. 342-349, 2002.

[43] H. Qi and M. Xu, "Stokes' first problem for a viscoelastic fluid with the generalized Oldroyd-B model," Acta Mechanica Sinica/Lixue Xuebao, vol. 23, no. 5, pp. 463-469, 2007.

[44] W. C. Tan and M. Y. Xu, "The impulsive motion of flat plate in a generalized second grade fluid," Mechanics Research Communications, vol. 29, no. 1, pp. 3-9, 2002.

Copyright of International Journal of Differential Equations is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.