# Effects of Heat Transfer and an Endoscope on Peristaltic Flow of a Fractional Maxwell Fluid in a Vertical Tube

Authors: H. Rachid
Journal: Abstract and Applied Analysis
 (2015)

---

## Abstract

We investigate the unsteady peristaltic transport of a viscoelastic fluid with fractional Maxwell model
through two coaxial vertical tubes. This analysis has been carried under low Reynolds number and long wavelength approximations.
Analytical solution of the problem is obtained by using a fractional calculus approach. The effects of Grashof number, heat
parameter, relaxation time, fractional time derivative parameter, amplitude ratio, and the radius ratio on the pressure gradient,
pressure rise, and the friction forces on the inner and outer tubes are graphically presented and discussed.


---

## Body

## 1. Introduction

Peristaltic transport is a mechanism of pumping fluids in tubes when progressive wave of area contraction or expansion propagates along the length on the boundary of a distensible tube containing fluid. In the physiological process, peristalsis or the blood flow is the major applications of this transport phenomenon. The contractions and relaxations of the walls propel the fluid. The first investigation to study the peristaltic pumping was by Latham [1]. The peristaltic transport of a Newtonian fluid at long wavelengths at low Reynolds number has been studied with particles suspended by Kaimal [2] or without particles by Shapiro et al. [3] and by Mishra and Rao [4]. The non-Newtonian effects of fluids on the peristaltic transport have been studied for a third-grade fluid by Ali et al. [5], for a Jeffrey fluid by Abd-Alla et al. [6], for a Maxwell model by Rachid and Ouazzani [7], for a Herschel-Bulkley fluid by Medhavi [8], or for a power-law fluid by Shukla and Gupta [9]. The effect of an endoscope on peristaltic flow for a Newtonian fluid through a nonuniform annulus has been analyzed by Mekheimer [10]. This effect has been also studied for a Johnson Segalman fluid by Akbar and Nadeem [11] or for a Maxwell fluid by Husseny et al. [12]. Recently, the peristaltic transport of non-Newtonian fluids has gained considerable interest because of its applications in industry and biology. Mathematically, the non-Newtonian fluids have a nonlinear relationship between the stress and the rate of strain. It is difficult to propose a single model which can exhibit all the properties of non-Newtonian fluids. To describe the viscoelastic properties of such fluids recently, constitutive equations with ordinary and fractional time derivatives have been introduced. Fractional calculus has proved to be very successful in the description of constitutive relations of viscoelastic fluids. The starting point of the fractional derivative model of viscoelastic fluids is usually a classical differential equation which is modified by replacing the time derivative of an integral order with fractional order and may be formulated both in the Riemann-Liouville and in the Caputo sense [13]. This generalization allows one to define precisely noninteger order integrals or derivatives and this plays an important role in studying the valuable tool of viscoelastic properties. In recent years and considering the relevance of fractional models of viscoelastic fluids, a number of articles have addressed unsteady flows of viscoelastic fluids with the fractional Maxwell model [14, 15], fractional second-grade fluid [16, 17], fractional Oldroyd-B model [18, 19], and fractional Burgers model [20]. For example, Tripathi [21] studied a mathematical model based on viscoelastic fractional Oldroyd-B model for the peristaltic flow of chyme in small intestine, which is assumed to be in the form of an inclined cylindrical tube. The peristaltic transport for a fractional second-grade model through a cylindrical tube has been investigated by Tripathi [22]. Rathod and Channakote [23] analyzed the interaction of heat transfer on peristaltic pumping of a fractional second-grade fluid through a vertical cylindrical tube by using Caputo’s definition. For the fractional Maxwell model, Tripathi et al. [24] studied the peristaltic flow through a channel where the homotopy perturbation method (HPM) and Adomian decomposition method (ADM) are used. For the fractional Burger’s model Tripathi et al. [25] studied the influence of slip condition on peristaltic transport by using a homotopy analysis method (HAM). This generalized model of fluid has been investigated for an infinite oscillating plate by Khan et al. [26], for an orthogonal rheometer by Siddiqui et al. [27] or between concentric cylinders by Shah [28]. Heat transfer analysis is one of the principle axes of research in chemical engineering, but nowadays due to its wide applications in biofluid mechanics it has great interaction with peristaltic motion. The effect of heat transfer on peristaltic transport of a Newtonian fluid through a porous medium in a vertical tube under the effect of a magnetic field has been studied by Vasudev et al. [29]. This effect on peristaltic flow with slip condition and variable viscosity in an asymmetric channel has been analyzed by Hayat and Abbasi [30]. Hayat et al. [31] investigated the effects of heat and mass transfer and wall properties on the peristalsis in a power-law fluid.The aim of this paper is to analytically study the effect of heat transfer on peristaltic flow of a viscoelastic fluid with fractional Maxwell model in the gap between two coaxial vertical tubes. The influence of different physical parameters on the pressure gradient, pressure rise, and the fictional forces has been graphically shown and discussed.


