Abstract

A mathematical study for investigating heat transfer for the peristaltic flow of a Newtonian fluid through a nonuniform inclined channel while considering surface roughness of the inner wall has been conducted. Special attention is focused to analyze the physical behaviour of various flow characteristics at different magnitudes of surface roughness parameter. The results obtained for temperature have been characterized by two boundary conditions: uniform surface temperature and uniform heat flux. Effect of inner wall surface roughness and angle of inclination on the heat transfer for the peristaltic flow through a nonuniform channel has been discussed in detail. The expression for pressure gradient, velocity, and the temperature was obtained. The effect of inner wall surface roughness and the angle of inclination, on pressure rise and temperature, has also been elaborated in detail. The problem formulation developed has been solved using the long wavelength and low Reynolds number assumption.

1. Introduction

Peristaltic transport has gained the attention of researchers during the last five decades due to its applications in industry and physiology. Peristaltic motion involves a continuous progressive wave of area shrinking and enlargement, along the length of an expandable duct or tube. Peristaltic motion is an involuntary neuromuscular function, which occurs in several body organs of living beings, like the movement of fodder in the esophagus, chyme transport, urine deliver to the bladder through the ureter, transport mechanism in the cervical canal, ova transportation from the ovary to the uterus, fluid flow in ductus efferent, blood movement from small blood arteries, and lymph transportation through the lymphatic vessels. The principle involved in this type of motion is used for designing various contemporary medical equipments such as roller-finger pumps and pumping of blood in the heart-lung machine. This mechanism is applied to the transport of toxic liquids in nuclear installations to prevent the mixing of radioactive elements in the environment. Apart from the above, this phenomenon finds applications in slurry transport, corrosive chemicals, and sanitary liquids from one point to the other.

Profound knowledge of the subject of peristaltic transport can be obtained from the vast literature available on the same. The first study in this area was conducted by Latham [1], which was further elaborated by Shapiro et al. [2] in detail. An effect of nonuniformity of the tube as well as the channel was observed in detail by Gupta and Seshadri [3] taking vas deferens as a reference. Heat transfer analysis of a toxic liquid subjected to peristaltic transport through a channel in nuclear operating power plants was carried out by Radhakrishnamacharya and Murty [4]. Excellent investigations have been done by El Koumy et al. [5] and Mekheimer et al. [6] to study the effect of constant magnetic field on peristaltic transport of electrically conducting Maxwell fluid, through a two-dimensional channel and microchannel, having porosity, with the consideration of Hall effect. An important study to find the effect of rotation on couple stress fluid flowing through a channel by peristaltic mechanism was done by Abd Elmaboud et al. [7]. Tripathi and Anwar Bég [8] studied peristaltic transport (chyme dynamics) through a diseased intestine which acts as a nonuniform porous medium. Two eminent investigations, for compressible, electrically conducting Jeffrey fluid and non-Newtonian nanofluid in a microfluidic and nonuniform channel considering the combined influence of the magnetic field and rheological properties with the consideration of slip velocity were conducted by Abdelsalam and Vafai [9] and Abdelsalam and Bhatti [10]. The collective impact of magnetic field, heat transfer, and porosity on the magneto-hydrodynamic peristaltic motion of the fluid through a channel has been investigated by Eldesoky et al. [11], and a decrement in temperature profile was found due to inclusion of fluid suspension. Victor and Shah [12, 13] found a solution for the energy equation of blood flow in a tube that is subjected to uniform wall temperature and heat flux in a fully developed region of velocity and temperature. The peristaltic mechanism exists in several body organs. The complications of blood flow using the peristaltic flow mechanism in a uniform and nonuniform channel as well as in tube were analyzed by Srivastava and Srivastava [14], supposing blood as a Casson fluid. Mekheimer [15] performed an analysis of blood flow during peristaltic motion through a nonuniform channel subjected to a magnetic field. An outstanding observation in the literature [16] is done to find the behavioural change in a small blood artery in which fluid with the suspended particle is transported by the peristaltic mechanism. Later, Abd Elmaboud et al. [17] conducted a study for the electromagnetic flow of two immiscible fluids, one permeable matrix soaked in a Newtonian fluid and the other pure fluid passes through the inclined channel, and it was found that electric field enhances velocity profile. Admirable investigations were done in studies by Abdelsalam et al. [18], Sohail et al. [19], and Abdelsalam and Sohail [20] to find the impact of heat and mass transfer during the transportation in the effect of the magnetic field through a two-dimensional channel and 3D flow. Sohail et al. [21] addressed the effect of momentum, heat, and mass transfer on the mixed convection boundary layer of Casson liquid for a porous medium. Sadaf and Abdelsalam [22] revealed the hostile effect of nanofluid in a wavy nonuniform annulus with flexible boundaries. Some other worthy research studies by Abdelsalam and Bhatti [23], Abdelsalam and Bhatti [24], Abdelsalam and Bhatti [25], and Eldesoky et al. [26] have been performed with magnetic field effect in a tube for blood flow. Biological organs in which biofluids flow have a small degree of roughness or uncertain (random) spatial variability. Park et al. [27] calculated the effective blood viscosity, with the observation of endothelium roughness of coronary arteries. The degree of roughness within the tissues is very helpful in recognizing the diseases. Evaluation of surface roughness of coronary arteries, in which blood flow takes place through the peristaltic mechanism, helps in developing a remedy for coronary heart disease. Burton and Espino [28] investigated the consequences of intense physiological loading on surface roughness; from the observations, it was revealed that initially, mechanical overloading has no significant effect on roughness; however, when a chemical processing was applied on the surface, an increase in surface roughness was noticed in the longitudinal and circumferential directions. This intense overloading results in damage to the arteries, which is notable after dehydration of the tissue.

