## Abstract

The present study deals with simultaneous effects of Joule heating and slip on peristaltic flow of a Bingham fluid in an inclined tapered porous channel with elastic walls. The closed form solutions for the stream function, the velocity and the temperature fields are obtained. The effects of the physical parameters on the flow characteristics are presented through graphs for both slip and no-slip cases. In addition, the performance of the temperature is studied with and without Joule heating effects. Moreover, the trapping phenomenon is analysed. The size of the trapped bolus increases with increasing values of the slip parameter and decreasing values of the magnetic, the permeability and the yield stress parameters. The present results are compared with the available results in the literature and our results agree well with the available results for some special cases.

## 1 Introduction

Peristaltic transport is a mechanism of a fluid flow produced by propagation of wave trains along the channel walls. This phenomenon has wide range of practical applications in physiology and biomedical engineering such as swallowing of foodstuff, blood movement in blood vessels, lymph drive in lymphatic vessels, urine transport through ureter, chyme movement in intestinal tract, ovum transport, bile flow in bile duct, etc. Initially Latham [1] and Shapiro et al. [2] investigated the mechanism of peristalsis. Later many investigators have tried to understand the behaviour of peristalsis by considering different fluids and geometries. Elshehawey and Mekhemier [3] studied the peristaltic flow of couple-stress fluids. Usha and Ramachandra Rao [4] investigated the impact of the peristalsis on two-layered power-law fluid flows. Misra and Pandey [5] discussed the peristaltic flow of a non-Newtonian fluid by considering a peripheral layer.

The study of non-Newtonian fluid flow problems with heat transfer has many applications in chemical engineering and related industries; also, due to its tremendous applications in biomedical engineering, researchers have paid considerable attention on peristaltic flow of convective non-Newtonian fluids. The effects of heat transfer on peristaltic transport of a Jeffrey fluid in a vertical channel with porous medium was analyzed by Vajravelu et al. [7]. Tripathi [8] presented a mathematical model for explaining the impact of heat transfer on swallowing of food bolus through the esophagus. Vajravelu et al. [9] made a theoretical model to study the peristaltic flow of a MHD phan-thien-tanner fluid in an asymmetric channel with heat transfer. Akram et al. [10] and Nadeem et al. [11] studied the impact of peristaltic flow of non-Newtonian fluids. Very recently the authors [12, 13, 14, 15, 16, 17, 18] have studied the peristalsis by considering different fluids and geometries.

Not much work has been reported on peristaltic flow by considering simultaneous effects of slip and magnetic field with heat transfer. Furthermore, the consideration of wall properties is essential to understand the behavior of physiological fluids. Srinivas et al. [27] examined the impact of slip conditions and heat transfer on MHD peristaltic transport. The effects of slip and wall properties on the peristaltic transport of a MHD Bingham fluid with heat transfer, was presented by Lakshminarayana et al. [28]. Satyanarayana et al. [29] presented a model to explain the effects of magnetohydrodynamics and heat transfer on peristaltic slip flow of Bingham fluid in porous channel with flexible walls. Moreover, the impact of wall properties is discussed by Srinivas and Kothandapani [30], Hayat et al. [31], Riaz et al. [32] and Sucharitha et al. [33, 34].

Present paper describes a mathematical model to investigate the impact of Joule heating and slip on magnetohydrodynamic peristaltic flow of a Bingham fluid in an inclined non-uniform porous channel with flexible walls. The expressions for the stream function, the velocity and the temperature fields are obtained. The effects of the physical parameters on the flow quantities are discussed in detail. The present study reveals many interesting results which could facilitate the further investigation in convective non-Newtonian fluid flow phenomenon.

## 2 Mathematical formulation

Consider the two dimensional flow of a MHD Bingham fluid in an inclined non-uniform channel with porous medium. Flow is due to sinusoidal wave trains propagating along the elastic walls of the channel with a constant speed *c* (see Fig. 1). The channel wall deformation is assumed as

where *d*(*x*) = *d*+* mx*,

*m*

Using long wavelength and small Reynolds number assumptions (see [28, 30, 34, 35] for details), the simplified non-dimensional governing equations and corresponding boundary conditions for the present study can be written as

Momentum and energy equations are

where *N*^{2} = *M*^{2}+*σ*^{2}

The corresponding boundary conditions (see Ref. [28, 30]) are

(Flexible boundary condition)

(Slip condition)

The non-dimensional parameters and quantities used in the above governing equations are

where *u*, *v*, *p* and *ψ* are the velocity components, pressure and stream functions respectively, *ρ* is the density, *μ* is the viscosity, *d* is the mean width of the channel, *a* is the amplitude, *λ* is the wavelength, *c* is the wave speed, *ξ* is the specific heat, *ν* is the kinematic viscosity, *k*_{0} is the thermal conductivity, *k* is the permeability, *σ*_{0} is the electrical conductivity, *B*_{0} is the magnetic field, *g* is the acceleration due to gravity, *T* is the temperature, *τ*_{0} is the yield stress, *Ec* is the Eckert number, *E*_{1}, *E*_{2} and *E*_{3} are the elasticity parameters, *m* is the non-uniform parameter, *σ* is the permeability parameter, *α* is the inclination angle, *β* is the slip parameter, *Br* is the Brinkman number, *M* is the magnetic parameter, *δ* is the wave number,*ε* is the amplitude ratio, *Pr* is the Prandtl number and *Re* is the Reynolds number.

## 3 Solution of the problem

By differentiating Eq. (2.2) with respect to y we obtain

By solving Eq. (3.1) with boundary conditions (2.5), (2.6) and (2.7) we obtain the stream functions and corresponding velocities in the plug and non-plug flow regions as

where

By using the Eqs. (3.4) and (2.8) in Eq. (2.4) we obtain the temperature as

where

The nusselt number at the wall is given by

## 4 Results and discussion

To find the impact of physical parameters, we have plotted the velocity, temperature and heat transfer coefficient profiles in figures 2 – 22 with the fixed values of *x* = 0.2, *t* = 0.1, *ε* = 0.24, *m* = 0.2, *M* = 1, *σ* = 1, *τ*_{0} = 0.4, *η* = 1, *α* = *π*/4, *E*_{1} = 0.2, *E*_{2} = 0.2, *E*_{3} = 0.1, *Br* = 0.2 (for temperature *E*_{1} = 0.005, *E*_{2} = 0.005, *E*_{3} = 0.001). It is observed that the velocity is higher in slip flow when compared with the nonslip flow whereas the temperature exhibits the opposite behavior. In fact, the growth in slip diminishes the friction between the wall and the fluid. This may be the cause to increase in the velocity and reduction in the internal heat production. The impacts of magnetic parameter *M* and permeability parameter *σ* on the velocity and the temperature fields are presented in figures 2 – 5. We found that an increase in *M* and *σ* reduces the velocity profiles. This is due to the influence of drag force (Lorentz force) which opposes the flow. We also noticed an enhancement in the temperature field. This may be due to the Joule heating impact. Similar behavior has been noticed in the study of Hayat et al. [23].

Velocity profiles for *M*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Velocity profiles for *σ*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *M*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *σ*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

From figures 6–9, it is observed that an increase in the yield stress *τ*_{0} is to reduce the velocity field and to increase the temperature field. We noticed that the higher values of the non-uniform parameter *m* increase both the velocity and the temperature profiles. The effects of the wall flexibility parameters are described trough the figures 10–13. We notice that the velocity increases with increasing *E*_{1} and *E*_{2} (due to the wall tension and mass charectrization property); but the velocity drop is noticed for increasing values of *E*_{3} (due to the damping force). The influence of the elasticity parameters on the temperature has similar effects as in the case of the velocity. This agrees with the results of Hayat et al. [21].

Velocity profiles for *τ*_{0}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Velocity profiles for *m*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *τ*_{0}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *m*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Velocity profiles for *E*_{1} and *E*_{2}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Velocity profiles for *E*_{3}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *E*_{1} and *E*_{2}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *E*_{3}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

From figures 14 and 15 we see that for increasing values of the inclination angle *al pha* both the velocity and temperature fields increase. The impact of the Brinkman number *Br* is shown in figure 16. An increase in *Br* enhances the temperature field.

