MSC: 76U05, 76N20, 76A02
The inability of most uids to meet the in-
dustrial thermophysical requirements promp-
ted Maxwell [1] to consider suspending solid
particles, mostly millimetre-sized, in uids.
The development improved the thermal and
electrical characteristics of the uid but the
solid particles, being heavy, settle at the wall
and give rise to challenges such as clogging and
pipe erosion. The suggestion of Choi and East-
man [2] to replace the millimetre-sized particle
with a nanometre-sized particle revolutionised
the study of uid ow. The idea resolved the
clogging and pipe erosion problem of Maxwell's
approach [3, 4]. Since the innovation, authors
have considered the suspension of more than
one type of nanoparticle. Lee et al. [5] showed
that the thermal conductivity of the base uid
increased by more than 20% when two nano-
Coriolis Effect on the Flow of Water Carrying CNT, Graphene, and
Alumina Nanoparticles over a Heated Moveable Non-porous Surface
OLUWASEUN BIODUN ONUOHA1, FAWWAZ BATAYNEH2,
ABAYOMI SAMUEL OKE3,*,a, MARIO RASO4
1Department of Mathematical Sciences, Adekunle Ajasin University, NIGERIA
2Mathematics Division, College of Engineering, Kuwait College of Science and Technology, Doha, KUWAIT
3Department of Mathematical Sciences, Adekunle Ajasin University, NIGERIA
4Department of Computer Science, Sapienza University of Rome, ITALY
*Corresponding Author
aORCiD: 0000-0003-3903-4112
Abstract: Heat transfer fluids, heat exchangers, and coolants in electronics are typical industrial applications
where improved fluids are required for optimal performance. Stemming from the increasing demand, this study
examines the effects of suction, heat source and stretch- ing ratio on the flow and heat transfer in a
magnetohydrodynamic ternary hybrid nanofluid across a moveable rotating surface. Carbon nanotubes,
graphene and alumina are considered as the nanoparticles with water as the base fluid. The governing equations
are transformed from partial to ordinary differential equations. The equilibrium point of the system was ob-
tained, the conditions for the system stability were established and the emerging parameters were chosen within
the acceptable interval. The equations are numerically solved using the MATLAB bvp4c solver. The effects of
the flow parameters on the velocity and temperature distribution are graphically illustrated. The analysis shows
that the stretching ratio reduces flow temperature and velocity but increases skin friction. Coriolis force
enhances the heat transfer rate and increases the primary skin friction. Heat source increases flow temperature
and secondary skin friction.
Keywords: Ternary hybrid nanofluid; Coriolis force; Heat source; Suction; Stretching ratio
Received: February 25, 2024. Revised: September 7, 2024. Accepted: October 8, 2024. Published: November 21, 2024.
1. Background Information
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
284
Volume 3, 2024
particles were suspended in the base uid and
Baghbanzadeh et al. [6] showed that the heat
transfer rate increased by 14%. Suresh et al.
[7] considered the hybrid nanouid of copper-
alumina hybrid nanouid, Hayat and Nadeem
[8] studied silver-copper(II)oxide in water and
Oke et al. [9] considered the ow of copper-
alumina in water. Other recent studies on hy-
brid nanouid can be found in [10, 11, 12]. Re-
cently, it has been shown that ternary hybrid
nanouids enhance the thermal and electrical
properties of the uids better than the nano-
uids or hybrid nanouids [13]. Mousavi et al.
[14] examined the impact of volume fraction
and temperature on the thermophysical proper-
ties of a water-based ternary hybrid nanouid
of CuO, MgO, and TiO. The results showed
an improvement in the thermal conductivity
of the ternary hybrid nanouid with increasing
volume fraction. Furthermore, Elnaqeeb et al.
[13] and Oke [15] considered the suspension of
carbon nanotubes, graphenes and alumina in
water. The study considered the eect of vari-
ous particle shapes and sizes on the thermal
properties of the uid. The results show that
temperature is enhanced at a small volume frac-
tion.
Power generation, cancer tumour treatment
and magnetic devices for cell separation are
all practical applications where the use of
magnetohydrodynamic (MHD) ows play vi-
tal roles. MHD ow becomes more applic-
able when the geometry of ow involves ro-
tation as seen in ([16, 17, 18, 19, 15]). Oke
et al. [20] explored the eect of rotating sur-
faces on the MHD ow of water and air. It
was discovered that increasing surface rota-
tion raises ow temperature and that increas-
ing magnetic eld strength counteracts the ef-
fects of increased Coriolis force on skin friction
and heat transfer rate. Oyem et al. [21] con-
sidered the eects of variable thermal conduct-
ivity on an MHD Blasius and Sakiadis ows.
The results show that velocity proles reduce
and the temperature increases with increasing
magnetic eld strength for both Sakiadis and
Blasius ows. The impact of magnetic eld on
the transport of micropolar uid was reported
by Fatunmbi et al. [22]. Ouru et al. [23] showed
that the magnetic eld strength also inhibits
the ow velocity and increases temperature pro-
les. The ow of a reactive tangent hyperbolic
uid over a non-linear stretching surface was
evaluated by Fatunmbi et al. [24]. The authors
pointed out that the momentum boundary layer
and uid velocity were reduced due to the drag-
like Lorentz force occasioned by the magnetic
eld. Insightful results on the MHD ow of
Eyring-Powell uid over a stretching sheet with
a convective boundary condition were shown in
Oke and Mutuku [25]. Williamson uid ow
over an inclined rotating surface is studied in
the presence of a magnetic eld by Juma et al.
[26].
Carbon is the fourth most abundant chemical
element in the universe, following the gases hy-
drogen, helium, and oxygen [27]. As the most
abundant solid element, carbon is an excel-
lent candidate for use in nanoparticles within
nanouids. Carbon exists in several allotropes
[28], including graphite, diamond, fullerenes,
graphene, and carbon nanotubes. These al-
lotropes have diverse industrial applications,
such as steel production, where graphite is used
for electrodes and lubricants, and diamonds for
cutting tools and abrasives [29, 30]. Given their
exceptional properties, this study focuses on
the use of carbon nanotubes and graphene as
nanoparticles in the nanouids. From the lit-
erature above, studies have been carried out to
consider the ow of ternary hybrid nanouid
over moveable surfaces. However, the studies
have been restricted to nanoparticles of dif-
ferent shapes, and non-rotating surfaces and
without considering heat source. These are the
research gaps lled with this study where the
ow of ternary hybrid nanouid. (made from
the suspension of carbon nanotube (CNT), alu-
mina and graphene nanoparticles in water) over
a rotating moveable surface under the inuence
of heat source. This study provides answers to
the following research questions;
1. How can the ow parameters be chosen
to sustain model stability?
2. what is the inuence of the stretching ra-
tio on the transport phenomenon of the
ow of ternary hybrid nanouid over a
moveable rotating surface?
3. to what extent does channel rotation af-
fect the heat transfer rate, skin friction
and other transport properties of the ow
of ternary-hybrid nanouid over a move-
able rotating surface?
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
285
Volume 3, 2024
4. how does the heat source aect the tem-
perature, velocity, skin friction and heat
transfer rate in the MHD ow of ternary
hybrid nanouid over a rotating surface?
This study models the transport processes of
a water-based ternary hybrid nanouid that
suspends three distinct nanoparticles of sim-
ilar shapes in a moveable rotating channel. The
ow is taken to be laminar and incompressible,
subject to a magnetic eld of constant mag-
netic eld strength
B0.
The surface produces
the imaginary Coriolis force due to rotation
with an angular velocity of
.
Assuming the
ow adheres to the boundary layer theory [31]
and the Sakiadis [32] theory on a moveable sur-
face, the physical conguration of the ow is
depicted in Figure (1). The uid layers adhered
to the surface obeyed the no-slip requirement,
and the surface is extended linearly in the
xy
-
plane. The
xy
- surface is assumed to move lin-
early while the uid layers glued to the surface
obey the no-slip condition. Hence, following
the work of Elnaqeeb et al. [13], Oke et al. [33]
and Oke [15], the ow can be modelled as the
following equations;
xu+yv+zw= 0,
(1)
(u∂x+vy+wz)u
=µtf
ρtf
zzu2Ωuσtf B2
0
ρtf
u,
(2)
(u∂x+vy+wz)v
=µtf
ρtf
zzv+ 2Ωvσtf B2
0
ρtf
v,
(3)
(u∂x+vy+wz)T
=αtf zzT+Q0
(ρCp)tf
(TwT) exp zra
νbf .
(4)
At the base of the channel
(z= 0) ,
the uid
layers on the wall admit the same velocity as
that of the stretching wall and thus, the no-slip
condition is retained as
on the
x
-axis
:
u=ax, w =zw(bf )1
2Zv
ay .
on the
y
-axis
:
v=ay, w =zw(bf )1
2Zu
ax,
on the
(x, y)
-plane
, T =Tw.
Furthermore, at the free stream
(z ),
the
free stream conditions are sustained as
u0,
on the
x
-axis
,
v0,
on the
y
-axis
,
T=T
on the
(x, y)
-plane
.
The characteristics of the ow that of practical
implication are the skin drags and heat transfer
rate in all directions which are dened as
Cfx =µtf ρbf a2x21u
z
z=0
,
(5)
Cgy =µtf ρbf a2y21v
z
z=0
,
(6)
Nux=tf (κbf (TwT))1T
z
z=0
,
(7)
Nuy=yκtf (κbf (TwT))1T
z
z=0
.
(8)
Since the properties of the resulting ternary
nanouid properties are inuenced by the prop-
erties of both the base uid and the nan-
oparticles, the following denitions extracted
from literature [33, 13] are adopted in this study
µtf =1
ϕ
3
X
j=1
µjϕj, κtf =1
ϕ
3
X
j=1
κjϕj,
(9)
ρtf = (1 ϕ)ρbf +
3
X
j=1
ρjϕjρtf =ρbf A1,
(10)
(ρcp)tf = (1 ϕ) (ρcp)bf +
3
X
j=1
(ρcp)jϕj
(ρcp)tf = (ρcp)bf A2
where
A1= 1 ϕ+1
ρbf
3
X
j=1
ρjϕj,
A2= 1 ϕ+1
(ρcp)bf
3
X
j=1
(ρcp)jϕj.
The thermophysical properties of the nano-
particles and the base uid are provided in table
(3). Also, the ternary nanouid is made up of
three dierent nanouids. Considering platelet-
shaped nanoparticles, then the viscosity and
2. Research Methodology
2.1 Mathematical Formulation
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
286
Volume 3, 2024
thermal conductivity of the nanouids are mod-
elled as
µj=1 + 37.1ϕj+ 612.6ϕ2
jµbf
κj
κbf
=κj+ 4.7κbf 4.7ϕ(κbf κj)
κj+ 4.7κbf +ϕ(κbf κj)
for
j= 1,2,3
and
ϕj>0.02.
Hence, with this
model, we have
µtf =1
ϕ
3
X
j=1
µjϕj
=µbf
ϕ
3
X
j=1 1 + 37.1ϕj+ 612.6ϕ2
jϕj
=µbf
ϕA3
(11)
κtf =1
ϕ
3
X
j=1
κjϕj
=κbf
ϕ
3
X
j=1
κj+ 4.7κbf 4.7ϕ(κbf κj)
κj+ 4.7κbf +ϕ(κbf κj)ϕj
=κbf
ϕA4,
(12)
where
A3=
3
X
j=1 1 + 37.1ϕj+ 612.6ϕ2
jϕj,
A4=
3
X
j=1
κj+ 4.7κbf 4.7ϕ(κbf κj)
κj+ 4.7κbf +ϕ(κbf κj)ϕj.
Employing the similarity variables
u=axf, v =acyg, w =(bf )1
2(f+cg)
(13)
T=T+ (TwT)θ, η =za
νbf 1
2
,
(14)
the equations are rendered dimensionless and
the number of parameters is reduced to become
A3
ϕA1
f′′′ + (f+cg)f′′ f2K+M
A1f
= 0,
(15)
A3
ϕA1
g′′′ + (f+cg)g′′ cg2+KM
A1g
= 0,
(16)
A4
ϕA2
Θ′′ +P rQ
A2
eη+P r (f+cg) Θ
(17)
with the boundary and initial conditions
at
η= 0; f=fw, f= 1, g =fw
c, g= 1,Θ = 1,
(18)
as
η ;f0, g0,Θ0.
(19)
A1= 1 ϕ+1
ρbf
3
X
j=1
ρjϕj,
A2= 1 ϕ+1
(ρcp)bf
3
X
j=1
(ρcp)jϕj,
A3=
3
X
j=1 1 + 37.1ϕj+ 612.6ϕ2
jϕj,
A4=
3
X
j=1
κj+ 4.7κbf 4.7ϕ(κbf κj)
κj+ 4.7κbf +ϕ(κbf κj)ϕj.
K=2Ω
a, M =σB2
0
bf
, Q =Q0
a(ρcp)bf
,
P r =νbf
αbf
, αbf =κbf
(ρcp)bf
.
The ow characteristics that are of practical in-
terest are also reduced to
Re
1
2
xCfx =ϕ1A1f′′ (0) , Re
1
2
yCgy =ϕ1A1g′′ (0)
(20)
Re1
2
xNux=Re1
2
yNuy=ϕ1A3Θ(0) .
(21)
The system of equations (15) - (17) has bound-
ary conditions (18) and (19) which makes it
dicult to adopt a direct numerical solution.
2.2 Similarity Transformation
2.3 Method of Solution
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
287
Volume 3, 2024
Hence, the shooting technique is required to
convert the boundary conditions to initial con-
ditions. To start with, equations (15) - (17) are
rewritten as a system of rst-order dierential
equations by setting
h1=f, h2=f, h3=f′′, h4=g, h5=g,
h6=g′′, h7= Θ, h8= Θ
to have
h
1=h2,
(22)
h
2=h3,
(23)
h
3=ϕA1
A3(h1+ch4)h3h2
2K+M
A1h2,
(24)
h
4=h5,
(25)
h
5=h6,
(26)
h
6=ϕA1
A3(h1+ch4)h6ch2
5+KM
A1h5,
(27)
h
7=h8,
(28)
h
8=ϕA2
A4
P r (h1+ch4)h8+Q
A2
exp (η).
(29)
However, we have a combination of the initial
and boundary conditions (18) and (19). The
shooting technique requires that all these condi-
tions be rewritten as initial conditions, by pla-
cing an arbitrary initial guess for the variables
whose conditions are not known. In this case,
we have the conditions as
h1=fw, h2= 1, h3=ε1, h4=fw
c,
(30)
h5= 1, h6=ε2, h7= 1, h8=ε3.
(31)
where
ε1, ε2, ε3
are arbitrary guesses that will
be updated after every iteration until the
boundary conditions
h2()=0, h5()=0, h7()=0
are satised up to a specied tolerance level.
The MATLAB bvp4c solver is used to solve
the system (See Oke [34] for some other semi-
analytical methods of solution). The bvp4c
solver [35] executes the nite dierence code
that implements the three-stage Lobatto IIIa
formula of the fourth order. The system is
solved with an absolute tolerance of
104
and
a relative tolerance of
104.
The stability analysis of the equations (22 - 29)
is carried out in this section. This becomes ne-
cessary to identify the range of values accept-
able for the parameter values. By setting right
hand side of the equations to
0,
we have the
equilibrium point as
(h1, h2, h3, h4, h5, h6, h7, h8)
=a1,0,0, a4,0,0, a7,Qexp (η)
A2(a1+ca4).
The characteristic equation obtained by substi-
tuting the equilibrium point into the Jacobian
of the system is
λ1 0 0 0 0 0 0
0λ1 0 0 0 0 0
0a32 a33 λ0 0 0 0 0
0 0 0 λ1000
0 0 0 0 λ100
0 0 0 0 0 a66 λ0 0
0 0 0 0 0 0 λ1
γ0 0 γc 000a88 λ
= 0,
where
γ=ϕA2
A4
P rh8, a32 =ϕA1
A3K+M
A1,
a33 =ϕA1
A3
(a1+ca4), a66 =ϕA1
A3
(a1+ca4),
a88 =ϕA2
A4
P r (a1+ca4).
The characteristic equation becomes
λ4(a88 λ) (a66 λ)λ2a33λa32= 0, λ1,2,3,4= 0,
λ5=a88 <0, λ6=a66 <0
provided
a1+ca4>0,
λ7,8<0,
provided
a33 <0, a32 <0.
Hence, the system is stable provided
a1+ca4>
0.
The values of the parameters will therefore
be chosen so that
a1+ca4>0
and since
ϕA1
A3K+M
A1>0,
then the values of
M
and
K
can be chosen arbitrarily.
The results obtained in this study are validated
with the results obtained using the MATLAB
bvp5c solver when
M0, K 0.
The res-
ults are shown in Table 4. There is an excel-
lent agreement between the results. Hence, the
bvp4c is used to numerically study the problem
under consideration.
2.4 Qualitative analysis
2.5 Validation of Results
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
288
Volume 3, 2024
In the ow under consideration, the ternary hy-
brid nanouid is made from the suspension of
carbon nanotubes, graphene nanoparticles and
Al2O3
nanoparticles in water base uid. By
rendering the governing equations dimension-
less, the ow becomes equipped with some di-
mensionless parameters. This section discusses
the dynamics of the ow as the dimensionless
parameters (stretching parameter
c,
suction ve-
locity
fw,
Coriolis force
K,
magnetic strength
M,
and heat source
Q
) vary during the ow.
The ow properties are visualised graphically
to illustrate the behaviour of the ow proper-
ties (the primary velocity, secondary velocity
and temperature) as the ow parameters are
varied. Figures (2a) and (2b) visually illustrate
the behaviour of the velocity components in the
x
- and
y
-directions. It is observed that the ve-
locity decreases in all directions as the stretch-
ing ratio increases. The stretching ratio is var-
ied over the interval
[0.1,1.0]
where
c= 0.1
indicates that stretching in the
x
-direction is
faster than stretching in the
y
-direction and
c= 1.0
indicates that stretching is happen-
ing at the same rate in both directions. Sim-
ilarly, temperature is shown in Figure (2c) to
reduce as the stretching ratio increases. The
inuence of suction velocity
fw
(chosen from
the interval
[0.1,0.7]
) on the ow velocity and
temperature is shown in Figures (3a), (3b) and
(3c). As observed in Figures (3a) and (3b), in-
creasing suction reduces velocity in all direc-
tions and Figure (3c) shows that temperature
reduces with increasing suction. The eects of
Coriolis force, measured by the dimensionless
parameter
K,
are studied in Figures (4a) and
(4b). By varying the Coriolis force over the
interval
K[0.01,1.00] ,
where
K= 0.01
in-
dicates the presence of nearly no rotation while
K= 1.00
represents a moderately rotating sur-
face. Coriolis force reduces the primary velo-
city decreases (as shown in Figure (4a)) but
increases the secondary velocity (as shown in
Figure (4b)). The ow is equipped with some
heat source and it is observed in Figure (5) that
the temperature of the ow increases as heat
source increases.
The skin friction in both
x
- and
y
-directions are
displayed in Figures (6a) and (6b) for varying
Coriolis force, stretching ratio and heat source.
Figure (6a) shows that skin friction in the
x
-
direction increases as either Coriolis force or
stretching ratio increases. More so, increas-
ing heat source only slightly decreases the skin
friction in the
x
-direction. Figure (6b) shows
that skin friction in the
y
-direction increases as
stretching ratio increases but decreases with in-
creasing Coriolis force. More so, increasing the
heat source only slightly increases the skin fric-
tion in the
y
-direction. Figure (7) shows the
behaviour of heat transfer rate. Increasing the
Coriolis force only contributes very slowly to an
increase in the heat transfer rate, a higher heat
source leads to a reduction in the heat transfer
rate and the stretching ratio decreases the heat
transfer rate.
A summary of the results is provided in Table
(1). The results are discussed in this section
and our ndings are compared with existing lit-
erature.
The stretching ratio described by Makhdoum
et al. [36] is the ratio of wall stretching in
one direction to wall stretching in another dir-
ection. In this case, the boundary condition
u(η= 0) = axf(η)
and
v(η= 0) = acyg(η)
gives a stretching ratio
c.
Establishing that
the primary direction is the
x
-direction, then
the ratio stretching ratio describes the ratio of
stretching in the
x
-direction to the stretching
in the
y
-direction. The values of
c
were restric-
ted to the interval
[0.1,1.0]
so that the stretch-
ing in the
y
-direction is less than or equal to
the stretching in the
x
-direction (where equal-
ity only holds when
c= 1
). According to the
study by Oke [15], the stretching ratio decreases
the primary velocity but increases the second-
ary velocity. However, our results show that
both the primary and the secondary velocities
are reduced by increasing the stretching ratio.
The conict in the eect of the stretching ratio
on the secondary velocity could be due to re-
strictions placed on the secondary wall so that
the secondary wall cannot move faster than the
primary wall. The stretching ratio increases the
skin friction in all directions but decreases the
heat transfer rate. This outcome agrees with
the study of Makhdoum et al. [36]. Suction
3. Analysis and Discussion of Results
3.1 Analysis of Results
3.2 Discussion of Results
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
289
Volume 3, 2024
measures the amount of partial vacuum cre-
ated when the uid is drawn towards a region
of lower pressure [37]. Makhdoum et al. [37]
showed that suction decreases the temperat-
ure and this is in agreement with our ndings.
The partial pressure created by suction leads
to a reduction in the thermal energy distribu-
tion in the ow and hence, increasing the suc-
tion leads to a reduction in the temperature.
Coriolis force as described by [12] is the iner-
tia force generated in a rotating system and is
responsible for the ctitious deection of the
path of moving objects in the rotating frame.
Coriolis force increases as angular velocity in-
creases. As angular velocity increases, the ow
experiences sideways lurching to the secondary
direction. The dynamic equilibrium is main-
tained by the reduction in ow in the primary
direction [12]. The heat source, measured by
the heat source parameter
Q,
represents the ex-
ternal heat source introduced to the ow. In-
creasing the heat source enhances the inherent
thermal energy in the ow and thereby raises
the temperature. This nding agrees with that
found in the ref [4].
u v T Cfx Cf y N u
c
fw
- - -
K
-
Q
- -
Table 1: Summary of analysis
This work analyses the ow of ternary hy-
brid nanouid obtained by the suspension of
graphene and alumina nanoparticles in water.
The ow is subjected to some heat source
while the wall of the channel is moveable and
rotating. The ow equations are developed,
rendered dimensionless and simulated to illus-
trate the eects of various parameters. The
main outcomes of the study include;
The stretching ratio reduces ow temper-
ature, velocity in all directions and heat
transfer rate but increases the skin fric-
tion in all directions.
Heat transfer rate is enhanced with in-
creasing Coriolis force.
Increasing rotation leads to an increase in
both the skin friction in the primary ow
direction and the heat transfer rate but a
reduction in the skin friction in the sec-
ondary ow direction.
More heat source increases ow temper-
ature and secondary skin friction but de-
creases primary skin friction and heat
transfer rate.
This study focuses on ows subject to a low to
moderately high rotation, excluding extremely
high rotation. In instances of high rotation, the
uid dynamics becomes turbulent. Examining
the behaviour of hybrid ternary uid ow over
a rapidly rotating surface could be of scientic
interest, oering insight into a domain char-
acterized by turbulent uid dynamics. This
study has emphasised on using similar nano-
particle shapes, and some other works of lit-
erature (see [38, 39, 13, 33]) have emphasised
dissimilar nanoparticles. It will be of benet
to the industry if comparisons between ows of
ternary hybrid nanouid with similar and dis-
similar nanoparticle shapes can be made.
[1] James Clerk Maxwell. A Treat-
ise on Electricity and Magnetism.
Nature
, 7:478480, 1873. doi: ht-
tps://doi.org/10.1038/007478a0.
[2] Stephen U. S. Choi and J. A. Eastman.
Enhancing Thermal Conductivity of Flu-
ids with Nanoparticles.
ASME Interna-
tional Mechanical Engineering Congress &
Exposition
, 1995.
[3] John Kinyanjui Kigio, Mutuku Wini-
fred Nduku, and Oke Abayomi Samuel.
Analysis of Volume Fraction and Con-
vective Heat Transfer on MHD Cas-
son Nanouid over a Vertical Plate.
Fluid Mechanics
, 7(1):18, 2021. doi:
10.11648/j.fm.20210701.11.
[4] A. S. Oke, I. L. Animasaun, W. N. Mu-
tuku, M. Kimathi, Nehad Ali Shah, and
S. Saleem. Signicance of Coriolis force,
volume fraction, and heat source/sink on
the dynamics of water conveying 47Â nm
4. Conclusion
References
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
290
Volume 3, 2024
alumina nanoparticles over a uniform sur-
face.
Chinese Journal of Physics
, 71:716
727, 2021. ISSN 0577-9073. doi: ht-
tps://doi.org/10.1016/j.cjph.2021.02.005.
[5] S. Lee, S. U.-S. Choi, S. Li, and J. A. East-
man. Measuring Thermal Conductivity
of Fluids Containing Oxide Nanoparticles.
Journal of Heat Transfer
, 121(2):280289,
may 1999. doi: 10.1115/1.2825978.
[6] Mohammadali Baghbanzadeh, Alimorad
Rashidi, Davood Rashtchian, Roghayeh
Lot, and Azadeh Amrollahi. Syn-
thesis of spherical silica/multiwall car-
bon nanotubes hybrid nanostructures
and investigation of thermal conductiv-
ity of related nanouids.
Thermochim-
ica Acta
, 549:8794, dec 2012. doi:
10.1016/j.tca.2012.09.006.
[7] S. Suresh, K. P. Venkitaraj, P. Selvak-
umar, and M. Chandrasekar. Eect
of Al2O3-Cu/water hybrid nanouid in
heat transfer.
Experimental Thermal and
Fluid Science
, 38:5460, apr 2012. doi:
10.1016/j.expthermusci.2011.11.007.
[8] Tanzila Hayat and S. Nadeem. Heat trans-
fer enhancement with Ag-CuO/water hy-
brid nanouid.
Results in Physics
, 7:
23172324, 2017. ISSN 2211-3797. doi:
10.1016/j.rinp.2017.06.034.
[9] A. S. Oke, B. C. Prasannakumara, W. N.
Mutuku, R. J. Punith Gowda, B. A. Ju-
maand R. Naveen Kumar, and O. I. Bada.
Exploration of the eects of Coriolis force
and thermal radiation on water-based hy-
brid nanouid ow over an exponentially
stretching plate.
Scientic Reports
, 12(1),
December 2022. ISSN 2045-2322. doi:
10.1038/s41598-022-21799-9.
[10] K. Thriveni and B. Mahanthesh. Sensit-
ivity analysis of nonlinear radiated heat
transport of hybrid nanoliquid in an annu-
lus subjected to the nonlinear Boussinesq
approximation.
Journal of Thermal Ana-
lysis and Calorimetry
, 143(3):27292748,
apr 2020. doi: 10.1007/s10973-020-09596-
w.
[11] K. Thriveni and B. Mahanthesh. Sens-
itivity computation of nonlinear con-
vective heat transfer in hybrid nanoma-
terial between two concentric cylinders
with irregular heat sources.
Interna-
tional Communications in Heat and Mass
Transfer
, 129:105677, dec 2021. doi:
10.1016/j.icheatmasstransfer.2021.105677.
[12] Abayomi Samuel Oke. Combined ef-
fects of Coriolis force and nanoparticle
properties on the dynamics of gold-
water nanouid across nonuniform sur-
face.
ZAMM - Journal of Applied
Mathematics and Mechanics / Zeitschrift
×r Angewandte Mathematik und Mech-
anik
, 0(0):e202100113, 2022. doi: ht-
tps://doi.org/10.1002/zamm.202100113.
[13] Thanaa Elnaqeeb, Isaac Lare Animasaun,
and Nehad Ali Shah. Ternary-hybrid
nanouids: signicance of suction and
dual-stretching on three-dimensional ow
of water conveying nanoparticles with vari-
ous shapes and densities.
Zeitschrift für
Naturforschung A
, 76(3):231243, 2021.
doi: doi:10.1515/zna-2020-0317. URL
https://doi.org/10.1515/zna-2020-0317
.
[14] S. M. Mousavi, F. Esmaeilzadeh, and X. P.
Wang. Eects of temperature and particles
volume concentration on the thermophys-
ical properties and the rheological beha-
vior of CuO/MgO/TiO2 aqueous ternary
hybrid nanouid.
Journal of Thermal Ana-
lysis and Calorimetry
, 137(3):879901, jan
2019. doi: 10.1007/s10973-019-08006-0.
[15] Abayomi Samuel Oke. Heat and Mass
Transfer in 3D MHD Flow of EG-Based
Ternary Hybrid Nanouid Over a Rotat-
ing Surface.
Arabian Journal for Science
and Engineering
, 2022. ISSN 2191-4281.
doi: 10.1007/s13369-022-06838-x. URL
https://doi.org/10.1007/s13369-022-06838-x
.
[16] Nadeem Abbas, S. Nadeem, Anber Saleem,
M. Y. Malik, Alibek Issakhov, and
Fahd M. Alharbi. Models base study
of inclined MHD of hybrid nanouid
ow over nonlinear stretching cylinder.
Chinese Journal of Physics
, 69:109117,
2021. ISSN 0577-9073. doi: ht-
tps://doi.org/10.1016/j.cjph.2020.11.019.
[17] Nur Syazana Anuar, Norfah Bachok,
and Ioan Pop. Cu-Al2O3/Water Hybrid
Nanouid Stagnation Point Flow Past
MHD Stretching/Shrinking Sheet in
Presence of Homogeneous-Heterogeneous
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
291
Volume 3, 2024
and Convective Boundary Conditions.
Mathematics
, 8(8), 2020. ISSN 2227-
7390. doi: 10.3390/math8081237. URL
https://www.mdpi.com/2227-7390/8/8/1237
.
[18] Nirmalendu Biswas, U. K. Sarkar, Ali J.
Chamkha, and Nirmal Kumar Manna.
Magneto-hydrodynamic thermal convec-
tion of Cu-Al2O3/water hybrid nano-
uid saturated with porous media subjec-
ted to half-sinusoidal nonuniform heating.
Journal of Thermal Analysis and Calor-
imetry
, 143(2):17271753, January 2021.
ISSN 1588-2926. doi: 10.1007/s10973-020-
10123-0.
[19] Iskander Tlili, Hossam A. Nabwey, M. Gir-
inath Reddy, N. Sandeep, and Maddileti
Pasupula. Eect of resistive heating on in-
cessantly poignant thin needle in magneto-
hydrodynamic Sakiadis hybrid nanouid.
Ain Shams Engineering Journal
, 12(1):
10251032, 2021. ISSN 2090-4479. doi: ht-
tps://doi.org/10.1016/j.asej.2020.09.009.
[20] Abayomi S. Oke, Winifred N. Mu-
tuku, Mark Kimathi, and Isaac Lare
Animasaun. Coriolis eects on MHD
newtonian ow over a rotating non-
uniform surface.
Proceedings of the
Institution of Mechanical Engineers,
Part C: Journal of Mechanical Engin-
eering Science
, 235(19):38753887, 2021.
doi: 10.1177/0954406220969730. URL
https://doi.org/10.1177/0954406220969730
.
[21] Anselm O. Oyem, Winifred N. Mutuku,
and Abayomi S. Oke. Variability eects
on magnetohydrodynamic for Blasius and
Sakiadis ows in the presence of Dufour
and Soret about a at plate.
Engineer-
ing Reports
, 2(10):e12249, aug 2020. doi:
10.1002/eng2.12249.
[22] Ephesus Olusoji Fatunmbi,
Hammed Abiodun Ogunseye, and Pre-
cious Sibanda. Magnetohydrodynamic
micropolar uid ow in a porous medium
with multiple slip conditions.
Interna-
tional Communications in Heat and Mass
Transfer
, 115:104577, jun 2020. doi:
10.1016/j.icheatmasstransfer.2020.104577.
[23] J. O. Ouru, W. N. Mutuku, and A. S. Oke.
Buoyancy-Induced MHD Stagnation Point
Flow of Williamson Fluid with Thermal
Radiation.
Journal of Engineering Re-
search and Reports
, pages 918, mar 2020.
doi: 10.9734/jerr/2020/v11i417065.
[24] Ephesus Olusoji Fatunmbi, Fazle Mabood,
Hedi Elmonser, and Iskander Tlili. Mag-
netohydrodynamic nonlinear mixed con-
vection ow of reactive tangent hyperbolic
nano uid passing a nonlinear stretchable
surface.
Physica Scripta
, 96(1):015204, nov
2020. doi: 10.1088/1402-4896/abc3e9.
[25] Abayomi S. Oke and Winifred N. Mu-
tuku. Signicance of viscous dissipation on
MHD EyringPowell ow past a convect-
ively heated stretching sheet.
Pramana
,
95(4), nov 2021. doi: 10.1007/s12043-021-
02237-3.
[26] Belindar A. Juma, Abayomi S. Oke, Wini-
fred N. Mutuku, Afolabi G. Ariwayo, and
Olum J. Ouru. Dynamics of William-
son uid over an inclined surface subject
to Coriolis and Lorentz forces.
Engineer-
ing and Applied Science Letter
, 5(1):37
46, March 2022. ISSN 2617-9709. doi:
10.30538/psrp-easl2022.0083.
[27] Chiaki Kobayashi, Amanda I. Karakas,
and Maria Lugaro. The Origin of Ele-
ments from Carbon to Uranium.
The As-
trophysical Journal
, 900(2):179, September
2020. ISSN 1538-4357. doi: 10.3847/1538-
4357/abae65.
[28] Xiaohui Ye, Ming Qi, Mengzhen Chen,
Lihui Zhang, and Jinying Zhang. Zero
to three dimension structure evolution
from carbon allotropes to phosphorus allo-
tropes.
Advanced Materials Interfaces
, 10
(5), December 2022. ISSN 2196-7350. doi:
10.1002/admi.202201941.
[29] Jiaqiang Li and Yu Han. Articial carbon
allotrope
γ
-graphyne: Synthesis, proper-
ties, and applications.
Giant
, 13:100140,
March 2023. ISSN 2666-5425. doi:
10.1016/j.giant.2023.100140.
[30] Leandro Nicolýs Sacco and Sten Vol-
lebregt. Overview of engineering car-
bon nanomaterials such as carbon nan-
otubes (cnts), carbon nanobers (cnfs),
graphene and nanodiamonds and other
carbon allotropes inside porous anodic alu-
mina (paa) templates.
Nanomaterials
, 13
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
292
Volume 3, 2024
(2):260, January 2023. ISSN 2079-4991.
doi: 10.3390/nano13020260.
[31] L. Prandtl. "Über Flüssigkeitsbewegung
bei sehr kleiner Reibung" translated to
"Motion of uids with very little viscos-
ity".
Int. Math.-Kongr. Heidelb.
, 8(13):1
8, 1904.
[32] B. C. Sakiadis. Boundary layer behavior
on continuous solidsurfaces: The boundary
layer on a continuous at surface.
Am.Inst.
Chem. Eng. (AIChE)
, 7:221225, 1961.
[33] Abayomi S. Oke, Ephesus O. Fatun-
mbi, Isaac L. Animasaun, and Belin-
dar A. Juma. Exploration of ternary-
hybrid nanouid experiencing Coriolis
and Lorentz forces: case of three-
dimensional ow of water conveying car-
bon nanotubes, graphene, and alumina
nanoparticles.
Waves in Random and
Complex Media
, pages 120, Septem-
ber 2022. ISSN 1745-5049. doi:
10.1080/17455030.2022.2123114.
[34] A. S. Oke. Convergence of dierential
transform method for ordinary dierential
equations.
Journal of Advances in Math-
ematics and Computer Science
, 24(6):1
17, 2017.
[35] L. F. Shampine, M. W. Reichelt, and
J. Kierzenka. Solving Boundary Value
Problems for Ordinary Dierential Equa-
tions in MATLAB with bvp4c. 2010.
[36] Basim M. Makhdoum, Zafar Mahmood,
Bandar M. Fadhl, Musaad S. Aldhabani,
Umar Khan, and Sayed M. Eldin. Sig-
nicance of entropy generation and nan-
oparticle aggregation on stagnation point
ow of nanouid over stretching sheet
with inclined Lorentz force.
Ara-
bian Journal of Chemistry
, 16(6):104787,
June 2023. ISSN 1878-5352. doi:
10.1016/j.arabjc.2023.104787.
[37] Basim M. Makhdoum, Zafar Mahmood,
Umar Khan, Bandar M. Fadhl, Ilyas Khan,
and Sayed M. Eldin. Impact of suc-
tion with nanoparticles aggregation and
Joule heating on unsteady MHD stag-
nation point ow of nanouids over ho-
rizontal cylinder.
Heliyon
, 9(4):e15012,
April 2023. ISSN 2405-8440. doi:
10.1016/j.heliyon.2023.e15012.
[38] I. L. Animasaun, A. S. Oke, Qasem M.
Al-Mdallal, and A. M. Zidan. Explor-
ation of water conveying carbon nan-
otubes, graphene, and copper nano-
particles on impermeable stagnant and
moveable walls experiencing variable tem-
perature: thermal analysis.
Journal of
Thermal Analysis and Calorimetry
, 148
(10):45134522, March 2023. ISSN 1588-
2926. doi: 10.1007/s10973-023-11997-6.
[39] Wenhao Cao, Animasaun I. L., Se-Jin
Yook, Oladipupo V. A., and Xianjun
Ji. Simulation of the dynamics of col-
loidal mixture of water with various nan-
oparticles at dierent levels of partial
slip: Ternary-hybrid nanouid.
Interna-
tional Communications in Heat and Mass
Transfer
, 135:106069, jun 2022. doi:
10.1016/j.icheatmasstransfer.2022.106069.
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
293
Volume 3, 2024
Appendix
Table 2:
Nomenclature
Dimensional quantities
x, y, z
Distance in three-dimensional space
(
L
)
σ
electrical conductivity
(
M1L3T3A2
)
u, v, w
Velocity component in the
x, y, z
-directions (
LT 1
)
α
thermal diusivity (
L2T1
)
T
Dimensional uid temperature (
K
)
B0
magnetic eld strength (
L1A
)
Tw
wall temperature (
K
)
κ
thermal conductivity (
K1
)
T
free stream temperature (
K
)
ρ
density (
ML3
)
cp
Specic heat capacity (
L2T2K1
)
angular velocity of the surface
(
T1
)
ν
kinematic viscosity (
L2T1
)
µ
dynamic viscosity (
ML1T1
)
Dimensionless quantities
f
primary velocity
K
rotation parameter
g
secondary velocity
P r
Prandtl number
Θ
temperature
fw
suction
M
Magnetic eld parameter
c
stretching ratio
Cfx
local skin friction in
x
-direction
ϕ
volume fraction
Q
Heat source parameter
Nux
Nusselt number in
x
-direction
Cfy
local skin friction in
y
-direction
Nuy
Nusselt number in
y
-direction
Subscripts
tf
ternary hybrid nanouid
bf
base uid
1,2,3
CNT, graphene,
Al2O3
nanoparticle
Z
Y
X
u=ax, w =zwbf Rv
ay
v=acy, w =zw(bf )1
2Ru
ax
¯
B¯
B¯
B
u
v
w
Figure 1: Flow conguration
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
294
Volume 3, 2024
Table 3: thermal and physical properties [20]
ρkgm3kW m1K1cpJkg1K1
carbon nanotubes 2100 3007.4 410
graphene 2200 5000 790
Al2O3
3970 40 765
water 997.1 0.613 4179
Table 4: Results Validation for
M=K= 0, fw= 0.3, P r = 6.2, η= 70
c Q f′′(0) g′′ (0) Θ(0)
bvp4c bvp5c bvp4c bvp5c bvp4c bvp5c
5.0 0.1 -0.6751 -0.6727 -0.9173 -0.9138 2.7848 -2.7685
5.5 0.1 -0.6943 -0.6918 -0.9568 -0.9531 2.8759 -2.858
6.0 0.1 -0.713 -0.7102 -0.9946 -0.9907 2.9636 -2.9441
0.3 0.5 -0.443 -0.4422 -0.3553 -0.3549 1.09 -1.0905
0.3 1.0 -0.443 -0.4422 -0.3553 -0.3549 0.4711 -0.4717
0.3 1.5 -0.443 -0.4422 -0.3553 -0.3549 -0.1478 0.1472
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
295
Volume 3, 2024
012345678910
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
f'( )
c = 0.1000
c = 0.3250
c = 0.5500
c = 0.7750
c = 1.0000
Horizontal velocity reduces
with increasing stretching
ratio c
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
(a) primary velocity with stretching ratio
0 2 4 6 8 10 12
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
g'( )
c = 0.1000
c = 0.3250
c = 0.5500
c = 0.7750
c = 1.0000
vertical velocity reduces
with increasing stretching
ratio c
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
(b) secondary velocity with stretching ratio
0123456
0
0.2
0.4
0.6
0.8
1
1.2
( )
c = 0.1000
c = 0.3250
c = 0.5500
c = 0.7750
c = 1.0000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
Temperature reduces
with increasing
stretching ratio c
(c) temperature with stretching ratio
Figure 2: stretching ratio eects
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
296
Volume 3, 2024
012345678910
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
f'( )
fw = 0.1000
fw = 0.2500
fw = 0.4000
fw = 0.5500
fw = 0.7000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
Horizontal velocity reduces
with increasing fw
(a) primary velocity with suction
012345678910
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
g'( )
fw = 0.1000
fw = 0.2500
fw = 0.4000
fw = 0.5500
fw = 0.7000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
vertical velocity reduces
with increasing fw
(b) secondary velocity with suction
0123456
0
0.2
0.4
0.6
0.8
1
1.2
( )
fw = 0.1000
fw = 0.2500
fw = 0.4000
fw = 0.5500
fw = 0.7000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
Temperature reduces
with increasing fw
(c) temperature with suction
Figure 3: suction eects
012345678910
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
f'( )
K = 0.0100
K = 0.2575
K = 0.5050
K = 0.7525
K = 1.0000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
horizontal velocity reduces
with increasing Coriolis force
(a) primary velocity with Coriolis force
0 5 10 15
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
g'( )
K = 0.0100
K = 0.2500
K = 0.5000
K = 0.7500
K = 1.0000
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
vertical velocity increases
with increasing Coriolis force
(b) secondary velocity with Coriolis force
Figure 4: Coriolis force eect
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
297
Volume 3, 2024
012345678
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
( )
Q = 0.1000
Q = 1.3250
Q = 2.5500
Q = 3.7750
Q = 5.0000
Temperature increases
with increasing heat
source
c = 0.5; K = 0.2;
M = 1; Pr = 6.2;
fw = 0.3; Q = 2;
Figure 5: temperature with heat source
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Parameters
3
3.1
3.2
3.3
3.4
3.5
3.6
3.7
3.8
Cfx
K
c
Q
(a) skin friction in
x
-direction
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Parameters
1.2
1.4
1.6
1.8
2
2.2
2.4
2.6
2.8
3
Cfy
K
c
Q
(b) skin friction in
y
-direction
Figure 6: skin friction
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Parameters
4
6
8
10
12
14
16
18
20
22
Cfy
K
c
Q
Figure 7: heat transfer rate
International Journal of Applied Sciences & Development
DOI: 10.37394/232029.2024.3.26
Oluwaseun Biodun Onuoha, Fawwaz Batayneh,
Abayomi Samuel Oke, Mario Raso
E-ISSN: 2945-0454
298
Volume 3, 2024