Effect of joule heating and MHD on peristaltic blood flow of Eyring Powell nanofluid in a nonuniform

Transcription

1 Journal of Taibah University for Science ISSN: Print 18-3 Online Journal homepage: Effect of joule heating and MHD on peristaltic blood flow of Eyring Powell nanofluid in a nonuniform channel S. K. Asha & G. Sunitha To cite this article: S. K. Asha & G. Sunitha 2019 Effect of joule heating and MHD on peristaltic blood flow of Eyring Powell nanofluid in a non-uniform channel, Journal of Taibah University for Science, 13:1, 1-18, DOI: / To link to this article: The Authors. Published by Informa UK Limited, trading as Taylor & rancis Group Published online: 30 Nov Submit your article to this journal Article views: 3 View Crossmark data ull Terms & Conditions of access and use can be found at

2 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 2018, VOL. 13, NO. 1, Effect of joule heating and MHD on peristaltic blood flow of Eyring Powell nanofluid in a non-uniform channel S. K. Asha and G. Sunitha Department of Mathematics, Karnatak University, Dharwad, India ABSTRACT This paper examines the effect of joule heating and MHD on peristaltic blood flow of a Eyring Powell nanofluid in a non-uniform channel. The transport equations involve the combined effects of Brownian and Thermophoresis diffusion of nanoparticles. The mathematical modelling is carried out by utilizing long wavelength and low Reynolds number assumptions. The energy equation is modelled by taking joule heating effect. The resulting nonlinear system of partial differential equations is solved for the stream function, velocity, pressure gradient, concentration and temperature distributions with the help of the Homotopy Analysis Method. The effect of various physical parameters on the flow characteristics is shown and discussed with the help of graphs. Pressure gradient gives opposite behaviour with increasing values of Eyring Powell fluid parameters A and B. Pressure gradient decreases by increasing Hartman number and joule heating. Nanoparticle concentration increases with increasing thermophoresis parameter, but the reverse trend is observed with the effect of Brownian motion parameter. ARTICLE HISTORY Received 18 September 2018 Revised 9 November 2018 Accepted 13 November 2018 KEYWORDS Peristalsis; Eyring Powell fluid; thermophoresis and Brownian diffusion; MHD; joule heating Nomenclature ũ, ṽ velocity components m/s x, ỹ Cartesian coordinates m t time s p pressure in fixed frame N/m 2 B 0 magnetic field β,c material parameters constant heat addition/absorption κ volume expansion coefficient T temperature of the fluid K C nanoparticle concentration D B Brownian diffusion coefficient D T Thermophoresis diffusion T m fluid mean temperature Re Reynolds number Pr Prandtl number Gr, Qr Grashof numbers due to local temperature and local nanoparticles concentration, respectively Ec Eckert number Br Brinkman number S stress tensor Greek symbols µ dynamic viscosity Ns/m 2 σ electrical conductivity of the fluid S/m δ wave number m 1 ρ f density of the fluid kg/m 3 k α λ thermal conductivity amplitude ratio m wavelength m 1. Introduction The peristaltic transport of Eyring Powell nanofluid in a non-uniform channel is an important problem in engineering, industrial and biological applications. It is useful for the description of the progressive wave contraction along a channel or tube. In particular, the mechanism involved in swallowing food through esophagus, urine transport from kidney to bladder through the ureter, movement of chime in the gastrointestinal tract, vasomotion of small blood vessels and the locomotion of some warms, etc. Peristaltic pump is also used in roller and finger pump. Such devices are used to pump the blood, in bypass surgery it is used to circulate the blood in heart lung machine, slurries is to prevent the transport of fluid from coming in contact with mechanical part of the pump. The word peristaltic comes from the Greek word as peristaltikos which means clutching and compressing. Latham [1] wasthe first who investigated the fluid motion in peristaltic pump. After his first investigation, a large number of researchers and scientist focus their attentions to study the peristaltic flows considering Newtonian and non- Newtonian fluids [2 ]. Bhatti et al. [8] have studied the study of variable magnetic field on the peristaltic CONTACT S. K. Asha as.kotnur2008@gmail.com 2018 The Authors. Published by Informa UK Limited, trading as Taylor & rancis Group This is an Open Access article distributed under the terms of the Creative Commons Attribution License which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

3 1 S. K. ASHA AND G. SUNITHA flow of Jeffrey fluid in a non-uniform rectangular duct having compliant walls. urther, nanofluid is a fluid containing nanometer-sized particles called nanoparticle. The nanoparticles used in nanofluids are typically made of metals Al, Cu, oxides ceramics, Al 2 O 3, CuO; nonmetal like graphite and carbon nanotubes, nanofibres, droplet, nanosheet, etc. Chio [9] was the first to introduce the word nanofluid that represents the fluid in which nanoscale particles diameter less than 100 nm are suspended in the base fluid. Recent articles on the nanofluid are cited in [10 13]. Powell and Eyring proposed a new fluid model in 19 known as Eyring Powell fluid model [1]. Viscoelastic fluids in physiology and industry have a prominent place. The Eyring Powell nanofluid is one such model which is advantageous to recover accurate results of viscous nanofluid at low and high shear rates. At low shear rates less than 100 s 1,thebest example of Eyring Powell nanofluid is human blood. Possibly the complex mathematical description of this model in flexible curved channel averts the study of nanoparticles subject to Eyring Powell fluid [1]. urther investigation of mixed convection peristaltic flow of Eyring Powell nanofluid with magnetic field in a non-uniform channel by Asha and Sunitha [1]. Nooren and Qasim [1] studied the peristaltic flow of MHD Eyring Powell fluid in a channel. Recent articles on the peristaltic transport of Eyring Powell fluid are cited in [18 21]. Peristalsis in connection with nanofluids has application in biomedicines, i.e. cancer treatment, radio therapy, etc. Peristaltic blood flow has gained considerable attention because of its major importance in physiological. The velocity of blood flow in the human body is different depending on what type of vessels it is flowing through. Velocity is indirectly related to the crosssectional area of the vascular bed, so when the crosssectional area is high i.e. in capillaries then the speed of blood flow decreases dramatically. Blood flow follows the laws of hydrodynamics where speed is proportional to pressure difference, fluidity, radius of the tube. When blood flows from heart to various organs it passes through aorta to artery to arterioles to capillaries. Capillaries have very small diameter but as arteriole is divided into many capillaries, the overall diameter increases, so blood volume decreases. So through flow of blood from aorta to capillaries the speed of flow gradually decreases and flow gradually become streamline. Recently, peristaltic flow with magnetic particles grabbed the attention of several researchers due to its emerging applications in magnetic drug targeting [22], pumping of blood, reduction of bleeding during surgery, continue casting process, hyperthermia, magnetic resonance imaging MRI and Magnetotherapy, etc. Moreover, in biomedical engineering, ferrofluid is beneficial for an experimental treatment of cancer disease known as targeted magnetic hyperthermia. urthermore, magneto-hydrodynamic MHD flow and heat transfer towards moving or fixed flat surface have become a vibrant problem in the current field and huge applications in the field of engineering problems, for instance, petroleum engineering, plasma studies, geothermal energy abstractions, aerodynamics and much more. Bhatti et al. [23] have studied the Combine effects of Magnetohydrodynamics MHD and partial slip on peristaltic Blood flow of Ree-Eyring fluid with wall properties. Particularly, to regulate the behaviour of boundary layer, numerous artificial approaches have been established in references [2 2]. Joule heating effect cannot be ignored when strong magnetic field is applied. There are many practical uses of joule heating such as conductors in electronics, electric stoves and other electric heaters, power lines, fuses. The study of heat and mass transfer with joule heating on MHD peristaltic blood flow under the influence of hall effect was given by Bhatti et al. [2]. Recent attempts on mixed convection peristaltic flow with joule effect can be seen in [28,29]. As far the best of author s knowledge, the effect of joule heating on MHD peristaltic blood flow of Eyring Powell nanofluid in a non-uniform channel has not been studied yet. However, Anum et al. [12] studied the mixed convection peristaltic flow of Eyring Powell nanofluid in a curved channel with compliant walls, but did not consider the effect of joule heating. The prime object of the present attempt is to explore the ideas of MHD and joule heating on the peristaltic blood flow by considering the Eyring Powell nanofluid. The present study has numerous applications in biomedical engineering. In most medical therapies, MHD with joule heating plays an important role to modulate the blood and to reduce the pain of human body like in laparoscopic treatment. This paper is organized in the following fashion. In Section 1, the mathematical formulation of the problem is constructed by employing the wave frame and then reduced by assumptions of long wavelength and low Reynolds number approximation. In Section 2, solutions for velocity, temperature and concentration are obtained by using the Homotopy Analysis Method HAM [30,31]. This method is a general approximate analytic solution, used to obtain series solutions of linear/nonlinear equations. This method is valid for the nonlinear problem, contains small or large physical parameters. The method provides great freedom in choosing initial approximation and auxiliary linear operators. By means of choosing such approximation and operators, any complicated nonlinear problems can be solved by transforming into linear subproblems. In the last section, the effects of various physical parameters on pressure gradient, frictional force, concentration and temperature are analysed through graphs.