Velocity profiles for *α*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *α*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for *Br*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Figures 17 and 18 are drawn to validate the present results with the published availableresults in the literature. It is perceived that for *α* = 0 and in the absence of Joule heating our results reduce to those of Satyanarayana et al. [19]. Further, we notice an increase in the values of the velocity and temperature fields when compared with the results in [19], due to the presence of Joule heating. The deviations in heat transfer coefficient are shown in figures 19–22. It can be seen that the absolute value of the Nusselt number is higher in the slip flow compared to the nonslip flow. Figures 19 and 20 reveal that the heat transfer coefficient increases for large values of *M* (due to the existence of Joule heating) whereas it is a decreasing function of *σ*. From figures 21 and 22, we observe an increase in the Nusselt number with the increasing *Br* and decreasing *α*.

Velocity profiles for validation.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Temperature profiles for validation.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Variation of Nusselt number for *M*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Variation of Nusselt number for *σ*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Variation of Nusselt number for *Br*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Variation of Nusselt number for *α*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

## 5 Trapping phenomenon

The process of the formation of fluid bolus by the closed streamlines in the fluid flow is called trapping phenomenon and the trapped bolus moves forward with the peristaltic wave. The behaviour of streamlines is presented in figures 23–26. The effect of the slip parameter *β* is shown in figure 23. An increase in *β* increases the size of the bolus. From figures 24–26, we notice that an increase in the magnetic parameter *M*, the permeability parameter *sigma* and the yield stress *τ*_{0} reduces the size of the trapped bolus.

Stream lines for *β*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Stream lines for *M*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Stream lines for *σ*.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

Stream lines for *τ*_{0}.

Citation: Applied Mathematics and Nonlinear Sciences 3, 1; 10.21042/AMNS.2018.1.00005

## 6 Conclusions

The effects of the wall slip and the Joule heating on MHD peristaltic flow of a Bingham fluid in a porous channel are studied. Long wavelength and small Reynolds number approximations are used to obtain the exact solutions of the problem. The results have applications in biomedical engineering and oil industries. Some of the interesting results are summarised as follows:

Velocity is an increasing function of the slip parameter

*β*whereas it is a decreasing function of the magnetic parameter*M*, the yield stress*τ*_{0}and the permeability parameter*σ*.The effect of the Joule heating is to increase the temperature.

Slip parameter

*β*reduces the temperature while it is an increasing function of the magnetic parameter*M*and the permeability parameter*σ*.Inclination angle

*α*and and the non-uniform parameter*m*increase both the velocity and the temperature fields.Increase in both the velocity and the temperature is identified due to increase in the wall parameters

*E*_{1}and*E*_{2}, while reduction is noticed in the case of other wall parameter*E*_{2}.Nusselt number is an increasing function of the slip parameter

*β*and the magnetic parameter*M*.Size of the trapped bolus is reduced for the increasing values of

*M*and*σ*while it increases with increasing*σ*.

Communicated by Juan L.G. Guirao

## References

- [1]↑
T. W. Latham, (1966), Fluid motions in a peristaltic pump, Doctoral dissertation, Massachusetts Institute of Technology.

- [2]↑
A. H. Shapiro, M. Y. Jaffrin and S. L. Weinberg, (1969), Peristaltic pumping with long wavelengths at low Reynolds number, Journal of Fluid Mechanics 37, No 4, 799-825.

- [3]↑
E. F. El Shehawey and K. S. Mekheimer, (1994), Couple-stresses in peristaltic transport of fluids, Journal of Physics D: Applied Physics 27, No 6, 1163-1170.

- [4]↑
S. Usha and A. Ramachandra Rao, (1997), Peristaltic Transport of Two-Layered Power-Law Fluids, Journal of Biomechanical Engineering 119, No 4, 483-488.

- [5]↑
J.C. Misra and S.K. Pandey, (1999), Peristaltic transport of a non-Newtonian fluid with a peripheral layer, International Journal of Engineering Science 37, No 14, 1841-1858.

