# Journal of Fluid Dynamics

### Article Outline

Review Article | VOLUME 3 | ISSUE 1 | DOI: 10.36959/717/659 OPEN ACCESS

## The Effect of Nanoparticles on MHD Blood Flow in Stretching Arterial Porous Vessel with the Influence of Thermal Radiation, Chemical Reaction and Heat Generation/Absorption

Binyam Zigta

• Binyam Zigta 1*
• Wachemo Univesity, College of Natural and Computational Science, Ethiopia

Zigta B (2022) The Effect of Nanoparticles on MHD Blood Flow in Stretching Arterial Porous Vessel with the Influence of Thermal Radiation, Chemical Reaction and Heat Generation/Absorption. J Fluid Dyn 3(1):16-25.

Accepted: April 28, 2022 | Published Online: April 30, 2022

## Abstract

Numerical and theoretical analysis of the effect of nanoparticles on MHD blood flow in stretching arterial porous vessel with the influence of thermal radiation, chemical reaction and heat generation/absorption has been examined. The governing non-linear partial differential equations of momentum, energy and nanoparticles concentration are converted into ordinary differential equations using similarity transformations which are solved numerically. The dimensionless governing equations are solved using Runge-Kutta-Fehlberg fourth-fifth order along with shooting method. The effect of physical parameters viz., Hartmann number, unsteadiness parameter, permeability parameter, Brownian motion parameter, thermophoresis parameter, thermal radiation parameter, heat source parameter, chemical reaction parameter and Lewis number on flow variables viz., velocity of nanofluid, temperature and nanoparticles concentration has been analyzed and discussed graphically. From the simulation study the following important results are obtained: velocity of nanofluid flow increases with an increment of unsteadiness and permeability parameter; velocity of nanofluid decreases with an increment of Hartmann number; temperature profile of the model problem increases as Brownian motion parameter, thermophoresis parameter, thermal radiation parameter and heat source parameter increases. An increment in Lewis number and chemical reaction parameter results in decrement of nanoparticles concentration of the nanofluid. As the value of thermophoresis parameter increases nanoparticles concentration of the nanofluid increases.

## Keywords

Nanoparticles, Brownian motion parameter, Thermophoresis parameter, Thermal radiation, Chemical reaction

## Introduction

Blood is considered as a multi component non-Newtonian fluid which is composed of red blood cells, white blood cells, platelets and plasma. Pulsatile blood flow in cardiovascular system employing the Navier-Stokes equation and detected an increase in the wall shear stress when the porosity of the medium has been increased was modeled in [1]. Mathematical model for pulsatile blood flow subjected to externally impose periodic body acceleration by considering blood as third grade fluid has been created in [2]. The pathological situation of blood vessel was modeled and numerical solution obtained by using finite difference in [3].

Nanotechnology becomes significant field of study and has very wide applications in biology, chemistry, physics and engineering. Nanofluid is a new class of nanotechnology based on heat transfer fluids that enhances thermal properties compared with that of classical fluid particles. Nanofluid is the interaction of nanoparticles in a fluid. Nanoparticles used in nanofluids are consisted of base fluids such as water, ethylene, glycol and oil which are conductive fluids; moreover, nanoparticles are made up of metals and nonmetals such as copper, mercury, oxides, carbides, nitrides, carbon nanotubes and graphite. Nanofluid has wide applications in convective heat transfer. Model of nanofluid flow based on the Brownian motion and thermophoresis was formulated in [4]. The boundary layer flow produced in a nanofluid due to a linearly stretching by using a convective heating boundary condition studied in [5]. Mixed convection flow of a power-law nanofluid over a non-linearly stretching surface examined in [6].

Recently, Magneto hydrodynamics has received more attention due to its applications viz., petroleum industry, space dynamics, plasma physics, nuclear science and engineering. The magneto-nanofluid is the combined effect of fluid flow and heat transfer influenced by external magnetic field. Many researchers have been studied the effect of magnetic field on nanofluid with different geometries. The influence of MHD and viscous dissipation past a non-linear stretching sheet studied in [7]. Numerical solution of MHD nanofluid across stretching sheet with thermal radiation and second order slip conditions have been studied in [8]. Heat generation/absorption and thermophoresis effects on MHD flow past an oscillatory stretching sheet presented in [9]. The joint effect of magnetic field and thermal radiation for nanofluid flow and heat transfer amidst two horizontal parallel plates by considering two component models have been studied in [10]. The double effect of magnetic field and thermal radiation past a vertical stretching sheet for 2D water based nanofluid flow has been examined in [11].

