# Effects of Wall Shear Stress on MHD Conjugate Flow over an Inclined Plate in a Porous Medium with Ramped Wall Temperature

Authors: Arshad Khan; Ilyas Khan; Farhad Ali; Sharidan Shafie
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This study investigates the effects of an arbitrary wall shear stress on unsteady magnetohydrodynamic (MHD) flow of a Newtonian
fluid with conjugate effects of heat and mass transfer. The fluid is considered in a porous medium over an inclined plate with ramped temperature. The influence of thermal radiation in the energy equations is also considered. The coupled partial differential equations governing the flow are solved by using the Laplace transform technique. Exact solutions for velocity and temperature in case of both ramped and constant wall temperature as well as for concentration fields are obtained. It is found that velocity solutions are more general and can produce a huge number of exact solutions correlative to various fluid motions. Graphical results are provided for various embedded flow parameters and discussed in detail.


---

## Body

## 1. Introduction

The study of the natural convection heat transfer from different geometries of the surfaces has received much attention in recent years both analytically and experimentally. Amongst them a very little attention has been given to the problem of natural convection over an inclined plate as shown in Figure1. Free convection heat transfer from an inclined surface has important role in making engineering devices and natural environment. Natural convection over an inclined plate was first studied experimentally by Rich [1]. A solution for the boundary layer on a horizontal plate showing that if the plate is heated and faces downwards or is cooled and faces upwards was presented by Stewartson [2]. Free convection heat transfer from an isothermal plate with arbitrary inclination was investigated by [3]. Chen et al. [4] have obtained a numerical solution for the problem of natural convection over an inclined plate with variable surface temperature. Ganesan and Ekambavanan [5] solved the unsteady natural convection boundary layer flow over a semi-infinite inclined plate with the wall temperature varying as the axial coordinate using an implicit finite difference scheme. The results of an experimental study on natural convection between inclined plates were presented by [6]. MHD natural convection from a nonisothermal inclined surface with multiple suction/injection slots embedded in a thermally stratified high porosity medium has been studied by [7]. A numerical solution of the transient free convection MHD flow of an incompressible viscous fluid past a semi-infinite inclined plate with variable surface heat and mass flux is presented by [8]. Unsteady free convection flow of water at its maximum density past a semi-infinite inclined plate with variable surface temperature of the plate was studied by [9]. An investigation deals with study of laminar natural convection flow of a viscous fluid over a semi-infinite flat plate inclined at a small angle to the horizontal with internal heat generation and variable viscosity was presented by [10]. Pressure effects on unsteady free convection and heat transfer flow of an incompressible fluid past a semi-infinite inclined plate with impulsive and uniformly accelerated motion were analyzed by [11]. Chemical reaction effects on MHD free convection flow through a porous medium bounded by an inclined surface were studied by [12]. An analysis to study the heat and mass transfer characteristics of natural convection flow along horizontal and inclined plates with variable surface temperature/concentration or heat/mass flux under the combined buoyancy effects of thermal and mass diffusion was presented by [13]. A steady two dimensional MHD free convection and mass transfer flow past an inclined semi-infinite surface in the presence of heat generation has been studied numerically in [14]. A two-dimensional steady MHD mixed convection and mass transfer flow over a semi-infinite porous inclined plate in the presence of thermal radiation with variable suction and thermophoresis has been analyzed numerically by [15]. Heat and mass transfer of an incompressible viscous fluid past a semi-infinite inclined surface with first-order homogeneous chemical reaction by Lie group analysis had been presented by [16]. Recently, the exact analysis of combined effects of radiation and chemical reaction on the MHD free convection flow of an electrically conducting incompressible viscous fluid over an inclined plate embedded in a porous medium was studied by [17]. The inherent irreversibility and thermal stability in a gravity driven temperature dependent variable viscosity thin liquid film along an inclined heated plate with convective cooling were investigated by Makinde [18]. The inherent irreversibility in hydromagnetic boundary layer flow of variable viscosity fluid over a semi-infinite flat plate under the influence of thermal radiation and Newtonian heating was analyzed by [19]. Saha et al. [20] investigated that a natural convection boundary layer adjacent to an inclined semi-infinite plate subject to a temperature boundary condition which follows a ramped function up until some specified time and then remains constant. Saha et al. [21] also performed a scaling analysis for the transient boundary layer established adjacent to an inclined flat plate following a ramp cooling boundary condition. Recently, Ismail et al. [22] conducted the combined effects of heat and mass transfer on unsteady MHD free convection flow in a porous medium past an infinite inclined plate with ramped wall temperature. Fetecau et al. [23] investigated free convection flow near a vertical plate that applies arbitrary shear stress to the fluid when the thermal radiation and porosity effects are taken into consideration. However, from the literature, it is found that no study has been presented to investigate the unsteady MHD conjugate flow of an incompressible viscous fluid in a porous medium past over an inclined plate with ramped temperature in the presence of radiation under the boundary condition of wall shear stress.Figure 1
Physical configuration of the problem.Therefore, the aim of the present investigation is to provide exact solutions for MHD conjugate flow of a Newtonian fluid past an infinite plate that applies arbitrary shear stress to the fluid. More exactly, we consider the inclined plate situated in the(
x
,
z
) plane of a Cartesian coordinate system O
x
y
z, the domain of the flow is the porous half-space y
>
0, and the arbitrary shear stress on the plate is given by f
(
t
)
/
μ, where f
(
t
)is an arbitrary function and μis the viscosity. Closed form solutions of the initial and boundary value problems that govern the flow are obtained by means of the Laplace transform. Some special cases are extracted from the general solutions together with some limiting solutions in the literature. The results for velocity, temperature, and concentration profiles are plotted graphically and discussed for the embedded flow parameters.


## 2. Mathematical Formulation

Let us consider the unsteady MHD free convection flow of an incompressible viscous fluid over an infinite inclined plate. The physical configuration of the problem is shown in Figure1. The x-axis is taken along the plate and the y-axis is taken normal to it. Initially, both the plate and fluid are at stationary conditions with the constant temperature T
∞ and concentration C
∞. After time t
=
0
+, the plate applies a time dependent shear stress f
(
t
) to the fluid along the x-axis. Meanwhile, the temperature of the plate is raised or lowered to T
∞
+
(
T
w
-
T
∞
)
(
t
/
t
0
)when t
≤
t
0, and, thereafter, for t
>
t
0, is maintained at constant temperature T
w and concentration is raised to C
w. The radiation term is also considered in the energy equation. However, the radiative heat flux is considered negligible in x-direction compared to y-direction. We assume that the flow is laminar and the fluid is grey absorbing-emitting radiation but no scattering medium. Using Boussinesq’s approximation and neglecting the viscous dissipation, the equations governing the flow are given by [22]
(1)
∂
u
∂
t
=
ν
∂
2
u
∂
y
2
+
g
β
T
(
T
-
T
∞
)
cos
⁡
(
α
)
+
g
β
C
(
C

-

C
∞
)
cos
⁡
(
α
)

-

ν
K
u

-