4 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 1 2. Mathematical formulation Let us consider the peristaltic blood flow of Eyring Powell nanofluid through a two-dimensional nonuniform channel with the sinusoidal wave propagating towards down its walls. The nanofluid is electrically conducting by an external magnetic fieldb 0 and the effect of joule heating. Here we consider the Cartesian coordinate system X, Ỹ suchthat X-axis is considered along the centre line in the direction of wave propagation and Ỹ-axis is transverse. The geometry of the wall surface can be written as 2π h X, t = a X + bsin λ X c t. 1 Here a X = a 20 + k X is the channel half-width, λ is the wavelength, t is the time and b represents the wave amplitude. The velocity components Ũ and Ṽ along the X and Ỹ directions, respectively, in the fixed frame, the velocity field V is taken as V = [Ũ X, Ỹ, t, Ṽ X, Ỹ, t, 0]. 2 The shear of non-newtonian fluid is considered to the study of Eyring Powell fluid model. The stress tensor of Eyring Powell fluid model is given by S = μ Ṽ β sinh 1 c Ṽ, 3 where μ is the coefficient of shear viscosity, β and c are the material constant of Eyring Powell fluid. The term sinh 1 is approximated using the second-order approximation of the hyperbolic sine function as follows: 1 sinh 1 c Ṽ 1 c c Ṽ, 1 c Ṽ 1. The governing equations for the Eyring Powell nanofluid can be written as follows [1]: The continuity equation: Ũ X + Ũ = 0. Ỹ The momentum equation: Ũ ρ f t + Ũ Ũ X + Ṽ Ũ Ỹ = p X + 1 2βc 3 μ Ũ βc X + 2 Ũ σ B 2 0 Ũ 2 Ỹ Ũ Ũ Ũ Ũ + +2 X Ỹ X Ỹ + ρ f gκ T T 0 +ρ f gκ C C 0, 2 Ũ X + 2 Ũ 2 Ỹ 2 Ṽ ρ f t + Ũ Ṽ X + Ṽ Ṽ Ỹ = p Ỹ + 1 2βc 3 μ Ṽ βc X + 2 Ṽ 2 Ỹ Ṽ Ṽ Ṽ Ṽ X Ỹ X Ỹ. 2 Ṽ X Ṽ Ỹ 2 The energy equation: T ρc f t + Ũ T X + Ṽ T Ỹ = k 2 T X + 2 T 2 Ỹ 2 T X + T Ỹ C + ρc p D B X + C Ỹ D T 2 T + ρc p T m X + 2 T 2 Ỹ σ B 0 2 Ũ 2. 8 The concentration equation: C t + Ũ C X + Ṽ C Ỹ 2 C = D B X + 2 C 2 Ỹ 2 + D T 2 T T m X + 2 T, 9 2 Ỹ 2 where p is the pressure, ρ f is the density of the fluid, σ is the electrical conductivity of the fluid, B 0 is the applied magnetic field, is the constant heat addition/absorption, g is the acceleration due to gravity, κ is the volume expansion coefficient, T is the temperature of the fluid, C is the nanoparticle concentration, ρc f is the heat capacity of the fluid, k is the thermal conductivity, ρc p is the effective heat capacity of the nanoparticle material, D B is the Brownian diffusion coefficient, D T is thermophoresis diffusion and T m is the fluid mean temperature. The relations between the laboratory and wave frame are introduced through x = X c t, ỹ = Ỹ, ũ x, ỹ = Ũ c, ṽ x, ỹ = Ṽ, 10 where ũ, ṽand x, ỹ indicate the velocity components and coordinates in the wave frame. Introducing the velocity field in terms of stream functions can be written as u = ψ ψ y, v = δ x. The corresponding boundary conditions for the above problem is

5 18 S. K. ASHA AND G. SUNITHA given by ψ = 0, ũ = ψ ỹ = 0, T = T 0, C = C 0 at ỹ = 0, ψ = q, ũ = ψ ỹ = c, T = T 1, C = C 1 at ỹ = h = a x + b sin 2π λ x. 11 Introducing the following non-dimensional quantities: ψ = ψ ca, B = 1 μβc, A = Bc2 X 2a 2, X = c 2 λ, x = x λ, Y = Ỹ a, y = ỹ a, t = c t λ, p = a2 p cλμ, v = ṽ c, δ = a λ, u = ũ c, Re = 2ρ f ca μ, M = Qr = ρ f gka 2 C 1 C 0 cμ σ μ B 0a, β = k,pr= μ ρc f β, = q ca Gr = ρ f gka 2 T 1 T 0, cμ, Nb = ρc pd B C 1 C 0, ρc f μ Nt = ρc pd T T 1 T 0 a 2, Q 0 =, ρc f T m μ T 1 T 0 ρc f E c = c 2 ρc p T 1 T 0, Br = Pr E c, = C C 0 C 1 C 0, θ = T T 0, h = h = 1 + λkx T 1 T 0 a 20 a 20 +α sin 2πx, α = b. a The non-dimensional symbols of the above-mentioned quantities are as follows: A and B are the nondimensional Eyring Powell fluid parameters, Pr is the Prandtl number, Gr and Qr are the Grashof numbers corresponding to the local temperature and local nanoparticles mass transfer respectively, Nb is the Brownian motion parameter, Nt is the Thermophoresis parameter, the Hartmann number M, Q 0 is the heat generation parameter, Br is the Brinkman number and Ec is the Eckert number, the temperature distribution θ, is the mass concentration, p is the dimensionless pressure, ψ is the stream function, x is the non-dimensional axial coordinate, y is the non-dimensional transverse coordinate, u, v is the velocity components, is the dimensionless average flux in the wave frame, µ is the dynamic viscosity and αis the amplitude ratio. The nonlinear terms in the momentum equation are determined to be zero Reδ 2, where Re = 2ρ f ca μ is the Reynolds number and δ = λ a is the wave number and using non-dimensional quantities. Introducing the velocity fields in terms of stream functions u = ψ ψ y, v = δ x with long wavelength and low Reynolds number approximation, the basic equations 1 12 reduce to p x = 1 + B 3 ψ y 3 A M 2 ψ y 2 ψ y ψ y 3 + Grθ + Br, 13 p y = 0, 1 2 θ θ + Pr Nb θ + Pr Nt y2 y y y 2 + Pr Q 0 + BrM 2 u 2 = 0, 1 2 y 2 + Nt 2 θ = 0. 1 Nb y2 The corresponding dimensionless boundary conditions can be written in the following form: ψ = 0, ψ =, 2 ψ y2 = 0, θ = 0, = 0 at y = 0, ψ y = 1, θ = 1, = 1 at y = h = 1 + λkx + α sin 2πx, a 20 1 where the dimensionless time mean flow rate is the wave frame is related to the dimensionless time mean flow rate in the laboratory frame as follows: = + 1, = h 0 ψ y dy, 18 here = ca Q and = q ca are the dimensionless time mean flow in a fixed and wave form respectively. 3. Method of solution Analysing the governing Equations 13 1 and corresponding boundary conditions 1, it is to find that the solution is an odd function 13 equation and even functions of 1 1. Besides, noticing that Equations 13 1 contain the term y, it is convenient to express uy, θy and y by power polynomials. So, we consider the auxiliary linear operators by using HAM. or HAM solutions, we calculate the initial guesses and auxiliary linear operators in the following form: ψ 0 y = yh3 y 3 + h hy 3, 19 θ 0 y = y h, 20 0 y = y h. 21

