Scholarly article on topic 'Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress'

Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress Academic research paper on "Mathematics"

0
0
Share paper
Academic journal
Bound Value Probl
OECD Field of science
Keywords
{""}

Academic research paper on topic "Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress"

0 Boundary Value Problems

a SpringerOpen Journal

RESEARCH

Open Access

Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress

Qamar Sultan1, Mudassar Nazar1*, Muhammad Imran2 and Usman Ali1

Correspondence: mudassar_666@yahoo.com 'Centre for Advanced Studies in Pure and Applied Mathematics, Bahauddin Zakariya University, Multan, Pakistan Full list of author information is available at the end of the article

Abstract

This paper presents the unsteady magnetohydrodynamic (MHD) flow of a generalized Burgers' fluid between two parallel side walls perpendicular to a plate. The plate applies a shear stress induced by rectified sine pulses to the fluid. The obtained solutions by means of the Laplace and Fourier cosine and sine transforms are presented as a sum of the corresponding Newtonian and non-Newtonian contributions. The effects of the magnetic field, permeability, and the period of the oscillation have been observed on the fluid motion. Moreover, the influence of the side walls on the fluid motion and the distance between the walls for which the velocity of the fluid in the middle of the channel is negligible are presented by graphical illustrations. MSC: 76A05; 76A10

Keywords: generalized Burgers fluid; rectified sine pulses shear stress; Fourier cosine and sine transforms

ringer

1 Introduction

Motion of non-Newtonian fluids on oscillating plates is not only of fundamental theoretical interest but it also occurs in many applied problems, e.g., clay rotation, heart pumping, artificial surfing etc. Erdogan [1] obtained a solution as a sum of steady and transient solutions for the flow of a viscous fluid produced by a plane boundary moving in its own plane with a sinusoidal variation of velocity. Exact solutions for unsteady flow of a generalized Burgers fluid due to a rigid plate between two infinite parallel plates, one of which is an oscillating and time-periodic plane Poiseuille flow, was established by Fetecau etal. [2]. Zheng et al. [3] established an exact solution for the unsteady flow of a generalized Maxwell fluid over a flat plate. The plate was set into oscillating motion induced by hyperbolic sine velocity. Some recent work involving oscillating flows has been presented in many studies [4-7].

MHD flow of fluids and motion of fluids through porous media occur in medicine, engineering problems and geophysics, e.g., cardiology, delivery of medicine to affected areas, regulation of skin, nuclear reactors and geomagnetic dynamo. Khan and Zeeshan [8], and Ghosh and Sana [9] investigated the MHD flow of an Oldroyd-B fluid through a porous space. The motions were generated in the fluid due to the velocity sawtooth pulses of the plate. Ghosh and Sana [10] discussed the unsteady motion of an Oldroyd-

© 2014 Sultan et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the originalworkis properly credited.

B fluid in a channel bounded by two infinite rigid parallel plates in the presence of an external magnetic field acting normal to the plates. The flow is generated from rest due to rectified sine pulses applied periodically on the upper plate with the lower plate held fixed.

In the above citations, the conditions on the boundary are given in terms of velocity. The stress at the boundary gives important information as regards the nature of dissipation at the boundary. Little work is available in the literature where the oscillating stress is given on the boundary. Vieru et al. [11] analyzed the unsteady motion of a second grade fluid between two parallel side walls induced by oscillating shear stress. Li et al. [12] presented an analysis for helical flows of a heated generalized Oldroyd-B fluid subject to a linear time-dependent shear stress in a porous medium, where the motion is induced by the longitudinal time-dependent shear stress and the oscillating velocity at the boundary. Jamil et al. [13] and Shahid et al. [14] determined the starting solutions for the motion of Oldroyd-B fluids induced by quadratic, and cosine and sine oscillating time-dependent shear stress, respectively. Sohail et al. [15] presented closed-form expressions for the starting solutions corresponding to the unsteady motion of a Maxwell fluid due to an infinite plate that applies oscillating shear stresses to the fluid. Rubbab et al. [16] derived the unsteady natural convection flow of an incompressible viscous fluid near a vertical plate that applies a shear stress which is of exponential order of time.

In spite of all these citations and work in this direction, no attempt is made for oscillations induced by rectified sine pulses stress in a bounded domain. The main objective of the present investigation is to study the MHD oscillatory flow of a generalized Burgers fluid through a porous medium between two parallel walls. The formulation of the governing problem is made using the modified Darcy law of a generalized Burgers fluid. The induced magnetic field is assumed to be small as compared with the applied magnetic field. Analytical expressions for the velocity field and the shear stress are determined by means of the Fourier cosine and sine transforms coupled with Laplace transform. Finally, a comprehensive study of some physical parameters involved is performed to illustrate the influence of these parameters on the velocity.

2 Governing equations

For the generalized Burgers fluid, the Cauchy stress tensor is given by

where -pi denotes the indeterminate spherical stress, S is the extra-stress tensor, A = L + LT is the first Rivlin-Ericksen tensor (L being the velocity gradient), j denotes the upper convective derivative, ¡x is the dynamic viscosity, X1 and X3 (< X^ are the relaxation and retardation times, X2 and X4 are the material parameters of the generalized Burgers fluid having the dimension of the square of time, and

We seek the velocity field V and the stress field S of the form

V = V(x, y, t) = w(x, y, t)k, S = S(x, y, t), (3)

where k is the unit vector along the z-direction. If the fluid is at rest up to the moment t = 0, then

V(x, y, 0) = 0, S(x, y, 0) = d S(^y,0) = 0. (4)

Equations (1), (2), and (4) give the trivial stresses Sxx = Sxy = Syy = 0 and the meaningful equations

/ d d2\ , x / d d2\ dw(x, y, t)

l1 + X1 dt + X2 Jf)T1(x,y, t) = H3 dt + X4 It?)-—aT~

for n(x, y, 0) = 0, (5)

.992\ / d d2\ d w(x,y, t)

+ X1 Ft + X2 d?)T2(x,y, t) = H1 + X3 dt + X4 ^^^^

for T2 (x, y, 0) = 0, (6)

where r1 = Sxz(x,y, t) and t2 = Syz(x,y, t) are the non-zero shear stresses. The Darcy resistance R in a generalized Burgers fluid satisfies the following expression:

/ 3 32 \ ¡$( 3 32 \

C1+X1 l^t + X2^) R = ~~\ 1+X3 ~3~t + X4d?)V (7)

where $ is the porosity and k is the permeability of the medium.

We assume that a uniform magnetic field of strength fio is applied to the fluid. We also assume that the direction of the magnetic field is perpendicular to the velocity field. Thus the Lorentz force due to the magnetic field becomes

J x B = -aB2V, (8)

where a is the electrical conductivity of the fluid.

The balance of linear momentum which governs the MHD flow through the porous medium becomes

p dV = V-t + J x B + R, (9)

here p is the density.

We consider the unsteady flow of an incompressible generalized Burgers fluid over an infinite flat plate between two parallel side walls separated by a distance a, perpendicular to the plate. At time t =0, the plate and the fluid are at rest. At time t = 0+, the plate applies a pulsating shear to the fluid induced by rectified sine pulses.

In view of Eqs. (5)-(9), the governing equation leads to

/ d d2\ dw(x, y, t) 1 + A.1— + ^2 —T —

V 1 d t d t2) d t

d d2 \ / d2 d2 \ , = v| 1 + À3 Ft + d^A dX2 + df) w(xy,t)

d d 2 \ id d 2 - ^ 1 + Ft + X2 dtp) W(x, y, t)-^ 1 + ^3 Ft + ^ ItP )W(X, y, t), (1°)

where v = ^ is the kinematic viscosity, ^ = is the magnetic parameter, and e = k is the porosity parameter. We use the following appropriate initial conditions:

d w(x, y,°) d2 w(x, y,°)

w(x,y,°) =---=-—=-= °, forx >°andy e [°, d], (11)

the boundary conditions

