particulate suspension slip flow induced by peristaltic ... · original article particulate...
TRANSCRIPT
Alexandria Engineering Journal (2016) xxx, xxx–xxx
HO ST E D BY
Alexandria University
Alexandria Engineering Journal
www.elsevier.com/locate/aejwww.sciencedirect.com
ORIGINAL ARTICLE
Particulate suspension slip flow induced by
peristaltic waves in a rectangular duct: Effect of
lateral walls
* Corresponding author.E-mail addresses: [email protected], safia_akram@yahoo.
com (S. Akram).
Peer review under responsibility of Faculty of Engineering, Alexandria
University.
http://dx.doi.org/10.1016/j.aej.2016.09.0111110-0168 � 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induced by peristaltic waves in a rectangular duct: Effect of lateral walls, AlEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
Safia Akrama,*, Kh.S. Mekheimer
b,c, Y. Abd Elmaboud
d,e
aDepartment of Basic Sciences, MCS, National University of Sciences and Technology, Islamabad 4400, PakistanbMathematics Department, Faculty of Science, Taif University Hawia, P.O. Box 888, Taif, Saudi ArabiacMathematics Department, Faculty of Science, Al-Azhar University, Nasr City, 11884 Cairo, EgyptdMathematics Department, Faculty of Science and Arts, Khulais, King Abdulaziz University (KAU), Saudi ArabiaeMathematics Department, Faculty of Science, Al-Azhar University(Assiut Branch), Assiut, Egypt
Received 29 March 2015; revised 19 September 2016; accepted 20 September 2016
KEYWORDS
Fluid suspension;
Lateral walls;
Peristaltic pumping;
Partial slip
Abstract This paper looks at the influence of lateral walls on peristaltic transport of a particle fluid
suspension model applied in a non-uniform rectangular duct with slip boundaries. The peristaltic
waves propagate on the horizontal sidewalls of a rectangular duct. The flow analysis has been devel-
oped for low Reynolds number and long wavelength approximation. Exact solutions have been
established for the axial velocity and stream function. The effects of aspect the ratio b (ratio of
height to width) and the volume fraction density of the particles C on the pumping characteristics
are discussed in detail. The expressions for the pressure rise and friction forces on the wall of a rect-
angular duct were computed numerically and were plotted with variation of the flow rate for differ-
ent values of the parameters. It is observed that in the peristaltic pumping ðDp > 0;Q > 0Þ and
retrograde pumping ðDp > 0;Q < 0Þ regions the pumping rate increases with an increase in M,
while in the copumping region ðDp < 0;Q > 0Þ the behavior is quite opposite. Furthermore it is also
observed that the pressure rise increases in the upper half of the channel and decreases in the lower
half of the channel with the increase in lslip parameter.� 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an
open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
1. Introduction
In the recent years peristaltic transport becomes an importantsubject for the researcher in biofluid. The peristaltic transportis very important because, it is a primary transport mechanism
inherent in many tubular organs of the human body such as
exandria
2 S. Akram et al.
the gastrointestinal tract, the urethra, the ureter and arterioles.The mechanism of peristaltic transport has been employed forindustrial applications such as sanitary fluid transport, blood
pumps in heart lung machine and transport of corrosive fluids.Peristaltic flow is generated in a channel (or a circular tube)when a progressive wave travels along the wall. Many studies
have been carried out for understanding the characteristics ofthe transport mechanism associated with peristaltic flow underthe assumption of low Reynolds number and infinitely long
wavelength such as Jaffrin and Shapiro [1] explained the basicprinciples of peristaltic pumping and brought out clearly thesignificance of the various parameters governing the flow. Anumber of analytical, numerical and experimental recent stud-
ies of peristaltic flows of different fluids have been reported[2–13].
The theory of the two phase fluid is very useful in under-
standing a number of diverse physical problems includingpowder technology, fluidization. The continuum theory ofmixtures is applicable to hydrodynamics of biological systems,
because it provides an improved understanding of diverse sub-jects such as diffusion of proteins, the rheology of blood,swimming of micro-organisms and particle deposition. The
model of a particulate suspension with peristaltic transport isinvestigated by many authors [14–18]. There are a recent fewstudies take into account the effect of the sidewalls in peri-staltic transport such as Reddy et al. [19], and Nadeem and
Akram [20]. None of the previous studies take into considera-tion slip boundaries and the sidewall effects especially in thenon-uniform channel when the fluid model is a two phase fluid.
So our motivation was to study the influence of lateral walls onperistaltic transport of a particle fluid suspension modelapplied in a non-uniform rectangular duct with slip
boundaries.
2. Mathematical formulation
Let us consider the two-phase (fluid particles) flow in a non-uniform duct of rectangular cross section. The duct walls areflexible, and an infinite train of sinusoidal waves propagate
with constant velocity c along the walls parallel to the XYplane in the axial direction. Cartesian coordinate system (X,Y,Z) is with X, Y, and Z axes corresponding to axial, lateral,and vertical directions, respectively, of a rectangular duct. We
assume that there is no flow in the lateral direction. So, thevelocity vector in this direction will be zero. The sinusoidalwaves propagating along the channel walls are of the following
forms:
Z0 ¼ H0ðX0; t0Þ ¼ �a� kx� b sin2pkðX0 � ct0Þ
� �; ð1Þ
where a is the half width of the channel, b is amplitude of thewave, k is the wavelength, d is the channel width, k(�1) is aconstant whose magnitude depends on the length of the chan-nel, c is the velocity of propagation, and t0 is the time. The
equations governing the conservation of mass and linearmomentum for both the fluid and particle phases using a con-tinuum approach are expressed as follows:
Fluid phase
@
@X0 ðð1� CÞU0fÞ þ
@
@Z0 ðð1� CÞV0fÞ ¼ 0; ð2Þ
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
ð1� CÞqf
@U0f
@t0þU0
f
@U0f
@X0 þ V0f
@U0f
@Z0
� �
¼ �ð1� CÞ @P0
@X0 þ ð1� CÞlsðCÞr2U0f þ CS U0
P �U0f
� �; ð3Þ
ð1� CÞqf
@V0f
@t0þU0
f
@V0f
@X0 þ V0f
@V0f
@Z0
� �
¼ �ð1� CÞ @P0
@Z0 þ ð1� CÞlsðCÞr2V0f þ CS V0
P � V0f
� �; ð4Þ
Particulate phase
@
@X0 CU0P
� þ @
@Z0 CV0P
� ¼ 0; ð5Þ
CqP
@U0P
@t0þU0
P
@U0P
@X0 þ V0P
@U0P
@Z0
� �
¼ �C@P0
@X0 þ CS U0f �U0
P
� �; ð6Þ
CqP
@V0P
@t0þU0
P
@V0P
@X0 þ V0P
@V0P
@Z0
� �
¼ �C@P0
@Z0 þ CS V0f � V0
P
� �; ð7Þ
where U0f; V0
f
� �and U0
P; V0P
� are the velocity components of
the fluid and the particle phase respectively, P0 is the pressure,C is the volume fraction density of the particles [17], ls(C) isthe mixture viscosity (effective or apparent viscosity of suspen-sion) and S is the drag coefficient of interaction for the force
exerted by one phase on the other. The expression for the dragcoefficient of interaction, S and the empirical relation for thevelocity of the suspension, ls for the present problem is
selected as done by Srivastava and Saxena [17]
S ¼ 9
2
l0
a20k0ðCÞ;
k0ðCÞ ¼ 4þ ½8C� 3C2�12 þ 3C
ð2� 3CÞ2 ; ð8Þ
ls ¼ lsðCÞ ¼l0
1�mC;
m ¼ 0:070 exp 2:49Cþ 1107
Texpð�1:69CÞ
� �; ð9Þ
where a0 is the radius of each solid particle suspended in thefluid, l0 is the constant fluid viscosity and T is measured in
absolute temperature (K0). The formula (9) has been testedby Charm and Kurland [21] by using a cone and a plateviscometer, and it has been proclaimed that it is reasonably
accurate up to C= 0.6.Let us define a wave frame ðx0; y0; z0Þ moving with the
velocity c away from the fixed frame ðX0;Y0;Z0Þ by thetransformation
x0 ¼ X0 � ct0; y0 ¼ Y0; z0 ¼ Z0; u0f;p ¼ U0f;p � c;
v0f;p ¼ V0f;p; p0ðx; zÞ ¼ P0ðX;Z; tÞ: ð10Þ
where u0f;p and v0f;p are the velocity components of fluid and par-
ticle respectively, p0 and P0 are the pressures in wave and fixed
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria
-0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
Q
p
M = 0.2M = 0.3
M = 0.4
M = 0.5
Δ
Figure 1 Variation of Dp with Q for different values of M at
/ ¼ 0:6, lslip ¼ 0:3, C ¼ 0:5, K ¼ 0:05, l ¼ 0:7 and b ¼ 1:2.
Particulate suspension slip flow 3
frame respectively. Consider the following non-dimensional
variables:
x ¼ x0
k; y ¼ y0
d; z ¼ z0
a; d ¼ a
k; up ¼
u0pc; uf ¼
u0fc; t ¼ ct0
k;
h ¼ H0
a; p ¼ a2p0
lck;
b ¼ a
d; / ¼ b
a; vf ¼
v0fcd
; vp ¼v0pcd
; Re ¼ acqf
lsð1� CÞ ;
M ¼ sa2
lsð1� CÞ ;
N ¼ sa2qf
qplsð1� CÞ ; �l ¼ ls
l; �q ¼ qp
qf
: ð11Þ
where b is an aspect ratio (aspect ratio b < 1 means thatheight is less compared to width, and b = 0 corresponds to atwo-dimensional channel. When b = 1, the rectangular duct
becomes a square duct and for b > 1, the height is more com-pared to width.). Re is a suspension Reynolds number, (M,N)are the suspension parameters, d is the wave number, and / is
the amplitude ratio. In moving frame, the nondimensionalequations of motion for both the fluid and particle phaseare
@
@xðð1� CÞufÞ þ d
@
@zðð1� CÞvfÞ ¼ 0; ð12Þ
Red�lð1� CÞ uf@uf@x
þ vf@uf@z
� �
¼ � @p
@xþ �l d2
@2uf@x2
þ b2 @2uf@y2
þ @2uf@z2
� �þMC�lðup � ufÞ;
ð13Þ
Red3�lð1� CÞ uf@vf@x
þ vf@vf@z
� �
¼ � @p
@zþ �ld2 d2
@2vf@x2
þ b2 @2vf@y2
þ @2vf@z2
� �þ Cd2�lMðvp � vfÞ;
ð14Þ
@
@xðCupÞ þ d
@
@zðCvpÞ ¼ 0; ð15Þ
�qRedð1� CÞ�l up@up@x
þ vp@up@z
� �
¼ � @p
@xþMð1� CÞ�lðuf � upÞ; ð16Þ
�qRed3ð1� CÞ�l up@vp@x
þ vp@vp@z
� �
¼ � @p
@zþMð1� CÞ�ld2ðvf � vpÞ; ð17Þ
under lubrication approach Eqs. (12)-(17) reduced to
@p
@x¼ �l b2 @
2uf@y2
þ @2uf@z2
� �þMC�lðup � ufÞ; ð18Þ
@p
@z¼ 0; ð19Þ
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
@p
@x¼ Mð1� CÞ�lðuf � upÞ; ð20Þ
The corresponding no-slip boundary conditions are
uf ¼ �lslip@uf@y
� 1 at y ¼ 1;
uf ¼ lslip@uf@y
� 1 at y ¼ �1;
uf ¼ �1 at z ¼ �hðxÞ ¼ �1� Kx� / sin 2px; ð21Þwhere 0 6 / 6 1, / ¼ 0 for straight duct and / ¼ 1 corre-sponds to total occlusion. Combining Eqs. (18) and (20) we get
b2 @2uf@y2
þ @2uf@z2
¼ 1
ð1� CÞ�l@p
@x: ð22Þ
3. Solution of the problem
Using the similar procedure as discussed by Reddy et al. [19]the solution of Eq. (22) satisfying the boundary conditions(21) can be directly written as
uf ¼�1þ 1
2ð1�CÞ�l
� dp
dxz2�h2ðxÞþ4
h
X1n¼1
ð�1Þn coshðany=bÞcosðanzÞa3nðcoshðan=bÞþ lslip
anb sinhðan=bÞÞ
( ):
ð23Þthe particulate velocity will be in the form
up ¼�1þ 1
2ð1�CÞ�l
� dp
dx
�2
Mþ z2�h2ðxÞþ4
h
X1n¼1
ð�1Þn coshðany=bÞcosðanzÞa3nðcoshðan=bÞþ lslip
anb sinhðan=bÞÞ
( ):
ð24Þ
where an ¼ ð2n� 1Þp=2hðxÞ.
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria
4 S. Akram et al.
The volumetric flow rate in the rectangular duct in the waveframe (in vertical half) is given by
qf ¼ð1�CÞZ h
0
Z 1
0
ðufÞdydz;
¼ð1�CÞ �hþ 1
2ð1�CÞ�ldp
dx
� �2
3h3þ4
bh
X1n¼1
ð�1Þn sinhðan=bÞa5nðcoshðan=bÞþ lslip
anb sinhðan=bÞÞ
!):
qp ¼C
Z h
0
Z 1
0
ðupÞdydz;
¼ C
ð1�CÞ� � h
M�ldp
dxþ 1
�ldp
dx
� �h3
3þ2b
h
X1n¼1
ð�1Þn sinhðan=bÞa5nðcoshðan=bÞþ lslip
anb sinhðan=bÞÞ
!)�hð1�CÞ
):
-0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
-1
-0.5
0
0.5
Q
p
lslip
= 0.1
lslip
= 0.3
lslip
= 0.5
lslip
= 0.7
Figure 2 Variation of Dp with Q for different values of lslip at
/ ¼ 0:5, M ¼ 1:3, C ¼ 0:3, K ¼ 0:005, l ¼ 0:7 and b ¼ 1:2.
0 0.2 0.4 0.6 0.8 1-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
Q
K = 0.0
K = 0.05
K = 0.1
K = 0.2
pΔ
Figure 3 Variation of Dp with Q for different values of K at
/ ¼ 0:6, M ¼ 0:5, C ¼ 0:5, lslip ¼ 0:3, l ¼ 0:9 and b ¼ 1:3.
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
where qf, qp are the volume flow rate for the fluid phase and
particulate phase.The instantaneous flux in the laboratory frame is
Qf ¼ ð1� CÞZ h
0
Z 1
0
ðuf þ 1Þdydz ¼ qf þ hð1� CÞ: ð25Þ
Qp ¼ C
Z h
0
Z 1
0
ðup þ 1Þdydz ¼ qp þ hC: ð26Þ
The average volume flow rate over one period T ¼ kc
� of
the peristaltic wave is defined as
�Qðx; tÞ ¼ 1
T
Z T
0
ðQf þQpÞdt ¼ qf þ qp þ 1: ð27Þ
0 0.2 0.4 0.6 0.8 1-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
Q
Δp
φ = 0.0
φ = 0.2
φ = 0.4
φ = 0.6
Figure 4 Variation of Dp with Q for different values of / at
K ¼ 0:005, M ¼ 0:8, C ¼ 0:5, lslip ¼ 0:3, l ¼ 0:9 and b ¼ 1:3.
0 0.2 0.4 0.6 0.8 1-40
-35
-30
-25
-20
-15
-10
-5
0
x
xd/pd
C = 0.01
C = 0.02
C = 0.03
C = 0.04
Figure 5 Variation of dp=dx with x for different values of C at
Q ¼ 2:5 , lslip ¼ 0:3, M ¼ 0:2, b ¼ 1, K ¼ 0:005, l ¼ 0:9 and
/ ¼ 0:6.
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria
0 0.2 0.4 0.6 0.8 1-5
-4.5
-4
-3.5
-3
-2.5
-2
-1.5
-1
-0.5
x
xd/pdM = 0.5
M = 1.0
M = 1.5
M = 2.0
Figure 6 Variation of dp=dx with x for different values of M at Q ¼ 2:0, lslip ¼ 0:4, C ¼ 0:5, b ¼ 1:3, K ¼ 0:05, l ¼ 0:9 and / ¼ 0:6.
0 0.2 0.4 0.6 0.8 1-1.3
-1.2
-1.1
-1
-0.9
-0.8
-0.7
-0.6
-0.5
x
xd/pd
lslip
= 0.1
lslip
= 0.3
lslip
= 0.5
lslip
= 0.7
Figure 7 Variation of dp=dx with x for different values of lslip at
Q ¼ 2:0, M ¼ 0:4, C ¼ 0:5, b ¼ 1:3, K ¼ 0:05, l ¼ 0:9 and
/ ¼ 0:6.
0 0.2 0.4 0.6 0.8 1-2
-1.5
-1
-0.5
x
xd/pd
Q = 1.0
Q = 1.2
Q = 1.4
Q = 1.6
Figure 8 Variation of dp=dx with x for different values of Q at
b ¼ 1, M ¼ 1, C ¼ 0:5, l ¼ 0:9, K ¼ 0:05, / ¼ 0:6 and lslip ¼ 0:6.
Particulate suspension slip flow 5
where �Q the mixture average flux of fluid phase and particulatephase.
From Eqs. (25) and (26) with (29), the pressure gradient isobtained as
dp
dx¼ 3ð�1þ CÞhMð �Q� 1þ hÞ�l
3Ch2 þ h4M� 6MbP1
n¼1sinhðan=bÞ
a5nðcoshðan=bÞþlslipanb sinhðan=bÞÞ
: ð28Þ
Integration of Eq. (22) over one wavelength yields
Dp¼10
dp
dxdx¼
Z 1
0
3ð�1þCÞhMð �Q�1þhÞ�l3Ch2þh4M�6Mb
P1n¼1
sinhðan=bÞa5nðcoshðan=bÞþlslip
anb sinhðan=bÞÞ
8<:
9=;dx:
ð29Þ
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
It is noticed that when C? 0, �l ! 0 and lslip ? 0 the solu-tions of Reddy et al. [19] are recovered as a special case of our
problem. Also, in the limit b ! 0 (keeping a fixed and d ! 1),the rectangular duct reduces to a two dimensional channel andresults of Shapiro [1] are recovered.
4. Numerical results and discussion
In this section, the graphical results of the problem under con-
sideration are discussed. The expression for the pressure riseand pressure gradient is calculated using a mathematics soft-ware Mathematica.
The variation of pressure rise with volume flow rate Q fordifferent values of M, lslip, K and / is shown in Figs. 1–4. Itis observed from Fig. 1 that in the peristaltic pumping
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria
Figure 9 Streamlines for different values of M. (a) for M ¼ 0:55 and (b) for M ¼ 0:6. The other parameters are Q ¼ 1, b ¼ 1:3,
lslip ¼ 0:3, K ¼ 0:04, C ¼ 0:5, l ¼ 0:9, y ¼ 0 and / ¼ 0:4.
Figure 10 Streamlines for different values of C. (a) for C ¼ 0:1 and (b) for C ¼ 0:15. The other parameters are Q ¼ 1, b ¼ 1:3, lslip ¼ 0:3,
K ¼ 0:04, M ¼ 0:5, l ¼ 0:9, y ¼ 0 and / ¼ 0:6.
Figure 11 Streamlines for different values of K. (a) for K ¼ 0:07 and (b) for K ¼ 0:01. The other parameters are Q ¼ 1, b ¼ 1:3;
lslip ¼ 0:3 C ¼ 0:5; M ¼ 0:5; l ¼ 0:9; y ¼ 0 and / ¼ 0:6.
6 S. Akram et al.
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induced by peristaltic waves in a rectangular duct: Effect of lateral walls, AlexandriaEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
Figure 12 Streamlines for different wave shape The other parameters are Q ¼ 1, b ¼ 1:3; lslip ¼ 0:3, C ¼ 0:1; K ¼ 0:04;M ¼ 0:5; l ¼ 0:9,
y ¼ 0 and / ¼ 0:6.
Particulate suspension slip flow 7
ðDp > 0;Q > 0Þ and retrograde pumping ðDp > 0;Q < 0Þregions the pumping rate increases with an increase inM, whilein the copumping region ðDp < 0;Q > 0Þ the behavior is quiteopposite, and here pumping rate decreases with an increase in
M. It is depicted from Fig. 2 that in the peristaltic pumpingðDp > 0;Q > 0Þ and retrograde pumping ðDp > 0;Q < 0Þregions the pressure rise remains constant in these regions,
while in the copumping region ðDp < 0;Q > 0Þ the pressurerise decreases with an increase in lslip. Figs. 3 and 4 show the
variation of pressure rise with volume flow rate Q. It isobserved from Figs. 3 and 4 that the pressure rise decreasesin all the regions with an increase in K and pressure rise
increases with an increase in /. The pressure gradient for dif-ferent values of C, M, lslip and Q are plotted in Figs. 5–8. It isdepicted from Figs. 5–8 that for x 2 ½0; 0:2� and x 2 ½0:8; 1�, thepressure gradient is small i.e., the flow can easily pass withoutimposition of a large pressure gradient, while in the regionx 2 ½0:2; 0:8�, the pressure gradient increases with an increase
in C, and decreases with an increase inM and Q, so much pres-sure gradient is required to maintain the flux to pass. It is alsoobserved that the pressure rise increases in the upper half of
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
the channel and decreases in the lower half of the channel withthe increase in lslip parameter. Streamlines for different valuesof M, C and K are shown in Figs. 9–11. It is observed fromFig. 9 that the size of the trapping bolus increases with an
increase in M. It is depicted from Fig. 10 that with the increasein C the size and number of the trapped bolus decrease. Fig. 11shows that the size of the trapped bolus decreases with the
decrease in K. Fig. 12 shows the streamlines for different waveform. Figure a is for sinusoidal wave, figure b for multisinu-soidal wave, figure c for triangular wave and figure d for trape-
zoidal wave.
References
[1] M.Y. Jaffrin, A.H. Shapiro, Peristaltic pumping, Annu. Rev.
Fluid Mech. 3 (1971) 13–36.
[2] V.P. Srivastava, M.A. Saxena, Two-fluid model of non-
Newtonian blood flow induced by peristaltic waves, Rheol.
Acta 34 (1995) 406–414.
[3] V.P. Srivastava, Effects of an inserted endoscope on chyme
movement in small intestine, Appl. Appl. Math. 2 (2007) 79–91.
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria
8 S. Akram et al.
[4] V.P. Srivastava, Particle-fluid suspension flow induced by
peristaltic waves in a circular cylindrical tube, Bull. Calcutta
Math. Soc. 94 (2002) 167–184.
[5] Kh.S. Mekheimer, Y. Abd Elmaboud, Peristaltic flow through a
porous medium in an annulus, Appl. An Endosc. 2 (1) (2008)
103.
[6] Kh.S. Mekheimer, Nonlinear peristaltic transport through a
porous medium in an inclined planner channel, J. Por. Media 6
(2003) 189–201.
[7] Kh.S. Mekheimer, S.Z.A. Husseny, Y. Abd Elmaboud, Effects
of heat transfer and space porosity on peristaltic flow in a
vertical asymmetric channel, Numer. Methods Partial Differ.
Eqs. (2009), http://dx.doi.org/10.1002/num.20451.
[8] T. Hayat, Y. Wang, A.M. Siddiqui, K. Hutter, Peristaltic
motion of a Segalman fluid in a planar channel, Math. Prob.
Eng. 1 (2003) 1–23.
[9] L.M. Srivastava, V.P. Srivastava, Peristaltic transport of a
particle-fluid suspension, J. Biomech. Eng. 111 (1989) 158.
[10] M.M. Bhattia, M. Ali Abbas, M.M. Rashidi, Combine effects of
Magnetohydrodynamics (MHD) and partial slip on peristaltic
blood flow of Ree-Eyring fluid with wall properties, Eng. Sci.
Technol., Int. J. (2016), http://dx.doi.org/10.1016/j.jestch.2016.
05.004.
[11] M.M. Bhatti, A. Zeeshan, N. Ijaz, Slip effects and endoscopy
analysis on blood flow of particle-fluid suspension induced by
peristaltic wave, J. Mol. Liq. 218 (2016) 240–245.
[12] S. Nadeem, A. Riaz, R. Ellahi, Noreen Sher Akbar, A. Zeeshan,
Heat and mass transfer analysis of peristaltic flow of nanofluid
Please cite this article in press as: S. Akram et al., Particulate suspension slip flow induEng. J. (2016), http://dx.doi.org/10.1016/j.aej.2016.09.011
in a vertical rectangular duct by using the optimized series
solution and genetic algorithm, J. Comput. Theor. Nanosci. 11
(4) (2014) 1133–1149.
[13] A. Riaz, S. Nadeem, R. Ellahi, A. Zeeshan, Exact solution for
peristaltic flow of Jeffrey fluid model in a three dimensional
rectangular duct having slip at the walls, Appl. Bion. Biomech.
11 (2014) 81–90.
[14] L.M. Srivastava, V.P. Srivastava, Peristaltic transport of a
particle-fluid suspension, J. Biomech. Eng. 111 (1989) 157.
[15] J.C. Misra, S.K. Pandey, Peristaltic transport of a particle-fluid
suspension in a cylindrical tube, Comput. Math. Appl. 28 (1994)
131.
[16] V.P. Srivastava, L.M. Srivastava, Influence of wall elasticity and
poiseuille flow on peristaltic induced flow of a particle-fluid
mixture, Int. J. Eng. Sci. 35 (1997) 1359.
[17] V.P. Srivastava, M. Saxena, Particulate suspension flow induced
by sinusoidal peristaltic waves, Jpn. J. Appl. Phys. 36 (1997) 385.
[18] Kh.S. Mekheimer, Y. Abd Elmaboud, Peristaltic transport of a
particle–fluid suspension through a uniform and non-uniform
annulus, Appl. Bion. Biomech. 5 (2) (2008) 47.
[19] M.V. Subba Reddy, M. Mishra, S. Sreenadh, A.R. Rao,
Influence of lateral walls on peristaltic flow in a rectangular
duct, J. Fluids Eng. 127 (2005) 824.
[20] S. Nadeem, Safia Akram, Peristaltic flow of a Jeffrey fluid in a
rectangular duct, Nonlin. Anal.: Real World Appl. 11 (2010)
4238.
[21] S.E. Charm, G.S. Kurland, Blood Flow and Micro-circulation,
John Wiley, New York, 1974.
ced by peristaltic waves in a rectangular duct: Effect of lateral walls, Alexandria