Peristaltic Transport of Prandtl-Eyring Liquid in a Convectively Heated Curved Channel.

Here peristaltic activity for flow of a Prandtl-Eyring material is modeled and analyzed for curved geometry. Heat transfer analysis is studied using more generalized convective conditions. The channel walls satisfy complaint walls properties. Viscous dissipation in the thermal equation accounted. Unlike the previous studies is for uniform magnetic field on this topic, the radial applied magnetic field has been utilized in the problems development. Solutions for stream function (ψ), velocity (u), and temperature (θ) for small parameter β have been derived. The salient features of heat transfer coefficient Z and trapping are also discussed for various parameters of interest including magnetic field, curvature, material parameters of fluid, Brinkman, Biot and compliant wall properties. Main observations of present communication have been included in the conclusion section.

Introduction

Peristaltic transport holds a considerable position in physiology and engineering. Extensive research has been addressed under different situations since the seminal works of Latham [1] and Shapiro et al. [2]. Locomotion of worm, gliding movement of some bacteria, corrosive and sanitary liquids transport, heart lung machine and roller and finger pumps also uses this mechanism for their working. Heat transfer in peristalsis is further significant in chemical and pharmaceutical industries, hemodialysis, oxygenation, tissue analysis, thermotherapy and human thermoregulatory process. Heat transfer is also quite prevalent in several peristaltic pumps. Gul et al. [3] discussed the effect of temperature dependent viscosity on the flow of third grade fluid over vertical belt. Gul et al. [4, 5] also study the heat transfer analysis by considering oscillating vertical and inclined belt. Besides this, advancement is also made about the interaction of magnetohydrodynamics in peristalsis, which finds great importance in connection with certain problems for motion of conductive fluids in physiology, for instance, the blood and blood pumps machines, hyperthermia, cancer therapy, drug delivery transport, magnetic resonance imaging (MRI) and theoretical research about operation of peristaltic magnetohydrodynamic (MHD) compressors. Motivated by all the aforementioned facts the recent investigators are also engaged in the analysis of peristalsis through diverse aspects. Few recent studies and several interesting references in this direction can be seen in the attempts [6-21].

All the aforementioned attempts and existing information on this topic witness that much attention has been given to the flows with peristalsis in a planer channel which seems inadequate in reality. It is because of the fact that most of the ducts in physiological and industrial applications are curved. Influence curvature on peristaltic transport liquid is discussed in only some studies. Sato et al. [22] initially examined the curvature effect in the peristalsis of viscous fluid. Ali et al. [23] reconsidered the problem of ref [24] in wave frame. Hayat et al. [25] study the peristaltic phenomenon for viscous fluid in curved geometry. The peristaltic flow of third order and Carreau-Yasuda materials in curved flow configuration has been discussed by Ali et al. [26] and Abbasi et al. [27] respectively. Heat/mass transport in peristalsis of pseudoplastic, Johnson-Segalman and third grade fluids are explored by Hina et al. [28-29] and Hayat et al. [30]. On the other hand it has also been noted that heat transfer in previous studies related to peristalsis has been dealt with either prescribing temperature or heat flux at the channel walls. Scarce information is available for peristalsis involving heat transfer through convective conditions (see [31-33]).

The facts of present attempt is to advance the theory of peristalsis of non-Newtonian materials via three important aspects i.e. curved channel, convective heat transfer condition and radial magnetic field. Hence we model here the governing flow problem employing constitutive relations of Prandtl-Eyring fluid. The flow formulation is completed through compliant properties of channel walls. Arising nonlinear analysis is computed for the series solutions. Arrangement of paper is as follows. Next section formulates the problems for flow and temperature. Section three includes solution expressions for the stream function, temperature and heat transfer coefficient. Discussion to graph of different physical quantities for various parameters is assigned in section four. section five consists of conclusions.

Formulation

