# Mixed Convection over an Inclined Wavy Surface in a Nanofluid Saturated Non-Darcy Porous Medium with Radiation Effect

Authors: Darbhasayanam Srinivasacharya; Poshala Vijay Kumar
Journal: International Journal of Chemical Engineering
 (2015)

---

## Abstract

The problem of mixed convection in a nanofluid along an inclined wavy surface embedded in a non-Darcy porous medium with radiation effect is analyzed. Coordinate transformation is employed to transform the complex wavy surface to a smooth surface. The governing equations are transformed into a set of ordinary differential equations using the appropriate similarity transformation and then solved using the successive linearization method. The present results are compared with previously published work and are found to be in very good agreement. The effect of pertinent parameters on the nondimensional velocity, temperature, nanoparticle volume fraction, heat, and nanoparticle mass transfer rates is studied and presented graphically.


---

## Body

## 1. Introduction

The research on heat transfer in nanofluids is receiving much attention from the past several years due to its unique characteristics in thermal engineering fields and broad range of applications such as microelectronics, microfluidics, solid-state lighting, transportation, biomedical, and manufacturing. To enhance the thermal conductivity and viscosity of the base fluid, colloidal suspensions of nanometer sized particles (of dimensions about 1–100 nm) are made to form nanofluid [1]. Choi et al. [2] showed that the addition of a small amount (less than 1% by volume) of nanoparticles to conventional heat transfer liquids increased the thermal conductivity of the fluid. The nanofluids are advantageous in the sense that they are stable and have an acceptable viscosity and better wetting, spreading, and dispersion properties. A detailed review of the literature on nanofluids and applications of nanofluids can be found in the book by Das et al. [3] and collection of papers by Kakaç and Pramuanjaroenkij [4] and Puliti et al. [5]. These reviews discuss in detail the work done on convective transport in nanofluids.Numerous models and methods have been proposed by different authors to study convective flows of nanofluids. The computational studies reported in this area include two main approaches:(
1
) a two-phase model, in which both liquid and solid heat transfer behaviors are solved in the flow fields [6, 7], and (
2
) a single-phase model, in which solid particles are considered to behave as fluids, because the nanoparticles are easily fluidized [8–10]. Buongiorno [11] showed that the high heat transfer coefficients in nanofluids cannot be explained satisfactorily by thermal dispersion phenomenon or increase in turbulence intensity promoted by the presence of nanoparticles or nanoparticle rotation as suggested in the literature. He noted that the nanoparticle absolute velocity can be viewed as the sum of the base fluid velocity and a relative velocity (which he calls the slip velocity). He developed a mathematical model to capture the nanoparticle/base fluid slip by treating nanofluid as a two-component mixture with Brownian diffusion and thermophoresis as the important factors in the convective transport process in a nanofluid. Kuznetsov and Nield [12] used the Buongiorno [11] model to study the natural convective flow of a nanofluid over a vertical plate. Buongiorno [11] model is the basis of the present study.Convective transport in porous media has been widely studied in recent years due to its wide range of applications in mechanical, chemical, and civil engineering. A review of convective heat transfer in porous medium is presented in the book by Nield and Bejan [13]. Considerable work has been reported in the literature focusing on the problem of mixed convection of nanofluid along an inclined surface embedded in porous medium. Ahmad and Pop [14] studied the steady mixed convection boundary layer flow past a vertical flat plate embedded in a porous medium filled with nanofluids using different types of nanoparticles and the model they used for the nanofluid incorporates only the nanoparticle volume fraction parameter. Chamkha and Aly [15] obtained numerical solution of steady natural convection boundary layer flow of a nanofluid consisting of a pure fluid with nanoparticles along a permeable vertical plate in the presence of magnetic field, heat generation or absorption, and suction or injection effects. Rosca et al. [16] numerically studied the problem of steady mixed convection boundary layer flow over an impermeable horizontal flat plate embedded in a porous medium saturated by a nanofluid. Rana et al. [17] obtained the numerical solution for steady mixed convection boundary layer flow of an incompressible nanofluid along an inclined plate in a porous medium. Kamali and Binesh [18] numerically studied the effects of nanoparticles size on convective heat transfer characteristics of metal oxide nanofluids considering thermal conductivity and viscosity of nanofluid as a function of nanoparticles size, volume fraction, and Brownian motion.Radiation effects on convection gained importance in the context of many industrial applications involving high temperatures such as gas turbines, nuclear power plant, and various propulsion engines for aircraft, missiles, satellites, and space technology. Chen et al. [19] studied the effects of thermal radiation on laminar forced and free convection along a wavy surface. Chamkha et al. [20] presented the nonsimilar solution of steady mixed convection of a nanofluid in the presence of thermal radiation. Hady et al. [21] studied the effect of radiation on viscous flow of a nanofluid and heat transfer over a nonlinearly stretching sheet. Abdel-Rahman [22] numerically investigated the problem of thermal radiation and unsteady MHD flow of a nanofluid in stretching porous medium.The prediction of heat and mass transfer from irregular surface has great importance in engineering applications. In many practical situations surfaces are roughened in order to enhance the rate of heat transfer. The presence of roughened surface not only alters the flow field but also alters the heat and mass transfer characteristics. Natural convection from wavy surfaces is used for transferring heat in several heat transfer devices, such as flat plate condensers in refrigerators, flat plate solar collectors, and cooling of electrical and nuclear components. Extensive studies of natural convection heat and mass transfer of a vertical wavy surface under boundary layer approximation have been undertaken by several authors. Mahdy and Ahmed [23] have analyzed laminar free convection over a vertical wavy surface embedded in a porous medium saturated with a nanofluid. Ahmed and Abd El-Aziz [24] studied the effect of local thermal nonequilibrium on unsteady heat transfer by natural convection of a nanofluid over a vertical wavy surface. On the other hand, convection along inclined surfaces and bluff bodies has been receiving the attention because of many industrial applications in areas such as electroplating, chemical processing of heavy metals, and ash or scrubber waste treatment. Ingham et al. [25] studied the natural convection from a semi-infinite flat plate inclined at a small angle to the horizontal in a saturated porous medium. Cheng [26, 27] studied the double diffusive natural convection along an inclined wavy surface in a porous medium and bidisperse porous medium.In this study, we aimed at exploring the effects of radiation and the angle of inclination of a wavy surface on mixed convection in a non-Darcy porous medium saturated with a nanofluid.