σ
B
0
2
ρ
u
;
y
,
t
>
0
,
(2)
ρ
C
p
∂
T
∂
t
=
k
∂
2
T
∂
y
2
-
∂
q
r
∂
y
y
,
t
>
0
,
(3)
∂
C
∂
t
=
D
∂
2
C
∂
y
2
y
,
t
>
0
,
where u, T, C, ν, ρ, g, β
T, β
C, K, σ, B
0, C
p, k, q
r, and D are the velocity of the fluid in x-direction, its temperature and concentration, the kinematic viscosity, the constant density, the gravitational acceleration, the heat transfer coefficient, the mass transfer coefficient, the permeability of the porous medium, the electric conductivity of the fluid, the applied magnetic field, the heat capacity at constant pressure, the thermal conductivity, the radiative heat flux, and mass diffusivity.The corresponding initial and boundary conditions are(4)
u
(
y
,
0
)
=
0
,
T
(
y
,
0
)
=
T
∞
,
C
(
y
,
0
)
=
C
∞
;
∀
y
⩾
0
,
∂
u
(
0
,
t
)
∂
y
=
f
(
t
)
μ
,
C
(
0
,
t
)
=
C
w
;
t
>
0
,
T
(
0
,
t
)
=
T
∞
+
(
T
w
-
T
∞
)
t
t
0
,
0
<
t
<
t
0
,
T
(
0
,
t
)
=
T
w
;
t
≥
t
0
,
u
(
∞
,
t
)
=
0
,
T
(
∞
,
t
)
=
T
∞
,
C
(
∞
,
t
)
=
C
∞
;
t
>
0
.
The radiation heat flux under Rosseland approximation for optically thick fluid [24–26] is given by
(5)
q
r
=
-
4
σ
*
3
k
R
∂
T
4
∂
y
,
where σ
* and k
R is the Stefan-Boltzman constant and the mean absorption coefficient, respectively. We can see from (5) that the radiation term is nonlinear. Recently David Maxim Gururaj and Anjali Devi [27] used nonlinear radiation effects and studied MHD boundary layer flow with forced convection past a nonlinearly stretching surface with variable temperature. Therefore, we follow David Maxim Gururaj and Anjali Devi [27] and assume that the temperature differences within the flow are sufficiently small; that is, the difference between the fluid temperature and the free stream temperature is negligible, so that (5) can be linearized by expanding T into the Taylor series about T
∞, which after neglecting higher order terms takes the form
(6)
T
4
≈
4
T
∞
3
T
-
3
T
∞
4
.
Substituting (6) into (5) and then putting the obtained result in (2), we get
(7)
Pr
∂
T
∂
t
=
ν
(
1
+
N
r
)
∂
2
T
∂
y
2
;
y
,
t
>
0
,
where Pr, ν, and N
r are defined by
(8)
Pr
=
μ
C
p
k
,
ν
=
μ
ρ
,
N
r
=
16
σ
T
∞
3
3
k
k
R
.
By introducing the following dimensionless variables
(9)
u
*
=
u
t
0
ν
,
T
*
=
T
-
T
∞
T
w
-
T
∞
,
C
*
=
C
-
C
∞
C
w
-
C
∞
,
y
*
=
y
ν
t
0
,
t
*
=
t
t
0
,
f
*
(
t
*
)
=
t
0
μ
f
(
t
0
t
*
)
,
into (1), (3) and (7) and eleminating the star notations, we get
(10)
∂
u
∂
t
=
∂
2
u
∂
y
2
+
Gr
T
cos
⁡
α
+
Gm
C
cos
⁡
α
-
K
p
u
-
M
u
,
(11)
P
r
eff
∂
T
∂
t
=
∂
2
T
∂
y
2
,
(12)
∂
C
∂
t
=
1
Sc
∂
2
C
∂
y
2
,
where P
r
eff
=
Pr
/
(
1
+
N
r
) is the effective Prandtl number [24, (
10
)] and
(13)
Gr
=
g
β
T
(
T
w
-
T
∞
)
ν
U
0
3
,
Gm
=
g
β
c
(
C
w
-
C
∞
)
ν
U
0
3
,
M
=
σ
B
0
2
t
0
ρ
,
Sc
=
ν
D
,
K
p
=
ν
t
0
K
,
t
0
=
ν
U
0
2
,
are the Grashof number, modified Grashof number, magnetic parameter, Schmidt number, the inverse permeability parameter for the porous medium, and the characteristic time, respectively.The corresponding dimensionless initial and boundary conditions are(14)
u
(
y
,
0
)
=
0
,
T
(
y
,
0
)
=
0
,
C
(
y
,
0
)
=
0
;
∀
y
≥
0
,
∂
u
∂
y
|
y
=
0
=
f
(
t
)
,
T
(
0
,
t
)
=
t
;
0
<
t
≤
1
,
T
(
0
,
t
)
=
1
;
t
>
1
,
C
(
0
,
t
)
=
1
,
C
(
∞
,
t
)
=
0
,
T
(
∞
,
t
)
=
0
,
u
(
∞
,
t
)
=
0
;
t
>
0
.


## 3. Solution of the Problem

In order to solve (10)–(12) under conditions (14), we use the Laplace transform technique and get the following differential equations:
(15)
q
u
-
(
y
,
q
)
=
∂
2
u
-
(
y
,
q
)
∂
y
2
+
Gr
T
-
(
y
,
q
)
cos
⁡
(
α
)
+
Gm
C
-
(
y
,
q
)
cos
⁡
(
α
)
-
K
p
u
-
(
y
,
q
)
-
M
u
-
(
y
,
q
)
,
(16)
T
-
(
y
,
q
)
=
1
P
r
eff
q
∂
2
T
-
(
y
,
q
)
∂
y
2
,
(17)
C
-
(
y
,
q
)
=
1
Sc
q
∂
2
C
-
(
y
,
q
)
∂
y
2
,
with boundary conditions
(18)
C
-
(
∞
,
q
)
=
0
,
C
-
(
0
,
q
)
=
1
q
,
T
-
(
∞
,
q
)
=
0
,
u
-
(
∞
,
q
)
=
0
,
∂
u
-
(
y
,
q
)
∂
y
|
y
=
0
=
F
(
q
)
,
T
(
0
,
q
)
=
1
-
e
-
q
q
2
.
Solving (16) in view of (18), we get
(19)
T
-
(
y
,
q
)
=
1
q
2
e
-
y
q
P
r
eff
-
e
-
q
q
2
e
-
y
q
P
r
eff
,
which upon inverse Laplace transform gives
(20)
T
(
y
,
t
)
=
f
(
y
,
t
)
-
f
(
y
,
t
-
1
)
H
(
t
-
1
)
,
where
(21)
f
(
y
,
t
)
=
(
P
r
eff
y
2
2
+
t
)
erf
⁡
c
(
P
r
eff
y
2
t
)
-
P
r
eff
t
π
y
exp
⁡
(
-
P
r
eff
y
2
4
t
)
,
(22)
∂
T
(
y
,
t
)
∂
y
|
y
=
0
=
2
P
r
eff
π
(
t
-
t
-
1
H
(
t
-
1
)
)
,
is the corresponding heat transfer rate also known as Nusselt number. Here erf
⁡
(
·
) and erf
⁡
c
(
·
)denote the error function and complementary error function of Gauss.Solution of (17) using boundary conditions from (18) yields
(23)
C
(
y
,
q
)
=
1
q
e
-
y
Sc
q
,
which upon inverse Laplace transform gives
(24)
C
(
y
,
t
)
=
erf
⁡
c
(
y
Sc
2
t
)
,
(25)
∂
C
(
y
,
t
)
∂
y
|
y
=
0
=
-
Sc
π
t
,
that is the corresponding mass transfer rate also known as Sherwood number.The solution of (15) under boundary conditions (18) gives
(26)
u
-
(
y
,
q
)
=
a
1
q
q
2
(
q
-
a
2
)
q
+
H
1
e
-
y
q
+
H
1
-
a
1
q
e
-
q
q
2
(
q

-

a
2
)
q

+

H
1
e
-
y
q
+
H
1

-

F
(
q
)
q

+

H
1
e
-
y
q
+
H
1
-
a
3
q
2
(
q
-
a
2
)
e
-
y
q
P
r
eff
+
a
3
e
-
q
q
2
(
q
-
a
2
)
e
-
y
q
P
r
eff
+
a
4
q
q
(
q

-

a
5
)
q

+

H
1
e
-
y
q
+
H
1

-

a
6
q
(
q
-
a
5
)
e
-
y
q
Sc
,
which upon inverse Laplace transform gives
(27)
u
(
y
,
t
)
=
u
c
(
y
,
t
)
+
u
m
(
y
,
t
)
,
where
(28)
u
c
(
y
,
t
)
=
a
1
∫
0
t
(
e
a
2
(
t
-
s
)
erf
⁡
(
a
2
(
t
-
s
)
)
(
a
2
)
3
/
2
-
2
t
-
s
π
a
2
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
[
a
1
a
2
π
∫
0
t
-
1
(
2
t
-
1
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
]
×
H
(
t
-
1
)
-
[
a
1
(
a
2
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
2
(
t
-
1
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
×
(
a
2
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
)
d
s
a
1
(
a
2
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
2
(
t
-
1
-
s
)
-
H
1
s
-
(
y
2
/
4
s
)
erf
⁡
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
5
(
t
-
s
)
erf
⁡
(
a
5
(
t
-
s
)
)
a
5
-
2
t
-
s
π
a
2
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
a
3
a
2
(
t
+
P
r
eff
y
2
2
)
erf
⁡
c
(
y
P
r
eff
2
t
)
-
a
3
a
2
y
P
r
eff
t
π
e
-
y
2
P
r
eff
/
4
t
+
a
3
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
)
-
a
3
e
a
2
t
+
y
P
r
eff
a
2
2
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
+
a
2
t
)
-
a
3
e
a
2
t
-
y
P
r
eff
a
2
2
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
-
a
2
t
)
-
a
3
a
2
(
(
t
-
1
)
+
P
r
eff
y
2
2
)
erf
⁡
c
(
y
P
r
eff
2
t
-
1
)
H
(
t
-
1
)
+
a
3
a
2
y
P
r
eff
t
-
1
π
e
-
y
2
P
r
eff
/
4
(
t
-
1
)
H
(
t
-
1
)
-
a
3
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
-
1
)
H
(
t
-
1
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
+
a
3
e
a
2
(
t
-
1
)
+
y
P
r
eff
a
2
2
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
-
1
+
a
2
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
3
e
a
2
(
t
-
1
)
-
y
P
r
eff
a
2
2
a
2
2
erf
⁡
c
(
y
P
r
eff
2
t
-
1
-
a
2
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
5
t
-
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
-
a
6
e
a
5
t
+
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
,
(29)
u
m
(
y
,
t
)
=
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
,
correspond to the convective and mechanical parts of velocity.It is noted from (20) and (28) that T
(
y
,
t
) is valid for all positive values of P
r
eff while the u
c
(
y
,
t
) is not valid for Pr
eff
=
1. Therefore, to get u
c
(
y
,
t
) when the effective Prandtl number is not equal to one, we make P
r
eff
=
1 into (11), use a similar procedure as discussed above, and obtain
(30)
u
-
(
y
,
q
)
=
-
a
14
q
3
/
2
q
+
H
1
e
-
y
q
+
H
1
+
a
14
e
-
q
q
3
/
2
q
+
H
1
e
-
y
q
+
H
1
-
F
(
q
)
q
+
H
1
e
-
y
q
+
H
1
+
a
14
q
2
e
-
y
q
-
a
14
e
-
q
q
2
e
-
y
q
+
a
4
q
q
(
q

-

a
5
)
q

+

H
1
e
-
y
q
+
H
1

-

a
6
q
(
q
-
a
5
)
e
-
y
q
Sc
.
By taking inverse Laplace transform, we find that
(31)
u
(
y
,
t
)
=
-
2
a
14
π
∫
0
t
‍
t
-
s
e
-
H
1
s
-
y
2
/
4
s
s
d
s
+
(
2
a
14
π
∫
0
t
-
1
‍
t
-
1
-
s
e
-
H
1
s
-
y
2
/
4
s
s
d
s
)
H
(
t
-
1
)
+
a
4
∫
0
t
(
e
a
5
(
t
-
s
)
erf
⁡
(
a
5
(
t
-
s
)
)
a
5
-
2
t
-
s
π
a
2
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
a
14
[
(
t
+
y
2
2
)
erf
⁡
c
(
y
2
t
)
-
y
t
π
e
-
y
2
/
4
t
]
-
a
6
e
a
5
t
-
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
-
a
14
[
(
t
-
1
+
y
2
2
)
erf
⁡
c
(
y
2
t
-
1
)
-
y
t
-
1
π
e
-
y
2
/
4
(
t
-
1
)
]
H
(
t
-
1
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
-
a
6
e
a
5
t
+
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
,
where
(32)
a
1
=
cos
⁡
α
Gr
P
r
eff
P
r
eff
-
1
,
a
2

=

H
1
P
r
eff

-

1
,
a
3

=

Grcos
α
P
r
eff

-

1
,
a
4
=
cos
⁡
α
Gm
Sc
Sc
-
1
,
a
5
=
H
1
Sc
-
1
,
a
6
=
cos
⁡
α
Gm
Sc
-
1
,
a
7
=
cos
⁡
α
Gr
Pr
Pr
-
1
,
a
8
=
H
1
Pr
-
1
,
a
9
=
cos
⁡
α
Gr
Pr
-
1
,
a
10
=
K
p
P
r
eff
-
1
,
a
11
=
K
p
Sc
-
1
,
a
12
=
M
P
r
eff
-
1
,
a
13
=
M
Sc
-
1
,
a
14
=
cos
⁡
α
Gr
H
1
,
H
1
=
K
p
+
M
.


## 4. Plate with Constant Temperature

Equations (20) and (27) give analytical expressions for the temperature and velocity near an inclined plate with ramped temperature. In order to highlight the effect of the ramped temperature distribution of the boundary on the flow, it is important to compare such a flow with the one near a plate with constant temperature. It can be shown that the temperature, rate of heat transfer, and velocity for the flow near an isothermal plate are
(33)
T
(
y
,
t
)
=
erf
⁡
c
(
y
P
r
eff
2
t
)
,
(34)
∂
T
(
0
,
t
)
∂
y
=
-
P
r
eff
π
t
,
(35)
u
c
(
y
,
t
)
=
a
1
π
a
2
∫
0
t
e
a
2
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
2
(
t
-
s
)
)
s
d
s
+
a
4
π
a
5
∫
0
t
e
a
5
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
5
(
t
-
s
)
)
s
d
s
-
a
3
2
a
2
e
a
2
t
+
y
a
2
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
+
a
2
t
)
-
a
3
2
a
2
e
a
2
t
-
y
a
2
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
-
a
2
t
)
-
a
6
2
a
5
e
a
5
t
-
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
+
a
3
a
2
erf
⁡
c
(
y
Pr
eff
2
t
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
5
e
a
5
t
+
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
,
(36)
u
m
(
y
,
t
)
=
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
.
As previously, (36) is not valid for P
r
eff
=
1. Therefore we calculate separately solution for velocity by taking P
r
eff
=
1 into (11) and finally get
(37)
u
(
y
,
t
)
=
a
14
erf
⁡
c
(
y
2
t
)
-
a
14
π
∫
0
t
‍
e
-
H
1
s
-
y
2
/
4
s
(
t
-
s
)
s
d
s
+
a
4
π
a
5
∫
0
t
‍
e
a
5
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
5
(
t
-
s
)
)
s
d
s
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
-
a
6
2
a
5
e
a
5
t
+
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
-
a
6
2
a
5
e
a
5
t
-
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
.


## 5. Limiting Cases

In this section, we discuss few limiting cases of our general solutions.


### 5.1. Solution in the Absence of Porous Effects for Ramped and Constant Wall Temperature (K
p
→
0)

Consider(38)
u
(
y
,
t
)
=
a
1
∫
0
t
‍
(
e
a
12
(
t
-
s
)
erf
⁡
(
a
12
(
t
-
s
)
)
(
a
12
)
3
/
2
-
2
t
-
s
π
a
12
)
×
e
-
M
s
-
y
2
/
4
s
π
s
d
s
+
[
a
1
a
12
π
∫
0
t
-
1
‍
(
2
t
-
1
-
s
)
e
-
M
s
-
y
2
/
4
s
s
d
s
]
H
(
t
-
1
)
-
[
a
1
(
a
12
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
12
(
t
-
1
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
×
(
a
12
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
)
d
s
a
1
(
a
12
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
12
(
t
-
1
-
s
)
-
M
s
-
(
y
2
/
4
s
)
erf
⁡
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
13
(
t
-
s
)
erf
⁡
(
a
13
(
t
-
s
)
)
a
13
-
2
t
-
s
π
a
12
)
×
e
-
M
s
-
y
2
/
4
s
π
s
d
s
+
a
3
a
12
(
t
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
M
s
-
y
2
/
4
s
s
d
s
-
a
3
a
12
y
Pr
eff
t
π
e
-
y
2
Pr
eff
/
4
t
+
a
3
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
)
-
a
3
e
a
12
t
+
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
+
a
12
t
)
-
a
3
e
a
12
t
-
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
a
12
t
)
-
a
3
a
12
(
(
t
-
1
)

+

Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t

-

1
)
+
a
3
a
12
y
Pr
eff
t
-
1
π
e
-
y
2
Pr
eff
/
4
(
t
-
1
)
H
(
t
-
1
)
-
a
6
e
a
13
t
-
y
a
13
Sc
2
a
13
erf
⁡
c
(
y
Sc
2
t

-

a
13
t
)
-
a
3
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t
-
1
)

+

a
6
a
13
erf
⁡
c
(
y
Sc
2
t
)
+
a
3
e
a
12
(
t
-
1
)
+
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

+

a
12
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
3
e
a
12
(
t
-
1
)
-
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
-
a
12
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
13
t
+
y
a
13
Sc
2
a
13
erf
⁡
c
(
y
Sc
2
t
+
a
13
t
)
,
u
(
y
,
t
)
=
a
1
π
a
12
∫
0
t
e
a
12
(
t
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
(
a
12
(
t
-
s
)
)
s
d
s
+
a
4
π
a
13
∫
0
t
e
a
13
(
t
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
(
a
13
(
t
-
s
)
)
s
d
s
-
a
6
2
a
13
e
a
13
t
-
y
a
13
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
13
t
)
-
a
3
2
a
12
e
a
12
t
+
y
a
12
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
+
a
12
t
)
-
a
3
2
a
12
e
a
12
t
-
y
a
12
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
-
a
12
t
)
+
a
3
a
12
erf
⁡
c
(
y
Pr
eff
2
t
)
+
a
6
a
13
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
13
e
a
13
t
+
y
a
13
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
13
t
)
.


