# Mathematical Analysis on Heat Transfer during Peristaltic Pumping of Fractional Second-Grade Fluid through a Nonuniform Permeable Tube

Authors: Siddharth Shankar Bhatt; Amit Medhavi; R. S. Gupta
Journal: Journal of Fluids
 (2016)

---

## Abstract

This mathematical study is related to heat transfer under peristaltic flow of fractional second-grade fluid through nonuniform cylindrical tube with permeable walls. The analysis is performed under low Reynolds number and long wavelength approximation. The analytical solution for pressure gradient, friction force, and temperature field is obtained. The effects of appropriate parameters such as Grashof number, nonuniformity of tube, permeability of tube wall, heat source/sink parameter, material constant, fractional time derivative parameter and amplitude ratio on pressure rise, friction force, and temperature distribution are discussed. It is found that an increase in amplitude ratio and material constant causes increase in pressure but increase in nonuniformity of the tube causes decrease in pressure. It is also observed that variation of friction force against flow rate shows opposite behavior to that of pressure. Increase in temperature is also observed due to increase in heat source/sink parameter at inlet as well as downstream.


---

## Body

## 1. Introduction

The word “peristalsis” originated from the Greek word “peristaltikos,” which means clasping and compressing. Peristalsis is a mechanism of fluid transport through deformable vessels with the aid of a progressive contraction/expansion wave along the vessel. It is an important mechanism of fluid transport in different parts of the entire physiological system. Peristaltic flow appears in urine transport from kidney to bladder, the movement of spermatozoa in the ductus efferentes of the male reproductive tract and in the vasomotion of small blood vessels such as arterioles, venules, and capillaries. Peristaltic motion finds application in industry such as heart lung machine and roller pump. Peristaltic motion was first studied clinically in an article given by Bayliss and Starling [1] and much later Latham [2] theoretically investigated peristalsis using fluid mechanics principles. The work of Jaffrin and Shapiro [3] throws light on various parameters involved in the analysis of peristaltic pumping. Kumar et al. [4] studied unsteady peristaltic pumping in finite length tube with permeable wall.Interaction of peristalsis with heat transfer plays an important role in biomedical science. Models of microvascular heat transfer are being increasingly used for optimizing thermal therapies such as hyperthermia treatment. Thermodynamic aspects of blood become significant in processes like oxygenation and hemodialysis. Victor and Shah [5] studied heat transfer to blood flowing in a tube. Srinivas and Kothandapani [6] investigated peristaltic transport in an asymmetric channel with heat transfer. Muthuraj and Srinivas [7] studied mixed convective heat and mass transfer in a vertical wavy channel with traveling thermal waves and porous medium. Lots of investigations have been done for uniform channel or tube; however, most physiological vessels, for example, ureters, esophagus, intestine, and ductus efferentes of the reproductive tract, possess nonuniform geometries. Some researcher have used nonuniform geometry for analysis such as Radhakrishnamacharya and Radhakrishna [8] who discussed heat transfer to peristaltic transport in a nonuniform channel. Ellahi et al. [9] discussed effect of heat and mass transfer on peristaltic flow in a nonuniform rectangular duct. However, the study of heat flow during peristalsis has not been given much consideration by investigators.Most physiological fluids exhibit both viscous and elastic properties. Viscoelastic models are derived by using classical mechanics laws, that is, Newton’s law for viscous liquids and Hooke’s law for elastic solids. In modern days, fractional calculus is a rapidly growing field of research in physics, biology, and medical engineering. By using fractional calculus, the viscoelastic behavior of fluid can be successfully explained. Fractional second-grade calculus operator is actually a generalization to deal with integrals and derivatives of noninteger order. Fractional second-grade model can be obtained by replacing ordinary time derivative to fractional time derivative. This model is applied to the study of movement of chyme through small intestine, esophagus, and so forth. Using the concept of fractional calculus some visoelastic models have been developed such as fractional Maxwell, fractional Zener, fractional anti-Zener or fractional Jeffrey, and fractional Burgers’ models. Number of researchers has studied unsteady flows of viscoelastic fluids using different models like fractional Maxwell model, fractional second-grade fluid model, fractional Burgers’ model and fractional generalized Burgers’ model, and fractional Oldroyed-B model in channels, annulus, or tubes. Qi and Xu [10] studied unsteady flow of viscoelastic fluid with fractional Maxwell model in a channel. Tripathi [11] explored the transportation of a viscoelastic fluid with fractional second-grade model by peristalsis through cylindrical tube under the assumptions of long wavelength and low Reynolds number. Tripathi and Bég [12] studied peristaltic propulsion of generalized Burgers’ fluid through a nonuniform porous medium with chyme dynamics through diseased intestine. Hameed et al. [13] analysed the heat transfer on peristaltic flow of the fractional second-grade fluid confined in a uniform cylindrical tube in presence of magnetic field. Rathod and Tuljappa [14] have studied peristaltic flow of fractional second-grade fluid through a cylindrical tube with heat transfer under the assumption of long wavelength and low Reynolds number assumption.The objective of this study is to investigate the effect of various concerned parameters on fractional second-grade fluid with heat transfer and peristaltic flow through nonuniform tube with permeable walls. The effects of these parameters have been studied under long wave length and the low Reynolds number approximation. The problem is solved analytically by the use of fractional calculus. The obtained expressions are utilized to discuss the influences of various physical parameters.


