# Multiple Solutions of Mixed Convective MHD Casson Fluid Flow in a Channel

Authors: Jawad Raza; Azizah Mohd Rohni; Zurni Omar
Journal: Journal of Applied Mathematics
 (2016)

---

## Abstract

A numerical investigation is made to determine the occurrence of the multiple solutions of MHD Casson fluid in a porous channel. Governing partial differential equation of the proposed problem converted into nonlinear ordinary differential equations by using similarity transformation. Numerical technique known as shooting method is used to investigate the existence of the multiple solutions for the variations of different parameters. Effects of physical parameters on velocity profile, temperature, concentration, and skin friction are presented in pictorial and tabulation representation.


---

## Body

## 1. Introduction

The fluid flow in a channel has abundant practical applications in industry from mathematical and engineering point of view [1–5]. In particular, those fluids in which the shear stresses are not linear proportional to the velocity gradient are characterized as non-Newtonian fluids which are of much interest among the researchers. Among the non-Newtonian fluids, Casson fluid has attracted more attention of researchers due to its applications in the fields of metallurgy, food processing, drilling operations, and bioengineering operations [6, 7]. Some more applications of Casson fluid can be seen in the manufacturing of pharmaceutical products, coal in water, china clay, paints, synthetic lubricants, and biological fluids such as synovial fluids, sewage sludge, jelly, tomato sauce, honey, soup, and blood due to its contents such as plasma, fibrinogen, and protein [8].Due to the novel application of MHD mixed convection flow in porous medium, in the field of industrial engineering, many researchers are attracted towards it. Design of MHD power generators, nuclear waste processing, and distribution of chemical waste control are some most conspicuous applications among all. The problem of laminar fully developed mixed convection flow in a channel was followed back to 1960 by Tao [9]. Recently, Fersadou et al. [10] investigated the problem of MHD mixed convection flow of nanofluid in a vertical porous channel numerically. Problem of mixed convection MHD flow of an Al2O3 water nanofluid in a channel with asymmetric heated walls was examined by Chen et al. [11]. Khan et al. [12] investigated numerically the problem of Blasius and Sakiadis flows of Casson fluid with viscous dissipation and convective boundary conditions. Umavathi and Sultana [13] talked about the mixed convective hydromagnetic flow of micropolar fluid in a vertical channel with related boundary states of third kind. Si et al. [14] utilized homotopy analysis method (HAM) to acquire the arrangement of the problem of micropolar fluid in a channel with heat and mass exchange impacts. The channel is thought to be permeable with extending/contracting dividers. Prakash and Muthtamilselvan [15] executed Crank-Nicolson limited contrast plan to take care of numerical issue of mixed convective MHD stream of micropolar fluid between two vertical permeable dividers with proper limit states of third kind. Recently, Raza et al. [16] examined the rotational effects of channel on the problem of nanofluid with shrinking channel walls.The Casson constitutive equation was derived by Casson [17] which shows that the rate of strain and stress relationship is nonlinear. Flow of Casson fluid between two rotating cylinders is studied by Eldabe and Salwa [18]. Attia and Sayed-Ahmed [19] considered the Couette flow of electrically conducting Casson fluid between parallel plates. Effect of mass transfer on MHD flow of Casson fluid was discussed analytically by Shehzad et al. [20]. Taylor’s series was employed in order to solve nonlinear differential equations. Explicit finite difference method of unsteady Casson fluid flow through parallel plates was investigated by Afikuzzaman et al. [21]. Recently, Reddy et al. [22] discussed the effects of Joule heating and Hall effects on free convection in an electrically conducting Casson fluid in a vertical channel in the presence of viscous dissipation. Analytical solutions were found with the help of homotopy analysis method (HAM) and compared with Adomian Decomposition Method (ADM). Walawander et al. [23], Batra and Jena [24], Sayed Ahmed and Attia [25], Kataria and Patel [26], and Das et al. [27] have reported the flow of Casson fluid under different flow regimes. None of the investigations cited above dealt with multiple solutions of Casson fluid in a channel.Motivated by the above-cited investigations, the aim of the present study is to investigate the multiple solutions of mixed convection flow of Casson fluid in a porous channel under the influence of the magnetic field. The effects of different parameters such as Reynold numberR, magnetic field M, Casson parameter β, and bouncy parameter λ on velocity and temperature profiles are discussed graphically and also in tabulation representation.


## 2. Formulation of the Problem

Consider the steady, incompressible MHD flow of Casson fluid in a channel. Thex
-axis is along the centerline of the channel, parallel to the channel surfaces and the y
-axis is perpendicular to it. Lower wall of the channel is located at y
=
-
H and upper wall is at y
=
H. The fluid is injected into the channel and extracted out at a uniform velocityV (V
>
0 suction and V
<
0 injection) from upper wall and lower wall, respectively. A uniform magnetic field of strength B
° is applied perpendicular to the velocity field. The induced magnetic field is negligible as compared with the imposed field.The constitutive equation for Casson fluid can be written as [17] (1)
τ
i
j
=
2
μ
B
+
τ
y
2
π
e
i
j
,
π
>
π
c
,
2
μ
B
+
τ
y
2
π
c
e
i
j
,
π
<
π
c
,where μ
B is the plastic dynamic viscosity of the non-Newtonian fluid, τ
y is the yield stress of the fluid, π is the product of the component of deformation rate with itself, and π
c is critical value of π based on non-Newtonian model. Under these assumptions the governing equations for MHD boundary layer flow of Casson fluid are expressed as the following equations:(2)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(3)
u
∂
u
∂
x
+
v
∂
u
∂
y
=
-
1
ρ
∂
p
∂
x
+
υ
1
+
1
β
∂
2
u
∂
y
2
-
σ
B
°
2
u
ρ
±
g
β
T
T
-
T
2
+
β
C
C
-
C
2
,
(4)
u
∂
v
∂
x
+
v
∂
v
∂
y
=
-
1
ρ
∂
p
∂
y
+
υ
1
+
1
β
∂
2
v
∂
x
2
,
(5)
u
∂
T
∂
x
+
v
∂
T
∂
y
=
k
ρ
C
p
∂
2
T
∂
y
2
,
(6)
u
∂
C
∂
x
+
v
∂
C
∂
y
=
D
∂
2
C
∂
y
2
-
κ
1
C
,where ρ is density, μ is dynamic viscosity, ν is kinematic viscosity, σ is electrical conductivity, β is Casson fluid parameter,T is temperature of the fluid,k is thermal conductivity, κ
1 is reaction rate,D is mass diffusion, andC is the concentration field.Along with boundary conditions,(7)
u
=
0
,
v
=
V
2
,
T
=
T
2
,
C
=
C
2
at
y
=
H
,
∂
u
∂
y
=
0
,
v
=
0
,
T
=
T
1
,
C
=
C
1
at
y
=
0
.Introducing stream function such that u
-
=
∂
ψ
/
∂
y, v
-
=
-
∂
ψ
/
∂
x and eliminating pressure term from (3) and (4) by introducing vorticity ω, we get(8)
u
∂
ω
∂
x
+
v
∂
ω
∂
y
=
υ
1
+
1
β
∂
2
ω
∂
x
2
+
∂
2
ω
∂
y
2
-
σ
B
°
2
u
′
ρ
±
∂
∂
y
g
β
T
T
-
T
2
+
β
C
C
-
C
2
,where ω
=
∂
v
-
/
∂
x
-
∂
u
-
/
∂
y.Define(9)
x
∗
=
x
H
,
y
∗
=
y
H
,
u
=
-
V
x
∗
f
′
y
∗
,
v
=
V
f
y
∗
,
θ
y
∗
=
T
-
T
2
T
1
-
T
2
,
ϕ
y
∗
=
C
-
C
2
C
1
-
C
2
.Then governing nonlinear momentum and energy equations of the proposed problem can be written as(10)
1
+
1
β
f
′
′
′
′
-
M
2
f
′
′
+
R
f
′
f
′
′
-
f
f
′
′
′
+
λ
θ
′
+
N
ϕ
′
=
0
,
θ
′
′
-
Pr
⁡
f
θ
′
=
0
,
φ
′
′
-
S
c
γ
φ
-
S
c
f
φ
′
=
0
.Boundary conditions (7) reduce into(11)
f
1
=
1
2
,
f
′
1
=
0
,
θ
1
=
0
,
φ
1
=
0
,
f
′
′
0
=
0
,
f
0
=
0
,
θ
0
=
1
,
φ
0
=
1
.Here, R
=
V
H
/
ν is Reynolds number (R
>
0 for suction R
<
0 and for injection), M
2
=
σ
B
°
2
H
2
/
μ is Hartman number, λ
=
G
r
x
/
R
2 is the thermal buoyancy parameter, G
r
x
=
V
H
4
g
β
T
T
1
-
T
2
/
x
ν
3 is Grashof number, N
=
β
C
C
1
-
C
2
/
(
T
1
-
T
2
) is concentration buoyancy parameter, P
r
=
ρ
C
P
H
V
/
κ is Prantl number, S
c
=
H
V
/
D is Smith number, and γ
=
κ
1
H
/
ν is chemical reaction rate.


## 3. Numerical Solution of the Problem via Shooting Method

We employ shooting method in order to find the numerical solution of the proposed problem. For this we have to reduce the governing boundary value problem (10) into initial value problem by assuming x
1
=
η
,
x
2
=
f
,
x
3
=
f
′
,
x
4
=
f
′
′
,
x
5
=
f
′
′
′
,
x
6
=
θ
′
,
x
7
=
φ
′; then the following system is obtained(12)
x
1
′
x
2
′
x
3
′
x
4
′
x
5
′
x
6
′
x
7
′
=
1
x
3
x
4
x
5
β
1
+
β
M
2
x
4
-
R
x
3
x
4
-
x
2
x
5
-
λ
x
6
+
N
x
7
Pr
⁡
x
2
x
6
S
c
γ
φ
+
S
c
x
2
x
7
,with initial condition(13)
x
1
1
x
2
1
x
3
1
x
4
1
x
5
1
x
6
1
x
7
1
=
1
1
2
0
α
1
α
2
α
3
α
4
.Here, α
1, α
2, α
3, α
4 are missing initial conditions. In a shooting method, the missing (unspecified) initial condition at the initial point of the interval is assumed, and the differential equation is then integrated numerically as an initial valued problem to the terminal point. The accuracy of the assumed missing initial condition is then checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. If a difference exists, another value of the missing initial condition must be assumed and the process is repeated. This process is continued until the agreement between the calculated and the given condition at the terminal point is within the specified degree of accuracy. It is important to notice that we have to shoot the values of α
1, α
2, α
3, α
4. Since these values are not given in the boundary condition (13), by trial and error, suitable guess values are made and integration is carried out. The details of shooting method with Maple implementationshoot have been described by Meade et al. [28].


## 4. Results and Discussions

