Peristaltic Blood Flow with Gold Nanoparticles on a Carreau Nanofluid through a Non-Darcian Porous Medium

Show more

1. Introduction

Cancer is a dangerous and deadly to most of its patients. Recent studies have shown that gold nanoparticles (GNP) can cure and overcome it because these particles have high atomic numbers which produce heat and leads to treatment of malignancy tumors. GNP has many properties that are very essential in cancer therapy. They are very small and can penetrate widely throughout a body. They have been considered as major area of research due to their applications in biomedical sciences [1] . Huang and El-Sayed [2] studied an important optical and photothermal property of gold nanoparticles in different shapes and structures and address some recent applications for cancer imaging, spectroscopic detection and photothermal therapy. The flow of GNP through two coaxial tubes has been investigated by some authors. Hamzehnezhad et al. [3] and Hatami et al. [4] studied the third-grade non-Newtonian fluid conveying GNP in a porous and hollow vessel by various analytical methods which is essential idea of the present paper.

The studying is that Darcian and non-Darcian flow in Newtonian and non-Newtonian fluid flow saturating porous medium. This because of their wide applications such as chemical reactors, building non-Darcian fluids in saturated porous media meets a great importance. The Dag insulation, packed bed, enhanced oil recovery, food technology and filtration processes. Eldabe and Abu Zeid [5] have studied a non-Darcian Couette flow through a porous medium of magnetohydrodynamic visco-elastic fluid with heat and mass transfer. They found the solutions of velocity, temperature and nanoparticles distribution by using the homotopy perturbation method. The effect of the boundary conditions on a peristaltic flow of a nanofluid in an asymmetric channel embedded with a non-Darcian porous medium was scrutinized by Rathod and Sanjeevkumar [6] . Wu et al. [7] studied the analysis of multiphase Non-Darcian flow in porous media.

In the literature, numerous studies regarding peristaltic flow of a Carreau nanofluid have been done for Non-Newtonian fluids. Akbar et al. [8] studied the numerical simulation of peristaltic flow of a Carreau nanofluid in an asymmetric channel. They found solutions for the distributions of stream function, pressure rise, and temperature and nanoparticle volume fraction by using the fourth and fifth order Runge-Kutta-Fehlberg. The effects of Nanofluid on peristaltic flow of a Carreau fluid model in an inclined magnetic field are proposed by Akram [9] . He found the solutions of the simplified coupled nonlinear equations using an analytical approach. Eldabe et al. [10] scrutinized the peristaltic transport of magnetohydrodynamic Carreau nanofluid with heat and mass transfer inside asymmetric channel.

Peristaltic flow of a nanofluid under the effect of the hall current and porous medium was studied by Nowar [11] . Abo-Eldahab et al. [12] introduced the hall currents and heat transfer effects on peristaltic transport in a vertical asymmetric channel through a porous medium. Steh and Ghosh [13] addressed the combined influence of rotation and Hall current in the presence of a transverse magnetic field, neglecting the induced magnetic field. Ghosh [14] discussed the effects of Hall current on an unsteady hydromagnetic flow in a rotating channel permeated by an inclined magnetic field in the presence of an oscillator.

Several applications in engineering occur at high temperature through variable thermal conductivity. Few such operations include nuclear power plants, in turbines, rockets, space vehicles and pumps operated at high temperatures (see Ref. [15] [16] ). In addition, studying velocity slip effects are adequate for certain flow problems at the walls of the pipes/channels.

The governing equations of motion become more complex to handle as supplemental nonlinear terms. So, the exact solutions of such problems are practically impossible, so we turn to find a new method such as the differential transform method (DTM). DTM is semi-analytical methods appear in the equations of motion, because of the flow behavior of Non-Newtonian nanofluids. This method was proposed by Zhou [17] , for solving ordinary, partial differential and integral equations. Different applications of DTM can be found in [18] . The Ms-DTM method accelerates the convergence of the series solution over large region and yields a series solution, this series will be truncated due to the required accuracy of solutions. This modified technique is verified through illustrative examples of non-chaotic or chaotic systems by Odibat et al. [19] .

To the best of our knowledge, the effect of the Hall current, variable thermal conductivity on peristaltic flow of Carreau nanofluid has not been discussed before. Motivated by the above discussions, the aim of the present paper is to examine the MHD peristaltic flow of Carreau nanofluid accompanying heat and mass transfer in the presence of viscous dissipation. In order to clarify the problem at hand, the mathematical modeling is affirmed in Section 2, solutions of the problem on using the multi-steps DTM are presented in Section 3. In Section 4 the influence of the various parameters, the pressure gradient, pressure rise, longitudinal velocity, temperature and concentration are explored through figures and tables. The main findings have been abbreviated in Section 5.

2. Formulation of the Problem

Consider a peristaltic flow of an incompressible Carreau nanofluid in a two dimensional asymmetric channel of width ${d}_{1}+{d}_{2}$. It is convenient to work with the Cartesian coordinates $\left(X,Y\right)$ where X-axis and Y-axis are the horizontal and vertical axes, respectively. The right and left walls have a constant temperature of ${T}_{1}$ and ${T}_{0}$, respectively. Also, they have nanoparticle volume fractions of ${C}_{1}$ and ${C}_{0}$. Asymmetry in the channel is produced by propagation of waves along the channel walls and traveling with different amplitudes, phases and of the constant speed. The sketch of the model is given in Figure 1. In the Cartesian coordinates system $\left(X,Y\right)$, the right wall $Y={H}_{1}$ and the left wall $Y={H}_{2}$. They are given from the following relations:

$Y={H}_{1}={d}_{1}+{a}_{1}\mathrm{cos}\left[\frac{2\text{\pi}}{\lambda}\left(X-ct\right)\right]$ Right hand side (1)

$Y={H}_{2}=-{d}_{2}-{b}_{1}\mathrm{cos}\left[\frac{2\text{\pi}}{\lambda}\left(X-ct\right)+\varphi \right]$ Left hand side (2)

