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

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 thatanincreaseinamplituderatioandmaterialconstantcausesincreaseinpressurebutincreaseinnonuniformityofthetubecausesdecreaseinpressure.Itisalsoobservedthatvariationoffrictionforceagainstflowrateshowsoppositebehaviortothatofpressure.Increaseintemperatureisalsoobservedduetoincreaseinheatsource/sinkparameteratinletaswellasdownstream.


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 2 Journal of Fluids 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 secondgrade 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.

Mathematical Modelling
Consider the flow of an incompressible fractional secondgrade fluid (as shown in Figure 1) due to peristaltic transport induced by sinusoidal wave trains propagating with constant speed . The temperature of walls of tube is 0 . The constitutive equation for viscoelastic fluid with fractional secondgrade model is given bỹ wherẽ,̃,, 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 The geometry of wall surface is given by Journal of Fluids   3 where , , , , and are radius of the tube at inlet, wave amplitude, wave length, wave propagation speed, and time, respectively. 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 by continuity equation: momentum equation: energy equation: We introduce nondimensional parameters where is density of fluid, 0 is the constant heat, and , , , , , , 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) The boundary conditions in dimensionless form are given as follows: = 0 at = 0, = 0 at = ℎ, where is the slip parameter including slip [16]. The nondimensional wall surface geometry, ℎ, is given as where = / , 0 = / is amplitude ratio.

Journal of Fluids
On solving (10) with boundary condition (11), the following is obtained: The volumetric flow rate is calculated as Using (13) in (15), the following is obtained: Equation (16) It is observed that as → 0 and → 0 (12) and (15) reduce to the corresponding result of Tripathi [11]. Pressure rise and friction force at wall are given by

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 geometry , amplitude ratio 0, material constant , Grashof number Gr, heat source/sink parameter , and slip parameter . To discuss the results for the above obtained quantities, the form of the where / is time average of flow over one period of the wave. The above form of ( , ) has been assumed in view of the fact that a constant value of ( , ) gives Δ negative, and consequently, there would no pumping action in the tube wall. Figures 2-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 on pressure. It can be seen from Figure 5  in Figure 6. It is clear from the Figure 6 that, for given flow rate, pressure increases with increases in 0. 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. Figures 8-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 ( ) and flow rate ( ) 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 . Frictional force is higher with increase in slip parameter at zero flow rate. The effect of amplitude ratio 0 on pressure rise is illustrated in Figure 12. It is noted that frictional force decreases with increase in 0 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. Figures 14-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 . It is clear that temperature remains unaffected with increase in at inlet but increases downstream. It is also observed that temperature is higher at inlet.

Conclusion
In this paper, viscoelastic fluid flow with fractional secondgrade 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. Journal of Fluids (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.