### 5.2. Solution in the Absence of Thermal Radiation (N
r
→
0)

In the absence of thermal radiation, the corresponding solutions for ramped and constant wall temperature are directly obtained from the general solutions (20), (22), (27), and (33)–(36) by taking N
r
→
0 and replacing Preff by Pr; that is,
(39)
u
(
y
,
t
)
=
a
7
∫
0
t
(
e
a
8
(
t
-
s
)
erf
⁡
(
a
8
(
t
-
s
)
)
(
a
8
)
3
/
2
-
2
t
-
s
π
a
8
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
[
a
7
π
a
8
∫
0
t
-
1
‍
2
t
-
1
-
s
e
-
H
1
s
-
y
2
/
4
s
s
d
s
]
×
H
(
t
-
1
)
-
[
a
7
(
a
8
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
×
e
a
8
(
t
-
1
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
)
d
s
a
7
(
a
8
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
5
(
t
-
s
)
erf
⁡
(
a
5
(
t
-
s
)
)
a
5
-
2
t
-
s
π
a
8
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
a
9
a
8
(
t
+
Pr
y
2
2
)
erf
⁡
c
(
y
Pr
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
-
a
9
a
8
y
Pr
t
π
e
-
y
2
Pr
/
4
t
+
a
9
a
8
2
erf
⁡
c
(
y
Pr
2
t
)
-
a
9
e
a
8
t
+
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
+
a
8
t
)
-
a
9
e
a
8
t
-
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
a
8
t
)
-
a
9
a
8
(
(
t
-
1
)
+
Pr
y
2
2
)
erf
⁡
c
(
y
Pr
2
t
-
1
)
H
(
t
-
1
)
-
a
9
a
8
erf
⁡
c
(
y
Pr
2
t
-
1
)
H
(
t
-
1
)
+
a
9
a
8
y
Pr
t
-
1
π
e
-
y
2
Pr
/
4
(
t
-
1
)
H
(
t
-
1
)
+
a
9
e
a
8
(
t
-
1
)
+
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
1
+
a
8
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
9
e
a
12
(
t
-
1
)
-
y
Pr
a
12
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
1
-
a
8
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
5
t
+
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
e
a
5
t
-
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
,
T
(
y
,
t
)
=
f
1
(
y
,
t
)
-
f
1
(
y
,
t
-
1
)
H
(
t
-
1
)
,
where
(40)
f
1
(
y
,
t
)
=
(
Pr
y
2
2
+
t
)
erf
⁡
c
(
Pr
y
2
t
)
-
Pr
t
π
y
exp
⁡
(
-
Pr
y
2
4
t
)
,
∂
T
(
y
,
t
)
∂
y
|
y
=
0
=
2
Pr
π
(
t
-
t
-
1
H
(
t
-
1
)
)
,
u
(
y
,
t
)
=
a
7
π
a
8
∫
0
t
‍
e
a
8
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
8
(
t
-
s
)
)
s
d
s
+
a
4
π
a
5
∫
0
t
‍
e
a
5
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
5
(
t
-
s
)
)
s
d
s
-
a
9
2
a
8
e
a
8
t
+
y
a
8
Pr
erf
⁡
c
(
y
Pr
2
t
+
a
8
t
)
-
a
9
2
a
8
e
a
8
t
-
y
a
8
Pr
erf
⁡
c
(
y
Pr
2
t
-
a
8
t
)
+
a
9
a
8
erf
⁡
c
(
y
Pr
2
t
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
5
e
a
5
t
+
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
-
a
6
2
a
5
e
a
5
t
-
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
,
T
(
y
,
t
)
=
erf
⁡
c
(
y
Pr
2
t
)
,
∂
T
(
0
,
t
)
∂
y
=
-
Pr
π
t
.


### 5.3. Solutions in the Absence of Free Convection

Let us assume that the flow is caused only due to bounding plate and the corresponding buoyancy forces are zero equivalently it shows the absence of free convection due to the differences in temperature and mass gradients that is, the termsGr and Gm are zero. This shows that the convective parts of velocities are zero in both cases of ramped wall and constant temperature and the flow is only governed by the mechanical part of velocitiesgiven by (29) and (36).


### 5.4. Solutions in the Absence of Mechanical Effects

In this case we assume that the infinite plate is in static position at every time; that is, the functionf
(
t
) is zero for all values of tand the mechanical parts for both ramped and constant wall temperature are equivalently zero. In such a situation, the motion in the fluid is induced only due to the free convection. Therefore, the velocities of the fluid in both cases of ramped and constant wall temperature are only represented by their convective parts given by (28) and (36).


### 5.5. Solution in the Absence of Magnetic Parameter (M
→
0)

It is clear from (20) and (24) that the temperature and concentration distributions are not influenced by the magnetic parameter M, and the velocities with M
=
0for both ramped and constant wall temperature are given by
(41)
u
(
y
,
t
)
=
a
1
∫
0
t
‍
(
e
a
10
(
t
-
s
)
erf
⁡
(
a
10
(
t
-
s
)
)
(
a
10
)
3
/
2
-
2
t
-
s
π
a
10
)
×
e
-
K
p
s
-
y
2
/
4
s
π
s
d
s
+
[
a
1
π
a
10
∫
0
t
-
1
‍
2
t
-
1
-
s
e
-
K
p
s
-
y
2
/
4
s
s
d
s
]
H
(
t
-
1
)
-
[
a
1
(
a
10
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
×
e
a
10
(
t
-
1
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
)
d
s
a
1
(
a
10
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
11
(
t
-
s
)
erf
⁡
(
a
11
(
t
-
s
)
)
a
11
-
2
t
-
s
π
a
10
)
×
e
-
K
p
s
-
y
2
/
4
s
π
s
d
s
+
a
3
a
10
(
t
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
)
-
a
3
a
10
y
Pr
eff
t
π
e
-
y
2
Pr
eff
/
4
t
+
a
3
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
K
p
s
-
y
2
/
4
s
s
d
s
-
a
3
e
a
10
t
+
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
+
a
10
t
)
-
a
3
e
a
10
t
-
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
a
10
t
)
-
a
3
a
10
(
(
t
-
1
)
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
×
H
(
t
-
1
)
+
a
3
a
10
y
Pr
eff
t
-
1
π
e
-
y
2
Pr
eff
/
4
(
t
-
1
)
H
(
t
-
1
)
-
a
6
e
a
11
t
-
y
a
11
Sc
2
a
11
erf
⁡
c
(
y
Sc
2
t
-
a
11
t
)
-
a
3
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t
-
1
)

+

a
6
a
11
erf
⁡
c
(
y
Sc
2
t
)
+
a
3
e
a
10
(
t
-
1
)
+
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

+

a
10
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
3
e
a
10
(
t
-
1
)
-
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

-

a
10
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
11
t
+
y
a
11
Sc
2
a
11
erf
⁡
c
(
y
Sc
2
t
+
a
11
t
)
,
u
(
y
,
t
)
=
a
1
π
a
10
∫
0
t
‍
e
a
10
(
t
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
10
(
t
-
s
)
)
s
d
s
+
a
4
π
a
11
∫
0
t
‍
e
a
11
(
t
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
11
(
t
-
s
)
)
s
d
s
-
a
3
2
a
10
e
a
10
t
+
y
a
10
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
+
a
10
t
)
-
a
3
2
a
10
e
a
10
t
-
y
a
10
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
-
a
10
t
)
+
a
3
a
10
erf
⁡
c
(
y
Pr
eff
2
t
)
+
a
6
a
11
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
11
e
a
11
t
+
y
a
11
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
11
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
K
p
s
-
y
2
/
4
s
s
d
s
-
a
6
2
a
11
e
a
11
t
-
y
a
11
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
11
t
)
.