where ${a}_{1}$ and ${b}_{1}$ are the amplitudes of the waves, $\lambda $ is the wave length, the phase difference $\varphi $ diverse in the range $\left(0\le \varphi \le \text{\pi}\right)$, $\varphi =0$, corresponds to

Figure 1. Geometry of physical model.

symmetric channel with waves out of phase and $\varphi =\text{\pi}$ the waves are in phase, and further ${a}_{1},{b}_{1},{d}_{1},{d}_{2}$ and satisfies the condition.

${a}_{1}^{2}+{b}_{1}^{2}+2{a}_{1}{b}_{1}\mathrm{cos}\varphi \le \left({d}_{1}^{2}+{d}_{2}^{2}\right)$ (3)

A strong uniform magnetic field with magnetic flux density $B=\left(0,0,{B}_{0}\right)$ is applied and the Hall effects are taken, also into account. Assuming that no applied polarization voltage, so that the electric field ignored. In case of the hall current, as given from generalization of the ohm [20] :

$\stackrel{\xaf}{J}=\sigma \left[E+\stackrel{\xaf}{V}\times \stackrel{\xaf}{B}-\frac{1}{e\left({n}_{e}\right)}\left(\stackrel{\xaf}{J}\times \stackrel{\xaf}{B}\right)\right]$, (4)

where $\sigma $ is the electrical conductivity of the fluid, $\stackrel{\xaf}{V}$ is the velocity vector, e is the electric charge of electrons, ${n}_{e}$ is the number density of the electrons. Equation (3) can be solved in $\stackrel{\xaf}{J}$ to yield the Lorentz force vector in the form:

$\stackrel{\xaf}{J}\times \stackrel{\xaf}{B}=\frac{\sigma {B}_{0}^{2}}{1+{m}^{2}}\left[\left(V+mU\right)\underset{\_}{i}+\left(mV-U\right)\underset{\_}{j}\right],$ (5)

where U and V are the X and Y components of the velocity vector $V=\left({V}_{x},{V}_{y},0\right)$, $m=\frac{\sigma {B}_{0}}{e\left({n}_{e}\right)}$ is the hall parameter.

The variable thermal conductivity, with dimensionless temperature $\theta $ is taken into account where,

$k\left(\theta \right)={k}_{0}\left(1+\u03f5\theta \right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\u03f5<1$ (6)

Here $\u03f5=\left({T}_{m}-{T}_{w}\right)$ is the thermal conductivity variation parameter, ${k}_{0}$ is the thermal conductivity constant parameter at $\u03f5=0$.

The constitutive equation for a Carreau fluid is given by [8] .

$\frac{\eta -{\eta}_{\infty}}{{\eta}_{0}-{\eta}_{\infty}}=\left[1+{\left(\Gamma \stackrel{\dot{}}{\gamma}\right)}^{2}\right]$, (7a)

and

${\tau}_{ij}={\eta}_{0}\left[1+\frac{n-1}{2}{\left(\Gamma \stackrel{\dot{}}{\gamma}\right)}^{2}\right]{\stackrel{\dot{}}{\gamma}}_{ij},$ (7b)

In which ${\tau}_{ij}$ is the extra stress tensor, ${\eta}_{\infty}$ is the infinite shear rate viscosity, ${\eta}_{0}$ is the zero shear rate viscosity, $\Gamma $ is the time constant, n is the power law index, and $\stackrel{\dot{}}{\gamma}$ is defined as:

$\stackrel{\dot{}}{\gamma}=\sqrt{\frac{1}{2}{\displaystyle {\sum}_{i}{\displaystyle {\sum}_{j}{\stackrel{\dot{}}{\gamma}}_{ij}{\stackrel{\dot{}}{\gamma}}_{ji}}}}=\sqrt{\frac{1}{2}{{\displaystyle \prod}}^{\text{}}.}$ (8)

where ${{\displaystyle \prod}}^{\text{}}.=trac{\left(gradV+{\left(gradV\right)}^{\text{T}}\right)}^{2}$ is the second invariant strain tensor. For considering the natural convection, the nanofluid density $\rho $ is defined by [2] as:

$\rho =\phi {\rho}_{p}+\left(1+\phi \right){\rho}_{{f}_{0}}\cong \phi {\rho}_{p}+\left(1+\phi \right)\left[{\rho}_{f}\left(1-{\beta}_{T}\left(\vartheta -{\theta}_{w}\right)\right)\right]$, (9)

where ${\rho}_{p}$ is the density of the gold nanoparticles, ${\rho}_{{f}_{0}}$ is the base fluid’s density, $\vartheta $ is the temperature distribution, ${\theta}_{w}$ is the reference temperature, ${\beta}_{T}$ is the volumetric coefficient of the expansion and ${\rho}_{f}$ is the nanofluid’s density at the reference temperature, $\phi $ is the concentration of nanoparticles.

The governing equations of peristaltic flow in the case of an incompressible flow of Carreau nanofluid are given by [4] [8] :

The continuity equation:

$\frac{\partial U}{\partial X}+\frac{\partial V}{\partial Y}=0,$ (10)

The conversation of momentum yields

$\begin{array}{l}{\rho}_{f}\left[\frac{\partial U}{\partial t}+U\frac{\partial U}{\partial X}+V\frac{\partial U}{\partial Y}\right]\\ =-\frac{\partial P}{\partial X}+\frac{\partial}{\partial X}\left({\tau}_{XX}\right)+\frac{\partial}{\partial Y}\left({\tau}_{XY}\right)-\frac{\sigma {B}_{0}^{2}}{1+{m}^{2}}\left(V+mU\right)+\phi {\rho}_{p}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\left(1+\phi \right)\left[{\rho}_{f}\left(1-{\beta}_{T}\left(\vartheta -{\theta}_{w}\right)\right)\right]g-\frac{{\eta}_{0}}{{k}_{1}}-\frac{{\rho}_{f}{C}_{b}}{\sqrt{{k}_{1}}}{U}^{2},\end{array}$ (11)

$\begin{array}{l}{\rho}_{f}\left[\frac{\partial V}{\partial t}+U\frac{\partial V}{\partial X}+V\frac{\partial V}{\partial Y}\right]\\ =-\frac{\partial P}{\partial Y}+\frac{\partial}{\partial X}\left({\tau}_{YX}\right)+\frac{\partial}{\partial Y}\left({\tau}_{YY}\right)+\frac{\sigma {B}_{0}^{2}}{1+{m}^{2}}\left(mV-U\right)-\frac{{\mu}_{f}}{{k}_{1}}V-\frac{{\rho}_{f}{C}_{b}}{\sqrt{{k}_{1}}}{V}^{2},\end{array}$ (12)

The energy equation gives

$\begin{array}{l}{\left(\rho c\right)}_{f}\left[\frac{\partial T}{\partial t}+U\frac{\partial T}{\partial X}+V\frac{\partial T}{\partial Y}\right]\\ =\frac{\partial}{\partial X}\left(k\left(T\right)\frac{\partial T}{\partial X}\right)+\frac{\partial}{\partial Y}\left(k\left(T\right)\frac{\partial T}{\partial y}\right)+\frac{\sigma {B}_{0}^{2}}{1+{m}^{2}}\left({U}^{2}+{V}^{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+\left[{\tau}_{XX}\frac{\partial U}{\partial X}+{\tau}_{XY}\frac{\partial U}{\partial Y}+{\tau}_{YX}\frac{\partial V}{\partial Y}+{\tau}_{YY}\frac{\partial V}{\partial Y}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+{\left(\rho c\right)}_{p}\left[\frac{{D}_{T}}{{T}_{0}}\left[\frac{{\partial}^{2}T}{\partial {X}^{2}}+\frac{{\partial}^{2}T}{\partial {Y}^{2}}\right]+{D}_{B}\left(\frac{\partial C}{\partial X}\frac{\partial T}{\partial X}+\frac{\partial C}{\partial Y}\frac{\partial T}{\partial Y}\right)\right],\end{array}$ (13)

The concentration equation gives

$\left[\frac{\partial C}{\partial t}+U\frac{\partial C}{\partial X}+V\frac{\partial C}{\partial Y}\right]={D}_{B}\left(\frac{{\partial}^{2}C}{\partial {X}^{2}}+\frac{{\partial}^{2}C}{\partial {Y}^{2}}\right)+\frac{{D}_{T}}{{T}_{0}}\left[\frac{{\partial}^{2}T}{\partial {X}^{2}}+\frac{{\partial}^{2}T}{\partial {Y}^{2}}\right].$ (14)

where P is the pressure, ${k}_{1}$ is the permeability of the porous medium, g is acceleration of the gravity, T is the temperature of the fluid and C is the nanoparticles concentration. Both the magnetic field and channel are inclined at angle, ${\left(\rho c\right)}_{f}$ is the heat capacity of the fluid, k is the thermal conductivity, (c) p is the effective heat capacity of the nanoparticles material, ${D}_{T}$ is the thermophoretic diffusion and ${D}_{B}$ is the Brownian diffusion coefficient. Using similar arguments as given by [2] [8] , we may assume the following transformations, $x=X-ct$, $y=Y$, $u=U-c$, $v=V$ and $p\left(x\right)=P\left(X,t\right)$.

To simplify the equations of motion, we may consider the following dimensionless quantities:

$\begin{array}{l}\stackrel{\xaf}{x}=\frac{x}{\lambda},\stackrel{\xaf}{y}=\frac{y}{{d}_{1}},\stackrel{\xaf}{u}=\frac{u}{c},\stackrel{\xaf}{v}=\frac{v}{c\delta},\stackrel{\xaf}{p}=\frac{{d}_{1}^{2}}{\mu c\lambda},{h}_{1}=\frac{{H}_{1}}{{d}_{1}},{h}_{2}=\frac{{H}_{2}}{{d}_{2}},d=\frac{{d}_{2}}{{d}_{1}},a=\frac{{a}_{1}}{{d}_{1}},\\ b=\frac{{b}_{1}}{{d}_{1}},{R}_{e}=\frac{{\rho}_{f}c{d}_{1}}{{\eta}_{0}},\stackrel{\xaf}{\psi}=\frac{\psi}{c{d}_{1}},{\stackrel{\xaf}{\tau}}_{xx}=\frac{\lambda}{{\eta}_{0}c}{\tau}_{xx},{\stackrel{\xaf}{\tau}}_{xy}=\frac{{d}_{1}}{{\eta}_{0}c}{\tau}_{xy},{\stackrel{\xaf}{\tau}}_{yy}=\frac{{d}_{1}}{{\eta}_{0}c}{\tau}_{yy},\\ {W}_{e}=\frac{\Gamma c}{{d}_{1}},\stackrel{\xaf}{\stackrel{\dot{}}{\gamma}}=\frac{{d}_{1}\stackrel{\dot{}}{\gamma}}{c},s=\frac{{k}_{1}}{{d}_{1}^{2}},\alpha =\frac{K}{{\left(\rho c\right)}_{f}},\theta =\frac{{\theta}^{\prime}-{\theta}_{w}}{{\theta}_{m}-{\theta}_{w}},\phi =\frac{{\phi}^{\prime}-{\phi}_{w}}{{\phi}_{m}-{\phi}_{w}},{P}_{r}=\frac{{\eta}_{0}}{\alpha {\rho}_{f}},\\ {G}_{r}=\frac{{\rho}_{f}{d}_{1}^{2}g\left({\theta}_{m}-{\theta}_{w}\right)\left(1-{\phi}_{m}\right)}{c{\eta}_{0}},{B}_{r}=\frac{{a}_{1}^{2}\rho \left({\rho}_{p}-{\rho}_{f}g\right)\left({\phi}_{m}-{\phi}_{w}\right)}{c{\eta}_{0}},\\ {N}_{b}={D}_{B}\left({\phi}_{m}-{\phi}_{w}\right),{N}_{t}=\frac{{D}_{T}\left({\theta}_{m}-{\theta}_{w}\right)}{{\theta}_{w}},{C}_{f}=\frac{{\rho}_{f}{C}_{b}{d}_{1}^{2}c}{{\eta}_{0}\sqrt{K}}.\end{array}$ (15)

where $\delta $ is the dimensionless wave number, ${R}_{e}$ is the Reynolds number, ${P}_{r}$ is the Prandtl number, ${C}_{f}$ is the Forchheimer number, ${G}_{r}$ is the local Grashof number, ${B}_{r}$ is the local nanoparticles Grash of number, ${N}_{b}$ is the Brownian motion parameter and ${N}_{t}$ is the thermophoresis parameter, ${W}_{e}$ is the Weissenberg number, M is the magnetic parameter. Here ${\theta}_{m}$ and ${\phi}_{m}$, denote the pipe temperature, fluid temperature and mass concentration, respectively. As usual in fluid mechanics, the stream function $\psi =\psi \left(x,y\right)$ may be taken as

$u=\frac{\partial \psi}{\partial y}$ and $v=-\frac{\partial \psi}{\partial x}$.

${\tau}_{xy}=\left[1+\frac{n-1}{2}{W}_{e}^{2}{\stackrel{\dot{}}{\gamma}}^{2}\right]\left(\frac{{\partial}^{2}\psi}{\partial {y}^{2}}-{\delta}^{2}\frac{{\partial}^{2}\psi}{\partial {x}^{2}}\right)={\tau}_{yx},$ (16)

The governing equations of the Carreau nanofluid are:

$\begin{array}{l}-\frac{\partial p}{\partial x}+\frac{\partial}{\partial y}\left[1+\frac{n-1}{2}{W}_{e}^{2}{\left(\frac{{\partial}^{2}\psi}{\partial {y}^{2}}\right)}^{2}\right]\frac{{\partial}^{2}\psi}{\partial {y}^{2}}-\left(\frac{{M}^{2}}{1+{m}^{2}}+2{C}_{f}\right)\frac{\partial \psi}{\partial y}\\ \text{\hspace{0.05em}}-{C}_{f}\frac{{\partial}^{2}\psi}{\partial {y}^{2}}+{G}_{r}\theta +{B}_{r}\phi =0,\end{array}$ (17)

$-\frac{\partial p}{\partial y}=0,$ (18)

$\begin{array}{l}\frac{1}{{P}_{r}}\left(\u03f5{\left(\frac{\partial \theta}{\partial y}\right)}^{2}+\left(1+\u03f5\frac{{\partial}^{2}\theta}{\partial {y}^{2}}\right)\right)+{E}_{c}\left[\frac{{\partial}^{2}\psi}{\partial {y}^{2}}+\frac{n-1}{2}{W}_{e}^{2}{\left(\frac{{\partial}^{2}\psi}{\partial {y}^{2}}\right)}^{4}\right]\\ +{E}_{c}\frac{{M}^{2}}{1+{m}^{2}}{\left(\frac{\partial \psi}{\partial y}+1\right)}^{2}+\frac{{N}_{b}}{{P}_{r}}\frac{\partial \theta}{\partial y}\frac{\partial \phi}{\partial y}+\frac{{N}_{t}}{{P}_{r}}{\left(\frac{\partial \theta}{\partial y}\right)}^{2}=0,\end{array}$ (19)

$\frac{{\partial}^{2}\phi}{\partial {y}^{2}}+\frac{{N}_{t}}{{N}_{b}}\frac{{\partial}^{2}\theta}{\partial {y}^{2}}=0$. (20)

The eliminating of the pressure from Equations (17) and (18) gives

$\begin{array}{l}\frac{{\partial}^{2}}{\partial {y}^{2}}\left[1+\frac{n-1}{2}{W}_{e}^{2}{\left(\frac{{\partial}^{2}\psi}{\partial {y}^{2}}\right)}^{2}\right]\frac{{\partial}^{2}\psi}{\partial {y}^{2}}-\left(\frac{{M}^{2}}{1+{m}^{2}}+2{C}_{f}\right)\frac{{\partial}^{2}\psi}{\partial {y}^{2}}\\ \text{\hspace{0.05em}}-{C}_{f}\frac{{\partial}^{3}\psi}{\partial {y}^{3}}+{G}_{r}\frac{\partial \theta}{\partial y}+{B}_{r}\frac{\partial \phi}{\partial y}=0\end{array}$ (21)

The appropriate boundary conditions may be represented as:

$\psi =\frac{q}{2},\theta =0,\phi =0,\frac{\partial \psi}{\partial y}-L\frac{{\partial}^{2}\psi}{\partial {y}^{2}}=-1,\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}y={h}_{1}=1+a\mathrm{cos}2\text{\pi}x,$ (22)

$\psi =\frac{-q}{2},\theta =0,\phi =0,\frac{\partial \psi}{\partial y}=-1,\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}y={h}_{2}=-d-b\mathrm{cos}\left(2\text{\pi}x+\varphi \right),$ (23)

The dimensional time mean flow rate Q in the laboratory frame is related to q through the relation [20] [21]

$Q=q+1+d$ (24)

Now, the system of nonlinear differential Equations (19)-(21) subjected to the appropriate boundary conditions (22)-(32) must be solved, thus we use the Ms-DTM [18] [19] .

3. The Method of Solution

The above system may be solved by the MS-DTM method as follows. Using similar arguments as given by [19] , the functions $\psi \left(y\right)$, $\theta \left(y\right)$ and $\phi \left(y\right)$ are given from the following Equations.

$\begin{array}{l}\left[k+1\right]\left[k+2\right]\left[k+3\right]\left[k+4\right]\Psi \left[k+4\right]\\ +\frac{1}{2}\left(-1+n\right){W}_{e}^{2}{\displaystyle {\sum}_{r=0}^{k}\left[\begin{array}{c}\left[r+1\right]\left[k-r+1\right]\left[k-r+2\right]\left[k-r+3\right]\\ \left[k-r+4\right]\Psi \left[r+1\right]\Psi \left[k-r+2\right]\Psi \left[k-r+4\right]\end{array}\right]}\\ \text{\hspace{0.05em}}-\left(\frac{{M}^{2}}{1+{m}^{2}}+\frac{1}{s}\right){\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\Psi \left[r+1\right]\Psi \left[k-r+1\right]}\\ \text{\hspace{0.05em}}+{G}_{r}\left[k+1\right]\Theta \left[k+1\right]+{B}_{r}\left[k+1\right]\Phi \left[k+1\right]\end{array}$

$\begin{array}{l}\text{\hspace{0.05em}}+3\left(-1+n\right){W}_{e}^{2}({\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\left[k-r+2\right]\left[k-r+3\right]}\\ \text{\hspace{0.05em}}\times \Psi \left[r+1\right]\Psi \left[k-r+2\right]\Psi \left[k-r+3\right])+{\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\left[k-r+2\right]}\\ \text{\hspace{0.05em}}\times \left[k-r+3\right]\left[k-r+4\right]\Psi \left[k-r+2\right]\Psi \left[k-r+4\right]\\ \text{\hspace{0.05em}}-2{C}_{f}\left[k+1\right]\left[k+2\right]\Psi \left[k+2\right]-{C}_{f}\left[k+1\right]\left[k+2\right]\left[k+3\right]\Psi \left[k+3\right],\end{array}$ (24)

$\begin{array}{l}\frac{\u03f5}{{P}_{r}}[{\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\Theta \left[r+1\right]\Theta \left[k-r+1\right]}\begin{array}{c}\\ \end{array}\\ +\frac{\u03f5}{{P}_{r}}{\displaystyle {\sum}_{r=0}^{k}\left[\begin{array}{c}\left[r+1\right]\left[k-r+1\right]\left[k-r+2\right]\\ \Theta \left[r+1\right]\Theta \left[k-r+2\right]\end{array}\right]}\\ +\frac{{N}_{b}}{{P}_{r}}{\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\left[k-r+2\right]\Phi \left[k-r+2\right]\Theta \left[k-r+1\right]}]\\ +\frac{{N}_{t}}{{P}_{r}}{\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\Theta \left[r+1\right]\Theta \left[k-r+1\right]}\end{array}$

$\begin{array}{l}\text{\hspace{0.05em}}+{E}_{c}\frac{{M}^{2}}{1+{m}^{2}}[1+2\left[k+1\right]\Psi \left[k+1\right]+{\displaystyle {\sum}_{r=0}^{k}\left[r+1\right]\left[k-r+1\right]\Psi \left[r+1\right]\Psi \left[k-r+1\right]}\\ \text{\hspace{0.05em}}+\left[k+1\right]\left[k+2\right]\Psi \left[k+2\right]]+\frac{1}{{P}_{r}}\left[k+1\right]\left[k+2\right]\Theta \left[k+2\right]\\ \text{\hspace{0.05em}}+{E}_{c}\frac{n-1}{2}{W}_{e}^{2}{\displaystyle {\sum}_{r=0}^{k}\left[k+1\right]\left[k+2\right]\Psi \left[k+2\right]+\left[r+1\right]\left[k-r+1\right]}\\ \text{\hspace{0.05em}}\times \left[k-r+2\right]\left[k-r+3\right]\left[k-r+4\right]\Psi \left[r+1\right]\Psi \left[k-r+4\right],\end{array}$ (25)

$\left[k+1\right]\left[k+2\right]\Phi \left[k+2\right]+\frac{{N}_{t}}{{N}_{b}}\left[k+1\right]\left[k+2\right]\Theta \left[k+2\right]=0$ (26)

where $\Psi \left[k\right]$, $\Theta \left[k\right]$ and $\Phi \left[k\right]$ are the differential transformation functions of $\psi \left(y\right)$, $\theta \left(y\right)$ and $\phi \left(y\right)$ respectively. The differential transform of the associated boundary conditions are given by:

${\Psi}_{0}=-\frac{q}{2},{\Psi}_{1}-2L{\Psi}_{2}=-1,{\Theta}_{0}=1,{\Phi}_{0}=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}y={h}_{2}\left(x\right)=-d-b\mathrm{cos}\left(2\text{\pi}x+\varphi \right),$ (27)

$\begin{array}{l}{\displaystyle {\sum}_{k=0}^{n}{\Psi}_{k}{\left[{h}_{1}\left(x\right)-{h}_{2}\left(x\right)\right]}^{k}}=\frac{q}{2},\\ {\displaystyle {\sum}_{k=1}^{n}k{\Psi}_{k}{\left[{h}_{1}\left(x\right)-{h}_{2}\left(x\right)\right]}^{k-1}}-L{\displaystyle {\sum}_{k=1}^{n}k\left(k-1\right){\Psi}_{k}{\left[{h}_{1}\left(x\right)-{h}_{2}\left(x\right)\right]}^{k-1}}=-1,\end{array}$

$\begin{array}{l}{\displaystyle {\sum}_{k=0}^{n}{\Theta}_{k}{\left[{h}_{1}\left(x\right)-{h}_{2}\left(x\right)\right]}^{k}}=0,\\ {\displaystyle {\sum}_{k=0}^{n}{\Phi}_{k}{\left[{h}_{1}\left(x\right)-{h}_{2}\left(x\right)\right]}^{k}}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}y={h}_{1}\left(x\right)=1+a\mathrm{cos}\left(2\text{\pi}\right).\end{array}$ (28)