Several analyses have been performed in the present area for mini and microchannels. Bhatti et al. [29] conducted a study on the swimming pattern of gyrotactic microorganisms in a blood-based nanofluid flowing through the anisotropically narrowing artery. Veatch [30] carried out experiments on eleven rectangular passages of different cross sections, sizes, and roughness; experiments were performed using four different fluids, and it was concluded that, for laminar flow, friction factor and Reynolds number are independent of roughness. Wall roughness has a significant impact on the heat transfer rate during fluid motion. Study of the impact of wall roughness while considering the roughness orientation, roughness structure, and shape of the roughness element has become a key area of research due to its vast applications in engineering. Excellent studies in the field of heat transfer for the roughed wall have been performed by Han et al. [31] and Taslim et al. [32]. Taylor et al. [33] found the influence of roughness on the fluid motion for surfaces having up to 14% relative roughness, and constricted flow diameters were not considered for analysis during the study; however, till that time, the studies were limited to 5% of relative roughness, which violated micro- and macrofluidic flow systems. Wagner and Kandlikar [34] and Dharaiya and Kandlikar [35] investigated experimentally the influence of sinusoidal structured roughness on heat transfer rate during liquid transport through micro- and macrochannels; experiments were performed for different pitches and roughness heights for a correct understanding of frictional losses in fully developed laminar flow. Kharati-Koopaee and Zare [36] and Attalla et al. [37] performed a study to find the heat transfer effect during fluid flow in microchannels having a rectangular cross section with aligned and offset structured sinusoidal roughness. From the above study, it was concluded that, to lower the pressure loss for the fluid under consideration, high roughness height or low channel height gave favourable results.

From the above, it is evident that, though many investigations have been undertaken for explaining the variation in heat transfer rate due to surface roughness during fluid flow, the variation in heat transfer rate due to wall roughness for peristalsis mechanism during transport of fluid has yet to be explored. The current analysis involves the study of heat transfer, while pondering the wall roughness for a two-dimensional nonuniform inclined channel when a Newtonian fluid flows through it under the influence of peristaltic waves.

2. Problem Formulation

Figure 1 shows the flow geometry, for the laboratory frame of reference. It is supposed that the walls of the nonuniform channel have sinusoidal-shaped roughness. The study has been carried out in the rectangular coordinate system .

The wall geometry iswithwhere represents the half channel width at entry, is the amplitude, is the wavelength, is the wave propagation velocity, is the nonuniformity parameter, is the pitch, is the time, is the axial variable, and is the roughness height.

Continuity equation is as follows:

Momentum equation is as follows:with boundary conditionswhere and are the velocity components for and directions, respectively, in the fixed frame of reference, is the pressure, is the density, is the kinematic coefficient of viscosity, is the specific heat at constant pressure, is the temperature of fluid, is the thermal conductivity, is the temperature at the centre of channel, is the temperature at the wall, and is the angle of inclination.

3. Analysis

The nondimensional parameters are as follows:where is the Reynolds Number, is the wavenumber, is the dimensionless temperature in the case of uniform heat flux, is the dimensionless temperature in the case of uniform surface temperature, is the Eckert number in the case of uniform heat flux, is the Eckert number in the case of uniform surface temperature, and is the Prandtl number.

Using equation (8) in equations (3)–(6) and dropping primes with the supposition of long wavelength and low Reynolds number,with boundary conditionswhere is the nonuniformity of the wall, is the amplitude ratio, is the roughness parameter, and is the dimensionless pitch.

On solving equations (9)–(12) with given boundary conditions equation (13), we get

The flow flux is found asor

As the roughness parameter and angle of inclination , the result obtained from equation (16) is the same as that found by Gupta and Seshadri [3].

The pressure rise across the channel of length in dimensionless forms is

4. Results and Discussion