## 5.1. Solution in the Absence of Porous Effects for Ramped and Constant Wall Temperature (K
p
→
0)

Consider(38)
u
(
y
,
t
)
=
a
1
∫
0
t
‍
(
e
a
12
(
t
-
s
)
erf
⁡
(
a
12
(
t
-
s
)
)
(
a
12
)
3
/
2
-
2
t
-
s
π
a
12
)
×
e
-
M
s
-
y
2
/
4
s
π
s
d
s
+
[
a
1
a
12
π
∫
0
t
-
1
‍
(
2
t
-
1
-
s
)
e
-
M
s
-
y
2
/
4
s
s
d
s
]
H
(
t
-
1
)
-
[
a
1
(
a
12
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
12
(
t
-
1
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
×
(
a
12
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
)
d
s
a
1
(
a
12
)
3
/
2
π
∫
0
t
-
1
(
(
e
a
12
(
t
-
1
-
s
)
-
M
s
-
(
y
2
/
4
s
)
erf
⁡
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
13
(
t
-
s
)
erf
⁡
(
a
13
(
t
-
s
)
)
a
13
-
2
t
-
s
π
a
12
)
×
e
-
M
s
-
y
2
/
4
s
π
s
d
s
+
a
3
a
12
(
t
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
M
s
-
y
2
/
4
s
s
d
s
-
a
3
a
12
y
Pr
eff
t
π
e
-
y
2
Pr
eff
/
4
t
+
a
3
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
)
-
a
3
e
a
12
t
+
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
+
a
12
t
)
-
a
3
e
a
12
t
-
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
a
12
t
)
-
a
3
a
12
(
(
t
-
1
)

+

Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t

-

1
)
+
a
3
a
12
y
Pr
eff
t
-
1
π
e
-
y
2
Pr
eff
/
4
(
t
-
1
)
H
(
t
-
1
)
-
a
6
e
a
13
t
-
y
a
13
Sc
2
a
13
erf
⁡
c
(
y
Sc
2
t

-

a
13
t
)
-
a
3
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t
-
1
)

+

a
6
a
13
erf
⁡
c
(
y
Sc
2
t
)
+
a
3
e
a
12
(
t
-
1
)
+
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

+

a
12
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
3
e
a
12
(
t
-
1
)
-
y
Pr
eff
a
12
2
a
12
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
-
a
12
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
13
t
+
y
a
13
Sc
2
a
13
erf
⁡
c
(
y
Sc
2
t
+
a
13
t
)
,
u
(
y
,
t
)
=
a
1
π
a
12
∫
0
t
e
a
12
(
t
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
(
a
12
(
t
-
s
)
)
s
d
s
+
a
4
π
a
13
∫
0
t
e
a
13
(
t
-
s
)
-
M
s
-
y
2
/
4
s
erf
⁡
(
a
13
(
t
-
s
)
)
s
d
s
-
a
6
2
a
13
e
a
13
t
-
y
a
13
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
13
t
)
-
a
3
2
a
12
e
a
12
t
+
y
a
12
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
+
a
12
t
)
-
a
3
2
a
12
e
a
12
t
-
y
a
12
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
-
a
12
t
)
+
a
3
a
12
erf
⁡
c
(
y
Pr
eff
2
t
)
+
a
6
a
13
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
13
e
a
13
t
+
y
a
13
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
13
t
)
.


## 5.2. Solution in the Absence of Thermal Radiation (N
r
→
0)

