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

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.

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 (𝑟,𝜃,𝑧) with 𝑟 in the radial direction, 𝑧 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 𝑟₁=𝑑2,𝑟₂=𝐻=𝑑+𝑏sin2𝜋𝜆(𝑧−𝑐𝑡), (1) in which 𝑑 signifies the radius of the outer tube, 𝑏 indicates the wave amplitude, 𝜆 is the wavelength, 𝑐 is the propagation velocity along 𝑧 direction, and 𝑡 is the time. The fluid is subjected to a constant magnetic field 𝐵₀ that acts in 𝜃 direction. Induced magnetic field and external electric field are neglected.

Figure 1 

The equations that govern the flow are the balance of mass: ∇⋅𝑉=0, (2) the equation of momentum: 𝜌𝑉⋅∇𝑉=−∇𝑃+∇⋅𝜏+𝐽×𝐵−𝜇𝐾𝑉, (3) Maxwell's equations: ∇⋅𝐵=0,∇×𝐵=𝜇_𝑚𝐽,∇×𝐸=−𝜕𝐵,𝐽𝜕𝑡=𝜎𝐸+𝑉×𝐵, (4a)(4b)(4c)(4d)the equation of energy: 𝜌𝐶_𝑝𝑉⋅∇𝑇=𝑘_𝑐∇²𝑇+Γ_𝑖𝑗𝜕𝑉_𝑖𝜕𝑥_𝑗+𝜌𝐷_𝑚𝐾_𝑇𝐶_𝑠∇²𝐶, (5) and the equation of concentration: 𝑉⋅∇𝐶=𝐷_𝑚∇²𝐶−𝐴𝐶−𝐶_∞^𝑚+𝐷_𝑚𝐾_𝑇𝑇_𝑚∇²𝑇. (6) The stress is defined by 𝜏=𝜇∇𝑉+1𝑎₂sinh⁻¹1𝑎₁∇𝑉, (7) where 𝜌 is the density, 𝑉 is the fluid velocity and 𝑃 are the pressure, 𝜏 is the stress tensor, 𝜇 is the plastic viscosity of the fluid, 𝑘_𝑝 is the permeability of the porous medium, 𝐵 is total magnetic field vector, 𝐽 is current density vector, 𝑇 is the fluid temperature, 𝐶 is the fluid concentration, 𝐶_𝑝 is the specific heat at constant pressure, 𝑘_𝑐 is the thermal conductivity of the fluid, 𝐷_𝑚 is the coefficient of mass diffusivity, 𝐾_𝑇 is the thermal diffusion ratio, 𝑇_𝑚 and 𝐶_𝑚 are the mean fluid temperature and concentration, 𝐶_𝑠 is the concentration susceptibility, 𝑞_𝑟 is the radiative heat flux, 𝐸 is the electric field vector, 𝜎 is the electrical conductivity of the fluid, 𝜇_𝑚 is the magnetic permeability, 𝐴 is the reaction rate constant, 𝑚 is the reaction order, and 𝑎₁ and 𝑎₂ are characteristics of Eyring-Powell model.