6 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 19 The corresponding auxiliary linear operators are with property L ψ = y, L θ = 2 y 2, L = 2 y 2 22 [ y 2 L ψ C 1 + C 2 y + C C y 3 + C y ] = 0, 2 L θ [C 1 + C 2 y] = 0, L [C 1 + C 2 y] = 0, 23 + h 1 + Nt y Nb h + h 3 Nt h θ Nb y Pr Nb h + Pr Nt h 2 + h 2 + Pr Q 0 1 +BrM y BrM y where C 1, C 2, C 3, C and C are arbitrary constants, the zeroth-order deformation equations are 1 ql ψ [ψy, q ψ 0 y, q] = qh ψ h ψ N[ψy, q], 2 1 ql θ [θy, q θ 0 y, q] = qh θ h θ N[θy, q], 2 1 ql [ y, q 0 y, q] = qh h N[ y, q], 2 where q [0, 1] is the embedding parameter, h ψ, h θ and h are auxiliary parameters, L is an auxiliary linear operator, H ψ, H θ and H are auxiliary functions, N is a nonlinear operator, ψy, q, θy, q and y, q are unknown functions, ψ 0 y, q, θ 0 y, q and 0 y, q are initial approximations ofψy, q, θy, q and y, q. rom Equations 2 2, we can easily find the system of equations with their relevant boundary conditions. According to the methodology of the HAM, the solutions are given by θy, q = θ 0 + h θ θ 1 + h θ 2 θ y, q = 0 + h 1 + h ψy, q = ψ 0 + h ψ ψ 1 + h 2 ψ ψ The solutions for velocity u = ψ y, temperature and concentration are evaluated using Equations 2 29 we get 1 uy, q = 3 + y 3 h ψ 2 a 11 + h ψ 3 a 12 + h ψ a 13 + h ψ a 1 + h θ h ψ 3 a 1 + h θ h 3 a 1, 30 θy, q = y h + 2h θ b 11 + h θ 2 b 12 + Pr Nb h θ 2 h h + h θ h ψ 2 b 1 + Pr Q 0 h θ y Nt b 13 + Pr Nth 3 θ b 1 Nb 2, 31 y, q = y h + h Nt y Nb h + h Nt y Nb h The volumetric flow rate in the wave frame is defined by Qx, t = h 0 uy, qdy. 33 The pressure gradient can be calculated from the volumetric flow rate and after some simplification they can be written as Q h + h 3 h 2 ψ c 11 h 3 ψ c 12 h ψ c 13 h ψ c 1 p x = h θ h ψ 3 c 1 h θ h 3 c 1 M 2 20 h h ψ 2 h3 h ψ 3 h3 3 h ψ 1 + B h3 h ψ c 1. 3 The pressure rise P L t and the frictional force L t at the wall in the non-uniform channel of length L and it can be defined in their non-dimensional form. L/λ p P L t = dx, 3 x L t = 0 L/λ. Results and discussion 0 h p dx. 3 x This section prepared to describe the results of our problem graphically using MATHEMATICA program. It consists of four parts: the first part illustrates the impact of some physical parameters on the pressure gradient, the second part illustrates the frictional force. However in the third part is the evaluation of temperature distribution with different parameters discussed and the fourth part illustrates the concentration with different parameters..1. Pressure gradient distribution Peristaltic pumps are used to pump the blood, in bypass surgery it is used to circulate the blood in heart lung machine. Peristaltic pumping mechanism is used to transport various biological fluids in diagnostic problems to analyse the pumping characteristic. igures 1 and 2 are plotted to show the behaviour