## 2. Formulation and Analysis

Consider the peristaltic flow of a viscoelastic fluid with fractional Maxwell model through the gap between two coaxial vertical tubes. The inner tube is rigid (endoscope) maintained at a temperatureT
1, and the outer tube has a sinusoidal wave traveling down its wall and it is exposed to temperature T
0 as seen in Figure 1.Figure 1
Geometry of the problem.In a cylindrical coordinate system(
R
¯
,
Z
¯
), the equations for the tube walls in their dimensional form are given by(1)
r
¯
1
=
a
1
,
r
¯
2
=
a
2
+
b
sin
⁡
2
π
λ
Z
¯
-
c
t
¯
,where a
1 is the radius of the endoscope, a
2 is the average radius, b is the amplitude of the wave, λ is the wavelength, and c is the wave speed of the outer tube.The equation of the fractional Maxwell fluid is given by(2)
1
+
λ
1
¯
α
1
D
~
t
¯
α
1
S
¯
=
μ
γ
˙with t
¯, S
¯, γ
˙, and μ being the time, shear stress, rate of shear strain, and the viscosity, respectively. λ
1
¯ is the relaxation time and α
1 is the fractional time derivative parameter such that 0
≤
α
1
≤
1. D
~
t
α
1 is the upper convected fractional derivative defined by(3)
D
~
t
¯
α
1
S
¯
=
D
t
¯
α
1
S
¯
+
V
¯
·
∇
S
¯
-
L
¯
S
¯
-
S
¯
L
¯
Tin which(4)
γ
˙
=
∇
V
¯
+
∇
V
¯
Tand D
t
¯
α
1
=
∂
t
¯
α
1 is the fractional differentiation operator of order α
1 with respect to t and may be defined as [32](5)
D
t
α
1
f
t
=
1
Γ
1
-
α
1
d
d
t
∫
0
t
f
ξ
t
-
ξ
α
1
d
ξ
,
0
≤
α
1
≤
1
.Here Γ
(
·
) denotes the Gamma function.We choose a cylindrical coordinate system(
r
¯
,
z
¯
) where z
¯-axis is the longitudinal direction and the r
¯-axis is transverse to it.The equations of motion of the flow for an incompressible fluid are given by(6)
ρ
∂
u
¯
∂
t
¯
+
u
¯
∂
u
¯
∂
r
¯
+
w
¯
∂
u
¯
∂
z
¯
=
-
∂
p
¯
∂
r
¯
+
1
r
¯
∂
∂
r
¯
r
¯
S
¯
r
¯
r
¯
+
∂
S
¯
r
¯
z
¯
∂
z
¯
,
ρ
∂
w
¯
∂
t
¯
+
u
¯
∂
w
¯
∂
r
¯
+
w
¯
∂
w
¯
∂
z
¯
=
-
∂
p
¯
∂
z
¯
+
1
r
¯
∂
∂
r
¯
r
¯
S
¯
r
¯
z
¯
+
∂
S
¯
z
¯
z
¯
∂
z
¯
+
ρ
g
α
T
¯
-
T
0
,
ρ
c
p
u
¯
∂
T
¯
∂
r
¯
+
w
¯
∂
T
¯
∂
z
¯
=
k
∂
2
T
¯
∂
r
¯
2
+
1
r
¯
∂
T
¯
∂
r
¯
+
∂
2
T
¯
∂
z
¯
2
+
Q
0
,
∂
u
¯
∂
r
¯
+
u
¯
r
¯
+
∂
w
¯
∂
z
¯
=
0
,where ρ is the fluid density, T
¯ is the temperature, α is the coefficient of linear thermal expansion of the fluid, k denotes thermal conductivity, c
p denotes specific heat at constant pressure, Q
0 is the constant heat addition/absorption, u
¯ and w
¯ are the velocity components in the wave frame, and p
¯ is the pressure.We assume that the extra stressS
¯ depends on r and t only. After using the initial condition S
¯
(
t
¯
=
0
), (2)–(4) yield S
¯
r
¯
r
¯
=
S
¯
z
¯
z
¯
=
S
¯
r
¯
θ
¯
=
S
¯
θ
¯
z
¯
=
0 and(7)
1
+
λ
1
¯
α
1
∂
α
1
∂
t
¯
α
1
S
¯
r
¯
z
¯
=
μ
∂
w
¯
∂
r
¯
,where S
¯
r
¯
z
¯ is the tangential stress.For carrying out further analysis, we introduce the following dimensionless parameters:(8)
z
=
z
¯
λ
;
Z
=
Z
¯
λ
;
r
1
=
r
¯
1
a
2
;
r
2
=
r
¯
2
a
2
;
R
=
R
¯
a
2
;
t
=
c
t
¯
λ
;
u
=
λ
u
¯
a
2
c
;
U
=
λ
U
¯
a
2
c
;
w
=
w
¯
c
;
W
=
W
¯
c
;
λ
1
=
c
λ
1
¯
λ
;
ɛ
=
a
1
a
2
,
p
=
a
2
2
p
¯
μ
λ
c
;
θ
=
T
¯
-
T
0
T
1
-
T
0
;
Q
=
Q
¯
π
c
a
2
2
;
δ
=
a
2
λ
;
Pr
⁡
=
μ
c
p
k
;
R
e
=
ρ
c
a
2
μ
0
;
β
=
a
2
2
Q
0
k
T
1
-
T
0
;
Gr
=
ρ
g
α
a
2
3
T
1
-
T
0
μ
,where ɛ is the radius ratio, δ is the wave number, R
e is the Reynolds number, Pr
⁡ is the Prandtl number, Gr is the Grashof number, β is the nondimensional heat source/sink parameter, and ϕ is the amplitude ratio with 0
<
ϕ
<
1
-
ɛ.Using the above nondimensional quantities and under the assumptions of long wavelength approximation (i.e.,δ
≪
1 or λ
≫
a
2) and the low Reynolds number (i.e., R
e
→
0), (6) yield(9)
∂
p
∂
r
=
0
,
(10)
1
+
λ
1
α
1
D
t
α
1
∂
p
∂
z
-
Gr
θ
=
∂
2
w
∂
r
2
+
1
r
∂
w
∂
r
,
(11)
∂
u
∂
r
+
u
r
+
∂
w
∂
z
=
0
,
(12)
∂
2
θ
∂
r
2
+
1
r
∂
θ
∂
r
+
β
=
0
.The boundary conditions are(13)
w
=
-
1
;
θ
=
1
at
r
=
r
1
=
ɛ
,
w
=
-
1
;
θ
=
0
at
r
=
r
2
,where r
2
=
1
+
ϕ
sin
⁡
2
π
z is the dimensionless equation of the outer tube radius in the wave frame.Integrating (12) with the boundary conditions (13), we find the dimensionless temperature as follows:(14)
θ
=
ln
⁡
r
/
r
2
ln
⁡
ɛ
/
r
2
+
β
4
ε
2
ln
⁡
r
/
r
2
-
r
2
2
ln
⁡
r
/
ɛ
ln
⁡
ɛ
/
r
2
-
β
4
r
2
.Use the transformations between the laboratory and the wave frames, in the dimensionless form, given by(15)
z
=
Z
-
t
;
r
=
R
;
u
=
U
;
w
=
W
-
1
,where U
¯ and W
¯ are the velocity components in the fixed frame.Substituting (14) in (10) and using the same boundary conditions, we obtain the velocity profile of the fluid as(16)
W
=
1
4
f
r
2
-
a
21
ln
⁡
r
+
a
25
∂
p
∂
z
+
1
64
f
Gr
β
r
4
+
4
β
a
24
ln
⁡
r
ɛ
-
a
41
ln
⁡
r
r
2
+
a
42
r
2
+
1
64
f
Gr
a
43
,where f
(
t
)
=
1
+
λ
1
α
1
D
t
α
1
=
f.Using the definition of the fractional differential operator (5) we find the expression of f as follows:(17)
f
=
1
+
λ
1
α
1
t
-
α
1
Γ
1
-
α
1
.The volume rate of flow in the fixed coordinate system (R, Z) is given as(18)
Q
Z
,
t
=
2
∫
ɛ
r
2
W
R
d
R
=
f
a
61
∂
p
∂
z
+
a
62
.Using the transformation (15) we find(19)
Q
Z
,
t
=
q
+
r
2
2
-
ε
2
,where q is the volume flow rate in the moving coordinate system given by(20)
q
=
2
∫
ɛ
r
2
w
r
d
r
.The time-averaged flow rate is defined by(21)
Q
=
∫
0
1
Q
Z
,
t
d
t
=
q
+
1
+
ϕ
2
2
-
ε
2
.Eliminating the flow rateq between (19) and (21) we find(22)
Q
Z
,
t
=
Q
+
r
2
2
-
1
-
ϕ
2
2
.From (18) and (22) we obtain the pressure gradient as follows:(23)
∂
p
∂
z
=
1
a
61
Q
+
r
2
2
-
1
-
ϕ
2
/
2
f
-
a
62
.