- [6]
S. Nadeem and Noreen Sher Akbar, (2009), Influence of heat transfer on a peristaltic transport of Herschel-Bulkley fluid in a non-uniform inclined tube, Communications in Nonlinear Science and Numerical Simulation 14, No 12, 4100-4113.

- [7]↑
K. Vajravelu, S. Sreenadh and P. Lakshminarayana, (2011), The influence of heat transfer on peristaltic transport of a Jeffrey fluid in a vertical porous stratum, Communications in Nonlinear Science and Numerical Simulation 16, No 8, 3107-3125.

- [8]↑
D. Tripathi, (2012), A mathematical model for swallowing of food bolus through the oesophagus under the influence of heat transfer, International Journal of Thermal Sciences 51, 91-101.

- [9]↑
K. Vajravelu, S. Sreenadh, S. Dhananjaya and P. Lakshminarayana, (2016), Peristaltic Flow and Heat Transfer of a Conducting Phan-Thien-Tanner Fluid in an Asymmetric Channel – Application to Chyme Movement in Small Intestine, International Journal of Applied Mechanics and Engineering 21, No 3, 713-736.

- [10]↑
S. Akram, S. Nadeem and M. Hanif, (2013), Numerical and analytical treatment on peristaltic flow of Williamson fluid in the occurrence of induced magnetic field, Journal of Magnetism and Magnetic Materials 346, 142-151.

- [11]↑
S. Nadeem, Arshad Riaz and R. Ellahi, (2013), Peristaltic flow of a Jeffrey fluid in a rectangular duct having compliant walls, Chemical Industry and Chemical Engineering Quarterly 19, 399-409.

- [12]↑
I. Shahzadi, S. Nadeem and F. Rabiei, (2017), Simultaneous effects of single wall carbon nanotube and effective variable viscosity for peristaltic flow through annulus having permeable walls, Results in Physics 7, 667-676.

- [13]↑
I. Shahzadi, H. Sadaf, S. Nadeem and A. Saleem, (2017), Bio-mathematical analysis for the peristaltic flow of single wall carbon nanotubes under the impact of variable viscosity and wall properties, Computer Methods and Programs in Biomedicine 139, 137-147.

- [14]↑
I. Shahzadi and S. Nadeem, (2017), Role of Inclined Magnetic Field and Copper Nanoparticles on Peristaltic Flow of Nanofluid through Inclined Annulus: Application of the Clot Model, Communications in Theoretical Physics 67, No 6, 704-716.

- [15]↑
S. Akram and S. Nadeem, (2017), Influence of Nanoparticles Phenomena on the Peristaltic Flow of Pseudoplastic Fluid in an Inclined Asymmetric Channel with Different Wave Forms, Iranian Journal of Chemistry and Chemical Engineering 36, No 2, 107-124.

- [16]↑
H. Sadaf and S. Nadeem, (2017), Analysis of Combined Convective and Viscous Dissipation Effects for Peristaltic Flow of Rabinowitsch Fluid Model, Journal of Bionic Engineering 14, No 1, 182-190.

- [17]↑
G. Sucharitha, S. Sreenadh, P. Lakshminarayana and K. Sushma, (2017), Brownian motion and thermophoresis effects on Peristaltic slip flow of a MHD nanofluid in a symmetric/asymmetric channel, IOP Conference Series: Materials Science and Engineering 263 062025.

- [18]↑
G. Sucharitha, P. Lakshminarayana and N. Sandeep, (2017), Dual solutions of cross diffusion effects on MHD Peristaltic flow in a conduit, IOP Conference Series: Materials Science and Engineering 263 062024.

- [19]↑
Y.S. Wu, K. Pruess, P.A. Witherspoon, (1990), Flow and Displacement of Bingham Non-Newtonian Fluids in Porous Media, Lawrence Berkeley National Laboratory. LBNL Report #: LBL-28340. Retrieved from https://escholarship.org/uc/item/7q14212w

- [20]
V. P. Rathod and D. Laxmi, (2014), Effects of heat transfer on the peristaltic MHD flow of a Bingham fluid through a porous medium in a channel, International Journal of Biomathematics 7, No 6, 1450060 (20 pages).