Present section is devoted to discuss the numerical results of our finding both in tabulation and in graphical form. Our main motive is to investigate the multiple solutions of the proposed problem. We have prepared Figures1–6 in order to evaluate the effects of different physical quantities on the velocity, temperature, and mass fraction. The effects include the existence of multiple solutions, Reynold number R, magnetic field M, Casson parameter β, Prandtl number P
r, Smith number S
c, and chemical reaction rate γ.Figure 1
Skin friction-
f
′
′
1 against the values of Reynold number R.Figure 2
Effect of Reynold numberR on velocity profile f
′
η.
(a)
(b)
(c)Figure 3
Effect of magnetic fieldM on velocity profile f
′
η.
(a)
(b)
(c)Figure 4
Effect of Casson parameterβ on velocity profile f
′
η.
(a)
(b)
(c)Figure 5
Effect of Prandtl numberP
r on temperature profile θ
η.
(a)
(b)
(c)Figure 6
Effect of Smith numberS
c on concentration profile ϕ
η for λ
=
0
.
(a)
(b)
(c)Based on the finding of multiple solutions it is concluded that there is only one solution in the case ofR
<
0, β
∈
0
,
∞ and β
∈
0,5, R
∈
0
,
∞. However, there exist multiple solutions for β
∈
5
,
∞ and R
∈
31.07
,
∞ for any value of magnetic number M
∈
0,2.0. So we can say that there is a critical value of Casson number β and suction parameter R such that β
c
r
i
t
i
c
a
l
=
5 and R
c
r
i
t
i
c
a
l
=
31.07. So it can be defended as there is no multiple solutions if β
<
β
c
r
i
t
i
c
a
l and <
R
c
r
i
t
i
c
a
l. The said phenomena can be observed from Figure 1. We plot the magnitude of the skin friction f
′
′
1 against the values of Reynold number R.Effects of Reynold numberR on velocity profile f
′
η for non-buoyant flow case λ
=
0 are presented in Figure 2. It is noticed that with enhanced values of Reynolds number R, velocity profile f
′
η decreases near the center of the channel for the first (1st) and third (3rd) solutions and increases near the walls of the channel. However, totally reverse behavior is observed for second (2nd) solution. Furthermore, triple solutions exist only for R
≥
31.07. Effect of magnetic field M on velocity profile f
′
η corresponding to forced convection or non-buoyant flow case λ
=
0 is shown in Figure 3. It is concluded that velocity profile f
′
η for 1st and 2nd solutions decreases near the center of the channel η
≈
0 and increases near the channel walls. Since magnetic field is applied perpendicular to the channel walls, effect of magnetic field can be seen clearly near the channel walls η
≈
1. Physically we can say that magnetic field enhances the viscosity of the fluid due to the chain deformation of the fluid particles. The chain-like structure retards the flow and decelerates the motion. This results in the fact that the fluid flow can be controlled by applying magnetic field which results in many control based applications including MHD power generation, casting of metals, blood flow in arteries, and many, many more. Effect of Casson parameter β on velocity profile f
′
η for the case of forced convection λ
=
0 is depicted in Figure 4. Velocity profile f
′
η decreases gradually for the 1st and 3rd solutions near the center of the channel and increases near the channel wall η
≈
1. However, velocity increases for the 2nd solution as the enhancement of the Casson parameter β is taking place in the neighborhood of η
≈
0. Furthermore, β
→
∞ corresponds to the Newtonian fluid. This is because of the fact that fluid is extracted from the walls of the channel with constant velocity that decreases the viscosity of the fluid particles which results in increase in the velocity near the wall of the channel. Figure 5 presented the effect of Prandtl number P
r on temperature profile θ
1 for non-buoyant λ
=
0 case. It is concluded from these profiles that temperature profile θ
1 increases strictly monotonically as the Prandtl number P
r increases for 1st and 2nd solutions. On the other hand, temperature profile θ
1 decreases as the Prandtl number P
r increases for 3rd solution. Effects of Smith number S
c on concentration profile ϕ
η for λ
=
0 (non-buoyant case) are plotted in Figure 6. Concentration profile ϕ
η decreases by the increase of the strength of Smith number S
c for all triple solutions. Physically it can be argued that an increase in Smith number S
c prompts the increment in the quantity of solute atoms experiencing substance reaction coming about the way that concentration of the fluid’s molecule diminishes monotonically. In Figure 7, we plot θ
′
1 against the values of Reynold number R for the case of forced convection λ
=
0. It is seen that 1st and 2nd solutions overlap each other as the Reynold number R increases.Figure 7
Effect of Reynold numberR on θ
′
1.In Tables1 and 2, we represented numerical values of skin friction f
′
′
1 and θ
′
1 for the variations of different physical parameters. Table 1 presented the numerical values of skin friction f
′
′
1 for the variation of buoyancy parameter λ by fixing R
=
36, M
=
N
=
0.5, β
=
5, P
r
=
1, S
c
=
1, γ
=
1.2. Magnitude of the skin friction increases for the 1st and 3rd solutions and decreases for 2nd solution in the case of opposing flow λ
<
0. Furthermore, same behavior is depicted for the case of assisting flow λ
>
0; therefore the fluid velocity near the channel walls η
≈
1 increases. From Table 2 it is concluded that magnitude of skin friction f
′
′
1 and θ
′
1 increases and decreases, respectively, as Reynold number increases for the case of assisting flow λ
>
0 only for 1st and 3rd solutions by setting M
=
N
=
0.5, P
r
=
1, S
c
=
2, γ
=
1.2. Numerical values of skin friction f
′
′
1 increase for 1st and 3rd solutions and decrease for 2nd solution by the variation of Reynolds number for the case of opposing flow λ
<
0 by setting M
=
N
=
0.5, P
r
=
1, S
c
=
2, γ
=
1.2.Table 1
Skin friction for different values of buoyancy parameters forR
=
36, M
=
N
=
0.5, β
=
5, P
r
=
1, S
c
=
1, γ
=
1.2.
λ
1st solution
2nd solution
3rd solution
f
′
′
1
f
′
′
1
f
′
′
1
−0.50
−3.29405873
−14.78794525
−27.20927012
−0.25
−4.19831393
−14.22940191
−27.20983823
0
−5.21362887
−13.55733409
−27.21039669
0.25
−6.42985179
−12.68188468
−27.21094599
0.50
−8.22263120
−11.22740570
−27.21148660Table 2
Skin friction and temperature gradient for different values of Reynold numberR.
R
λ
1st solution
2nd solution
3rd solution
1st solution
2nd solution
3rd solution
f
′
′
1
f
′
′
1
f
′
′
1
θ
′
1
θ
′
1
θ
′
1
31.07
−0.25
−3.56529871
−16.7538394
−23.89501535
−1.20630237
−1.12415316
−1.04645775
35
−4.06923819
−14.50163765
−26.71544153
−1.20492033
−1.14909835
−1.01703703
40
−4.70569981
−13.72648291
−28.87613992
−1.20354537
−1.16169799
−0.99886321
31.07
0
−4.22815985
−16.33685245
−23.92995344
−1.20259758
−1.12723885
−1.045737
35
−5.00282319
−13.89608196
−26.71780315
−1.20017552
−1.15257708
−1.01679981
40
−6.10401248
−12.73103850
−28.87353590
−1.19720323
−1.16644810
−0.99874016
31.07
0.25
−4.93400210
−15.87676698
−23.96319996
−1.19862360
−1.13057525
−1.04504186
35
−6.08426152
−13.14036724
−26.72012735
−1.19463415
−1.15685896
−1.01656453
40
−8.28517017
−10.94971752
−28.87095120
−1.18721976
−1.17484905
−0.99861750


## 5. Conclusion

Multiple solutions of MHD Casson fluid flow in a channel with heat and mass transfer are analyzed. Here are some important observations have been engendered in the light of numerical investigation:(i)
There is only single solution for the case of injectionR
<
0 for any value of Casson number β or Hartman number M.
(ii)
There exist multiple solutions only forβ
∈
5
,
∞ and R
∈
31.07
,
∞ for any value of magnetic number M
∈
0,2.0.
(iii)
Effect of Casson numberβ
≥
5 and Reynolds number R
>
31.06 on velocity profile f
′
η increases for 1st and 3rd solutions near the channel wall η
≈
1.
(iv)
Effect of Reynold numberR on skin friction f
′
′
1 and θ
′
1 is the same for assisting flow λ
>
0 and opposing flow λ
<
0.



---

*Source: 7535793-2016-09-04.xml*