## 3. The Pumping Characteristics

The pressure riseΔ
p and the frictional force F
λ at the walls of the inner and outer tubes, in the nondimensional form, are given by(24)
Δ
p
=
∫
0
1
∂
p
∂
z
d
z
,
F
i
=
∫
0
1
ε
2
-
∂
p
d
z
d
z
,
F
o
=
∫
0
1
r
2
2
-
∂
p
d
z
d
z
.


## 4. Results and Discussions

We note that the fractional Maxwell model reduces to a Newtonian fluid whenλ
1
¯
=
0.The analytical expressions for the pressure gradient∂
p
/
∂
z, the pressure rise Δ
p, and the frictional forces on the inner and outer tubes, respectively, F
(
i
) and F
(
o
), are derived in the previous sections. In order to compute these physical quantities with respect to interest parameters of problem we observe that the integrals in (24) are not integrable in the closed form. They are evaluated numerically using a mathematics software.Pressure gradient∂
p
/
∂
z against the axial distance z on one wavelength for different given values of the Grashof number Gr, the heat parameter β, the relaxation time λ
1, the fractional time derivative parameter α
1, the amplitude ratio ϕ, and the radius ratio ɛ is plotted in Figures 2–4. These figures show that, in the wider part of the tube (
z
∈
[
0,0.5
]
), ∂
p
/
∂
z is relatively small, where the flow can easily pass without giving any large pressure gradient. However, in the narrow part (
z
∈
[
0.5,1
]
) a much larger pressure gradient is required to maintain the same flux to pass it, especially for the narrowest position near z
=
0.75. Moreover the maximum of the pressure gradient increases with increasing Gr, ϕ, ɛ, and α
1 while it decreases with the increase in λ
1 and β.Figure 2
Pressure gradient∂
p
/
∂
z versus axial distance z for different values of the relaxation time λ
1 and radius ratio ɛ with Q
=
-
2, Gr
=
3, α
1
=
0.4, β
=
10, ϕ
=
0.4, and t
=
1.Figure 3
Pressure gradient∂
p
/
∂
z versus axial distance z for different values of the parameter α
1 and amplitude ratio ϕ with Q
=
-
2, Gr
=
3, ɛ
=
0.35, β
=
10, λ
1
=
1, and t
=
1.Figure 4
Pressure gradient∂
p
/
∂
z versus axial distance z for different values of the heat parameter β and the Grashof number Gr with Q
=
-
2, λ
1
=
1, α
1
=
0.4, ɛ
=
0.35, ϕ
=
0.4, and t
=
1.The effects of various parameters on the pressure riseΔ
p versus the time-averaged flow rate Q are shown in Figures 5–7. It is known that we have the following peristaltic regions, pumping region (Δ
p
>
0), free-pumping region (Δ
p
=
0), and copumping region (Δ
p
<
0). These figures show a linear relationship between Δ
p and Q in all the regions; that is, the pressure rise decreases with increasing the time-averaged flow rate. From Figure 5 we observe that the time-averaged flow rate Q increases with increasing Gr while it decreases with the increase in β in all the regions. In addition, Δ
p and the flow rate Q and the pumping increase with increasing Gr while they decrease with the increase in β. Figures 6-7 show that for different values of ɛ, λ
1, ϕ, and α
1 the curves are intersecting in the pumping region. In addition, for a given value of the pressure rise Δ
p we observe that in both the free-pumping and copumping regions the time-averaged flow rate Q increases with increasing λ
1 while it decreases with the increase in α
1. The same figures show that Q decreases in the free-pumping and copumping regions with increasing ɛ and ϕ.Figure 5
Pressure riseΔ
p versus the time-averaged flow rate Q for different values of the heat parameter β and the Grashof number Gr with λ
1
=
1, α
1
=
0.4, ɛ
=
0.35, ϕ
=
0.4, and t
=
1.Figure 6
Pressure riseΔ
p versus the time-averaged flow rate Q for different values of the relaxation time λ
1 and radius ratio ɛ with Gr
=
3, α
1
=
0.4, β
=
10, ϕ
=
0.4, and t
=
1.Figure 7
Pressure riseΔ
p versus the time-averaged flow rate Q for different values of the parameter α
1 and amplitude ratio ϕ with Gr
=
3, ɛ
=
0.35, β
=
10, λ
1
=
1, and t
=
1.In Figures8–13 we plotted the frictional forces F
(
i
) and F
(
o
) on the inner and outer tubes, respectively. These figures show that F
(
i
) and F
(
o
) have an opposite behavior compared to the pressure rise Δ
p versus the physical parameters. In addition the frictional force on the outer tube F
(
o
) is greater than the frictional force on the outer tube F
(
i
).Figure 8
Frictional force on the inner tubeF
(
i
) versus the time-averaged flow rate Q for different values of heat parameter β and the Grashof number Gr with λ
1
=
1, α
1
=
0.4, ɛ
=
0.35, ϕ
=
0.4, and t
=
1.Figure 9
Frictional force on the inner tubeF
(
i
) versus the time-averaged flow rate Q for different values of the relaxation time λ
1 and radius ratio ɛ with Gr
=
3, α
1
=
0.4, β
=
10, ϕ
=
0.4, and t
=
1.Figure 10
Frictional force on the inner tubeF
(
i
) versus the time-averaged flow rate Q for different values of the parameter α
1 and amplitude ratio ϕ with Gr
=
3, ɛ
=
0.35, β
=
10, λ
1
=
1, and t
=
1.Figure 11
Frictional force on the outer tubeF
(
o
) versus the time-averaged flow rate Q for different values of heat parameter β and the Grashof number Gr with λ
1
=
1, α
1
=
0.4, ɛ
=
0.35, ϕ
=
0.4, and t
=
1.Figure 12
Frictional force on the outer tubeF
(
o
) versus the time-averaged flow rate Q for different values of the relaxation time λ
1 and radius ratio ɛ with Gr
=
3, α
1
=
0.4, β
=
10, ϕ
=
0.4, and t
=
1.Figure 13
Frictional force on the outer tubeF
(
o
) versus the time-averaged flow rate Q for different values of the parameter α
1 and amplitude ratio ϕ with Gr
=
3, ɛ
=
0.35, β
=
10, λ
1
=
1, and t
=
1.


