Abstract

The existing analysis deals with heat transfer occurrence on peristaltic transport of a Carreau fluid in a rectangular duct. Flow is scrutinized in a wave frame of reference moving with velocity away from a fixed frame. A peristaltic wave propagating on the horizontal side walls of a rectangular duct is discussed under lubrication approximation. In order to carry out the analytical solution of velocity, temperature, and pressure gradient, the homotopy perturbation method is employed. Graphical results are displayed to see the impact of various emerging parameters of the Carreau fluid and power law index. Trapping effects of peristaltic transport is also discussed and observed that number of trapping bolus decreases with an increase in aspect ratio .

1. Introduction

The applications of peristaltic flows in medical and engineering sciences have attracted the attention of a number of researchers. Applications of peristalsis occur in swallowing food through the esophagus, urine transport from the kidney to the bladder through the ureter, transport of the spermatozoa in the efferent ducts of the male reproductive tract, movement of the ovum in the fallopian tube, movement of the chyme in the gastrointestinal tract, the transport of lymph in the lymphatic vessels, and vasomotion in the small blood vessels such as the arterioles, veins, and capillaries. The peristaltic phenomenon was first discussed by Latham [1]. Based on his experimental theory, numerous researchers have inspected the phenomenon of peristaltic transport under many conjectures [210].

Another fascinating area in connection with peristaltic motion is the heat transfer which has industrial applications like sanitary fluid transport, blood pumps in the heart-lungs machine and transport of corrosive fluids where the contact of fluid with the machinery parts are prohibited. Only a limited attention has been focused to the study of peristaltic flows with heat transfer [1115].

An immense amount of literature is presented on two-dimensional peristaltic flow problems. The study of peristaltic phenomenon in a rectangular channel was first examined by [16]. Based on the theory of [16], several researchers have studied the phenomenon of peristaltic transport in a rectangular duct under various approximations [1723]. In the papers cited above, the phenomena of heat transfer are not taken into account. Keeping in mind the present information, the heat transfer phenomena on the peristaltic flow of non-Newtonian fluid have not been discussed in a three-dimensional channel. So, the aim of the present problem is to discuss the effects of heat transfer on peristaltic flow of a non-Newtonian fluid in a rectangular duct with different wave forms. The governing equations for the three-dimensional rectangular channel are first modeled for Carreau fluid and then simplified under the long wavelength and low Reynolds number approximation. Homotopy perturbation technique is carried out to calculate the analytical solution of the highly nonlinear partial differential equations. The expressions for velocity, pressure rise, pressure gradient, and temperature have been computed and discussed through graphs. Three-dimensional graphical representations for the velocity field are also discussed through graphs. In the end, different wave forms are also used.

2. Mathematical Formulation of the Problem

Let us consider the peristaltic flow of an incompressible Carreau fluid in a duct of rectangular cross section having the channel width 2d and height 2a. We are considering the Cartesian coordinates system in such a way that the is taken along the axial direction, is taken along the lateral direction, and is along the vertical direction of a rectangular duct.

The peristaltic waves on the walls are represented as where and are the amplitudes of the waves, is the wave length, is the velocity of propagation, is the time, and is the direction of wave propagation. The walls parallel to the plane remain undisturbed and are not subjected to any peristaltic wave motion. We assume that the lateral velocity is zero as there is no change in the lateral direction of the duct cross section. Let be the velocity for a rectangular duct. The governing equations for the flow problem are in which is the density, is the pressure, is the time, and is the stress tensor for the Carreau fluid; is the specific heat; and is the temperature. The stress tensor for the Carreau fluid is defined by [17].

Let us define a wave frame moving with the velocity away from the fixed frame by the transformation

Defining the following nondimensional quantities,

Using the above nondimensional quantities in Equations (2), (3), (4), (5) and (7), the resulting equations after dropping the bars can be written as where

Under the assumption of long wavelength and low Reynolds number , Equations (10), (11), (12), (13), (14) and (21) take the form

The corresponding boundary conditions are where , for straight duct and corresponds to total occlusion.

3. Solution of the Problem

3.1. Homotopy Perturbation Method

The solution of the nonlinear partial differential Equations (22) and (23) has been calculated using homotopy perturbation technique. The homotopy perturbation method for Equations (22) and (23) can be defined as [2430]. or

Here, is an embedding parameter which has the range , under the condition that for , we get the initial solution, and for , we seek the final solution. Here, £ is the linear operator which is taken here as . We define the initial guess as

Let us define substituting Equation (30) into Equations (26) and (27), and then comparing the like powers of , one obtains the following problems with the corresponding boundary conditions.

3.2. Zeroth Order System

3.3. First-Order System

3.4. Second-Order System

The resulting series solutions after three iterations are determined using Equation (30) as (when ) and are evaluated as The volumetric flow rate is given by

The instantaneous flux is defined as

The average volume flow rate over one period of the peristaltic wave is defined as

The pressure gradient is obtained from Equations (35) and (37) as

