## Abstract

In this paper, the peristaltic wave propagation of a Non-Newtonian Casson liquid in a non-uniform (flexible)channel with wall properties and heat transfer is analyzed. Long wavelength and low Reynolds number approximations are considered. Analytical solution for velocity, stream function and temperature in terms of various physical parameters is obtained. The impact of yield stress, elasticity, slip and non-uniformity parameters on the peristaltic flow of Casson liquidare observed through graphs and discussed. The important outcome is that an increase in rigidity, stiffness and viscous damping force of the wall results in the enhancement of the size and number of bolus formed in the flow pattern.

## 1 Introduction

Peristaltic wave proliferation is a system of transporting liquid from lower pressure to higher pressure. This peristaltic wave spread discovers its applications in biology, medical and engineering field. Further, in view of this guideline, modern peristaltic pumps are likewise outlined. Numerous examinations on peristaltic stream wonder have been performed in tube and channels. A large portion of the examinations in industry and science demonstrate that the liquid conduct is non-Newtonian. Subsequently a few specialists are focusing on the stream of non-Newtonian liquids through peristalsis in tubes and channels (Scott Blair [1], Vajravelu et al. [2]- [3]. A large portion of the examinations in the literature has influenced a way to deal with concentrate to and comprehend the urine move through ureter and flow of blood in large arteries where shear rate is high, yet it neglects to clarify the complex rheological conduct of blood in thin veins where the shear rates are low. The examination on non-Newtonian nature of bloodstream has been of most significance to scientists lately because of their application in exploring the conduct of blood in narrow arteries. Casson liquid is one such non-Newtonian liquid which shows yield pressure and fits the streaming blood when the shear rates are low (See Casson [4]). Ref. [1] watched that at low shear rates, Casson model was more precise in anticipating the physiological practices of blood. Srivastava and Srivastava [5] assumed blood as immiscible fluids and studied on the peristaltic pumping of blood by considering Casson fluid. Recently, Vajravelu et al. [6] investigated the flow of Casson liquid under peristalsis and elasticity. Numerous researchers did the investigation of Casson model under various geometric conditions as of late (Nagarani [7]; Prasad et al. [8]- [9]; Vajravelu et al. [10]; Prasad et al. [11]). In the event that the liquid is moving with the impact of peristalsis, it is extremely intriguing reality to think about the elastic nature of the channel. Numerous examinations have been conveyed with a stream of non-Newtonian liquids in elastic tubes (Vajravelu et al. [12], Nadeem and Ijaz [13], Siva et. al. [14], Badari et. al. [15]- [16], Rajashekhar et al. [17]).

Heat transfer in a biological framework is a natural phenomenon. Further, bio-heat transfer show is considered keeping in mind the end goal to outline the impacts of blood perfusion and metabolic heat age in living tissues. Heat is a type of energy that is transferred across a boundary because of the temperature distinction. The standards of variety in temperature in designing structures can be associated with the human body to choose how the body trades heat. Heat is created in the body by the constant absorption of nutrient supplements which offers vitality to the systems of the body. Blood moving through the vessels goes about as a convective liquid and keeps in charge any advancement of heat inside the tissues of the body. The heat conveyed by the blood is administered by the temperature of the neighbouring tissues, the measurement of the veins, the thickness of the liquid, liquid speed and the heat trade coefficient of the blood. The investigation of heat transfer impacts alongside slip conditions on peristalsis has procured a tremendous measure of enthusiasm among the scientists in the previous decades. The examination of heat transfer has discovered its application in the field of biofluid mechanics, substance designing and prescription(See Prasad et al. [18], [19], [20]). Numerous analysts analysed the communication amongst peristalsis and heat transfer in various geometries along with/without slip conditions.Hayat et al [21], [22] employed analytical method to examine the Influence of slip and heat transferon the peristaltic transport in a channel.Radhakrishnamacharya et al. [23] observed the flow of Newtonian liquid in a channel with wall effects and heat transfer.Hayat et al. [24] continued the work of Ref. [23] by considering Power law fluid. Further, Lakshminarayana et al. [25] examined the heat transfer analysis on the MHD peristaltic flow of a Bingham fluid in a channel with wall properties.Recently, Nabil et al. [26] concentrated on the peristaltic motion of couple stress liquid in a porous channel with heat transfer. In a flexible channel the peristaltic motion of dusty fluid with MHD and heat transfer was investigated by Hayat andJaved [27].

In perspective of this, the analysis concentrates on the impact of elasticity, boundary slip and heat transfer on the flow of Casson liquid in a channel. An expression for velocity, stream function, and temperature has been calculated analytically under the assumptions of long wavelength and small Reynolds number approximations. The liquid flow depends on many physical expressions such as wall properties, slip parameter, non-uniformity parameter and yield stress. These effects of parameters are discussed in detail through graphs by using MATLAB. Since Casson model closely describes blood flow in physiological systems, the results obtained have important applications in the cardiovascular system.

## 2 Mathematical Formulation

Consider a peristaltic flow of a Casson liquid in a channel with heat transfer and wall effects (see Fig. 1) on which sinusoidal waves of moderate amplitude are imposed. The walls are taken like stretched membranes. The geometry of the channel wall is given by

where *D*(*x*) = *d*+*ω*′*X*, *ω*′ ≪ 1, *a* is the amplitude, *λ* is the wave length, *d* is the mean half width of the channel, *ω*′is the dimensional non-uniformity of the channel.

The equations governing the motion for the present problem are

Where Ω and *V* are the components of velocity along *x* and y directions respectively, *ρ* is the density, *μ*is the coefficient of viscosity of the liquid, *p* is the pressure,*d* is the mean half width of the channel, *a* is the amplitude, *λ* is the wave length, *c* is the phase speed of the wave, *m*′ is the dimensional non-uniformity of the channel, *ξ* is the specific heat at constant volume, *υ* is the kinematic viscosity of the liquid, *k* is the thermal conductivity of the liquid, *Θ* is the temperature of the liquid.

The governing equations of motion of the flexible wall may be expressed as

where Γ^{∗} is an operator, which is used to represent the motion of stretched membrane with viscosity damping forces such that

Here *α* is the elastic tension in the membrane,*β* is the mass per unit area, *γ* is the coefficient of viscous damping forces, *p*_{0} is the pressure on the outside surface of the wall due to the tension in the muscles. Continuity of stress at *y* = *η* and using *x* – momentum equation yields

where
*p*_{0} = 0 and introduce the non-dimensional stream function and non-dimensional quantities as

After dropping primes, we obtain non-dimensional governing equations

Non-dimensional boundary conditions are

where
*ζ* is the Knudsen number (slip parameter).

## 3 Solution of the Problem

Using the long wavelength and low Reynolds number approximations, one can find from equations (12) to (16) that

Equation (19) shows that *p* is not a function of *y*

On differentiating Eq. (18) with respect to *y*, we get

From Eq. (16) we get

The closed form solution for equation (21) using the boundary conditions (15), (17) and (22) can be obtained as

We find the upper limit of plug flow region using the boundary condition that Ψ_{yy} = 0 at*y* = *y*_{0}. It is given by

Taking *y* = *y*_{0} in equation (23) and using the relation (24), we get the velocity in the plug flow region as

By using Equations (23) and (25), we get

where

By solving equation (20) with the help of equation (26) and (17), an expression for temperature field is obtained

## 4 Results and Discussions

The velocity, stream function and temperature are calculated analytically via MATLAB under the assumptions of long wave length and low Reynolds number. The effect of physical parameters like elastic parameters *A*_{1}, *A*_{2} and *A*_{3}, non-uniform parameter *ω*, yield stress *τ*_{0} and Brinkman number *Br* on velocity, temperature and stream function are discussed graphically from Fig 2 to 16.

**Effect of** *A*_{1} **on the velocity distribution for fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω* = 0.1 *τ*_{0} = 0.1,*A*_{2} = 0.3,*A*_{3} = 0.5,*ζ* = 0.1

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

### 4.1 Velocity Field and Temperature distribution

Figures 2 to 4 shows the profiles for velocity Ω versus *y* with different elastic wall parameters namely *A*_{1}, *A*_{2} and *A*_{3}. An increase in the rigidity parameter *A*_{1} gives rise to an increment in velocity Ω. A similar trend may be observed in the case of elastic parameter *A*_{2} whereas the parameter *A*_{3} exhibits opposite behavior.From this it is clear that the fluid flow is more if the tension in the membrane and the mass per unit area of the elastic wall is more and fluid flow will be less if the viscous damping forces in the fluid are more. The change in the velocity Ω for different values of non-uniformity parameter *ω* is elucidated in the Fig. 5. Physically, for the convergent channel (*ω* < 0), uniform channel (*ω* = 0) and divergent channel (*ω* > 0). It is noted from these profiles that the velocity Ω decreases with increase in non-uniformity parameter and in order to obtain a better flow of liquid the channel must be convergent. Fig. 6 demonstrates that the velocity is a decreasing function of yield stress *τ*_{0} .The effect of slip parameter on velocity is observed in Fig. 7. An increase in slip parameter results in the enhancement of the velocity of the liquid.

**Effect of** *A*_{2} **on the velocity distribution for Fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω* = 0.1 *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{3} = 0.5, *ζ* = 0.1

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

**Effect of** *A*_{3} **on the velocity distribution for Fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω* = 0.1 *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *ζ* = 0.1

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

**Effect of** *ω* **on the velocity distribution for Fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ζ* = 0.1 *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

**Velocity Distribution for different** *τ*_{0} **for fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω* = 0.1 *ζ* = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

**Velocity Distribution for different** *ζ* **for fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω* = 0.1 *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

The effect of various parameters on temperature is illustrated in Figures 8 to 11. From Fig. 8 it can be noticed that an increase in *Br* increments in the temperature field. Further, it is noted from the Fig. 9 that temperature increase for large values of *ω*. Fig. 10 depicts that an increase in the value of the yield stress decreases the magnitude of temperature. This behaviour is expected due to the presence of *τ* (minimum amount of energy required to begin the flow) in the Casson model. We see from Fig.11 that as increasing values of *A*_{1}and *A*_{2} increases the temperature. Physically, the tension in the membrane and mass per unit area is more, then we have high temperature.

**Effect of** *Br* **on the Temperature for fixed values of** *x* = 0.2,*t* = 0.l, *ε* = 0.1, *ω* = 0.1, *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

**Effect of** *ω* **on the temperature for fixed values of** *x* = 0.2,*t* = 0.1,*ε* = 0.1, *Br* = 1, *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

**Temperature profiles for different** *τ*_{0} **for fixed values of** *x* = 0.2,*t* = 0.1, *ε* = 0.1, *ω* = 0.1, *Br* = 1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5

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

**Effect of elastic parameters on Temperature profiles for fixed values of** *x* = 0.2,*t* = 0.1,*ε* = 0.1, *ω* = 0.1, *Br* = 1, *τ*_{0} = 0.1

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

### 4.2 Trapping Phenomenon

The most important concept to be noted in peristalsis is trapping, a closed bolus of liquid that moves along with the peristaltic wave. Fig. 12 shows that the size of the bolus increases with increase in rigidity *A*_{1}. From Fig. 13 we observe that as stiffness parameter (*A*_{2}) increases the number of bolus increases. Further increment in *A*_{3} increases the size of the trapped bolus which is noticed in Fig. 14. From Fig. 15 we conclude that the pattern and size of bolus based on the change of non-uniformity parameter. From these figures it is noticed that there is a symmetric behaviour in uniform channel and non-symmetric behaviour in convergent and divergent channel. Also, the pattern of the formation of bolus is opposite in nature for convergent and divergent channels. The effect of yield stress parameter on the trapping is illustrated in Fig. 16. It can be concluded that the size of the trapped bolus decreases with increase in *τ*_{0}. Thus, we have seen the effect of associated parameters *A*_{1}, *A*_{2},*A*_{3},*ω* *and* *τ*_{0} on the progress of the trapping phenomena.These qualitative results may have some significance in understanding the transport of blood in the small blood vessels.

**Effect of** *A*_{1} **on Trapping (I)** *A*_{1} = 0.5 **(II)** *A*_{1} = 0.6 **(III)** *A*_{1} = 0.8 **for fixed values of** *x* = 0.01,*t* = 0.4,*ε* = 0.2, *ω* = 0.1, *τ*_{0} = 0.1,*A*_{2} = 0.3, *A*_{3} = 0.5,*ζ* = 0.1

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

**Effect of** *A*_{2} **on Trapping (I)** *A*_{2} = 0.2 **(II)** *A*_{2} = 0.4 **(III)** *A*_{2} = 0.5 **for fixed values of** *x* = 0.01,*t* = 0.4,*ε* = 0.2, *ω* = 0.1, *τ*_{0} = 0.1,*A*_{1} = 0.1, *A*_{3} = 0.5,*ζ* = 0.1

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

**Effect of** *A*_{3} **on Trapping (I)** *A*_{3} = 0.1 **(II)** *A*_{3} = 0.3 **(III)** *A*_{3} = 0.5 **for fixed values of** *x* = 0.01, *t* = 0.4, *ε* = 0.2, *ω* = 0.1, *τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *ζ* = 0.1

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

**Effect of** *ω***on Trapping (I)** *ω*=–0.1 **(II)** *ω*=0 **(III)** *ω*=0.1 **for fixed values of** *x* = 0.01,*t* = 0.4,*ε* = 0.2,*τ*_{0} = 0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5, *ζ* = 0.1

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

**Effect of** *τ*_{0} **on Trapping (I)** *τ*_{0} = 0.001 **(II)** *τ*_{0} = 0.01 **(III)** *τ*_{0} = 0.1 **for fixed values of** *x* = 0.01,*t* = 0.4, *ε* = 0.2, *ω*=0.1, *A*_{1} = 0.1, *A*_{2} = 0.3, *A*_{3} = 0.5, *ζ* = 0.1

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

## 5 Concluding remarks

The main findings are listed below.

As the elastic effect of the channel increases the velocity of the liquid increases.

The increase in yield stress results out in the decrease in velocity of the liquid flow.

Increase in rigidity

*A*_{1}, stiffness*A*_{2}, viscous damping force*A*_{3}of the wall increases the size and number of bolus formed in the flow pattern.The size of the tapered bolus decreases as the yield stress increases.

If the channel is uniform, then the bolus is symmetric and opposite behavior is observed if the channel is either convergent or divergent.

Communicated by Juan L.G. Guirao

The authors appreciate the constructive comments of the reviewers which led to definite improvement in the paper.

## References

- [1]↑
G. W. Scott Blair,(1959), An equation for the flow of blood, plasma and serum through glass capillaries, Nature, 183, 613-614.

- [2]↑
K. Vajravelu, S. Sreenadh and V. Ramesh Babu, (2005), Peristaltic transport of a Herschel-Bulkley fluid in a channel, Appl. Math and Comput, 169, 726-735.

- [3]↑
K. Vajravelu, S. Sreenadh and V. Ramesh Babu, (2005), Peristaltic transport of a Herschel- Bulkley fluid in an inclined tube, Int. J. Nonlinear Mech. 40, 83-90.

- [4]↑
Casson, N., (1959), A flow equation for pigment oil suspensions of printing ink type Rheology of Dispersed System 81-102, Pergamon press, London.

- [5]↑
L.M. Srivastava, V.P. Srivastava, (1984), Peristaltic transport of blood Casson Model II, J.Biomechanics, 17,821-829.

- [6]↑
K.Vajravelu, S. Sreenadh, P.Devaki, K.V.Prasad, (2016), Peristaltic transport of Casson fluid in an elastic tube, Journal of Applied Fluid Mechanics, 9, No 4, 1897-1905.

- [7]↑
P. Nagarani, (2010), Peristaltic transport of Casson fluid in an inclined channel, Korea-Australia Rheology Journal, 22, 105-111.

- [8]↑
K. V. Prasad, K. Vajravelu, H. Vaidya, I. S. Shivakumara and N. Z. Basha, (2016), Flow and heat transfer of a Casson Nanofluid over a nonlinear stretching sheet, Journal of Nanofluids, 5, 743-752.

- [9]↑
K. V. Prasad, K. Vajravelu, H. Vaidya, (2016), MHD Casson Nanofluid Flow and Heat Transfer at a Stretching Sheet with Variable Thickness, Journal of Nanofluids, 5, No 3, 423-435.

- [10]↑
K. Vajravelu, K. V. Prasad, H. Vaidya, N. Z. Basha and Chiu-On-Ng, (2017), Mixed convective flow of a Casson fluid over a stretching sheet, International Journal of Applied and Computational Mathematics, 3, 1619-1638.

- [11]↑
K. V. Prasad, K. Vajravelu, Hanumesh Vaidya, M. M. Rashidi, and Neelufer, Z.Basha, (2018), Flow and Heat Transfer of a Casson Liquid over a Vertical Stretching Surface: Optimal Solution, American Journal of Heat and Mass Transfer, 5, No 1, 1-22.

- [12]↑
K.Vajravelu, S. Sreenadh, P. Devaki, K.V.Prasad, (2011), Mathematical model for a Herschel-Bulkley fluid flow in an Elastic tube, Central European Journal of Physics, 9, No 5, 1357-1365.

- [13]↑
S. Nadeem and S. Ijaz, (2014), Nanoparticles analysis on the blood flow through a tapered catheterized elastic artery with overlapping stenosis, Eur. Phys. J. Plus, 129: 249.

- [14]↑
P. Siva, A. Govindarajan and M. Vidhya, (2015), Analysis on the effects of wall properties on MHD peristaltic flow of a Dusty fluid through a porous medium, 102, No 2, 247-263.

- [15]↑
C. H. Badari Narayana, P. Devaki and S. Sreenadh, (2016), Effect of Elasticity on Herschel-Bulkley Fluid Flow in a Tube, International Journal of Advanced Research in Innovative Discoveries in Engineering and Applications, 1, No. 2, 9-16.

- [16]↑
C. H. Badari Narayana, P. Devaki and S. Sreenadh, (2017), Effect of Elasticity and inclination on Herschel-Bulkley Fluid Flow in a Tube, International Journal of Advanced Information Science and Technology, 6 No 2, 6-12.

- [17]↑
C. Rajashekhar, G. Manjunatha, K. V. Prasad, B. B. Divya and Hanumesh Vaidya, (2018), Peristaltic transport of two-layered blood flow using Herschel鈦揃ulkley Model, Cogent Engineering, 5, No 1, 1495592.

- [18]↑
K. V. Prasad, K. Vajravelu and Hanumesh Vaidya, (2016), Hall Effect on MHD Flow and Heat Transfer over a Stretching Sheet with Variable Thickness, International Journal for Computational Methods in Engineering Science and Mechanics, 17, No 4, 288-297.

- [19]↑
K. V. Prasad, H. Vaidya, K. Vajravelu and M. M. Rashidi, (2017), Effects of Variable Fluid Properties on MHD Flow and Heat Transfer over a Stretching Sheet with Variable Thickness, 33, No 4, 501-512.

- [20]↑
K. V. Prasad, K. Vajravelu, H. Vaidya and Robert A. Van Gorder, (2017), MHD flow and heat transfer in a nanofluid over a slender elastic sheet with variable thickness, Results in Physics, 7, 1462-1474.

- [21]↑
T. Hayat, S. Hinaand N. Ali, (2010), Simultaneous effects of slip and heat transfer on the peristaltic flow, Communications in Nonlinear Science and Numerical Simulation, 15, 1526-1537.

- [22]↑
T. Hayat, Q. Hussain, U. Qureshi, N. Ali and A. A. Hendi, (2011), Influence of slip condition on the peristaltic transport in an asymmetric channel with heat transfer: an exact solution, International Journal for Numerical Methods in Fluids, 67, 1944-1959.

- [23]↑
G. Radhakrishnamacharya, Ch. Srinivasulu and C. R. Mecanique, (2007), Influence of wall properties on peristaltic transport with heat transfer, ComptesRendusM茅canique, 335, No 7, 369-373.

- [24]↑
T. Hayat, M. Javed, S. Asghar and A. A. Hendi, (2013), Wall properties and heat transfer analysis of the peristaltic motion in a power-law fluid, International journal for Numerical methods in fluids, 71, No. 1, 65-79.

- [25]↑
P. Lakshminarayana, S. Sreenadh and G. Sucharitha, (2015), Influence of Slip, Wall Properties on Peristaltic Transport of a Conducting Bingham fluid with heat transfer, IOP Conf. Series: Materials Science and Engineering, Procedia Engineering, 127, 1087-1094.

- [26]↑
Nabil T. M. Eldabe, Mohamed A. Hassan and Mohamed Y. Abou-Zeid,(2016), Wall properties effect on the peristaltic motion of a couple stress fluid with heat and mass transfer through a porous medium, Journal of engineering mechanics, 142, No 3.

- [27]↑
T. Hayat and M. Javed, (2017), Effects of heat transfer on MHD peristaltic transport of dusty fluid in a flexible channel, Proceedings of 2017 14th International Bhurban Conference on Applied Sciences & Technology (IBCAST).