Open AccessArticle

Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface

Khuram Rafique

Muhammad Imran Anwar

^1,2,3,

Masnita Misiran

¹,

Ilyas Khan

^4,*

Asiful H. Seikh

⁵

El-Sayed M. Sherif

^5,6

and

Kottakkaran Sooppy Nisar

⁷

School of Quantitative Sciences, Universiti Utara Malaysia, Sintok 06010, Kedah, Malaysia

Department of Mathematics, Faculty of Science, University of Sargodha, Sargodha 40100, Pakistan

Higher Education Department (HED) Punjab, Lahore 54000, Pakistan

⁴

Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam

⁵

Center of Excellence for Research in Engineering Materials (CEREM), King Saud University, P.O. Box 800, Al-Riyadh 11421, Saudi Arabia

⁶

Electrochemistry and Corrosion Laboratory, Department of Physical Chemistry, National Research Centre, El-Behoth St. 33, Dokki, Cairo 12622, Egypt

⁷

Department of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Al-Dawaser 11991, Saudi Arabia

Author to whom correspondence should be addressed.

Symmetry 2019, 11(11), 1379; https://doi.org/10.3390/sym11111379

Submission received: 22 September 2019 / Revised: 17 October 2019 / Accepted: 26 October 2019 / Published: 6 November 2019

(This article belongs to the Special Issue Fluid Mechanics Physical Problems and Symmetry)

Download

Browse Figures

Versions Notes

Abstract

The prime aim of this paper is to probe the flow of micropolar nanofluid towards an inclined stretching surface adjacent to the stagnation region with Brownian motion and thermophoretic impacts. The chemical reaction and heat generation or absorption are also taken into account. The energy and mass transport of the micropolar nanofluid flow towards an inclined surface are discussed. The numerical solution is elucidated for the converted non-linear ordinary differential equation from the set of partial nonlinear differential equations via compatible similarity transformations. A converted system of ordinary differential equations is solved via the Keller-box scheme. The stretching velocity and external velocity are supposed to change linearly by the distance from the stagnation point. The impacts of involved parameters on the concerned physical quantities such as skin friction, Sherwood number, and energy exchange are discussed. These results are drawn through the graphs and presented in the tables. The energy and mass exchange rates show a direct relation with the stagnation point. In the same vein, skin friction diminishes with the growth of the stagnation factor. Heat and mass fluxes show an inverse correspondence with the inclination factor.

Keywords:

chemical reaction; micropolar nanofluid; MHD; stagnation point; heat generation or absorption; inclination effect

1. Introduction

During the last few years, the necessity to model the liquid which incorporates the gyrating micro-components has brought about presenting the hypothesis of micropolar liquids. Micropolar liquids are liquids which connect the atomic rotating motion and the microscopic velocity field. These liquids are comprised of rigid particles which are adjourned in a thick medium. A few instances of micropolar liquids are bubbly fluids, Ferro-liquids and animal blood et cetera [1]. A couple of examples of noteworthy uses of these fluids are organic structures, lubricant liquids and polymer solutions [2]. Numerous specialists examined the hypothesis of micropolar fluids across the world. In the light of aforementioned uses, at first, Eringen [3] demonstrated the idea of micropolar fluids. Rehman et al. [4] studied the flow of micropolar fluid towards an inclined surface. Moreover, Uddin [5] investigated the heat exchange of micropolar liquid towards a porous inclined surface. The numerical simulation of micropolar fluid flow towards an inclined surface has been reviewed by Shamshuddin and Thumma [6]. The latest study on micropolar nanofluid flow towards a slanted surface has been done by Rafique et al. [7]. For further detail on the flow of micropolar fluid flow with different geometries see [8,9,10,11,12].

Nanofluid has served in a wide range of applications for a more stable, more sustainable, and more efficient energy delivery, for example, porous materials [13], petroleum engineering, fuel cell industry, and medical treatment, etc., due to the significant increase in a number of properties of the fluid, such as heat transfer rate, permeability, chemical stability, etc., compared to conventional engineering fluids [14,15,16]. Nanofluids are a combination of nano-meter sized particles with the conventional liquids. By adding nanoparticles Choi [17] named these liquids as nanofluids. He scrutinized that these particles have better thermal conductivity than the regular liquids, which reasons the reduction pumping cost in the case of the heat exchange method many times.

Boundary layer flow generated via a steady stretching surface expanded significant curiosity in the previous couple of years because of its extensive uses in various engineering forms. Some important instances of valuable applications of moving stretching sheets are: solidifying of molten stones, sheet and paper construction, fiber and glass generation, hot rolling, common stock use in the kitchen, and wire outlines [18]. Sakiadis [19] initiated the boundary layer flow towards an enlarging surface the first time. Beginning now and in the probable future, Crane [20] researched a flow because of a sheet issuing with a consistent velocity. By virtue of different unique and mechanical uses, Crane’s effort has taken into account many investigators under different physical views and specific sheets. In the present day, some useful flows that move towards a stretching sheet with different extents are power law, exponential, oscillatory, and quadratic which are under idea in the current age [21,22,23,24,25,26,27].