In the absence of thermal radiation, the corresponding solutions for ramped and constant wall temperature are directly obtained from the general solutions (20), (22), (27), and (33)–(36) by taking N
r
→
0 and replacing Preff by Pr; that is,
(39)
u
(
y
,
t
)
=
a
7
∫
0
t
(
e
a
8
(
t
-
s
)
erf
⁡
(
a
8
(
t
-
s
)
)
(
a
8
)
3
/
2
-
2
t
-
s
π
a
8
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
[
a
7
π
a
8
∫
0
t
-
1
‍
2
t
-
1
-
s
e
-
H
1
s
-
y
2
/
4
s
s
d
s
]
×
H
(
t
-
1
)
-
[
a
7
(
a
8
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
×
e
a
8
(
t
-
1
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
)
d
s
a
7
(
a
8
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
8
(
t
-
1
-
s
)
)
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
5
(
t
-
s
)
erf
⁡
(
a
5
(
t
-
s
)
)
a
5
-
2
t
-
s
π
a
8
)
×
e
-
H
1
s
-
y
2
/
4
s
π
s
d
s
+
a
9
a
8
(
t
+
Pr
y
2
2
)
erf
⁡
c
(
y
Pr
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
-
a
9
a
8
y
Pr
t
π
e
-
y
2
Pr
/
4
t
+
a
9
a
8
2
erf
⁡
c
(
y
Pr
2
t
)
-
a
9
e
a
8
t
+
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
+
a
8
t
)
-
a
9
e
a
8
t
-
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
a
8
t
)
-
a
9
a
8
(
(
t
-
1
)
+
Pr
y
2
2
)
erf
⁡
c
(
y
Pr
2
t
-
1
)
H
(
t
-
1
)
-
a
9
a
8
erf
⁡
c
(
y
Pr
2
t
-
1
)
H
(
t
-
1
)
+
a
9
a
8
y
Pr
t
-
1
π
e
-
y
2
Pr
/
4
(
t
-
1
)
H
(
t
-
1
)
+
a
9
e
a
8
(
t
-
1
)
+
y
Pr
a
8
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
1
+
a
8
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
9
e
a
12
(
t
-
1
)
-
y
Pr
a
12
2
a
8
2
erf
⁡
c
(
y
Pr
2
t
-
1
-
a
8
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
5
t
+
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
e
a
5
t
-
y
a
5
Sc
2
a
5
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
,
T
(
y
,
t
)
=
f
1
(
y
,
t
)
-
f
1
(
y
,
t
-
1
)
H
(
t
-
1
)
,
where
(40)
f
1
(
y
,
t
)
=
(
Pr
y
2
2
+
t
)
erf
⁡
c
(
Pr
y
2
t
)
-
Pr
t
π
y
exp
⁡
(
-
Pr
y
2
4
t
)
,
∂
T
(
y
,
t
)
∂
y
|
y
=
0
=
2
Pr
π
(
t
-
t
-
1
H
(
t
-
1
)
)
,
u
(
y
,
t
)
=
a
7
π
a
8
∫
0
t
‍
e
a
8
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
8
(
t
-
s
)
)
s
d
s
+
a
4
π
a
5
∫
0
t
‍
e
a
5
(
t
-
s
)
-
H
1
s
-
y
2
/
4
s
erf
⁡
(
a
5
(
t
-
s
)
)
s
d
s
-
a
9
2
a
8
e
a
8
t
+
y
a
8
Pr
erf
⁡
c
(
y
Pr
2
t
+
a
8
t
)
-
a
9
2
a
8
e
a
8
t
-
y
a
8
Pr
erf
⁡
c
(
y
Pr
2
t
-
a
8
t
)
+
a
9
a
8
erf
⁡
c
(
y
Pr
2
t
)
+
a
6
a
5
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
5
e
a
5
t
+
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
5
t
)
-
a
6
2
a
5
e
a
5
t
-
y
a
5
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
5
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
H
1
s
-
y
2
/
4
s
s
d
s
,
T
(
y
,
t
)
=
erf
⁡
c
(
y
Pr
2
t
)
,
∂
T
(
0
,
t
)
∂
y
=
-
Pr
π
t
.


## 5.3. Solutions in the Absence of Free Convection

Let us assume that the flow is caused only due to bounding plate and the corresponding buoyancy forces are zero equivalently it shows the absence of free convection due to the differences in temperature and mass gradients that is, the termsGr and Gm are zero. This shows that the convective parts of velocities are zero in both cases of ramped wall and constant temperature and the flow is only governed by the mechanical part of velocitiesgiven by (29) and (36).


## 5.4. Solutions in the Absence of Mechanical Effects

In this case we assume that the infinite plate is in static position at every time; that is, the functionf
(
t
) is zero for all values of tand the mechanical parts for both ramped and constant wall temperature are equivalently zero. In such a situation, the motion in the fluid is induced only due to the free convection. Therefore, the velocities of the fluid in both cases of ramped and constant wall temperature are only represented by their convective parts given by (28) and (36).


## 5.5. Solution in the Absence of Magnetic Parameter (M
→
0)

It is clear from (20) and (24) that the temperature and concentration distributions are not influenced by the magnetic parameter M, and the velocities with M
=
0for both ramped and constant wall temperature are given by
(41)
u
(
y
,
t
)
=
a
1
∫
0
t
‍
(
e
a
10
(
t
-
s
)
erf
⁡
(
a
10
(
t
-
s
)
)
(
a
10
)
3
/
2
-
2
t
-
s
π
a
10
)
×
e
-
K
p
s
-
y
2
/
4
s
π
s
d
s
+
[
a
1
π
a
10
∫
0
t
-
1
‍
2
t
-
1
-
s
e
-
K
p
s
-
y
2
/
4
s
s
d
s
]
H
(
t
-
1
)
-
[
a
1
(
a
10
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
×
e
a
10
(
t
-
1
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
)
×
(
s
)
-
1
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
)
d
s
a
1
(
a
10
)
3
/
2
π
∫
0
t
-
1
(
(
erf
⁡
(
a
10
(
t
-
1
-
s
)
)
]
×
H
(
t
-
1
)
+
a
4
∫
0
t
‍
(
e
a
11
(
t
-
s
)
erf
⁡
(
a
11
(
t
-
s
)
)
a
11
-
2
t
-
s
π
a
10
)
×
e
-
K
p
s
-
y
2
/
4
s
π
s
d
s
+
a
3
a
10
(
t
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
)
-
a
3
a
10
y
Pr
eff
t
π
e
-
y
2
Pr
eff
/
4
t
+
a
3
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
K
p
s
-
y
2
/
4
s
s
d
s
-
a
3
e
a
10
t
+
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
+
a
10
t
)
-
a
3
e
a
10
t
-
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
a
10
t
)
-
a
3
a
10
(
(
t
-
1
)
+
Pr
eff
y
2
2
)
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
×
H
(
t
-
1
)
+
a
3
a
10
y
Pr
eff
t
-
1
π
e
-
y
2
Pr
eff
/
4
(
t
-
1
)
H
(
t
-
1
)
-
a
6
e
a
11
t
-
y
a
11
Sc
2
a
11
erf
⁡
c
(
y
Sc
2
t
-
a
11
t
)
-
a
3
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1
)
H
(
t
-
1
)

+

a
6
a
11
erf
⁡
c
(
y
Sc
2
t
)
+
a
3
e
a
10
(
t
-
1
)
+
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

+

a
10
(
t
-
1
)
)
×
H
(
t
-
1
)
+
a
3
e
a
10
(
t
-
1
)
-
y
Pr
eff
a
10
2
a
10
2
erf
⁡
c
(
y
Pr
eff
2
t
-
1

-

a
10
(
t
-
1
)
)
×
H
(
t
-
1
)
-
a
6
e
a
11
t
+
y
a
11
Sc
2
a
11
erf
⁡
c
(
y
Sc
2
t
+
a
11
t
)
,
u
(
y
,
t
)
=
a
1
π
a
10
∫
0
t
‍
e
a
10
(
t
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
10
(
t
-
s
)
)
s
d
s
+
a
4
π
a
11
∫
0
t
‍
e
a
11
(
t
-
s
)
-
K
p
s
-
y
2
/
4
s
erf
⁡
(
a
11
(
t
-
s
)
)
s
d
s
-
a
3
2
a
10
e
a
10
t
+
y
a
10
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
+
a
10
t
)
-
a
3
2
a
10
e
a
10
t
-
y
a
10
Pr
eff
erf
⁡
c
(
y
Pr
eff
2
t
-
a
10
t
)
+
a
3
a
10
erf
⁡
c
(
y
Pr
eff
2
t
)
+
a
6
a
11
erf
⁡
c
(
y
Sc
2
t
)
-
a
6
2
a
11
e
a
11
t
+
y
a
11
Sc
erf
⁡
c
(
y
Sc
2
t
+
a
11
t
)
-
1
π
∫
0
t
‍
f
(
t
-
s
)
e
-
K
p
s
-
y
2
/
4
s
s
d
s
-
a
6
2
a
11
e
a
11
t
-
y
a
11
Sc
erf
⁡
c
(
y
Sc
2
t
-
a
11
t
)
.