( d d2\ , ( d d2\ dw(x, y, t)

[1 + hFt + "2 Jfi)T1(x,y,t)u=° = H1 + "3dt + Jfi)^^ x=° = Uf(t)

for y e (°, d) and t > °,

w(x, °, t) = w(x, d, t) = ° for x, t > °,

and the natural conditions

. . d w(x, y, t)

w(x,y, t) =--► ° as x ^œ,y e [°, d], t > °.

(12) (13)

According to the nature of the applied stress, we assume that the mathematical form of the function/(t) is [8]

/(t) = sin(TtW) + 2;£(-1)psin(T(t-pT)W(t), p > 0 and T >0, (15)

\ / p=i \ /

where H(■) is the Heaviside unit step function of period T and is defined as

HpT (t) = 0 fort < pT and HpT (t) = 1 fort > pT.

In order to solve the problem, we use the Laplace transform technique and Fourier cosine and sine transforms in this order.

3 Calculation of velocity field

Applying the Laplace transform to Eq. (10), we obtain the following problem: (1 + k\q + k2q2) qw (x, y, q)

, 2\/ 'd2 w (x, y, q) d 2w (x, y, q)

= v( 1 + X3q + ^q2) I -—-+

dx2 dy2

- ^(1 + k\q + k2q2)w(x,y,q) - e( 1 + X3q + k4q2)w(x,y,q).

The Laplace transform w(x, y, q) of the function w(x, y, t) has to satisfy the conditions

(l + Xiq + X2q2) ti(x, y, q)|x=o = [i(i l + X3q + X4q2)

^ d w (x, y, q)

= U—| l + 2 ^(-1)p exp(-pTq) j for y e [0, d],

q2 + —

q + T2

w(x, 0, q) = w(x, d, q) = 0 for x >0,

. d w (x, y, q) w(x,y, q) =--► 0 as x ^œ,y e (0, d).

(18) (19)

Multiplying both sides of Eq. (16) by ^cos(M x) sin(Xny), where Xn = d, integrating with respect to x and y from 0 to œ and 0 to d, respectively, and bearing in mind the conditions (17)-(19), we find that

wvn(H, q) =

2 U ((-1)n -1)

n p Xn

x 1/( (X2q3 + (X1 + vX4(M2 + X2n) + X2a + X4c) q2 + (1 + vX3(f2 + X2n) + X^ + X3f)q + v(M2 + X2n) + fi + c))

q2 + (T )2)

^1 + ^(-1)P exP(-pTq )j,

_ /2 /"œ fd_ wn(M, q) = J— / I w (x, y, q) cos(M x) sin(Xny) dydx; n = 1,2,3,. V n J0 J0

Equation (20) can also be written as

_ /2 U ((-1)n -1)

wvn(H, q) = J-

n p Xn -1 (

q + q4,n(f )

x a (q1 ,n (M ) q2,n (M ))(q1,n(M ) q3,n (M ))(q1,n(M ) + q4,n (M))(q q1,n (M ))

(q2,n(M ) -q1,n(M ))(q2,n(M ) - q3,n(M ))(q2,n(M ) + q4,n(M ))(q - q2,n(M ))

<fen(M )

(q3,n(M ) -q1,n(M ))(q3,n(M ) - q2,n (M ))(q3,n (M ) + q4,n (M ))(q - q3,n(M ))

<&4,n(M )

(q1,n(M ) + q4,n(M ))(q2,n(M ) + q4,n(M ))(q3,n(M ) + q4,n(M ))(q + q4,n (M ))

q^(f ) = v (f 2 + xI) + a + c,

Xi + vk4(f 2+ X2n) + k2tt +

qi,n(f ) = Mf )--—2-,

4>Uf) = k2q\n (f) + X + vX4(f2 + X2n) + + A.4c) q2n(f) (23)

+ (vX3(f 2+ X2n) + ki& + qi,n(f), i = 1,2,3, and 04,n (f) = -A.2 q4,n(f) + X + vk4(f2 + xn) + + k4c)q\n(f)

- (vX3(f 2+ Xl) + Xia + X3c)qUn(f). In the above relations

( hn(f) | (Pl,n(f ))2 ("i(f jFV/3

5i,n(f H-—+ V-^"ST"

+ ( - ^nM - y(^i,n(f ))2 + Kn"(iFy/3, (24)

A,n(f) I (A,n(f))2 (ai,n(IF\i/3

s2n(f ) = z--+ w-+-

2,nvs/ > 2 V 4 27

+ .^-.i^ + ^2F)1/3, (25)

„2/ ^i,n(f) /QMf))2 Kn(fFV/3

«3,n(f )=Z2--— + \/--- + '

2 4 27

+z( - finm -1(^ + Y/3, (26)

... i + vX3(f 2 + xn) + xia + X3C (Xi + vX4(f 2+ X2n) + X2a + X4C)2

ai,n(f ) = -----TT"2-, (27)

X2 3X2

2 2 2 2 3

v(f2 + X2n) + a + c (i + vX3(f2 + x2n) + Xifi + X3C)

fiin(f ) = —---+ 2 3

X2 27X2

(Xi + vX4(f2 + x2n) + X2a + X4c)(i + vX3(f2 + x2n) + Xia + X3C)

3 X 22

n to n

-i + i43 , ,

z =—. (()

To solve Eq. (22), let us take

Hn(q) =12 n ^ T-T&\\ + ^(-i)^, 2 n ^ T-exp(-pTq), (30)

(q2 + (f )2)(q - qi,n(f)) p=i (q2 + (f )2)(q - qi,n(f))

Gn(q) =

(q2 + (T )2)(q - qi,n(H))'

we can prove that

exp(q1,n(f )t)

£n(t)=£-1( Gn(q)) = n T

n 2 + T 2q2,n (f)

n 2 + T2q2,n (f) T2q1,n(f) sini T ) + n T cos

The inverse Laplace transform of Eq. (30), using Eq. (32), is given by

Hnf, t) = n T H (t)

n 2 + T 2q2(f) 1

n 2 + T2 q2n(f)

T 2q1(f) sin( ^ + n T cos( nT^^H(t) + 2£(-1)p \ / \ / /

exp(q1,n(f )(t - pT))__1

n 2 + T2q2n(f) n 2 + T 2q2,n (f)

T 2q1,n(f) sin

(n (t - pT)

+ n T cos

n (t - pT)

HpT (t).

Inversion of Eq. (22) by means of the Laplace formula and Fourier cosine and sine transforms and using Eq. (33) result in

2 ^ Pi C ^

w(x, y, t) = a ^Mj Bnf) cos(f x)

\ exp[-q4,n(f )t] n ( n2 + T2q\,n(f)

x J"« + ^(-1)P ex^-q4,n(f )(t - PT)) Hpt (t)^

n 2 + T2 q4,n(f)

(Tq4,n(f )F1(t) - nF2(t))

2 ^ /2

sin(kny) An(f) cos(fx)

n T| ^1,n(f )

exp[q1,n(f )t] n 2 + T2q2,n(f)

x H(t) + 2^(-1)p exp((t - pT)q1,n (f ))HpT(t)

f 2,n (f (h (t) + 2 ¿(-1)p exp((t - pT )q2,n(f ))HpT «

n + T q2,n(f) \

+ f 3,n (f) 11 (h (t) + 2 £(-1)p exp((t - pT )q3,n (f ))HpT (t)

- f4,n(f )

n 2 + T 2q2n(f)'

exp[ q4,n (f)t]

n 2 + T 2q4,n(f)

x (^(t) + 2^(-1)pexP(-q4,n(M)(t-pT))HPT(t)jj

T-J l.\rj,2l q1,n (M )fi,n(M ) q2,n (M )^2,n(M )

1() U 2 + T2q2,n (M ) + n 2 + T 2q2,n(M )

q3,n (M )f%n(M ) q4,n (M )^4,n(M )

+ n 2 + T2 q2,n (M ) n 2 + T2q2,n(M )

- Î2(t)n T

^2,n(M )

n 2 + T 2q2,n(M ) n 2 + T2 q2,n (M )

^3,n(M )

^4,n(M )

n 2 + T2 q2,n (M ) n 2 + T2q2,n (M )

MM ) =

^2,n(M ) =

f3,n(M ) = ^4,n(M ) =

An(M )

(q1 ,n (M) q2,n (M ))(q1,n(M ) q3,n (M ))(q1,n(M ) + q4,n (M))

_fo,n(M )_

(q2 ,n (M)

q1,n (M ))(q2,n(M ) q3,n (M))(q2,n(M) + q4,n (M)) _fo,n(M )_

(q3 ,n (M)

q1,n (M))(q3,n(M) q2,n (M))(q3,n(M) + q4,n (M)) _^4,n(M )_

(q1,n (M ) + q4,n(M ))(q2,n(M ) + q4,n(M ))(q3,n(M ) + q4,n(M )),

ii(t) = sin( yV (t) + 2^(- 1)p sin( n (t -pT)\ HpT (t),

\ / p=1 \ /

B,№ ,= /! u mi->,

\ n p Xn ( n t

Î2(t)=cosf ^V (t) + 2^(- 1)p cos( n (t-pT )\HpT (t),

\ / p=1 \ /

A„(M )=/! U M.

V n p XnX2

The first part of Eq. (34) gives the corresponding solution for a Newtonian fluid, while the second part gives the corresponding non-Newtonian contribution. The transient part of velocity for Re(q1,n(M )), Re(q2,n(M )), Re(q3,n(M )) < 0, Re(q4,n(M )) > 0 is

wt (x, y, t)

2 œ F2 f œ

= sin(Xny) Bn(M ) cos(Mx)

d Z^ V n J0

n 2 + T2 q4,n(M )

H (t) + 2^( - 1)pHpT (t)

2 2 œ sin(Xny) V n J0

An(M ) cos(M x)

r, fi,n(f) f2,n(f)

n T{——_„ n ,,, +

n 2 + T2qln(f) n 2 + T 2q2,n(f)

^4,n(f )

n 2 + T 2q2,n(f) n 2 + T2q4,n(f) while the steady state part is given by

2 F2 f to

Ws(x, y, t) = ^Y] \ 2 sin(Xny) Bn(f) cos(f x)

H (t) + 2j2(-i)PHpT (t) p=i )

df, (38)

n 2 + T 2q2,n(f)

(Tq 4,n(f )Fi(t)-n F2(t))

2 f° VJ- sin(Xny) .. 1 ' n ./o

An(f ) C0s(f x)

n=i ^i,n(f )

n 2 + T 2q2,n(f)

^2,n(f )

(TFi(t)qin(f ) + n F2(t))

n 2 + T 2q2,n (f)

^3,n(f )

n 2 + T 2q2,n (f)

(TFi(t)q2(f ) + n F2(t)) (TFi(t)q3,n(f ) + n F2(t))

^4,n(f2\j,, (TFi(t)q4,n(f)-nF2(t))

n 2 + T 2q4,n (f)

4 Calculation of tangential stresses

To obtain the expressions for the shear stresses ri(x,y, t) and r2(x,y, t), applying the Laplace transform to Eqs. (5) and (6), we have the expressions

ti(x, y, q) = ix T2(x, y, q) = —

(i + X3q + X4q2) dw(x,y, q) (i + Xiq + X2q2) d x

(i + X3q + X4q2) dw(x,y, q)

(i + Xiq + X2q2) dy From Eq. (20), with inverse Fourier cosine and sine transforms, we have

2 -TO /2 f

w(x,y,q) = -> J—An(f)sin(Xny) cosfx

d vn io

(q - qi,n (f ))(q - q2,n (f ))(q - q3,n (f))

q2 + (T )2)

f to \

1 + 2^(-1)p exp(-pTq) df. V p=i /

Using Eq. (4i) in the set of Eq. (40), we have

2— /2 Ti(x, y, q) = —— V A -An(f) sin(Xny) ^ , V n

f sin f x(i + X3q + X4q2)

o (i + Xiq + X2q2)(q - qi,n (f))(q - q2,n (f))(q - q3,n (f))

n TO \

y + ^(-i)p exP(-pTq)J df, (42)

T2(x, y, q) = 21 V J - XnAn(f) cos(Xny) d V n

cos fx(i + X3q + X4q2)

o (i + X3q + X4q2)(q - qi,n (f))(q - q2,n (f))(q - q3,n (f))

n / TO \

i + 2^(-i)p exp(-pTq) df. (43)

(q2 + (T )2K p=i /

Let us take

A (q) = i + X3^ + X^, (4)

i + Xiq + X2q2 which can also be written in the form

- , , q + ai bi

A (q)=fl2+ a.3---—-2+ a.4-.-^—(5)

(q + ai)2 - b2 (q + ai)2 - b2

Xi X4 X2X3 - XiX4

ai = -, «2 = -, 0-3 =

2X2 X2 X

t(M, N, t) = (A * B)(t) = i A(t - q)B(M,N, q) dq. Jo

2X2(X2 - X4) - Xi(X2X3 - XiX4) JX^- 4X2 «4 = -, , , -, bi = ---,

X2V X2-4X2 2X2

where Xj2 - 4X2 > o. Applying the inverse Laplace transform to Eq. (45), we obtain

A(t) = a2 + a3 cosh(bit) exp(-ait) + a4 sinh (bit) exp(-ait). (47)

Employing the methodology as for the velocity field, the inverse Laplace transform of Eqs. (42) and (43) results in

2— [2 fm Ct

Ti(x,y, t) = —— .An(f) sin(Xny) I f sin fxl («2 + («3cosh(bi(t- q))

d V n Jo Jo

^ , V n

+ a4 sinh (bi(t - q))) exp(-ai(t - q)))

n ^ "i;(f ^^ fa) + 21 (-1)p exp(-pTqi,n(f ))HpT(q)^

(f) exp[q2,n(f )q] + n 2 + T 2q2,n(f)

H (q) + 2 £](-i)p exp(-prq2,„ (f ))HpT (q)

(f) exp[q3,«(f )q]

H (q) + 2 £](-i)p exp(-pTq 3,n(f H (q)

* 2 + T 2q2,n(f) y p=i

-Fi(q)T2| ——„0 „ + ——„„ „ +

* 2 + T 2q2,n(f) * 2 + T2q2n(f) * 2 + T 2q2(f)

- ^2(q)* T

ni,n(f)

%,n(f )

* 2 + T 2q2,n (f) * 2 + T 2q2,n(f)

n3,n(f)

* 2 + T 2q2,n(f)

<iq<if,

2^ /2 ^to

r2(x,y,t) =— —X.nAn(f)cos(kny) I cosfxl (a2 + (a3cosh^^t- q))

« n=A * Jo Jo

+ fl4 sinh {bi(t - q))) exp {-ai(t - q))) (f) exp[qi,n(f )q]

*2 + t 2q2,n (f)

x H (q) + 2 £](-i)p exp( -pTqi,n(f )) HpT (q)

+ n2*f )eXT2[q22,fq^H(q) + 2 ^(-i)p exp(-pTq,n(f ))HpT(q)

* + T q2,n(f) \ p=i )

+ Wf )eTP2[q23,nf)qW H (q) + 2 £(-i)p exp(-pTq 3,n(f ))HpT (q)

* + T q3,n(f) \

^2/ qi,n(f)ni,n(f) q2,n(f)n2,n(f) q3(f)n3,n(f) -Fi(q)T2| ——„„ „ + ——„„ „ +

* 2 + T 2q2,n(f) * 2 + T2q2n(f) * 2 + T 2q§(f)

- F2(q)* T

n3,n(f)

ni,n(f)

n2,n(f )

* 2 + T 2q2,n(f)

* 2 + T2q2,n (f) * 2 + T 2q2,n(f)

<iq<if,

ni,n(f) = n2,n(f) = n3,n(f) =

(qi ,n (f ) q2,n (f ))(qi,n(f) q3,n (f)) i

(q2 ,n (f ) qi,n (f ))(q2,n(f) q3,n (f)) i

(q3 ,n (f) qi,n (f))(q3,n(f) q2,n (f ))

the shear stresses for the generalized Burgers fluid.

5 Limiting cases 5.1 Burgers fluid

Letting k4 = 0, we obtain the velocity field and the associated shear stresses corresponding to a Burgers fluid performing the same motion.

5.2 Oldroyd-B fluid

Letting k2 = X4 = 0 into Eq. (20), and following the same way as before, the velocity and tangential stresses expressions for Oldroyd-B fluid take the form

T [ exp[-g4,„(f )t]

2 C /2 fœ wo(x, y, t) = ^J-sin(kny) Bn(f ) cos(f x)

\n 2 + T 2qljf )

x H (t) + 2^(-l)p exp(-q4,„(f )(t - pT)) Hpt (t)

V p=i )

n 2 + T 2q\n (f )

(Tq 4,„(f )Fi(t)-n Fz(t))

2 W2 r

-sm(kny) V n Jo

Bn(f ) cos(f x)

n T| f5,„(f )

exp[q5,«(f )t]

n 2 + T 2q5,„(f ) x ( H(t) + 2^(-l)pexp((t-pT)qSn(f))HpT(t)

+ ^6,«(f ) - Mf )

exp[q6,n(f )t] n2 + T2q2„(f )

exp[ q4,n (f )t] n 2 + T 2q4,n(f )

H (t) + 2 £](-l)p exp( (t - pT)q6,n (f )) HpT (t)

V p=i y

x ^H (t) + 2 E(-l)P ex^-q4,n(f )(t - pT)) HpT (t)j J

2 q5,n (f )Mf ) q6,n (f )Kn(f ) q4,n (f )^7,n(f )

- Fl(t)T

- F2(t)n T

n 2 + T2q5,n(f ) n 2 + T 2q6,n(f ) n 2 + T2q4,n(f )

^5,n(f ) n2 + T2q5,n(f )

^6,n(f )

^7,n(f )

n2 + T2q2n (f ) n2 + T2q2n(f )

2^ /2 /"C

tlo(x,y, t) = —An(f ) sin(kny) f sinfx d in Jo

(f ) exp[q5,n(f )q] n 2 + T2q2n(f )

H (q) + 2 £](-l)p exp(-pTq5,n(f ))HpT (q)

(M ) exp[q6,n(M )q]

n 2 + T 2q6,n (M )

(H(q) + 2£(-1)pexp(-pTq6,n(M))HpT(q) )

n4,n(M ) exp[-£■]

n2 + ( XT1 )2

' H(q) + 2 X>1)p ex^pTjHpT(q)^

r/ ^2/ q5,n(M)n5,n(M) q6,n(M)n6,n(M) 1 n4,n(M)

- FXq)T21 ——„„ „ +

n 2 + T 2q2,n (M ) n 2 + T 2q6,n(M ) X1 n 2+ (T.)2

- F2(q)n T

V5,n(M )

n6,n(M )

_ _ n4,n(M )

n2 + T2q2,n(M ) + n 2 + T 2q6,n (M ) + n2 + (^)2

dM, (52)

2^ 2 fœ

T20(x,y, t) = — /— XnAn(M) cos(Xny) cos Mx

d n=i V n ^0

(M ) exp[q5,n(M )q] n 2 + T 2q2,n(M )

x ^H(q) + ^(-1)p ex^-pTq5,n(M))HpT(q)j

+ n6n|+xpq|M /h (q) + 2 £ (-1)p exp(-pTq 6,n(M ))HpT (q))

n + T q6,n(M ) \ p=1 )

+ iinM+lil (hO + 2 £ (-1)p exp^ )h,t (q))]

w ^2/ q5,n(M)n5,n(M) q6,n(M)n6,n(M) 1 n4,n(M)

- F1(q) T21 ——„„ „ +

n 2 + T 2q2,n (M ) n 2 + T 2q6,n(M ) X1 n 2+ (T1)2

- F2(q)n T

n5,n(M )

n6,n(M )

_ _ n4,n(M )

n2 + T2q2,n(M ) + n 2 + T 2q6,n (M ) + n2 + (T1)2

dM, (53)

X1q2,n(M ) + (vX3(M2 + X2n) + X1fi + X3C )q5,n(M )

^5,n(M ) =

(q5 ,n (M) q6,n (M ))(q5,n(M ) + q4,n (M))

, ^ X1q2n(M ) + (vX3(M 2+ Xn) + X1fi + X3c)q6,n(M ) ,C/I,

^6,n(M ) = -n-777-¡7777-7-7-7-77-, (54)

(q6,n(M) - q5,n(M))(q6,n(M) + q4,n(M)) (M) = X1q2,n(M) - (vX3(M2 + X2n) + X1fi + X3c)q4,n(M) 7,n (q4,n(M ) + q5,n(M ))(q4,n(M ) + q6,n(M )) ,

q5,n(M ), q6,n(M )

= (-(1 + vX3(M 2+ X2n ) + X1fi + X3C)

(1 + vX3(M2 + xn) + X1fi + X3C)2 -4X1(vX3(M2 + xn) + X1fi + X3C^/(2X1),

^(f) = _k1- k3_,

4,n (1 + Xiq5,„(f))(1 + Xiq6,„(f)),

n5,n(f) = l6,n(f) =

1 + ^3q5,«(f)

(1 + Xiq5,„(f ))(q5,n(f) q6,n (f))

_1 + ^3q6,«(f)_

(1 + Xiq6,„(f ))(q6,n(f) q5,n (f))'

5.3 Maxwell fluid

Letting k3 = 0 in the set of equations for Oldroyd-B fluid, we obtain the corresponding expressions for the Maxwell fluid:

\ exp[-q4,„(f)t] n ( n 2 + T2 ql n (f)

2 C /2

(x,y, t) = - V J — sin(kny) Bn(f) cos(f x) d in Jo

x tf (t) + 2^(-1)p exp(-q4,„(f )(t - pT)) tfr (t)

V p=1 )

n 2 + T 2q\,n(f)

(Tq4,n(f )Fi(t) - ni2(t))

2 C /2

sin(kny) Bn(f) cos(f x)

d n=1 V n •'o

n T| f8,n(f )

exp[q7,n(f )t] n 2 + T2 q2n(f)

x H (t) + 2^(-1)p exp((t - pT )q7,n(f (t)

+ Mf) e2xp[Sf (2f 11 (H (t) + 2f(-1)p exp((t - pT )q8,n(f ))tfpT (t)

n + T q8,n(f) \

- ^10,n(f )

exp[ q4,n (f )t]

n 2 + T2 q2,n(f)

x [H (t) + 2^(-1)p exp(-q4,n (f )(t - pT)) HpT (t)

V p=1 )

- F1(t)T

al q7,n(f)f8,n(f) , q&.n(f)fg,n(f) q4,n(f)^1o,n(f)

n 2 + T2 q7,n(f) n 2 + T 2q8,n(f) n 2 + T2q4,n(f)

J-. /,\ ^8,n(f) ^g,n(f)

-F2(t)nT ——_ „ ,., +

n 2 + T2q2n(f) n 2 + T 2q2,n(f)

^10,n(f )

n 2 + T 2q4,n(f)

2^ /2 /"C

tlm (x, y, t) = —2An (f) sin(kny) f sin f x d in Jo

(f) exp[q7,n(f )q]

n 2 + T 2q2,n (f)

x (H(q) + 2 £](-1)p exp(-pTq7,n(f))HpT(q)

+ n9;f+eTp^)q] H(q) + 2£(-1)p exp(-pTq8,n(M ))HpT (q)

n7,n(M ) exp[-£]

■ n 2 + (t )2 ^H(q) + ^(-1)pex^pX^)HpT(q)

N-^/ q7,n(M)n8,n(M) q8,n(M)n9,n(M) 1 ^7,n(M) -i1(q)T2| —-„„ „ + ■

n2 + T2q2,n(M) ' n2 + T2q2,n(M) X1 n2 + (£)2

- f2(q)n t

n9,n(M )

n8,n(M )

_ _ l7,n(M )

n2 + T2q2,n (M ) + n 2 + T 2q8,n (M ) + n2 + (^)2

dM, (59)

T2M (x,.

2^ /2 fœ

', y, t) = 22 XnAn(M ) cos(Xny) / cos M x

d vn J0

(M ) exp[q7,n(M )q]

n 2 + T 2q2,n (M )

x H (q) + 2 £](-1)p exp( -pTq7,n(M )) HpT (q)

V p=1 )

+ ^(M ^f^f H(q) + 2 £(-1)p exp(-pTq8,n(M ))HpT (q) n + T q8,n(M ) \ p=1 y

n7,n(M ) exp[-X^]

n2 + (XT )2

H (q) + ^(-1)p exp( pT)HpT (q)

N^2/ q7,n(M)n8,n(M) q8,n(M)n9,n(M) 1 l7,n(M) - i1(q)T2| —-„„ „ +

n 2 + T 2q2,n (M ) n 2 + T 2q2,n(M ) X1 n 2+ (X1)2

- F2(q)n T

n9,n(M )

n8,n(M )

_ _ l7,n(M )

n2 + T2q2,n (M ) + n 2 + T 2q2,n (M ) + n2 + (^)2

dM, (60)

^8,n(M ) =

f9,n(M ) = ; fw,n(M, ) =

X1q2n (M ) + X1fiq7,n(M )

(q7 ,n (M) q8,n (M))(q7,n(M) + q4,n (M))

X1q8,n(M ) + X1fiq8,n(M )

(q8,n (M) q7,n (M ))(q8,n(M ) + q4,n (M))

q7,n(M ^ q8,n(M ) = n7,n(M ) = n8,n(M ) = n9,n(M ) =

X1q4,n(M) - X1^q4,n (M)

(M) + q7,n (M ))(q4,n(M ) + q8,n(M ))' -(1 + X1fi) (1 + X1fi)2-4X2fi

(1 + X1q7,n(M ))(1 + X1q8,n(M ))' 1

(1 + X1q7,n(M ))(q7,n (M ) q8,n (M))

(1 + X1q8,n (M ))(q8,n(M ) q7,n (M))

Figure 1 Velocity profile for generalized Burgers fluid for different values of x and d. Other parameters and values are taken as p = 975, x = 0.01, y = 2, f = 3.9, T = f, U =15, Xi =3, X2 =2, X3 = 0.5, X4 = 2, € = 1.7, and ft = 3.

Figure 2 Velocity profile for a generalized Burgers fluid for different values of x and y. Other parameters and values are taken as p = 975, d = 0.2, f = 3.9, T = 4, U =15, Xi =3, X2 = 2, X3 =0.5, X4 = 2, € = 0.3, and ft = 0.5.

Figure 3 Diagrams of the velocity for different values of the parameter ft. Comparison between five models.

6 Results and discussion

The present problem is concerned with unsteady motion of the generalized Burgers fluid generated from rest induced by rectified sine pluses shear stress. The Laplace transform technique and Fourier cosine and sine transforms have been used as mathematical tools in this order. The obtained expression for the velocity field has been written as the sum of Newtonian and non-Newtonian contributions.

By using the numerical calculations and graphical illustrations, the following physical aspects of the fluid behavior have been analyzed: (a) Influence of side walls on the velocity field.

In order to study the influence of the side walls on the fluid velocity we prepared the diagrams contained in Figures 1 and 2. These diagrams present the velocity field w(x, y, t) for three values of the distance x at the bottom plate. It is seen that if d = 0.1 (small distance between walls), the amplitude of oscillation of velocity in the middle of the channel is

ol_I_!_I_!_

0 2 4 6 S 10

2i-1-1-

Figure 4 Diagrams of the velocity for different values of the parameter €. Comparison between five models.

small. By increasing the value of d, the amplitude of oscillation of velocity increases. But it exists up to a value of 0.5 (in our case) after which the velocity remains the same (velocity profiles are almost identical) if d increases. Therefore, in the considered case, after d = 0.5 the influence of the side walls on the velocity is negligible. In Figure 2 have been sketched the velocity profiles for different positions in the channel, starting from the side walls till the middle of the channel, i.e. for different values of the variable y, for a fixed distance between the side walls. The amplitude of oscillations decreases far from the plate. It is also observed that, if the bottom plate is set into oscillation, the velocity increases with respect to the y-coordinate, from zero to a maximum in the middle of the channel. (b) Comparative study of various models.

Figure 5 Diagrams of the velocity for different values of the parameter T. Comparison between five models.

Figures 3, 4, and 5 are sketched in order to compare the various fluid types. Also, the influence of the magnetic field, permeability and the period of oscillation of sawtooth pulses on the velocity field can be observed from these diagrams. In these figures, we have considered the following values of the parameters: U = 15, X1 = 2.9, X2 = 4, X3 = 0.35, X4 = 1.5, x = 0.05, d = 0.1, y = d/2, v = 0.004, ^ = 3.9.

Figure 3 shows the diagrams of velocity w(x, y, t), versus t, for the porosity parameter 6 = 0.1, the pulse period T = 4, and for three values of the magnetic parameter, Q = 0.2,1.7,3.2. It is observed that there is a time interval in which the velocity has an oscillating behavior for all kinds of fluids. After this moment, the velocity tends to a common value (the differences between velocities of different fluids are insignificant). For low values of the magnetic field strength, the amplitudes of the velocity oscillations are smaller for the generalized Burgers, Oldroyd-B, and Newtonian fluids, and much larger for Maxwell and Burgers fluids (see diagrams for Q = 0.2). If the magnetic field is stronger, the ve-

Figure 6 The influence of the parameter X1 on the velocity field for x = 0.25, d = 0.50, y = 0.25 and

X2 = 2, X3 = 0.5, X4 = 2.

locity amplitudes of Maxwell and Burgers fluids decrease while, the velocity amplitudes of generalized Burgers and Newtonian fluids increase (case Q = 3.2). It is important to note that velocities of the fluids tend to a common value in shorter time if the magnetic field is stronger. Figure 4 is plotted for Q = 0.2, T = 2 and for three values of the porosity parameter e = 0.1,0.5,0.g. In this case, amplitudes of the velocity oscillations of the generalized Burgers, Oldroyd-B, and Newtonian fluids are lower than those corresponding to the Burgers and Maxwell fluids.

Increasing permeability leads to increase the velocity amplitudes. The effect of the pulse period T on the velocity field is shown in Figure 5 for Q = 0.4, e = 0.1 and three values of the parameter T = 0.75,1.25,1.75. The Burgers fluid oscillates with larger amplitude, while other fluids have oscillations with amplitude close as order of magnitude. For low values of the parameter T, velocities tend to a common value in a shorter time than in the case of large values of the parameter T.

(c) Influence of parameters k1, k2, k3, k4 on the velocity field.

In order to study the behavior of the fluid for various values of the material parameters k1, k2, k3, k4, the diagrams of the velocity field from Figures 6, 7, 8 and g were plotted. In these figureswe usedthe numericalvalues U = 25, x = 0.25, d = 0.5, y = 2, T = 8, P = 830, fi = 3.78, e = 2, and three values for the magnetic parameter, namely Q = 0.05,0.5,1.5. Therefore, the influence of the magnetic parameter on the velocity field is also analyzed and, similar conclusions to those from the case b) were obtained. In Figure 6 the param-

^lxio 3 jxf.

£ 5xl0"4 3 0 -5xl0"4 - lxlO"3

Figure 7 The influence of the parameter X2 on the velocity field for x = 0.25, d = 0.50, y = 0.25 and X1 =2,^3 =0.5,^4 = 2.

eter X1 is variable and parameters X2, X3, X4 are constant. It can be seen that, if the values of the parameter X1 increase, the fluid flows more slowly. Also, be noted that for the same value of the parameter X1, the increasing of magnetic parameter values result in decreasing velocity of fluid flow (the velocity amplitudes decrease if the values of the magnetic parameter increase). Figure 7 corresponds to the variation of the parameter X2. Unlike the previous case when X1 is variable, in this case the velocity amplitudes are higher in the early period of the movement. Another difference appears in the behavior of the fluid, namely, velocity amplitudes increase if the parameter X2 increases. In this case, for a constant value of the parameter X2, the increasing of the magnetic field strength leads to attenuated fluid motion. The diagrams of Figure 8 correspond to the variable parameter X3. It is clear that the fluid behavior is similar to the case of the variation of the parameter X1. For increasing values of the parameter X3, the fluid velocity amplitudes decrease and, also, if X3 remains constant and the values of the magnetic parameter increase, then the fluid flows more slowly. The influence of the parameter X4 on the fluid velocity is shown in Figure 9. The fluid behavior is similar like in the cases of variations of parameters X1 and X3.

7 Conclusions

The purpose of this work is to provide exact solutions for the velocity field as well as shear stresses corresponding to the oscillating flow of a generalized Burgers fluid between two parallel side walls over a plate. The oscillation is induced by rectified sine

Figure 8 The influence of the parameter X3 on the velocity field for x = 0.25, d = 0.50, y = 0.25 and

X1 =2, X2 = 2, X4 = 2.

pulses shear stress applied to the bottom plane. These solutions, presented as a sum of the Newtonian and non-Newtonian contributions, are obtained by using Fourier cosine and sine transforms, and the Laplace transform. The main findings are summarized as follows:

• The amplitude of oscillation of velocity in the middle of the channel is small. With the increase of distance between the walls, the amplitude of oscillation of velocity increases.

• Amplitude of oscillations decreases far from the plate. As the bottom plate is set into oscillation, with respect to y-coordinate, the velocity increases from zero to a maximum in the middle of the channel.

• Velocities of all the fluids tend to a common value in shorter time if the magnetic field is stronger.

• Increasing permeability leads to the increasing magnitudes of amplitudes of velocity.

• For low values of the time period T, velocities tend to a common value in a shorter time as compared to large values of the parameter T.

• As the values of the parameters k1, k3, and k4 increase, the fluid flows more slowly whereas behavior of k2 is opposite.

Competing interests

The authors declare that they have no competing interests. Authors' contributions

The authors' contribution for this article is as follows: QS 40%, MN 30%, MI 15% and UA 15%. All authors read and approved the final manuscript.

Author details

1 Centre for Advanced Studies in Pure and Applied Mathematics, Bahauddin Zakariya University, Multan, Pakistan. 2Department of Mathematics, GC University, Faisalabad, Pakistan.

Acknowledgements

The authors would like to express their gratitude to the reviewers for fruitful remarks and suggestions. Also we are grateful to the Higher Education Commission (HEC) of Pakistan for supporting this research work.

Received: 19 February 2014 Accepted: 5 June 2014 Published online: 23 September 2014

References

1. Erdogan, ME: A note on an unsteady flow of a viscous fluid due to an oscillating plane wall. Int. J. Non-Linear Mech. 35,1-6 (2000)

2. Fetecau, C, Hayat,T, Fetecau, C: Steady-state solutions for some simple flows of generalized Burgers fluids. Int. J. Non-Linear Mech. 41, 880-887 (2006)

3. Zheng, L, Wang, KN, Gao, YT: Unsteady flow and heat transfer of a generalized Maxwell fluid due to a hyperbolic sine accelerating plate. Comput. Math. Appl. 61,2209-2212 (2011)

4. Fetecau, C, Hayat,T, Khan, M, Fetecau, C: A note on longitudinal oscillations of a generalized Burgers fluid in cylindrical domains. J. Non-Newton. Fluid Mech. 165,350-351 (2010)

5. Fetecau, C, Hayat,T, Fetecau, C: Starting solutions for oscillating motions of Oldroyd-B fluids. J. Non-Newton. Fluid Mech. 153,191-201 (2008)

6. Anjum, A, Ayub, M, Khan, M: Starting solutions for oscillating motions of an Oldroyd-B fluid over a plane wall. Commun. Nonlinear Sci. Numer. Simul. 17,472-482 (2012)

7. Bose, D, Basu, U: Unsteady incompressible flow of a generalised Oldroyed-B fluid between two infinite parallel plates. World J. Mech. 3,146-151 (2013)

8. Khan, M, Zeeshan: MHD flow of an Oldroyd-B fluid through a porous space induced by sawtooth pulses. Chin. Phys. Lett. 28,84701-84704(2011)

9. Ghosh, AK, Sana, P: On hydromagnetic flow of an Oldroyd-B fluid near a pulsating plate. Acta Astronaut. 64, 272-280 (2009)

10. Ghosh, AK, Sana, P: On hydromagnetic channelflow of an Oldroyd-B fluid induced by rectified sine pulses. Comput. Appl. Math. 28, 365-395 (2009)

11. Vieru, D, Fetecau, C, Rana, M: Starting solutions for the flow of second grade fluids in a rectangular channeldue to an oscillating shear stress. In: AlPConf. Proc.,vol. 1450, pp. 45-54 (2012). doi:10.1063/1.4724116

12. Li, C, Zheng, L, Zhang, Y, Ma, L, Zhang, X: Helicalflows of a heated generalized Oldroyd-B fluid subject to a time-dependent shear stress in porous medium. Commun. Nonlinear Sci. Numer. Simul. 17, 5026-5041 (2012)

13. Jamil, M, Fetecau, C, Imran, M: Unsteady helicalflows of Oldroyd-B fluids. Commun. Nonlinear Sci. Numer. Simul. 16, 1378-1386 (2011)

14. Shahid, N, Rana, M, Siddique, I: Exact solution for motion of an Oldroyd-B fluid over an infinite flat plate that applies an oscillating shear stress to the fluid. Bound. Value Probl. (2012). doi:10.1186/1687-2770-2012-48

15. Sohail, A, Vieru, D, Imran, MA: Influence of side walls on the oscillating motion of a Maxwellfluid over an infinite plate. Mechanika 19(3), 269-276 (2013)

16. Rubbab, Q, Vieru, D, Fetecau, C, Fetecau, C: Naturalconvection flow near a verticalplate that applies a shear stress to a viscous fluid. PLoSONE 8(11), e78352 (2013). doi:10.1371/journal.pone.0078352

doi:10.1186/s13661-014-0152-0

Cite this article as: Sultan et al.: Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress. BoundaryValue Problems 2014 2014:152.

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Immediate publication on acceptance

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com