The expansions of the distributions of $\psi \left(y\right)$, $\theta \left(y\right)$ and $\phi \left(y\right)$ up to the tenth-order may be listed as follows:

$\begin{array}{c}\psi \left(y\right)=0.75-1.00429y-0.0428583{y}^{2}+0.196714{y}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+0.09676565239842302{y}^{4}+0.001193861687515796{y}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.0006033224257142588{y}^{6}-0.0022550902628023194{y}^{7}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+0.0004527902505426113{y}^{8}+0.0003671545327072134{y}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.0000736874408943607{y}^{10},\end{array}$

$\begin{array}{c}\theta \left(y\right)=0.807854y-2.5917363287244672{y}^{2}+2.2974739187465896{y}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-1.5353449243714945{y}^{4}+0.7069612997113166{y}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.27258261651934795{y}^{6}+0.11458745801237043{y}^{7}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.05626101276801457{y}^{8}+0.027201065873438676{y}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.0124587997837748{y}^{10},\end{array}$

And

$\begin{array}{c}\phi \left(y\right)=-0.00870575y+4.535538575267818{y}^{2}-4.020579357806532{y}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+2.686853617650115{y}^{4}-1.237182274494804{y}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+0.47701957890885893{y}^{6}-0.20052805152164827{y}^{7}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+0.0984567723440255{y}^{8}-0.04760186527851768{y}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+0.0218028996216059{y}^{10}.\end{array}$

4. Graphical Results and Discussions

The solutions acquired by the Ms-DTM are displayed through the following numerical calculation. The aim of these calculations is to obtain the influences of various parameters on the distributions for each of pressure gradient and pressure rise as well as velocity, temperature and concentration as the follows:

Figures 2(a)-(e) show the influences of the different parameters on the

distribution of pressure gradient $\frac{\text{d}P}{\text{d}x}$. It is found that the pressure gradient

increases with the increase of $M,{C}_{f}$ and ${N}_{b}$. Therefore larger pressure gradient is needed for the motivation the flow through the channel. It is a usual sense. Meanwhile, it decreases by increasing ${N}_{t}$. As seen in Figure 2(e) the pressure gradient decreases in $x\in \left[0.0,0.2989\right]$ and $x\in \left[0.544,1.0\right]$. Meanwhile it increases in $x\in \left[0.342,0.489\right]$ with the increase of ${G}_{r}$, which is named by phenomenon dual role.

Figures 3(a)-(f) show the influences of the different parameters on the distribution of pressure rise $\Delta P$ with volume flow rate Q for different values of ${N}_{t},{B}_{r},{N}_{b},L$ and ${G}_{r}$ It is seen in Figure 3(a), Figure 3(b) the pressure rise decreases in all pumping regions with an increase in ${N}_{t}$ and ${B}_{r}$. Meanwhile the converse behavior was observed with the increment in ${N}_{b}$ through Figure 3(c). Figure 3(d), Figure 3(e) show that pressure rise increases in the retrograde pumping $\left(\Delta P>0,Q<0\right)$, free pumping $\left(Q={Q}_{0}>0;\Delta P=0\right)$ and peristaltic pumping $\left(\Delta P>0,Q>0\right)$ regions with an increase in L and ${C}_{f}$. Meanwhile in the augmented pumping region $\left(\Delta P<0,Q>0\right)$ the pressure rise increases

Figure 2. (a)-(e). Effects of embedded parameters on the pressure gradient, (a) ${G}_{r}=1$, ${C}_{f}=0.5$, ${N}_{b}=0.7$, ${N}_{t}=0.8$; (b) ${G}_{r}=1$, ${N}_{b}=0.5$, $M=1$, ${N}_{t}=0.8$; (c) ${G}_{r}=1$, ${C}_{f}=0.5$, $M=1$, ${N}_{t}=0.8$; (d) ${G}_{r}=1$, ${C}_{f}=0.5$, $M=1$, ${N}_{b}=0.8$; (e) ${N}_{b}=0.7$, ${C}_{f}=0.5$, $M=1$, ${N}_{t}=0.8$, the other parameters are chosen as $x=0$, ${B}_{r}=1$, ${P}_{r}=2$, $q=1$, $\vartheta =\frac{\text{\pi}}{4}$, $a=0.4$, $b=0.5$, $d=1.5$, ${W}_{e}=0.5$, $n=1$, ${E}_{c}=0.2$, $m=1$, $\u03f5=0.1$, $L=0.1$, $\alpha =0.2$.

Figure 3. (a)-(f). Effects of embedded parameters on the pressure rise, (a) ${B}_{r}=1$, ${G}_{r}=1$, ${C}_{f}=0.5$, $L=0.1$, ${N}_{b}=0.7$; (b) ${G}_{r}=1$, ${C}_{f}=0.5$, $L=0.1$, ${N}_{b}=0.7$, ${N}_{t}=0.8$; (c) ${B}_{r}=1$, ${G}_{r}=1$, $L=0.1$, ${C}_{f}=0.7$, ${N}_{t}=0.8$; (d) ${B}_{r}=1$, ${G}_{r}=1$, $L=0.1$, ${N}_{b}=0.7$, ${N}_{t}=0.8$; (e) ${B}_{r}=1$, ${G}_{r}=1$, ${C}_{f}=0.5$, ${N}_{b}=0.7$, ${N}_{t}=0.8$; (f) ${B}_{r}=1$, ${C}_{f}=0.5$, $L=0.1$, ${N}_{b}=0.7$, ${N}_{t}=0.8$, the other parameters are chosen as $x=0$, ${P}_{r}=2$, $q=1$, $\vartheta =\frac{\text{\pi}}{4}$, $a=0.4$, $b=0.5$, $d=1.5$, ${W}_{e}=0.5$, $n=1$, ${E}_{c}=0.2$, $m=1$, $\u03f5=0.1$, $\alpha =0.2$.

till a certain value Q then it decreases (case of increasing in ${C}_{f}$ ) , also in case of increasing in L the pressure rise increases till a certain value ${Q}^{*}$ then the pressure rise decreases in the augmented pumping region. Figure 3(f) portrayed that the pressure rise decreases in all pumping regions till a certain value ${Q}^{*}$ were lies in the augmented pumping region then the pressure rise increases.

Figures 4(a)-(e) show the influences of the different parameters on the distribution of velocity profile. It is seen in Figure 4(a), Figure 4(b) the velocity profile increases at the central of the channel in $y\in \left[1.201,0.21\right]$. Meanwhile it decreases at the left and right hand sides of the channel in $y\in \left[0.22,1.28\right]$ and $y\in \left[1.85,1.31\right]$ with the increase of in M and ${C}_{f}$. Therefore, these Figs show that the gold nanoparticles are active in drug carrying and drug delivery systems, because they can encapsulate large quantities of therapeutic molecules. Figure 4(c), Figure 4(d) depicted that the velocity increases in the region $y\in \left[1.8769,0.6923\right]$, Meanwhile it decreases at the right hand side in a region $y\in \left[0.72,0.764\right]$ of the channel with the increase of ${W}_{e}$ and ${N}_{b}$. However the converse behavior can be observed through Figure 4(e) with the increase of ${N}_{t}$. This shows the dual role of the influence of the parameters on temperature.

Figures 5(a)-(f) show the influences of the different parameters on the distribution of temperature $\theta $ It is seen in Figures 5(a)-(d) the temperature distribution increases with the increase of ${W}_{e},{N}_{b},{N}_{t}$ and M. Also it is noted that increases of the temperature distribution in the case of gold. Nanoparticles can be of help the treatment of cancer. Figure 5(e), Figure 5(f) portrayed that by increasing variable thermal conductivity parameter $\u03f5$ and Forchheimer number ${C}_{f}$ the temperature distribution decreases.

Figures 6(a)-(d) show the influences of the different parameters on the distribution of nanoparticles concentration. It is seen in Figures 6(a)-(c) the nanoparticles concentration decreases with the increase of ${W}_{e},M$ and ${N}_{t}$. Figure 6(d) portrayed that the nanoparticles concentration decreases with the increase of ${N}_{b}$.

5. Conclusions

In this investigation, which behavior like a mathematical model describing the heat and mass transfer of blood flow as a Carreau nanofluid with the peristaltic flow conveying gold nanoparticles in an asymmetric channel. The coupled nonlinear governing equations have been modeled and then simplified using long wavelength approximations. The governing equations of motion are analytically solved through Ms-DTM to obtain the distribution of velocity, temperature and concentration. Ms-DTM is applicable to nonlinear models such as Carreau nanofluid models which is more complicated and has a higher degree of non-linearity, in a direct way without using linearization or restrictive assumptions. The main outcomes of the present study are concisely summarized as:

・ Thermophoresis and Brownian motion parameters have an opposite effects on pressure gradient and pressure rise.

Figure 4. (a)-(e). Effects of embedded parameters on the velocity profile, (a) ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (b) $M=1$, ${W}_{e}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (c) $M=1$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (d) $M=1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$; (e) $M=1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{b}=0.7$, the other parameters are chosen as $x=0$, ${P}_{r}=2$, $q=1$, $\vartheta =\frac{\text{\pi}}{4}$, $a=0.4$, $b=0.5$, $d=1.5$, ${B}_{r}=1$, $n=1$, ${E}_{c}=0.2$, $m=1$, $\u03f5=0.1$, $\alpha =0.2$, ${G}_{r}=1$, $L=0.1$.

Figure 5. (a)-(f). Effects of embedded parameters on the temperature profile, (a) $M=1$, $\u03f5=0.1$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (b) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$; (c) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{b}=0.7$; (d) $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (e) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (f) $M=1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$, the other parameters are chosen as $x=0$, ${P}_{r}=2$, $q=1$, $\vartheta =\frac{\text{\pi}}{4}$, $a=0.4$, $b=0.5$, $d=1.5$, ${B}_{r}=1$, $n=1$, ${E}_{c}=0.2$, $m=1$, $\alpha =0.2$, ${G}_{r}=1$, $L=0.1$.

Figure 6. (a)-(f). Effects of embedded parameters on the concentration profile, (a) $M=1$, $\u03f5=0.1$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (b) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$; (c) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{b}=0.7$; (d) $\u03f5=0.1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (e) $M=1$, $\u03f5=0.1$, ${W}_{e}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$; (f) $M=1$, ${W}_{e}=0.5$, ${C}_{f}=0.5$, ${N}_{t}=0.8$, ${N}_{b}=0.7$, the other parameters are chosen as $x=0$, ${P}_{r}=2$, $q=1$, $\vartheta =\frac{\text{\pi}}{4}$, $a=0.4$, $b=0.5$, $d=1.5$, ${B}_{r}=1$, $n=1$, ${E}_{c}=0.2$, $m=1$, $\alpha =0.2$, ${G}_{r}=1$, $L=0.1$.

・ The behavior of Weissenberg and Forchheimer number on velocity profile are too similar.

・ The behavior of pressure gradient with the increase of local Grash of number appears as the dual role phenomenon.

・ The gold particles produce the heat which is helpful for treatment of the tumor glands, because it possesses a large atomic number.

・ Increases in ${N}_{t}$ and ${N}_{b}$ cause to increase in the temperature profile.

・ On the nanoparticles concentration the behavior of increasing in ${N}_{t}$ and M are similar.

Excellent agreement are founded between our results and published results by Akbar et al. [8] through Table 1.

Table 1. Comparison of velocity, temperature and concentration solutions with exacting published paper Akbar, et al. [8] at the same value of parameter and by taking ${C}_{f}=0$, $\u03f5=0$, $L=0$ and $M=0$.

References

[1] Khan, A.K., Rashid, R., Murtaza, G. and Zahra, A. (2014) Gold Nanoparticles: Synthesis and Applications in Drug Delivery. Tropical Journal of Pharmaceutical Research, 13, 1169-1177.

https://doi.org/10.4314/tjpr.v13i7.23

[2] Huang, X. and El-Sayed, M.A. (2010) Gold Nanoparticles: Optical Properties and Implementations in Cancer Diagnosis and Photothermal Therapy. Journal of Advanced Research, 1, 13-28.

https://doi.org/10.1016/j.jare.2010.02.002

[3] Hamzehnezhad, A., Fakour, M., Ganji, D.D. and Rahbari, A. (2017) Heat Transfer and Fluid Flow of Blood Flow Containing Nanoparticles through Porous Blood Vessels with Magnetic Field. Mathematical Biosciences, 283, 38-47.

https://doi.org/10.1016/j.mbs.2016.11.009

[4] Hatami, M., Hatami, J. and Ganji, D.D. (2014) Computer Simulation of MHD Blood Conveying Gold Nanoparticles as a Third Grade Non-Newtonian Nanofluid in a Hollow Porous Vessel. Computer Methods and Programs in Biomedicine, 113, 632-641.

https://doi.org/10.1016/j.cmpb.2013.11.001

[5] Eldabe, N.T. and Abu Zeid, M. (2014) Magnetohydrodynamic Flow with Heat and Mass Transfer of Non-Newtonian Fluid Past a Vertical Heated Plate Embedded in Non-Darcy Porous Medium with Variable Porosity. Innovative Systems Design and Engineering, 5, 60-74.

[6] Rathod, V.P. and Sanjeevkumar, D. (2015) Effect of General Boundary Conditions on the Peristaltic Flow of a Nano Fluid in an Asymmetric Channel Embedded with a Non-Darcy Porous Medium. Journal of Chemical, Biological and Physical Sciences, 6, 72-95.

[7] Wu, Y.S., Bitao, L., Jennifer, L. and Miskimins, M. (2011) Analysis of Multiphase Non-Darcy Flow in Porous Media. Transport Porous Medium, 88, 205-223.

https://doi.org/10.1007/s11242-011-9735-8

[8] Akbar, N.S., Nadeem, S. and Khan, Z.H. (2014) Numerical Simulation of Peristaltic Flow of a Carreau Nanofluid in an Asymmetric Channel. Alexandria Engineering Journal, 53, 191-197.

https://doi.org/10.1016/j.aej.2013.10.003

[9] Akram, S. (2014) Effects of Nanofluid on Peristaltic Flow of a Carreau Fluid Model in an Inclined Magnetic Field. Heat Transfer—Asian Research, 43, 368-383.

https://doi.org/10.1002/htj.21082

[10] Eldabe, N.T.M., Abo-Seida, O.M., Abo-Seliem, A.S., ElShekhipy, A.A. and Hegazy, N. (2017) Peristaltic Transport of Magnetohydrodynamic Carreau Nanofluid with Heat and Mass Transfer inside Asymmetric Channel. American Journal of Computational Mathematics, 7, 1-20.

https://doi.org/10.4236/ajcm.2017.71001

[11] Nowar, K. (2014) Peristaltic Flow of a Nanofluid under the Effect of Hall Current and Porous Medium. Mathematical Problems in Engineering, 2014, Article ID: 389581.

https://doi.org/10.1155/2014/389581

[12] Abo-Eldahab, E., Barakat, E. and Nowar, K. (2012) Hall Currents and Heat Transfer Effects on Peristaltic Transport in a Vertical Asymmetric Channel through a Porous Medium. Mathematical Problems in Engineering, 2012, Article ID: 840203.

https://doi.org/10.1155/2012/840203

[13] Seth, G.S. and Ghosh, S.K. (1986) Combined Influence of Rotation and Hall Current in the Presence of a Transverse Magnetic Field. Indian Journal of Pure and Applied Mathematics, 17, 819-826.

[14] Ghosh, S.K. (1999) Hall Effect on Unsteady Hydromagnetic Flow in a Rotating Channel Permeated by an Inclined Magnetic Field in the Presence of an Oscillator. Czechoslovak Journal of Physics, 49, 465-472.

https://doi.org/10.1023/A:1022812824266

[15] Mekheimer, K.S. and Abd Elmaboud, Y. (2014) Simultaneous Effects of Variable Viscosity and Thermal Conductivity on Peristaltic Flow in a Vertical Asymmetric Channel. Canadian Journal of Physics, 92, 1541-1555.

https://doi.org/10.1139/cjp-2013-0465

[16] Thakur, P.M. and Hazarika, G.C. (2014) Effects of Variable Viscosity and Thermal Conductivity on Flow and Heat Transfer over a Stretching Surface with Variable Heat Flux in Micropolar Fluid in Presence Magnetic Field. International Journal of Scientific and Innovative Mathematical Research, 2, 554-566.

[17] Zhou, J.K. (1986) Differential Transformation and Its Applications for Electrical Circuits. Huazhong University Press, Wuhan, China.

[18] Erturk, V., Momani, S. and Odibat, Z. (2008) Application of Generalized Differential Transform Method to Multi-Order Fractional Differential Equations. Communication Nonlinear Science and Numerical Simulation, 13, 1642-1654.

https://doi.org/10.1016/j.cnsns.2007.02.006

[19] Odibata, Z.M., Bertelle, C., Aziz-Alaoui, M.A. and Duchamp, G.H.E. (2010) A Multi-Step Differential Transform Method and Application to Non-Chaotic or Chaotic Systems. Computers and Mathematics with Applications, 59, 1462-1472.

https://doi.org/10.1016/j.camwa.2009.11.005

[20] Bharali, A. and Borkakati, A.K. (1982) The Effect of Hall Currents on MHD Flow and Heat Transfer between Two Parallel Porous Plates. Applied Scientific Research, 39, 155-165.

https://doi.org/10.1007/BF00457017

[21] Akbar, N.S., Nadeem, S., Hayat, T. and Hendi, A.A. (2012) Peristaltic Flow of a Nanofluid with Slip Effects. Meccanica, 47, 1283-1294.

https://doi.org/10.1007/s11012-011-9512-3