# Effects of Heat and Mass Transfer on MHD Peristaltic Flow of a Non-Newtonian Fluid through a Porous Medium between Two Coaxial Cylinders

Authors: Abeer A. Shaaban; Mohamed Y. Abou-zeid
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We investigated the influence of heat and mass transfer on the peristaltic flow of magnetohydrodynamic Eyring-Powell fluid under low Reynolds number and long-wavelength approximation. The fluid flows between two infinite cylinders; the inner tube is uniform, rigid, and rest, while the outer flexible tube has a sinusoidal wave traveling down its wall. The governing equations are solved numerically using finite-difference technique. The velocity, temperature, and concentration distribution are obtained. The features of flow characteristics are analyzed by plotting graphs and discussed in detail.


---

## Body

## 1. Introduction

Expansion and contraction of an extensible tube in a fluid generate progressive waves which propagate along the length of the tube, mixing and transporting the fluid in the direction of wave propagation. This phenomenon is known as peristalsis. It is an inherent property of many tubular organs of the human body. In some biomedical instruments, such as heart-lung machines, peristaltic motion is used to pump blood and other biological fluids [1]. In addition, peristaltic pumping occurs in many practical applications involving biomedical systems. It has now been accepted that most of the physiological fluids behave like a non-Newtonian fluids. Raju and Devanathan [2] were probably the first to consider this aspect. The analysis of the mechanisms responsible for peristaltic transport for non-Newtonian fluid has been studied by many authors.Mernone et al. [3] considered the peristaltic flow of rheologically complex physiological fluids when modelled by a non-Newtonian Casson fluid  in a two-dimensional channel. Haroun [4] studied the effect of a third-order fluid on the peristaltic transport in an asymmetric channel. In his study, the wavelength of the peristaltic waves is assumed to be large compared to the varying channel width, whereas the wave amplitudes need not be small compared to the varying channel width. Eldabe et al. [5] analyzed the incompressible flow of electrically conducting biviscosity fluid through an axisymmetric nonuniform tube with a sinusoidal wave under the considerations of long wavelength and low Reynolds number. The effect of rotation, porous medium, and magnetic field on peristaltic transport of a Jeffrey fluid in tube is studied by Mahmoud [6]. Another work dealing with non-Newtonian peristaltic flow is mentioned in [7–11].In all the abovementioned studies, no porous medium has been taken into account. But it is well known that flow through a porous medium has practical applications especially in geophysical fluid dynamics. Examples of natural porous media are beach sand, sandstone, limestone, rye bread, wood, the human lung, bile duct, gall bladder with stones and in small blood vessels. In the arterial systems of humans or animals, it is quite common to find localized narrowings, commonly caused by in travascular plaques [12]. Afifi [13] was the first who studied the peristaltic flow through a porous medium.  The peristaltic  motion  of  a generalized  Newtonian  fluid  through  a  porous  medium is studied by Elshehawey et al. [14]. Mekheimer [15] studied the nonlinear  peristaltic  transport  through  a porous medium in an inclined planar channel taking the gravity effect on pumping characteristics. Such analysis is of great value in medical research. Hayat et al.  [16]   investigated   the   peristaltic   transport   of   electrically  conducting Maxwell fluid through a  porous medium in a  planar  channel.Magnetohydrodynamics (MHD) is the science which deals with the motion of a highly conducting fluids in the presence of a magnetic field. The motion of the conducting fluid across the magnetic field generates electric currents which change the magnetic field, and the action of the magnetic field on these currents gives rise to mechanical forces which modify the flow of the fluid [17]. The study of peristaltic flow with heat and mass transfer for an electrically conducting  fluid past a porous plate under the influence of a magnetic field has attracted the interest of many investigators in view of  its  applications in many engineering problems such as MHD generator, plasma studies, nuclear reactors, oil exploration, geothermal energy extractions, and the boundary layer control in the field of aerodynamics [18]. Vajravelu et al. [19] examined the flow through vertical porous tube with peristalsis and heat transfer. Srinivas and Kothandapani [20] have studied the influence of heat and mass transfer on MHD peristaltic flow through a porous space with compliant walls. The peristaltic transport of Newtonian fluid in a vertical asymmetric channel with heat transfer and porous medium is studied by Srinivas and Gayathri [21]. For the benefit of the readers in this direction you can take these additional papers ([22–25]) into account.With the above discussion in mind, we propose to study the fluid mechanics effects of MHD peristaltic transport through a porous media in a gap between two coaxial tubes, filled with an incompressible non-Newtonian (Eyring-Powell) fluid. The effects of body temperature and concentration are taken into consideration. The inner tube is rigid and the outer one has a sinusoidal wave traveling down its wall. The flow analysis is developed in a wave frame of reference moving with the velocity of the wave. Our problem can be considered as a mathematical model for the blood flow with catheterization. The problem is solved numerically for the velocity, temperature, and concentration. The effects of various emerging parameters on the flow, temperature, and concentration distributions are shown and discussed with the help of graphs.