## 6. Special Cases

We noted that the solutions for velocity obtained in Section3 are more general. Therefore, we want to discuss some special cases of the present solutions together with some limiting solutions in order to know more about the physical insight of the problem. Hence, we discuss the following important special cases in the case of ramped wall temperature whose technical relevance is well known in the literature. Similarly we can discuss some special cases of constant wall temperature solutions.


### 6.1. Case-I:f
(
t
)
=
f
H
(
t
)

In this first case we take the arbitrary functionf
(
t
)
=
f
H
(
t
)
,where f
(
·
) is a dimensionless constant and H
(
·
)denotes the unit step function. After time t
=
0
,the infinite inclined plate applies a constant shear stress to the fluid. The convective part of the velocity remains unchanged while the mechanical part takes the following form:
(42)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
e
-
y
2
/
4
s
-
H
1
s
s
d
s
,
equivalently
(43)
u
m
(
y
,
t
)
=
-
f
H
1
e
-
y
H
1
+
2
f
π
∫
t
∞
‍
e
-
y
2
/
4
z
2
-
H
1
z
2
d
z
,
for K
p
≠
0, M
≠
0. Moreover, if we take M
=
0, (42) reduces to the form
(44)
u
m
(
y
,
t
)
=
-
f
K
p
e
-
y
K
p
+
2
f
π
∫
t
∞
e
-
y
2
/
4
z
2
-
K
p
z
2
d
z
,
which is equivalent to [28, Equation (
28
)] with the correction of K
p.Furthermore, in the absence of bothK
p
=
0and M
=
0, (42) is identical with [23, Equation (
23
)]
(45)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
e
-
y
2
/
4
s
s
d
s
.


### 6.2. Case-II:f
(
t
)
=
f
sin
(
ω
t
)

In the second case, we take the arbitrary function of the formf
(
t
)
=
f
sin
(
ω
t
)in which the plate applies an oscillating shear stress to the fluid. Here ω denotes the dimensionless frequency of the shear stress. As previously, the convective part of velocity remains the same whereas the mechanical part takes the form
(46)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
It can be further written as a sum of the steady-state and transient solutions
(47)
u
m
(
y
,
t
)
=
u
m
s
(
y
,
t
)
+
u
m
t
(
y
,
t
)
,
where
(48)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
‍
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
,
u
m
t
(
y
,
t
)
=
f
π
∫
t
∞
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
By taking M
=
0, the steady-state component reduces to [28, Equation (
35
)]
(49)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
K
p
s
s
d
s
.
In addition when K
p
=
0
,physically it corresponds to the absence of porous effects and (49) results in
(50)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
‍
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
s
d
s
,
which can be written in simplified form as
(51)
u
m
s
(
y
,
t
)
=
f
ω
exp
⁡
(
-
y
ω
2
)
cos
⁡
(
ω
t
-
y
ω
2
+
π
4
)
,
equivalent to [23, Equation (
33
)].


### 6.3. Case-III:f
(
t
)
=
f
t
a
(
a
>
0
)

In the final case, we takef
(
t
)
=
f
t
a
,in which the plate applies an accelerating shear stress to the fluid where the mechanical part takes the following form:
(52)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
The corresponding solution for M
=
0, namely,
(53)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
-
K
p
s
s
d
s
,
is identical with [28, Equation (
32
)].Additionally, if we takeK
p
=
0
,(53) yields
(54)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
s
d
s
.


## 6.1. Case-I:f
(
t
)
=
f
H
(
t
)

In this first case we take the arbitrary functionf
(
t
)
=
f
H
(
t
)
,where f
(
·
) is a dimensionless constant and H
(
·
)denotes the unit step function. After time t
=
0
,the infinite inclined plate applies a constant shear stress to the fluid. The convective part of the velocity remains unchanged while the mechanical part takes the following form:
(42)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
e
-
y
2
/
4
s
-
H
1
s
s
d
s
,
equivalently
(43)
u
m
(
y
,
t
)
=
-
f
H
1
e
-
y
H
1
+
2
f
π
∫
t
∞
‍
e
-
y
2
/
4
z
2
-
H
1
z
2
d
z
,
for K
p
≠
0, M
≠
0. Moreover, if we take M
=
0, (42) reduces to the form
(44)
u
m
(
y
,
t
)
=
-
f
K
p
e
-
y
K
p
+
2
f
π
∫
t
∞
e
-
y
2
/
4
z
2
-
K
p
z
2
d
z
,
which is equivalent to [28, Equation (
28
)] with the correction of K
p.Furthermore, in the absence of bothK
p
=
0and M
=
0, (42) is identical with [23, Equation (
23
)]
(45)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
e
-
y
2
/
4
s
s
d
s
.


## 6.2. Case-II:f
(
t
)
=
f
sin
(
ω
t
)

In the second case, we take the arbitrary function of the formf
(
t
)
=
f
sin
(
ω
t
)in which the plate applies an oscillating shear stress to the fluid. Here ω denotes the dimensionless frequency of the shear stress. As previously, the convective part of velocity remains the same whereas the mechanical part takes the form
(46)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
It can be further written as a sum of the steady-state and transient solutions
(47)
u
m
(
y
,
t
)
=
u
m
s
(
y
,
t
)
+
u
m
t
(
y
,
t
)
,
where
(48)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
‍
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
,
u
m
t
(
y
,
t
)
=
f
π
∫
t
∞
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
By taking M
=
0, the steady-state component reduces to [28, Equation (
35
)]
(49)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
-
K
p
s
s
d
s
.
In addition when K
p
=
0
,physically it corresponds to the absence of porous effects and (49) results in
(50)
u
m
s
(
y
,
t
)
=
-
f
π
∫
0
t
‍
sin
(
ω
t
-
ω
s
)
e
-
y
2
/
4
s
s
d
s
,
which can be written in simplified form as
(51)
u
m
s
(
y
,
t
)
=
f
ω
exp
⁡
(
-
y
ω
2
)
cos
⁡
(
ω
t
-
y
ω
2
+
π
4
)
,
equivalent to [23, Equation (
33
)].


## 6.3. Case-III:f
(
t
)
=
f
t
a
(
a
>
0
)

In the final case, we takef
(
t
)
=
f
t
a
,in which the plate applies an accelerating shear stress to the fluid where the mechanical part takes the following form:
(52)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
-
H
1
s
s
d
s
.
The corresponding solution for M
=
0, namely,
(53)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
-
K
p
s
s
d
s
,
is identical with [28, Equation (
32
)].Additionally, if we takeK
p
=
0
,(53) yields
(54)
u
m
(
y
,
t
)
=
-
f
π
∫
0
t
‍
(
t
-
s
)
a
e
-
y
2
/
4
s
s
d
s
.


## 7. Results and Discussion