Since the flow parameters are independent of the azimuthal coordinate 𝜃; hence, the velocity is given by 𝑉=(𝑢,0,𝑤). A uniform magnetic field 𝐵=(0,𝐵₀,0) is applied in a direction perpendicular to the fluid flow. The governing equations (2)–(6) become 𝜕𝑢+𝑢𝜕𝑟𝑟+𝜕𝑤𝜌𝑢𝜕𝑧=0,𝜕𝑢𝜕𝑟+𝑤𝜕𝑢𝜕𝑧=−𝜕𝑃+𝜕𝑟𝜕𝜏_𝑟𝑟+𝜏𝜕𝑟𝑟^𝑟𝑟+𝜕𝜏_𝑟𝑧−𝜕𝑧𝜎𝐵²₀+𝜇𝐾𝜌𝑢𝑢,𝜕𝑤𝜕𝑟+𝑤𝜕𝑤𝜕𝑧=−𝜕𝑃+𝜕𝑧𝜕𝜏_𝑟𝑧+𝜏𝜕𝑟𝑟^𝑟𝑧+𝜕𝜏_𝑧𝑧−𝜕𝑧𝜎𝐵²₀+𝜇𝐾𝑤,𝜌𝐶_𝑝𝑢𝜕𝑇𝜕𝑟+𝑤𝜕𝑇𝜕𝑧=𝑘_𝑐𝜕²𝑇𝜕𝑟²+1𝑟𝜕𝑇+𝜕𝜕𝑟²𝑇𝜕𝑧²+𝜏_𝑟𝑟𝜕𝑢𝜕𝑟+𝜏_𝑟𝑧𝜕𝑢𝜕𝑧+𝜏_𝑧𝑟𝜕𝑤𝜕𝑟+𝜏_𝑧𝑧𝜕𝑤+𝜕𝑧𝜌𝐷_𝑚𝐾_𝑇𝐶_𝑠𝜕²𝐶𝜕𝑟²+1𝑟𝜕𝐶+𝜕𝜕𝑟²𝐶𝜕𝑧²𝑢𝜕𝐶𝜕𝑟+𝑤𝜕𝐶𝜕𝑧=𝐷_𝑚𝜕²𝐶𝜕𝑟²+1𝑟𝜕𝐶+𝜕𝜕𝑟²𝐶𝜕𝑧²+𝐷_𝑚𝐾_𝑇𝑇_𝑚𝜕²𝑇𝜕𝑟²+1𝑟𝜕𝑇+𝜕𝜕𝑟²𝑇𝜕𝑧²−𝐴𝐶−𝐶_𝑚^𝑚. (8) The boundary conditions are given by 𝑢=0,𝑤=0,𝑇=𝐶=0at𝑟=𝑟₁.𝑢=−𝑐𝜕𝐻𝜕𝑧,𝑤=−𝑐,𝑇=𝑇_𝑤,𝐶=𝐶_𝑤at𝑟=𝑟₂, (9) where 𝑇_𝑤 and 𝐶_𝑤 are the uniform temperature and concentration at the outer tube.

