**Keller-Box Analysis of Buongiorno Model with Brownian and Thermophoretic Di**ff**usion for Casson Nanofluid over an Inclined Surface**

**Khuram Rafique 1, Muhammad Imran Anwar 1,2,3, Masnita Misiran 1, Ilyas Khan 4,\*, Sayer O. Alharbi 5, Phatiphat Thounthong <sup>6</sup> and Kottakkaran Sooppy Nisar <sup>7</sup>**


Received: 10 August 2019; Accepted: 17 September 2019; Published: 5 November 2019

**Abstract:** The key objective of the study under concern is to probe the impacts of Brownian motion and thermophoresis diffusion on Casson nanofluid boundary layer flow over a nonlinear inclined stretching sheet, with the effect of convective boundaries and thermal radiations. Nonlinear ordinary differential equations are obtained from governing nonlinear partial differential equations by using compatible similarity transformations. The quantities associated with engineering aspects, such as skin friction, Sherwood number, and heat exchange along with various impacts of material factors on the momentum, temperature, and concentration, are elucidated and clarified with diagrams. The numerical solution of the present study is obtained via the Keller-box technique and in limiting sense are reduced to the published results for accuracy purpose.

**Keywords:** Keller-box technique; Casson nanofluid; MHD; Power law fluid; Convective boundaries; Radiation effect; Inclined surface

#### **1. Introduction**

Brownian motion and thermophoresis diffusions are the key notions of abnormal improvement in thermal conductivity by using binary fluids (base fluid along with nanoparticles). The influence of Brownian motion and thermophoresis is focused in the Buongiorno model. This model supports engineers and scholars through its utilization in the field of science and technology. It is also pointed out that nanoparticles occupying Brownian motion and thermophoresis effects cause improvement of thermal conductivity. The Brownian motion principle along thermophoresis particle installation supports, in manufacturing, germanium dioxide optical fibers and, in communication engineering, silicon. The impacts of Brownian motion and thermophoresis diffusion on Casson nanofluid flow on a stretching sheet were discussed by Anwar et al. [1]. Afify [2] scrutinized the Brownian movement and thermophosis impact on Casson nanofluid flow with convective boundaries. The impacts of radiations on Casson nanofluid flow with Brownian motion and thermophoresis influence were studied by Souayeh et al. [3]. Rashidi et al. [4] discussed heat exchange and particle motion by considering the discrete phase model (DPM). Bhatti et al. [5] examined electro-magnetohydrodynmaic (MHD) flow with heat exchange by incorporating the thermal radiations effect. Ellahi et al. [6] investigated a shiny thin film with metallic tactile covering nanoparticles through a rotating disk. Numerous scholars [7–12] considered the Buongiorno model to investigate flow characteristics.

The most significant concerns of a creator and craftsman in the construction of different items, in the pursuit of excellence, is the lessening of expenses and time. The role of heat and fluids in industries is undeniable. Discovering approaches towards the advancement of procedures and the quantity of energy exchange has consistently been a concern for researchers and specialists from the earlier times to the current era. The discovery of nanoparticles and the advancement in nanotechnology is viewed as a tremendous change in innovation and science. Choi [13] was the initiator of the nanofluid concept. A mixture of a base fluid (water, ethylene glycol and so on) with nano-scale particles called nanoparticles is termed as a nanofluid. Nanofluids have a higher thermal conductivity as compared to base fluids, due to which the energy exchange procedure is enhanced. The radiation effects on Casson nanofluid flow on a nonlinear slanted sheet were investigated by Ghadikolaei et al. [14]. The effect of a magnetic field on the flow of a nanofluid over an inclined sheet was studied by Suriyakumar and Devi [15]. Khan et al. [16] examined the flow of a Jeffery nanofluid over a slanted sheet. Thumma et al. [17] discussed the flow of a nanofluid over a nonlinear inclined stretching sheet. Parkash et al. [18] discussed the nanofluid flow through a channel analytically. Zeeshan et al. [19] examined the flow of titanium dioxide-water base nanofluid because of entropy generation. Shehzad et al. [20] calculated the silver-water base nanofluid flow in a porous medium because of entropy generation. Hussain et al. [21] studied multiphase flow synthesis with nano-size hafnium particles with the effect of electro-hydrodynamic effect. Ellahi et al. [22] discussed the thermally charged MHD bi-phase flow coating with non-Newtonian nanofluid along slipper walls. Recently, many scholars discussed nanofluid flow by incorporating different impacts [23–26].

Non-Newtonian fluids have gained considerable attention from scientists and engineers because of their key role in the field of industry and engineering. The study under concern has direct noteworthy use in association with non-Newtonian fluids, such as Casson fluids (honey, human blood etc.), power law fluids and nanofluids etc. Reddy [27] investigated the movement of Casson liquid over a slanted sheet. Hakeem et al. [28] investigated the inclined Lorentz force on Casson fluid flow on an extended sheet. Rawi et al. [29] studied the unsteady flow of Casson fluid through a slanted sheet. Casson fluid flow on an inclined sheet with multiple impacts was discussed by Jain and Parmar [30]. Ellahi et al. [31] discussed the two-phase Couette flow of couple stress by incorporating the magnetic field impacts. Ellahi et al. [32] investigated the blood flow of couple stress fluid with chemical reaction effects. For further detailed literature related to non-Newtonian fluid flow on inclined sheets, see [33–40].

Heat exchange due to thermal radiations has become an active area of research due to its vast range of applications in the field of nuclear power plants, missiles, satellites and in nanotechnology. Moreover, it is significant that thermal radiation is not suitable for the engineering of thermal tools with large variations in temperature [41]. The thermal radiation impact on flow and heat exchange is a key factor to design advanced energy conversion systems [42]. Recently, Ghadikolaei et al. [43] investigated the flow of Casson nanofluid on a porous inclined sheet numerically. Saidulu [44] discussed the radiation impacts on the flow of a nanofluid over an exponential inclined surface.

To the best of the authors' knowledge, no study on Casson nanofluid flow over an inclined nonlinear stretching sheet along with radiation effects and convective boundaries has been reported yet. Besides, the fact that a lot of work has already been done on non-Newtonian fluids with different geometries, but due to the growing applications of non-Newtonian fluids in the field of industry, the authors choose this study on an inclined sheet. The non-Newtonian fluid flow on an inclined sheet plays a vital role in MHD generators, gas turbines, and extrusion of plastic sheets. The numerical solution of the current problem is obtained using the Keller-box method.

#### **2. Problem Formulation**

Suppose two-dimensional incompressible Casson nanofluid flow over a nonlinear inclined stretching sheet slanted at <sup>γ</sup>, where *uw*(*x*) = *ax<sup>m</sup>* is the extending speed and *<sup>u</sup>*∞(*x*) = 0 is free stream speed, in which *x* is the coordinate stated towards the extending sheet and '*a*' is considered as constant. The transverse magnetic field '*B*0' is taken as normal to the track of flow. The Brownian motion and thermophoresis effects are considered. The temperature *T* and nanoparticle fraction *C* take the values *Tw* and *Cw* at the wall. The thermal radiation impact is incorporated with a convective heating procedure considered by the temperature *Tf* and heat exchange factor *hf* , which is proportional to *x*<sup>−</sup>1. Meanwhile, the encompassing structures for nanofluid temperature and mass divisions *T*<sup>∞</sup> and *C*<sup>∞</sup> are achieved as y keeps an eye on infinity, as displayed in Figure 1.

**Figure 1.** Physical geometry with coordinate system.

The flow equations for this study are given by

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

$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = \nu \left(1 + \frac{1}{\beta}\right)\frac{\partial^2 u}{\partial y^2} + g\left[\beta\_l (T - T\_{\text{co}}) + \beta\_c (\text{C} - \text{C}\_{\text{co}})\right] \text{cos}\gamma - \frac{\sigma B\_0^2(x)}{\rho} u\_\prime \tag{2}$$

$$
\mu \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = \left. \alpha \frac{\partial^2 T}{\partial y^2} - \frac{1}{\left(\rho c\right)\_f} \frac{\partial q\_f}{\partial y} + \tau \left[ D\_\mathcal{B} \frac{\partial \mathcal{C}}{\partial y} \frac{\partial T}{\partial y} + \frac{D\_T}{T\_\mathcal{co}} \left( \frac{\partial T}{\partial y} \right)^2 \right] \tag{3}
$$

$$
\mu \frac{\partial \mathbf{C}}{\partial \mathbf{x}} + v \frac{\partial \mathbf{C}}{\partial y} = D\_B \frac{\partial^2 \mathbf{C}}{\partial y^2} + \frac{D\_T}{T\_{\infty}} \frac{\partial^2 T}{\partial y^2}. \tag{4}
$$

Here, the Rosseland estimation (for radiation flux) is characterized as

$$q\_{l'} = -\frac{4\sigma^\*}{3k^\*} \left. \frac{\partial T^4}{\partial y} \right|\_{} \tag{5}$$

where the Stephen-Boltzmann coefficient is given by σ∗ and the mean absorption constant is represented by *k*∗ . Meanwhile, the temperature changes between the local temperature *T* and free steam *T*<sup>∞</sup> are very small, by ignoring higher order terms in the expansion of *<sup>T</sup>*<sup>4</sup> in Taylor succession about *<sup>T</sup>*<sup>∞</sup> for:

$$T^4 \cong 4T\_{\infty}^3 T - \Im T\_{\infty}^4. \tag{6}$$

By using Equations (5) and (6), the Equation (3) is converted into

$$
\mu \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = \left( \alpha + \frac{16 \sigma^\* T\_{\infty}^3}{3k^\*(\delta c)\_f} \right) \frac{\partial^2 T}{\partial y^2} + \tau \left[ D\_B \frac{\partial C}{\partial y} \frac{\partial T}{\partial y} + \frac{D\_T}{T\_{\infty}} \left( \frac{\partial T}{\partial y} \right)^2 \right]. \tag{7}
$$

where in the directions *x* and *y*, the velocity constituents are *u* and *v*, individually, *g* is the gravitational acceleration, the strength of the magnetic field is defined by *B*0, σ is the electrical conductivity, viscosity is given by μ, the density of conventional fluid is given by ρ*<sup>f</sup>* , the density of the nanoparticle is given by ρ*p*, β is the Casson parameter, the thermal expansion factor is denoted by β*t*, the concentration expansion constant is given by β*c*, *DB* denotes the Brownian dissemination factor and *DT* represents the thermophoresis dispersion factor. The thermal conductivity is given by *k*, the heat capacity of the nanoparticles symbolically is given as (ρ*c*)*p*, the heat capacity of the conventional liquid is given by (ρ*c*)*<sup>f</sup>* , <sup>α</sup> = *<sup>k</sup>* (ρ*c*)*<sup>f</sup>* denotes the thermal diffusivity parameter, and the symbolic representation of the relation among the current heat capacity of the nanoparticle and the liquid is <sup>τ</sup> <sup>=</sup> (ρ*c*)*<sup>p</sup>* .

In this problem, boundary conditions are considered as

$$\begin{aligned} \mu = \mu\_{\overline{\boldsymbol{w}}}(\mathbf{x}) &= a\mathbf{x}^{\mathfrak{m}}, \; \boldsymbol{\upsilon} = \boldsymbol{0} \; \text{ } -k\frac{\partial \boldsymbol{T}}{\partial \boldsymbol{y}} = h\_{f}(\boldsymbol{T}\_{f} - \boldsymbol{T}) \; \text{ } \mathbb{C} = \mathbb{C}\_{\mathfrak{w}} \; \text{at } \boldsymbol{y} = \boldsymbol{0}, \\\ \boldsymbol{u} &\to \boldsymbol{u}\_{\boldsymbol{\alpha}}(\boldsymbol{x}) = \boldsymbol{0}, \; \boldsymbol{v} \to \boldsymbol{0} \; \text{ } \boldsymbol{T} &\to \boldsymbol{T}\_{\boldsymbol{\alpha}\boldsymbol{\alpha}} \; \text{ } \mathbb{C} \to \mathbb{C}\_{\boldsymbol{\alpha}\boldsymbol{\alpha}} \; \text{at } \boldsymbol{y} \to \boldsymbol{\infty} \end{aligned} \tag{8}$$

For the conversion of the Equations (2), (4) and (7) into ordinary differential equations, we use ψ = ψ(*x*, *y*), called the stream function, characterized as

$$u = \frac{\partial \psi}{\partial y'}, \ v = -\frac{\partial \psi}{\partial x}. \tag{9}$$

(ρ*c*)*<sup>f</sup>*

The similarity transformations are considered as

$$\psi = \sqrt{\frac{2\nu\alpha x^{m+1}}{m+1}} f(\eta), \; \theta(\eta) = \frac{T - T\_{\infty}}{T\_w - T\_{\infty}} \phi(\eta) = \frac{\mathbb{C} - \mathbb{C}\_{\infty}}{\mathbb{C}\_w - \mathbb{C}\_{\infty}}, \; \eta = y \sqrt{\frac{(m+1)\alpha x^{m-1}}{2\nu}} \tag{10}$$

By using Equations (9) and (10), Equation (1) is fulfilled indistinguishably. Besides, Equations (2), (4) and (7) are transformed to the following

$$\left(1+\frac{1}{\beta}\right)f'''' + ff'' - \left(\frac{2m}{m+1}\right)f'^2 + \frac{2}{m+1}(\lambda\theta + \delta\phi)\cos\gamma - \left(\frac{2M}{m+1}\right)f' = 0\tag{11}$$

$$Pr\_N\theta'^\prime + f\theta' + Nb\phi'\theta' + Nt\theta'^2 = 0\tag{12}$$

$$
\phi'' + L\varepsilon f \phi' + Nt\_b \theta'' = 0 \tag{13}
$$

where

$$\begin{array}{c} \lambda = \frac{Gr\_{\rm x}}{Rc^{2}}, \,\delta = \frac{Gr\_{\rm l}}{Rc^{2}}, \,\mathcal{M} = \frac{a\mathbb{B}\_{0}^{-2}\langle x\rangle}{a\rho}, \,\, L\boldsymbol{x} = \frac{\nu}{D\_{b}}\mathcal{P}\boldsymbol{r} = \frac{\nu}{a}\,\,\mathcal{N}\_{b} = \frac{\tau\mathcal{D}\_{B}(\mathbb{C}\_{w}-\mathbb{C}\_{w})}{\nu},\\ \mathcal{N}\_{b} = \frac{\tau\mathcal{D}\_{l}(\mathbb{T}\_{w}-\mathbb{T}\_{w})}{\nu\mathcal{T}\_{w}}, \,\,\mathcal{N}t\_{b} = \frac{\mathcal{N}\_{b}}{\mathcal{N}\_{b}}, \,\,\mathcal{G}\boldsymbol{r}\_{\rm x} = \frac{\mathcal{G}^{\otimes}(\mathbb{T}\_{w}-\mathbb{T}\_{w})\mathcal{X}^{\otimes}}{\nu^{2}}, \,\,\mathcal{R}\boldsymbol{c}\_{\rm x} = \frac{\mu\_{\rm X}\mathcal{X}}{\nu}, \,\,\mathcal{G}\boldsymbol{c}\_{\rm x} = \frac{\mathcal{G}\mathbb{B}(\mathbb{C}\_{w}-\mathbb{C}\_{w})\mathcal{X}^{\otimes}}{\nu^{2}},\\ \mathcal{P}\boldsymbol{r}\_{\rm N} = \frac{1}{\mathcal{N}}(\mathbb{1} + \frac{\mathcal{A}}{\mathcal{N}}\mathcal{N}), \,\,\mathcal{N} = \frac{4\nu^{\prime}\mathcal{T}\_{w}^{3}}{ab^{\prime}}. \end{array} \tag{14}$$

Here, primes signify the differentiation concerning η, λ is the buoyancy parameter, δ is the solutal buoyancy parameter, the magnetic constraint is given by *M*, ν denotes the kinematic viscosity of the liquid, the Prandtl number is given as *Pr*, the Lewis number is given by *Le*, *Nb* denotes the Brownian motion parameter, *Nt* indicates the thermophoresis factor and *N* denotes the radiation factor.

The resultant boundary settings are

$$\begin{array}{c} f(\eta) = 0, \ f'(\eta) = 1, \ \theta'(0) = -\gamma\_1(1 - \theta(0)), \ \phi(\eta) = 1 \ at \ \eta = 0, \\\ f'(\eta) \to 0, \ \theta(\eta) \to 0, \ \phi(\eta) \to 0 \ as \ \eta \to \infty. \end{array} \tag{15}$$

Here, γ<sup>1</sup> = *<sup>n</sup> k* √ *Rex* is the convective parameter termed as Biot number. The skin friction, Sherwood number and Nusselt number for the current study are regarded as

$$\text{Nu}\_{\text{x}} = \frac{\text{x}q\_{\text{w}}}{k(T\_{\text{w}} - T\_{\text{ov}})}, \text{ Sh}\_{\text{x}} = \frac{\text{x}q\_{\text{m}}}{D\_{\text{B}}(\text{C}\_{\text{w}} - \text{C}\_{\text{ov}})}, \text{ C}\_{f} = \frac{t\_{\text{w}}}{\mu\_{\text{w}}\,^{2}\rho\_{f}},\tag{16}$$

where

$$q\_{\rm uv} = - \left[ k + \frac{4\sigma^\* T\_{\rm ox}^3}{3k^\*} \right] \frac{\partial T}{\partial y'}, \quad q\_{\rm uv} = -D\_B \frac{\partial \mathcal{C}}{\partial y'}, \quad \tau\_w = \mu \left( 1 + \frac{1}{\beta} \right) \frac{\partial u}{\partial y'}, \text{ at } y = 0.$$

The related terms of dimensionless reduced Nusselt number −θ (0), reduced Sherwood number −φ (0) and skin friction coefficient *Cf x* = 1 + <sup>1</sup> β *f* (0) are defined as

$$-\theta'(0) = \frac{N\mu\_x}{\left(1 + \frac{4}{3}N\right)\sqrt{\frac{m+1}{2}}\text{Re}\_x},\\ -\phi'(0) = \frac{S\hbar\_x}{\sqrt{\frac{m+1}{2}\text{Re}\_x}},\\ \mathbb{C}\_{f\ge} = \mathbb{C}\_f\sqrt{\frac{2}{m+1}\text{Re}\_{\ge}}\tag{17}$$

where *Rex* = *uwx* <sup>ν</sup> , is the local Reynolds number.

The converted nonlinear differential Equations (11)–(13) with the boundary settings (15) are elucidated by a Keller-box scheme consisting of the steps as finite-differences scheme, Newton's technique and block elimination process, clearly explained by Anwar et al. [7]. The Keller-box technique has been widely applied because it is the 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 study, the numerical outcomes of the converted nonlinear ordinary differential Equations (11)–(13) with boundary settings (15) are elucidated by the Keller-box method. For the numerical results of physical parameters of our concern, namely, Brownian motion denoted by *Nb*, thermophoresis given by *Nt*, magnetic factor *M*, buoyancy factor λ, solutal buoyancy constraint δ, inclination factor γ, Prandtl number *Pr*, Lewis number *Le*, radiation factor *N*, Casson fluid parameter β, Biot number γ<sup>1</sup> and parameter *m*, several figures and tables are prepared. In Table 1, in the deficiency of λ, δ, *M*, and *N*, and taking factor *m* = 1, with γ = 90◦ and β → ∞, the outcomes of −θ (0), −φ (0) (reduced Nusselt number, reduced Sherwood number) are equated with the results of Khan and Pop [45]. The magnitudes are established as brilliant settlement. The effects on −θ (0), −φ (0) and *Cf x*(0) against several values of involved physical parameters *Nb*, β, *Nt*, *M*, *N*, λ, δ, γ, *Pr*, *Le*, γ<sup>1</sup> and *m* are presented in Table 2. It is noted that −θ (0) drops when increasing the values of *Nb*, *Pr*, β, *N*, *Le*, *m*, and γ1, and it increased by enhancing the numerical values of γ, λ, δ, *M* and *Nt*. Moreover, it is perceived that −φ (0) is enhanced with larger values of *Nb*, *Pr*, *N*, *Le*, λ, and δ, and drops for bigger values of *m*, *M*, β, γ, γ<sup>1</sup> and *Nt*. Physically, by enhancing the Brownian motion impact, the thermal boundary layer thickness increases, and it effects a large amount of the fluid. Moreover, the Sherwood number increases and the Nusselt number decreases as we boost the thermophoresis effect; this is due to the fact that the thermal boundary layer turns thicker due to deeper diffusion penetration into

the fluid. On the other hand, *Cf x*(0) rises with the growing values of *Nb*, *Pr*, *Le*, β, *M*, *N*,γ, and *m*, and drops with the higher values of *Nt*, λ, δ, γ<sup>1</sup> and *Pr*.


**Table 1.** Contrast of the reduced Nusselt number −θ (0) and the reduced Sherwood number −φ (0) against γ = 90◦, γ<sup>1</sup> → ∞, β → ∞, *M*, *N*, λ, δ = 0, with *m* = 1, and *Pr* = *Le* = 10.

**Table 2.** Values of the reduced Nusselt number −θ (0), the reduced Sherwood number −φ (0) and the skin friction coefficient *Cf x*(0).


Figure 2 demonstrates the velocity profile against the magnetic effect. It is observed that the magnetic parameter produces Lorentz force, due to which the velocity of the fluid retards and the velocity profile drops for higher values of the magnetic parameter. Moreover, Figure 3 shows that the temperature contour is enhanced by upgrading the magnetic parameter, and the reason behind this is that the Lorentz force boosts the temperature. Consequently, the thickness of the boundary layer upturns with the increasing of the magnetic parameter. Besides, a different effect of magnetic field on concentration is noticed in Figure 4.

The influence of the nonlinear parameter on the velocity profile is shown in Figure 5. It is noted that the velocity field is not much pronounced in the case of a linear or nonlinear stretching sheet as compared to a uniformly moving surface. Similar behavior is shown in Figure 6 for the temperature profile. Moreover, an opposite effect is shown in Figure 7 for the concentration profile. Figure 8 represents the Casson effect suppressed the velocity of the fluid. It is meaningful because β reduces the yield stress in the Casson fluid. Physically, an enhancement in the Casson parameter tends to reduce the yield stress, which implies that the plastic dynamic viscosity of the liquid is enhanced and the momentum boundary layer becomes thicker [46].

**Figure 2.** Variations in velocity profile for several values of *M*.

**Figure 3.** Variations in temperature profile for several values of *M*.

Figure 9 indicates that the buoyancy force parameter λ has a directly proportional relation with the velocity profile. Physically, an increment in the buoyancy force causes a reduction of the viscous force, due to which the fluid particles move faster. In summary, the enhancement in buoyancy force tends to enhance the velocity profile. Figure 10 reveals the effect of solutal buoyancy impact on the velocity profile. The concentration difference, length and viscosity of the fluid affected the solutal buoyancy parameter. Therefore, as we enhance the solutal buoyancy parameter, the viscosity declines and the concentration increases, due to which the velocity of the fluid increases [47]. Figure 11 reflects the impact of inclination factor γ on the velocity profile. It is observed in Figure 11 that the velocity contour runs down by improving the values of γ. Moreover, the conditions specify that the maximum gravitational force acts on the flow in the case of γ = 0, because in this state the sheet will be vertical. On the other hand, for γ = 90◦ , the sheet will be horizontal, which causes a drop in velocity profile as the power of the bouncy forces drop. Figure 12 represents the impact of radiations

on the temperature profile. It reveals that the temperature profile increases with large values of thermal radiation parameter; the reason behind this is that the heat exchange is enhanced and the boundary layer thickness declines [48]. Figure 13 indicates that the temperature profile is enhanced near the boundary layer by improving the values of the Biot number. A similar behavior in the case of concentration outline against higher values of the Biot number is seen in Figure 14.

**Figure 4.** Variations in concentration profile for several values of *M*.

**Figure 5.** Variations in velocity profile for several values of *m*.

**Figure 6.** Variations in temperature profile for several values of *m*.

**Figure 7.** Variations in concentration profile for several values of *m*.

**Figure 8.** Variations in velocity profile for several values of β.

**Figure 9.** Variations in velocity profile for several values of λ.

**Figure 10.** Variations in velocity profile for several values of δ.

**Figure 11.** Variations in velocity profile for several values of γ.

**Figure 12.** Variations in temperature profile for several values of *N*.

**Figure 13.** Variations in temperature profile for several values of γ1.

**Figure 14.** Variations in concentration profile for several values of γ1.

Figure 15 represents that the temperature profile drops by improving the parameter *Pr*. This is because the bigger values of *Pr* cause improvement in viscosity and decline in the thermal boundary layer thickness. Figure 16 shows the result of Lewis number *Le* on the concentration profile. The boundary layer viscosity reduces by enhancing the values of Lewis number *Le*.

**Figure 15.** Variations in temperature profile for several values of *Pr*.

**Figure 16.** Variations in concentration profile for several values of *Le*.

Figures 17 and 18 show the effect of Brownian motion on the temperature and concentration profiles, respectively. The temperature profile upturns with improving *Nb*; on the other hand, concentration distribution has the opposite impact. Physically, the boundary layer heats up due to the development in Brownian motion, which accelerates the nanoparticles from the extending sheet to the stationary fluid. Therefore, the concentration of nanoparticles reduces. Figures 19 and 20 reveal the impact of *Nt* on temperature and concentration profile for altered values. It is observed that both temperature and concentration contours are directly proportional to the *Nt*.

**Figure 17.** Variations in temperature profile for several values of *Nb*.

**Figure 18.** Variations in concentration profile for several values of *Nb*.

**Figure 19.** Variations in temperature profile for several values of *Nt*.

**Figure 20.** Variations in concentration profile for several values of *Nt*.

#### **4. Conclusions**

In the article under study, we investigate the heat and mass transfer of Casson nanofluid flow over an inclined sheet, with convective boundaries and thermal radiation effects taken in account. The numerical results are elucidated with the Keller-box method. We found an excellent agreement between the current outcomes and already published results. The core findings of the problem under concern are the following:


**Author Contributions:** M.I.A. and K.R. formulated the problem. M.M.; S.O.A. and P.T. elucidated the numerical results. K.S.N. and I.K. plotted and discussed the results.

**Funding:** This research received no external funding.

**Acknowledgments:** The authors would like to thank Deanship of Scientific Research, Majmaah University for supporting this work.

**Conflicts of Interest:** The authors declare that they have no conflict of interest.

#### **References**


© 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/).

*Article*
