删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Effectiveness of heat and mass transfer on mixed convective peristaltic motion of nanofluid with irr

本站小编 Free考研考试/2022-01-02

Y Akbar1, F M Abbasi1, U M Zahid1, S A Shehzad,21Department of Mathematics, COMSATS University Islamabad, Islamabad, Pakistan
2Department of Mathematics, COMSATS University Islamabad, Sahiwal, Pakistan

Received:2021-03-25Revised:2021-07-2Accepted:2021-07-2Online:2021-08-31


Abstract
Entropy generation is one of the key features in analysis as it exhibits irreversibility of the system. Therefore, the present study investigates the entropy generation rate in a mixed convective peristaltic motion of a reactive nanofluid through an asymmetrical divergent channel with heat and mass transfer characteristics. The endorsed nanofluid model holds thermophoresis and Brownian diffusions. Mathematical modeling is configured under the effects of mixed convection, heat generation/absorption and viscous dissipation. A chemical reaction is also introduced for the description of mass transportation. The resulting system of differential equations is numerically tackled by employing the Shooting method. The findings reveal that entropy generation rises by improving the Brownian motion and thermophoresis parameters. The temperature of the nanofluid decreases due to rising buoyancy forces caused by the concentration gradient. The concentration profile increases by increasing the chemical reaction parameter. The velocity increases by enhancing the Brownian motion parameter.
Keywords: peristaltic flow;thermophoresis and Brownian motion;mixed convection;chemical reaction;irreversibility;nanofluid


PDF (1502KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Y Akbar, F M Abbasi, U M Zahid, S A Shehzad. Effectiveness of heat and mass transfer on mixed convective peristaltic motion of nanofluid with irreversibility rate. Communications in Theoretical Physics, 2021, 73(10): 105004- doi:10.1088/1572-9494/ac10bf

1. Introduction

The mechanism of nanofluids is becoming a hot topic for research owing to their heat conductive properties. Typically, a nanofluid is a dilute aggregation of fibers and particles of nanometer scale distributed in a carrier liquid. Nanomaterials are composed of metals, their oxides and nanotubes, whereas base liquids include engine oils, lubricants, grease, water and ethylene glycol etc. Some nanofluids possess idiosyncratic properties. They have had a broad array of applications in many biomedical, industrial and engineering systems. Some of them include domestic refrigerators, fission and fusion reactions, fuel cells, cancer treatments, chemotherapy, hybrid-powered engines, milk sterilization, and electronic cooling systems. In 1995, innovative ideas of nanofluids were reported by Choi [1]. Later, Buongiorno [2] stated that nanofluids are the colloidal suspensions of small particles in liquids. In this study, the author also designed a model to assess the flow of nanoliquids using Brownian motion and thermophoresis effects. Prasad et al [3] reported a new concept to deliver nanomedicines. The main aim of this research was to treat rheumatoid arthritis. Kleinstreuer and Feng [4] conducted an experimental and theoretical study on optimizing nanofluids' thermal conductivity. A water-based nanoliquid moving over wavy surfaces in a porous space of spherical packing beds was highlighted by Hassan et al [5]. Vallejo et al [6] studied Ag–graphene-based hybrid nanoliquids for direct solar absorption. Recently, Shahrestani et al [7] numerically reviewed the fluid flow and thermal analysis of ${\mathrm{Al}}_{2}{{\rm{O}}}_{3}/{{\rm{H}}}_{2}{\rm{O}}$ nanoliquid inside a microchannel of inner diameter 0.1 mm.

Peristalsis is a fluid flow mechanism that occurs due to the successive series of muscle relaxation and contraction by which fluids are transferred from one place to another in a physiological system. Peristaltic flows are encountered in cilia motion, movement of spermatozoa, heart–lung and dialysis machines, medical infusion pumps and in rocket chamber fuel control technologies. In 1966, Latham [8] proposed a mechanism of liquid motion via peristaltic pump. Shapiro et al [9] made a considerable development on the mathematical modeling of peristaltic motion [8] and explored this mechanism by applying ‘long wavelength and small Reynolds number assumptions.' This research employed both theoretical and numerical methods. Akbar and Nadeem [10] discussed the peristaltic movement of a nanoliquid with endoscopic effects. The process of peristalsis has been studied by many researchers both theoretically and experimentally [1119]. Riaz et al [20] explored the peristaltic propulsion of a Jeffrey nanoliquid through a porous rectangular duct.

The aspects of Brownian motion and thermophoresis are very important in the fields of science and technology. These effects are of great importance in the production of germanium oxide and silicon, etc. Beg and Tripathi [21] presented a theoretical analysis to explore the phenomenon of peristalsis with double-diffusive convection in nanoliquids via a deformable channel. For this review, the authors utilized Buongiorno's model to investigate the behavior of a nanoliquid. Peristaltic transportation of a magnetohydrodynamic (MHD) nanoliquid under the influence of Brownian motion and thermophoresis were highlighted by Hayat et al [22]. A numerical study on peristaltic motion of nanoliquids between the gap of two coaxial tubes with distinct structures and configurations was conducted by Ellahi et al [23]. Khan et al [24] analyzed the movement of non-Newtonian material subject to Dufour and Soret impacts over a porous space. The effects of wall flexibility with convective conditions on MHD peristaltic motion of an Eyring–Powell nanoliquid were highlighted by Nisar et al [25].

In thermal design and heat transfer, the analysis of the second law of thermodynamics is a widely used and useful tool. The study of entropy production is indeed one of the main aspects in analysis because it illustrates the disorder of the system. When any engineering system operates, it encounters irreversible energy losses that affect the efficiency of the machine. Therefore, circumstances are needed that minimize the production of entropy. The innovative idea of entropy generation (second law) was explained by Bejan [26]. In this analysis, the author outlined the key ways to minimize entropy at the system level of the structure. Nadeem and Sadia [27] studied entropy generation for a Johnson–Segalman fluid by assuming variable viscosity. Akbar et al [28] studied entropy generation for peristalsis of $\mathrm{Cu}+{{\rm{H}}}_{2}{\rm{O}}$ nanoliquid. Abbas et al [29] studied irreversible analysis for peristaltic movement of nanoliquids flowing in a two-dimensional (2D) channel with compliant walls. Farooq et al [30] explored the role of entropy production in peristalsis of multi- and single-walled carbon nanotubes. Consideration of entropy production for peristaltic transportation of $\mathrm{Cu}-{{\rm{H}}}_{2}{\rm{O}}$ nanoliquid with radiative heat flux was provided by Akbar et al [31]. MHD peristaltic transportation of a Sutterby nanoliquid with Hall current and mixed convection was highlighted by Hayat et al [32]. Some important studies related to irreversibility analysis for peristaltic transport of a nanoliquid can also be viewed in references [3339].

In view of these investigations, the main aim of the present study is to address the irreversibility rate, heat and mass transfer for peristaltic flow of a reactive nanofluid flowing through a divergent channel. The effects of thermophoresis, chemical reaction and Brownian motion are taken into consideration. Lubrication theory is used in the mathematical formulation. The Shooting method is employed to obtain the solution of differential equations. Entropy generation, Bejan number, mass and heat transfers, and skin friction are addressed under the influence of different pertinent parameters. Further, comparisons of heat, mass transfer rate and skin friction at the channel wall in light of two different numeric methods are also provided. The present study finds its application in transport phenomena in chemical engineering and energy systems to increase the heat capacity of machines.

2. Problem formulation

A 2D mixed convective peristaltic motion of reactive nanofluid through a divergent channel is considered. The wave trains are propagating with velocity c on the channel walls. Further, the channel is assumed to be asymmetric. A physical model of the problem is shown in figure 1. Mathematically, deformation in the channel walls is defined as [16]:$\begin{eqnarray}\overline{{H}_{2}}(\overline{X},\overline{t})=-{d}_{2}-\bar{m}\bar{X}-{a}_{2}\sin \left(\displaystyle \frac{2\pi }{\lambda }(\overline{X}-c\overline{t})\right),\end{eqnarray}$$\begin{eqnarray}\overline{{H}_{1}}(\overline{X},\overline{t})={d}_{1}+\bar{m}\bar{X}+{a}_{1}\sin \left(\displaystyle \frac{2\pi }{\lambda }(\overline{X}-c\overline{t})+\gamma \right),\end{eqnarray}$where ${d}_{1}+{d}_{2}$, a1, a2, c, t, m, γ and λ are the channel width, amplitudes, wave speed, time, nonuniform parameter, phase difference and wavelength, respectively. The governing equations for the present model using chemical reaction, thermophoresis, mixed convection, viscous dissipation and Brownian motion are given as [24, 33, 38]:$\begin{eqnarray}\displaystyle \frac{\partial \overline{U}}{\partial \overline{X}}+\displaystyle \frac{\partial \overline{V}}{\partial \overline{Y}}=0,\end{eqnarray}$$\begin{eqnarray}\begin{array}{l}{\rho }_{f}\left(\displaystyle \frac{\partial }{\partial \overline{t}}+\overline{U}\displaystyle \frac{\partial }{\partial \overline{X}}+\overline{V}\displaystyle \frac{\partial }{\partial \overline{Y}}\right)\overline{U}=-\displaystyle \frac{\partial \overline{P}}{\partial \overline{X}}\\ \quad +\,{\mu }_{f}\left(\displaystyle \frac{{\partial }^{2}}{\partial {\overline{X}}^{2}}+\displaystyle \frac{{\partial }^{2}}{\partial {\overline{Y}}^{2}}\right)\overline{U}+\alpha {\rho }_{f}g(T-{T}_{m})\\ \quad +\,{\rho }_{f}g{\alpha }^{* }(C-{C}_{m}),\end{array}\end{eqnarray}$$\begin{eqnarray}\begin{array}{l}{\rho }_{f}\left(\displaystyle \frac{\partial }{\partial \overline{t}}+\overline{U}\displaystyle \frac{\partial }{\partial \overline{X}}+\overline{V}\displaystyle \frac{\partial }{\partial \overline{Y}}\right)\overline{V}\\ \quad =-\displaystyle \frac{\partial \overline{P}}{\partial \overline{Y}}+{\mu }_{f}\left(\displaystyle \frac{{\partial }^{2}}{\partial {\overline{X}}^{2}}+\displaystyle \frac{{\partial }^{2}}{\partial {\overline{Y}}^{2}}\right)\overline{V},\end{array}\end{eqnarray}$$\begin{eqnarray}\begin{array}{l}{\rho }_{f}{C}_{f}\left(\displaystyle \frac{\partial }{\partial \overline{t}}+\overline{U}\displaystyle \frac{\partial }{\partial \overline{X}}+\overline{V}\displaystyle \frac{\partial }{\partial \overline{Y}}\right)T\\ =\,{K}_{f}\left(\displaystyle \frac{{\partial }^{2}T}{\partial {\overline{Y}}^{2}}+\displaystyle \frac{{\partial }^{2}T}{\partial {\overline{X}}^{2}}\right)+{\rm{\Phi }}+{\mu }_{f}\left[2\left({\left(\displaystyle \frac{\partial \overline{V}}{\partial \overline{Y}}\right)}^{2}\right.\right.\\ \left.\left.+\,{\left(\displaystyle \frac{\partial \overline{U}}{\partial \overline{X}}\right)}^{2}\right)+{\left(\displaystyle \frac{\partial \overline{V}}{\partial \overline{X}}+\displaystyle \frac{\partial \overline{U}}{\partial \overline{Y}}\right)}^{2}\right]+\tau {\rho }_{f}{C}_{f}\left[{D}_{B}\right.\\ \times \,\left(\displaystyle \frac{\partial C}{\partial \overline{Y}}\displaystyle \frac{\partial T}{\partial \overline{Y}}+\displaystyle \frac{\partial C}{\partial \overline{X}}\displaystyle \frac{\partial T}{\partial \overline{X}}\right)\\ \left.+\,\displaystyle \frac{{D}_{T}}{{T}_{m}}{\left(\displaystyle \frac{\partial T}{\partial \overline{X}}+\displaystyle \frac{\partial T}{\partial \overline{Y}}\right)}^{2}\right],\end{array}\end{eqnarray}$$\begin{eqnarray}\begin{array}{l}\left(\displaystyle \frac{\partial }{\partial \overline{t}}+\overline{U}\displaystyle \frac{\partial }{\partial \overline{X}}+\overline{V}\displaystyle \frac{\partial }{\partial \overline{Y}}\right)\overline{C}=\displaystyle \frac{{D}_{T}}{{T}_{m}}\left(\displaystyle \frac{{\partial }^{2}T}{\partial {\overline{Y}}^{2}}+\displaystyle \frac{{\partial }^{2}T}{\partial {\overline{X}}^{2}}\right)\\ +\,{D}_{B}\left(\displaystyle \frac{{\partial }^{2}C}{\partial {\overline{Y}}^{2}}+\displaystyle \frac{{\partial }^{2}C}{\partial {\overline{X}}^{2}}\right)-\hat{{\gamma }_{c}}(C-{C}_{m}).\end{array}\end{eqnarray}$In the above equations P, Cf, $\hat{{\gamma }_{c}}$, DB, DT, C, α and ${\alpha }^{* }$ represent pressure, specific heat of fluid, chemical reaction rate, mass diffusivity, thermal diffusion, concentration, and thermal and concentration expansion coefficients, respectively.

Figure 1.

New window|Download| PPT slide
Figure 1.Geometry of the problem.


Invoking the following dimensionless quantities:$\begin{eqnarray}\begin{array}{rcl}x & = & \displaystyle \frac{\overline{X}}{\lambda },\,\,u=\displaystyle \frac{\overline{U}}{c},\,\,y=\displaystyle \frac{\overline{Y}}{{d}_{1}},\,\,v=\displaystyle \frac{\overline{V}}{c\delta },\\ \delta & = & \displaystyle \frac{{d}_{1}}{\lambda },\,\,\,\,{h}_{1}=\displaystyle \frac{\overline{{H}_{1}}}{{d}_{1}},\,\,d=\displaystyle \frac{{d}_{2}}{{d}_{1}},\,\,\,{h}_{2}=\displaystyle \frac{\overline{{H}_{2}}}{{d}_{1}},\\ b & = & \displaystyle \frac{{a}_{2}}{{d}_{1}},\,\,\theta =\displaystyle \frac{T-{T}_{m}}{{T}_{1}-{T}_{0}},\,\,p=\displaystyle \frac{{d}_{1}^{2}\overline{p}}{c\lambda {\mu }_{f}},\\ E & = & \displaystyle \frac{{c}^{2}}{{C}_{f}({T}_{1}-{T}_{0})},\,\,{\Pr }=\displaystyle \frac{{\mu }_{f}{C}_{f}}{{K}_{f}},\,\,{Br}={\Pr }\,\cdot \,E,\\ u & = & {\psi }_{y},\,\,\,{Re}\,=\,\displaystyle \frac{\rho {cd}}{{\mu }_{f}},\,\,\,v=-{\psi }_{x},\\ {N}_{b} & = & \tau \displaystyle \frac{{D}_{B}({C}_{1}-{C}_{0})}{\nu },\,\,\,{Sr}=\displaystyle \frac{\rho {D}_{B}{K}_{T}(T-{T}_{0})}{{\mu }_{f}{T}_{m}(C-{C}_{0})},\\ {N}_{t} & = & \tau \displaystyle \frac{{D}_{T}({T}_{1}-{T}_{0})}{\nu {T}_{m}},\,\phi =\displaystyle \frac{C-{C}_{m}}{{C}_{1}-{C}_{0}},\,{Sc}=\displaystyle \frac{{\mu }_{f}}{\rho {D}_{B}},\\ {Sr} & = & \displaystyle \frac{{D}_{B}{K}_{T}(C-{C}_{0})}{{C}_{f}{C}_{s}{\mu }_{f}(T-{T}_{0})},{\gamma }_{c}=\displaystyle \frac{\hat{{\gamma }_{c}}{d}_{1}^{2}}{\nu },\\ {m}_{1} & = & \displaystyle \frac{\lambda \overline{m}}{{d}_{1}},\,\,\,\,{G}_{t}=\displaystyle \frac{{\rho }_{f}g\alpha ({T}_{1}-{T}_{0}){d}_{1}^{2}}{{\mu }_{f}c},\\ a & = & \displaystyle \frac{{a}_{1}}{{d}_{1}},\quad {G}_{c}=\displaystyle \frac{{\rho }_{f}g{\alpha }^{* }({C}_{1}-{C}_{0}){d}_{1}^{2}}{{\mu }_{f}c},\end{array}\end{eqnarray}$and applying ‘low Reynolds number and long wavelength' approximations, equations (4)–(7) take the following form:$\begin{eqnarray}0=-\displaystyle \frac{\partial p}{\partial x}+\displaystyle \frac{{\partial }^{3}\psi }{\partial {y}^{3}}+{G}_{t}\theta +{G}_{c}\phi ,\end{eqnarray}$$\begin{eqnarray}0=-\displaystyle \frac{\partial p}{\partial y},\end{eqnarray}$$\begin{eqnarray}0=\displaystyle \frac{{\partial }^{2}\theta }{\partial {y}^{2}}+{Br}{\left(\displaystyle \frac{{\partial }^{2}\psi }{\partial {y}^{2}}\right)}^{2}+{N}_{b}{\Pr }\displaystyle \frac{\partial \phi }{\partial y}\displaystyle \frac{\partial \theta }{\partial y}+{N}_{t}{\Pr }{\left(\displaystyle \frac{\partial \theta }{\partial y}\right)}^{2}+\varepsilon ,\end{eqnarray}$$\begin{eqnarray}0=\displaystyle \frac{{\partial }^{2}\phi }{\partial {y}^{2}}+\displaystyle \frac{{N}_{t}}{{N}_{b}}\displaystyle \frac{{\partial }^{2}\theta }{\partial {y}^{2}}-{Sc}{\gamma }_{c}\phi .\end{eqnarray}$Here ψ, Gt, Br, Gc, Sc, Nb, ${\gamma }_{c}$, Nt, and Pr are the dimensionless stream function, temperature Grashof number, Brinkman number, concentration Grashof number, Schmidt number, Brownian motion parameter, chemical reaction parameter, thermophoresis parameter and Prandtl number, respectively. Determining the nondimensional flow rates η and F, defined by [16]:$\begin{eqnarray}F=\eta +b\sin 2\pi (x-t)+a\sin 2\pi [(x-t)+\gamma ],\end{eqnarray}$the dimensionless boundary conditions are given as [16]:$\begin{eqnarray}\begin{array}{rcl}\psi & = & \displaystyle \frac{F}{2},\quad \displaystyle \frac{\partial \psi }{\partial y}=0,\quad \theta =-\displaystyle \frac{1}{2},\\ \phi & = & -\displaystyle \frac{1}{2},\quad {at}\,\,y={h}_{1},\\ \psi & = & -\displaystyle \frac{F}{2},\quad \displaystyle \frac{\partial \psi }{\partial y}=0,\quad \theta =\displaystyle \frac{1}{2},\\ \phi & = & \displaystyle \frac{1}{2},\quad {at}\,\,y={h}_{2},\quad \end{array}\end{eqnarray}$where$\begin{eqnarray}\begin{array}{rcl}{h}_{1} & = & 1+{m}_{1}x+a\sin 2\pi [(x-t)+\gamma ],\\ {h}_{2} & = & -d-{m}_{1}x-b\sin 2\pi (x-t).\end{array}\end{eqnarray}$

3. Entropy generation

Mathematically, the dimensional entropy generation rate is defined as [33, 38]:$\begin{eqnarray}\begin{array}{l}{S}_{G}=\displaystyle \frac{{K}_{f}}{{\left({T}_{1}-{T}_{0}\right)}^{2}}\left[{\left(\displaystyle \frac{\partial T}{\partial X}\right)}^{2}+{\left(\displaystyle \frac{\partial T}{\partial Y}\right)}^{2}\right]+\displaystyle \frac{\mu }{({T}_{1}-{T}_{0})}\\ \left[2\left({\left(\displaystyle \frac{\partial \overline{V}}{\partial \overline{Y}}\right)}^{2}+{\left(\displaystyle \frac{\partial \overline{U}}{\partial \overline{X}}\right)}^{2}\right)+{\left(\displaystyle \frac{\partial \overline{V}}{\partial \overline{X}}+\displaystyle \frac{\partial \overline{U}}{\partial \overline{Y}}\right)}^{2}\right]\\ +\,\displaystyle \frac{\tau {\rho }_{f}{C}_{f}}{({T}_{1}-{T}_{0})}\left[{D}_{B}\left(\displaystyle \frac{\partial T}{\partial \overline{Y}}\displaystyle \frac{\partial C}{\partial \overline{Y}}+\displaystyle \frac{\partial C}{\partial \overline{X}}\displaystyle \frac{\partial T}{\partial \overline{X}}\right)\right.\\ \left.+\,\displaystyle \frac{{D}_{T}}{{T}_{m}}{\left(\displaystyle \frac{\partial T}{\partial \overline{X}}+\displaystyle \frac{\partial T}{\partial \overline{Y}}\right)}^{2}\right]+\displaystyle \frac{{\rm{\Phi }}}{({T}_{1}-{T}_{0})}.\end{array}\end{eqnarray}$In dimensionless form, the entropy generation number $({N}_{s})$ is presented as [33]:$\begin{eqnarray}\begin{array}{rcl}{N}_{s}=\displaystyle \frac{{S}_{G}}{{S}_{c}} & = & {\left(\displaystyle \frac{\partial \theta }{\partial y}\right)}^{2}+{Br}{\left(\displaystyle \frac{{\partial }^{2}\psi }{\partial {y}^{2}}\right)}^{2}\\ & & +{{PrN}}_{b}\displaystyle \frac{\partial \phi }{\partial y}\displaystyle \frac{\partial \theta }{\partial y}+{{PrN}}_{t}{\left(\displaystyle \frac{\partial \theta }{\partial y}\right)}^{2}+\varepsilon ,\end{array}\end{eqnarray}$where ${S}_{c}=\tfrac{{K}_{f}}{{d}_{1}^{2}}$.

The Bejan number $({B}_{e})$ is the ratio of entropy produced by conduction to the total entropy. The mathematical interpretation of the Bejan number $({B}_{e})$ is expressed as follows [33]:$\begin{eqnarray}{B}_{e}=\displaystyle \frac{{\left(\tfrac{\partial \theta }{\partial y}\right)}^{2}}{{\left(\tfrac{\partial \theta }{\partial y}\right)}^{2}+{Br}{\left(\tfrac{{\partial }^{2}\psi }{\partial {y}^{2}}\right)}^{2}+{{PrN}}_{b}\tfrac{\partial \phi }{\partial y}\tfrac{\partial \theta }{\partial y}+{{PrN}}_{t}{\left(\tfrac{\partial \theta }{\partial y}\right)}^{2}+\varepsilon }.\end{eqnarray}$

4. Numerical solution

For numerical solution we intend to use the Shooting technique with RK-4 on differential equations (9)–(12) subject to boundary conditions (14). In this connection the boundary conditions (14) for the system of differential equations (9)–(12) are simplified into initial conditions. Therefore, the prescribed problem to be addressed is as follows:$\begin{eqnarray}0=\displaystyle \frac{{\partial }^{4}\psi }{\partial {y}^{4}}+{G}_{t}\displaystyle \frac{\partial \theta }{\partial y}+{G}_{c}\displaystyle \frac{\partial \phi }{\partial y},\end{eqnarray}$$\begin{eqnarray}0=\displaystyle \frac{{\partial }^{2}\theta }{\partial {y}^{2}}+{Br}{\left(\displaystyle \frac{{\partial }^{2}\psi }{\partial {y}^{2}}\right)}^{2}+{N}_{b}{\Pr }\displaystyle \frac{\partial \phi }{\partial y}\displaystyle \frac{\partial \theta }{\partial y}+{N}_{t}{\Pr }{\left(\displaystyle \frac{\partial \theta }{\partial y}\right)}^{2}+\varepsilon ,\end{eqnarray}$$\begin{eqnarray}0=\displaystyle \frac{{\partial }^{2}\phi }{\partial {y}^{2}}+\displaystyle \frac{{N}_{t}}{{N}_{b}}\displaystyle \frac{{\partial }^{2}\theta }{\partial {y}^{2}}-{Sc}{\gamma }_{c}\phi .\end{eqnarray}$The initial conditions are written as:$\begin{eqnarray*}\psi =\displaystyle \frac{F}{2},\,\displaystyle \frac{\partial \psi }{\partial y}=0,\,\displaystyle \frac{\partial \theta }{\partial y}={\xi }_{1},\,\displaystyle \frac{\partial \phi }{\partial y}={\xi }_{2},\,\mathrm{at}\,y={h}_{2},\end{eqnarray*}$$\begin{eqnarray}\displaystyle \frac{{\partial }^{2}\psi }{\partial {y}^{2}}={\xi }_{3},\,\displaystyle \frac{{\partial }^{3}\psi }{\partial {y}^{3}}={\xi }_{4},\,\theta =\displaystyle \frac{1}{2},\,\phi =\displaystyle \frac{1}{2},\,\mathrm{at}\,y={h}_{2}.\end{eqnarray}$We partially differentiate equations (19)–(21) with respect to ${\xi }_{1}$, ${\xi }_{2}$, ${\xi }_{3}$ and ${\xi }_{4}$ separately by letting:$\begin{eqnarray*}\displaystyle \frac{\partial \psi }{\partial {\xi }_{1}}={\psi }_{1},\,\displaystyle \frac{\partial \theta }{\partial {\xi }_{1}}={\theta }_{1},\,\displaystyle \frac{\partial \phi }{\partial {\xi }_{1}}={\phi }_{1},\end{eqnarray*}$$\begin{eqnarray*}\displaystyle \frac{\partial \psi }{\partial {\xi }_{2}}={\psi }_{2},\,\displaystyle \frac{\partial \theta }{\partial {\xi }_{2}}={\theta }_{2},\,\displaystyle \frac{\partial \phi }{\partial {\xi }_{2}}={\phi }_{2},\end{eqnarray*}$$\begin{eqnarray*}\displaystyle \frac{\partial \psi }{\partial {\xi }_{3}}={\psi }_{3},\,\displaystyle \frac{\partial \theta }{\partial {\xi }_{3}}={\theta }_{3},\,\displaystyle \frac{\partial \phi }{\partial {\xi }_{3}}={\phi }_{3},\end{eqnarray*}$$\begin{eqnarray}\displaystyle \frac{\partial \psi }{\partial {\xi }_{4}}={\psi }_{4},\,\displaystyle \frac{\partial \theta }{\partial {\xi }_{4}}={\theta }_{4},\,\displaystyle \frac{\partial \phi }{\partial {\xi }_{4}}={\phi }_{4}.\end{eqnarray}$Hence, the stated problem is converted into four systems of equations. These equations are further simplified to a system of first-order differential equations. The system of differential equations are then solved numerically by using the Shooting technique with the RK-4 method. Appropriate values of unknown initial conditions ${\xi }_{1}$, ${\xi }_{2}$, ${\xi }_{3}$ and ${\xi }_{4}$ are approximated and optimized through Newton's technique until the boundary conditions are met. The step scale is considered as ${\rm{\Delta }}h=0.05$ and convergence criterion is taken as 10−6. All necessary calculations are performed in Mathematica software.

5. Results and discussion

This section provides the graphs of different flow quantities for several parameters of interest. Plots for heat transfer, irreversible rate, mass transfer and velocity profile are presented and analyzed.

5.1. Heat transfer analysis

Analysis of temperature profile is presented through figures 25. The maximum value of temperature profile is observed near the center of the channel. Figure 2 indicates that the temperature profile decreases by increasing buoyancy forces caused by concentration difference (Gc). The displayed result is due to an increment in viscosity which leads to a drop in temperature. However, figure 3 shows that temperature increases by improving the temperature Grashof number (Gt). Physically, it refers to the condition where improved buoyancy forces are strengthened due to the temperature gradient. Figure 4 shows that temperature increases for a rise in Nb. The reason is that kinetic energy is augmented as Nb is improved which incorporates it into internal energy and thus raises the temperature within the channel. An increase in temperature is obtained by increasing the amount of the thermophoresis parameter (Nt) (see figure 5). For higher Nt, nanoparticles migrate from warm regions to cool regions and hence the temperature increases. Here, the numerical outcomes presented via figures 4 and 5 are found to be well matched with the results reported by Hayat et al [32] and Ellahi et al [23].

Figure 2.

New window|Download| PPT slide
Figure 2.Variation in θ with Gc.


Figure 3.

New window|Download| PPT slide
Figure 3.Variation in θ with Gt.


Figure 4.

New window|Download| PPT slide
Figure 4.Variation in θ with Nb.


Figure 5.

New window|Download| PPT slide
Figure 5.Variation in θ with Nt.


Numerical values of the heat transfer rate at the channel wall ($-\theta ^{\prime} ({h}_{1})$) for variations in Gt, Gc, Nt, and Nb are shown in table 1. It is seen that ($-\theta ^{\prime} ({h}_{1})$) increases when Gt is assigned higher values. This fact shows that mixed convection improves heat transfer phenomena between the solid boundary and the base fluid. A similar trend is also seen for Nb and Nt. However, ($-\theta ^{\prime} ({h}_{1})$) decreases for large Gc. Moreover, it is interesting to note that results calculated by both methods, i.e. the Shooting method and NDSolve, are in good agreement.


Table 1.
Table 1.Comparison of $-\theta ^{\prime} ({h}_{1})$ by using the Shooting method and the built-in function NDSolve in Mathematica.
GtGcNtNbShooting methodNDSolveDifference
0.54.058 594.05860.000 01
1.04.065 214.065 220.000 01
1.54.084 884.084 890.000 01
0.04.08514.085 110.000 01
0.54.058 594.05860.000 01
1.04.037 634.037 640.000 01
0.23.819 073.819 080.000 01
0.43.976 663.976 670.000 01
0.64.142 694.142 710.000 02
0.23.930 513.930 520.000 01
0.44.016 464.016 470.000 01
0.64.100 284.10030.000 02

New window|CSV

5.2. Entropy generation

Figures 69 demonstrate the change in entropy generation (Ns) for different values of Gc, Gt, Nt and Nb. Figure 6 shows that through enhancing the values of Gc, entropy generation decreases. This is mainly due to a rise in buoyancy forces caused by the concentration difference, which lowers the temperature and contributes to a decrease in entropy generation. However, figure 7 indicates that entropy generation is enhanced by increasing the values of Gt. This occurs due to mixed convection caused by temperature differences, as mixed convection helps the nanoliquid in achieving higher temperatures, resulting in increased entropy production. It is visualized in figures 8 and 9 that increases in Nt and Nb enhance entropy generation. This is because kinetic energy increases with increases in Nt and Nb, which generates heat in the system, resulting in increased irreversibility in the system.

Figure 6.

New window|Download| PPT slide
Figure 6.Variation in Ns with Gc.


Figure 7.

New window|Download| PPT slide
Figure 7.Variation in Ns with Gt.


Figure 8.

New window|Download| PPT slide
Figure 8.Variation in Ns with Nt.


Figure 9.

New window|Download| PPT slide
Figure 9.Variation in Ns with Nb.


Figures 10 and 11 are outlined to show the effects of Nt and Gt on entropy generation along the channel's length. These graphs show oscillating behavior with the highest value near a wider portion of the channel. Figure 10 shows that entropy increases by improving Nt. A similar trend is seen for Gt (see figure 11). By improving buoyancy forces, the randomness in a system rises and hence the irreversibility rate increases.

Figure 10.

New window|Download| PPT slide
Figure 10.Effects of Nt on Ns versus axial distance.


Figure 11.

New window|Download| PPT slide
Figure 11.Effects of Gt on Ns versus axial distance.


Figures 12 and 13 provide the distribution of Bejan number $({B}_{e})$ for different values of Gt and Gc. These graphs show that the Bejan number reaches its minimum value near the center of the channel, while the maximum value appears near the channel walls. The Bejan number has an increasing effect by increasing the values of Gt (see figure 12). This is due to an increase in heat transfer irreversibility over total irreversibility. However, the Bejan number has a decreasing behavior for enhanced values of Gc (see figure 13).

Figure 12.

New window|Download| PPT slide
Figure 12.Variation in Be with Gt.


Figure 13.

New window|Download| PPT slide
Figure 13.Variation in Be with Gc.


Figures 14 and 15 are outlined to show the effects of Nt and Gt on the Bejan number along the channel's length. These plots exhibit sinusoidal behavior with the highest value near the wider portion of the channel. Figure 14 shows that the Bejan number increases by increasing Gt. A similar trend is also seen for Nt (see figure 15).

Figure 14.

New window|Download| PPT slide
Figure 14.Effects of Nt on Be versus axial distance.


Figure 15.

New window|Download| PPT slide
Figure 15.Effects of Gt on Be versus axial distance.


5.3. Mass transfer analysis

Figures 1619 display the consequences of Sc, ${\gamma }_{c}$, Nb, and Nt on concentration profile $(\phi )$. Figure 16 indicates that φ increases when Sc is increased. This is primarily due to an increase in momentum diffusivity. A similar trend is noted for ${\gamma }_{c}$ (see figure 17). The concentration increases when ${\gamma }_{c}$ is assigned higher values. As the chemical reaction enhances the interfacial mass transfer rate, thus the local concentration decreases, resulting in an increased concentration flux. It is inferred from figures 18 and 19 that φ increases with an enhancement in Nb and decreases for larger Nt. When Nt increases, there is a decrease in viscosity which leads to a decrease in φ.

Figure 16.

New window|Download| PPT slide
Figure 16.Variation in φ with Sc.


Figure 17.

New window|Download| PPT slide
Figure 17.Variation in φ with ${\gamma }_{c}$.


Figure 18.

New window|Download| PPT slide
Figure 18.Variation in φ with Nb.


Figure 19.

New window|Download| PPT slide
Figure 19.Variation in φ with Nt.


Table 2 shows the effects of Sc, ${\gamma }_{c}$, Nt and Nb on mass transfer rate at the channel wall ($\phi ^{\prime} ({h}_{1})$). Table 2 shows that ‘$\phi ^{\prime} ({h}_{1})$' augments by increasing the values of Nt while it decreases by increasing Sc, ${\gamma }_{c}$ and Nb. Moreover, it is interesting to note that the results calculated by both methods are in good agreement.


Table 2.
Table 2.Comparison of $\phi ^{\prime} ({h}_{1})$ by using the Shooting method and the built-in function NDSolve in Mathematica.
Sc${\gamma }_{c}$NtNbShooting methodNDSolveDifference
0.22.923 022.923030.000 01
0.42.751 892.751880.000 01
0.62.611 372.611380.000 01
0.13.115 393.115400.000 01
0.23.091 713.091720.000 01
0.33.068 673.068680.000 01
0.20.843 570.843580.000 01
0.42.197 922.197930.000 01
0.63.681 783.681790.000 01
0.27.593 697.593700.000 01
0.43.701 533.701540.000 01
0.62.403 822.403830.000 01

New window|CSV

5.4. Flow behavior analysis

Figures 2023 are designed to compute the variations in Gt, Gc, Nb and Nt on velocity profile. These plots show that an extreme value of velocity occurs near the channel center. From figure 20 it is noted that velocity is an increasing function of Gt. This is because an increase in Gr leads to a decrease in viscosity. This activity causes less resistance to flow and therefore increases velocity. In contrast to the previous case, velocity decreases with an increase in Gc (see figure 21). This is because Gc with its increasing values improves the concentration of the nanofluid and results in a decrease in velocity. Figure 22 shows an increase in velocity as Nb is enhanced. It is found that Nb lowers the viscosity of the fluid due to its heat absorbing properties and thus tends to increase the fluid velocity. A decreasing effect on velocity is reported for increasing Nt as it makes the nanofluid denser (see figure 23).

Figure 20.

New window|Download| PPT slide
Figure 20.Variation in u with Gt.


Figure 21.

New window|Download| PPT slide
Figure 21.Variation in u with Gc.


Figure 22.

New window|Download| PPT slide
Figure 22.Variation in u with Nb.


Figure 23.

New window|Download| PPT slide
Figure 23.Variation in u with Nt.


It is important to look at the nature of streamlines under the influence of various parameters that affect the flow characteristics. The most important characteristic of peristaltic flow is associated with the circulation of streamline, characterized as trapping that leads to the formation of the trapped volume of liquids, which is called a bolus. This captured bolus is carried along peristaltic waves. This tendency is extremely beneficial for transporting fluids through the body in a proper manner. This phenomenon is examined through figures 24 and 25. From these plots, it is noted that the bolus formation and streamline pattern slightly changes by enhancing Gc and Nt.

Figure 24.

New window|Download| PPT slide
Figure 24.Impacts of Gc on ψ.


Figure 25.

New window|Download| PPT slide
Figure 25.Impacts of Nt on ψ.


Table 3 displays the numerical results of Gt, Gc, Nb and Nt on skin friction at the wall. Table 3 shows that skin friction at the wall consistently augments with an increase in Gt and Nt whereas it reduces with Gc and Nb.


Table 3.
Table 3.Comparison of skin friction by using the Shooting method and the built-in function NDSolve in Mathematica.
GtGcNtNbShooting methodNDSolveDifference
0.5−1.06618−1.066170.000 01
1.0−0.73846−0.738450.000 01
1.5−0.40594−0.405930.000 01
0.0−1.01375−1.013740.000 01
0.5−1.06618−1.066170.000 01
1.0−1.11807−1.118060.000 01
0.2−0.96284−0.962830.000 01
0.4−1.03219−1.032180.000 01
0.6−1.09966−1.099680.000 01
0.2−1.32612−1.326110.000 01
0.4−1.10898−1.108970.000 01
0.6−1.03813−1.038120.000 01

New window|CSV

6. Conclusions

Entropy generation rate, heat and mass transfer analysis for mixed convective peristaltic transportation of a reactive nanofluid under the effects of thermophoresis and Brownian motion are investigated. The major findings produced from the present study are: The temperature of a nanofluid increases by increasing the Brownian motion parameter.
Entropy generation enhances by improving the Brownian motion parameter.
An increase in the thermophoresis parameter produces more entropy.
The Bejan number has an increasing effect for increasing temperature Grashof number.
Concentration of nanoparticles improves when the Brownian motion parameter is increased, and decreases by enhancing the thermophoresis parameter.
An increase in the chemical reaction parameter improves the concentration characteristics of nanoparticles.
The mass transfer rate at the wall enhances by increasing the thermophoresis parameter. On the other hand, it decreases by enhancing the Schmidt number and chemical reaction parameter.
An augmentation in axial velocity is seen for a large temperature Grashof number while it decreases by increasing the concentration Grashof number.
Velocity increases by enhancing the Brownian motion parameter while it decreases for a large thermophoresis parameter.


Reference By original order
By published year
By cited within times
By Impact factor

Choi S U Eastman J A 1995 Enhancing Thermal Conductivity of Fluids with Nanoparticles (No. ANL/MSD/CP-84938; CONF-951135-29) Lemont, IL Argonne National Lab
[Cited within: 1]

Buongiorno J 2006 Convective Transport in Nanofluids 128 240 250
DOI:10.1115/1.2150834 [Cited within: 1]

Prasad L K O'Mary H Cui Z 2015 Nanomedicine delivers promising treatments for rheumatoid arthritis
Nanomedicine 10 2063 2074

DOI:10.2217/nnm.15.45 [Cited within: 1]

Kleinstreuer C Feng Y 2011 Experimental and theoretical studies of nanofluid thermal conductivity enhancement: a review
Nanoscale Res. Lett. 6 1 13

DOI:10.1186/1556-276X-6-229 [Cited within: 1]

Hassan M Marin M Alsharif A Ellahi R 2018 Convective heat transfer flow of nanofluid in a porous medium over wavy surface
Phys. Lett. A 382 2749 2753

DOI:10.1016/j.physleta.2018.06.026 [Cited within: 1]

Vallejo J P Sani E Zyla G Lugo L 2019 Tailored silver/graphene nanoplatelet hybrid nanofluids for solar applications
J. Mol. Liq. 296 112007

DOI:10.1016/j.molliq.2019.112007 [Cited within: 1]

Irandoost Shahrestani M Maleki A Safdari Shadloo M Tlili I 2020 Numerical investigation of forced convective heat transfer and performance evaluation criterion of ${\mathrm{Al}}_{2}{{\rm{O}}}_{3}/\mathrm{water}$ nanofluid flow inside an axisymmetric microchannel
Symmetry 12 120

DOI:10.3390/sym12010120 [Cited within: 1]

Latham T W 1966 Fluid motion in a peristaltic pump
[MS Thesis]Cambridge M A

[Cited within: 2]

Shapiro A H Jaffrin M Y Weinberg S L 1969 Peristaltic pumping with long wavelengths at low Reynolds number
J. Fluid Mech. 37 799 825

DOI:10.1017/S0022112069000899 [Cited within: 1]

Akbar N S Nadeem S 2011 Endoscopic effects on peristaltic flow of a nanofluid
Commun. Theor. Phys. 56 761

DOI:10.1088/0253-6102/56/4/28 [Cited within: 1]

Noreen S 2013 Mixed convection peristaltic flow of third order nanofluid with an induced magnetic field
PLoS One 8 e78770

DOI:10.1371/journal.pone.0078770 [Cited within: 1]

Abbasi F M Hayat T Ahmad B 2015 Peristaltic transport of copper water nanofluid saturating porous medium
Physica E 67 47 53

DOI:10.1016/j.physe.2014.11.002

Gnaneswara Reddy M Kumara B C Makinde O D 2017 Cross diffusion impacts on hydromagnetic radiative peristaltic Carreau–Casson nanofluids flow in an irregular channel
Defect and Diffusion Forum vol 377 Bäch, Switzerland Trans Tech Publications Ltd 62 83

DOI:10.4028/www.scientific.net/DDF.377.62

Makinde O D Reddy M G Reddy K V 2017 Effects of thermal radiation on MHD peristaltic motion of Walters-B fluid with heat source and slip conditions
J. Appl. Fluid Mech. 10 1105 1112

DOI:10.18869/acadpub.jafm.73.241.27082

Reddy K V Makinde O D Reddy M G 2018 Thermal analysis of MHD electro-osmotic peristaltic pumping of Casson fluid through a rotating asymmetric micro-channel
Indian J. Phys. 92 1439 1448

DOI:10.1007/s12648-018-1209-1

Kothandapani M Prakash J 2015 Effects of thermal radiation parameter and magnetic field on the peristaltic motion of Williamson nanofluids in a tapered asymmetric channel International
J. Heat Mass Transfer 81 234 245

DOI:10.1016/j.ijheatmasstransfer.2014.09.062 [Cited within: 3]

Mekheimer K S Hasona W M Abo-Elkhair R E Zaher A Z 2018 Peristaltic blood flow with gold nanoparticles as a third grade nanofluid in catheter: applications of cancer therapy
Phys. Lett. A 382 85 93

DOI:10.1016/j.physleta.2017.10.042

Prakash J Siva E P Tripathi D Kothandapani M 2019 Nanofluids flow driven by peristaltic pumping in occurrence of magnetohydrodynamics and thermal radiation
Mater. Sci. Semicond. Process. 100 290 300

DOI:10.1016/j.mssp.2019.05.017

Zeeshan A Bhatti M M Muhammad T Zhang L 2020 Magnetized peristaltic particle–fluid propulsion with Hall and ion slip effects through a permeable channel
Physica A 550 123999

DOI:10.1016/j.physa.2019.123999 [Cited within: 1]

Riaz A Zeeshan A Bhatti M M Ellahi R 2020 Peristaltic propulsion of Jeffrey nano-liquid and heat transfer through a symmetrical duct with moving walls in a porous medium
Physica A 545 123788

DOI:10.1016/j.physa.2019.123788 [Cited within: 1]

Beg O A Tripathi D 2011 Mathematica simulation of peristaltic pumping with double-diffusive convection in nanofluids: a bio-nano-engineering model
Proc. Inst. Mech. Eng., Part N 225 99 114

DOI:10.1177/1740349912437087 [Cited within: 1]

Hayat T Abbasi F M Al-Yami M Monaquel S 2014 Slip and Joule heating effects in mixed convection peristaltic transport of nanofluid with Soret and Dufour effects
J. Mol. Liq. 194 93 99

DOI:10.1016/j.molliq.2014.01.021 [Cited within: 1]

Ellahi R Zeeshan A Hussain F Asadollahi A 2019 Peristaltic blood flow of couple stress fluid suspended with nanoparticles under the influence of chemical reaction and activation energy
Symmetry 11 276

DOI:10.3390/sym11020276 [Cited within: 2]

Khan M I Hayat T Afzal S Khan M I Alsaedi A 2020 Theoretical and numerical investigation of Carreau-Yasuda fluid flow subject to Soret and Dufour effects
Comput. Methods Programs Biomed. 186 105145

DOI:10.1016/j.cmpb.2019.105145 [Cited within: 2]

Hayat T Nawaz S Alsaedi A 2020 Entropy analysis for the peristalsis flow with homogeneous–heterogeneous reaction
Eur. Phys. J. Plus 135 1 16

DOI:10.1140/epjp/s13360-020-00293-z [Cited within: 1]

Bejan A 1979 A study of entropy generation in fundamental convective heat transfer

DOI:10.1115/1.3451063 [Cited within: 1]

Nadeem S Sadia H 2013 Influence of temperature dependent viscosity and entropy generation on the flow of a Johnson–Segalman fluid
Walailak J. Sci. Technol. 10 553 579

[Cited within: 1]

Akram S Zafar M Nadeem S 2018 Peristaltic transport of a Jeffrey fluid with double-diffusive convection in nanofluids in the presence of inclined magnetic field
Int. J. Geom. Meth. Mod. Phys. 15 1850181

DOI:10.1142/S0219887818501815 [Cited within: 1]

Abbas M A Bai Y Rashidi M M Bhatti M M 2016 Analysis of entropy generation in the flow of peristaltic nanofluids in channels with compliant walls
Entropy 18 90

DOI:10.3390/e18030090 [Cited within: 1]

Farooq S Hayat T Alsaedi A Asghar S 2018 Mixed convection peristalsis of carbon nanotubes with thermal radiation and entropy generation
J. Mol. Liq. 250 451 467

DOI:10.1016/j.molliq.2017.11.179 [Cited within: 1]

Akbar Y Abbasi F M Shehzad S A 2020 Thermal radiation and Hall effects in mixed convective peristaltic transport of nanofluid with entropy generation
Appl. Nanoscience 10 5421 5433

DOI:10.1007/s13204-020-01446-3 [Cited within: 1]

Hayat T Nisar Z Alsaedi A Ahmad B 2021 Analysis of activation energy and entropy generation in mixed convective peristaltic transport of Sutterby nanofluid
J. Therm. Anal. Calorim. 143 1867 1880

DOI:10.1007/s10973-020-09969-1 [Cited within: 2]

Akbar Y Abbasi F M 2020 Impact of variable viscosity on peristaltic motion with entropy generation
Int. Commun. Heat Mass Transfer 118 104826

DOI:10.1016/j.icheatmasstransfer.2020.104826 [Cited within: 5]

Akbar Y Shanakhat I Abbasi F M Shehzad S A 2020 Entropy generation analysis for radiative peristaltic motion of silver–water nanomaterial with temperature dependent heat sink/source
Phys. Scr. 95 115201

DOI:10.1088/1402-4896/abbbce

Akbar Y Abbasi F M Shehzad S A 2020 Thermodynamical analysis for mixed convective peristaltic motion of silver–water nanoliquid having temperature dependent electrical conductivity
Appl. Nanoscience1 12

DOI:10.1007/s13204-020-01569-7

Akbar Y Abbasi F M Shehzad S A 2020 Effectiveness of Hall current and ion slip on hydromagnetic biologically inspired flow of $\mathrm{Cu}-{\mathrm{Fe}}_{3}{{\rm{O}}}_{4}/{{\rm{H}}}_{2}{\rm{O}}$ hybrid nanomaterial
Phys. Scr. 96 025210

DOI:10.1088/1402-4896/abcff1

Riaz A Bhatti M M Ellahi R Zeeshan A M Sait S 2020 Mathematical analysis on an asymmetrical wavy motion of blood under the influence entropy generation with convective boundary conditions
Symmetry 12 102

DOI:10.3390/sym12010102

Hayat T Nawaz S Alsaedi A Fardoun H M 2020 Entropy optimization for peristalsis of Rabinowitsch nanomaterial
Appl. Nanoscience 10 4177 4190

DOI:10.1007/s13204-020-01535-3 [Cited within: 2]

Akbar Y Abbasi F MIrreversibility analysis of nanofluid flow induced by peristaltic waves in the presence of concentration-dependent viscosity
Heat Transfer1 18

DOI:10.1002/htj.22133 [Cited within: 1]

相关话题/Effectiveness transfer mixed