In chemical engineering there is chemical reaction between foreign mass and the fluid. This process takes place in industrial application viz., food processing, polymer production and manufacturing of ceramics. Chemical reaction, heat generation/absorption effects on MHD boundary layer flow over a permeable stretching surface have been investigated in [12]. Unsteady MHD flow and heat transfer of nanofluid over a permeable shrinking sheet with thermal radiation and chemical reaction studied in [13]. The effect of MHD on free convection oscillatory Couette flow when the temperature and concentration oscillate with time in the presence of the thermal radiation and chemical reaction has been studied in [14]. The effect of thermal radiation, chemical reaction and viscous dissipation on MHD flow has been studied in [15]. Thermal radiation, chemical reaction, viscous and Joule dissipation effects on MHD flow embedded in a porous medium investigated in [16]. Mixed convection on MHD flow with thermal radiation, chemical reaction and viscous dissipation embedded in a porous medium examined in [17]. Effect of thermal radiation and chemical reaction on MHD flow of blood in stretching permeable vessel studied in [18].

Chemical reaction effect on MHD boundary layer flow of two phase nanofluid model over an exponentially stretching with a heat generation studied in [19]. In this study the effect of nanoparticles with the influence of time dependent magnetic field and time dependent velocity has not been considered. This omission encourages studying the effect of nanoparticles on MHD blood flow in stretching arterial porous vessel with the influence of thermal radiation, chemical reaction and heat generation/absorption.

The main objective of this study focuses on the effect of nanoparticles on MHD blood flow in stretching arterial porous vessel with the influence of thermal radiation, chemical reaction and heat generation/absorption with no velocity, thermal and concentration slip in the arterial porous vessel.

## Mathematical Analysis of Nanofluid Model

Consider two dimensional unsteady, incompressible, mixed convection, viscous, electrically conducting nanofluid and permeable stretching arterial vessel with the influence of time dependent magnetic field, thermal radiation, chemical reaction and heat generation/absorption. The x axis is taken along the stretching arterial porous vessel in the direction of motion and the y axis is normal to it as shown in Figure 1. The nanofluid is moving due to permeable stretching arterial vessel caused by time dependent magnetic field applied in transverse direction of the flow. Furthermore, the fluid is considered to be a gray in color, has radiation absorbing emitting nature but non-scattering medium in the optically thick limit. Rossel and approximation is used to describe the radiative heat flux applied perpendicular to the stretching arterial porous vessel. The ambient values of temperature and concentration are given by T∞ and C∞ respectively.

Consider MHD nanofluid flow in the arterial porous vessel is influenced by time dependent magnetic field which is applied in transverse direction of nanofluid flow. The arterial porous vessel is assumed to be stretched with velocity $Uw\left(x,t\right);Tw\left(x,t\right)$ and is arterial surface temperature and nanoparticles concentration of the blood are taken along the x axis respectively; and time dependent chemical reaction parameter ${k}_{r}\left(t\right)$ takes the following form

Where a,b,c and ${k}_{o}$ are the constants such that $a>0,b\ge 0,c\ge 0$ and $t<{c}_{}^{-1}$. Let us choose $B\left(t\right)=\frac{{B}_{o}}{\sqrt{\left(1-ct\right)}},{B}_{o}$ denotes the magnetic field strength at t = 0 and ${K}_{r}\left(t\right)$ is chemical reaction rate of the fluid.

The expressions, and are valid if $ct<1$ unless $c=0$ . Furthermore, the velocity describes the stretching arterial porous vessel is stretched in the positive direction and the rate of stretching increases with time. The arterial surface temperature ${T}_{w}\left(x,t\right)$ of the stretching arterial porous vessel describes the situation in which the temperature increases if $b>0$ and decreases if $b<0.$ The nanoparticles concentration of blood describes the situation in which it increases if $b>0$ and decreases if $b<0$.

To derive the governing equations the following model assumptions are made

1. MHD nanofluid is an electrically conducting fluid. Blood can be considered as MHD nanofluid.

2. Time dependent magnetic field is applied in the direction transverse to the flow of nanofluid.

3. The magnetic Reynolds numbers are small so that the induced magnetic field of the fluid is negligible.

4. Time dependent permeability of arterial porous vessel is considered.

5. The external electric field is supposed to be zero.

6. The electric field due to polarization of charges is negligible.

7. Velocity, thermal and concentration slip are negligible.

8. Thermal radiation, chemical reaction and heat absorption/generation are taken in to consideration.

Based on model assumptions, the governing equations of conservation of mass, momentum, energy and nanoparticles concentration of equation of MHD nanofluid can be expressed as follows:

Here in Eqs. (2)-(5), u and v denotes the velocity components along x and y axis respectively. Thermal diffusivity of the nanofluid is represented as and the ratio of heat capacity of nanoparticles to the fluid is denoted by. Density of the fluid is denoted by; density of the nanoparticles is represented as ; denotes electric conductivity of the fluid; is the specific heat capacity at constant pressure; represents time dependent permeability of arterial porous vessel ; denotes effective heat capacity of the nanoparticles; represents effective heat capacity of the fluid; represents Brownian diffusion coefficient; denotes thermophoresis diffusion coefficient; is the volumetric rate of heat generation/absorption within the fluid; is the radiative heat flux; represents time dependent first order chemical reaction rate, destructive chemical reaction; generative chemical reaction and no chemical reaction T and C are temperature and nanoparticles concentration respectively.

u → 0, T → T, C → C as y → ∞

It is assumed that the boundary layer is optically thick and the Roseland approximation for radiation is valid [Cortel.R, 2008]. Hence the radiative heat flux for an optically thick boundary layer [Brewster, 1992] is given by

In (8), the parameters ${\sigma }_{}^{*}$ and ${k}_{s}$ represent the Stefan Boltzmann constant and the Roseland mean absorption coefficient, respectively.

Taking the nanoparticles flow temperature as very small, in Eq. (8) can be expanded as a linear function of ,

After neglecting higher order terms in Eq. (9), ${T}_{}^{4}$ can be expanded as

Using Eq. (10), the radiative heat flux can be expressed as

Using Eq. (11) in Eq. (4) the following equation is obtained

## Non-Dimensionalization of the Model

In order to find similarity solution of the problem the following non dimensional variables are introduced [Ludlow, Clarkson and Bassom, 2000]

The continuity equation (2) satisfies the stream function defined by

Where $\psi \left(x,y,t\right)$ the stream function is defined by which satisfies continuity equation (2).

Substitution of Eq. (13) into equations (2)-(5), the following similarity equations are obtained:

In equation (16)-(18), A is unsteadiness parameter; Nb is Brownian motion parameter; is thermophoresis parameter; is Lewis number; M is the Hartmann number; R is the radiation parameter; k3 is permeability parameter; denotes the ratio of heat capacity of nanoparticles to the fluid; Pr is the Prandtl number; is chemical reaction parameter and ω are heat generation/ absorption represented respectively by

$A=c/a;Nb=\left(\frac{{D}_{B}}{v}\right)\left(C-{C}_{\infty }\right);Nt=\left(\frac{{D}_{T}}{{T}_{\infty }v}\right)\left(T-{T}_{\infty }\right);Le=v/{D}_{B}\text{ }M={B}_{0}\sqrt{\sigma /{\rho }_{f}a;}$

$\omega =\frac{{Q}_{o}}{k}$ and

The corresponding boundary conditions (6)-(7) reduces to