Integration of Equation (38) over one wavelength yields

It is noticed here that the limit (keeping it fixed and ), the rectangular duct reduces to a two-dimensional channel. It is also noticed that when the rectangular duct becomes a square duct.

4. Expressions for Different Wave Shapes

The nondimensional expressions for five considered wave forms are given by [11]. The expression for the triangular, square, and trapezoidal waves are derived from the Fourier series. (1)Sinusoidal wave (2)Multisinusoidal wave (3)Triangular wave (4)Trapezoidal wave (5)Square wave

5. Numerical Results and Discussion

In this segment, the graphical representation of the proposed problem is discussed. With the help of mathematics software Mathematica, the expression for pressure rise and pressure gradient is calculated. In order to see the behavior of pressure rise with volume flow rate for different values of aspect ratio , power law index , amplitude ratio , and Weissenberg number We. It is observed from Figure 1 that in the peristaltic pumping and retrograde pumping regions, the pumping rate decreases with an increase in aspect ratio β, while in the copumping region , the behavior is quite opposite; here, the pumping rate increases with an increase in aspect ratio . Figures 2 and 3 show the variation of pressure rise with volume flow rate for different values of and . It is observed from Figure 2 that in the retrograde pumping and peristaltic pumping regions, the pumping rate increases with an increase in and , while in the copumping region , the pressure rise decreases with an increase in values of and . In order to see the behavior of pressure rise with volume flow rate for different values of , Figure 4 is plotted. It is depicted from Figure 4 that in the peristaltic pumping and retrograde pumping regions, the pumping rate decreases with an increase in , while in the copumping region , the behavior is quite opposite; here, the pumping rate increases with an increase in . The pressure gradient for different values of aspect ratio β, , amplitude ratio , and We against space variable is plotted in Figures 58. It is depicted from Figures 58 that for and , the pressure gradient is small, i.e., the flow can easily pass without the imposition of a large pressure gradient, while in the region that the pressure gradient decreases with an increase in aspect ratio β and and increases with an increase in amplitude ratio and so much pressure gradient is required to maintain the flux to pass. Figures 912 show the stream lines for different wave forms. It is observed from Figures 912 that the pressure gradient retains the same shape of the waves as we consider.

Figures 1316 show the temperature profile for different values of , , Pr, and . It is depicted from Figure 13 that the magnitude value of the temperature profile increases with an increase in . It is observed from Figure 14 that the magnitude value of the temperature profile decreases with an increase in Figures 15 and 16 show the temperature profile for different values of and . It is observed from Figures 15 and 16 that temperature profile increases with an increase in and . Velocity profile for different values of , , and is plotted in Figures 1719. It is observed from Figures 17 and 18 that the magnitude value of the velocity profile decreases with an increase in volume flow rate . It is depicted from Figure 19 that the magnitude value of the velocity profile increases with an increase in . Stream lines for different values of aspect ratio amplitude ratio , power law index , and Weissenberg number is plotted in Figures 2023. It is observed from Figure 20 that the number of the trapping bolus decreases with an increase in aspect ratio . It is observed from Figures 21 and 22 that the size of the trapping bolus increases with an increase in power law index and amplitude ratio Figure 23 shows the stream lines for different values of the Weissenberg number . It is depicted from Figure 23 that the number of the trapping bolus increases in the lower half of the channel with an increase in , while in the upper half of the channel the size of the trapping bolus increases. Figure 24 shows the stream lines for different wave forms.

6. Conclusion

We have discussed the impact of heat transfer on peristaltic flow of a Carreau fluid in a rectangular duct with different wave forms. The governing equations for a three-dimensional rectangular channel are the first model for Carreau fluid and then simplified under long wave length and low Reynolds number approximation. The analytical solutions of highly nonlinear equations are calculated. The results are discussed through graphs. It is observed that pumping rate decreases with an increase in aspect ratio and in the peristaltic and retrograde pumping regions, while in the copumping region , the behavior is quite opposite. It is also noted that in the retrograde pumping and peristaltic pumping regions, the pumping rate increases with an increase in and . It is observed that the profile of temperature increases with an increase in , , and , while temperature profile shows the opposite behavior in the case of . It is also observed that number of trapping bolus increases in the lower half of the channel with an increase in , while in the upper half of the channel, size of the trapping bolus increases.

Nomenclature

,:Amplitudes of waves
:Velocity of propagation
:Direction of wave propagation
:Specific heat
:Temperature
:Dimensionless wave number
:Aspect ratio
:Weissenberg number
:Power law index
:Temperature of fluid in dimensionless form
:Wavelength
:Time
:Density
:Pressure
:Reynolds number
:Volume flow rate
:Prandtl number
:Amplitude ratio
:Eckert number
:Stress tensor for Carreau fluid.

Data Availability

Data used to support the findings of this work are included in the article.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

The authors would like to acknowledge the support provided by the Deanship of Scientific Research (DSR) at the Prince Mohammad Bin Fahd University (PMU) for funding this work.