# Flow of a Burger’s Fluid in a Channel Induced by Peristaltic Compliant Walls

Authors: I. Ahmad; N. Ali; A. Abbasi; W. Aziz; M. Hussain; M. Ahmad; M. Taj; Q. Zaman
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

A theoretical analysis is presented for the peristaltic motion of a magneto-hydrodynamic (MHD) non-Newtonian fluid in channel with complaint walls. The fluid obeys viscoelastic non-Newtonian model with Burger’s constitutive equation. The relevant equations are first developed and then solved using perturbation technique. Expressions of stream function and velocity components are constructed under the assumption thatδ (characteristic ratio of transversal and axial scales of peristaltic motion) is small. The results indicate the strong effects of Burger’s fluid parameter, Hartman number, Reynolds number, and complaint wall parameters on the velocity field and stream function. The obtained solutions are shown graphically for the different values of involved parameters.


---

## Body

## 1. Introduction

Peristaltic pumping is a mean of fluid transport caused by waves of contraction along the boundary of a channel/tube. To analyze the properties of peristalsis theoretically, the flow caused by sinusoidal traveling waves of contraction along the boundary has been considered in an infinite channel/tube. In general, flow in a peristaltic channel is driven by two independent mechanisms: the movement of the walls and due to the mean pressure gradient along the channel. After the pioneering work of Shapiro et al. [1] the fluid flow due to prescribed wall motion has been described theoretically by Ebaid [2]. The approach of Fung and Yih [3] was used to analyze the peristaltic flow by Hayat et al. [4] and Haroun [5]. In the above mentioned studies, peristaltic motion was investigated under various assumptions. One of the important aspects in such studies is the nature of the fluid. Since, it is a fact that most of the biological and industrial fluids are non-Newtonian in nature, therefore, peristaltic flow of these fluids has also been studied by Hayat et al. [6] considering non-Newtonian Maxwell fluid in a channel. Srinivas et al. [7] studied the MHD peristaltic motion with slip and heat characteristics. The effects of wall flexibility on the peristaltic motion have also attracted the attention of various researchers in the field. Mittra and Prasad [8] initiated the study of peristaltic transport of viscous fluid with compliant walls. Their work has been generalized by Muthu et al. [9] for viscous incompressible fluid in circular tube and Abd Elnaby and Haroun [10] for two-dimensional channel. The work in [10] has been extended for Maxwell fluid by Ali et al. [11]. Peristaltic motion of Jeffrey fluid in a channel having compliant walls was also studied by Hayat et al. [12]. Recently the heat and mass transfer effects on the peristaltic flow of an Oldroyd-B fluid in a channel with complaint walls were investigated by Hayat et al. [13]. Motivated by the above mentioned investigations, the aim of the present paper is to extend the work of [13] by considering the fluid which obeys the constitutive equation of Burger’s fluid model. The Burger’s fluid is a rate type viscoelastic fluid model [14] which has been found quite successful in predicting the behavior of viscoelastic material, food products such as cheese and soil. The viscoelastic behavior of fine-grained polycrystalline olivine can also be predicted by Burger’s fluid model [15]. Moreover, this model has also been used for calculation of transient creep proportion of the earth’s mantle [16] and in the interpretation of postglacial uplift [17].The paper is arranged in the following manner. Flow geometry is explained in Section2 while governing equations are developed in Section 3. Series solution is provided in Section 4. Results and discussion presented in Section 5. Some conclusions are made at the end of article in Section 6.


## 2. Flow Geometry

Let us consider the flow of magneto hydrodynamic (MHD) Burger’s fluid due to sinusoidal motion of the compliant walls in a channel. A constant magnetic field of strengthB
0 is applied in the direction normal to the flow. The deformation of the walls of the channel is represented by the expression
(1)
y
=
±
η
(
x
,
t
)
=
±
{
d
1
+
a
sin
2
π
λ
(
x
-
c
t
)
}
,
where a is the amplitude of the peristaltic wave, λ is the wavelength, t is the time, c is the speed of the wave, and d
1 is the half-width of the channel.


## 3. Governing Equations