## 2. Basic Definitions

Definition 1.
The Riemann-Liouville fractional integral operator of orderα>0 of a function fx: (0,∞)→R is given by [15](1)fx=1Γα∫0xx-ξα-1fξdξ,α>0,x>0.Definition 2.
The fractional derivative of orderα>0 of a continuous function fx:(0,∞)→R is given by [15](2)Dαfx=1Γm-αddxm∫0xx-ξα-1fξdξ,form-1<α≤m,m∈N,x>0,f∈C-1m, where m=α+1, provided that right-hand side is point-wise defined on (0,∞).Remark 3.
For example,fx=xβ; we quote for β>-1; in (2) one can get(3)Dαxβ=Γβ+1Γβ-α+1xβ-α,giving in particular Dαxβ-n=0, n=1,2,3,…,N, where N is the smallest integer greater than or equal to α.


## 3. Mathematical Modelling

Consider the flow of an incompressible fractional second-grade fluid (as shown in Figure1) due to peristaltic transport induced by sinusoidal wave trains propagating with constant speed c. The temperature of walls of tube is T0. The constitutive equation for viscoelastic fluid with fractional second-grade model is given by(4)s~=μ1+λ~1α∂α∂t~αγ˙,where t~, s~, γ˙, and λ~1 is time, shear stress, rate of shear strain, and material constant, respectively, μ is viscosity, and α is fractional time derivative parameters such that 0<α≤1. This model reduces to second-grade models when α=1 and classical Navier-Stokes model is obtained by substituting λ~1=0.Figure 1
Geometry of the problem.The geometry of wall surface is given by(5)H=a+B′X+bsin⁡2πηX-ct, where a, b, η, c, and t are radius of the tube at inlet, wave amplitude, wave length, wave propagation speed, and time, respectively. B′ is a constant whose magnitude depends on the length of tube.The governing equations of the motion of viscoelastic fluid with fractional second-grade model through inclined tube for axisymmetric flow are given bycontinuity equation:(6)1R∂RV∂R+∂U∂X=0,momentum equation:(7)ρ∂U∂t+U∂U∂X+V∂U∂R=-∂p∂X+μ1+λ~1α∂α∂t~α1R∂R∂U/∂R∂R+∂2U∂X2+ρgα1T-T0,ρ∂V∂t+U∂V∂X+V∂V∂R=-∂p∂R+μ1+λ~1α∂α∂t~α1R∂RV∂R+∂2V∂X2,energy equation:(8)ρcp∂T∂t+U∂T∂X+V∂T∂R=K∂2T∂R2+1R∂T∂R+∂2T∂X2+Q0.We introduce nondimensional parameters(9)V′=Vcδ,U′=Uc,X′=Xη,R′=Ra,λ′=λ~1cη,δ=aη,p′=pa2μcη,Re=ρacδμ,θ=T-T0T0,Pr=μcpK,Gr=ρgα1a2T0μc,t′=ctη,β=a2Q0KT0,h=Ha,where ρ is density of fluid, Q0 is the constant heat, and p, U, V, R, μ, K, α1, Pr, β, and Gr stand for pressure, axial velocity, radial velocity, radial coordinate, coefficient of viscosity, thermal conductivity, coefficient of expansion, Prandtl number, source/sink parameter, Grashof number, respectively.Using nondimensional parameters given in (9) and applying long wavelength and low Reynolds number approximation, (7) and (8) reduce to (after dropping primes) (10)∂p∂X=1+λα∂α∂tα1R∂R∂U/∂R∂R+θGr,∂p∂R=0,∂2θ∂R2+1R∂θ∂R+β=0.The boundary conditions in dimensionless form are given as follows:(11)∂θ∂R=0atR=0,θ=0atR=h,∂U∂R=0atR=0,U=-k∂U∂RatR=h,where k is the slip parameter including slip [16].The nondimensional wall surface geometry,h, is given as(12)h=1+BX+∅sin⁡2πX-t,where B=B′η/a, ∅=b/a is amplitude ratio.On solving (10) with boundary condition (11), the following is obtained:(13)1+λα∂α∂tαU=R24-h24-kh2∂p∂X-Grβ644R2h2-R4-3h4-4kh3,(14)θ=β4h2-R2.The volumetric flow rate is calculated as(15)Q-=∫0h2πRUdR.Using (13) in (15), the following is obtained:(16)1+λα∂α∂tαQ-π=-h48-kh32∂p∂X+4Grβh6192+Grβkh516.Equation (16) gives (17)∂p∂X=1+λα∂α/∂tαQ-/π-2Grβh6/96-Grβkh5/16-h4/8-kh3/2.It is observed that as β→0 and k→0 (12) and (15) reduce to the corresponding result of Tripathi [11].Pressure rise and friction force at wall are given by(18)Δp=∫01dpdXdX,F=∫01h2-dpdXdX.