In the present study, graphs have been drawn using MATHEMATICA software. The instantaneous flow rate q(x, t), periodic in (x − t), is as follows [3]:and magnitude of various parameters is taken as [14, 15], and [28].

Figure 2 depicts changes in the pressure rise as a function of time over one period of the wave, at different surface roughness parameters and channel inclinations, at a particular value of all other parameters like . It is obvious from Figure 2 that, during one period of the wave, an increase in the angle of inclination enhances the pressure rise. The effect of surface roughness is most significant at the time , which corresponds to the point where maximum occlusion occurs from the entry of the channel. It is seen that pressure rise increases with an increase in the roughness parameter for a particular choice of all other parameters. This can be interpreted that, with increasing surface roughness at any particular section, the cross section area will reduce due to which the pressure rise increases with increasing surface roughness at that section.

Figure 3 displays the changes in pressure rise against flow rate for various surface roughness parameters at different angles of inclinations, at the particular value of all other parameters . The magnitude of the pressure rise decreases linearly with an increase in flow rate for all magnitude of the surface roughness parameter and angle of inclination. In a noninclined channel, the maximum flow rate occurs up to , and increase in inclination from to enhances the flow rate. At any particular value of flow rate in-between and , pressure rise increases continuously; when the surface roughness parameter increases from to , reverse impact on the value of pressure rise is seen with increasing surface roughness parameter. The pressure rise increases also with an increase in the angle of inclination from to . The pressure rise is maximum at no flow condition, and this value rises with the surface roughness parameter and angle of inclination. It can be observed that, in the peristaltic pumping zone, the maximum flow rate is obtained at zero pressure rise. From Figure 3, it is also observed that the magnitude of the maximum flow rate can be enhanced by increasing the angle of inclination for a particular choice of all other parameters.

The amplitude ratio has a significant impact on the pressure rise. Figure 4 shows that, at higher magnitudes of amplitude ratio , a rapid augment is visualised in the magnitude of pressure rise. With an increase in inclination from to , the pressure rise increases at any particular value of amplitude ratio.

Figures 5 and 6 depict the variation in pressure rise against flow rate at different values of and , respectively, for different surface roughness parameters. It can be noticed from Figures 5 and 6 that the pressure rise increases with an increase in surface roughness parameter at zero flow rate, and this difference goes on reducing as flow rate increases.

From Figure 5, it is also observed that, for different values of surface roughness parameter, the magnitude of maximum pressure rise and maximum flow rate increases with an increase in . Again from Figure 6, it is seen that the value of maximum pressure rise and maximum flow rate decreases with an increase in for different values of surface roughness parameter. Variation in flow rate and pressure rise with has opposite nature to that of .

Figures 7(a) and 7(b) and Table 1 describe the effect of the angle of inclination on the temperature distribution for different values of surface roughness parameter; similarly, Figures 8(a) and 8(b), Table 2, Figures 9(a) and 9(b), and Table 3 describe the effect of and on the temperature distribution at different values of surface roughness parameter. The entire analysis performed above has been carried out under two conditions, uniform heat flux and uniform surface temperature condition.

Figures 7(a) and 7(b) and Table 1 describe the effect of the angle of inclination on the temperature distribution for different values of surface roughness parameter; similarly, Figures 8(a) and 8(b), Table 2, Figures 9(a) and 9(b), and Table 3 describe the effect of and on the temperature distribution at different values of surface roughness parameter. The entire analysis performed above has been carried out under two conditions: uniform heat flux and uniform surface temperature condition.

Figure 7(a) shows the variation in against for several values of channel inclination and surface roughness. As increases, increases from a minimum value at the centre of the channel to a maximum value at the channel wall, for a channel having zero inclination. For a noninclined channel, it is visible from Table 1 that increases and decreases at any particular value of , as the surface roughness parameter is increased. It is also concluded from Table 1 that, for a particular value of surface roughness parameter in a noninclined channel, when increases, increases and decreases.

With the variation in against , for channels inclined at and inclination, initially increases rapidly to a maximum value after which it starts decreasing; finally, it coincides with the maximum value of for a channel having zero inclination, at the wall of the channel. Furthermore, for channels inclined at and , from Table 1, it is found that the value and decreases at any particular value of , as the value of surface roughness parameter increases. For any particular value of surface roughness, with an increase in , increases and decreases for a channel having inclination. It is also clear from Table 1 that, for any particular value of surface roughness, increases as increases from 0 to 0.8, after which it starts decreasing and increases as increases from 0 to 0.4, after which it starts decreasing for a channel having inclination.

Figure 7(b) shows the variation in against for different values of channel inclination and surface roughness. The value of decreases with increase in and at ; is maximum, for channels having zero and inclination. But in the case of channels having inclination, with increase in , increases up to a maximum value, and then, it decreases to zero at the channel wall.