Letu and v denote the velocity components in the axial and transverse direction, respectively, then we can write
(2)
V
=
[
u
(
x
,
y
,
t
)
,
v
(
x
,
y
,
t
)
,
0
]
.
The law of conservation of mass and momentum for the flow of an incompressible fluid are
(3)
div
⁡
V
=
0
,
ρ
d
V
d
t
=
-
∇
p
+
div
⁡
S
,
where V is the velocity, ρ is the density, p is the pressure, and S is the extra stress tensor, which satisfies the following relation for a Burger’s fluid:(4a)
S
+
λ
1
D
S
D
t
+
β
D
2
S
D
t
2
=
μ
(
1
+
λ
2
D
D
t
)
A
1
.
In the above equation μ is the dynamic viscosity, A
1 is the first Rivlin-Ericksen tensor, λ
1 and λ
2 are the relaxation and retardation times, respectively, β is Burger’s fluid parameter, L is the gradient of velocity, and D
/
D
t is the upper convective derivative and defined as
(4b)
D
S
D
t
=
∂
S
∂
t
+
(
V
·
∇
)
S
-
L
S
-
S
L
T
.Using the definition of velocity and D
/
D
t in (3) we get the following equations:
(5)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(6)
ρ
(
∂
∂
t
+
u
∂
∂
x
+
v
∂
∂
y
)
u
=
-
∂
p
∂
x
+
∂
S
x
x
∂
x
+
∂
S
x
y
∂
y
-
σ
B
0
2
u
,
(7)
ρ
(
∂
∂
t
+
u
∂
∂
x
+
v
∂
∂
y
)
v
=
-
∂
p
∂
y
+
∂
S
x
y
∂
x
+
∂
S
y
y
∂
y
,
(8)
ρ
d
d
t
(
∂
u
∂
y
+
∂
v
∂
x
)
=
(
∂
2
∂
y
2
-
∂
2
∂
x
2
)
S
x
y
+
∂
2
∂
x
∂
y
(
S
x
x
-
S
y
y
)
-
σ
B
0
2
∂
u
∂
y
,
where the components of extra stress appearing in the above equations satisfy three equations that can be obtained directly from (4a).The relevant boundary conditions are [8]
(9)
u
=
0
,
[
-
τ
∂
3
∂
x
3
+
m
∂
3
∂
x
∂
t
2
+
d
∂
3
∂
x
∂
t
]
η
=
∂
S
x
x
∂
x
+
∂
S
x
y
∂
y
-
ρ
d
u
d
t
-
σ
B
0
2
u
,
∂
S
x
x
∂
x
+
∂
S
x
y
∂
y
-
ρ
i
i
d
u
d
t
d
u
d
t
at
y
=
±
η
,
where τ is the tension in the wall, m is the mass per unit area, d is the coefficient of viscous damping, and σ is the electrical conductivity of fluid. Defining the stream function ψ such that u
=
∂
ψ
/
∂
y, v
=
-
∂
ψ
/
∂
x, and the following dimensionless variables:
(10)
ψ
*
=
ψ
c
d
1
,
x
*
=
x
λ
,
y
*
=
y
d
1
,
t
*
=
c
t
λ
,
η
*
=
η
d
1
,
p
*
=
d
1
2
p
c
λ
μ
,
S
i
j
*
=
d
1
S
i
j
c
μ
,
λ
1
*
=
λ
1
c
d
1
,
λ
2
*
=
λ
2
c
d
1
,
we find that the continuity equation (5) is identically satisfied and (8)-(9), after dropping the primes, reduce to
(11)
R
δ
(
∂
∂
t
+
∂
ψ
∂
y
∂
∂
x
-
∂
ψ
∂
x
∂
∂
y
)
(
δ
2
∂
2
ψ
∂
x
2
+
∂
2
ψ
∂
y
2
)
=
(
∂
2
∂
y
2
-
δ
2
∂
2
∂
x
2
)
S
x
y
-
M
2
∂
2
ψ
∂
y
2
+
δ
∂
2
∂
x
∂
y
(
S
x
x
-
S
y
y
)
,
(12)
[
E
1
∂
3
∂
x
3
+
E
2
∂
3
∂
x
∂
t
2
+
E
3
∂
3
∂
x
∂
t
]
η
=
δ
∂
S
x
x
∂
x
-
M
2
∂
ψ
∂
y
+
∂
S
x
y
∂
y
-
R
δ
(
∂
∂
t
+
∂
ψ
∂
y
∂
∂
x
-
∂
ψ
∂
x
∂
ψ
∂
y
)
∂
∂
y
,
ψ
y
=
0
at
y
=
±
η
.
In the above equations η
=
±
(
1
+
ε
sin
2
π
(
x
-
t
)
), ε
=
a
/
d
1 is the dimensionless amplitude, δ
=
d
1
/
λ denotes the characteristic ratio of transversal and axial scales of peristaltic motion, R
=
c
d
1
/
ν is the Reynolds number, M
=
σ
/
μ
B
0
d
1 denotes the Hartman number, E
1
=
-
τ
d
1
3
/
λ
3
μ
c, E
2
=
m
c
d
1
3
/
λ
3
μ, E
3
=
d
d
1
3
/
λ
3
μ are the nondimensional elasticity parameters, and β
*
=
β
c
2
/
d
1
2 is the dimensionless Burger’s fluid parameter.


## 4. Solution of the Problem

Due to nonlinearity of (11) it seems impossible to get a closed form solution. Therefore, some plausible assumptions must be made to put these equations in tractable form. We would like to mention here that without such assumptions one cannot obtain solution of these equations even for Newtonian fluid. Following [1], we assume that the δ is small and proceed to write
(13)
ψ
=
∑
k
=
0
∞
δ
k
ψ
k
,
S
i
j
=
∑
k
=
0
∞
δ
k
S
i
j
k
,
i
,
j
=
1,2
.
Substituting expressions (13) into (11) and (12) and collecting the terms of equal powers of δ, we get the following systems.


### 4.1. Zero-Order System

(14)
∂
2
S
0
x
y
∂
y
2
-
M
2
∂
2
ψ
0
∂
y
2
=
0
,
S
0
y
y
=
0
,
S
0
x
x
-
2
λ
1
∂
2
ψ
0
∂
y
2
S
0
x
y
-
2
λ
2
(
∂
2
ψ
0
∂
y
2
)
2
,
S
0
x
y
-
λ
1
∂
2
ψ
0
∂
y
2
S
0
y
y
=
∂
2
ψ
0
∂
y
2
,
[
E
1
∂
3
∂
x
3
+
E
2
∂
3
∂
x
∂
t
2
+
E
3
∂
2
∂
x
∂
t
]
η
=
-
M
2
∂
ψ
0
∂
y
+
∂
S
0
x
y
∂
y
,
∂
ψ
0
∂
y
=
0
,
at
y
=
±
η
.


### 4.2. First-Order System

(15)
R
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
+
M
2
∂
2
ψ
1
∂
y
2
=
∂
2
S
1
x
y
∂
y
2
+
∂
2
∂
x
∂
y
(
S
0
x
x
-
S
0
y
y
)
,
S
1
x
x
+
λ
1
[
2
∂
2
ψ
0
∂
y
2
S
1
x
y
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
x
x
-
2
∂
2
ψ
0
∂
x
∂
y
S
0
x
x
-
2
∂
2
ψ
1
∂
y
2
S
0
x
y
-
2
∂
2
ψ
0
∂
y
2
S
1
x
y
]
+
β
[
(
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
)
2
{
4
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
0
∂
y
2
HHHHHHH
-
2
(
∂
ψ
0
∂
y
∂
3
ψ
0
∂
x
∂
y
2
-
∂
ψ
0
∂
x
∂
3
ψ
0
∂
y
3
)
}
S
0
x
y
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
S
1
y
y
×
{
-
4
∂
2
ψ
0
∂
y
2
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
}
S
0
x
y
-
2
(
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
)
2
S
0
y
y
]
=
2
∂
2
ψ
0
∂
x
∂
y
-
4
λ
2
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
,
S
1
y
y
+
λ
1
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
y
y
+
2
∂
2
ψ
0
∂
x
∂
y
S
0
y
y
]
=
-
2
∂
2
ψ
0
∂
x
∂
y
,
S
1
x
y
+
λ
1
[
∂
2
ψ
0
∂
y
2
S
1
x
y
-
∂
2
ψ
1
∂
y
2
S
0
y
y
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
x
y
-
∂
2
ψ
0
∂
y
2
S
1
x
y
-
∂
2
ψ
1
∂
y
2
S
0
y
y
]
=
∂
2
ψ
1
∂
y
2
+
λ
2
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
=
∂
2
ψ
1
∂
y
2
+
λ
2
+
2
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
0
∂
y
2
]
,
∂
S
0
x
x
∂
x
+
∂
S
1
x
y
∂
y
-
R
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
ψ
0
∂
y
-
M
2
∂
ψ
1
∂
y
=
0
,
ψ
1
y
=
0
,
at
y
=
±
η
.


### 4.3. Second-Order System

(16)
R
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
1
∂
y
2
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
]
=
∂
2
S
2
x
y
∂
y
2
-
∂
2
S
0
x
y
∂
x
2
-
M
2
∂
2
ψ
2
∂
y
2
+
∂
2
∂
x
∂
y
(
S
1
x
x
-
S
1
y
y
)
,
S
2
x
y
+
λ
1
[
∂
2
ψ
1
∂
y
2
S
1
y
y
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
S
0
x
y
+
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
1
x
y
+
∂
2
ψ
0
∂
x
2
S
0
x
x
-
∂
2
ψ
0
∂
y
2
S
2
y
y
-
∂
2
ψ
1
∂
y
2
S
1
y
y
-
∂
2
ψ
2
∂
y
2
S
0
y
y
]
+
β
{
∂
2
∂
t
2
+
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
×
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
2
-
4
∂
2
ψ
0
∂
y
2
∂
ψ
0
∂
x
2
}
S
0
x
y
+
{
-
2
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
-
2
∂
2
ψ
0
∂
x
∂
y
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
-
(
∂
ψ
0
∂
y
∂
3
ψ
0
∂
x
2
∂
y
-
∂
ψ
0
∂
x
∂
3
ψ
0
∂
x
∂
y
2
)
-
2
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
}
S
0
y
y
=
(
∂
2
ψ
2
∂
y
2
-
∂
2
ψ
0
∂
x
2
)
+
λ
2
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
1
∂
y
2
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
+
2
∂
2
ψ
1
∂
x
∂
y
∂
2
ψ
0
∂
y
2
+
2
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
1
∂
y
2
]
,
S
2
y
y
+
λ
1
[
∂
2
ψ
0
∂
x
2
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
i
y
y
-
2
∂
2
ψ
0
∂
x
∂
y
S
1
y
y
-
2
∂
2
ψ
0
∂
x
2
S
0
x
y
]
+
β
[
{
(
∂
2
ψ
0
∂
x
∂
y
)
2
∂
∂
t
2
+
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
2
-
2
∂
3
ψ
0
∂
x
∂
y
∂
t
+
4
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
∂
∂
t
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
x
∂
y
+
4
(
∂
2
ψ
0
∂
x
∂
y
)
2
+
∂
2
ψ
0
∂
x
∂
y
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
2
∂
2
ψ
0
∂
y
2
∂
2
ψ
0
∂
x
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
}
S
0
y
y
]
=
-
2
∂
2
ψ
1
∂
x
∂
y
-
2
λ
2
[
(
∂
2
ψ
0
∂
x
∂
y
)
2
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
x
∂
y
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
+
∂
2
ψ
0
∂
x
2
∂
2
ψ
0
∂
y
2
]
,
∂
S
1
x
x
∂
x
+
∂
S
2
x
y
∂
y
-
M
2
∂
ψ
2
∂
y
-
R
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
ψ
1
∂
y
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
ψ
0
∂
y
]
=
0
,
∂
ψ
2
∂
y
=
0
,
at
y
=
±
η
.
The solution of these systems yields
(17)
u
=
u
0
+
δ
u
1
+
δ
2
u
2
,
where
(18)
u
0
=
L
(
sech
⁡
M
η
cosh
⁡
M
y
-
1
)
,
u
1
=
A
2
y
+
A
4
sinh
⁡
M
y
+
L
2
y
cosh
⁡
M
y
+
L
3
y
2
sinh
⁡
M
y
and value of u
2 is given in the appendix.The value of constants involved in (18) is given in the appendix.


## 4.1. Zero-Order System

(14)
∂
2
S
0
x
y
∂
y
2
-
M
2
∂
2
ψ
0
∂
y
2
=
0
,
S
0
y
y
=
0
,
S
0
x
x
-
2
λ
1
∂
2
ψ
0
∂
y
2
S
0
x
y
-
2
λ
2
(
∂
2
ψ
0
∂
y
2
)
2
,
S
0
x
y
-
λ
1
∂
2
ψ
0
∂
y
2
S
0
y
y
=
∂
2
ψ
0
∂
y
2
,
[
E
1
∂
3
∂
x
3
+
E
2
∂
3
∂
x
∂
t
2
+
E
3
∂
2
∂
x
∂
t
]
η
=
-
M
2
∂
ψ
0
∂
y
+
∂
S
0
x
y
∂
y
,
∂
ψ
0
∂
y
=
0
,
at
y
=
±
η
.


## 4.2. First-Order System

(15)
R
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
+
M
2
∂
2
ψ
1
∂
y
2
=
∂
2
S
1
x
y
∂
y
2
+
∂
2
∂
x
∂
y
(
S
0
x
x
-
S
0
y
y
)
,
S
1
x
x
+
λ
1
[
2
∂
2
ψ
0
∂
y
2
S
1
x
y
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
x
x
-
2
∂
2
ψ
0
∂
x
∂
y
S
0
x
x
-
2
∂
2
ψ
1
∂
y
2
S
0
x
y
-
2
∂
2
ψ
0
∂
y
2
S
1
x
y
]
+
β
[
(
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
)
2
{
4
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
0
∂
y
2
HHHHHHH
-
2
(
∂
ψ
0
∂
y
∂
3
ψ
0
∂
x
∂
y
2
-
∂
ψ
0
∂
x
∂
3
ψ
0
∂
y
3
)
}
S
0
x
y
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
S
1
y
y
×
{
-
4
∂
2
ψ
0
∂
y
2
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
}
S
0
x
y
-
2
(
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
)
2
S
0
y
y
]
=
2
∂
2
ψ
0
∂
x
∂
y
-
4
λ
2
∂
2
ψ
0
∂
y
2
∂
2
ψ
1
∂
y
2
,
S
1
y
y
+
λ
1
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
y
y
+
2
∂
2
ψ
0
∂
x
∂
y
S
0
y
y
]
=
-
2
∂
2
ψ
0
∂
x
∂
y
,
S
1
x
y
+
λ
1
[
∂
2
ψ
0
∂
y
2
S
1
x
y
-
∂
2
ψ
1
∂
y
2
S
0
y
y
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
0
x
y
-
∂
2
ψ
0
∂
y
2
S
1
x
y
-
∂
2
ψ
1
∂
y
2
S
0
y
y
]
=
∂
2
ψ
1
∂
y
2
+
λ
2
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
=
∂
2
ψ
1
∂
y
2
+
λ
2
+
2
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
0
∂
y
2
]
,
∂
S
0
x
x
∂
x
+
∂
S
1
x
y
∂
y
-
R
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
ψ
0
∂
y
-
M
2
∂
ψ
1
∂
y
=
0
,
ψ
1
y
=
0
,
at
y
=
±
η
.


## 4.3. Second-Order System