## 2. Mathematical Formulation

Consider the steady laminar incompressible two-dimensional boundary layer mixed convection flow along a semi-infinite inclined wavy surface embedded in a nanofluid saturated non-Darcy porous medium. The coordinate system is shown in Figure1. The wavy surface is described by (1)
y
=
δ
x
=
a
sin
⁡
π
x
l
,where a is the amplitude of the wavy surface and 2
l is the characteristic length of the wavy surface. The wavy surface is maintained at uniform wall temperature T
w and nanoparticle concentration ϕ
w. These values are assumed to be greater than the ambient temperature T
∞ and the nanoparticle concentration ϕ
∞ at any arbitrary reference point in the medium.Figure 1
Physical model.The porous medium is considered to be homogeneous and isotropic and is saturated with a fluid which is in local thermodynamic equilibrium with the solid matrix. We assume the fluid properties to be constant except the density in the buoyancy term of the balance of momentum equation. The fluid is considered to be a gray, absorbing emitting radiation but nonscattering medium and the Rosseland approximation [28] is used to describe the radiative heat flux in the energy equation.The governing equations for this problem under the laminar boundary layer flow assumptions, Boussinesq approximation, and using the non-Darcy flow through a homogeneous porous medium near the inclined wavy surface are given by [12, 29, 30] (2)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(3)
1
+
K
~
ν
u
2
+
v
2
∂
u
∂
y
-
∂
v
∂
x
+
K
~
ν
u
2
+
v
2
u
2
∂
u
∂
y
+
u
v
∂
v
∂
y
-
∂
u
∂
x
-
v
2
∂
v
∂
x
=
1
-
ϕ
∞
ρ
f
∞
β
K
g
μ
∂
T
∂
y
sin
⁡
A
-
∂
T
∂
x
cos
⁡
A
-
ρ
p
-
ρ
f
∞
K
g
μ
∂
ϕ
∂
y
sin
⁡
A
-
∂
ϕ
∂
x
cos
⁡
A
,
(4)
u
∂
T
∂
x
+
v
∂
T
∂
y
=
α
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
+
γ
D
B
∂
ϕ
∂
x
∂
T
∂
x
+
∂
ϕ
∂
y
∂
T
∂
y
+
D
T
T
∞
∂
T
∂
x
2
+
∂
T
∂
y
2
+
16
σ
T
∞
3
3
K
e
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
,
(5)
u
∂
ϕ
∂
x
+
v
∂
ϕ
∂
y
=
D
B
∂
2
ϕ
∂
x
2
+
∂
2
ϕ
∂
y
2
+
D
T
T
∞
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
,where u and v are the velocity components in the x and y directions, respectively, T is the temperature, ϕ is the nanoparticle concentration, g is the acceleration due to gravity, K is the permeability, ρ
f is the the density of the base fluid, ρ
p is the density of the particles, ρ
C
p is the heat capacitance of the nanoparticles, ν is the kinematic viscosity of the fluid, α
=
k
/
(
ρ
c
)
f is the effective thermal diffusivity, β is the volumetric thermal expansion coefficient of the nanofluid, μ is the dynamic viscosity of the fluid, D
B is the Brownian diffusion coefficient, D
T is the thermophoretic diffusion coefficient, γ is the ratio between the effective heat capacity of the nanoparticle material and heat capacity of the fluid i
.
e
.
,
γ
=
ɛ
(
ρ
c
)
p
/
(
ρ
c
)
f, K
e and σ are the mean absorption coefficient and Stefan-Boltzmann constant, and K
~ is a material parameter which is a measure of inertia impedance of the matrix to account for non-Darcian inertial effects.The boundary conditions are(6a)
v
=
0
,
T
=
T
w
,
ϕ
=
ϕ
w
,
a
t
y
=
δ
x
,
(6b)
u
=
U
∞
,
T
⟶
T
∞
,
ϕ
⟶
ϕ
∞
a
s
y
⟶
∞
.Introducing the stream functionψ by(7)
u
=
∂
ψ
∂
y
,
v
=
-
∂
ψ
∂
xand the following nondimensional variables(8)
x
^
,
y
^
,
a
^
,
δ
^
=
l
-
1
x
,
y
,
a
,
δ
,
θ
=
T
-
T
∞
T
w
-
T
∞
,
ψ
^
=
ψ
U
∞
l
,
s
=
ϕ
-
ϕ
∞
ϕ
w
-
ϕ
∞in (3)–(5), we get the following system of nondimensional equations:(9)
1
+
F
c
Q
∂
2
ψ
^
∂
x
^
2
+
∂
2
ψ
^
∂
y
^
2
+
F
c
Q
∂
ψ
^
∂
y
^
2
∂
2
ψ
^
∂
y
^
2
+
2
∂
ψ
^
∂
x
^
∂
ψ
^
∂
y
^
∂
2
ψ
^
∂
x
^
∂
y
^
+
∂
ψ
^
∂
x
^
2
∂
2
ψ
^
∂
x
^
2
=
Δ
∂
θ
∂
y
^
sin
⁡
A
-
∂
θ
∂
x
^
cos
⁡
A
-
N
r
∂
s
∂
y
^
sin
⁡
A
-
∂
s
∂
x
^
cos
⁡
A
,
∂
ψ
^
∂
y
^
∂
θ
∂
x
^
-
∂
ψ
^
∂
x
^
∂
θ
∂
y
^
=
1
P
e
1
+
4
R
3
∂
2
θ
∂
x
^
2
+
∂
2
θ
∂
y
^
2
+
N
b
∂
s
∂
x
^
∂
θ
∂
x
^
+
∂
s
∂
y
^
∂
θ
∂
y
^
+
N
t
∂
θ
∂
x
^
2
+
∂
θ
∂
y
^
2
,
∂
ψ
^
∂
y
^
∂
s
∂
x
^
-
∂
ψ
^
∂
x
^
∂
s
∂
y
^
=
1
L
e
P
e
∂
2
s
∂
x
^
2
+
∂
2
s
∂
y
^
2
+
N
t
N
b
∂
2
θ
∂
x
^
2
+
∂
2
θ
∂
y
^
2
,where Q
=
∂
ψ
^
/
∂
x
^
2
+
∂
ψ
^
/
∂
y
^
2
1
/
2 is the nondimensional velocity, R
a
=
(
1
-
ϕ
∞
)
ρ
f
∞
β
K
g
(
T
w
-
T
∞
)
l
/
μ
α is the Rayleigh number, P
e
=
U
∞
l
/
α is the Peclet number, Δ
=
R
a
/
P
e is the mixed convection number, F
c
=
K
~
U
∞
/
ν is the non-Darcy parameter, N
r
=
(
ρ
p
-
ρ
f
∞
)
(
ϕ
w
-
ϕ
∞
)
/
ρ
f
∞
β
(
T
w
-
T
∞
)
(
1
-
ϕ
∞
) is the buoyancy ratio, R
=
4
σ
T
∞
3
/
K
K
e is the Radiation parameter, N
b
=
γ
D
B
(
ϕ
w
-
ϕ
∞
)
/
α is the Brownian motion parameter, N
t
=
γ
D
T
(
T
w
-
T
∞
)
/
α
T
∞ is the thermophoresis parameter, and L
e
=
α
/
D
B is the Lewis number.The effect of the wavy surface can be transferred from the boundary conditions into the governing equations by the coordinate transformation(10)
ξ
=
x
^
,
η
=
y
^
-
δ
P
e
1
/
2
ξ
1
/
2
1
+
δ
˙
2
,
ψ
^
=
P
e
1
/
2
ξ
1
/
2
f
ξ
,
η
.Substituting (10) into (9) and letting R
a
→
∞ (i.e., boundary layer approximation), we obtain the following equations:(11)
f
′
′
+
2
F
c
1
+
δ
˙
2
-
1
/
2
f
′
f
′
′
=
Δ
sin
⁡
A
+
δ
˙
cos
⁡
A
θ
′
-
N
r
s
′
,
1
+
4
R
3
θ
′
′
+
1
2
f
θ
′
+
N
b
s
′
θ
′
+
N
t
θ
′
2
=
0
,
s
′
′
+
1
2
L
e
f
s
′
+
N
t
N
b
θ
′
′
=
0
.The associated boundary conditions are(12a)
f
=
0
,
θ
=
1
,
s
=
1
,
a
t
η
=
0
,
(12b)
f
′
⟶
1
,
θ
⟶
0
,
s
⟶
0
,
a
s
η
⟶
∞
.The primary objective of this study is to estimate the parameters of engineering interest in fluid flow, heat, and nanoparticle mass transport problems, namely, the Nusselt numberN
u
x and nanoparticle Sherwood number N
S
h
x. These parameters characterize the wall heat and nanoparticle mass transfer rates, respectively.The local heat and nanoparticle mass fluxes from the wavy plate can be obtained from(13)
q
w
=
-
k
n
·
∇
T
+
n
·
q
r
,
q
n
p
=
-
D
B
n
·
∇
ϕ
,where n
=
-
δ
˙
/
1
+
δ
˙
2
,
1
/
1
+
δ
˙
2 is the unit normal vector to the wavy plate and q
r is the radiative heat flux given by q
r
=
(
-
4
σ
/
3
K
e
)
∇
T
4.The dimensionless local Nusselt numberN
u
x
=
x
q
w
/
k
(
T
w
-
T
∞
) and the nanoparticle Sherwood number N
S
h
x
=
x
q
n
p
/
D
B
(
ϕ
w
-
ϕ
∞
) are given by(14)
N
u
x
P
e
x
=
-
1
+
4
R
3
1
1
+
δ
˙
2
∂
θ
∂
η
η
=
0
,
N
S
h
x
P
e
x
=
-
1
1
+
δ
˙
2
∂
s
∂
η
η
=
0
.