The appropriate nondimensional variables for the flow are defined as 𝑟^∗=𝑟𝑑,𝑧^∗=𝑧𝜆,𝑢^∗=𝜆𝑐𝑑𝑢,𝑤^∗=𝑤𝑐,𝑃^∗=𝑑²𝜆𝑐𝜇𝑃,𝑡^∗=𝑐𝜆𝑡,𝑇^∗=𝑇−𝑇_𝑚𝑇_𝑤−𝑇_𝑚,𝐶^∗=𝐶−𝐶_𝑚𝐶_𝑤−𝐶_𝑚𝐻,ℎ=𝑑𝑏,𝜀=𝑑𝑑,𝛿=𝜆,1𝛼=𝜇𝑎₁𝑎₂𝑐,𝛾=𝑎₁𝑑²,√𝛾=𝛾𝛼. (10) In terms of these variables and dropping the star mark for simplicity, (8) becomes 𝜕𝑢+𝑢𝜕𝑟𝑟+𝜕𝑤𝜕𝑧=0,Re𝛿³𝑢𝜕𝑢𝜕𝑟+𝑤𝜕𝑢𝜕𝑧=−𝜕𝑃−𝛿𝜕𝑟³𝑢Da+𝛿²𝜕²𝑢𝜕𝑟²+1𝑟𝜕𝑢−𝑢𝜕𝑟𝑟²+𝛿²𝜕²𝑢𝜕𝑧²×11+𝛼𝑎₁²2𝛿²𝜕𝑢𝜕𝑟²+2𝛿²𝑢𝑟²+𝛿𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟²+2𝛿²𝜕𝑤𝜕𝑧²+1^−1/2,Re𝛿𝑢𝜕𝑤𝜕𝑟+𝑤𝜕𝑤𝜕𝑧=−𝜕𝑃−1𝜕𝑧𝑤+𝜕Da+𝑀₂𝑤𝜕𝑟²+1𝑟𝜕𝑤𝜕𝑟+𝛿²𝜕²𝑤𝜕𝑧²×11+𝛼𝑎₁²2𝛿²𝜕𝑢𝜕𝑟²+2𝛿²𝑢𝑟²+𝛿𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟²+2𝛿²𝜕𝑤𝜕𝑧²+1^−1/2,Re𝛿𝑢𝜕𝑇𝜕𝑟+𝑤𝜕𝑇=1𝜕𝑧𝜕Pr²𝑇𝜕𝑟²+1𝑟𝜕𝑇𝜕𝑟+𝛿²𝜕²𝑇𝜕𝑧²+Ec2𝛿²𝜕𝑢𝜕𝑟²+2𝛿²𝑢𝑟²+𝛿𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟²+2𝛿²𝜕𝑤𝜕𝑧²+Ec𝛾𝛿𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟sinh⁻¹√𝛾𝛿𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟+𝛿𝜕𝑢𝜕𝑟sinh⁻¹𝛿√𝛾𝜕𝑢𝑢𝜕𝑟+𝛿𝑟sinh⁻¹𝛿√𝛾𝑢𝑟+𝛿𝜕𝑤𝜕𝑧sinh⁻¹𝛿√𝛾𝜕𝑤𝜕𝜕𝑧+Df²𝐶𝜕𝑟²+1𝑟𝜕𝐶𝜕𝑟+𝛿²𝜕²𝐶𝜕𝑧²,Re𝛿𝑣𝜕𝐶𝜕𝑟+𝑢𝜕𝐶=1𝜕𝑧𝜕Sc²𝐶𝜕𝑟²+1𝑟𝜕𝐶𝜕𝑟+𝛿²𝜕²𝐶𝜕𝑧²𝜕+Sr²𝑇𝜕𝑟²+1𝑟𝜕𝑇𝜕𝑟+𝛿²𝜕²𝑇𝜕𝑧²−𝐴𝐶^𝑚. (11) Here, Re=𝜌𝑐𝑑/𝜇 is Reynolds number, Da=𝐾/𝑑² is Darcy number, 𝑀=𝜎𝐵²₀𝑑²/𝜇 is the magnetic field parameter, Pr=𝜇𝐶_𝑃/𝑘_𝑐 is Prandtl number, Ec=𝑐²/𝐶_𝑃(𝑇_𝑤−𝑇_𝑚) is Eckert number, Sc=𝜇/𝜌𝐷_𝑚 is Schmidt number, Sr=𝐷_𝑚𝑘_𝑇𝜌(𝑇_𝑤−𝑇_𝑚)/𝜇𝑇_𝑚(𝐶_𝑤−𝐶_𝑚) is Soret number, 𝐴=𝐴𝑟_𝑒(𝐶_𝑤−𝐶_𝑚)^𝑚−1𝑑²𝜌/𝜇 is the reaction rate parameter, Df=𝐷_𝑚𝑘_𝑇(𝐶_𝑤−𝐶_𝑚)/2𝜈𝐶_𝑠𝐶_𝑝(𝑇_𝑤−𝑇_𝑚) is Dufour number, and 𝛼=1/𝜇𝑎₁𝑎₂ is Eyring Powell parameter.

Thus, the boundary conditions (9) in their dimensionless form read 𝑢=0,𝑤=0,𝑇=𝐶=0at𝑟=𝑟₁=12,𝑢=−𝜕ℎ𝜕𝑧,𝑤=−1,𝑇=𝐶=1at𝑟=ℎ=1+𝜀sin2𝜋𝑧. (12)