- [21]↑
J. C. Misra, S. Maiti and G. C. Shit, (2008), Peristaltic transport of a physiological fluid in an asymmetric porous channel in the presence of an external magnetic field, Journal of Mechanics in Medicine and Biology 8, No 4, 507-525.

- [22]
S. K. Pandey and M. K. Chaube, (2011), Effect of magnetic field on peristaltic transport of couple stress fluids through a porous medium, Journal of Biological Systems 19, No 2, 251-262.

- [23]↑
O. Anwar Bég, M. M. Rashidi, T. A. Bég and M. Asadi, (2012), Homotopy analysis of transient mgneto-bio-fluid dynamics of micropolar squeeze film in a porous medium: A model for mgneto-bio-rheological lubrication, Journal of Mechanics in Medicine and Biology 12, No 3, 1250051 (21 pages).

- [24]
R. Ellahi, M. M. Bhatti, A. Riaz and M. Sheikholeslami, (2014), Effects of magnetohydrodynamics on peristaltic flow of Jeffrey fluid in a rectangular duct through a porous medium, Journal of Porous Media 17, No 2, 143-157.

- [25]
K. Chakradhar, K. Nandagopal, S. Sreenadh, P. Lakshminarayana and G. Sucharitha, (2015), Peristaltic pumping of a micropolar fluid in a tube with permeable wall, International Journal of Mathematical Archive 6, No 3, 173-183.

- [26]
K. Vajravelu, S. Sreenadh, P. Lakshminarayana and G. Sucharitha, (2016), The effect of heat transfer on the nonlinear peristaltic transport of a Jeffrey fluid through a finite vertical porous channel, International Journal of Biomathematics 9, No 2, 1650023 (2016) (24 pages).

- [27]↑
S. Srinivas, R. Gayathri and M. Kothandapani, (2009), The influence of slip conditions, wall properties and heat transfer on MHD peristaltic transport, Computer Physics Communications 180, No 11, 2115-2122.

- [28]↑
P. Lakshminarayana, S. Sreenadh and G. Sucharitha, (2015), The Influence of Slip, Wall Properties on the Peristaltic Transport of a Conducting Bingham Fluid with Heat Transfer, Procedia Engineering 127, 1087-1094.

- [29]↑
K. V. V. Satyanarayana, S. Sreenadh, G. Sucharitha and P. Lakshminarayana, (2016), The Effect of Wall Properties on the Convective Peristaltic Transport of a Conducting Bingham Fluid through Porous Medium, Indian Journal of Science & Technology 9, No 42, 1-9.

- [30]↑
S. Srinivas and M. Kothandapani, (2009), The influence of heat and mass transfer on MHD peristaltic flow through a porous space with compliant walls, Applied Mathematics and Computation 213, No 1, 197-208.

- [31]↑
T. Hayat, Z. Nisar, H. Yasmin and A. Alsaedi, (2016), Peristaltic transport of nanofluid in a compliant wall channel with convective conditions and thermal radiation, Journal of Molecular Liquids 220, 448-453.

- [32]↑
A. Riaz, S. Nadeem, R. Ellahi and N. S. Akbarc, (2014), The influence of wall flexibility on unsteady peristaltic flow of Prandtl fluid in a three dimensional rectangular duct, Applied Mathematics and Computation 241, 389-400.

- [33]↑
G. Sucharitha, K. Vajravelu, S. Sreenadh and P. Lakshminarayana, (2017), Peristaltic flow and heat transfer of a Herschel-Bulkley fluid in an inclined non-uniform channel with wall properties, IOP Conference Series: Materials Science and Engineering 263 062026.

- [34]↑
G. Sucharitha, P. Lakshminarayana and N. Sandeep, (2017), Joule heating and wall flexibility effects on the peristaltic flow of magnetohydrodynamic nanofluid, International Journal of Mechanical Sciences 131-132, 52-62.

- [35]↑
T. Hayat, F.M. Abbasi, M. Al-Yami and S. Monaquel, (2014), Slip and Joule heating effects in mixed convection peristaltic transport of nanofluid with Soret and Dufour effects, Journal of Molecular Liquids 194, 93-99.