## 3. Numerical Solution

Equations (11) along with the boundary conditions (12a) and (12b) were solved numerically using the successive linearization method (SLM) [31–33]. Using this method the nonlinear boundary layer equations reduce to a system of linear differential equations. The Chebyshev pseudo spectral method is then used to transform the iterative sequence of linearized differential equations into a system of linear algebraic equations which are converted into a matrix system. In this method we assume that the independent variables f
(
η
), θ
(
η
), and s
(
η
) can be expressed as(15)
f
η
=
f
i
η
+
∑
n
=
0
i
-
1
f
n
η
,
θ
η
=
θ
i
η
+
∑
n
=
0
i
-
1
θ
n
η
,
s
η
=
s
i
η
+
∑
n
=
0
i
-
1
s
n
η
,where f
i, θ
i, and s
i, (
i
=
1,2
,
3
,
…
) are unknown functions and f
n, θ
n, and s
n are the approximations which are obtained by recursively solving the linear part of the equation system that results from substituting (15) in (11) and neglecting the nonlinear terms in f
i, θ
i, and s
i
(16)
a
1
,
i
-
1
f
i
′
′
+
a
2
,
i
-
1
f
i
′
+
a
3
,
i
-
1
θ
i
′
+
a
4
,
i
-
1
s
i
′
=
r
1
,
i
-
1
,
b
1
,
i
-
1
θ
i
′
′
+
b
2
,
i
-
1
θ
i
′
+
b
3
,
i
-
1
s
i
′
+
b
4
,
i
-
1
f
i
=
r
2
,
i
-
1
,
s
i
′
′
+
c
1
,
i
-
1
s
i
′
+
c
2
,
i
-
1
θ
i
′
′
+
c
3
,
i
-
1
f
i
=
r
3
,
i
-
1
,where(17)
a
1
,
i
-
1
=
1
+
2
F
c
1
+
δ
˙
2
-
1
/
2
∑
n
=
0
i
-
1
f
n
′
,
a
2
,
i
-
1
=
2
F
c
1
+
δ
˙
2
-
1
/
2
∑
n
=
0
i
-
1
f
n
′
′
,
a
3
,
i
-
1
=
-
Δ
sin
⁡
A
+
δ
˙
cos
⁡
A
,
a
4
,
i
-
1
=
N
r
Δ
sin
⁡
A
+
δ
˙
cos
⁡
A
,
b
1
,
i
-
1
=
1
+
4
R
3
,
b
2
,
i
-
1
=
1
2
∑
n
=
0
i
-
1
f
n
+
N
b
∑
n
=
0
i
-
1
s
n
′
+
2
N
t
∑
n
=
0
i
-
1
θ
n
′
,
b
3
,
i
-
1
=
N
b
∑
n
=
0
i
-
1
θ
n
′
,
b
4
,
i
-
1
=
1
2
∑
n
=
0
i
-
1
θ
n
′
,
c
1
,
i
-
1
=
1
2
L
e
∑
n
=
0
i
-
1
f
n
,
c
2
,
i
-
1
=
N
t
N
b
,
c
3
,
i
-
1
=
1
2
L
e
∑
n
=
0
i
-
1
s
n
′
,
r
1
,
i
-
1
=
-
∑
n
=
0
i
-
1
f
n
′
′
-
2
F
c
1
+
δ
˙
2
-
1
/
2
∑
n
=
0
i
-
1
f
n
′
∑
n
=
0
i
-
1
f
n
′
′
+
Δ
sin
⁡
A
+
δ
˙
cos
⁡
A
∑
n
=
0
i
-
1
θ
n
′
-
Δ
sin
⁡
A
+
δ
˙
cos
⁡
A
N
r
∑
n
=
0
i
-
1
s
n
′
,
r
2
,
i
-
1
=
-
1
+
4
R
3
∑
n
=
0
i
-
1
θ
n
′
′
-
1
2
∑
n
=
0
i
-
1
f
n
∑
n
=
0
i
-
1
θ
n
′
-
N
b
∑
n
=
0
i
-
1
s
n
′
∑
n
=
0
i
-
1
θ
n
′
-
N
t
∑
n
=
0
i
-
1
θ
n
′
∑
n
=
0
i
-
1
θ
n
′
,
r
3
,
i
-
1
=
-
∑
n
=
0
i
-
1
s
n
′
′
-
1
2
∑
n
=
0
i
-
1
f
n
∑
n
=
0
i
-
1
s
n
-
N
t
N
b
∑
n
=
0
i
-
1
θ
n
′
′
.The boundary conditions reduce to(18)
f
i
0
=
f
i
′
∞
=
0
,
θ
i
0
=
θ
i
∞
=
s
i
0
=
s
i
∞
=
0
.The initial approximationsf
0
(
η
), θ
0
(
η
), and s
0
(
η
) are chosen such that they satisfy the boundary conditions (12a) and (12b) and are taken as f
0
(
η
)
=
1
+
η
-
e
-
η, θ
0
(
η
)
=
e
-
η, and s
0
(
η
)
=
e
-
η.The solutionsf
i, θ
i, and s
i (i
≥
1) are obtained by iteratively solving (16). The approximate solutions for f
(
η
), θ
(
η
), and s
(
η
) are then obtained as(19)
f
η
≈
∑
m
=
0
M
f
m
η
,
θ
η
≈
∑
m
=
0
M
θ
m
η
,
s
η
≈
∑
m
=
0
M
s
m
η
,where M is the order of SLM approximation. The linearized equations (16) are solved using the Chebyshev spectral collocation method [34]. In this method the unknown functions are approximated by the Chebyshev interpolating polynomials in such a way that they are collocated at the Gauss-Lobatto points defined as(20)
ξ
j
=
cos
⁡
π
j
N
,
j
=
0,1
,
2
,
…
,
N
,where N is the number of collocation points used. The physical region 0
,
∞ is transformed into the region [
-
1
,
1
] using the domain truncation technique in which the problem is solved on the interval [
0
,
L
] instead of 0
,
∞. This leads to the mapping(21)
η
L
=
ξ
+
1
2
,
-
1
≤
ξ
≤
1
,where L is a scaling parameter used to invoke the boundary condition at infinity. The functions f
i, θ
i, and S
i are approximated at the collocation points by(22)
f
i
ξ
=
∑
k
=
0
N
f
i
ξ
k
T
k
ξ
j
,
θ
i
ξ
=
∑
k
=
0
N
θ
i
ξ
k
T
k
ξ
j
,
s
i
ξ
=
∑
k
=
0
N
s
i
ξ
k
T
k
ξ
j
,
j
=
0,1
,
2
,
…
,
N
,where T
k is the k
t
h Chebyshev polynomial defined by(23)
T
k
ξ
=
c
o
s
k
cos
-
1
⁡
ξ
.The derivatives of the variables at the collocation points are represented as(24)
d
a
f
i
d
η
a
=
∑
k
=
0
N
D
k
j
a
f
i
ξ
k
,
d
a
θ
i
d
η
a
=
∑
k
=
0
N
D
k
j
a
θ
i
ξ
k
,
d
a
s
i
d
η
a
=
∑
k
=
0
N
D
k
j
a
s
i
ξ
k
,
j
=
0,1
,
2
,
…
,
N
,where a is the order of differentiation and D
=
(
2
/
L
)
D with D being the Chebyshev spectral differentiation matrix. Substituting (21)–(24) into (16) leads to the matrix equation(25)
A
i
-
1
X
i
=
R
i
-
1
,subject to the boundary conditions(26)
f
i
ξ
N
=
0
,
∑
k
=
0
N
D
0
k
f
i
ξ
k
=
0
,
θ
i
ξ
N
=
θ
i
ξ
0
=
s
i
ξ
N
=
s
i
ξ
0
=
0
.In (25), A
i
-
1 is a (
3
N
+
3
)
×
(
3
N
+
3
) square matrix and X
i and R
i
-
1 are (
3
N
+
3
)
×
1 column vectors defined by(27)
A
i
-
1
=
A
11
A
12
A
13
A
21
A
22
A
23
A
31
A
32
A
33
,
X
i
=
F
i
Θ
i
Φ
i
,
R
i
-
1
=
r
1
,
i
-
1
r
2
,
i
-
1
r
3
,
i
-
1
,where(28)
F
i
=
f
i
ξ
0
,
f
i
ξ
1
,
…
,
f
i
ξ
N
-
1
,
f
i
ξ
N
T
,
Θ
i
=
θ
i
ξ
0
,
θ
i
ξ
1
,
…
,
θ
i
ξ
N
-
1
,
θ
i
ξ
N
T
,
Φ
i
=
s
i
ξ
0
,
s
i
ξ
1
,
…
,
s
i
ξ
N
-
1
,
s
i
ξ
N
T
,
r
1
,
i
-
1
=
r
1
,
i
-
1
ξ
0
,
r
1
,
i
-
1
ξ
1
,
…
,
r
1
,
i
-
1
ξ
N
-
1
,
r
1
,
i
-
1
ξ
N
T
,
r
2
,
i
-
1
=
r
2
,
i
-
1
ξ
0
,
r
2
,
i
-
1
ξ
1
,
…
,
r
2
,
i
-
1
ξ
N
-
1
,
r
2
,
i
-
1
ξ
N
T
,
r
3
,
i
-
1
=
r
3
,
i
-
1
ξ
0
,
r
3
,
i
-
1
ξ
1
,
…
,
r
3
,
i
-
1
ξ
N
-
1
,
r
3
,
i
-
1
ξ
N
T
,
A
11
=
a
1
,
i
-
1
D
2
+
a
2
,
i
-
1
D
,
A
12
=
a
3
,
i
-
1
D
,
A
13
=
a
4
,
i
-
1
D
,
A
21
=
b
4
,
i
-
1
I
,
A
22
=
b
1
,
i
-
1
D
2
+
b
2
,
i
-
1
D
,
A
23
=
b
3
,
i
-
1
D
,
A
31
=
c
3
,
i
-
1
I
,
A
32
=
c
2
,
i
-
1
D
2
,
A
33
=
D
2
+
c
1
,
i
-
1
D
.Herea
k
,
i
-
1, b
l
,
i
-
1, and c
l
,
i
-
1
(
k
=
1,2 and l
=
1,2
,
3
) are diagonal matrices of size (
N
+
1
)
×
(
N
+
1
) and I is an identity matrix of size (
N
+
1
)
×
(
N
+
1
). After modifying the matrix system (25) to incorporate boundary conditions (26), the solution is obtained as(29)
X
i
=
A
i
-
1
-
1
R
i
-
1
.