After applying the long-wavelength approximation to our analysis as described in [5], then the field (11) now gives 𝜕𝑢+𝑢𝜕𝑟𝑟+𝜕𝑤𝜕𝑧=0,𝜕𝑃𝜕𝜕𝑟=0,₂𝑤𝜕𝑟²+1𝑟𝜕𝑤×𝛼𝜕𝑟1+1/𝑎₁²(𝜕𝑤/𝜕𝑟)²−1+1Da+𝑀𝑤=𝜕𝑃,1𝜕𝑧𝜕Pr²𝑇𝜕𝑟²+1𝑟𝜕𝑇𝜕𝑟+Ec𝜕𝑤𝜕𝑟²+Ec𝛾𝜕𝑤𝜕𝑟×sinh⁻¹√𝛾𝜕𝑤𝜕𝜕𝑟+Df²𝐶𝜕𝑟²+1𝑟𝜕𝐶𝜕𝜕𝑟=0,₂𝐶𝜕𝑟²+1𝑟𝜕𝐶𝜕𝜕𝑟+Sc⋅Sr²𝑇𝜕𝑟²+1𝑟𝜕𝑇𝜕𝑟−𝐴⋅Sc⋅𝐶^𝑚=0. (13)(14)(15)(16)(17) Equation (14) indicates that 𝑃 is not function of 𝑟. Hence 𝑃 is function of 𝑧 and 𝑡 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 𝑤_𝑖+1−2𝑤_𝑖+𝑤_𝑖−1ℎ²+1𝑤𝑖ℎ^𝑖+1−𝑤_𝑖−1×𝛼2ℎ1+1/𝑎²𝑤_𝑖+1−𝑤_𝑖−1/2ℎ²−1+1𝑤Da+𝑀_𝑖=𝜕𝑃,1𝜕𝑧𝑇Pr^𝑖+1−2𝑇_𝑖+𝑇_𝑖−1ℎ²+1𝑇𝑖ℎ^𝑖+1−𝑇_𝑖−1𝑤2ℎ+Ec^𝑖+1−𝑤_𝑖−12ℎ²+Ec𝛾𝑤_𝑖+1−𝑤_𝑖−12ℎ×sinh⁻¹√𝛾𝑤_𝑖+1−𝑤_𝑖−1𝐶2ℎ+Df^𝑖+1−2𝐶_𝑖+𝐶_𝑖−1ℎ²+1𝐶𝑖ℎ^𝑖+1−𝐶_𝑖−1𝐶2ℎ=0,_𝑖+1−2𝐶_𝑖+𝐶_𝑖−1ℎ²+1𝐶𝑖ℎ^𝑖+1−𝐶_𝑖−1𝑇2ℎ+Sc⋅Sr^𝑖+1−2𝑇_𝑖+𝑇_𝑖−1ℎ²+1𝑇𝑖ℎ^𝑖+1−𝑇_𝑖−1𝐶2ℎ−𝐴⋅Sc⋅_𝑖^𝑚=0. (18) The boundary conditions (12) transform into 𝑤_𝑖=0,𝑇_𝑖=𝐶_𝑖=0at𝑟=𝑟₁=12,𝑤_𝑖=−1,𝑇_𝑖=𝐶_𝑖=1at𝑟=ℎ=1+𝜀sin2𝜋𝑧. (19)

4. Results and Discussion

To discuss the effects of various parameters involved in the problem such as Prandtl number Pr, magnetic field parameter 𝑀, Eckert number Ec, Darcy number Da, Dufour number Df, reaction rate parameter 𝐴, Soret number Sr, Schmidt number Sc, Eyring Powell parameter 𝛼, and the pressure gradient 𝑑𝑃/𝑑𝑧, 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 𝑤, the temperature 𝑇, and the concentration 𝐶 are illustrated in Figures 2–25.

Figure 2 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, 𝑀=2, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 3 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 4 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 5 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 6 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 7 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, Df=0.07, and 𝑚=2.

Figure 8 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, 𝐴=0.1, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, Df=0.07, and 𝑚=2.