7 10 S. K. ASHA AND G. SUNITHA igure 1. Pressure gradient p x versus x when t = 0.1, x = 0.2, Pr =.9,ϕ = 0., Br = 2.0, M = 0.,Q 0 = 2.0,Q = 0.2, λ = 10, k = 0.1, a 20 = 2.0. a B = 2.0, Nt = 0., Nb = 0., Gr = 0.3, Qr = 0.3. b A = 0.001, Nt = 0., Nb = 0., Gr = 0.3,Qr = 0.3. c A = 0.001, B = 2.0, Nt = 0., Nb = 0., Qr = 0.3. d A = 0.001, B = 2.0, Nt = 0., Nb = 0., Gr = 0.3. e A = 0.001, B = 2.0, Nb = 0., Gr = 0.3, Qr = 0.3. f A = 0.001, B = 2.0, Nt = 0., Gr = 0.3, Qr = 0.3. of different parameters on the pressure gradient distribution. igure 1a and b are plotted to show the behaviour of fluid parameters A and B on pressure gradient. It is observed that the pressure gradient has the opposite behaviour on these two fluid parameters. The pressure gradient decreases with an increase in the Eyring Powell fluid parameter A, which appears with the nonlinear term of the governing momentum, in the wider part of channel is relatively small x [0, 0.] and x [0.8, 1.0] the pressure gradient is relatively small and the flow can easily pass without the imposition of large pressure gradient. However, in the narrow part of the channel x [0., 0.8] as much larger pressure gradient is required to maintain the same flux to pass through it. Moreover, it is seen that pressure gradient > 0 appears to increase when fluid parameters B is increased. igure 1c and d show the effect of Grashof number Gr and local nanoparticle Grashof number Qr on pressure gradient. Here it is observed that the pressure gradient increases as Grashof number Gr and local nanoparticle Grashof number Qr increases. igure 1e shows the effect thermophoresis parameter Nt on pressure gradient. The pressure gradient appears to increase when the thermophoresis parameter Nt increased. Whereas, in igure 1f, the pressure gradient decreases when Brownian motion parameter Nb increased. igure 2g captures the effects of Hartman number M on the pressure gradient. Here it is observed that the pressure gradient is decreasing with the increase of M

8 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 11 igure 2. Pressure gradient p x versus x when t = 0.1, x = 0.2, λ = 10, k = 0.1, a 20 = 2.0, Pr =.9, A = 0.001, B = 2.0, Nt = 0., Nb = 0., Gr = 0.3, Qr = 0.3, Br = 2.0, Q 0 = 2.0; g ϕ = 0., Q = 0.2. h M = 0.,ϕ = 0.. i M = 0., Q = 0.2. j ϕ = 0., Q = 0, Q = 0.2. which is very helpful during heart surgery. igure 2h includes flow rate Q the pressure gradient appears to decrease when Q is increased and it is observed that the pressure is maximum at zero flow rate. In igure 2i, the effect of amplitude ratio ϕ on pressure gradient is plotted. The pressure gradient appears to increase when the amplitude ratio ϕ increased. In igure 2j, the effect of flow rate Q and Hartman number M on pressure gradient is plotted. The pressure gradient decreases when the flow rate decreases and Hartman number also decreases..2. rictional force distribution igure 3 is plotted to show the behaviour of different parameters on the frictional force distribution. igure 3a shows that the frictional force decreases with an increase in the Eyring Powell fluid parameter A. igure 3b shows that as Eyring Powell fluid parameter B increases the frictional force increases. igure 3c captures the effects of Hartman number M on the frictional force. Here it is observed that frictional force is a decreasing function of M..3. Temperature distribution igures are plotted to show the behaviour of different parameters on the temperature distribution θ. igure a shows that the temperature decreases with an increase in the Eyring Powell fluid parameter A. igure b shows the effect of Eyring Powell fluid parameter B on frictional force. It is observed that the temperature appears to increase when B is increased. igure c shows that the temperature increases with increasing of Hartman number. In igure d is the effect joule heating Q 0. The temperature appears to increase when joule heating effect Q 0 increased. igure e shows that the temperature increases with an increase in the Brownian motion parameter Nb. igure f captures the effect of thermophoresis parameter Nt on the temperature, the temperature increases when the thermophoresis parameter increases. igure g shows that as Grashof number Gr increases the temperature increases. igure h includes the effect of local nanoparticle Grashof number Qr. Here temperature appears to decrease when Qr is increased. igure i capturesthe effects of Brinkman number Br on the temperature. Here it is observed that temperature increases when Brinkman number Br increased. igure j shows that the temperature increases with an increase in the Prandtl number Pr... Concentration distribution igure describes the variation of the concentration profile for several values of the Hartman number M, joule heating effect Q 0, thermophoresis parameter Nt and the Brownian motion parameter Nb, Brinkman number Br and Prandtl number Pr. In igure a, we observe that the concentration profile increases with

9 12 S. K. ASHA AND G. SUNITHA igure 3. rictional force L t versus Q when t = 0.1, x = 0.2, λ = 10, k = 0.1, a 20 = 2.0, Pr =.9, Nt = 0., Nb = 0., Gr = 0.3, Qr = 0.3, Br = 2.0, Q 0 = 2.0, ϕ = 0., Q = 0.2. a B = 2.0, M = 0.. b A = 0.001,M = 0.. c A = 0.001, B = 2.0. igure. Temperature θversus y when t = 0.1, x = 0.2, λ = 10, k = 0.1, a 20 = 2.0, Pr =.9, Nt = 0., Nb = 0., Gr = 0.3, Qr = 0.3, Br = 2.0, ϕ = 0., Q = 0.2. a B = 2.0, M = 0., Q 0 = 2.0. b A = 0.001, M = 0.,Q 0 = 2.0. c A = 0.001, B = 2.0, Q 0 = 2.0. d A = 0.001, B = 2.0, M = 0..

10 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 13 igure. Temperature θ versus y when t = 0.1, x = 0.2, λ = 10, k = 0.1, a 20 = 2.0, A = 0.001, Q 0 = 2.0, ϕ = 0., Q = 0.2. e Pr =.9, Nt = 0., Gr = 0.3, Qr = 0.3, Br = 2.0. f Pr =.9, Nb = 0., Gr = 0.3, Qr = 0.3, Br = 2.0. g Pr =.9, Nb = 0., Nt = 0., Qr = 0.3, Br = 2.0. h Pr =.9, Nb = 0., Nt = 0., Gr = 0.3, Br = 2.0. i Pr =.9, Nb = 0., Nt = 0., Gr = 0.3, Qr = 0.3. j Nb = 0., Nt = 0., Gr = 0.3, Qr = 0.3, Br = 2.0. an increase in Hartman number M. igure b captures the effect of joule heating Q 0 on the concentration profile, the concentration increases when the joule heating effect Q 0 increases. igure c shows that the concentration profile decreases with an increase in Nb. igure d captures the effect of thermophoresis parameter Nt on the concentration profile, the concentration increases when the thermophoresis parameter increased. igure e includes the effect of Brinkman number Br on the concentration. Here it is observed that the concentration increases when Brinkman number Br increased. igure f shows that the concentration increases with an increase in the Prandtl number Pr.. Concluding remarks Here we analysed the effect of joule heating on peristaltic blood flow of Eyring Powell nanofluid in a nonuniform channel, is modelled in the presence of MHD, is studied under the assumptions of long wavelength and low Reynolds number. The results are discussed through graphs. We have the following main observations of the present analysis. Both pressure gradient and temperature distribution give opposite results with increasing values of Eyring Powell fluid parameters A and B.

11 1 S. K. ASHA AND G. SUNITHA igure. Concentration versus y when t = 0.1, x = 0.2, ϕ = 0., λ = 10, k = 0.1, a 20 = 2.0. a Pr =.9, Nb = 0., Nt = 0., Br = 2.0,Q 0 = 0.. b M = 0., Pr =.9, Nb = 0., Nt = 0., Br = 2.0. c Pr =.9, M = 0., Nt = 0., Br = 2.0,Q 0 = 0.. d Pr =.9, M = 0., Nb = 0., Br = 2.0, Q 0 = 0.. e Pr =.9, M = 0., Nb = 0., Nt = 0., Q 0 = 0.. f M = 0., Nb = 0., Nt = 0., Br = 2.0, Q 0 = 0.. Similar behaviour of pressure gradient and temperature profiles is observed with an increase in thermal Grashof numbergr. Also opposite behaviour of pressure gradient and temperature in the presence of local nanoparticle Grashof number Qr. Both pressure gradient and concentration profiles give opposite results with increasing values of Nt and Nb. Also the temperature increases with an increase in Brownian motion parameter Nb and thermophoresis parameter Nt. The pressure gradient decreases with increasing Hartman number M and flow rate Q. Pressure gradient decreases by increasing amplitude ratio ϕ. rictional force gives the opposite behaviour with increasing values of Eyring Powell fluid parameters A and B, also frictional force decreases with an increasing Hartman number M. Similar behaviour of temperature and concentration profiles increases in the presence of joule heating effect Q 0, Hartman number M, Brinkman number Br and Prandtl number Pr. Acknowledgements We would like to thank the editor and the reviewer s for their constructive comments that significantly improved this manuscript.