## 4. Results and Discussion

In order to validate the numerical procedure generated, the results of the present problem have been compared with work of Ping [35] as a special case by taking A
=
π
/
2, a
=
0, ξ
=
0, N
r
=
0, N
t
=
0, N
b
→
0, L
e
=
0, F
c
=
0, and R
=
0 and we found that they are in good agreement, as shown in Table 1. Numerical solutions for the dimensionless velocity, temperature and nanoparticle volume fraction functions, and heat and nanoparticle mass transfer rates are computed and presented graphically in Figures 2–11. The effects of radiation R, angle of inclination A, Brownian motion parameter N
b, thermophoretic parameter N
t, non-Darcy parameter F
c, and amplitude a of the wavy surface have been discussed.Table 1
Comparison of values of-
θ
′
(
0
) for aiding and opposing flow by the present method and Ping [35] for fixed values of A
=
π
/
2, a
=
0, ξ
=
0, N
r
=
0, N
t
=
0, N
b
→
0, Le
=
0, and R
=
0.
Aiding flow
Opposing flow
Δ
Ping [35]
Present
Δ
Ping [35]
Present
0
0.5641
0.56415775
−0.2
0.5269
0.52691089
0.5
0.6473
0.64736510
−0.4
0.4865
0.48653284
1
0.7205
0.72055401
−0.6
0.442
0.44202064
3
0.9574
0.95744512
−0.8
0.3916
0.39166292
10
1.516
1.51623967
−1.0
0.332
0.33202116
20
2.066
2.066Figure 2
Variation of (a) velocity, (b) temperature, and (c) nanoparticle volume fraction profiles with the wave amplitude (a).
(a)
(b)
(c)Figure 3
Variation of (a) velocity, (b) temperature, and (c) nanoparticle volume fraction profiles with angle of inclination (A).
(a)
(b)
(c)Figure 4
Variation of (a) velocity, (b) temperature, and (c) nanoparticle volume fraction profiles with non-Darcy parameter (F
c).
(a)
(b)
(c)Figure 5
Variation of (a) velocity, (b) temperature, and (c) nanoparticle volume fraction profiles with radiation parameter (R).
(a)
(b)
(c)Figure 6
Effect of the wave amplitude (a) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure 7
Effect of the angle of inclination (A) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure 8
Effect of the Brownian motion parameter (N
b) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure 9
Effect of the themophoresis parameter (N
t) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure 10
Effect of the non-Darcy parameter (F
c) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure 11
Effect of the radiation parameter (R) on the heat and nanoparticle mass transfer rates.
(a)
(b)Figure2 shows the effect of the amplitude of the wavy surface on velocity, temperature, and nanoparticle volume fraction distributions. It is observed that as a increases, velocity increases near the plate and decreases away from the plate, whereas the temperature and nanoparticle volume fraction decrease for aiding flow and the reverse trend is observed in the case of opposing flow.The effect of the angle of inclination of the wavy surface on velocity, temperature, and nanoparticle volume fraction is plotted in Figure3. The similarity equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting A
=
0
° and A
=
9
0
°, respectively. It is noted from Figures 3(a)–3(c) that as A increases, the velocity increases near the plate and decreases away from the plate for aiding flow whereas the velocity decreases near the plate and increases away from the plate for opposing flow. But the temperature and nanoparticle volume fraction decrease within the boundary layer region for the aiding flow and increase for opposing flow. When the surface is vertical, the smallest temperature and nanoparticle volume fraction distributions are observed for aiding flow whereas largest temperature and nanoparticle volume fraction distributions are observed for opposing flow. While for the horizontal surface, largest temperature and nanoparticle volume fraction distributions are observed for aiding flow whereas smallest temperature and nanoparticle volume fraction distributions are observed for opposing flow.Figure4 displays the effect of non-Darcy parameter on the dimensionless velocity, temperature, and nanoparticle fraction. Figure 4(a) shows that an increase in non-Darcy parameter reduces velocity for aiding flow. Since F
c represent the inertial drag, thus an increase in the non-Darcy parameter increases the resistance to the flow and so a decrease in the fluid velocity is ensured. Here F
c
=
0 represents the case where the flow is Darcian. The velocity is maximum in this case due to the total absence of inertial drag. In case of opposing flow the velocity is increasing with the increase in the non-Darcy parameter. An increase in F
c increases temperature values for aiding flow, since as the fluid is decelerated, energy is dissipated as heat and serves to increase temperatures. As such the temperature is minimised for the lowest value of F
c and maximised for the highest value of F
c as shown in Figure 4(b). Further it is noticed that increase in F
c decreases the temperature for opposing flow. From Figure 4(c), it is observed that the nanoparticle volume fraction increases for aiding flow and decreases in the case of opposing flow.The effect of radiationR on velocity, temperature, and nanoparticle volume fraction is shown in Figure 5. Figure 5(a) reveals that an increase in the radiation parameter increases the velocity for aiding flow and reduces the same for the case of opposing flow. Increase in radiation parameter enhances the temperature for both aiding and opposing flows as shown in Figure 5(b). We observe a negligible effect of radiation on nanoparticle volume fraction as shown in Figure 5(c). These results can be explained by the fact that an increase in the radiation parameter (R
=
4
σ
T
∞
3
/
K
K
e) for given K and T
∞ means a decrease in the Rosseland radiation absorptivity K
e. Hence, the divergence of radiative heat flux q
r increases as K
e decreases. Therefore, the rate of radiative heat transferred to the fluid increases, and consequently the fluid temperature and simultaneously the velocity of the fluid also increase.Figure6 shows the effect of wave amplitude on the Nusselt and nanoparticle Sherwood number. It reveals that an increase in the wave amplitude reduces the local heat and nanoparticle mass transfer for both aiding and opposing flows. In general, we conclude that increasing the wave amplitude makes the surface more roughened.The variation of heat and nanoparticle mass transfer rates for various values of the angle of inclinationA is displayed in Figure 7. This figure shows that increasing the angle of inclination increases the buoyancy force and assists the flow, leading to an increase in the heat and nanoparticle mass transfer rates for aiding flow whereas reverse trend is observed in the case of opposing flow. The minimum values of the dimensionless heat and nanoparticle mass transfer rates are observed when the surface is horizontal; in which case, the buoyancy force is at its maximum. Therefore, the heat and nanoparticle mass transfer rates are at a lower level when this effect is considered (A
≠
0).The effect of Brownian motion parameterN
b on the heat and nanoparticle mass transfer rates is presented in Figure 8. Figure 8(a) depicts that the dimensionless heat transfer rate decreases with the increase in the Brownian motion parameter for both aiding and opposing flows. An increase in the value of Brownian motion parameter enhances the nanoparticle volume fraction transfer rate for both aiding and opposing flows, as shown in Figure 8(b). Figure 9 depicts the streamwise distribution of Nusselt and Sherwood numbers for different values of thermophoretic parameter N
t. It is revealed that the heat and nanoparticle mass transfer rate reduces with increase in the thermophoresis parameter for both aiding and opposing flows. Brownian motion is proportional to the volumetric fraction of nanoparticles in the direction from high to low concentration, whereas the thermophoresis is proportional to the temperature gradient from hot to cold. Hence, we conclude that the effect of the combination of Brownian motion and thermophoresis is to reduce the value of Nusselt number.Figure10 displays the effect of non-Darcy parameter F
c on the Nusselt and nanoparticle Sherwood number. It is observed that, in the case of aiding flow, the heat and nanoparticle mass transfer reduces with increase in the value of non-Darcy parameter and a reverse trend is observed in the case of opposing flow. The effect of radiation R on the heat and nanoparticle mass transfer is illustrated in Figure 11. It is noticed from Figure 11(a) that increase in the radiation parameter increases the heat transfer rates for both aiding and opposing flows. And a negligible effect of radiation on nanoparticle mass transfer is shown in Figure 11(b).