Figure 9 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, 𝐴=0.1, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sr=2, Df=0.07, and 𝑚=2.

Figure 10 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 11 

Profiles of the dimensionless temperature 𝑇(𝑟) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 12 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, 𝑀=2, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 13 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 14 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 15 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 16 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 17 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, Sr=2, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, Df=0.07, and 𝑚=2.

Figure 18 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, 𝐴=0.1, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sc=0.5, Df=0.07, and 𝑚=2.

Figure 19 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Ec=0.3, 𝐴=0.1, Da=0.5, 𝛼=0.3, 𝑎=0.4, Sr=2, Df=0.07, and 𝑚=2.

Figure 20 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 21 

Profiles of the dimensionless concentration 𝐶(𝑟) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 22 

Profiles of the dimensionless velocity 𝑤(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 23 

Profiles of the dimensionless velocity 𝑤(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, 𝑀=2, Sr=2, Df=0.07, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 24 

Profiles of the dimensionless velocity 𝑤(𝑟) for a system have the particulars 𝑑𝑃/𝑑𝑧=−30, Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

Figure 25 

Profiles of the dimensionless velocity 𝑤(𝑟) for a system have the particulars Pr=3, Da=0.5, Sr=2, Df=0.07, 𝑀=2, 𝛼=0.3, 𝑎=0.4, Ec=0.3, Sc=0.5, 𝐴=0.1, and 𝑚=2.

The values of 𝑇 are plotted versus the radial coordinate in Figures 2–11. In Figure 2, the effect of Pr on the temperature distribution 𝑇 is indicated, and this figure shows that the temperature increases with the increasing of Pr. Also the relation between 𝑇 and 𝑟 is approximately linear for small values of Pr (<0.7), but for large values of Pr (>0.7), 𝑇 increases with increasing 𝑟 till a finite value 𝑟 (maximum value), after which 𝑇 decreases with increasing 𝑟. The effect of 𝑀 on the temperature distribution 𝑇 is shown in Figure 3; it is clear from this figure that the temperature decreases with increasing 𝑀. It is also noted that the temperature increases with increasing 𝑟 and reaches a maximum value (at a finite value of 𝑟: 𝑟=0.85) after which it decreases. Figure 4 gives the effect of Ec on the temperature distribution 𝑇; from this figure, we observe that 𝑇 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 𝑇, which were taken at Ec = 0.2, 0.3, 0.4, and 0.5 when the other parameters are constant, that is, Ec affects 𝑇 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 𝑇 disappears at the range 0.96<𝑟<1. Also, it is clear that the relation between 𝑇 and 𝑟 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 𝑇 occurs at 0.85<𝑟<0.9, and it increases as 𝑟 increases. In Figure 7, the relation between 𝑇 and 𝑟 for various values of 𝐴 is obtained. From this figure, it is clear that the temperature increases with increasing 𝐴, and, from this figure, we can see that although there is a big difference between the values of 𝐴, there is a very small difference between the four curves of 𝑇, which were taken at 𝐴 = 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 𝑇 disappears near the boundaries of the tubes and all curves have a maximum value of 𝑇 in a finite range of 𝑟 namely 0.88<𝑟<0.9. Figures 9 and 10 show that the effect of both Sc and 𝛼 on the temperature distribution 𝑇 is similar to the effect of Sr on 𝑇 given in Figure 8. The effect of 𝑑𝑃/𝑑𝑧 on the temperature distribution is indicated in Figure 11. It is seen from this figure that the temperature decreases with increasing 𝑑𝑃/𝑑𝑧, and all carves have a maximum value as in Figure 8.

Figures 12–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 𝐶 and 𝑟, and this relation is approximately linear at the values 0.7 and 1 of Pr, but, for large values of Pr (≥3), 𝐶 decreases with increasing 𝑟 till a finite value of 𝑟 (minimum values), after which 𝐶 increases with increasing 𝑟. Figure 13 shows the effect of 𝑀 on the concentration distribution, it is observed that the concentration increases with increasing 𝑀 and all curves have a minimum value at the range 0.98<𝑟<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 𝐶. Figure 15 describes the relation between 𝐶 and 𝑟 for various values of Da, and it is clear that the concentration decreases with increasing Da but the effect of Da on 𝐶 disappears at the range 0.96<𝑟<1. Also, it is clear that the relation between 𝐶 and 𝑟 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 𝐶 decreases with increasing 𝑟 in the range from 0 to 0.7 but after which it increases with increasing 𝑟. Figure 17 describes the relation between 𝐶 and 𝑟 for various values of 𝐴, it is obvious that the effect of 𝐴 on 𝐶 is similar to the effect of Df on 𝐶, 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 𝐶 decreases with increasing Sr, and, from this figure, we can see that the effect of Sr on 𝐶 disappears in the range 0.98<𝑟<1. Also, the relation between 𝐶 and 𝑟 is approximately linear at the value namely, 1.5 of Sr. All the curves of 𝐶 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 𝐶 decreases with increasing Sc; also Sc affects the relation between 𝐶 and 𝑟. This relation is approximately linear for small values of Sc, but, for large values of Sc, 𝐶 decreases with increasing 𝑟 till a finite value of 𝑟 (minimum value), after which it increases. Figure 20 describes the relation between 𝐶 and 𝑟 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 𝐶 disappears at the range 0.5<𝑟<0.55 and 0.96<𝑟<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 𝑑𝑃/𝑑𝑧 on the concentration distribution. The calculations were obtained at 𝑑𝑃/𝑑𝑧 = −35, −30, −25, and −20 when the other parameters are constant. 𝐶 increases with increasing 𝑑𝑃/𝑑𝑧. It is also noted from all curves in this figure that the concentration decreases by increasing 𝑟 at the range 0.5<𝑟<0.65 but after that it increases with increasing 𝑟.

Figures 22–25 show the effect of the various parameters on the velocity distribution. Figure 22 gives the effect of 𝑀 on the velocity distribution. This figure shows that the velocity decreases by increasing 𝑀 at the range 0.5<𝑟<0.84 and after that the velocity increases by increasing 𝑀. Also, it is clear that, in all curves in this figure, the maximum value of 𝑤 occurs at 0.65<𝑟<0.7. Figure 23 illustrates the effect of Da on the velocity distribution. It is clear that 𝑤 increases by increasing Da in the range 0.5<𝑟<0.84 and 𝑤 decreases by increasing Da in the range 0.84<𝑟<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 𝑤 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 𝑤 and 𝑟 for various values of 𝛼. From this figure, it is clear that the velocity decreases with increasing 𝛼, but the effect of 𝛼 on 𝑤 disappears in the range 0.96<𝑟<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 𝑑𝑃/𝑑𝑧 on the velocity distribution. The calculation was obtained at 𝑑𝑃/𝑑𝑧 = −35, −30, −25, and −20, and the other parameters are constant. The velocity decreases with increasing 𝑑𝑃/𝑑𝑧 and in all curves the maximum value of 𝑤 occurs in the range 0.65<𝑟<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 velocity 𝑤 rise decreases with the increase of each of 𝑀, 𝛼 and 𝑑𝑃/𝑑𝑧, whereas it increases as Da increases. But, all the other parameters such as, Pr, Ec, Df, 𝐴, Sr, and Sc do not affect on the velocity (the flow of fluid).(2)The temperature 𝑇 rise increases with the increase of each of Pr, Ec, Da, Df, 𝐴, Sr, Sc, and 𝛼, whereas it decreases as 𝑀 and 𝑑𝑃/𝑑𝑧 increase.(3)The concentration 𝐶 rise increases with the increase each of 𝑀 and 𝑑𝑃/𝑑𝑧, whereas it decreases as Pr, Ec, Da, Df, 𝐴, 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.