## 5. Conclusions

We have analytically studied the effect of heat transfer on peristaltic flow of fractional Maxwell model in the gap between two vertical coaxial tubes. The problem is simplified under the assumptions of long wavelength approximation and low Reynolds number. Analytical solution of problem is obtained by using a fractional calculus approach. The pressure gradient, pressure rise, and frictional forces on the endoscope and on the outer tube are discussed with the physical parameters, Grashof numberGr, the heat parameter β, the relaxation time λ
1, the fractional time derivative parameter α
1, the amplitude ratio ϕ, and the radius ratio ɛ. The graphical solutions have shown the following:(1)
The maximum of the pressure gradient∂
p
/
∂
z increases with increasing Gr, ϕ, ɛ, and α
1 while it decreases with the increase in λ
1 and β.
(2)
The time-averaged flow rateQ increases with increasing Gr while it decreases with the increase in β in all the regions of pumping for a given value of the pressure rise Δ
p.
(3)
Q increases with the increase in λ
1 in both the free-pumping and co-pumping regions while we observe an opposite behavior versus α
1 in these two regions.
(4)
Q decreases with increasing ϕ and ɛ in the free-pumping and copumping regions.
(5)
The frictional forcesF
(
i
) and F
(
o
) have opposite behavior compared to the pressure rise Δ
p versus the physical parameters.



---

*Source: 360918-2015-10-22.xml*