$f\left(\eta \right)=0,f\text{'}\left(\eta \right)=1,\theta \left(\eta \right)=1$ and $\varphi \left(\eta \right)=1$ at $\eta =0,$      (19)

$f\text{'}\left(\eta \right)=0,\theta \left(\eta \right)=0$ and $\varphi \left(\eta \right)=0$ as      (20)

## Numerical Solution to the Problem

The set of coupled non linear ordinary differential equations (16)-(18) with boundary conditions (19)-(20) are solved numerically using Runge-Kutta-Fehlberg fourth-fifth order along with shooting method. The values of boundary conditions at and are satisfied. In order to choose the appropriate value of we start with some initial guess value and solve the boundary value problem to obtain and at. Using computer language MATLAB code for different values of and step sizes it can be observed that there is no change in velocity, temperature and concentration for values of Therefore in the present study the value ofis selected to vary from 2 to 10 depending on the physical parameters and step size are taken. The coupled non linear ordinary differential equations are third order in and second order in and which have been converted in to seven simultaneous equations. There are four initial conditions at $\eta =0$ and three boundary conditions To find the solution of the problem there are three more initial conditions at, the values of $f",\theta \text{'}$ and $\varphi \text{'}$ these conditions are determined by shooting technique. Then the solutions of the problem are solved using Runge-Kutta-Fehlberg fourth-fifth order.

First convert higher order non linear differential equations to system of first order ordinary differential by using

${y}_{1}=f,{y}_{2}=f\text{'},{y}_{3}=f",{y}_{4}=\theta ,{y}_{5}=\theta \text{'},{y}_{6}=\varphi ,{y}_{7}=\varphi \text{'}$

${y}_{2}^{\text{'}}={y}_{2}$

${y}_{2}^{\text{'}}={y}_{3}$

${y}_{3}^{\text{'}}={y}_{2}{}_{}^{2}-{y}_{1}{y}_{3}+A\left({y}_{2}+\eta /2{y}_{3}\right)+{y}_{2}\left({M}_{}^{2}+1/{k}_{3}\right)$

${y}_{4}^{\text{'}}={y}_{5}$

${y}_{5}^{\text{'}}=\mathrm{Pr}/\left(1+R\right)\left\{\Gamma \left(Nb{y}_{7}{y}_{5}+Nt{y}_{5}^{2}\right)-\omega {y}_{4}+A\left(2{y}_{4}+\eta /2{y}_{5}\right)+{y}_{2}{y}_{4}-{y}_{1}{y}_{5}\right\}$

${y}_{6}^{\text{'}}={y}_{7}$

${y}_{7}^{\text{'}}=Le\left[-Nt{y}_{5}^{2}+A\left(\eta /2{y}_{7}+2{y}_{6}\right)+\left({y}_{2}{y}_{6}-{y}_{1}{y}_{7}\right)\right]+{k}_{r}^{\text{'}}{y}_{6}$

With corresponding boundary conditions

${y}_{1}\left(\eta \right)=0,\text{ }{y}_{2}\left(\eta \right)=1,\text{ }{y}_{4}\left(\eta \right)=1,\text{ }{y}_{6}\left(\eta \right)=1,\text{ }{y}_{3}\left(\eta \right)=m$

${y}_{5}\left(\eta \right)=n,\text{ }{y}_{7}\left(\eta \right)=l$

Where and $w$ are unknown to be computed using shooting technique to solve numerically.

## Simulation Study of the Model

In this section the effect of nanoparticles on MHD blood flow in permeable stretching arterial vessel with the influence of time dependent magnetic field, thermal radiation, chemical reaction and heat generation/ absorption has been studied. The effect of physical parameters viz., Hartmann number, unsteadiness parameter, permeability parameter, Brownian motion parameter, thermophoresis parameter, thermal radiation parameter, heat source parameter, chemical reaction parameter and Lewis number on flow variables viz., velocity, temperature and nanoparticles concentration has been analyzed and discussed graphically.

The governing nonlinear partial differential equations are converted in to simulataneous set of ordinary differential equations using similarity transformations which can be solved numerically using Runge-Kutta-Fehlberg fourth-fifth order along with shooting technique.

Graphical represntations of velocity profile for different values of Hartmann number, unsteadiness parameter and permeability parameters; Temperature profile for different values of Brownian motion parameter, thermophoresis parameter, thermal radiation parameter and heat source parameter. Furthermore, concentration profile for different values of thermophoresis parameter, chemical reaction parameter and Lewis number has been presented.

Figure 2 describes the effect of Hartmann number on velocity profile of the model. The graph is plotted versus for different values of dimensionless velocity against dimensionless distance respectively. From the figure it can be observed that as the values of Hartmann number increases the velocity of nanofluid decreases this occurs due to the Lorentz force which tends to resist the fluid flow and hence reduce the velocity of nanofluid flow.

Blood containing nanoparticles are considered as non-Newtonian nanofluid. Physiologically, this can be expressed velocity of nanoparticles in the blood vessel decreases due to red blood cells which contain high concentration of hemoglobin molecules. This results in an increment of concentration of red blood cells hence internal blood viscosity increases thus velocity of blood flow decreases. This is due to attribution of magnetic field which causes red blood cells to flow in the direction parallel to the magnetic field.

Figure 3 depicts the effect of unsteadiness parameter on velocity profile of the model. The graph is plotted versus with dimensionless velocity in the vertical axis and dimensionless distance in horizontal axis respectively. For particular value of unsteadiness parameter the value of nanofluid velocity starts from zero and increases as the value of dimensionless distance increases. From the simulation results it can be deduced that as unsteadiness parameter increases the velocity flow of nanofluid increases this causes increment of momentum boundary layer. Physiologically, this can be interpreted as the flow of nanofluid with slip or no slip boundary layer of arterial vessel increases unsteadiness parameter which results in an increment of velocity of the nanofluid.

Figure 4 describes the effect of permeability parameter on velocity of nanofluid flow. The graph is drawn versus for different values of dimensionless velocity against dimensionless distance. As the values of permeability parameter increases, the velocity profile of nanofluid flow increases. Physiologically, this shows the relationship between velocity of nanoparticles on blood flow and permeability of arterial porous medium parameter i.e., volume of nanoparticles that passes through porous medium dominates the total volume. Hence, the nanoparticles move from one place to another place easily in the vessel.

Figure 5 describes the effect of Brownian motion parameter on temperature profile of the model. The graph is plotted versus denoting dimensionless temperature in the vertical axis and dimensionless distance in the horizontal axis respectively. From the graph it can be concluded that as Brownian motion parameter increases temperature profile of the model increases and consequently thermal boundary layer thickness increases. Physiologically, this can be explained as temperature increases in the arterial vessel the Brownian motion of the blood increases. This is due to the fact that Brownian motion parameter is directly proportional to Brownian diffusion coefficient associated with the nanofluid which results in an increment of mass transfer rate in the arterial vessel.

Figure 6 depicts the effect of thermophores is parameter on temperature profile of the model. The graph is plotted versus denoting dimensionless temperature in the vertical axis and dimensionless distance in the horizontal axis respectively. From the simulation study it can be observed that as thermophores is parameter increases temperature profile of the model also increases. Physiologically, this can be interpreted as temperature in the arterial vessel increases this is due to the fact that thermophores is parameter is directly proportional to the heat transfer coefficient associated with the nanofluid.

Figure 7 presents temperature profile of the model for different values of thermal radiation parameter. The graph is drawn versus denoting dimensionless temperature in the vertical axis and dimensionless distance in the horizontal axis respectively. From the simulation study it can be observed that an increment in the values of thermal radiation parameter results in an increment of temperature profile of the model. Physiologically, this can be interpreted as thermal radiation increases nanoparticles flow in arterial vessel increases thus temperature of the boundary layer increases. Hence this reduces the effect of thermal buoyancy of body force which increases temperature distribution of the nanofluid. Therefore, an increment in thermal radiation results in an increment of temperature profile of the model.

Figure 8 presents the effect of heat source parameter on temperature profile of the model. The graph is drawn versus with dimensionless temperature along the vertical axis and dimensionless distance along the horizontal axis respectively. From the simulation study it can be observed that an increment of heat source parameter results in an increment of temperature profile of the model. Heat source causes to increase thermal boundary layer and the nanofluid becomes hotter. Physiologically, this can be described heat transfer is larger along the arterial vessel this results in an increment of temperature and hence this enhances thermal radiation in the arterial vessel.

Figure 9 presents the influence of Lewis number on concentration of nanoparticles of the model. The graph is drawn versus for different values of dimensionless concentration on the vertical axis and dimensionless distance on the horizontal axis. From the graphical results it can be observed that as Lewis number increases concentration of nanoparticles decreases and hence concentration boundary layer thickness. Physiologically, this shows an increment of Lewis number results in decrement of Brownian diffusion coefficient. This decrement results in decrement of nanoparticles concentration.

Figure 10 shows concentration of nanoparticles for different values of chemical reaction parameter. The graph is drawn versus for different values of dimensionless concentration on the vertical axis and dimensionless distance on the horizontal axis. From the simulation study it can be concluded that as chemical reaction parameter increases the concentration of nanofluid decreases. Physiologically, this shows that diffusion rates of nanoparticles are changed due to the influence endothermic chemical reaction. Chemical reaction is said to be endothermic if heat is absorbed. Consequently, an increment of chemical reaction parameter results in decrement of concentration.

Figure 11 depicts variation of concentration profile of the model in response to a change in thermophoresis parameter. The graph is plotted versus for different values of dimensionless concentration on the vertical axis and dimensionless distance on the horizontal axis. From the simulation study it can be observed that as thermophoresis parameter increases concentration of nanoparticles of the model increases. Consequently concentration boundary layer thickness increases. Physiologically, this can be interpreted as nanoparticles can diffuse quickly in the arterial vessel.

## Conclusion

Theoretical and numerical analysis have been investigated on the effect of nanoparticles on MHD blood flow in porous arterial vessel with the influence of thermal radiation, chemical reaction and heat generation/ absorption. The non linear partial differential equations are converted to ordinary differential equations using similarity transformations which are solved numerically.

The following important results are obtained from the simulation study:

i. An increment in unsteadiness and permeability parameter results in an increment of velocity profile of the model problem.

ii. As Hartmann number increases velocity profile of the model decreases.

iii. As Brownian motion parameter, thermophoresis parameter, thermal radiation parameter and heat source parameter increases temperature profile of the model problem increases.

iv. An increment in Lewis number and chemical reaction parameter results in decrement of nanoparticles concentration of the nanofluid.

v. As the values of thermophoresis parameter increases nanoparticles concentration of the nanofluid increases.

## Application of the Model Problem in Medical Science

The numerical results of the present study are compared with Ferdows, et al. (2012) and Bidin and Nazar, et al. (2009) in order to check the validity of the present work. Thus the results are in good agreement with mentioned literature reviews.

Blood is considered as non-Newtonian fluid, the small sized nanoparticles have very wide application in biomedical application.

i. Nanoparticles like gold can be applied to activate or prevent growth of blood vein (Hatami M, HatamiJ,Ganji DD[2014]).

ii. To enhance or reduce blood capillary growth in some specified diseases some drugs are applied but efficient for short time. Recently researchers revealed that nanoparticles are effective in drug carrying and delivery vehicles because large quantities of therapeutic molecules are encapsulated in them (Rahbari A, Fakour M, Hamzehezhad A, Vakilabadi MA, Ganji DD [2017]).

iii. In this study the effect of nanoparticles on MHD blood flows of arterial porous vessel with the influence of thermal radiation, chemical reaction and heat generation/absorption.

## Nomenclature

velocity of the nanofluid flow along x axis of the arterial vessel

transverse velocity of nanofluid flow in arterial porous vessel

temperature of blood at any point in the arterial porous vessel

concentration of nanoparticles of the blood at any point in the arterial porous vessel

streching velocity of arterial porous vessel of blood flow

nanoparticles concentration of blood along the x axis

temperature of arterial porous vessel wall

ambient temperature of the medium

ambient nanoparticles concentration of the solute

effective heat capacity of the nanoparticles

effective heat capacity of the fluid

time dependent magnetic field intensity

Brownian diffusion coefficient

thermophoresis diffusion coefficient

Lewis number

thermophoresis parameter

Brownian motion parameter

time dependent permeability of arterial porous vessel

constant permeability of arterial porous medium

non-dimensional permeability of arterial porous medium parameter

$M-$ Hartmann number of blood

prandtl number

time dependent chemical reaction of nanoparticles of arterial blood vessel

thermal conductivity of nanoparticles of arterial blood vessel

stream function

heat generation/absorption of the arterial vessel

dimensionless temperature

the ratio of heat capacity of nanoparticles to the fluid

electrical conductivity of fluid

density of fluid

density of nanoparticles

$v-$ kinematic viscosity of blood

dynamic viscosity of blood

dimensionless similarity variable

dynamic viscosity

## Abstract

Numerical and theoretical analysis of the effect of nanoparticles on MHD blood flow in stretching arterial porous vessel with the influence of thermal radiation, chemical reaction and heat generation/absorption has been examined. The governing non-linear partial differential equations of momentum, energy and nanoparticles concentration are converted into ordinary differential equations using similarity transformations which are solved numerically. The dimensionless governing equations are solved using Runge-Kutta-Fehlberg fourth-fifth order along with shooting method. The effect of physical parameters viz., Hartmann number, unsteadiness parameter, permeability parameter, Brownian motion parameter, thermophoresis parameter, thermal radiation parameter, heat source parameter, chemical reaction parameter and Lewis number on flow variables viz., velocity of nanofluid, temperature and nanoparticles concentration has been analyzed and discussed graphically. From the simulation study the following important results are obtained: velocity of nanofluid flow increases with an increment of unsteadiness and permeability parameter; velocity of nanofluid decreases with an increment of Hartmann number; temperature profile of the model problem increases as Brownian motion parameter, thermophoresis parameter, thermal radiation parameter and heat source parameter increases. An increment in Lewis number and chemical reaction parameter results in decrement of nanoparticles concentration of the nanofluid. As the value of thermophoresis parameter increases nanoparticles concentration of the nanofluid increases.