## 4. Result and Discussion

MATHEMATICA package is used to see quantitative effects of various parameters involved in the result on the pumping characteristic and heat transfer. Extensive computation has been performed to reveal the influence of geometry, hydrodynamical parameters on pressure distributation, frictional force, and heat transfer. Parameters analyzed are frictional parameters, nonuniformity of the geometryB, amplitude ratio ∅, material constant λ, Grashof number Gr, heat source/sink parameter β, and slip parameter k. To discuss the results for the above obtained quantities, the form of the instantaneous volume flow rate Q(X,t), periodic in (X-t), is assumed as follows [17]:(19)Q-X,tπ=Qπ-ϕ22+2B′ηXaϕsin⁡2πX-t+2ϕsin⁡2πX-t+ϕ2sin⁡2πX-t2,where Q/π is time average of flow over one period of the wave. The above form of Q-(X,t) has been assumed in view of the fact that a constant value of Q-(X,t) gives Δp negative, and consequently, there would no pumping action in the tube wall.Figures2–7 are plotted to see the variation of pressure distribution for different physical parameters. Figure 2 shows that pressure increases with increase in fractional parameter α for given flow rate. Maximum pressure is obtained at zero flow rate. It is also observed that pressure decreases with increases in slip parameter. Figure 3 reveals the fact that pressure increases with increase in source/sink parameter β for fixed flow rate. For given pressure, flow rate increases with increase in β in pumping as well as copumping region. Again maximum pressure is obtained for zero flow rate which is a strong trend of peristaltic fluid dynamics. It can be also seen that pressure decreases with increases in slip parameter. Figure 4 shows the effects of Grashof Number Gr. It is clear from Figure 4 that for given pressure, flow rate increases with increase in Grashof Number Gr in pumping as well as co-pumping region. It can be also seen that pressure increases with increase in Gr for given flow rate. Again, it is observed that pressure decreases with increases in slip parameter. The nonuniformity parameter will exert an important influence on pressure distribution in peristaltic regime. Figure 5 shows the effect of change in nonuniformity parameter B on pressure. It can be seen from Figure 5 that, on increasing the value of B, pressure decreases for a given value of flow rate. Pressure decreases with increases in slip parameter at zero flow rate. The variation of pressure against flow rate for various values of amplitude ratio ∅ is presented in Figure 6. It is clear from the Figure 6 that, for given flow rate, pressure increases with increases in ∅. Figure 7 shows the variation of pressure with change in material constant λ. It is evident from the Figure 7 that, on increasing material constant, pressure also increases. Pressure decreases with increases in slip parameter at zero flow rate.Figure 2
Pressure versus averaged flow rate for various values ofα at λ=1, B=0.1, t=0.5, ϕ=0.4, Gr=3, and β=5.Figure 3
Pressure versus averaged flow rate for various values ofβ at λ=1, B=0.1, t=0.5, ϕ=0.4, Gr=3, and α=0.2.Figure 4
Pressure versus averaged flow rate for various values of Gr atλ=1, B=0.1, t=0.5, ϕ=0.4, β=5, and α=0.2.Figure 5
Pressure versus averaged flow rate for various values ofB at λ=1, Gr=3, t=0.5, ϕ=0.4, β=5, and α=0.2.Figure 6
Pressure versus averaged flow rate for various values ofϕ at λ=1, Gr=3, t=0.5, B=0.1, β=5, and α=0.2.Figure 7
Pressure versus averaged flow rate for various values ofϕλ at ϕ=0.4, Gr=3, t=0.5, B=0.1, β=5, and α=0.2.Figures8–13 are plotted to examine the variation of friction force against flow rate for different physical parameters. It is evident from Figures 8–13 that there exist direct proportionality between frictional force (F) and flow rate (Q) irrespective of the parameter varied. Figure 8 shows the variation of frictional force against flow rate with change in fractional parameter α. It is clear from the Figure 8 that frictional force decreases with increase in α and it is observed that frictional force is higher with increase in slip parameter at zero flow rate. Figure 9 depicts the variation of frictional force against flow rate for different values of source/sink parameter β. It can be examined from Figure 9 that frictional force decreases with increase in β for given flow rate. The greatest value of friction force is associated with lowest value of β. Frictional force is higher with increase in slip parameter at zero flow rate. Figure 10 portrays the variation of frictional force against flow rate for different values of Grashof Number Gr. It is observed from Figure 10 that frictional force decreases with increase in Gr for given flow rate. Frictional force is higher with increase in slip parameter at zero flow rate. It is clear from the Figure 11 that frictional force decreases with increase in B. Frictional force is higher with increase in slip parameter at zero flow rate. The effect of amplitude ratio ∅ on pressure rise is illustrated in Figure 12. It is noted that frictional force decreases with increase in ∅ for given flow rate. Frictional force is higher with increase in slip parameter at zero flow rate. The effect of material constant λ on pressure rise is illustrated in Figure 13. It is evident from Figure 13 that frictional force decreases with increase in λ. For low flow rate, frictional force is negative but becomes positive for higher flow rate. Frictional force is higher with increase in slip parameter at zero flow rate.Figure 8
Friction force versus averaged flow rate for various values ofα at λ=1, t=0.5, ϕ=0.4, B=0.1, Gr=3, and β=5.Figure 9
Friction force versus averaged flow rate for various values ofαβ at λ=1, t=0.5, ϕ=0.4, B=0.1, Gr=3, and α=0.2.Figure 10
Friction force versus averaged flow rate for various values of Gr atλ=1, t=0.5, ϕ=0.4, B=0.1, α=0.2, and β=5.Figure 11
Friction force versus averaged flow rate for various values ofB at λ=1, t=0.5, ϕ=0.4, Gr=3, α=0.2, and β=5.Figure 12
Friction force versus averaged flow rate for various values ofϕ at λ=1, t=0.5, B=0.1, Gr=3, α=0.2, and β=5.Figure 13
Friction force versus averaged flow rate for various values ofλ at ϕ=0.4, t=0.5, B=0.1, Gr=3, α=0.2, and β=5.Figures14–16 show the effects of various parameters on temperature. Temperature increases with increase in source/sink parameter β at inlet as well as downstream as portrayed in Figure 14. Temperature is higher at inlet but as we move downstream temperature gradually decreases. Figure 15 illustrates the variation of temperature with amplitude ratio at inlet and downstream. It is clear that temperature decreases with increase in amplitude ratio at inlet and remains constant at downstream. Figure 16 depicts the variation of temperature with nonuniformity parameter B. It is clear that temperature remains unaffected with increase in B at inlet but increases downstream. It is also observed that temperature is higher at inlet.Figure 14
Effect ofβ on temperature at ϕ=0.4, B=0.1, and t=0.5 at inlet as well as downstream.Figure 15
Effect ofϕ on temperature at β=5, B=0.1, and t=0.5 at inlet as well as downstream.Figure 16
Effect ofB on temperature at β=5, ϕ=0.4, and t=0.5 at inlet as well as downstream.


## 5. Conclusion

In this paper, viscoelastic fluid flow with fractional second-grade model for heat effect with peristaltic transport through nonuniform inclined tube with permeable wall. Exact expression for pressure gradient and temperature are obtained under long wavelength and low Reynolds number. Effects of fractional parameter, nonuniformity of tube, permeability of wall, material constant, source/sink parameter, amplitude ratio, and Grashof number are studied. Fractional calculus theory is used to find pressure gradient expression. The following conclusion can be summarized.(1)
The linear relation is found between pressure and flow.(2)
The linear relation is found between friction force and flow.(3)
The pressure function increases with increasing value of fractional parameter, source/sink parameter, Grashof Number, and amplitude ratio.(4)
The variation of friction force against flow rate shows opposite behavior to that of pressure.(5)
Temperature increases with increasing value of source/sink parameter at inlet as well as downstream of the tube.(6)
Temperature remains constant with increasing value of nonuniformity parameter and amplitude ratio at inlet but may increase with increase in nonuniformity and decreases with increases in amplitude ratio at downstream.



---

*Source: 7506953-2016-06-14.xml*