12 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 1 Disclosure statement No potential conflict of interest was reported by the authors. unding This work was supported by the University Grant Commission IN: [Grant Number NST-KAR-0121 Dated: 01 April 201]; University Grants Commission: [Grant Number 10/3/DRS-III/201SAP-I Dated:29-eb-201]. ORCID S. K. Asha G. Sunitha References [1] Latham TW. luid motion in a peristaltic pump. MS Thesis, M.I.T. Cambridge. 19. [2] Jaffrin MY, Shapiro AH. Peristaltic pumping. Ann Rev luids Mech. 191;3:13 3. [3] Hayat T, Tanveer A, Yasmin H, et al. Effects of convective conditions and chemical reaction on peristaltic flow of Eyring Powell fluid. Appl Bionics Biomech. 201;11: [] Khan AA, Masood, Ellahi R, et al. Mass transport on chemicalized fourth-grade propagating peristaltically through a curved channel with magnetic effects. J Mol Liq. 2018;28: [] Bhatti MM, Zeenshan A, Ellahi R, et al. Mathematical modelling of heat and mass transfer effects on MHD peristaltic propulsion of two-phase flow througha Darcy Brinckman orcheimer porous medium. Adv Powder Technol. 2018;29: [] Zeeshan A, Ijaz N, Abbas T, et al. The sustainable characteristic of BIO-Bi-phase flow of peristaltic transport of MHD Jeffrey fluid in the human body. Sustainability. 2018;10:1 1. [] Asha SK, Sunita G. Effect of couple stress in peristaltic transport of blood flow by homotopy analysis method. AJST. 201;12:98 9. [8] Bhatti MM, Ellahi R, Zeeshan A. Study of variable magnetic field on the peristaltic flow of Jeffrey fluid in a nonuniform rectangular duct having compliant walls. J Mol Liq. 201;222: [9] Choi SUS. Enhancing thermal conductivity of the fluids with nanoparticles. ASME luids Eng Div. 199;231: [10] Bhatti MM, Zeeshan A, Ellahi R. Simultaneous effects of coagulation and variable magnetic field on peristaltically induced motion of Jeffrey nanofluid containing gyrotactic microorganism. Microvasc Res. 201;110:32 2. [11] Ellahi R, Alamri SZ, Basit A, et al. Structural impact of kerosene-al 2 O 3 nanoliquid on MHD Poiseuille flow with variable thermal conductivity. application of cooling process. J Mol Liq. 2018;2:0 1. [12] Tanveer A, Hayat T, Alsaadi, et al. Mixed convection peristaltic flow of Eyring Powell nanofluid in a curved channel with compliant walls. Comput Biol Med. 201;82:1 9. [13] Prakash J, Sharma A, Tripathi D. Thermal radiation effects on electroosmosis modulated of ionic nanoliquids in biomicrofluidics channel. J Mol Liq. 2018;29:83 8. [1] Powell RE, Eyring H. Mechanism for the relaxation theory of viscosity. Nature. 19;1:2 28. [1] Abbasi M, Alsaedi A, Hayat T. Peristaltic transport of Eyring Powell fluid in a curved channel. J Aerosp Eng 201;2:1. [1] Asha SK, Sunitha G. Mixed convection peristaltic flow of a Eyring Powell nanofluid with magnetic field in a nonuniform channel. JAMS. 2018;28: [1] Noreen S, Qasim M. Peristaltic flow of MHD Eyring Powell fluid in a channel. Eur Phys J Plus. 2013;12891:1 10. [18] Hayat T, Irfan Shah S, Ahmad B, et al. Effect of slip on peristaltic flow of Powell Eyring fluid in a symmetric channel. Appl Bionics Biomech. 201;11:9 9. [19] Abbasi M, Hayat T, Alsaedi A. Numerical analysis for peristaltic motion of MHD Eyring Prandtl fluid in an inclined symmetric cannel with inclined magnetic field. J Appl luid Mech. 201;9: [20] Hina S, Mustafa M, Hayat T, et al. Peristaltic flow of Powell Eyring in curved channel with heat transfer: a useful application in biomedicine. Comput Methods Programs Biomed. 201;13: [21] Hayat A, Abdulhadi A. Influence of magnetic field on peristaltic transport for Eyring Powell fluid in a symmetric channel during a porous medium. Math Theor Model. 201;9:9 22. [22] Ali Abbas M, Bai YQ, Rashidi MM, et al. Application of drug delivery in magnetohydrodynamics peristaltic blood flow of nanofluid in a non-uniform channel. J Mech Med Biol. 201;10:1002. [23] Bhatti MM, Ali Abbas M, Rashidi MM. Combine effects of magnetohydrodynamics MHD and partial slip on peristaltic blood flow of Ree Eyring fluid with wall properties. JESTECH. 201;193: [2] Zeeshan A, Bhatti MM, Akbar NS, et al. Hydromagnetic blood flow of sisko fluid in a non-uniform channel induced by peristaltic wave. CTP. 201;81: [2] Ranjit NK, Shit GC, Tripathi D. Joule heating and zeta potential effect on peristaltic blood flow through porous micro-vessels altered by elctrohydrodynamics. Microvasc Res. 2018;11: 89. [2] Akbar NS, Nadeem S. Characteristics of heating scheme and mass transfer on the peristaltic flow for an Eyring Powell fluid in an endoscope. Int Commun Heat Mass. 2012;: [2] Bhatti MM, Rashidi MM. Study of heat and mass transfer with joule heating on magnetohydrodynamic MHD peristaltic blood flow under the influence of hall effect. Propul Power Res. 201;3:1 18. [28] Abbasi M, Hayat T. Effects of inclined magnetic field and joule heating in mixed convective peristaltic transport of non-newtonian fluids. Bull Pol Ac Tech. 201;33: [29] Hayat T, Aslam N, Rafiq M, et al. Hall and joule heating effects on peristaltic flow of Powell Eyring liquid in an inclined symmetric channel. Results Phys. 201;: [30] Turkyilmazoglu M. A note on the homotopy analysis method. Appl Math Lett. 2010;23: [31] Gupta VG, Gupta S. Application of homotopy analysis method for solving nonlinear Cauchy problem. Surv Math Appl. 2012;: Appendix The values of a 11, a 12, a 1, a 13, a 1, a 1 in Equation 30 are given as 1 a 11 = M 2 y Py2 2