Consider a channel in curved shape with width 2d looped in a circle with radius R* and center O (see Fig 1). Coordinate axis are selected in such a way that x-axis lies along the length of the channel and r-axis lies normal to it. An incompressible electrically conducting Prandtl-Eyring fluid fills the channel. Flow in the channel is generated by propagating peristaltic waves along the channel walls in the axial direction with constant speed c. A radial magnetic field is applied. Induced magnetic field is neglected for small magnetic Reynolds number assumption. Electric field is further absent.

Fig 1

Geometry of the problem.

The wave form at the walls are here a depicts the wave amplitude, λ the wavelength and t the time. The constitutive equations are: where σ is the electrical conductivity of fluid, S (i, j = r, x) the components of extra stress tensor, C denotes the specific heat, ρ the density, κ stands for thermal conductivity and T the temperature of fluid.

Extra stress tensor (S) for Prandtl-Eyring fluid is written as: where L = (grad V) and A/C are material constants of Prandtl-Eyring fluid model.

The boundary conditions have been assumed in the form

Here T is the temperature at both upper and lower walls of the channel, h the heat transfer coefficient at upper/lower walls whereas S, S, S and S are the elements of S. Moreover where m, τ and d are the mass per unit area, longitudinal tension and the viscous damping coefficient respectively. Having

Eqs (3)–(11) can be reduced as follows: with the following non-dimensional boundary conditions with

Taking the stream function ψ (x, r, t) by

Eq (2) is identically satisfied and Eqs (13)–(19) after utilizing lubrication approximation take the forms: where ϵ (= a/d), δ (= d/λ) and K represent the dimensionless amplitude ratio, wave number and curvature parameter respectively. , , the non-dimensional elasticity parameters, Re (= cd /ν) depicts the Reynolds number, Pr (= μC/κ) the Prandtl number, We (= mc/d) the Weissenberg number, the Hartman number, Ec (= c/CT) the Eckert number, Br (= EcPr) the Brinkman number and Bi (= h/d) the Biot number. Further the dimensionless form of extra stress tensor after invoking long wavelength and low Reynolds number approximation becomes with and . Heat transfer coefficient is given by

Solution Methodology

The governing equations are highly non-linear and exact solution seems impossible. Therefore, perturbation method for small parameter β is used to find the solution. Thus we expand ψ, S, θ and Z as follows:

Zeroth order system and solution

Using eqs (29)–(32) into eqs (21)–(28) and comparing the coefficients of β we have with

Solving the above system we get and heat transfer coefficient is given by

First order system and solution

with

The results corresponding to first order system are and heat transfer coefficient is given by

Here the algebraic values of C1 → C6 and B1 → B12 can be evaluated using MATHEMATICA.

Results and Discussion

This portion analyzes the impact of several parameters of interest on the velocity (u), temperature distribution (θ), heat transfer coefficient (Z) and stream function (ψ).

Figs 2–6 are prepared to study velocity profile. It is shown by these figures that velocity profile is parabolic in nature. Also, maximum value is observed at the center of channel. Fig 2 shows decrease in u by enhancing H. This is due to the fact that when the magnetic field is applied in the transverse direction, it provides resistance to the flow which in turn decreases the velocity. Fig 3 shows increase in the axial velocity u increases for larger curvature k near upper half of the channel. Whereas opposite behavior is seen near lower wall. Figs 4 and 5 indicate that the axial velocity acts like an increasing function of Prandtl Eyring fluid parameters α and β. Fig 6 illustrates that with an increase in E and E the velocity enhances. It is due to the fact that less resistance is offered to the flow because of the wall elastance and thus velocity increases. However reverse effect is observed for E. This is because of the fact that larger E more resistive force due to damping and thus velocity decreases. Here we observe that the results obtained are in good agreement with the one get by Hina et al. [27].

Fig 2

Variation of H on u when ϵ = 0.2, x = 0.2, t = 0.1, E1 = 0.2, E2 = 0.01, E3 = 0.1, k = 3.5, α = 1.5 and β = 0.2.

Fig 6