Figures 8(a) and 8(b) show the temperature profile for different values of and surface roughness parameters for both boundary conditions mentioned above, i.e., uniform heat flux and uniform surface temperature condition. The above is also described numerically through data presented in Table 2. Figures 9(a) and 9(b) show the temperature profile for different values of and surface roughness parameters, and the same is described numerically through the data presented in Table 3. From the above two figures, it is clear that temperature variation at different values of Re and Fr follows a similar pattern to that depicted in Figures 7(a) and 7(b) except in the case of a channel having inclination, wherein initially increases and then it decreases. The variation in for different values of Re and is similar in nature; decreases with increasing for all values of Re and . From Table 2, it is clear that, at Re = 0.1 and 0.2, for any particular value of , goes on increasing, while goes on decreasing upon increasing surface roughness parameter. As the value of the surface roughness parameter is incremented, at Re = 0.3, for a particular value of , decreases for all values of from 0 to 0.8, and also decreases for all values of . Excluding the value of at , the variation in and , for any particular value of , with increasing values of surface roughness parameter at Re = 0.4 depicts a similar nature to that observed for Re = 0.3. For all values of Re at any given value of surface roughness parameter, as increases, increases while decreases.

From Table 3 it is found that, except at for all other values of , for any particular value of , increases and decreases upon incrementing the value of surface roughness parameter. At  = 0.1, for any particular value of , with an increase in surface roughness, both and decrease. For a particular choice of surface roughness parameter, as increases, increases and decreases for all values of Fr.

It can be seen that temperature profiles show a reverse behaviour at uniform surface temperature and uniform heat flux condition. For better understanding, the effect of surface roughness on temperature for various values of channel inclination Re and is shown below in the tabular form.

5. Conclusion

The objective of this analysis is to explain the variation occurring in heat transfer due to the inner wall roughness of a nonlinear inclined channel, for transportation of Newtonian Fluid through the peristaltic mechanism. From the literature survey conducted during the study, it has been found that, so far, no attempt has been made by researchers to address the above problem. The mathematical results obtained in the present study may apply to the transportation of liquids used in various industries. This study may be useful for analysis of heat conduction in many body organs in which physiological fluid transported. The conclusion drawn from the analysis is as follows:(1)An increase in the maximum magnitude of pressure rise is observed when the surface roughness parameter increased.(2)The effect of surface roughness on pressure rise is majorly seen at the time , which corresponds to the point where maximum occlusion occurs.(3)The maximum pressure rise value increases, as the angle of inclination is increased from to .(4)Pressure rise increases at all points for the period of a wave as the inclination of the channel is increased from to .(5)For a particular choice of pressure rise value, the flow rate increases with angle of inclination from to .(6)The significant effect of the angle of the inclination of the channel has been observed in the region 0.6 <  < 0.94.(7)Maximum flow rate increases with an inclination of the channel from to .(8)The value of maximum pressure rise and flow rate rises when Reynolds number is augmented.(9)For any value of Reynolds number, the maximum pressure attains higher values when the surface roughness parameter is increased.(10)The maximum pressure rise and maximum flow rate decrease as the Froude number is increased.(11)For any particular Froude number, maximum pressure rise attains a higher magnitude at a higher surface roughness parameter.(12)The temperature profiles show opposite behaviours at uniform heat flux conditions to that of uniform surface temperature conditions.(13)Excluding the value of at zero inclination, the value of and decreases with an increase in surface roughness, at all angles of inclination for any particular value of .

Nomenclature

:Vertical displacement
:Amplitude of the wave
:Inlet half channel width
:Nonuniformity parameter
:Wave velocity
:Time
:Wavelength
:Roughness height
:Pitch
:Axial variable
:Transverse variable
:Axial velocity
:Transverse velocity
:Density
:Pressure
:Coefficient of viscosity
:Gravitational acceleration
:Angle of inclination of channel
:Kinematic coefficient of viscosity
:Specific heat at constant pressure
:Thermal conductivity
:Temperature of fluid
:Temperature at the wall
:Temperature at centre of the channel
:Dimensionless temperature in case of uniform heat flux
:Dimensionless temperature in case of uniform surface temperature
:Eckert number in case of uniform heat flux
:Eckert number in case of uniform surface temperature
:Prandtl number
:Dimensionless vertical displacement
:Dimensionless time
:Roughness parameter
:Amplitude ratio
:Dimensionless axial variable
:Dimensionless transverse variable
:Dimensionless axial velocity
:Dimensionless transverse velocity
:Dimensionless pressure
:Reynolds number
:Froude number
:Wavenumber
:Pressure rise
:Friction force at the wall of channel
:Length of channel
:Instantaneous flow rate
:Time mean flow rate.

Data Availability

The data used to support the findings of this study are available from the first author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.