13 1 S. K. ASHA AND G. SUNITHA B 1 A 1 a 12 = 1 + B 1 + B 1 A 2h Qr y 3 2h y M 2 y Py 2 1 y. 2h Qr 2h y 3 a 1 = APM y AM2 P B 2h Qr y 2h 1 + A 2 PM AM B 2h Qr y 3 2h AP 1 + B 2h Qr y 2h 1 + 2AM B 2h Qr y 1 2h + 288AP 1 + B 2h Qr y 1 2h 31 + A 1 + B 2h Qr 2 y 1 2h + 1A2 M y 1 + 2A 2 M 2 P + 1 y 8 + PA y A 2 M B 2h Qr y 9 1 2h A2 P B 2h Qr y 10 2h A B 2h Qr y h A3 M y A3 2 y A B 2h Qr y 1 1 2h A 2 y a 13 = 1 + B M B 3 + y y P1 + B + M B M h Qr 2h B 1 + B 3 + y 3 + PM2 y M B 2h Qr y A1 + B 2h 1 + M 2 A y 9 3. a 1 = Gr y3 h + Gr Pr Nb Nt + Pr 2h2 2h 2 + Pr Q BrM y 12 BrGr 1 y 80. a 1 = Qr 1 + Nt y 3 Nb h h Qr 2h y The values of c 11, c 12, c 13, c 1, c 1, c 1, c 1 in Equation 32 are given as 1 c 11 = h M h B 2 1 2h Qr h 2h 1 2A h 8. 1 c 12 = 1 + B h M h B h 1 2A h 8. 1 c 13 = 1 + B h M B 3 + h M B 3 + h 3 M B 1 + B 2h Qr 2h c 1 = AM 1 + M B 1 2A1 + B 1 + 2M 2 A 2h Qr h 2h 30 h 8 h Qr h 2h + 2A Qr 2h h 2 h 1 + B + 288A 2 + AM B 1 + B h 2 + AM2 2h 1 2h 2h h

14 JOURNAL O TAIBAH UNIVERSITY OR SCIENCE 1 Qr h 2h A2 M A B 1 c 1 = AM 1 h A3 2 h h 1 + AM h Qr h 2h 2 + AM2 1 + B h 10 2h Qr 2h h B 1 + 2A 2h Qr h 2h + 288A 1 + B 2h Qr h 1 2h A2 M h A B 2h Qr 2h 1 1 h c 1 = Gr b3 2 Gr 1 1 h A3 Pr Nb Nt + Pr 2h2 2h 2 + Pr Q BrM2 h Br h c 1 = AM 3 + Qr h 2h + 12AM B 2h Qr h 2h + 2A Qr 2 h 8 2h + 18A2 M A 2 M B A B B h 1 + B 1 2h Qr 2h A 3 2 2h Qr 2h 1 1 h h h h B 2h Qr h 1 2h A 3 h c 1 = Qr 1 + Nt h 3 Nb 2. The values of b 11, b 12, b 13, b 1, b 1 in Equation 33 are given as b 11 = y y2 y2 + Pr Nb + Pr Nt h 2h2 2h 2 + Pr Q y BrM y BrM2 y 8 1 BrM y 28. b 12 = Pr Nb y2 y2 + Pr Nt 2h2 2h 2 + Pr Q BrM y BrM2 1 BrM b 13 = y2 y3 y3 + Pr Nb + Pr Nt 2h2 h2 1 + BrM y 3 + 8BrM2 y y 20. h 2 + Pr Q 0 y BrM y 120. y 8 y 9 31 b 1 = y2 Nb 2h 2 + Pr h 2 + Nt 1 h 2 + Q 0 + BrM y 9 31 y 1 12 BrM y BrM 2 2 y h Nb Pr h 2 + Nt h 2 1 BrM Q 0 + BrM y 3 2 h 1 y BrM h y 9 3 Nb + Pr h 2 + Nt 1 h 2 + Q 0 + BrM BrM y 8 Nb 1 Pr BrM 2 1 h 2 + Nt h 2 + Q 0 y BrM b 1 = 8A y B h 3 + h

15 18 S. K. ASHA AND G. SUNITHA + M2 1 1 h 3 + h A2 1 + B 1 h 3 + h Gr 120h + Qr 2 y h 9 1 h 3 + y 18 h A h 3 + h + M 2 1 h 3 + h + Gr 120h + Qr y 120h A2 1 h 3 + h y A1 + B 181 h 3 + h + M2 1 + Gr 120h + Qr 1 120h h 3 + h y. 1 h 3 + h