Variation of complaint wall parameters on u when ϵ = 0.2, x = 0.2, t = 0.1, k = 3.5, α = 1.5, H = 2.8 and β = 0.2.

Fig 3

Variation of k on u when ϵ = 0.2, x = 0.2, t = 0.1, E1 = 0.02, E2 = 0.01, E3 = 0.3, H = 0.2, α = 1.5 and β = 0.2.

Fig 4

Variation of α on u when ϵ = 0.2, x = 0.2, t = 0.1, E1 = 0.04, E2 = 0.03, E3 = 0.1, k = 3.5, H = 2.5 and β = 0.2.

Fig 5

Variation of β on u when ϵ = 0.2, x = 0.2, t = 0.1, E1 = 0.04, E2 = 0.03, E3 = 0.3, k = 3.5, α = 1.5 and H = 2.0.

Figs 7–13 indicate the effect of significant parameters involved in the temperature distribution θ. Fig 7 reveals that θ decreases when Hartman number H is increased. Fig 8 shows that increasing curvature K of the channel, we get opposite results for θ in upper/lower half of the channel. Figs 9 and 10 show that the temperature profile increases for larger Prandtl Eyring fluid parameters α and β. It is shown in Fig 11 that temperature increases via E and E and it decreases through E. Fig 12 illustrate that the temperature enhances when Brinkman number Br is increased. Fig 13 discloses that by increasing the Bi the temperature decreases. Here we have considered the values of Biot number much larger than 0.1 due to non-uniform temperature fields within the fluid. Temperature is the average kinetic energy of the molecules. Increase/decrease in temperature directly effects the velocity. Therefore, we get almost similar qualitative behavior for velocity and temperature profiles.

Fig 7

Variation of H on θ when ϵ = 0.2, x = 0.3, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.02, α = 1.5, β = 0.7, Bi1 = 10 and k = 4.

Fig 13

Variation of Bi on θ when ϵ = 0.2, x = 0.3, t = 0.1, E1 = 0.04, E2 = 0.03, E3 = 0.02, α = 1.5, H = 2.0, β = 0.7 and k = 3.5.

Fig 8

Variation of k on θ when ϵ = 0.2, x = 0.3, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.02, α = 1.5, β = 0.7, Bi1 = 10 and k = 0.2.

Fig 9

Variation of α on θ when ϵ = 0.2, x = 0.3, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.01, β = 0.7, H = 2.0, Bi1 = 10 and k = 3.5.

Fig 10

Variation of β on θ when ϵ = 0.2, x = 0.3, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.01, α = 1.5, H = 2.0, Bi1 = 10 and k = 3.5.

Fig 11

Variation of complaint wall parameters on θ when ϵ = 0.2, x = 0.3, t = 0.1, Br = 2, β = 0.7, α = 1.5, H = 2.0, Bi1 = 10 and k = 3.5.

Fig 12

Variation of Br on θ when ϵ = 0.2, x = 0.3, t = 0.1, β = 0.7, E1 = 0.04, E2 = 0.03, E1 = 0.01, α = 1.5, H = 2.0, Bi1 = 10 and k = 3.5.

In Figs 14–20 show the impact of various values of emerging parameters of Z(x). Fig 14 portrays that magnitude of Z(x) decreases when Hartman number H is increased. Fig 15 shows that the Z(x) increases when curvature parameter K is increased. The magnitude of Z(x) increases for larger Prandtl Eyring fluid parameters α and β (Figs 16 and 17). Fig 18 depicts that absolute value of Z(x) increases when there is an increase in E and E. However heat transfer coefficient decreases for E. Fig 19 illustrates that absolute value of heat transfer coefficient enhances by increasing Br. Further Z(x) is increasing function of Bi (Fig 20).

Fig 14

Variation of H on Z when ϵ = 0.2, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.02, α = 1.5, β = 0.7, Bi1 = 10 and k = 4.

Fig 20