For a significantly long era, several investigators have given noteworthy attention to the movement of non-Newtonian fluids near the stagnation area (where fluid velocity tends to zero) of a rigid body. It is due to its multidisciplinary uses in the submarine, that of thermal oil salvage, designing the missiles, et cetera. Heimenz [28] investigated the steady flow near the stagnation area. Additionally, the stagnation point flow of nanofluid with heat exchange by incorporating the chemical reaction was investigated by Gupta et al. [29]. Moreover, Rehman et al. [30] examined the flow of nanomaterial liquid adjacent to the stagnation point to an exponential sheet. Recently, Anwar et al. [31] investigated the flow of micropolar nanofluids in the stagnation point numerically. For different investigators who focused on the flow near the stagnation area for different geometries see [32,33,34,35].

Nevertheless, there are numerous liquid designing gadgets where the temperature distinction among the surface and the boundary layer fluid performs a significant role. The temperature difference is the reason to yield heat generation or absorption impacts, which have huge implications on heat exchange features, for instance, in the procedures where the working liquid experiences exothermic or endothermic chemical reactions, and in production of metal waste gained as a by-product from utilized atomic energy. In this manner, the investigation, including heat generation or absorption, has turned into a primary point of fascination for the specialists who are keen on researching fluid dynamics problems. The effect of heat generation or absorption on the flow of nanofluid towards a slanted surface was discussed by Vasanthakumari and Pondy [36]. Saeed et al. [37] investigated the heat generation or absorption of Casson fluid flow towards an inclined disk. For detailed literature on the heat generation or absorption effects on different geometries, see [38,39,40,41].

Motivated by the above-stated literature review, and due to the growing needs of non-Newtonian nanofluid flows in industry and engineering areas, we intended to investigate the flow of micropolar nanofluid adjacent to the stagnation section generated by the linear extending inclined sheet. To the best of our knowledge, the solution of micropolar nanofluid flow over a linear inclined stretching surface by incorporating the chemical reaction along with heat generation or absorption effects with the Keller-Box method has not been reported yet. The model under consideration is newly developed from Khan and Pop [42], and outcomes recovered from the current study are novel. This study is very useful in nuclear reactors, MHD generators, and geothermal energy.

2. Mathematical Formulation

Here we have an attention to scrutinize the chemical reaction and heat generation or absorption effects on the flow of micropolar type nanofluid adjacent to the stagnation region. Flow is produced by the linear extension in the inclined sheet with a stretching rate ‘

a

’. The inclination of the stretching sheet with the vertical direction is

ζ

. Brownian motion and thermophoretic impacts are considered. The wall temperature is

T_{w}

;

C_{w}

denotes nanoparticle concentration; and

u_{w}

dignifies the velocity on the wall, respectively. Moreover, the free-stream temperature is denoted by

T_{\infty}

; and

C_{\infty}

exhibits the nanoparticle concentration as y tends to infinity.

The flow under consideration is designated with the equations of Khan and Pop [42]:

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0,

(1)

\begin{array}{l} u \frac{\partial u}{\partial x} + v \frac{\partial v}{\partial y} = u_{\infty} \frac{d u_{\infty}}{d x} + (\frac{μ + k_{1}^{*}}{ρ}) \frac{\partial^{2} u}{\partial y^{2}} + (\frac{k_{1}^{*}}{ρ}) \frac{\partial N^{*}}{\partial y} + \frac{σ B_{0}^{2} (x)}{ρ} (u_{\infty} - u) \\ + g [β_{t} (T - T_{\infty}) + β_{c} (C - C_{\infty})] \cos ζ, \end{array}

(2)

u \frac{\partial N^{*}}{\partial x} + v \frac{\partial N^{*}}{\partial y} = (\frac{γ^{*}}{j^{*} ρ}) \frac{\partial^{2} N^{*}}{\partial y^{2}} - (\frac{k_{1}^{*}}{j^{*} ρ}) (2 N^{*} + \frac{\partial u}{\partial y}),

(3)

u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = α \frac{\partial^{2} T}{\partial y^{2}} + τ [D_{B} \frac{\partial C}{\partial y} \frac{\partial T}{\partial y} + \frac{D_{T}}{T_{\infty}} {(\frac{\partial T}{\partial y})}^{2}] + \frac{Q_{0}}{ρ c_{p}} (T - T_{\infty}),

(4)

u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial y} = D_{B} \frac{\partial^{2} C}{\partial y^{2}} + \frac{D_{T}}{T_{\infty}} \frac{\partial^{2} T}{\partial y^{2}} - R^{*} (C - C_{\infty}) .

(5)

where in the directions

x

and

y

, the velocity constituents are

u

and

v

individually,

g

is the gravitational acceleration, strength of magnetic field is defined by

B_{0}

σ

is the electrical conductivity, viscosity is given by

μ

, density of conventional fluid is given by

ρ_{f}

, density of the nanoparticle is given by

ρ_{p}

, thermal expansion factor is denoted by

β_{t}

, concentration expansion constant is given by

β_{c}

D_{B}

denotes the Brownian dissemination factor, and

D_{T}

represents the thermophoresis dispersion factor. The thermal conductivity given by

k

, the heat capacity of the nanoparticles symbolically is given as

{(ρ c)}_{p}

, heat capacity of the conventional liquid is given by

{(ρ c)}_{f}

α = \frac{k}{{(ρ c)}_{f}}

denotes thermal diffusivity parameter, and the symbolic representation of the relation among current heat capacity of the nanoparticle and the liquid is

τ = \frac{{(ρ c)}_{p}}{{(ρ c)}_{f}}

The boundary settings are:

\begin{array}{l} u = u_{w} (x) = a x, v = 0, N^{*} = - m_{0} \frac{\partial u}{\partial y}, T = T_{w}, C = C_{w}, at y = 0, \\ u \to u_{\infty} (x) = d x, v \to 0, N^{*} \to 0, T \to T_{\infty}, C \to C_{\infty} as y \to \infty . \end{array}

(6)

here, the stream function

ψ = ψ (x, y)

is demarcated as:

u = \frac{\partial ψ}{\partial y}, v = - \frac{\partial ψ}{\partial x},

(7)

where, equation of continuity in Equation (1) is fulfilled. The similarity transformations are defined as:

\begin{array}{l} u = a x f^{'} (η), v = - \sqrt{a v} f (η), η = y \sqrt{a / v}, \\ N^{*} = a x (\sqrt{a / v}) h (η), θ (η) = \frac{T - T_{\infty}}{T_{w} - T_{\infty}}, ϕ (η) = \frac{C - C_{\infty}}{C_{w} - C_{\infty}} . \end{array}

(8)

Equation (8) converted the system of Equations (2)–(5) into:

(1 + K) f^{‴} + f f^{″} - f {^{'}}^{2} + K h^{'} + (λ θ + δ ϕ) \cos ζ + γ^{2} + M (γ - f^{'}) = 0,

(9)

(1 + \frac{K}{2}) h^{″} + f h^{'} - f^{'} h - K (2 h + f^{″}) = 0,

(10)

(\frac{1}{\Pr}) θ^{″} + f θ^{'} + λ_{1} θ + N b ϕ^{'} θ^{'} + N t θ {^{'}}^{2} = 0,

(11)

ϕ^{″} + L e f ϕ^{'} + N t_{b} θ^{″} - L e R ϕ = 0 .

(12)

where,

\begin{matrix} λ = \frac{G r_{x}}{{Re}^{2}_{x}}, δ = \frac{G c_{x}}{{Re}^{2}_{x}}, M = \frac{σ B_{0}^{2} (x)}{a ρ}, γ = \frac{d}{a}, L e = \frac{ν}{D_{B}}, \Pr = \frac{ν}{α}, \\ N b = \frac{τ D_{B} (C_{w} - C_{\infty})}{ν}, N t = \frac{τ D_{T} (T_{w} - T_{\infty})}{ν T_{\infty}}, G r_{x} = \frac{g β_{t} (T_{w} - T_{\infty}) x^{3}}{ν^{2}}, \\ {Re}_{x} = \frac{u_{w} (x) x}{ν}, G c_{x} = \frac{g β_{c} (C_{w} - C_{\infty}) x^{3}}{ν^{2}}, N t_{b} = \frac{N t}{N b}, λ_{1} = \frac{Q_{0}}{a ρ c_{p}}, R = \frac{R^{*}}{a} . \end{matrix}

(13)

here, primes denote the differentiation with respect to

η

M

denotes the magnetic factor,

ν

denotes the kinematic viscosity,

\Pr

denotes the Prandtl number, Lewis number is denoted by

L e

, the material parameter is represented by

K

λ_{1}

denotes the heat generation or absorption parameter, and

R

indicates chemical reaction factor.

The corresponding boundary settings are changed to:

\begin{matrix} f (η) = 0, f^{'} (η) = 1, h (η) = 0, θ (η) = 1, ϕ (η) = 1, at η = 0, \\ f^{'} (η) \to γ, h (η) \to 0, θ (η) \to 0, ϕ (η) \to 0, at η \to \infty . \end{matrix}

(14)

The concerned physical quantities

C_{f}

, Nu and Sh (skin friction, Nusselt number, and Sherwood number) are demarcated as:

C_{f} = \frac{τ_{w}}{p u_{w}^{2}}, N u = \frac{x q_{w}}{k (T_{w} - T_{\infty})}, S h = \frac{x q_{m}}{D_{B} (C_{w} - C_{\infty})},

(15)

where, at y = 0 mass and heat fluxes, as well as the shear stress are expressed by

q_{m} = - D_{B} \frac{\partial C}{\partial y}

q_{w} = - k \frac{\partial T}{\partial y}

, and

τ_{w} = (μ + k_{1}^{*}) \frac{\partial u}{\partial y} + k_{1}^{*} N^{*}

, respectively. The associated expressions of the skin friction coefficient i.e.,

C_{f x} (0) = (1 + K) f^{″} (0)

- θ^{'} (0)

(reduced Nusselt number), and

- ϕ^{'} (0)

(reduced Sherwood number) are defined as:

- θ^{'} (0) = \frac{N u_{x}}{\sqrt{{Re}_{x}}}, - ϕ^{'} (0) = \frac{S h_{x}}{\sqrt{{Re}_{x}}}, C_{f x} = C_{f} \sqrt{{Re}_{x}},

(16)

where

{Re}_{x} = u_{w} (x) x / v

, is the local Reynolds number based on the stretching velocity. The converted nonlinear differential Equations (9)–(12) with the boundary settings of Equation (14) are elucidated by Keller-box scheme consisting of the steps as, finite-differences scheme, Newton’s technique, and block elimination process clearly explained by Anwar et al. [31]. The Keller-box technique has been widely applied because it is most flexible as compared to other approaches. It is informal to practice, much quicker, friendly to program, and effective.

3. Results and Discussion

In this part of the article, the numerical outcomes of physical parameters of our concern containing Brownian motion factor

N b

, thermophoresis assumed by Nt, magnetic factor

M

, buoyancy factor

λ

, solutal buoyancy constraint

δ

, inclination factor

ζ

, Prandtl number i.e.,

\Pr

, chemical reaction effect

R

, heat generation or absorption factor

λ_{1}

, stagnation parameter

γ

, Lewis number

L e

, and material factor

K

are arranged in different graphs and tables. Table 1 is supposed to be presented as a verification of the proposed model in the deficiency of

γ

λ, δ, λ_{1}, R, M, K

, and taking factor

\Pr = L e = 10

with

ζ = 90^{0}

. The magnitudes are established a brilliant settlement. The effects on

- θ^{'} (0)

- ϕ^{'} (0)

and C_fx(0) against changed magnitudes of concerned physical factors are portrayed in Table 2. From Table 2, it can be found that

- θ^{'} (0)

diminishes with the growth of

N b, N t, \Pr, L e, M, R, λ_{1}

and

ζ

, whereas it boosts on improving the magnitudes of

K, λ, δ

and

γ

. Moreover,

- ϕ^{'} (0)

increases on the growing of

N b, N t, \Pr, L e, K, λ_{1}, λ, δ, R

and

γ

, while it decreases against the cumulative magnitudes of

M

and

ζ

. Physically, by increasing the Brownian motion effect, the thermal boundary layer thickness increases and it affects a large amount of the liquid. Besides, the Sherwood number increases and the Nusselt number declines as we increase the thermophoresis effect. It is due to the fact that the thermal boundary layer turns thicker due to deeper dispersal penetration into the liquid. In sum,

C_{f x} (0)

improves on the growth of

N b, L e, M, K, R, λ_{1}

and

ζ

. On the other hand,

C_{f x} (0)

decreases with the increment in

N t, \Pr, λ, δ

and

γ

. It is observed that the skin friction reduces on improving the values of stagnation parameter, and its negative values identify the presence of drag force (employs the stretching sheet) on the motion of the micropolar nanofluid. It is not shocking, because the boundary layer is developed due to the stretching.

3.1. Velocity Profile

The behavior of the velocity field corresponding to the involved impacts are portrayed in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. For the magnetic field effect on the velocity profile for

γ < 1

and

γ > 1

see Figure 1. It demonstrates that

f^{'} (η)

(dimensionless velocity profile) reduces for

γ < 1

on improving the magnetic field strength and increases in the case of

γ > 1

. Moreover,

f^{'} (η)

enhances with the growth of

γ

for both

γ < 1

and

γ > 1

shown in Figure 2. The reason behind this is when

γ > 1

i.e., the free stream velocity greater than the stretching velocity, the flow develops a boundary layer. Physically, the fluid motion enhances near the stagnation point due to which the acceleration of the external stream increases, and in return, the boundary layer thickness declines with the enhancement of

γ

. On the other hand, when the stream velocity is less than the stretching velocity i.e.,

γ < 1

, the reversed boundary layer develops, but no boundary layer develops with

γ = 1

because in this circumstance, both velocities are equal. Figure 3 reveals the inclination impact on velocity. It demonstrated that

f^{'} (η)

decreases on increasing the inclination. Physically, the strength of the buoyancy force reduces at

ζ = 90^{0}

, leading to the decline in the velocity outline.

The velocity outline displays a direct relation with the buoyancy force, see Figure 4. Similarly, in Figure 5,

f^{'} (η)

follows the same pattern on the growth of the solutal buoyancy parameter as the buoyancy force impact in Figure 4. The effect of the magnetic field impact on

h^{'} (η)

is presented in Figure 6. The angular velocity upturns when

γ < 1

whereas, it reduces in the case of

γ > 1

for a growing magnetic field effect. Figure 7 demonstrates the effect of the material factor on the angular velocity against

γ < 1

and

γ > 1

. It is found that the angular velocity for the altered values of

K

against

γ < 1

increases but shows an inverse correspondence in the case

γ > 1

for growing the magnitude of

K

3.2. Temperature Profile

Figure 8, Figure 9 and Figure 10 exhibit the behavior of the temperature profile via altered values of incorporated factors. For the effect of the Brownian motion factor on the temperature profile for

γ < 1

and

γ > 1

drawn, see Figure 8. The temperature profile increases against the higher values of

N b

because Brownian motion is the irregular movement of the particles, which warm up the boundary layer in return and boost the temperature of the fluid. Thermophoretic impacts on the temperature profile against

γ < 1

and

γ > 1

are presented in Figure 9. The thermophoresis effect shows a direct relation with the temperature field because the variation in reference temperature and wall temperature improved on the growth of the thermophoretic impacts. Figure 10 demonstrates the behavior of temperature profile corresponds

λ_{1}

. Temperature profile increases with the enhancement in the heat generation factor; because of this, the heat generates in the flow field, and in return, the temperature enhances within the thermal boundary layer.

3.3. Concentration Profile

Figure 11, Figure 12 and Figure 13 signify the performance of the concentration profiles with under concerned parameters. An enhancement in

N b

against

γ < 1

and

γ > 1

decline the boundary layer thickness, which leads to a reduction in concentration profile (see Figure 11). Figure 12 specifies the thermophoretic effect on

ϕ^{'} (η)

by considering

γ < 1

and

γ > 1

. It is found from the sketch portrayed against the concentration for altered values of

N t

reduced. In sum, the concentration field against altered magnitudes of

R

, i.e. chemical reaction parameter, is portrayed in Figure 13. This figure discerned that

ϕ^{'} (η)

shows a reverse relation with R. Noteworthy, with the increment in the chemical reaction factor, the boundary layer thickness, as well as concentration profile, diminishes.

4. Conclusions

Heat and mass exchange of micropolar nanofluid flow adjacent to the stagnation point with magnetic effects as well as chemical reaction and heat generation or absorption reported in this article. This article fills the gap of the flow of micropolar nanofluid generated by the linear stretching inclined sheet, which plays a vital role in the heat transfer processes in industry and cooling systems. The primary outcomes of this study are as follows:

➢: Energy and mass exchange enhance with the growth of the stagnation parameter.
➢: The chemical reaction diminishes the concentration field with higher values.
➢: The stagnation parameter shows direct correspondence with the velocity profile.
➢: The heat generation or absorption factor declines the energy transport rate, whereas it improves the mass flux rate.
➢: The velocity profile shows an opposite behavior for $γ < 1$ , and $γ > 1$ .
➢: The velocity profile shows direct relation against growing magnitudes of bouncy impacts for $γ < 1$ , and $γ > 1$ .
➢: This study can be utilized in the building envelop applications because of the heat transfer.

Author Contributions

Formulation done by K.R. and M.M. Problem solved by K.R. and M.I.A. Results computed by K.R.; I.K. and A.H.S. Results discussed by K.R.; E.-S.M.S. and K.S.N. All authors contributed equally in writing manuscript.

Funding

This research was funded by Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia.

Acknowledgments

Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations Table

$C$	Fluid Concentration	$R^{*}$	Dimensionless Reaction Rate	$λ_{1}$	Heat Generation Parameter
$C_{f}$	Skin friction coefficient	$R$	Chemical reaction parameter	${Re}_{x}$	Reynolds number
$C_{\infty}$	Ambient nanoparticle volume fraction	$L e$	Lewis number	$S h$	Sherwood number
$C_{w}$	Surface volume fraction	$N b$	Brownian motion parameter	$T$	Fluid temperature
$c_{p}$	Specific heat at constant pressure	$N t$	Thermophoretic parameter	$T_{w}$	Wall temperature
$D_{B}$	Brownian diffusion coefficient	$N u$	Nusselt number	$T_{\infty}$	Ambient temperature
$D_{T}$	Thermophoretic diffusion coefficient	$\Pr$	Prandtle number	$U$	Composite velocity
$f$	Similarity function for velocity	$Q_{0}$	Dimensionless heat generation	$U_{\infty}$	Wall velocity
$ρ c_{p}$	Volume heat capacity	$μ$	Kinematic viscosity	$ν$	Dynamic viscosity
$ϕ$	Dimensionless solid volume fraction	$w$	Condition at the wall	$\infty$	Ambient condition
$δ$	Solutal buoyancy parameter	$β_{t}$	Thermal expansion coefficient	$β_{c}$	Concentration expansion coefficient
$σ$	Electric conductivity	$γ^{*}$	Spin gradient viscosity	$k_{1}^{*}$	Vertex viscosity
$j^{*}$	Micro inertia per unit mass	$ζ$	Inclination parameter	$^{^{'}}$	Differentiation with respect to $η$
$u$	Velocity in $x$ direction	$v$	Velocity in $y$ direction	$x$	Cartesian coordinate
$θ$	Dimensionless temperature	$γ$	Velocity ratio parameter	$k$	Thermal conductivity
$ρ$	Fluid density	$λ$	Bouncy parameter	$B_{0}$	Uniform magnetic field strength
$K$	Material parameter	$η$	Similarity independent variable	$α$	Thermal diffusivity
$N^{*}$	Non-dimensional angular velocity	$g$	Gravitational acceleration	$- θ^{'} (0)$	Reduced Nusselt number
$- ϕ^{'} (0)$	Reduced Sherwood number

References

Mishra, S.; Khan, I.; Al-Mdallal, Q.; Asifa, T. Free convective micropolar fluid flow and heat transfer over a shrinking sheet with heat source. Case Stud. Therm. Eng. 2018, 11, 113–119. [Google Scholar] [CrossRef]
Salleh, M.Z.; Rohni, A.M.; Amin, N. Boundary Layer Flow Due to a Moving Flat Plate in Micropolar Fluid. J. Teknol. 2005, 43, 67–83. [Google Scholar] [CrossRef]
Eringen, A.C. Simple microfluids. Int. J. Eng. Sci. 1964, 2, 205–217. [Google Scholar] [CrossRef]
Rahman, M.M.; Aziz, A.; Al-Lawatia, M.A. Heat transfer in micropolar fluid along an inclined permeable plate with variable fluid properties. Int. J. Therm. Sci. 2010, 49, 993–1002. [Google Scholar] [CrossRef]
Uddin, M.J. Convective Flow of Micropolar Fluids along an Inclined Flat Plate with Variable Electric Conductivity and Uniform Surface Heat Flux. Daffodil Int. Univ. J. Sci. Technol. 2011, 6, 69–79. [Google Scholar] [CrossRef]
Md, S.; Thumma, T.; Shamshuddin, M. Numerical study of a dissipative micropolar fluid flow past an inclined porous plate with heat source/sink. Propuls. Power Res. 2019, 8, 56–68. [Google Scholar] [CrossRef]
Rafique, K.; Anwar, M.I.; Misiran, M. Numerical Study on Micropolar Nanofluid Flow over an Inclined Surface by Means of Keller-Box. Asian J. Probab. Stat. 2019, 4, 1–21. [Google Scholar] [CrossRef]
Ali, V.; Gul, T.; Afridi, S.; Ali, F.; Alharbi, S.O.; Khan, I. Thin Film Flow of Micropolar Fluid in a Permeable Medium. Coatings 2019, 9, 98. [Google Scholar] [CrossRef]
Matta, A.; Gajjela, N. Order of chemical reaction and convective boundary condition effects on micropolar fluid flow over a stretching sheet. AIP Adv. 2018, 8, 115212. [Google Scholar] [CrossRef] [Green Version]
Jaiswal, S.; Yadav, P.K. A micropolar-Newtonian blood flow model through a porous layered artery in the presence of a magnetic field. Phys. Fluids 2019, 31, 071901. [Google Scholar] [CrossRef]
Rauf, A.; Shehzad, S.A.; Abbas, Z.; Hayat, T. Unsteady three-dimensional MHD flow of the micropolar fluid over an oscillatory disk with Cattaneo-Christov double diffusion. Appl. Math. Mech. 2019, 40, 1471–1486. [Google Scholar] [CrossRef]
Soid, S.K.; Ishak, A.; Pop, I. MHD Stagnation-Point Flow over a Stretching/Shrinking Sheet in a Micropolar Fluid with a Slip Boundary. Sains Malays. 2018, 47, 2907–2916. [Google Scholar] [CrossRef]
Xiao, B.; Wang, W.; Zhang, X.; Long, G.; Fan, J.; Chen, H.; Deng, L. A novel fractal solution for permeability and Kozeny-Carman constant of fibrous porous media made up of solid particles and porous fibers. Powder Technol. 2019, 349, 92–98. [Google Scholar] [CrossRef]
Duncan, T.V. Applications of nanotechnology in food packaging and food safety: Barrier materials, antimicrobials and sensors. J. Colloid Interface Sci. 2011, 363, 1–24. [Google Scholar] [CrossRef]
Qu, X.; Alvarez, P.J.; Li, Q. Applications of nanotechnology in water and wastewater treatment. Water Res. 2013, 47, 3931–3946. [Google Scholar] [CrossRef]
Krajnik, P.; Pusavec, F.; Rashid, A. Nanofluids: Properties, Applications and Sustainability Aspects in Materials Processing Technologies. In Advances in Sustainable Manufacturing; Springer Science and Business Media LLC: Berlin, Germany, 2011; pp. 107–113. [Google Scholar]
Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; No. ANL/MSD/CP-84938; CONF-951135-29; Argonne National Lab.: DuPage, IL, USA, 1995. [Google Scholar]
Abbas, Z.; Hasnain, J.; Sajid, M. Effects of Slip on MHD Flow of a Dusty Fluid over a Stretching Sheet through Porous Space. J. Eng. Thermophys. 2019, 28, 84–102. [Google Scholar] [CrossRef]
Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7, 26–28. [Google Scholar] [CrossRef]
Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. ZAMP 1970, 21, 645–647. [Google Scholar] [CrossRef]
Ali, M.; Alim, M.A.; Alam, M.S. Heat transfer boundary layer flow past an inclined stretching sheet in the presence of magnetic field. Int. J. Adv. Res. Technol. 2014, 3, 34–40. [Google Scholar]
Khan, I.; Shafquatullah; Malik, M.; Hussain, A.; Khan, M. Magnetohydrodynamics Carreau nanofluid flow over an inclined convective heated stretching cylinder with Joule heating. Results Phys. 2017, 7, 4001–4012. [Google Scholar] [CrossRef]
Rashad, A.M. Unsteady nanofluid flow over an inclined stretching surface with convective boundary condition and anisotropic slip impact. Int. J. Heat Technol. 2017, 35, 82–90. [Google Scholar] [CrossRef]
Reddy, P.B.A.; Bareddy, P. Magnetohydrodynamic flow of a Casson fluid over an exponentially inclined permeable stretching surface with thermal radiation and chemical reaction. Ain Shams Eng. J. 2016, 7, 593–602. [Google Scholar] [CrossRef] [Green Version]
Samanta, A. Effect of electric field on an oscillatory film flow. Phys. Fluids 2019, 31, 034109. [Google Scholar] [CrossRef]
Ramreddy, C.; Naveen, P.; Srinivasacharya, D. Department of Mathematics National Institute of Technology Warangal-; India Quadratic Convective Flow of a Micropolar Fluid along an Inclined Plate in a Non-Darcy Porous Medium with Convective Boundary Condition. Nonlinear Eng. 2017, 6, 139–151. [Google Scholar] [CrossRef]
Hayat, T.; Asad, S.; Mustafa, M.; Alsaedi, A. Radiation Effects on the Flow of Powell-Eyring Fluid Past an Unsteady Inclined Stretching Sheet with Non-Uniform Heat Source/Sink. PLoS ONE 2014, 9, e103214. [Google Scholar] [CrossRef]
Hiemenz, K. Die Grenzschicht an einem in den gleichformigen Flussigkeitsstrom eingetauchten geraden Kreiszylinder. Dinglers Polytech. J. 1911, 326, 321–324. [Google Scholar]
Gupta, S.; Kumar, D.; Singh, J. MHD mixed convective stagnation point flow and heat transfer of an incompressible nanofluid over an inclined stretching sheet with chemical reaction and radiation. Int. J. Heat Mass Transf. 2018, 118, 378–387. [Google Scholar] [CrossRef]
Rehman, F.U.; Nadeem, S.; Rehman, H.U.; Haq, R.U. Thermophysical analysis for three-dimensional MHD stagnation-point flow of nano-material influenced by an exponential stretching surface. Results Phys. 2018, 8, 316–323. [Google Scholar] [CrossRef]
Anwar, M.I.; Shafie, S.; Hayat, T.; Shehzad, S.A.; Salleh, M.Z. Numerical study for MHD stagnation-point flow of a micropolar nanofluid towards a stretching sheet. J. Braz. Soc. Mech. Sci. Eng. 2017, 39, 89–100. [Google Scholar] [CrossRef]
Dholey, S. Unsteady separated stagnation-point flow over a permeable surface. Z. Angew. Math. Phys. 2019, 70, 10. [Google Scholar] [CrossRef]
Kamal, F.; Zaimi, K.; Ishak, A.; Pop, I. Stability Analysis of MHD Stagnation-point Flow towards a Permeable Stretching/Shrinking Sheet in a Nanofluid with Chemical Reactions Effect. Sains Malays. 2019, 48, 243–250. [Google Scholar] [CrossRef]
Ghasemian, A.; Dinarvand, S.; Adamian, A.; Sheremet, M.A. Unsteady General Three-Dimensional Stagnation Point Flow of a Maxwell/Buongiorno Non-Newtonian Nanofluid. J. Nanofluids 2019, 8, 1544–1559. [Google Scholar] [CrossRef]
Nadeem, S.; Khan, M.R.; Khan, A.U. MHD oblique stagnation point flow of nanofluid over an oscillatory stretching/shrinking sheet: Existence of dual solutions. Phys. Scr. 2019, 94, 075204. [Google Scholar] [CrossRef]
Vasanthakumari, R.; Pondy, P. Mixed convection of silver and titanium dioxide nanofluids along inclined stretching sheet in presence of MHD with heat generation and suction effect. Math. Model. Eng. Probl. 2018, 5, 123–129. [Google Scholar] [CrossRef]
Saeed, A.; Shah, Z.; Islam, S.; Jawad, M.; Ullah, A.; Gul, T.; Kumam, P. Three-Dimensional Casson Nanofluid Thin Film Flow over an Inclined Rotating Disk with the Impact of Heat Generation/Consumption and Thermal Radiation. Coatings 2019, 9, 248. [Google Scholar] [CrossRef]
Elgazery, N.S. Nanofluids flow over a permeable unsteady stretching surface with non-uniform heat source/sink in the presence of inclined magnetic field. J. Egypt. Math. Soc. 2019, 27, 9. [Google Scholar]
Mjankwi, M.A.; Masanja, V.G.; Mureithi, E.W.; James, M.N. Unsteady MHD Flow of Nanofluid with Variable Properties over a Stretching Sheet in the Presence of Thermal Radiation and Chemical Reaction. Int. J. Math. Math. Sci. 2019, 2019, 1–14. [Google Scholar] [CrossRef]
Rafique, K.; Anwar, M.I.; Misiran, M. Keller-box Study on Casson Nano Fluid Flow over a Slanted Permeable Surface with Chemical Reaction. Asian Res. J. Math. 2019, 14, 1–17. [Google Scholar] [CrossRef]
Mishra, A.; Kumar, M. Ohmic–Viscous Dissipation and Heat Generation/Absorption Effects on MHD Nanofluid Flow Over a Stretching Cylinder with Suction/Injection. In Advanced Computing and Communication Technologies; Springer: Singapore, 2019; pp. 45–55. [Google Scholar]
Khan, W.; Pop, I.; Khan, W. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53, 2477–2483. [Google Scholar] [CrossRef]

Figure 1. Velocity profile for

M

Figure 1. Velocity profile for

M

Figure 2. Velocity profile for

γ

Figure 2. Velocity profile for

γ

Figure 3. Velocity profile for

ζ

Figure 3. Velocity profile for

ζ

Figure 4. Velocity profile for

λ

Figure 4. Velocity profile for

λ

Figure 5. Velocity profile for

δ

Figure 5. Velocity profile for

δ

Figure 6. Angular velocity profile for

M

Figure 6. Angular velocity profile for

M

Figure 7. Angular velocity for

K

Figure 7. Angular velocity for

K

Figure 8. Temperature profile for

N b

Figure 8. Temperature profile for

N b

Figure 9. Temperature profile for

N t

Figure 9. Temperature profile for

N t

Figure 10. Temperature profile for

λ_{1}

Figure 10. Temperature profile for

λ_{1}

Figure 11. Concentration profile for

N b

Figure 11. Concentration profile for

N b

Figure 12. Concentration profile for

N t

Figure 12. Concentration profile for

N t

Figure 13. Concentration profile for

R

Figure 13. Concentration profile for

R

Table 1. Contrast of

- θ^{'} (0)

and

- ϕ^{'} (0)

against

M, R, λ_{1}, γ, λ, δ = 0

with

L e, \Pr = 10

and

ζ = 90^{0}

Table 1. Contrast of

- θ^{'} (0)

and

- ϕ^{'} (0)

against

M, R, λ_{1}, γ, λ, δ = 0

with

L e, \Pr = 10

and

ζ = 90^{0}

$N b$	$N t$	Khan and Pop [42]		Present Results
$N b$	$N t$	$- θ^{'} (0)$	$- ϕ^{'} (0)$	$- θ^{'} (0)$	$- ϕ^{'} (0)$
0.1	0.1	0.9524	2.1294	0.9524	2.1294
0.2	0.2	0.3654	2.5152	0.3654	2.5152
0.3	0.3	0.1355	2.6088	0.1355	2.6088
0.4	0.4	0.0495	2.6038	0.0495	2.6038
0.5	0.5	0.0179	2.5731	0.0179	2.5731

Table 2. Values of

- θ^{'} (0)

- ϕ^{'} (0)

and

C_{f x} (0)

Table 2. Values of

- θ^{'} (0)

- ϕ^{'} (0)

and

C_{f x} (0)

$N b$	$N t$	$\Pr$	$L e$	$M$	$K$	$R$	$λ_{1}$	$λ$	$δ$	$γ$	$ζ$	$- θ^{'} (0)$	$- ϕ^{'} (0)$	$C_{f x} (0)$
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.8966	2.8563	0.6967
0.5	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.0640	2.8806	0.7234
0.1	0.3	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.4931	3.6809	0.6509
0.1	0.1	10.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.8322	2.9741	0.6955
0.1	0.1	7.0	10.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.8100	4.1661	0.7506
0.1	0.1	7.0	5.0	0.5	1.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.8885	2.8524	0.8200
0.1	0.1	7.0	5.0	0.1	3.0	1.0	0.1	0.1	0.9	0.5	45⁰	0.9120	2.8642	0.9514
0.1	0.1	7.0	5.0	0.1	1.0	3.0	0.1	0.1	0.9	0.5	45⁰	0.8122	4.4944	0.7482
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.5	0.1	0.9	0.5	45⁰	0.1890	3.3436	0.7051
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.5	0.9	0.5	45⁰	0.9007	2.8579	0.6064
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	2.0	0.5	45⁰	0.9071	2.8607	0.4742
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	1.5	45⁰	1.0243	2.9339	−1.4412
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	60⁰	0.8937	2.8552	0.7570
0.1	0.1	7.0	5.0	0.1	1.0	1.0	0.1	0.1	0.9	0.5	90⁰	0.8867	2.8523	0.9033

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Rafique, K.; Anwar, M.I.; Misiran, M.; Khan, I.; Seikh, A.H.; Sherif, E.-S.M.; Nisar, K.S. Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface. Symmetry 2019, 11, 1379. https://doi.org/10.3390/sym11111379

AMA Style

Rafique K, Anwar MI, Misiran M, Khan I, Seikh AH, Sherif E-SM, Nisar KS. Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface. Symmetry. 2019; 11(11):1379. https://doi.org/10.3390/sym11111379

Chicago/Turabian Style

Rafique, Khuram, Muhammad Imran Anwar, Masnita Misiran, Ilyas Khan, Asiful H. Seikh, El-Sayed M. Sherif, and Kottakkaran Sooppy Nisar. 2019. "Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface" Symmetry 11, no. 11: 1379. https://doi.org/10.3390/sym11111379

APA Style

Rafique, K., Anwar, M. I., Misiran, M., Khan, I., Seikh, A. H., Sherif, E. -S. M., & Nisar, K. S. (2019). Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface. Symmetry, 11(11), 1379. https://doi.org/10.3390/sym11111379

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Numerical Analysis with Keller-Box Scheme for Stagnation Point Effect on Flow of Micropolar Nanofluid over an Inclined Surface

Abstract

1. Introduction

2. Mathematical Formulation

3. Results and Discussion

3.1. Velocity Profile

3.2. Temperature Profile

3.3. Concentration Profile

4. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

Abbreviations Table

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI