Open AccessArticle

Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions

Abdul Samad Khan

^1,*

Yufeng Nie

¹,

Zahir Shah

Abdullah Dawar

Waris Khan

⁴

and

Saeed Islam

Department of Applied Mathematics, School of Science, Northwestern Polytechnical University, Dongxiang Road, Chang’an District, Xi’an 710129, China

Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan

Department of Mathematics, Qurtuba University of Science and Information Technology Peshawar, Peshawar 25000, Pakistan

⁴

Department of Mathematics, Kohat University of Science and technology, Kohat 26000, Pakistan

Author to whom correspondence should be addressed.

Appl. Sci. 2018, 8(11), 2244; https://doi.org/10.3390/app8112244

Submission received: 23 October 2018 / Revised: 7 November 2018 / Accepted: 10 November 2018 / Published: 14 November 2018

(This article belongs to the Section Nanotechnology and Applied Nanosciences)

Download

Browse Figures

Versions Notes

Abstract

In this study, we analyzed the three-dimensional flow of Williamson (pseudoplastic) fluids upon a linear porous stretching sheet. The thermal radiation impact was taken into account. The transformed non-linear equations were solved by the homotopy analysis method (HAM). The influence of the embedded parameters tretching parameter, Williamson parameter, porosity parameter, thermal radiation parameter, thermophoresis parameter, Brownian motion parameter, Prandtl number and Biot number are presented on velocity, temperature and concentration functions in the graphs and explained in detail. The velocity function along the x-direction reduces with the impact of the stretching, porosity and Williamson parameters. Velocity along the y-direction increases with the stretching parameter, while it reduces with the porosity and Williamson parameters. The effect of Skin friction, heat transfer and mass transfer are shown numerically. The numerical values of surface drag force and the impact of different parameters are calculated and it is observed that increasing the stretching parameter and the porosity parameter reduces the surface drag force, while increasing the Williamson parameter augments the surface drag force. Higher values of the stretching parameter, the Prandtl number and the radiation parameter enhance the heat transfer rate, while the augmented value of the thermophoresis and Brownian motion parameters reduces the heat transfer rate, where higher values of the stretching parameter, thermophoresis and Brownian motion parameters enhance the mass transfer rate.

Keywords:

nanofluid; heat transfer; mass transfer; Williamson fluid; convective boundary condition; HAM

1. Introduction

The non-Newtonian fluids have many applications in the field of fluid mechanics. For instance, geophysics, chemical industries, biological sciences and petroleum industries are some common applications in various industries. Due to these applications, research in the area of non-Newtonian fluids has received attention. Such flows appear in food processing, plastic manufacturing, polymer processing, biological fluids, ice and magma flows. Cord depiction, extrusion and copper spiraling are significant aspects of fluid flow upon a stretching surface. There are various studies that examined fluid flow over a stretching surface in two dimensions since the works of Sakadais and Crane [1,2]. In recent times, researchers drew attention to the field by studying fluid flow over a stretching surface in three dimensions (3-D). In 1929, Williamson [3] presented the theory of pseudo plastic fluids and found the practical significance of these fluids. Using the homotopy perturbation technique, the steady and laminar flow past a surface in three dimensions was presented by Ariel [4]. Using the Adomian decomposition method, the magnetohydrodynamic and boundary layer flow of Casson fluid past an exponentially shrinking sheet has been analyzed by Nadeem et al. [5]. In another article, Nadeem et al. [6] analytically probed the two-dimensional peristaltic flow of the Williamson model. The analytical solution of an incompressible micro polar fluid over an extending sheet has been analyzed by Pop et al. [7]. The analysis of viscous fluid flow in 3-D, past over extending surface, has been investigated by Sakiadis [8]. The two cases of heat transfer for recommended heat flux, and surface temperature with unsteady fluid flow over a continuous moving surface temperature, have been investigated by Tsou et al. [9]. They concluded that the velocity escalates with the escalation of the stretching parameter. Recently, premeditated three-dimensional Williamson fluid flow past a stretching surface has been studied by Malik et al. [10]. The process of blowing and suction combined with heat and mass transference rate past a stretching sheet has been introduced by Gupta and Gupta [11]. They extended the work of Erickson et al. [12] by replacing the constant surface speed with the linear surface speed. Aziz [13] perceived the heat radiation effect by studying the flow on an unsteady stretched sheet. In the presence of a porous medium, the impact of thermal radiation on a vertically stretched surface was studied by Mukhopadyay [14]. Shateyi and Motsa [15] numerically studied the changes in mass transfer and concentration rates over a horizontal stretched sheet. The effect of heat and momentum of incompressible fluid flow past a linear stretched surface has been investigated by Aziz [16]. The viscous fluid flow over a non-linear stretching surface under the impact of thermal radiation has been analyzed by Cortell [17]. In the case of uniform and non-uniform heat fluxes, the numerical study of heat transfer has been analyzed by Ishak et al. [18]. They claimed that the non-Newtonian model (Micro polar fluid) has a higher coefficient of convective heat transfer than viscous fluid flow. Eyring−Powell fluid flow under the impact of thermal radiation over a porous medium was investigated by Dawar et al. [19]. Magnetohydrodynamic (MHD) flow carbon nanotubes Casson nanofluid in rotating channels was studied by Dawar et al. [20]. Analytically-studied MHD Jeffery fluid under the impact of thermal radiation was studied by Abro et al. [21]. The recent analysis by Sheikholeslami on nanofluids can be seen in [22,23,24]. The recent relevant investigation on Casson fluids and nanofluids with different effects and their modern application can be studied in [25,26,27,28,29,30,31,32]. Bhatti et al. [33] and Darbari et al. [34] have investigated entropy generation in nanofluids. Selecting dynamic working nanofluids has become a popular subject in the recent decade [35,36,37,38,39,40,41,42,43,44].

Brownian motion is the random motion of nanoparticles in a fluid flow. It is understood that Brownian motion is the significant factor in thermal performance of nanoparticles in nanofluids [45]. The knowledge of heat transmission in nanofluids flow through diverse geometries is significant for heat exchanger design, transpiration, fiber coating etc. Currently, nanomaterials are among the well-known tackles for refining the low thermal conductivity of working liquids. Naturally, nanoparticles move randomly within the base fluid. The detailed study of heat transmission in nanofluids and their importance can be seen in [46,47,48].

Keeping in view the above mentioned literature, to the best of our knowledge, there is no work found similar to the study of three-dimensional Williamson fluid over a linear porous stretching sheet under the impact of thermal radiation. The transformed non-linear and coupled equations have been solved by the homotopy analysis method [49,50]. The impact of the stretching parameter, Williamson parameter, porosity parameter, thermal radiation parameter, thermophoresis parameter, Brownian motion parameter, Prandtl number and Biot number are presented in graphs and discussed in detail. The surface drag force, heat transfer rate and mass transfer rate are presented in tables.

2. Problem Formulation

Suppose the unsteady three-dimensional incompressible Williamson fluid flows past a linear stretching sheet. Also, it is assumed that along the

x y -direction

, the sheet is stretched. The motion of the fluid causes stretching of the sheet and the fluid occupies the space

z > 0

. The stretching sheet has been kept porous. Using the boundary layer approximation, the equations in Williamson fluid flow for continuity and momentum is given as [26,27,28,29,30]

\nabla . \vec{V} = 0,

(1)

ρ \frac{d \vec{V}}{d t} = d i v \vec{S} + ρ b .

(2)

where

ρ, V, S, b, d / d t

represent the density, velocity vector, Cauchy stress tensor, specific body force vector and material derivative, respectively.

The mathematical equations for Williamson fluid are given as [10,28]

\vec{S} = P \vec{I} + τ,

(3)

τ = [\frac{μ_{0} - μ_{\infty}}{1 - Γ γ^{*}}] {\vec{A}}_{1} .

(4)

where

τ

represents the extra stress tensor,

μ_{0}

indicates the limiting viscosity at zero,

μ_{\infty}

shows the limiting viscosity at infinite shear rate,

Γ > 0

is the constant of time,

{\vec{A}}_{1} = \nabla v + {(\nabla v)}^{T}

and

{\vec{A}}_{1}

represents the Rivlin-Erickson tensor and

γ^{*}

is given as

γ^{*} = \sqrt{\frac{1}{2} π},

(5)

where

π = t r a c e ({\vec{A}}_{1}^{2}) .

(6)

π

is the second invariant strain tensor.

The mathematical equations for the stated problem are as follows [10,28,32]

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

(7)

u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} = ν \frac{\partial^{2} u}{\partial z^{2}} + \sqrt{2} ν Γ \frac{\partial u}{\partial z} \frac{\partial^{2} u}{\partial z^{2}} - \frac{ν}{k} u,

(8)

u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} + w \frac{\partial v}{\partial z} = ν \frac{\partial^{2} v}{\partial z^{2}} + \sqrt{2} ν Γ \frac{\partial v}{\partial z} \frac{\partial^{2} v}{\partial z^{2}} - \frac{ν}{k} v,

(9)

u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} + w \frac{\partial T}{\partial z} = \frac{k}{ρ c_{p}} \frac{\partial^{2} T}{\partial z^{2}} - \frac{\partial q_{r}}{\partial z} + τ_{1} [D_{B} \frac{\partial C}{\partial z} \frac{\partial T}{\partial z} + \frac{D_{T}}{T_{\infty}} {(\frac{\partial T}{\partial z})}^{2}],

(10)

u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial y} + w \frac{\partial C}{\partial z} = D_{B} \frac{\partial^{2} C}{\partial z^{2}} + \frac{D_{T}}{T_{\infty}} (\frac{\partial^{2} T}{\partial z^{2}}) .

(11)

In Equations (7)–(9),

u

indicates the component of velocity in the x-direction,

v

is in the y-direction and

w

is in the z-direction,

Γ

represents Williamson fluid and

ν

indicates kinematic viscosity. In Equations (10) and (11),

k

denotes thermal conductivity,

T

is the temperature,

ρ

is the dynamic viscosity,

c_{p}

is the heat capacity,

τ_{1}

represents the ratio of the heat capacity,

D_{B}

is the Brownian motion coefficient,

D_{T}

indicates the thermophoretic coefficient and

q_{r}

is the radiative heat flux, which is defined as

q_{r} = - \frac{4 σ^{*}}{3 k^{*}} \frac{\partial T^{4}}{\partial z},

(12)

where

σ^{*}, k

are symbols used for the absorption coefficient and the Stefan Boltzmann constant, respectively. As

T^{4} = 4 T_{\infty}^{3} T - 3 T_{\infty}^{4}

, implementing this in Equation (10), it reduces to

u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} + w \frac{\partial T}{\partial z} = \frac{k}{ρ c_{p}} (1 + \frac{4 σ^{*} T_{\infty}^{3}}{k^{*}}) \frac{\partial^{2} T}{\partial z^{2}} + τ_{1} [D_{B} \frac{\partial C}{\partial z} \frac{\partial T}{\partial z} + \frac{D_{T}}{T_{\infty}} {(\frac{\partial T}{\partial z})}^{2}] .

(13)

The equivalent boundary conditions for the said problem are given here as

\begin{array}{l} \begin{array}{l} u = U_{w} = a x, v = V_{w} = b y \\ - k \frac{\partial T}{\partial z} = h_{f} (T_{f} - T), - D_{B} \frac{\partial C}{\partial z} = h_{s} (C_{s} - C) a t η = 0 \end{array}} \\ u \to 0, v \to 0, w \to 0, T \to T_{\infty}, C \to C_{\infty} a s η \to \infty . \end{array}

(14)

In the preceding equation,

U_{w}

and

V_{w}

are the stretching velocities,

a

and

b

are positive constants, where

c = \frac{a}{b}

h_{f}, h_{s}

are the convective heat and mass transmission coefficients,

T_{f}

is the convective fluid temperature and

C

is the concentration below the moving sheet.

In view of the above equations, the similarity variables are defined as

\begin{array}{l} u = a x f^{'} (η), v = b y g^{'} (η), w = - {(a ν)}^{\frac{1}{2}} (f (η) + c g (η)), \\ θ (η) = \frac{T - T_{\infty}}{T_{f} - T_{\infty}}, ϕ (η) = \frac{C - C_{\infty}}{C_{f} - C_{\infty}}, η = z \sqrt{\frac{a}{ν}} . \end{array}

(15)

In view of Equation (15), Equation (7) is satisfied, Equations (8), (9), (11), and (13) are reduced to

f^{‴} - {(f^{'})}^{2} + (f + c g) f^{″} + W e f^{″} f^{‴} + g f^{″} - γ f^{'} = 0,

(16)

g^{‴} - {(g^{'})}^{2} + (f + c g) g^{″} + W e g^{″} g^{‴} + f g^{″} - γ g^{'} = 0,

(17)

\frac{1}{\Pr} (1 + \frac{4}{3} R d) θ^{″} + (f + c g) θ^{'} + N b ϕ^{'} θ^{'} + N t {(θ^{'})}^{2} = 0,

(18)

ϕ^{″} (η) + L e (f + c g) ϕ^{'} (η) + (\frac{N t}{N b}) θ^{″} (η) = 0.

(19)

With boundary conditions

\begin{array}{l} f (0) = 0, f^{'} (0) = 1, g^{'} (0) = c, \\ g^{'} (0) = 0, f^{'} (\infty) = 0, g^{'} (\infty) = 0 \\ θ^{'} (0) = - B i_{1} (1 - θ (0)), ϕ^{'} (0) = - B i_{2} (1 - ϕ (0)), \\ θ (\infty) = ϕ (\infty) = 0. \end{array}

(20)

In the above equations,

W e = Γ x \sqrt{2 a^{3} / ν}

represents Williamson parameter,

γ = ν / a k

indicates the porosity parameter,

\Pr = ρ ν c_{p} / k

indicates the Prandtl number,

R d = 4 σ T_{\infty}^{3} / k k^{*}

represents the thermal radiation parameter,

N t = τ_{1} D_{T} (T_{f} - T_{\infty}) / ν T_{\infty}, N b = τ_{1} D_{T} (C_{f} - C_{\infty}) / ν

represent the thermophoresis and Brownian motion parameters.

B i_{1} = \frac{h_{f}}{k} / \sqrt{ν / a}, B i_{2} = \frac{h_{s}}{D_{B}} / \sqrt{ν / a}

, where

B i_{1}

and

B i_{2}

the Biot numbers. Skin friction, local Nusselt number and Sherwood number in dimensionless form are defined as

\begin{array}{l} \frac{1}{\sqrt{2}} C_{f} \sqrt{{Re}_{x}} = (f^{″} (0) + W e {(f^{″})}^{2} (0)), \\ \frac{1}{\sqrt{2}} C_{f} \sqrt{{Re}_{x}} = (g^{″} (0) + W e {(g^{″})}^{2} (0)), \end{array}

(21)

Ν u = - (1 + \frac{4}{3} R d) Θ^{'} (0), S h = - Φ^{'} (0) .

(22)

3. Solution by Homotopy Analysis Method

For the solution of obtained nonlinear model Equations (16)–(19) with modeled boundary Conditions (20), we have used the Homotopy Analysis Method (HAM). The preliminary suppositions are chosen as follows:

f_{0} (η) = 1 - e^{- η}, g_{0} (η) = c (1 - e^{- η}), θ_{0} (η) = (\frac{B i_{1}}{1 + B i_{1}}) e^{- η}, ϕ_{0} (η) = (\frac{B i_{2}}{1 + B i_{2}}) e^{- η} .

(23)

The

L_{f} {, L}_{g} {, L}_{θ} {and L}_{ϕ}

are linear operators, which are taken as

L_{f} (f) = f^{‴} - f^{'}, L_{g} (g) = g^{‴} - g^{'}, L_{θ} (θ) = θ^{″} - θ, L_{ϕ} (ϕ) = ϕ^{″} - ϕ,

(24)

which have the succeeding properties:

\begin{array}{l} L_{f} (λ_{1} + λ_{2} e^{- η} + λ_{3} e^{η}) = 0, L_{g} (λ_{4} + λ_{5} e^{- η} + λ_{6} e^{η}) = 0, \\ L_{θ} (λ_{7} e^{- η} + λ_{8} e^{η}) = 0, L_{ϕ} (λ_{9} e^{- η} + λ_{10} e^{η}) = 0. \end{array}

(25)

where

λ_{1}, λ_{2}, λ_{3}, λ_{4}, λ_{5}, λ_{6}, λ_{7}, λ_{8}, λ_{9}, λ_{10}

are constant coefficients.

The non-linear operators

N_{f}, N_{g}, N_{θ}

and

N_{ϕ}

are indicated as

\begin{array}{l} N_{f} [f (η; ℜ), g (η; ℜ)] = \frac{\partial^{3} f (η; ℜ)}{\partial η^{3}} - {(\frac{\partial f (η; ℜ)}{\partial η})}^{2} + (f (η; ℜ) + g (η; ℜ)) \frac{\partial^{2} f (η; ℜ)}{\partial η^{2}} \\ + (W e) \frac{\partial^{2} f (η; ℜ)}{\partial η^{2}} \frac{\partial^{3} f (η; ℜ)}{\partial η^{3}} + g (η; ℜ) \frac{\partial^{2} f (η; ℜ)}{\partial η^{2}} - (γ) \frac{\partial f (η; ℜ)}{\partial η}, \end{array}

(26)

\begin{array}{l} N_{g} [f (η; ℜ), g (η; ℜ)] = \frac{\partial^{3} g (η; ℜ)}{\partial η^{3}} - {(\frac{\partial g (η; ℜ)}{\partial η})}^{2} \\ + (f (η; ℜ) + g (η; ℜ)) \frac{\partial^{2} g (η; ℜ)}{\partial η^{2}} + (W e) \frac{\partial^{2} g (η; ℜ)}{\partial η^{2}} \frac{\partial^{3} g (η; ℜ)}{\partial η^{3}} \\ + f (η; ℜ) \frac{\partial^{2} g (η; ℜ)}{\partial η^{2}} - (γ) \frac{\partial g (η; ℜ)}{\partial η}, \end{array}

(27)

\begin{array}{l} N_{θ} [θ (η; ℜ), ϕ (η; ℜ)] = \frac{1}{\Pr} (1 + \frac{4}{3} R d) \frac{\partial^{2} θ (η; ℜ)}{\partial η^{2}} \\ + (f + c g) \frac{\partial θ (η; ℜ)}{\partial η} + (N b) \frac{\partial ϕ (η; ℜ)}{\partial η} \frac{\partial θ (η; ℜ)}{\partial η} + (N t) {(\frac{\partial θ (η; ℜ)}{\partial η})}^{2}, \end{array}

(28)

N_{ϕ} [θ (η; ℜ), ϕ (η; ℜ)] = \frac{\partial^{2} ϕ (η; ℜ)}{\partial η^{2}} + L e (f + c g) (\frac{\partial ϕ (η; ℜ)}{\partial η}) + \frac{N t}{N b} (\frac{\partial^{2} ϕ (η; ℜ)}{\partial η^{2}}) .

(29)

the zeroth-order problems are defined as

(1 - ℜ) L_{f} [f (η; ℜ) - f_{0} (η)] = ℜ ℏ_{f} N_{f} [f (η; ℜ), g (η; ℜ)],

(30)

(1 - ℜ) L_{g} [g (η; ℜ) - g_{0} (η)] = ℜ ℏ_{g} N_{g} [f (η; ℜ), g (η; ℜ)],

(31)

(1 - ℜ) L_{θ} [θ (η; ℜ) - θ_{0} (η)] = ℜ ℏ_{θ} N_{θ} [f (η; ℜ), θ (η; ℜ)],

(32)

(1 - ℜ) L_{ϕ} [ϕ (η; ℜ) - ϕ_{0} (η)] = ℜ ℏ_{ϕ} N_{ϕ} [θ (η; ℜ), ϕ (η; ℜ)] .

(33)

The equivalent boundary conditions are

\begin{array}{l} {f (η; ℜ) |}_{η = 0} = 0, {\frac{\partial f (η; ℜ)}{\partial η} |}_{η = 0} = 0, {\frac{\partial f (η; ℜ)}{\partial η} |}_{η \to \infty} = 0, \\ {g (η; ℜ) |}_{η = 0} = c, {\frac{\partial g (η; ℜ)}{\partial η} |}_{η = 0} = 0, {g (η; ℜ) |}_{η \to \infty} = 0, \\ {\frac{\partial θ (η; ℜ)}{\partial η} |}_{η = 0} = - B i_{1} (1 - θ), {θ (η; ℜ) |}_{η \to \infty} = 0, \\ {ϕ (η; ℜ) |}_{η = 0} = - B i_{2} (1 - ϕ), {ϕ (η; ℜ) |}_{η \to \infty} = 0. \end{array}

(34)

where

ℜ \in [0, 1]

is the implanting parameter,

ℏ_{f}, ℏ_{g}, ℏ_{θ}, ℏ_{ϕ}

are the parameters which were used to regulate the convergence of the solution. When

ℜ = 0 and ℜ = 1

, we have

\begin{array}{l} f (η; 0) = f_{0} (η), f (η; 1) = f (η), \\ g (η; 0) = g_{0} (η), g (η; 1) = g (η), \\ θ (η; 0) = θ_{0} (η), θ (η; 1) = θ (η), \\ ϕ (η; 0) = ϕ_{0} (η), ϕ (η; 1) = ϕ (η) . \end{array}

(35)

Expanding

f (η; ℜ), g (η; ℜ), θ (η; ℜ) and ϕ (η; ℜ)

by Taylor’s series

\begin{array}{l} f (η; ℜ) = f_{0} (η) + \sum_{q = 1}^{\infty} f_{q} (η) ℜ^{q}, \\ g (η; ℜ) = g_{0} (η) + \sum_{q = 1}^{\infty} g_{q} (η) ℜ^{q}, \\ θ (η; ℜ) = θ_{0} (η) + \sum_{q = 1}^{\infty} θ_{q} (η) ℜ^{q}, \\ ϕ (η; ℜ) = ϕ_{0} (η) + \sum_{q = 1}^{\infty} ϕ_{q} (η) ℜ^{q} . \end{array}

(36)

where

\begin{array}{l} f_{q} (η) = {\frac{1}{q!} \frac{\partial f (η; ℜ)}{\partial η} |}_{ℜ = 0}, g_{q} (η) = {\frac{1}{q!} \frac{\partial g (η; ℜ)}{\partial η} |}_{ℜ = 0}, \\ θ_{q} (η) = {\frac{1}{q!} \frac{\partial θ (η; ℜ)}{\partial η} |}_{ℜ = 0} and ϕ_{q} (η) = {\frac{1}{q!} \frac{\partial ϕ (η; ℜ)}{\partial η} |}_{ℜ = 0} . \end{array}

(37)

The secondary constraints

ℏ_{f}, ℏ_{g}, ℏ_{θ} and ℏ_{ϕ}

are nominated in such a way that the Series (35) converges at

ℜ = 1

. Changing

ℜ = 1

in (35), we get

\begin{array}{l} f (η) = f_{0} (η) + \sum_{q = 1}^{\infty} f_{q} (η), \\ g (η) = g_{0} (η) + \sum_{q = 1}^{\infty} g_{q} (η), \\ θ (η) = θ_{0} (η) + \sum_{q = 1}^{\infty} θ_{q} (η), \\ ϕ (η) = ϕ_{0} (η) + \sum_{q = 1}^{\infty} ϕ_{q} (η) . \end{array}

(38)

The

q^{t h}

-order problem satiates the following

\begin{array}{l} L_{f} [f_{q} (η) - χ_{q} f_{q - 1} (η)] = ℏ_{f} U_{q}^{f} (η), \\ L_{g} [g_{q} (η) - χ_{q} g_{q - 1} (η)] = ℏ_{g} U_{q}^{g} (η), \\ L_{θ} [θ_{q} (η) - χ_{q} θ_{q - 1} (η)] = ℏ_{θ} U_{q}^{θ} (η), \\ L_{ϕ} [ϕ_{q} (η) - χ_{q} ϕ_{q - 1} (η)] = ℏ_{ϕ} U_{q}^{ϕ} (η), \end{array}

(39)

The boundary conditions are

\begin{array}{l} f (0) = f^{'} (0) = f^{'} (\infty) = 0, \\ g^{'} (0) - c = g^{'} (0) = g^{'} (\infty) = 0, \\ θ^{'} (0) + B i_{1} (1 - θ (0)) = θ (\infty) = 0, \\ ϕ^{'} (0) + B i_{2} (1 - ϕ (0)) = ϕ (\infty) = 0. \end{array}

(40)

Here

U_{q}^{f} (η) = {f^{‴}}_{q - 1} + {({f^{'}}_{q - 1})}^{2} + (f + c g) {f^{″}}_{q - 1} + W e (\sum_{k = 0}^{q - 1} {f^{″}}_{q - 1} {f^{‴}}_{k}) + g {f^{″}}_{q - 1} - (γ) {f^{'}}_{q - 1},

(41)

U_{q}^{g} (η) = {g^{‴}}_{q - 1} + {({g^{'}}_{q - 1})}^{2} + (f + c g) {g^{″}}_{q - 1} + W e (\sum_{k = 0}^{q - 1} {g^{″}}_{q - 1} {g^{‴}}_{k}) + f {g^{″}}_{q - 1} - γ ({g^{'}}_{q - 1}),

(42)

U_{q}^{θ} (η) = \frac{1}{\Pr} (1 + \frac{4}{3} R d) {θ^{″}}_{q - 1} + (f + c g) {θ^{'}}_{q - 1} + N b (\sum_{k = 0}^{q - 1} {ϕ^{'}}_{q - 1 - k} {θ^{'}}_{k}) + N t {({θ^{'}}_{q - 1})}^{2},

(43)

U_{q}^{ϕ} (η) = {θ^{″}}_{q - 1} + (f + c g) {ϕ^{'}}_{q - 1} + N b (\sum_{k = 0}^{q - 1} {ϕ^{'}}_{q - 1 - k} {θ^{'}}_{k}) + \frac{N t}{N b} ({ϕ^{″}}_{q - 1}) .

(44)

where

χ_{q} = {\begin{cases} 0, if ℜ \leq 1 \\ 1, if ℜ > 1 \end{cases} .

(45)

4. Discussion

The impacts of emerging parameters have been presented in this section. The emerging parameters are stretching parameter

(c)

, Williamson parameter

(W e)

, Porosity parameter

(γ)

, thermal radiation parameter

(R d)

, Prandtl number

(\Pr)

, thermophoresis parameter

(N t)

, Brownian motion parameter

(N b)

and Biot numbers

(B i_{1}, B i_{2})

on velocities profiles

(f^{'} (η), g^{'} (η))

, temperature profile

θ (η)

and concentration profile

ϕ (η)

.Figure 1 shoe the physical sketch of the flow. Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 are plotted to observe the impact of these emerging parameters. Figure 2a is presented to see the impression of stretching parameter

c

f^{'} (η), g^{'} (η)

θ (η)

and

ϕ (η)

. Physically the stretching parameter is the ratio of the fluid velocity along the x-direction to the fluid velocity along the y-direction. So this parameter is directly related to the velocity in the y-direction and related indirectly to the velocity in the x-direction. Because of this reason, the stretching parameter reduces

f^{'} (η)

in the x-direction, as shown in Figure 2a, and the stretching parameter enhances

g^{'} (η)

in the y-direction, as shown in Figure 2b. The influences of stretching parameter

c

θ (η)

and

ϕ (η)

are shown in Figure 2c,d. The increasing values of stretching parameter

c

show a reduction in

θ (η)

and

ϕ (η)

. Figure 3a,b shows the impact of

(W e)

f^{'} (η), g^{'} (η)

in the x- and y-directions respectively. This parameter is the ratio of the relaxation to retardation time. The relaxation time is enhanced by increasing the Williamson parameter. Because of this, the fluid viscosity increases and as a result the velocity of fluid particles falls. Figure 4a,b shows the impact of

(γ)

f^{'} (η), g^{'} (η)

. Physically, the porous medium influences the boundary layer, which leads to acceleration of the fluid. Due to this fact, the velocity profile reduces with the enhancement of the porosity parameter. Figure 5 shows the impact of

(R d)

θ (η)

. We see from the figure that the enhancement of

(R d)

upsurges

θ (η)

(R d)

plays a significant role in the heat transmission process. When we increase thermal radiation, the temperature of the boundary layer flow enhances. Figure 6 depicts the impact of

(\Pr)

θ (η)

. Physically, with the increasing in the Prandtl number, the thermal diffusion diminishes, and also a small amplitude of oscillation in temperature occurs. Therefore,

θ (η)

reduces with increasing

(\Pr)

. Figure 7a, b depicts the impact of

(N t)

θ (η)

and

ϕ (η)

. Here we see that the enhancement of

(N t)

reduces

θ (η)

and

ϕ (η)

. In addition, thermophoresis is high in Newtonian fluids as compared to non-Newtonian fluids. Figure 8a,b depicts the impact of

(N b)

θ (η)

and

ϕ (η)

. In Figure 8a, we see that

θ (η)

enhances with the enhancement of

(N b)

, while the opposite behavior of

(N b)

ϕ (η)

can be seen in Figure 8b. Figure 9a–c depicts the impact of

(B i_{1}, B i_{2})

θ (η)

and

ϕ (η)

. From the figure, we can see that the increase of Biot numbers

(B i_{1}, B i_{2})

causes an increase in the temperature and concentration profiles. The numerical values of the surface drag force and impact of the different parameters;

W e

γ

and

c

, are shown in Table 1. It is observed that increasing

c

and

γ

reduces the surface drag force

f^{″} (0)

and

g^{″} (0)

, while increasing

W e

augments the surface drag force

f^{″} (0)

and

g^{″} (0)

. Table 2 displays the numerical values of the heat transfer rate for the different embedded parameters;

c, \Pr, R d, N b

and

N t

. It is concluded from Table 2, that the higher values of

c, \Pr

and

R d

enhance the heat transfer rate, while the augmented value of

N b

and

N t

decrease the mass transfer rate. Table 2 presents the numerical values of the mass transfer rate for the different embedded parameters;

c, L e, N b

and

N t

. Here we observed that the higher values of

c, N b

and

N t

enhance the mass transfer rate, while the augmented value of

L e

reduces the heat transfer rate. Table 1, Table 2 and Table 3 display the validation of the current results with the previous published results. Both the current and previous results are compared and an excellent outcome is found. This verifies the validity of the current results along with the accurateness of the analytical method (HAM) we used in this research work.

5. Conclusions

In this article we have analyzed three-dimensional Williamson fluid flow past a linear porous and stretching sheet with the effect of thermal radiation. The system of equations has been solved by the homotopy analysis method. The impact of implanted parameters has been deliberated through graphs and discussed in detail. The concluding observations are as follows:

The velocity function

f^{'} (η)

reduces with enhancement of

f^{'} (η)

(W e)

and

(γ)

The velocity function $g^{'} (η)$ improves with enhancement of the stretching parameter, while reduces with an increase in $(γ)$ and $(W e)$ .
The temperature function $θ (η)$ improves with an increase in $(R d)$ , $(N t)$ , $(N b)$ and $(B i_{1}, B i_{2})$ , while reduces with an increase in $(c)$ and $(\Pr)$ .
The concentration function $ϕ (η)$ enhances with an increase in $(N t)$ and $(B i_{2})$ , while reduces with an increase in $(c)$ and $(N b)$ .
It is observed that increasing $c$ and $γ$ reduces the surface drag forces $f^{″} (0)$ and $g^{″} (0)$ , while increasing $W e$ augmented the surface drag forces $f^{″} (0)$ and $g^{″} (0)$ .
$N t$ . It is concluded from Table 2, that the higher values of $c, \Pr$ and $R d$ enhance the heat transfer rate, while the augmented values of $N b$ and $N t$ decline the heat transfer rate.

Author Contributions

A.S.K., Y.N. and Z.S. modeled and solved the problem and wrote the manuscript. S.I. and A.A.D. thoroughly checked the mathematical modeling and English corrections. W.K. contributed in the results and discussions and solution. All the corresponding authors finalized the manuscript after its internal evaluation.

Funding

This research was funded by National Natural Science Foundation of China No. [11471262].

Conflicts of Interest

The authors declare no conflict of interest.

References

Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. Am. Inst. Chem. Eng. 1961, 7, 26–28. [Google Scholar] [CrossRef]
Crane, L.J. Flow past a stretching plate. Int. J. Appl. Mech. 1970, 21, 645–647. [Google Scholar] [CrossRef]
Williamson, R.V. The flow of pseudoplastic materials. Int. J. Ind. Eng. Chem. 1929, 21, 1108–1111. [Google Scholar] [CrossRef]
Ariel, P.D. The three-dimensional flow past a stretching sheet and the homotopy perturbation method. Int. J. Appl. Math. 2007, 54, 920–925. [Google Scholar] [CrossRef]
Nadeem, S.; Haq, R.U.; Lee, C. MHD flow of a Casson fluid over an exponentially shrinking sheet. Sci. Iran. 2012, 19, 1550–1553. [Google Scholar] [CrossRef]
Nadeem, S.; Akram, S. Peristaltic flow of a Williamson fluid in an asymmetric channel. Int. J. Nonlinear Sci. Numer. Simul. 2010, 15, 1705–1712. [Google Scholar] [CrossRef]
Pop, I.; Na, T.Y. Unsteady flow past a stretching sheet. Int. J. Acta Mech. 1996, 4, 413–422. [Google Scholar] [CrossRef]
Sakiadis, B.C. Boundary layer behavior on continuous solid surfaces and boundary layer on a continuous flat surface. J. Am. Inst. Chem. Eng. 1961, 7, 221–225. [Google Scholar] [CrossRef]
Tsou, F.K.; Sparrow, E.M.; Goldstein, R.J. Flow and heat transfer in the boundary layer on a continuous moving surface. Int. J. Heat Mass Transf. 1967, 10, 219–235. [Google Scholar] [CrossRef]
Malik, M.Y.; Bilal, S.; Salahuddin, T.; Khalil, U.R. Three-Dimensional Williamson Fluid Flow over a Linear Stretching Surface. Math. Sci. Lett. 2017, 6, 53–61. [Google Scholar] [CrossRef]
Gupta, P.S.; Gupta, A.S. Heat and mass transfer on a stretching sheet with suction or blowing. Can. J. Chem. Eng. 1977, 55, 744–746. [Google Scholar] [CrossRef]
Ericksod, L.E.; Fan, L.T.; Fox, V.G.; McCormack, P.D.; Crane, L. Phvsical fluid dynamics. Ind. Eng. Chem. 1966, 6, 19. [Google Scholar]
Abdei-Aziz, M. Effect on the flow and heat transfer over an unsteady stretching sheet. Int. Commun. Heat Mass Transf. 2009, 36, 521–524. [Google Scholar] [CrossRef]
Mukhopadyay, S. Effect of thermal radiation on unsteady mixed convection flow and heat treansfer over a porous stretching surface in porous medium. Int. J. Heat Mass Transf. 2009, 52, 3261–3265. [Google Scholar] [CrossRef]
Shateyi, S.; Motsa, S.S. Thermal radiation effects on heat and mass transfer over an unsteady stretching surface. Math. Probl. Eng. 2009, 13, 965603. [Google Scholar] [CrossRef]
Abdei-Aziz, M. Thermal-diffusion and diffusion-thermo effects on combined heat and mass transfer by hydromagnetic three-dimensional free convection over a permeable stretching surface with radiation. Phys. Lett. 2007, 372, 263–272. [Google Scholar] [CrossRef]
Cortell, R. Effects of viscous dissipation and radiation on thermal boundary layer over a non-linearly stretching sheet. Phys. Lett. A 2008, 372, 631–636. [Google Scholar] [CrossRef]
Ishak, A.; Nazar, R.; Pop, I. Heat transfer over a stretching surface with variable heat flux in viscous fluid. Phys. Lett. A 2008, 372, 559–561. [Google Scholar] [CrossRef]
Dawar, A.; Shah, Z.; Islam, S.; Idress, M.; Khan, W.; Islam, S.; Gul, T. Impact of Thermal Radiation and Heat Source/Sink on Eyring–Powell Fluid Flow over an Unsteady Oscillatory Porous Stretching Surface. Math. Comput. Appl. 2018, 23, 20. [Google Scholar] [CrossRef]
Dawar, A.; Shah, Z.; Islam, S.; Idress, M.; Khan, W. Magnetohydrodynamic CNTs Casson Nanofl uid and Radiative heat transfer in a Rotating Channels. J. Phys. Res. Appl. 2018, 1, 017–032. [Google Scholar] [CrossRef]
Kashif, A.A.; Irfan, A.A.; Sikandar, M.A.; Ilyas, K. On the thermal analysis of magnetohydrodynamic Jeffery fluid via modern non integer order derivative. J. King Saud Univ. Sci. 2018. [Google Scholar] [CrossRef]
Sheikholeslami, M. Numerical investigation of nanofluid free convection under the influence of electric field in a porous enclosure. J. Mol. Liq. 2018, 249, 1212–1221. [Google Scholar] [CrossRef]
Sheikholeslami, M. Fe₃O₄-Ethylene glycol nanofluid forced convection inside a porous enclosure in existence of Coulomb force. Mol. Liq. 2018, 249, 429–437. [Google Scholar] [CrossRef]
Sheikholeslami, M. Magnetic field influence on nanofluid thermal radiation in a cavity with tilted elliptic inner cylinder. J. Mol. Liq. 2017, 229, 137–147. [Google Scholar] [CrossRef]
Sheikholeslami, M. Magnetohydrodynamic nanofluid forced convection in a porous lid driven cubic cavity using Lattice Boltzmann method. J. Mol. Liq. 2017, 231, 555–565. [Google Scholar] [CrossRef]
Nadeem, S.; Ul Haq, R.; Noreen, S.A.; Khan, Z.H. MHD three-dimensional Casson fluid flow past a porous linearly stretching sheet. Alex. Eng. J. 2013, 52, 577–582. [Google Scholar] [CrossRef]
Shah, Z.; Islam, S.; Ayaz, H.; Khan, S. Radiative Heat and Mass Transfer Analysis of Micropolar Nanofluid Flow of Casson Fluid between Two Rotating Parallel Plates with Effects of Hall Current. ASME J. Heat Transf. 2018. [Google Scholar] [CrossRef]
Shah, Z.; Bonyah, E.; Islam, S.; Khan, W.; Ishaq, M. Radiative MHD thin film flow of Williamson fluid over an unsteady permeable stretching. Heliyon 2018, 4, e00825. [Google Scholar] [CrossRef] [PubMed]
Shah, Z.; Dawar, A.; Islam, S.; Khan, I.; Ching, D.L.C. Darcy-Forchheimer Flow of Radiative Carbon Nanotubes with Microstructure and Inertial Characteristics in the Rotating Frame. Case Stud. Therm. Eng. 2018. [Google Scholar] [CrossRef]
Khan, A.; Shah, Z.; Islam, S.; Khan, S.; Khan, W.; Khan, Z.A. Darcy–Forchheimer flow of micropolar nanofluid between two plates in the rotating frame with non-uniform heat generation/absorption. Adv. Mech. Eng. 2018, 10, 1–16. [Google Scholar] [CrossRef]
Jawad, M.; Shah, Z.; Islam, S.; Islam, S.; Bonyah, E.; Khan, Z.A. Darcy-Forchheimer flow of MHD nanofluid thin film flow with Joule dissipation and Navier’s partial slip. J. Phys. Commun. 2018. [Google Scholar] [CrossRef]
Khan, N.; Zuhra, S.; Shah, Z.; Bonyah, E.; Khan, W.; Islam, S. Slip flow of Eyring-Powell nanoliquid film containing graphene nanoparticles. AIP Adv. 2018, 8, 115302. [Google Scholar] [CrossRef]
Bhatti, M.A.; Tehseen, R.; Mohammad, A.; Mohamed, Y.Z. Entropy generation on MHD Eyring-Powell nanofluid through a permeable stretching surface. Entropy 2016, 18, 224. [Google Scholar] [CrossRef]
Darbari, B.; Rashidi, S.; Abolfazli Esfahani, J. Sensitivity analysis of entropy generation in nanofluid flow inside a channel by response surface methodology. Entropy 2016, 18, 52. [Google Scholar] [CrossRef]
Hassan, M.; Marin, M.; Abdullah, A.; Ellahi, R. Convection heat transfer flow of nanofluid in a porous medium over wavy surface. Phys. Lett. A 2018, 382, 2749–2753. [Google Scholar] [CrossRef]
Ahmed, Z.; Nouman, I.; Tehseen, A.; Rahmat, E. The sustainable characteristic of Bio-bi-phase flow of peristaltic transport of MHD Jeffery fluid in human body. Sustainability 2018, 10, 2671. [Google Scholar]
Ghofrane, S.; Sébastien, P. Further investigation on laminar forced convection of nanofluid flows in a uniformly heated pipe using direct numerical simulations. Appl. Sci. 2016, 6, 332. [Google Scholar]
Hussain, F.; Ellahi, R.; Zeeshan, A. Mathematical models of electro magnetohydrodynamic Multiphase Flows Synthesis with Nanosized Hafnium Particles. Appl. Sci. 2018, 8, 275. [Google Scholar] [CrossRef]
Pieters, G.J.M.; Duijn, C.J.V. Transient growth in linearly stable gravity-driven flow in porous media. Eur. J. Mech. 2006, 25, 83–94. [Google Scholar] [CrossRef] [Green Version]
Machrafi, H.; Rednikov, A.; Colinet, P.; Dauby, P.C. Bénard instabilities in a binary-liquid layer evaporating into an inert gas: Stability of quasi-stationary and time-dependent reference profiles. Eur. Phys. J. Spec. Top. 2011, 192, 71–81. [Google Scholar] [CrossRef] [Green Version]
Ozen, O.; Narayanan, R. The physics of evaporative and convective instabilities in bilayer systems. Linear Theory Phys. Fluids 2004, 16, 4644. [Google Scholar] [CrossRef]
Lapwood, E.R. Convection of a fluid in a porous medium. Math. Proc. Camb. Philos. Soc. 1948, 44, 508–521. [Google Scholar] [CrossRef]
Sulochana, C.; Ashwinkumar, G.P.; Sandeep, N. Similarity solution of 3D Casson nanofluid flow over a stretching sheet with convective boundary conditions. J. Niger. Math. Soc. 2016, 35, 128–141. [Google Scholar] [CrossRef]
Taza, G.; Abdul, S.K.; Saeed, I.; Aisha, M.A.; Ilyas, K.; Ali, S.A.; Abdullah, K.A. Muradullah. Heat transfer investigation of the unsteady thin film flow of Williamson fluid past an inclined and oscillating moving plate. Appl. Sci. 2017, 7, 369. [Google Scholar]
Jang, S.P.; Choi, S.U.S. Role of Brownian motion in the enhanced thermal conductivity of nanofluids. Appl. Phys. Lett. 2004, 84, 4316. [Google Scholar] [CrossRef]
Machrafi, H.; Lebon, G. The role of several heat transfer mechanisms on the enhancement of thermal conductivity in nanofluids. Contin. Mech. Thermodyn. 2016, 28, 1461–1475. [Google Scholar] [CrossRef] [Green Version]
Anoop, K.B.; Kabelac, S.; Sundararajan, T.; Das, S.K. Rheological and flow characteristics of nanofluids: Influence of electroviscous effects and particle agglomeration. J. Appl. Phys. 2009, 106, 034909. [Google Scholar] [CrossRef]
Azizian, R.; Doroodchi, D.; Moghtaderi, B. Effect of Nanoconvection Caused by Brownian Motion on the Enhancement of Thermal Conductivity in Nanofluids. Ind. Eng. Chem. Res. 2012, 51, 1782–1789. [Google Scholar] [CrossRef]
Hammed, K.; Haneef, M.; Shah, Z.; Islam, I.; Khan, W.; Asif, S.M. The Combined Magneto hydrodynamic and electric field effect on an unsteady Maxwell nanofluid Flow over a Stretching Surface under the Influence of Variable Heat and Thermal Radiation. Appl. Sci. 2018, 8, 160. [Google Scholar] [CrossRef]
Muhammad, S.; Ali, G.; Shah, Z.; Islam, S.; Hussain, A. The Rotating Flow of Magneto Hydrodynamic Carbon Nanotubes over a Stretching Sheet with the Impact of Non-Linear Thermal Radiation and Heat Generation/Absorption. Appl. Sci. 2018, 8, 482. [Google Scholar] [CrossRef]

Figure 1. Physical representation of the problem.

Figure 2. (a–d) Effect of the stretching parameter on velocity, temperature, and concentration profiles. (a) Effect of

c

f^{'} (η)

, when

W e = 0.1, γ = 0.2.

(b) Effect of

c

g^{'} (η)

, when

W e = 0.1, γ = 0.2.

c

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, R d = 0.8.

(d) Effect of

c

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N b = 0.4, N t = 0.5

Figure 2. (a–d) Effect of the stretching parameter on velocity, temperature, and concentration profiles. (a) Effect of

c

f^{'} (η)

, when

W e = 0.1, γ = 0.2.

(b) Effect of

c

g^{'} (η)

, when

W e = 0.1, γ = 0.2.

c

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, R d = 0.8.

(d) Effect of

c

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N b = 0.4, N t = 0.5

Figure 3. (a,b) Influence of the Williamson parameter on velocity profiles. (a) Effect of

W e

f^{'} (η)

, when

c = 0.1, γ = 0.2.

(b) Effect of

W e

g^{'} (η)

, when

c = 0.1, γ = 0.2.

Figure 3. (a,b) Influence of the Williamson parameter on velocity profiles. (a) Effect of

W e

f^{'} (η)

, when

c = 0.1, γ = 0.2.

(b) Effect of

W e

g^{'} (η)

, when

c = 0.1, γ = 0.2.

Figure 4. (a,b) Effect of

(γ)

f^{'} (η), g^{'} (η)

. (a) Influence of

γ

f^{'} (η)

, when

c = 0.1, W e = 0.2.

(b) Influence of

γ

g^{'} (η)

, when

c = 0.1, W e = 0.2

Figure 4. (a,b) Effect of

(γ)

f^{'} (η), g^{'} (η)

. (a) Influence of

γ

f^{'} (η)

, when

c = 0.1, W e = 0.2.

(b) Influence of

γ

g^{'} (η)

, when

c = 0.1, W e = 0.2

Figure 5. Effect of the thermal radiation parameter

R d

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8

Figure 5. Effect of the thermal radiation parameter

R d

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8

Figure 6. Effect of Prandtl number

\Pr

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, R d = 0.6, L e = 0.7, c = 0.8

Figure 6. Effect of Prandtl number

\Pr

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, N t = 0.5, R d = 0.6, L e = 0.7, c = 0.8

Figure 7. (a,b) Effect of the thermophoresis parameter on temperature and concentration profiles. (a) Effect of

N t

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, R d = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

N t

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N b = 0.4, L e = 0.7, c = 0.8

Figure 7. (a,b) Effect of the thermophoresis parameter on temperature and concentration profiles. (a) Effect of

N t

θ (η)

, when

B i_{1} = 0.1, N b = 0.4, R d = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

N t

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N b = 0.4, L e = 0.7, c = 0.8

Figure 8. (a,b) Effect of the Brownian motion parameter on temperature and concentration profiles. (a) Effect of

N b

θ (η)

, when

B i_{1} = 0.1, R d = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

N b

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N t = 0.4, L e = 0.7, c = 0.8

Figure 8. (a,b) Effect of the Brownian motion parameter on temperature and concentration profiles. (a) Effect of

N b

θ (η)

, when

B i_{1} = 0.1, R d = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

N b

ϕ (η)

, when

B i_{1} = 0.1, B i_{2} = 0.2, N t = 0.4, L e = 0.7, c = 0.8

Figure 9. (a–c) Effect of the Biot numbers on temperature and concentration profiles. (a) Effect of

B i_{1}

θ (η)

, when

R d = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

B i_{1}

ϕ (η)

, when

N b = 0.1, B i_{2} = 0.2, N t = 0.4, L e = 0.7, c = 0.8.

B i_{2}

ϕ (η)

, when

N b = 0.1, B i_{1} = 0.2, N t = 0.4, L e = 0.7, c = 0.8

Figure 9. (a–c) Effect of the Biot numbers on temperature and concentration profiles. (a) Effect of

B i_{1}

θ (η)

, when

R d = 0.1, N b = 0.4, N t = 0.5, \Pr = 0.6, L e = 0.7, c = 0.8.

(b) Effect of

B i_{1}

ϕ (η)

, when

N b = 0.1, B i_{2} = 0.2, N t = 0.4, L e = 0.7, c = 0.8.

B i_{2}

ϕ (η)

, when

N b = 0.1, B i_{1} = 0.2, N t = 0.4, L e = 0.7, c = 0.8

Table 1. Numerical values of surface drag force for different embedded parameters.

$c$	$W e$	$γ$	$C_{f}$ $for f^{″} (0)$ Previous Results [43]	$C_{f}$ $for f^{″} (0) Present Results$	$C_{f}$ $for g^{″} (0)$ Previous Results [43]	$C_{f}$ $for g^{″} (0) Present Results$
1.0	0.1	0.5	−1.489532	−1.23942	−0.695215	−2.04702
1.1			−1.564022	−1.25192	−1.564022	−2.23996
1.2			−1.682498	−1.26417	−2.647314	−2.38904
	0.2		−1.731420	−1.06979	−1.892121	−1.22441
	0.3		−1.544321	−0.83279	−1.744322	−0.42052
	0.4		−1.324151	−0.38181	−1.212141	0.35420
		0.6	−1.829575	−1.28570	−1.829575	−1.82788
		0.7	−2.083805	−1.32900	−2.083805	−1.60457
		0.8	−2.310564	−1.36945	−2.310564	−1.38981

Table 2. Numerical values of heat transfer rate for different embedded parameters.

$c$	$\Pr$	$R d$	$N t$	$N b$	$N u$ Previous Results [43]	$N u Present Results$
1.0	0.4	0.3	0.01	0.01	0.261896	0.119212
1.1					0.318272	0.120401
1.2					0.379145	0.121608
	0.5				-	0.127300
	0.6				-	0.134579
	0.7					0.140705
		0.4			-	0.155214
		0.5			-	0.192457
		0.6			-	0.241744
			0.02		0.293770	0.155120
			0.06		0.274299	0.154747
			0.08		0.255629	0.154467
				0.03	0.285647	0.155044
				0.05	0.253739	0.154874
				0.07	0.222739	0.154704

Table 3. Numerical values of mass transfer rate for different embedded parameters.

$c$	Le	$N t$	$N b$	$S h$ Previous Results [43]	$S h Present Results$
1.0	0.1	0.01	0.01	0.223924	0.291024
1.1				0.255507	0.293757
1.2				0.289266	0.296524
	0.2			-	0.291022
	0.3			-	0.291021
	0.4			-	0.291019
		0.02		0.213208	0.200865
		0.03		0.144515	0.290705
		0.04		0.089302	0.290545
			0.02	0.272275	0.290847
			0.03	0.286051	0.290671
			0.04	0.290593	0.290494

© 2018 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

Khan, A.S.; Nie, Y.; Shah, Z.; Dawar, A.; Khan, W.; Islam, S. Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions. Appl. Sci. 2018, 8, 2244. https://doi.org/10.3390/app8112244

AMA Style

Khan AS, Nie Y, Shah Z, Dawar A, Khan W, Islam S. Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions. Applied Sciences. 2018; 8(11):2244. https://doi.org/10.3390/app8112244

Chicago/Turabian Style

Khan, Abdul Samad, Yufeng Nie, Zahir Shah, Abdullah Dawar, Waris Khan, and Saeed Islam. 2018. "Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions" Applied Sciences 8, no. 11: 2244. https://doi.org/10.3390/app8112244

APA Style

Khan, A. S., Nie, Y., Shah, Z., Dawar, A., Khan, W., & Islam, S. (2018). Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions. Applied Sciences, 8(11), 2244. https://doi.org/10.3390/app8112244

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

Article Menu

Three-Dimensional Nanofluid Flow with Heat and Mass Transfer Analysis over a Linear Stretching Surface with Convective Boundary Conditions

Abstract

1. Introduction

2. Problem Formulation

3. Solution by Homotopy Analysis Method

4. Discussion

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI