Convective heat transfer for Peristaltic flow of SWCNT inside a sinusoidal elliptic duct

Abstract

A mathematical model is presented to analyse the flow characteristics and heat transfer aspects of a heated Newtonian viscous fluid with single wall carbon nanotubes inside a vertical duct having elliptic cross section and sinusoidally fluctuating walls. Exact mathematical computations are performed to get temperature, velocity and pressure gradient expressions. A polynomial solution technique is utilized to obtain these mathematical solutions. Finally, these computational results are presented graphically and different characteristics of peristaltic flow phenomenon are examined in detail through these graphs. The velocity declines as the volume fraction of carbon nanotubes increases in the base fluid. Since the velocity of fluid is dependent on its temperature in this study case and temperature decreases with increasing volumetric fraction of carbon nanotubes. Thus velocity also declines for increasing volumetric fraction of nanoparticles.

Keywords

Peristaltic flow elliptic duct carbon nanotubes Grashof number

Introduction

Peristalsis is the study of transportation of fluid inside ducts having sinusoidally deformable walls. The applications of this peristaltic flow phenomenon include many physiological flow problems (i.e. movement of food inside digestive system, urine transport to bladder and blood circulation in small blood vessels etc). Moreover, there are plenty of engineering applications like transportation of noxious and corrosive fluids, aggressive chemicals and slurries.¹ Many researchers have been interpreting the peristaltic flow problems for different geometries (tubes, rectangular duct, asymmetric channels etc) due to their variety of applications. Pozrikidis² had provided a theoretical investigation for two-dimensional flow in channels having sinusoidally deformable walls. Yin and Fung³ had analysed the peristaltic flow of a viscous Newtonian fluid inside cylindrical tubes. Srinivas and Kothandapani⁴ had studied the combined effect of peristalsis and heat transfer inside an asymmetric geometry. Nadeem and Akram⁵ had presented a mathematical investigation to study the peristaltic flow inside a duct having rectangular cross section. Further, some recent research articles that present the study of peristaltic flow inside different geometries is provided.^6–9 Moreover, a huge literature is available on study of peristaltic flow inside rectangular ducts, cylindrical and asymmetric channels. Recently, the analysis of fluid transport along with heat transfer inside ducts having elliptic cross-section and sinusoidally deformable walls has gained the interest of researchers. These elliptic duct peristaltic flow problems also have many applications in the field like the peristaltic flow problems for rectangular and circular ducts.

When compared to a circular duct, the elliptic duct is more practical to use if a better cooling effect is required for heat exchanger and the design has space as an important factor. The reason behind this advantage is that if an elliptic duct and a circular duct both have equal cross sectional area then the elliptic duct has a longer circumferential length as compared to circular duct.¹⁰ The flow of heated Newtonian fluid inside a geometry with elliptic cross section was mathematically interpreted by Richardson.¹¹ Abdel-Wahed et al.¹² had presented the experimental analysis for heated and fully developed laminar fluid flow inside an elliptic duct. Some recent researches that include the flow in elliptic ducts with heat transfer are Sakalis et al.,¹³ Shariat et al.,¹⁴ and Ragueb and Mansouri.¹⁵

Some fluids have low thermal conductivity and it restricts the analysis of heat transfer in many engineering problems. This issue can be resolved by adding a specific quantity of nano-sized particles in the base fluid. Since the addition of these nanoparticles play an important role to enhance thermal conductivity.¹⁶ Thus, the analysis of fluid transport with nanoparticles has its role in different geometries used for fluid flow problems. Akbar et al.¹⁷ had mathematically examined the peristaltic flow of nanofluids inside a vertical cylindrical duct. The fluid transport with single wall carbon nanotubes inside a curved cylindrical channel is mathematically modelled by Shahzadi et al.¹⁸ Akbar¹⁹ had provided the mathematical computations for the flow of carbon nanotubes across an asymmetrical channel. Ellahi et al.²⁰ had interpreted mathematically the peristaltic flow of nanofluid inside a duct with rectangular cross section. Some researchers have also studied the nanofluid flow inside ducts having elliptic cross section.^21,22 Further, some of the recent researches that interpret the nanofluid flow with distinct applications of practical importance are referred as Ellahi et al.,²³ Maskeen et al.,²⁴ and Zeeshan et al.^25,26

All this literature review has revealed that the peristaltic flow of nanoparticles across a duct having elliptic cross section is not computed mathematically yet. We have mathematically computed the peristaltic flow of a heated Newtonian viscous fluid with single wall carbon nanotubes across a vertical duct having elliptic cross-section. In this study, water is considered as base fluid. The heat transfer analysis is explained by considering the constant heat absorption for base fluid. Exact mathematical computations are done to obtain velocity as well as temperature profile solutions. We have added some graphical results to explain the flow behaviour with heat transfer in detail. Further, different aspects of peristaltic flow are discussed with these graphical results. Streamlines depict the side view of this vertical elliptic duct peristaltic flow problem.

Mathematical model

The peristaltic flow of a heated Newtonian viscous fluid inside a vertical elliptical duct is mathematically interpreted. The Cartesian coordinates $(\bar{X,} \bar{Y}, \bar{Z})$ are used to interpret this problem mathematicall (See Figure 1).

The geometry of deformable walls is provided by the following sinusoidal equations

\begin{matrix} \bar{a} (\bar{Z}, \bar{t}) = a_{0} + dSin (\frac{2 π}{λ} (\bar{Z} - c \bar{t})), \\ \bar{b} (\bar{Z}, \bar{t}) = b_{0} + dSin (\frac{2 π}{λ} (\bar{Z} - c \bar{t})), \end{matrix}

(1)

The mathematical equations governing the flow of this incompressible, Newtonian fluid are provided as²⁷

\frac{\partial \bar{U}}{\partial \bar{X}} + \frac{\partial \bar{V}}{\partial \bar{Y}} + \frac{\partial \bar{W}}{\partial \bar{Z}} = 0,

(2)

ρ_{nf} (\frac{\partial \bar{U}}{\partial \bar{t}} + \bar{U} \frac{\partial \bar{U}}{\partial \bar{X}} + \bar{V} \frac{\partial \bar{U}}{\partial \bar{Y}} + \bar{W} \frac{\partial \bar{U}}{\partial \bar{Z}}) = - \frac{\partial \bar{P}}{\partial \bar{X}} + μ_{nf} (\frac{\partial^{2} \bar{U}}{\partial {\bar{X}}^{2}} + \frac{\partial^{2} \bar{U}}{\partial {\bar{Y}}^{2}} + \frac{\partial^{2} \bar{U}}{\partial {\bar{Z}}^{2}}),

(3)

ρ_{nf} (\frac{\partial \bar{V}}{\partial \bar{t}} + \bar{U} \frac{\partial \bar{V}}{\partial \bar{X}} + \bar{V} \frac{\partial \bar{V}}{\partial \bar{Y}} + \bar{W} \frac{\partial \bar{V}}{\partial \bar{Z}}) = - \frac{\partial \bar{P}}{\partial \bar{Y}} + μ_{nf} (\frac{\partial^{2} \bar{V}}{\partial {\bar{X}}^{2}} + \frac{\partial^{2} \bar{V}}{\partial {\bar{Y}}^{2}} + \frac{\partial^{2} \bar{V}}{\partial {\bar{Z}}^{2}}),

(4)

\begin{matrix} ρ_{nf} (\frac{\partial \bar{W}}{\partial \bar{t}} + \bar{U} \frac{\partial \bar{W}}{\partial \bar{X}} + \bar{V} \frac{\partial \bar{W}}{\partial \bar{Y}} + \bar{W} \frac{\partial \bar{W}}{\partial \bar{Z}}) = - \frac{\partial \bar{P}}{\partial \bar{Z}} \\ + μ_{nf} (\frac{\partial^{2} \bar{W}}{\partial {\bar{X}}^{2}} + \frac{\partial^{2} \bar{W}}{\partial {\bar{Y}}^{2}} + \frac{\partial^{2} \bar{W}}{\partial {\bar{Z}}^{2}}) + {(ρ β)}_{nf} g (\bar{T} - {\bar{T}}_{w}), \end{matrix}

(5)

{(ρ C_{p})}_{nf} (\frac{\partial \bar{T}}{\partial \bar{t}} + \bar{U} \frac{\partial \bar{T}}{\partial \bar{X}} + \bar{V} \frac{\partial \bar{T}}{\partial \bar{Y}} + \bar{W} \frac{\partial \bar{T}}{\partial \bar{Z}}) = k_{nf} (\frac{\partial^{2} \bar{T}}{\partial {\bar{X}}^{2}} + \frac{\partial^{2} \bar{T}}{\partial {\bar{Y}}^{2}} + \frac{\partial^{2} \bar{T}}{\partial {\bar{Z}}^{2}}) + Q_{0},

(6)

The dimensional form of the relevant boundary conditions is given as

\bar{W} = 0, \bar{T} = {\bar{T}}_{w}, for \frac{{\bar{x}}^{2}}{{\bar{a}}^{2}} + \frac{{\bar{y}}^{2}}{{\bar{b}}^{2}} = 1,

(7)

The two frames of reference (i.e. fixed frame and moving frame) are mathematically related by the following equations

\bar{x} = \bar{X}, \bar{y} = \bar{Y}, \bar{z} = \bar{Z} - c \bar{t}, \bar{p} = \bar{P}, \bar{u} = \bar{U}, \bar{v} = \bar{V}, \bar{w} = \bar{W} - c,

(8)

The useful dimensionless variables are provided by

\begin{matrix} x = \frac{\bar{x}}{D_{h}}, y = \frac{\bar{y}}{D_{h}}, z = \frac{\bar{z}}{λ}, t = \frac{c \bar{t}}{λ}, w = \frac{\bar{w}}{c}, p = \frac{D_{h}^{2} \bar{p}}{μ λ c}, θ = \frac{\bar{T} - {\bar{T}}_{w}}{{\bar{T}}_{b} - {\bar{T}}_{w}}, δ = \frac{b_{0}}{a_{0}}, ϕ = \frac{d}{b_{0}}, \\ u = \frac{λ \bar{u}}{D_{h} c}, v = \frac{λ \bar{v}}{D_{h} c}, a = \frac{\bar{a}}{D_{h}}, b = \frac{\bar{b}}{D_{h}}, G_{r} = \frac{g ({\bar{T}}_{b} - {\bar{T}}_{w}) D_{h}^{2} {(ρ β)}_{f}}{μ_{f} c}, γ = \frac{Q_{0} D_{h}^{2}}{({\bar{T}}_{b} - {\bar{T}}_{w}) k_{f}}, \end{matrix}

(9)

The hydraulic diameter of ellipse is written as

D_{h} = \frac{π b_{0}}{E (e)},

(10)

The eccentricity of ellipse is $e = \sqrt{1 - δ^{2}},$ and the second kind elliptical integral $E (e)$ is written as²⁸

Table 1 and 2 represents the numerical values and experimental formulas for thermo physical features of hybrid nanofluid respectively.

E (e) = \int_{0}^{\frac{π}{2}} \sqrt{1 - e^{2} S i n^{2} α} d α,

(11)

The equations (3)–(6) are obtained in their non-dimensional form after using equations (8) and (9), then final simplified form of non-dimensional equations is obtained by applying lubrication theory, (i.e. λ $\to \infty$ ).

\frac{\partial p}{\partial x} = 0,

(12)

\frac{\partial p}{\partial y} = 0,

(13)

\frac{dp}{dz} = \frac{μ_{nf}}{μ_{f}} (\frac{\partial^{2} w}{\partial x^{2}} + \frac{\partial^{2} w}{\partial y^{2}}) + \frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}} G_{r} θ,

(14)

\frac{k_{nf}}{k_{f}} (\frac{\partial^{2} θ}{\partial x^{2}} + \frac{\partial^{2} θ}{\partial y^{2}}) + γ = 0,

(15)

With non-dimensional boundary conditions

w = - 1, for \frac{x^{2}}{a^{2}} + \frac{y^{2}}{b^{2}} = 1,

(16)

θ = 0, for \frac{x^{2}}{a^{2}} + \frac{y^{2}}{b^{2}} = 1,

(17)

Where $a = \frac{E (e)}{π} [\frac{1}{δ} + ϕ Sin (2 π z)],$ and $b = \frac{E (e)}{π} [1 + ϕ Sin (2 π z)] .$

Table 1.

Thermo-physical properties of base fluid and CNT referred as Akbar et al.²⁹

Physical properties	Water $(H_{2} O)$	SWCNT
$C_{p}$	4179	425
$k$	0.613	6600
$ρ$	997.1	2600
$β \times 10^{- 5}$	21	1.5

Table 2.

Carbon nanotube model referred as Nadeem.³⁰

Density	$ρ_{nf} = (1 - ψ) ρ_{f} + ψ ρ_{CNT}$ ,
Heat capacity	${(ρ C_{p})}_{nf} = (1 - ψ) {(ρ C_{p})}_{f} + ψ {(ρ C_{p})}_{CNT}$ ,
Viscosity	$\frac{μ_{nf}}{μ_{f}} = \frac{1}{{(1 - ψ)}^{2.5}},$
Thermal conductivity	$\frac{k_{nf}}{k_{f}} = \frac{1 - ψ + 2 ψ (\frac{k_{CNT}}{k_{CNT} - k_{f}}) Log (\frac{k_{CNT} + k_{f}}{2 k_{f}})}{1 - ψ + 2 ψ (\frac{k_{f}}{k_{CNT} - k_{f}}) Log (\frac{k_{CNT} + k_{f}}{2 k_{f}})},$
Thermal expansion coefficient and density relation	${(ρ β)}_{nf} = (1 - ψ) {(ρ β)}_{f} + ψ {(ρ β)}_{CNT}$ .

Exact solution

Assume a polynomial form solution of temperature equation (15) given as

θ (x, y) = C_{1} x^{4} + C_{2} y^{4} + C_{3} x^{2} y^{2} + C_{4} x^{2} + C_{5} y^{2} + C_{6},

(18)

By using equation (18) in equation (15) and comparing coefficients of $x^{2}, y^{2}, x^{0}, y^{0}$ , we get these three equations given as

12 C_{1} + 2 C_{3} = 0,

(i)

2 C_{3} + 12 C_{2} = 0,

(ii)

2 C_{4} + 2 C_{5} = - \frac{γ}{\frac{k_{n f}}{k_{f}}},

(iii)

Now by using equation (18) in the relevant boundary condition for temperature given in equation (17) and by comparing coefficients of $x^{4}, x^{2}, x^{0}$ , we get these three equations

C_{1} a^{4} + C_{2} b^{4} - C_{3} a^{2} b^{2} = 0,

(iv)

- 2 C_{2} b^{4} + C_{3} a^{2} b^{2} + C_{4} a^{2} - C_{5} b^{2} = 0,

(v)

C_{2} b^{4} + C_{5} b^{2} + C_{6} = 0,

(vi)

By solving equations (i)–(vi), we get values of constants $C_{1}, C_{2}, C_{3}, C_{4}, C_{5}, C_{6}$ given in equation (19) as follows

\begin{matrix} C_{1} = 0, C_{2} = 0, C_{3} = 0, C_{4} = - \frac{b^{2} γ}{2 (a^{2} + b^{2}) (\frac{k_{nf}}{k_{f}})}, \\ C_{5} = - \frac{a^{2} γ}{2 (a^{2} + b^{2}) (\frac{k_{nf}}{k_{f}})}, C_{6} = \frac{a^{2} b^{2} γ}{2 (a^{2} + b^{2}) (\frac{k_{nf}}{k_{f}})}, \end{matrix}

(19)

Now utilizing values of these constants in equation (18), we get an exact solution for temperature profile that satisfies both equation and boundary conditions exactly and is given as

θ (x, y) = \frac{- a^{2} b^{2} (\frac{x^{2}}{a^{2}} + \frac{y^{2}}{b^{2}} - 1) γ}{2 (a^{2} + b^{2}) (\frac{k_{nf}}{k_{f}})},

(20)

The same procedure that is used to obtain temperature solution is used again to get an exact solution for velocity profile that satisfies both the momentum equation (14) and boundary condition given in equation (16). The exact solution for $w (x, y)$ is given as

\begin{matrix} w (x, y) = \frac{1}{24 {(a^{2} + b^{2})}^{2} (a^{4} + 6 a^{2} b^{2} + b^{4}) (\frac{k_{nf}}{k_{f}}) (\frac{μ_{nf}}{μ_{f}})} \\ [b^{8} {12 (\frac{k_{nf}}{k_{f}}) (\frac{dp}{dz} x^{2} - 2 (\frac{μ_{nf}}{μ_{f}})) + G_{r} x^{4} γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})}) \\ + 6 a^{2} b^{6} {2 (\frac{k_{nf}}{k_{f}}) (\frac{dp}{dz} (- b^{2} + 7 x^{2} + y^{2}) - 16 (\frac{μ_{nf}}{μ_{f}})) + G_{r} x^{2} (- b^{2} + x^{2} + y^{2}) γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})} \\ + a^{8} {- 12 (\frac{k_{nf}}{k_{f}}) (\frac{dp}{dz} (b^{2} - y^{2}) + 2 (\frac{μ_{nf}}{μ_{f}})) + G_{r} (5 b^{4} - 6 b^{2} y^{2} + y^{4}) γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})} \\ + a^{4} b^{4} {- 84 (\frac{k_{nf}}{k_{f}}) (b^{2} \frac{dp}{dz} - \frac{dp}{dz} (x^{2} + y^{2}) + 4 (\frac{μ_{nf}}{μ_{f}})) \\ + G_{r} (5 b^{4} + 5 x^{4} + 12 x^{2} y^{2} + 5 y^{4} - 2 b^{2} (16 x^{2} + 5 y^{2})) γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})} \\ + 2 a^{6} b^{2} {6 (\frac{k_{nf}}{k_{f}}) (\frac{dp}{dz} (- 7 b^{2} + x^{2} + 7 y^{2}) - 16 (\frac{μ_{nf}}{μ_{f}})) \\ + G_{r} (13 b^{4} + 3 y^{2} (x^{2} + y^{2}) - b^{2} (5 x^{2} + 16 y^{2})) γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})}] \end{matrix}

(21)

The mathematical result for dimensionless flow rate $q (z)$ is obtained by integrating the velocity solution $w (x, y)$ over the elliptic cross-section and given as

q (z) = \frac{- 3 ab (a^{2} + b^{2}) (\frac{k_{nf}}{k_{f}}) π (a^{2} b^{2} \frac{dp}{dz} + 4 (a^{2} + b^{2}) (\frac{μ_{nf}}{μ_{f}})) + a^{5} b^{5} G_{r} π γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})}{12 {(a^{2} + b^{2})}^{2} (\frac{k_{nf}}{k_{f}}) (\frac{μ_{nf}}{μ_{f}})}

(22)

The mathematical result of pressure gradient is evaluated from equation (22) as follows

\frac{dp}{dz} = - \frac{4 (a^{2} + b^{2}) (\frac{μ_{nf}}{μ_{f}}) [- \int_{0}^{1} abdz + ab π + Q - \frac{a^{5} b^{5} G_{r} π γ (\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}})}{12 {(a^{2} + b^{2})}^{2} (\frac{k_{nf}}{k_{f}}) (\frac{μ_{nf}}{μ_{f}})}]}{a^{3} b^{3} π},

(23)

The pressure rise expression over a single wavelength is calculated by

Δ P = \int_{0}^{1} \frac{\partial p}{\partial z} dz,

(24)

Result and discussion

This segment provides a graphical assessment of the mathematical results that are obtained in the previous section. We have plotted some 2-dimensional as well as 3-dimensional graphs for both velocity and temperature profiles. The velocity of fluid is plotted for increasing value of different parameter as displayed in Figures 2 to 5. All these velocity graphs reveal that the velocity gains its highest value in the centre of elliptic duct and it declines toward the walls of duct. Also, an axial symmetry flow behaviour is revealed by these velocity graphs. Thus, we can conclude that the peristaltic flow in this vertical elliptic duct is an axially symmetric flow. Figure 2(a) gives a 2-dimensional graphical result of velocity for incrementing value of $G_{r}$ and it is revealed that velocity gains higher value with increasing values of $G_{r}$ . Figure 2(b) provides a 3-dimensional graph of velocity for incrementing value of $G_{r}$ and as a result velocity increases. The increasing value of $G_{r}$ implies the stronger buoyancy forces and as a result the flow is accelerated toward the axial direction. There is increase in the velocity for increasing values of $γ$ , displayed in Figure 3(a). Figure 3(b) depicts a 3-dimensional graph for rising value of $γ$ and the velocity of fluid increases with this rising value of $γ$ . Figure 4(a) depicts a decline in the velocity of fluid with increasing $ψ$ . A 3-dimensional graph of velocity for increasing value of $ψ$ is presented in Figure 4(b). The velocity decreases as the volume fraction of carbon nanotubes increase in the base fluid. Since the velocity of fluid is dependent on its temperature in this study case and temperature decreases with increasing value of volume fraction of carbon nanotubes. Thus velocity also declines for increasing volumetric fraction of nanoparticles. There is increase in velocity for incrementing value of flow rate as revealed in Figure 5(a). Figure 5(b) shows a 3-dimensional graphical result of velocity for increasing value of $Q$ and velocity enhances with incrementing values of flow rate. The temperature of this heated fluid is plotted against distinct physical parameters as shown in Figures 6 and 7. These temperature graphs show that the temperature has maximum value in the centre but declines toward the walls of duct. Also an axially symmetric behaviour is seen in these temperature graphs. Figure 6(a) depicts an increase in the temperature of fluid for increasing value of $γ$ . A 3-dimensional graph of temperature for increasing value of $γ$ is presented in Figure 6(b). The temperature of fluid increases as the value of heat absorption parameter increases. Figure 7(a) reveals a decline in temperature as the value of $ψ$ increases. Figure 7(b) depicts a 3-dimensional graph of temperature profile for increasing values of $ψ$ . The temperature of fluid declines with increasing the volume fraction of carbon nanotubes in the base fluid. The high value of thermal conductivity of carbon nanotubes is the main reason behind quick decline in the temperature of fluid. In Figure 8(a) to (d), the pressure gradient is plotted for different physical parameters against the axial coordinate of elliptic duct. The value of $\frac{dp}{dz}$ increases for increasing values of $G_{r}$ , as shown in Figure 8(a). Figure 8(b) shows an increasing behaviour in value of $\frac{dp}{dz}$ for increasing values of $γ$ . The value of $\frac{dp}{dz}$ decreases with increasing volume fraction $ψ$ of carbon nanotubes, displayed in Figure 8(c). Figure 8(d) depicts a decline in the value of $\frac{dp}{dz}$ for increasing values of $Q$ . In Figure 9(a) to (c), the pressure rise graphical result $Δ P$ is plotted against the volumetric flow rate $Q$ . The value of $Δ P$ increases for increasing value of $G_{r}$ , depicted in Figure 9(a). Figure 9(b) also shows increase in $Δ P$ for increasing value of $γ$ . There is a decline in the value of $Δ P$ with increasing value of $ψ$ . The streamlines for this peristaltic flow problem in a vertical duct having elliptic cross-section are provided in Figure 10(a) to (d). The graphs are plotted for enhancing values of flow rate. The trapping slightly increases in size with enhancing value of $Q$ . Further, the sinusoidally moving walls can be clearly seen in the streamline graph.

Conclusions

The flow of a heated Newtonian viscous fluid with single wall carbon nanotubes inside a vertical elliptical duct with sinusoidally deformable walls is mathematically interpreted. The main advantage of using carbon nanotubes in the base fluid is enhancement in the thermal conductivity. The major findings that appear in this work are given below

The peristaltic flow in this vertical elliptic duct is an axially symmetric flow.

An increment in the values of $G_{r}$ implies stronger buoyancy forces and as a result the flow is accelerated in the axial direction.

An axially symmetric behaviour is seen in these temperature graphs.

The temperature of fluid decreases with increasing the volume fraction of carbon nanotubes in the base fluid. The high value of thermal conductivity of carbon nanotubes is the main reason behind quick decline in the temperature of fluid.

The trapping slightly increases in size with increasing values of flow rate $Q$ . Further, the sinusoidally moving walls can be clearly seen in these streamline graphs.

Footnotes

Appendix

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The author Salman Saleem extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through research groups program under Grant No. RGP.2/38/42

ORCID iDs

Sohail Nadeem

Alibek Issakhov

Author biographies

Salman Akhtar is a PhD research fellow at Quaid-i-Azam University, Islamabad, Pakistan. His field of research is applied mathematics and computational fluid mechanics. He has published many articles under the supervision of Prof. Sohail Nadeem.

Luthais B McCash is an Honorary Fellow at the University of Leicester. He was previously the Lead Data Scientist and Head of the Mathematical Modelling, Simulation & Computation at a Scottish Energy company. Although he is interested in Applied and Industrial Mathematics in general, Luthais’ current research focuses on quantum fluids, and applied PDEs. Currently he collaborates with colleagues in both Leicester and externally.

Sohail Nadeem is a Professor and chairman of department of Mathematics at Quaid-i-Azam University, Islamabad, Pakistan. His field of research is applied mathematics and computational fluid mechanics. His contribution to the field is recognized at national and international levels. He has produced a good number of PhD students. Also, he is serving as editor of various well reputed international journals.

Salman Saleem earned his PhD degree in mathematics and presently working as an Assistant professor at King Khalid University, Saudi Arabia. His area of research is computational fluid mechanics. He has published a good number of articles and serving as editor of Hindawi journals.

Alibek Issakhov is currently working Head of department at Alfarabi Kazkh National University, Kazakhstan. He is a scholar, researcher, and teacher per excellence. He has authored many research articles which show his excellence in the field.

References

Jaffrin

Shapiro

. Peristaltic pumping. Annu Rev Fluid Mech 1971; 3(1): 13–37.

Pozrikidis

. A study of peristaltic flow. J Fluid Mech 1987; 180: 515–527.

Yin

Fung

. Peristaltic waves in circular cylindrical tubes. J Appl Mech 1969; 36(3): 579–587.

Srinivas

Kothandapani

. Peristaltic transport in an asymmetric channel with heat transfer—a note. Int Commun Heat Mass Transfer 2008; 35(4): 514–522.

Nadeem

Akram

. Peristaltic flow of a Jeffrey fluid in a rectangular duct. Nonlinear Anal Real World Appl 2010; 11(5): 4238–4247.

Ali

Saleem

Mumraiz

Saleem

Awais

Khan Marwat

. Investigation on TiO2-Cu/H2O hybrid nanofluid with slip conditions in MHD peristaltic flow of Jeffery material. J Ther Anal Calori, 2021; 143: 1985–1996.

Saleem

Akhtar

Nadeem

Saleem

Ghalambaz

Issakhov

. Mathematical study of Electroosmotically driven peristaltic flow of Casson fluid inside a tube having systematically contracting and relaxing sinusoidal heated walls. Chin J Phys 2021; 71: 300–311.

Akbar

Nadeem

Khan

. Numerical simulation of peristaltic flow of a Carreau nanofluid in an asymmetric channel. Alex Eng J 2014; 53(1): 191–197.

Nadeem

Akram

. Magnetohydrodynamic peristaltic flow of a hyperbolic tangent fluid in a vertical asymmetric channel with heat transfer. Acta Mech Sin 2011; 27(2): 237–250.

10.

Huang

. Study of the fluid flow in the elliptical duct by the method of characteristics. J Press Vessel Technol 1993; 115(1): 80–84.

11.

Richardson

. Leveque solution for flow in an elliptical duct. Lett Heat Mass Transfer 1980; 7(5): 353–362.

12.

Abdel-Wahed

Attia

Hifni

. Experiments on laminar flow and heat transfer in an elliptical duct. Int J Heat Mass Transfer 1984; 27(12): 2397–2413.

13.

Sakalis

Hatzikonstantinou

Kafousias

. Thermally developing flow in elliptic ducts with axially variable wall temperature distribution. Int J Heat Mass Transfer 2002; 45(1): 25–35.

14.

Shariat

Akbarinia

Nezhad

, et al. Numerical study of two phase laminar mixed convection nanofluid in elliptic ducts. Appl Therm Eng 2011; 31(14–15): 2348–2359.

15.

Ragueb

Mansouri

. A numerical study of viscous dissipation effect on non-Newtonian fluid flow inside elliptical duct. Energy Convers Manag 2013; 68: 124–132.

16.

Choi

Eastman

. Enhancing thermal conductivity of fluids with nanoparticles. No. ANL/MSD/CP-84938; CONF-951135-29, Argonne National Lab, Argonne, IL, United States, 1995.

17.

Akbar

Nadeem

Hayat

, et al. Peristaltic flow of a nanofluid in a non-uniform tube. Heat Mass Transfer 2012; 48(3): 451–459.

18.

Shahzadi

Sadaf

Nadeem

, et al. Bio-mathematical analysis for the peristaltic flow of single wall carbon nanotubes under the impact of variable viscosity and wall properties. Comput Methods Programs Biomed 2017; 139: 137–147.

19.

Akbar

. MHD peristaltic flow with carbon nanotubes in an asymmetric channel. J Comput Theor Nanosci 2014; 11(5): 1323–1329.

20.

Ellahi

Riaz

Nadeem

. A theoretical study of Prandtl nanofluid in a rectangular duct through peristaltic transport. Appl Nanosci 2014; 4(6): 753–760.

21.

Huminic

. The influence of hybrid nanofluids on the performances of elliptical tube: recent research and numerical study. Int J Heat Mass Transfer 2019; 129: 132–143.

22.

Shariat

Moghari

Sajjadi

, et al. Numerical investigation of Al₂O₃/water nanofluid in horizontal elliptic ducts using two phase mixture model. J Comput Theor Nanosci 2013; 10(1): 199–207.

23.

Ellahi

Zeeshan

Waheed

, et al. Natural convection nanofluid flow with heat transfer analysis of carbon nanotubes–water nanofluid inside a vertical truncated wavy cone. Math Methods Appl Sci. Epub ahead of print 22 February 2021. DOI: 10.1002/mma.7281.

24.

Maskeen

Zeeshan

Mehmood

, et al. Heat transfer enhancement in hydromagnetic alumina–copper/water hybrid nanofluid flow over a stretching cylinder. J Therm Anal Calorim 2019; 138(2): 1127–1136.

25.

Zeeshan

Ellahi

Mabood

, et al. Numerical study on bi-phase coupled stress fluid in the presence of Hafnium and metallic nanoparticles over an inclined plane. Int J Numer Methods Heat Fluid Flow 2019; 29(8): 2854–2869.

26.

Zeeshan

Bhatti

Ijaz

, et al. Biologically-inspired transport of solid spherical nanoparticles in an electrically-conducting viscoelastic fluid with heat transfer. Therm Sci 2020; 24(2(B)): 1251–1260.

27.

Nadeem

Akhtar

Saleem

. Peristaltic flow of a heated Jeffrey fluid inside an elliptic duct: streamline analysis. Appl Math Mech 2021; 42(4): 583–592.

28.

Yang

Chu

Zhang

. Monotonicity of the ratio for the complete elliptic integral and Stolarsky mean. J Inequal Appl 2016; 2016(1): 1–10.

29.

Akbar

Raza

Ellahi

. Influence of induced magnetic field and heat flux with the suspension of carbon nanotubes for the peristaltic flow in a permeable channel. J Magn Magn Mater 2015; 381: 405–415.

30.

Nadeem

. Single wall carbon nanotube (SWCNT) analysis on peristaltic flow in an inclined tube with permeable walls. Int J Heat Mass Transfer 2016; 97: 794–802.