In order to understand the physical aspects of the problem, the numerical results for velocity, temperature, and concentration are computed and plotted for various parameters of interest such as magnetic parameterM, porosity parameter K
p, effective Prandtl number P
r
eff
,Grashof number Gr, modified Grashof number Gm, dimensionless time t, Schmidt number Sc, and shear stress f. The graphs for velocity are shown in Figures 2–17, where t
=
1.2 corresponds to isothermal velocity and t
=
0.9 is for ramped velocity. Figures 18–21 are plotted to show the temperature variations for two types of boundary conditions, namely, ramped and constant wall temperatures. Furthermore, Figure 22 is displayed to show variations in fluid concentration. Figures 2 and 3 illustrate the influence of Grashof number Gr on the velocity. It is observed that velocity increases with increasing Gr. This implies that thermal buoyancy force tends to accelerate velocity for both ramped temperature and isothermal plates. In Figures 4 and 5, the velocity profiles for different values modified Grashof number Gm are shown. It is found that velocity increases on increasing Gm for both ramped temperature and isothermal plate. Further, it can be observed that the velocity and boundary layer thickness decrease along y with increasing distance from the leading edge. Moreover, from Figures 4 and 5, we observed that the amplitude of velocity in case of isothermal plate is greater and converges slowly as compared to ramped velocity. In Figures 6 and 7, the velocity profiles are shown for different values of Schmidt number Sc. Here the values of Sc are chosen 0.22, 0.60, and 0.96. to represent the presence of species by hydrogen, water vapor, and carbon dioxide respectively. It is observed that the velocity decreases with increasing Schmidt number. Physically, this refers to the phenomenon that increasing Schmidt number implies the dominance of the viscous forces over the diffusional effects. As a result, the flow will be therefore decelerated with a rise in Schmidt number. The velocity profiles for different values of magnetic parameter M are shown in Figures 8 and 9. The range of magnetic field is taken from 0 to 2. It is found that the velocity is decreasing with increasing values of M in both cases of ramped and isothermal plates. Physically, it is true due to the fact that increasing values of M causes the frictional force to increase which tends to resist the fluid flow, thus reducing its velocity. It is further observed that when the magnetic field imposed on the flow is zero (M
=
0), the MHD effect vanishes and the flow is termed as hydrodynamic flow. Physically, it is true due to the fact that increasing values of M causes the frictional force to increase which tends to resist the fluid flow, thus reducing its velocity. Figures 10 and 11 are plotted to see the difference between the ramped and isothermal plate velocities. The values of t
<
1 correspond to ramp velocity whereas t
>
1 is for isothermal plate. It is found that ramp velocity is less than isothermal plate and converges faster. Further velocity in both cases increases with increasing time. The effects of inverse permeability parameter K
p on the velocity profiles are presented in Figures 12 and 13. It is found that velocity decreases with increasing K
p in both cases of ramp and isothermal plate. Physically, it is due to the fact that increasing permeability of the porous medium increases the resistance and consequently velocity decreases. This observation is an excellent agreement with the previous study [28, Figure 3]. The effects of the shear stress f induced by the bounding plate on the nondimensional velocity profiles are shown in Figures 14 and 15. The velocity of fluid is found to decrease with increasing f in both cases of ramped velocity and isothermal plate. Graphical results to show the influence of the effective Prandtl number P
r
eff on velocity profiles are presented in Figures 16 and 17. It is observed that the velocity is a decreasing function with respect to P
r
eff. These graphical results are in accordance with [28, Figure 2]. The temperature variations against y for various values of effective Prandtl number are highlighted in Figures 18 and 19. The significant decrease of the temperature is found as a result of an increase of the effective Prandtl number. The fluid temperature decreases from maximum at the boundary to a minimum value as far from the plate in both cases of ramped and constant temperature. In Figures 20 and 21, we have shown the temperature variations for two types of boundary conditions ramped and constant plate temperatures. It is noted that the fluid temperature is greater in the case of isothermal plate than in the case of ramped temperature at the plate. This should be expected since, in the latter case, the heating of the fluid takes place more gradually than in the isothermal case [29]. Moreover, with increasing time, the temperature is found to increase in both cases of ramped and constant wall temperature. The concentration profiles for different values of Schmidt number Sc are shown in Figure 22. It is clear from this figure that the concentration profiles and the concentration boundary layer thickness decrease with increasing values of Sc. Physically, it is true, since increase of Sc means decrease of molecular diffusivity which results in a decrease of concentration boundary layer.Figure 2
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7,   t
=
1.2 and different values of Gr when the plate applies a constant shear stress f
=
-
1.Figure 3
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, t
=
0.9 and different values of Gr when the plate applies a constant shear stress f
=
-
1.Figure 4
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, t
=
1.2 and different values of Gm when the plate applies a constant shear stress f
=
-
1.Figure 5
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
0.9 and different values of Gm when the plate applies a constant shear stress f
=
-
1.Figure 6
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
1.2 and different values of Sc when the plate applies a constant shear stress f
=
-
1.Figure 7
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
0.9 and different values of Sc when the plate applies a constant shear stress f
=
-
1.Figure 8
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
1.2 and different values of M when the plate applies a constant shear stress f
=
-
1.Figure 9
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
0.9 and different values of M when the plate applies a constant shear stress f
=
-
1.Figure 10
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and K
p
=
0.7 and different values of t when the plate applies a constant shear stress f
=
-
1.Figure 11
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and K
p
=
0.7 and different values of t when the plate applies a constant shear stress f
=
-
1.Figure 12
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and t
=
1.2 and different values of K
p when the plate applies a constant shear stress f
=
-
1.Figure 13
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and t
=
0.9 and different values of K
p when the plate applies a constant shear stress f
=
-
1.Figure 14
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
1.2 and different values of constant shear stress f.Figure 15
Velocity profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
), K
p
=
0.7, and t
=
0.9 and different values of constant shear stress f.Figure 16
Velocity profiles forK
p
=
0.7 and t
=
1.2 and different values of P
r
eff when the plate applies a constant shear stress f
=
-
1.Figure 17
Velocity profiles forK
p
=
0.7 and t
=
0.9 and different values of P
r
eff when the plate applies a constant shear stress f
=
-
1.Figure 18
Temperature profile fort
=
1.2 and different values of P
r
eff.Figure 19
Temperature profile fort
=
0.9 and different values of P
r
eff.Figure 20
Temperature profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and different values of t.Figure 21
Temperature profiles forP
r
eff
=
0.350
(
N
r
=
1
,
Pr
=
0.7
) and different values of t.Figure 22
Concentration profiles fort
=
1.2 and different values of Sc.


## 8. Conclusions

The purpose of this work was to analyze the unsteady MHD free convection flow of an incompressible viscous fluid over an infinite inclined plate with ramped wall temperature and applies an arbitrary shear stress to the fluid. Exact solutions for velocity, temperature (for both cases of ramped and constant wall temperature), and concentration are obtained using the Laplace transform technique and expressed in terms of the complementary error function. They satisfy all imposed initial and boundary conditions. These solutions are plotted in various figures for different parameters of interest. The following conclusions are extracted from this study.(i)
The velocity of the fluidu
(
y
,
t
) can be written as a sum of its convective and thermal components.
(ii)
For the velocity solution in which the plate applies an oscillating shear stress to the fluidf
(
t
)
=
f
sin
(
ω
t
)
,the mechanical part can be further written as a sum of the steady-state and transient solutions.
(iii)
The concentration boundary layer thickness decreases with increasing values ofSc.
(iv)
The thermal boundary layer thickness in case of ramped wall temperature is less than isothermal wall temperature.
(v)
The magnetic parameterM develops shear resistance which retards the fluid whereas the inverse permeability parameter K
p enhances the fluid motion.



---

*Source: 861708-2014-04-02.xml*