## 5. Conclusions

In this paper, the mixed convection in a nanofluid along an inclined wavy surface embedded in a non-Darcy porous medium with radiation effect is analyzed. The governing equations are transformed into a set of ordinary differential equations using the appropriate similarity transformation and then are linearized using the successive linearization method. These linearized equations are solved using Chebyshev spectral collocation method. The main findings are summarized as follows:(1)
An increase in the wave amplitude enhances the velocity for aiding flow but reduces the same for opposing flow whereas the temperature and nanoparticle volume fraction reduce for aiding flow and enhances the same for opposing flow. The local heat and nanoparticle mass transfer rates are increased for both aiding and opposing flows.
(2)
The effect of angle of inclination is to increase the velocity, local heat, and nanoparticle mass transfer for aiding flow and a reverse trend is observed in the case of opposing flow. Moreover, the temperature and nanoparticle volume fraction reduce for aiding flow and a reverse trend is seen for opposing flow.
(3)
The effect of non-Darcy parameter is to reduce the velocity, local heat, and nanoparticle mass transfer for aiding flow and to increase in the case of opposing flow while the temperature and nanoparticle volume fraction increase for aiding flow and reduce for opposing flow.
(4)
The influence of radiation is to enhance the velocity and temperature for aiding flow and reduce the same in the case of opposing flow while the local heat transfer enhances for both aiding and opposing flows and a negligible effect on nanoparticle volume fraction and nanoparticle mass transfer rates is observed.



---

*Source: 927508-2015-09-07.xml*