Y Akbar¹, F M Abbasi¹, U M Zahid¹, S A Shehzad^,21Department of Mathematics, COMSATS University Islamabad, Islamabad, Pakistan
²Department 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 [11–19]. 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 [33–39].

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]:(1)$\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}$(2)$\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}$, a₁, a₂, 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]:(3)$\begin{eqnarray}\displaystyle \frac{\partial \overline{U}}{\partial \overline{X}}+\displaystyle \frac{\partial \overline{V}}{\partial \overline{Y}}=0,\end{eqnarray}$(4)$\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}$(5)$\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}$(6)$\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}$(7)$\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, C_f, $\hat{{\gamma }_{c}}$, D_B, D_T, 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:(8)$\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:(9)$\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}$(10)$\begin{eqnarray}0=-\displaystyle \frac{\partial p}{\partial y},\end{eqnarray}$(11)$\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}$(12)$\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 ψ, G_t, Br, G_c, Sc, N_b, ${\gamma }_{c}$, N_t, 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]:(13)$\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]:(14)$\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(15)$\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]:(16)$\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]:(17)$\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]:(18)$\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:(19)$\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}$(20)$\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}$(21)$\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*}$(22)$\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*}$(23)$\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⁻⁶. 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 2–5. 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 (G_c). 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 (G_t). 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 N_b. The reason is that kinetic energy is augmented as N_b 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 (N_t) (see figure 5). For higher N_t, 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 G_c.

Figure 3.

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

Figure 4.

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

Figure 5.

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


Numerical values of the heat transfer rate at the channel wall ($-\theta ^{\prime} ({h}_{1})$) for variations in G_t, G_c, N_t, and N_b are shown in table 1. It is seen that ($-\theta ^{\prime} ({h}_{1})$) increases when G_t 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 N_b and N_t. However, ($-\theta ^{\prime} ({h}_{1})$) decreases for large G_c. 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.

Gt	Gc	Nt	Nb	Shooting method	NDSolve	Difference
0.5				4.058 59	4.0586	0.000 01
1.0				4.065 21	4.065 22	0.000 01
1.5				4.084 88	4.084 89	0.000 01
	0.0			4.0851	4.085 11	0.000 01
	0.5			4.058 59	4.0586	0.000 01
	1.0			4.037 63	4.037 64	0.000 01
		0.2		3.819 07	3.819 08	0.000 01
		0.4		3.976 66	3.976 67	0.000 01
		0.6		4.142 69	4.142 71	0.000 02
			0.2	3.930 51	3.930 52	0.000 01
			0.4	4.016 46	4.016 47	0.000 01
			0.6	4.100 28	4.1003	0.000 02

New window|CSV

5.2. Entropy generation

Figures 6–9 demonstrate the change in entropy generation (N_s) for different values of G_c, G_t, N_t and N_b. Figure 6 shows that through enhancing the values of G_c, 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 G_t. 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 N_t and N_b enhance entropy generation. This is because kinetic energy increases with increases in N_t and N_b, which generates heat in the system, resulting in increased irreversibility in the system.

Figure 6.

New window|Download| PPT slide
Figure 6.Variation in N_s with G_c.

Figure 7.

New window|Download| PPT slide
Figure 7.Variation in N_s with G_t.

Figure 8.

New window|Download| PPT slide
Figure 8.Variation in N_s with N_t.

Figure 9.

New window|Download| PPT slide
Figure 9.Variation in N_s with N_b.


Figures 10 and 11 are outlined to show the effects of N_t and G_t 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 N_t. A similar trend is seen for G_t (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 N_t on N_s versus axial distance.

Figure 11.

New window|Download| PPT slide
Figure 11.Effects of G_t on N_s versus axial distance.


Figures 12 and 13 provide the distribution of Bejan number $({B}_{e})$ for different values of G_t and G_c. 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 G_t (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 G_c (see figure 13).

Figure 12.

New window|Download| PPT slide
Figure 12.Variation in B_e with G_t.

Figure 13.

New window|Download| PPT slide
Figure 13.Variation in B_e with G_c.


Figures 14 and 15 are outlined to show the effects of N_t and G_t 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 G_t. A similar trend is also seen for N_t (see figure 15).

Figure 14.

New window|Download| PPT slide
Figure 14.Effects of N_t on B_e versus axial distance.

Figure 15.

New window|Download| PPT slide
Figure 15.Effects of G_t on B_e versus axial distance.

5.3. Mass transfer analysis

Figures 16–19 display the consequences of S_c, ${\gamma }_{c}$, N_b, and N_t on concentration profile $(\phi )$. Figure 16 indicates that φ increases when S_c 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 N_b and decreases for larger N_t. When N_t 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 S_c.

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 N_b.

Figure 19.

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


Table 2 shows the effects of S_c, ${\gamma }_{c}$, N_t and N_b 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 N_t while it decreases by increasing S_c, ${\gamma }_{c}$ and N_b. 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}$	Nt	Nb	Shooting method	NDSolve	Difference
0.2				2.923 02	2.92303	0.000 01
0.4				2.751 89	2.75188	0.000 01
0.6				2.611 37	2.61138	0.000 01
	0.1			3.115 39	3.11540	0.000 01
	0.2			3.091 71	3.09172	0.000 01
	0.3			3.068 67	3.06868	0.000 01
		0.2		0.843 57	0.84358	0.000 01
		0.4		2.197 92	2.19793	0.000 01
		0.6		3.681 78	3.68179	0.000 01
			0.2	7.593 69	7.59370	0.000 01
			0.4	3.701 53	3.70154	0.000 01
			0.6	2.403 82	2.40383	0.000 01

New window|CSV

5.4. Flow behavior analysis

Figures 20–23 are designed to compute the variations in G_t, G_c, N_b and N_t 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 G_t. This is because an increase in G_r 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 G_c (see figure 21). This is because G_c 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 N_b is enhanced. It is found that N_b 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 N_t as it makes the nanofluid denser (see figure 23).

Figure 20.

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

Figure 21.

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

Figure 22.

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

Figure 23.

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


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 G_c and N_t.

Figure 24.

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

Figure 25.

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


Table 3 displays the numerical results of G_t, G_c, N_b and N_t on skin friction at the wall. Table 3 shows that skin friction at the wall consistently augments with an increase in G_t and N_t whereas it reduces with G_c and N_b.


Table 3.
Table 3.Comparison of skin friction by using the Shooting method and the built-in function NDSolve in Mathematica.

Gt	Gc	Nt	Nb	Shooting method	NDSolve	Difference
0.5				−1.06618	−1.06617	0.000 01
1.0				−0.73846	−0.73845	0.000 01
1.5				−0.40594	−0.40593	0.000 01
	0.0			−1.01375	−1.01374	0.000 01
	0.5			−1.06618	−1.06617	0.000 01
	1.0			−1.11807	−1.11806	0.000 01
		0.2		−0.96284	−0.96283	0.000 01
		0.4		−1.03219	−1.03218	0.000 01
		0.6		−1.09966	−1.09968	0.000 01
			0.2	−1.32612	−1.32611	0.000 01
			0.4	−1.10898	−1.10897	0.000 01
			0.6	−1.03813	−1.03812	0.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 View Option By original order
By published year
By cited within times
By Impact factor

[1]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]

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

[3]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]

[4]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]

[5]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]

[6]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]

[7]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]

[8]Latham T W 1966 Fluid motion in a peristaltic pump
[MS Thesis]Cambridge M A
[Cited within: 2]

[9]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]

[10]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]

[11]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]

[12]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

[13]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

[14]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

[15]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

[16]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]

[17]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

[18]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

[19]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]

[20]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]

[21]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]

[22]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]

[23]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]

[24]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]

[25]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]

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

DOI:10.1115/1.3451063 [Cited within: 1]

[27]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]

[28]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]

[29]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]

[30]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]

[31]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]

[32]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]

[33]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]

[34]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

[35]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

[36]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

[37]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

[38]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]

[39]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]