(16)
R
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
1
∂
y
2
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
]
=
∂
2
S
2
x
y
∂
y
2
-
∂
2
S
0
x
y
∂
x
2
-
M
2
∂
2
ψ
2
∂
y
2
+
∂
2
∂
x
∂
y
(
S
1
x
x
-
S
1
y
y
)
,
S
2
x
y
+
λ
1
[
∂
2
ψ
1
∂
y
2
S
1
y
y
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
S
0
x
y
+
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
1
x
y
+
∂
2
ψ
0
∂
x
2
S
0
x
x
-
∂
2
ψ
0
∂
y
2
S
2
y
y
-
∂
2
ψ
1
∂
y
2
S
1
y
y
-
∂
2
ψ
2
∂
y
2
S
0
y
y
]
+
β
{
∂
2
∂
t
2
+
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
×
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
2
-
4
∂
2
ψ
0
∂
y
2
∂
ψ
0
∂
x
2
}
S
0
x
y
+
{
-
2
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
-
2
∂
2
ψ
0
∂
x
∂
y
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
-
(
∂
ψ
0
∂
y
∂
3
ψ
0
∂
x
2
∂
y
-
∂
ψ
0
∂
x
∂
3
ψ
0
∂
x
∂
y
2
)
-
2
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
}
S
0
y
y
=
(
∂
2
ψ
2
∂
y
2
-
∂
2
ψ
0
∂
x
2
)
+
λ
2
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
1
∂
y
2
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
y
2
+
2
∂
2
ψ
1
∂
x
∂
y
∂
2
ψ
0
∂
y
2
+
2
∂
2
ψ
0
∂
x
∂
y
∂
2
ψ
1
∂
y
2
]
,
S
2
y
y
+
λ
1
[
∂
2
ψ
0
∂
x
2
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
S
i
y
y
-
2
∂
2
ψ
0
∂
x
∂
y
S
1
y
y
-
2
∂
2
ψ
0
∂
x
2
S
0
x
y
]
+
β
[
{
(
∂
2
ψ
0
∂
x
∂
y
)
2
∂
∂
t
2
+
∂
∂
t
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
2
-
2
∂
3
ψ
0
∂
x
∂
y
∂
t
+
4
∂
2
ψ
0
∂
x
∂
y
∂
∂
t
∂
∂
t
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
x
∂
y
+
4
(
∂
2
ψ
0
∂
x
∂
y
)
2
+
∂
2
ψ
0
∂
x
∂
y
(
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
+
2
∂
2
ψ
0
∂
y
2
∂
2
ψ
0
∂
x
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
}
S
0
y
y
]
=
-
2
∂
2
ψ
1
∂
x
∂
y
-
2
λ
2
[
(
∂
2
ψ
0
∂
x
∂
y
)
2
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
2
ψ
0
∂
x
∂
y
-
2
(
∂
2
ψ
0
∂
x
∂
y
)
2
+
∂
2
ψ
0
∂
x
2
∂
2
ψ
0
∂
y
2
]
,
∂
S
1
x
x
∂
x
+
∂
S
2
x
y
∂
y
-
M
2
∂
ψ
2
∂
y
-
R
[
(
∂
∂
t
+
∂
ψ
0
∂
y
∂
∂
x
-
∂
ψ
0
∂
x
∂
∂
y
)
∂
ψ
1
∂
y
+
(
∂
ψ
1
∂
y
∂
∂
x
-
∂
ψ
1
∂
x
∂
∂
y
)
∂
ψ
0
∂
y
]
=
0
,
∂
ψ
2
∂
y
=
0
,
at
y
=
±
η
.
The solution of these systems yields
(17)
u
=
u
0
+
δ
u
1
+
δ
2
u
2
,
where
(18)
u
0
=
L
(
sech
⁡
M
η
cosh
⁡
M
y
-
1
)
,
u
1
=
A
2
y
+
A
4
sinh
⁡
M
y
+
L
2
y
cosh
⁡
M
y
+
L
3
y
2
sinh
⁡
M
y
and value of u
2 is given in the appendix.The value of constants involved in (18) is given in the appendix.


## 5. Results and Discussion

The purpose of this section is to provide the graphical illustration of the solution obtained in previous section. Emphasis has been given to examine the influence of emerging parameters such asδ
,
R
,
M, and the rheological parameter of the Burger’s fluid β on the axial velocity and the stream function.In fact the novelty of the present analysis is based on the effects of Burger’s fluid parameterβ on various flow characteristics. This is because of the fact that the results in the absence of β are available already in the literature [13] for comparison. Particularly the velocity u is plotted for different values of the involved parameter including the Burger’s fluid parameter β. Figure 1 interprets the effects of Hartman number M and β on u. It is noted from this figure that u is an increasing function of β, while it decreases with M. The decrease of u with M is a result of the fact that magnetic force acts in the transverse direction and hence resists the flow due to peristalsis.Variation onu when E
1
=
1, E
2
=
0.2, E
3
=
0.5, ε
=
0.2, δ
=
0.01, t
=
0.2, R
=
0.01, λ
1
=
0.1, λ
2
=
0.05, and x
=
0.3 plotted for (a) for different values of M and (b) different Burger’s fluid parameter β.
(a)
(b)The effects ofR and δ on u can be seen through Figure 2. This figure shows that u increases by increasing R and δ. The Reynolds number represents the relative magnitude of inertial to the viscous force and with the increase of Reynolds number, the viscous forces decrease resulting in increases in velocity. Similarly with an increase in δ the wavelength of wave decreases which causes the increase in the velocity. It is worth mentioning that these results are in accordance with those obtained in [11] for a Maxwell fluid. A maximum in the velocity at center of the channel and an increase of velocity with R and δ are also observed from Figure 2.Variation onu when E
1
=
1, E
2
=
0.2, E
3
=
0.5, ε
=
0.2, β
=
0.05, t
=
0.2, M
=
2.5, λ
1
=
0.1, λ
2
=
0.0, and x
=
0.3 plotted for (a) for different values of R and (b) different values of δ.
(a)
(b)Figure3 is prepared to see the effects of complaint wall parameters E
1
,
E
2, and E
3 on longitudinal distribution of velocity u. This figure reveals that the velocity increases when both wall parameters E
1 and E
2 increase because of the increase in flexibility and elasticity of the wall. However, velocity decreases with E
3 due to the increase in rigidity of the wall. Note that we recover the results reported in [11] for λ
2
=
β
=
0.Variation onu when R
=
0.01, δ
=
0.01, ε
=
0.2, β
=
0.01, t
=
0.2, M
=
2.6, λ
1
=
0.1, λ
2
=
0.0, and x
=
0.3 plotted for (a) for different values of E
1, (b) for different values of E
2, and (c) for different values of E
3.
(a)
(b)
(c)Figure4 illustrates the effects of λ
1 and λ
2 on u. It is observed that when we increase the relaxation parameter the velocity u gets increased due to the increment in λ
1. The effects of retardation parameter λ
2 on the velocity are quite opposite to that of λ
1.Variation onu when E
1
=
1, E
2
=
0.2, E
3
=
0.5, ε
=
0.2, β
=
0.05, t
=
0.2, M
=
1, δ
=
0.01, R
=
0.01, and x
=
0.3 plotted for (a) for different values of λ
1 and (b) for different values of λ
2.
(a)
(b)The purpose of Figures5, 6, 7, and 8 is to capture the variation in size and circulation of the trapped bolus, in particular, the effects of E
1
,
E
2
,
β, and E
3. It is observed from these figures that the size of bolus decreases by increasing E
1 and E
2 while it increases with E
3 and β.Streamlines whenE
1
=
0.2, E
2
=
0.1, E
3
=
0.1, ε
=
0.1, δ
=
0.01, t
=
0.0, R
=
0.01, M
=
1, λ
1
=
0.1, and λ
2
=
0.05. (a) β
=
0.01 and (b) β
=
0.05.
(a)
(b)Streamlines whenβ
=
0.05, E
2
=
0.1, E
3
=
0.1, ε
=
0.1, δ
=
0.01, t
=
0.0, R
=
0.01, M
=
1, λ
1
=
0.1, and λ
2
=
0.05. (a) E
1
=
0.1 and (b) E
1
=
0.3.
(a)
(b)Streamlines whenβ
=
0.05, E
1
=
0.1, E
3
=
0.1, ε
=
0.1, δ
=
0.01, t
=
0.0, R
=
0.01, M
=
1, λ
1
=
0.1, and λ
2
=
0.05. (a) E
2
=
0.1 and (b) E
2
=
0.3.
(a)
(b)Streamlines whenβ
=
0.05, E
1
=
0.1, E
2
=
0.1, ε
=
0.1, δ
=
0.01, t
=
0.0, R
=
0.01, M
=
1, λ
1
=
0.1, and λ
2
=
0.05. (a) E
3
=
0.1 and (b) E
3
=
0.5.
(a)
(b)


## 6. Conclusions

Peristaltic motion of a (MHD) non-Newtonian fluid with Burger’s constitutive equation in channel with complaint walls is taken into account. The results for velocity and stream function are obtained and compared with the literature. It is observed that Reynolds number andδ cause an increase in the velocity while a decrease in velocity is noted by increasing Hartman number. In previous mathematical model [13], Burger’s parameter is absent and therefore such models cannot accurately describe the peristaltic motion of food stuff like cheese and so forth. The physical impact of this parameter is quite interesting. Our mathematical model predicts an increase in velocity by increasing Burger’s fluid parameter. The model further reveals that the size and circulation of the trapped bolus increase with Burger’s fluid parameter.



---

*Source: 236483-2014-10-28.xml*