Variation of Bi on Z when ϵ = 0.2, t = 0.1, H = 2, E1 = 0.04, E2 = 0.03, E3 = 0.1, α = 1.5, β = 0.7, Br = 2, H = 2 and k = 3.5.

Fig 15

Variation of k on Z when ϵ = 0.2, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.1, α = 1.5, β = 0.7, Bi1 = 5, H = 2.

Fig 16

Variation of α on Z when ϵ = 0.2, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.1, H = 1.5, β = 0.7, Bi1 = 5 and k = 3.5.

Fig 17

Variation of β on Z when ϵ = 0.2, t = 0.1, Br = 2, E1 = 0.04, E2 = 0.03, E3 = 0.1, α = 1.5, H = 2, Bi1 = 5 and k = 3.5.

Fig 18

Variation of complaint wall parameters on Z when ϵ = 0.2, t = 0.1, Br = 2, H = 02, α = 1.5, β = 0.7, Bi1 = 10 and k = 0.5.

Fig 19

Variation of Br on Z when ϵ = 0.2, t = 0.1, H = 2, E1 = 0.04, E2 = 0.03, E3 = 0.1, α = 1.5, β = 0.7, Bi1 = 10 and k = 3.5.

Figs (21–25) display the streamline pattern for various values of invoked parameters. Fig 21a and 21b discusses the impact of Hartman number H on streamlines. Decrease in the size is noticed for increased H. Fig 22a and 22b show the effect of curvature parameter k on the streamlines. These figures show that the bolus size enhances for larger K. Fig 23a and 23b illustrate the fact that the trapped bolus size increases when fluid parameter α is enhanced. Number of streamlines are more. Fig 24a and 24b show that size of trapped bolus decreases when we increase the values of fluid parameter β. We have analyzed through Fig 25a–25d that the streamlines increases through the increase in elastic parameters E and E while increase in E has no show significant effect. Moreover, by taking α = β = H = Bi = 0 results can be obtained for viscous case [24].

Fig 21

Variation of H on ψ for E1 = 0.02, E2 = 0.01, E3 = 0.03, α = 0.2, β = 0.02, ϵ = 0.2, t = 0.0, k = 3.5 when (a): H = 0.8 and (b): H = 1.1.

Fig 25

Variation of wall properties on ψ for H = 5, α = 1.5, k = 3.5, ϵ = 0.2, t = 0.0, β = 0.4 when (a): E1 = 0.03, E2 = 0.02, E3 = 0.1 (b): E1 = 0.04, E2 = 0.02, E3 = 0.1 (c): E1 = 0.03, E2 = 0.04, E3 = 0.1 and (d): E1 = 0.03, E2 = 0.02, E3 = 0.2.

Fig 22

Variation of k on ψ for E1 = 0.02, E2 = 0.15, E3 = 0.05, α = 0.2, β = 0.02, ϵ = 0.2, t = 0.0, H = 0.8 when (a): k = 3.8 and (b): k = 5.

Fig 23

Variation of α on ψ for E1 = 0.02, E2 = 0.15, E3 = 0.05, k = 3.5, β = 0.02, ϵ = 0.2, t = 0.0, H = 2.0 when (a): α = 0.5 and (b): α = 1.1.

Fig 24

Variation of β on ψ for E1 = 0.02, E2 = 0.15, E3 = 0.05, α = 0.7, k = 3.5, ϵ = 0.2, t = 0.0, H = 2.0 when (a): β = 0.3 and (b): β = 1.5.

Conclusions

Peristaltic motion of MHD Prandtl Eyring fluid flowing through curved geometry is discussed. Walls of channel are chosen to be compliant. Heat transfer phenomenon is also analyzed. Main findings of this study are:

Symmetry of velocity profile about the centre line is disturbed for the flow in curved channel.

The velocity profile has decreasing behavior for increasing values of Hartman number H.

Temperature profile is a decreasing function of Biot number Bi.

The absolute value of heat transfer coefficient for planer channel is higher than the curved one.

Size of trapped bolus increases for α but it decreases via β.