## 2. Formulation of the Problem

Let us consider the MHD flow with heat and mass transfer of an incompressible electrically conducting non-Newtonian obeying Eyring-Powell model between two coaxial tubes through a porous medium. The inner tube is rigid and uniform, while the outer tube has a sinusoidal wave traveling down its wall. We use cylindrical coordinate system(r,θ,z) with r in the radial direction, z lies along the center line of the inner and outer tubes, the geometry of the wall surfaces is described in (Figure 1), and the equations for the radii are
(1)r1=d2,r2=H=d+bsin2πλ(z-ct),
in which d signifies the radius of the outer tube, b indicates the wave amplitude, λ is the wavelength, c is the propagation velocity along z direction, and t is the time. The fluid is subjected to a constant magnetic field B0 that acts in θ direction. Induced magnetic field and external electric field are neglected.Figure 1The equations that govern the flow are the balance of mass:(2)∇·V_=0,
the equation of momentum:
(3)ρ(V_·∇V_)=-∇P+∇·τ_+J_×B_-μKV_,
Maxwell's equations: (4a)∇·B_=0,(4b)∇×B_=μmJ_,(4c)∇×E_=-∂B_∂t,(4d)J_=σ(E_+V_×B_),the equation of energy:
(5)ρCp(V_·∇T)=kc∇2T+Γij∂Vi∂xj+ρDmKTCs∇2C,
and the equation of concentration:
(6)V_·∇C=Dm∇2C-A¯(C-C∞)m+DmKTTm∇2T.
The stress  is defined by
(7)τ_=μ∇V_+1a2sinh-1(1a1∇V_),
where ρ is the density,  V_ is the fluid velocity and P are the pressure, τ_ is the stress tensor, μ is the plastic viscosity of the fluid, kp is the permeability of the porous medium, B_ is total magnetic field vector, J_ is current density vector, T is the fluid temperature, C is the fluid concentration, Cp is the specific heat at constant pressure, kc is the thermal conductivity of the fluid, Dm is the coefficient of mass diffusivity, KT is the thermal diffusion ratio, Tm and Cm are the mean fluid temperature and concentration, Cs is the concentration susceptibility, q_r is the radiative heat flux, E_ is the electric field vector, σ is the electrical conductivity of the fluid, μm is the magnetic permeability, A¯ is the reaction rate constant, m is the reaction order, and a1 and a2 are characteristics of Eyring-Powell model.Since the flow parameters are independent of the azimuthal coordinateθ; hence, the velocity is given by V_=(u,0,w). A uniform magnetic field B_=(0,B0,0) is applied in a direction perpendicular to the fluid flow. The governing equations (2)–(6) become
(8)∂u∂r+ur+∂w∂z=0,ρ(u∂u∂r+w∂u∂z)=-∂P∂r+∂τrr∂r+τrrr+∂τrz∂z-(σB02+μK)u,ρ(u∂w∂r+w∂w∂z)=-∂P∂z+∂τrz∂r+τrzr+∂τzz∂z-(σB02+μK)w,ρCp(u∂T∂r+w∂T∂z)=kc(∂2T∂r2+1r∂T∂r+∂2T∂z2)+τrr∂u∂r+τrz∂u∂z+τzr∂w∂r+τzz∂w∂z+ρDmKTCs(∂2C∂r2+1r∂C∂r+∂2C∂z2)u∂C∂r+w∂C∂z=Dm(∂2C∂r2+1r∂C∂r+∂2C∂z2)+DmKTTm(∂2T∂r2+1r∂T∂r+∂2T∂z2)-A¯(C-Cm)m.
The boundary conditions are given by
(9)u=0,w=0,T=C=0atr=r1.u=-c∂H∂z,w=-c,T=Tw,C=Cwatr=r2,
where Tw and Cw are the uniform temperature and concentration at the outer tube.The appropriate nondimensional variables for the flow are defined as(10)r*=rd,z*=zλ,u*=λcdu,w*=wc,P*=d2λcμP,t*=cλt,T*=T-TmTw-Tm,C*=C-CmCw-Cm,h=Hd,ε=bd,δ=dλ,α=1μa1a2,γ=(ca1d)2,γ¯=γα.
In terms of these variables and dropping the star mark for simplicity, (8) becomes
(11)∂u∂r+ur+∂w∂z=0,Reδ3(u∂u∂r+w∂u∂z)=-∂P∂r-δ3Dau=.+δ2(∂2u∂r2+1r∂u∂r-ur2+δ2∂2u∂z2)=.×(2δ2(∂w∂z)2)+1)-1/21+α(1a12(2δ2(∂u∂r)2+2δ2(ur)2=.+(δ∂u∂z+∂w∂r)2.+2δ2(∂w∂z)2)+11a12(2δ2(∂u∂r)2+2δ2(ur)2)-1/21+α(1a12(2δ2(∂u∂r)2+2δ2(ur)2),Reδ(u∂w∂r+w∂w∂z)=-∂P∂z-(1Da+M)w+(∂2w∂r2+1r∂w∂r+δ2∂2w∂z2)×(2δ2(∂w∂z)2)+1)-1/21+α(1a12(2δ2(∂u∂r)2+2δ2(ur)2+(δ∂u∂z+∂w∂r)2Z.+2δ2(∂w∂z)2)+11a12(2δ2(∂u∂r)2+2δ2(ur)2)-1/2),Reδ(u∂T∂r+w∂T∂z)=1Pr(∂2T∂r2+1r∂T∂r+δ2∂2T∂z2).+Ec(2δ2(∂u∂r)2+2δ2(ur)2+(δ∂u∂z+∂w∂r)2=.+2δ2(∂w∂z)2)+Ecγ¯[(δ∂u∂z+∂w∂r)sinh-1(γ(δ∂u∂z+∂w∂r))Z+δ(∂u∂r)sinh-1(δγ∂u∂r)Z+δ(ur)sinh-1(δγur)Z+δ(∂w∂z)sinh-1(δγ∂w∂z)(γ(δ∂u∂z+∂w∂r))]+Df(∂2C∂r2+1r∂C∂r+δ2∂2C∂z2),Reδ(v∂C∂r+u∂C∂z)=1Sc(∂2C∂r2+1r∂C∂r+δ2∂2C∂z2)+Sr(∂2T∂r2+1r∂T∂r+δ2∂2T∂z2)-ACm.
Here, Re=ρcd/μ is Reynolds number, Da=K/d2 is Darcy number, M=σB02d2/μ is the magnetic field parameter, Pr=μCP/kc is Prandtl number, Ec=c2/CP(Tw-Tm) is Eckert number, Sc=μ/ρDm is Schmidt number, Sr=DmkTρ(Tw-Tm)/μTm(Cw-Cm) is Soret number, A=A¯re(Cw-Cm)m-1d2ρ/μ is the reaction rate parameter, Df=DmkT(Cw-Cm)/2νCsCp(Tw-Tm) is Dufour number, and α=1/μa1a2 is Eyring Powell parameter.Thus, the boundary conditions (9) in their dimensionless form read
(12)u=0,w=0,T=C=0atr=r1=12,u=-∂h∂z,w=-1,T=C=1atr=h=1+εsin2πz.After applying the long-wavelength approximation to our analysis as described in [5], then the field (11) now gives
(13)∂u∂r+ur+∂w∂z=0,(14)∂P∂r=0,(15)(∂2w∂r2+1r∂w∂r)×(1+α(1/a12)(∂w/∂r)2+1)-(1Da+M)w=∂P∂z,(16)1Pr(∂2T∂r2+1r∂T∂r)+Ec(∂w∂r)2+Ecγ¯(∂w∂r)×sinh-1(γ∂w∂r)+Df(∂2C∂r2+1r∂C∂r)=0,(17)∂2C∂r2+1r∂C∂r+Sc·Sr(∂2T∂r2+1r∂T∂r)-A·Sc·Cm=0.
Equation (14) indicates that P is not function of r. Hence P is function of z and t only.


## 3. Method of Solution

We shall solve the system of non-linear ordinary differential equations numerically using the finite difference technique which is discussed by Soundalgekar and Ganesan [26], and so the system of (15)–(17) yields
(18)(wi+1-2wi+wi-1h2+1ihwi+1-wi-12h)×(1+α(1/a2)((wi+1-wi-1)/2h)2+1)-(1Da+M)wi=∂P∂z,1Pr(Ti+1-2Ti+Ti-1h2+1ihTi+1-Ti-12h)+Ec(wi+1-wi-12h)2+Ecγ¯(wi+1-wi-12h)×sinh-1(γ(wi+1-wi-12h))+Df(Ci+1-2Ci+Ci-1h2+1ihCi+1-Ci-12h)=0,(Ci+1-2Ci+Ci-1h2+1ihCi+1-Ci-12h)+Sc·Sr(Ti+1-2Ti+Ti-1h2+1ihTi+1-Ti-12h)-A·Sc·(Ci)m=0.
The boundary conditions (12) transform into
(19)wi=0,Ti=Ci=0atr=r1=12,wi=-1,Ti=Ci=1atr=h=1+εsin2πz.


## 4. Results and Discussion

To discuss the effects of various parameters involved in the problem such as Prandtl numberPr, magnetic field parameter M, Eckert number Ec, Darcy number Da, Dufour number Df, reaction rate parameter A, Soret number Sr, Schmidt number Sc, Eyring Powell parameter α, and the pressure gradient dP/dz, on the solution of the considered problem, numerical results are calculated by solving (18) with boundary conditions (19) using finite difference technique. The distributions of axial velocity w, the temperature T, and the concentration C are illustrated in Figures 2–25.Figure 2
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, M=2, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 3
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 4
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 5
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Df=0.07, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 6
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Da=0.5, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 7
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Da=0.5, α=0.3, a=0.4, Sc=0.5, Df=0.07, and m=2.Figure 8
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, A=0.1, Da=0.5, α=0.3, a=0.4, Sc=0.5, Df=0.07, and m=2.Figure 9
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, A=0.1, Da=0.5, α=0.3, a=0.4, Sr=2, Df=0.07, and m=2.Figure 10
Profiles of the dimensionless temperatureT(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 11
Profiles of the dimensionless temperatureT(r) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 12
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, M=2, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 13
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 14
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 15
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Df=0.07, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 16
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Da=0.5, α=0.3, a=0.4, Sc=0.5, A=0.1, and m=2.Figure 17
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, Sr=2, Da=0.5, α=0.3, a=0.4, Sc=0.5, Df=0.07, and m=2.Figure 18
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, A=0.1, Da=0.5, α=0.3, a=0.4, Sc=0.5, Df=0.07, and m=2.Figure 19
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Ec=0.3, A=0.1, Da=0.5, α=0.3, a=0.4, Sr=2, Df=0.07, and m=2.Figure 20
Profiles of the dimensionless concentrationC(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 21
Profiles of the dimensionless concentrationC(r) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 22
Profiles of the dimensionless velocityw(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 23
Profiles of the dimensionless velocityw(r) for a system have the particulars dP/dz=-30, Pr=3, M=2, Sr=2, Df=0.07, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 24
Profiles of the dimensionless velocityw(r) for a system have the particulars dP/dz=-30, Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.Figure 25
Profiles of the dimensionless velocityw(r) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, M=2, α=0.3, a=0.4, Ec=0.3, Sc=0.5, A=0.1, and m=2.The values ofT are plotted versus the radial coordinate in Figures 2–11. In Figure 2, the effect of Pr on the temperature distribution T is indicated, and this figure shows that the temperature increases with the increasing of Pr. Also the relation between T and r is approximately linear for small values of Pr (<0.7), but for large values of Pr (>0.7), T increases with increasing r till a finite value r (maximum value), after which T decreases with increasing r. The effect of M on the temperature distribution T is shown in Figure 3; it is clear from this figure that the temperature decreases with increasing M. It is also noted that the temperature increases with increasing r and reaches a maximum value (at a finite value of r: r=0.85) after which it decreases. Figure 4 gives the effect of Ec on the temperature distribution T; from this figure, we observe that T increases with increasing Ec, and although there is a small difference between the values of Ec, there is a remarkable difference between the four curves of T, which were taken at Ec = 0.2, 0.3, 0.4, and 0.5 when the other parameters are constant, that is, Ec affects T strongly.  The effect of Da on the temperature distribution is described in Figure 5, it is clear that the temperature increases with increasing Da but the effect of Da on T disappears at the range 0.96<r<1. Also, it is clear that the relation between T and r is approximately linear for small values of Da but for large values of Da there is a normal curves. Figure 6 shows the effect of Df on the temperature distribution, and it is seen from this figure that the temperature increases with increasing Df, and, in all curves, the maximum value of T occurs at 0.85<r<0.9, and it increases as r increases. In Figure 7, the relation between T and r for various values of A is obtained. From this figure, it is clear that the temperature increases with increasing A, and, from this figure, we can see that although there is a big difference between the values of A, there is a very small difference between the four curves of T, which were taken at A = 0.1, 100, 200, and 300 when the other parameters are constant. Figure 8 shows the effect of Sr on the temperature distribution; it is noted that the temperature increases with increasing Sr. It is also noted that the effect of Sr on T disappears near the boundaries of the tubes and all curves have a maximum value of T in a finite range of r namely 0.88<r<0.9. Figures 9 and 10 show that the effect of both Sc and α on the temperature distribution T is similar to the effect of Sr on T given in Figure 8. The effect of dP/dz on the temperature distribution is indicated in Figure 11. It is seen from this figure that the temperature decreases with increasing dP/dz, and all carves have a maximum value as in Figure 8.Figures12–21 show the effect of the various parameters on the concentration distribution. Figure 12 gives the effect of Pr on the concentration distribution. This figure shows that the concentration decreases with the increase of Pr. Also Pr affects the relation between C and r, and this relation is approximately linear at the values 0.7 and 1 of Pr, but, for large values of Pr (≥3), C decreases with increasing r till a finite value of r (minimum values), after which C increases with increasing r. Figure 13 shows the effect of M on the concentration distribution, it is observed that the concentration increases with increasing M and all curves have a minimum value at the range 0.98<r<1. The effect of Ec on the concentration distribution is illustrated in Figure 14, and we can see from this figure that although there is a small difference between the values of Ec, there is a remarkable difference between the four curves of C. Figure 15 describes the relation between C and r for various values of Da, and it is clear that the concentration decreases with increasing Da but the effect of Da on C disappears at the range 0.96<r<1. Also, it is clear that the relation between C and r is approximately linear for small values of Da, but for large values of Da, there are normal curves. Figure 16 shows the effect of Df on the concentration distribution. It is seen from this figure that the concentration decreases with increasing Df. Also, the concentration C decreases with increasing r in the range from 0 to 0.7 but after which it increases with increasing r. Figure 17 describes the relation between C and r for various values of A, it is obvious that the effect of A on C is similar to the effect of Df on C, with the only difference that the obtained curve is closer to those obtained in Figure 16. The effect of Sr on the concentration distribution is shown in Figure 18; from this figure, we can see that C decreases with increasing Sr, and, from this figure, we can see that the effect of Sr on C disappears in the range 0.98<r<1. Also,  the relation between C and r is approximately linear at the value namely, 1.5 of Sr. All the curves of C were taken at Sr = 1.5, 1.8, 2, and 2.5 when the other parameters are constant. Figure 19 gives the effect of Sc on the concentration distribution, and it is clear that C decreases with increasing Sc; also Sc affects the relation between C and r. This relation is approximately linear for small values of Sc, but, for large values of Sc, C decreases with increasing r till a finite value of r (minimum value), after which it increases. Figure 20 describes the relation between C and r for various values of α. It is clear from this figure that the concentration decreases with increasing α. From all curves in this figure, we obtained that the effect of α on C disappears at the range 0.5<r<0.55 and 0.96<r<1. All curves were taken at α = 0.1, 0.3, 0.5, and 0.9 when the other parameters are constant. Figure 21 shows the effect of dP/dz on the concentration distribution. The calculations were obtained at dP/dz = −35, −30, −25, and −20 when the other parameters are constant. C increases with increasing dP/dz. It is also noted from all curves in this figure that the concentration decreases by increasing r at the range 0.5<r<0.65 but after that it increases with increasing r.Figures22–25 show the effect of the various parameters on the velocity distribution. Figure 22 gives the effect of M on the velocity distribution. This figure shows that the velocity decreases by increasing M at the range 0.5<r<0.84 and after that the velocity increases by increasing M. Also, it is clear that, in all curves in this figure, the maximum value of w occurs at 0.65<r<0.7. Figure 23 illustrates the effect of Da on the velocity distribution. It is clear that w increases by increasing Da in the range 0.5<r<0.84 and w decreases by increasing Da in the range 0.84<r<1. Also, from this figure, we can see that although there is a big difference between the values of Da, there is a small difference between the four curves of w which were taken at Da = 0.005, 0.01, 0.09, and 0.5, when the other parameters are constant. Figure 24 describes the relation between w and r for various values of α. From this figure, it is clear that the velocity decreases with increasing α, but the effect of α on w disappears in the range 0.96<r<1. All curves in this figure were taken at α = 0.1, 0.3, 0.5, and 0.9 when the other parameters are constant. Figure 25 shows the effect of dP/dz on the velocity distribution. The calculation was obtained at dP/dz = −35, −30, −25, and −20, and the other parameters are constant. The velocity decreases with increasing dP/dz and in all curves the maximum value of w occurs in the range 0.65<r<0.7.


## 5. Conclusion

In this paper, we studied the effects of the physical parameters of the considered problem on peristaltic transport in the gap between two coaxial tubes, filled with an incompressible non-Newtonian (Eyring-Powell) fluid, and considered the effects of body temperature and concentration. The inner tube is rigid and the outer one has a sinusoidal wave traveling down its wall. The present analysis can serve in understanding the mechanism of many physiological flows. Physically, our model corresponds to the transport of the gastric juice in the small intestine when an endoscope is inserted through it. The system is solved numerically by finite difference technique. The effects of various emerging parameters on the flow, the temperature, and the concentration distributions are shown and discussed with the help of graphs. The main findings can be summarized as follows.(1)
The velocityw rise decreases with the increase of each of M, α and dP/dz, whereas it increases as Da increases. But, all the other parameters such as, Pr, Ec, Df, A, Sr, and Sc do not affect on the velocity (the flow of fluid).(2)
The temperatureT rise increases with the increase of each of Pr, Ec, Da, Df, A, Sr, Sc, and α, whereas it decreases as M and dP/dz increase.(3)
The concentrationC rise increases with the increase each of M and dP/dz, whereas it decreases as Pr, Ec, Da, Df, A, Sr, Sc, and α increase.(4)
The temperature has the opposite behavior compared to the concentration rise behavior at all different parameters in this problem.


## 6. Applications

(
1
) Gastrointestinal motility is the function of gastrointestinal smooth muscle. It is controlled by both the intrinsic and extrinsic nerves of the gastrointestinal tract and, to a lesser degree, the gastrointestinal hormones. Therefore, any abnormality of the above factors, theoretically, can cause gastrointestinal dysmotility. In a clinical situation, commonly seen is gastrointestinal dysmotility caused by either smooth muscle or intrinsic and extrinsic nerves dysfunction.Diseases that cause smooth muscle dysfunction include familial visceral myopathies, nonfamilial visceral myopathies, collagen disease, muscular dystrophies, amyloidosis, and thyroid disease. Diseases that cause enteric nerve dysfunction include familial visceral neuropathies, nonfamilial visceral neuropathies, diabetes mellitus, Chagas’ disease, ganglioneuromatosis of the intestine, visceral neuropathy of carcinomatosis, Parkinson’s disease.(
2
) Gastroparesis syndrome is a recognized complication of longstanding diabetes mellitus and is attributed to reduced gastric contractility due to “autovagotomy.” However, motor abnormalities associated with this syndrome may not be limited to the stomach.  Abnormal intestinal manometric patterns were observed in twelve out of the fourteen patients.Our findings show that the small intestine is frequently affected in patients with diabetic gastroparesis, and that the motility disorder both in the stomach and the small bowel is not invariably of a “paretic” type. The occurrence of incoordinated intestinal long bursts and continuous antral activity suggests that disturbed sympathetic innervations participate in the aetiopathogenesis of their upper-gut dysfunction.



---

*Source: 819683-2013-06-26.xml*