Next Article in Journal
A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method
Next Article in Special Issue
Experimental Analysis of the Thermal Performance Enhancement of a Vertical Helical Coil Heat Exchanger Using Copper Oxide-Graphene (80-20%) Hybrid Nanofluid
Previous Article in Journal
Preliminary Study on the Influence of the Polyphenols of Different Groups on the Digestibility of Wheat Starch, Measured by the Content of Resistant Starch
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives

1
Department of Industrial Engineering, Jazan University, Jazan 82822, Saudi Arabia
2
Department of Mathematics, University of Engineering and Technology, Lahore 54890, Pakistan
3
Department of Mathematics, COMSATS University Islamabad, Sahiwal 57000, Pakistan
4
Department of Mathematics and Statistics, Riphah International University I-14, Islamabad 44000, Pakistan
5
Department of Mechanical Engineering, Lebanese American University, Beirut 1102 2801, Lebanon
6
Mechanical Engineering Department, College of Engineering and Islamic Architecture, Umm Al-Qura University, P.O. Box 5555, Makkah 21955, Saudi Arabia
7
Department of Mechanical Engineering, College of Engineering, University of Ha’il, Ha’il 81451, Saudi Arabia
8
Laboratory of Thermal and Energy Systems Studies, National School of Engineering of Monastir, University of Monastir, Monastir 5000, Tunisia
9
Center of Research, Faculty of Engineering, Future University in Egypt, New Cairo 11835, Egypt
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(21), 10860; https://doi.org/10.3390/app122110860
Submission received: 10 September 2022 / Revised: 19 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Engineering Applications of Nanofluids)

Abstract

:
A fractional model was developed for presenting the thermal assessment of nanoparticles in an inclined moving surface. Water was used as a base fluid, while the nanofluid utilized copper oxide and silver nanoparticles. The modification of the thermal model was further supported by mixed convection, magnetic force, and porous saturated space. Slip effects to the porous surface were also introduced. The fluctuation in temperature at different times was assumed by following the ramped thermal constraints. The fractional computations for the set of flow problems were performed with implementations of the Atangana–Baleanu (AB) and Caputo–Fabrizio (CF) analytical techniques. The integration process for such computations was achieved using the Laplace transformation. The comparative velocity and thermal analysis for the water and kerosene-oil-based nanofluid model is presented. The declining change in the velocity was observed due to the increase in the volume fraction of nanoparticles. It was observed that the increment in the temperature profile was more progressive for the kerosene oil and silver nanoparticle suspension.

1. Introduction

Because of the rapid advancement of technology in all areas of life, global energy consumption continually increases. On the other hand, the use of fossil fuels has been reduced, and alternative energy resources are still being researched in order to improve their efficiency. Furthermore, fossil fuels have contributed to environmental damage and global warming. Nanotechnology’s revolution and its distinctive qualities in comparison to the huge scale of its originality have received a lot of attention. The multiple applications of nanotechnology in numerous spheres of life, such as the medical, engineering, industrial, and agricultural sectors, fueled this rapid rise. Nanotechnology is the study of the properties of tiny materials as a scientific field. Small nanoparticles of some solid materials with strong heat conductivity, such as titanium oxide and alumina, could be used to create nanotechnology-based approaches. When these nanometer-sized particles are suspended with a poor thermal conductivity base fluid, the heat transfer characteristics of the base fluids are likely to improve. This innovative fluid is known as a nanofluid, and it has a novel heat transfer property due to nanotechnology. Saving energy is, after all, the same as reducing the size of heat transfer equipment. Nanofluids have unique features that could make them beneficial in a variety of heat transfer applications, such as domestic refrigerators, pharmaceutical processes, machining, fuel cells, engine cooling, grinding, and so on. Choi [1] was the pioneer in investigating nanofluids as coolants and cooling agents. Recently, Mahdavi et al. [2] numerically studied the impact of the Poiseuille nanofluid flow between two flat plates and concluded that Brownian motion plays a crucial role in nanoparticle migration, especially when thermophoresis is involved. Dawar et al. [3] explored the effects of radiation on the heat transfer of a 3D MHD nanofluid across a spinning disk and found out that the velocity profiles decreased as the magnetic parameter increased. Using a nonuniform dynamical model, Uddin et al. [4] investigated the convective nanofluid flow inside a square cavity. Ghasemi et al. [5] deliberated the effects of the MHD and nonlinear radiation on the nanofluid flow across a stretched surface at the stagnation point. Entropy production phenomena were numerically examined using the activation energy and s heat source/sink by Li et al. [6]. Usman et al. [7] theoretically studied the heat and mass transport phenomena in a nanofluid flow upon a radially spinning disk filled with microorganisms. Mohsenian et al. [8] addressed the convective heat transfer in an MHD nanofluid flow inside two tubes. Klazly and Bognar [9] enhanced the heat transfer in a nanofluid flow in a channel-shaped geometry. Acharya et al. [10] analyzed the time-dependent MHD flow of a nanofluid upon a spinning disk. Muthtamilselvan et al. [11] reported nanofluid properties with microorganisms. Al-Mdallal et al. [12] addressed the copper nanofluid thermal impact for the spinning of disks.
Fractional calculus is the interesting branch of mathematics for which work was initiated after 1695 for suggesting the extension of integer derivatives [13]. Fractional models often result due to the formulation of many physical and scientific problems, such as quantum physics, plasma, chemistry, biological sciences, engineering, and industrial problems. The main motivations for addressing fractional simulations are the presence of singularities and local and nonlocal kernels. Different definitions regarding the fractional framework have been suggested by researchers in fractional calculus. The Caputo and Fabrizio technique [14] is one such mathematical fractional tool that can help to access the solutions of fractional problems. The definition of a CF operator is based on the exponential decay law. The novel aspect of the CF definition is the presence of a nonsingular kernel. On other hand, the Atangana–Baleanu (AB) derivative is another useful fractional tool that was recently implemented by many researchers [15]. The expression and implementation of the AB model is based on the definition of the Mittag–Leffler hypothesis. Investigators in recent decades [16,17,18,19,20,21,22,23,24,25] have focused a number of studies on some new techniques.
After presenting the above-mentioned research survey, the aim of the current continuation is to develop a fractional mathematical model for the nanofluid problem in the presence of base fluids. Kerosene oil and water base fluids are used with suspensions of cooper oxide C u O and silver A g nanoparticles. The motivations for enhancing the thermal efficiency of water and kerosene oil materials are subject to various engineering and industrial applications where such materials are comprehensively used as fuel. Kerosene oil is a distillate that is commonly used as a solvent and fuel in many industrial processes such as coal, oil shale, and petroleum refining, etc. An inclined surface was subjected to a magnetic force, and a porous medium was used for thermal inspection. The fluctuation in temperature with ramped conditions with flow time constraints is presented. The analytical simulations are presented with help of two fractional approaches: the Atangana–Baleanu (AB) and Caputo time-fractional (CF) derivatives. Both techniques are supported with recent definitions, and computations were performed via MATLAB software. The thermal dynamic the of parameters is physically addressed.

2. Problem Formulation

The thermal features of a nanofluid in the presence of kerosene oil and water base particles were studied. An inclined surface flow containing a uniform suspension of copper oxide and silver nanoparticles was considered. The magnetic force was normally implemented with an angle of inclination, θ . The impact of a porous medium was subjected to the flat surface (Figure 1). The suspension of nanoparticles and base materials with a uniform (50-50%) ratio was taken into consideration. The surface was further supported with slip constraints to control the fluid velocity. The flow in the stationary plate was induced after time t > 0 . The governing model is presented in terms of the following mathematical expressions [26,27,28]:
ρ n f w ξ , t t = μ n f 1 + α 1 t 2 w ξ , t ξ 2 + g ρ β T n f T ξ , t T C o s δ μ n f φ K w ξ , t σ n f B o 2 S i n θ w ξ , t ;             ξ , t > 0
ρ C p n f T ξ , t t = q ξ       ;                   ξ , t > 0
q ξ , t = k n f T ξ , t ξ
The boundary conditions for the fractional model are [26,27,28]:
w ξ , 0 = 0 ,                 T ξ , 0 = T               ;             ξ > 0
w 0 , t b w ξ , t ξ ξ = 0 = g t μ ,                 T 0 , t = T + T w T t t o ,                     0 t t o   T w ,                                                                                         t > t o
w ξ , t 0 ,                 T ξ , t T       ;             ξ ,   t > 0
The problem was entertained with following new variables [26,27,28]:
ξ * = ξ ν o υ f   ,         w * = w U o   ,           t * = ν o 2 t υ f     ,             ϑ * = T T T w T b * = h k b ,             q * = q q o ,           q o = k n f T w T ν o υ f           ,   g * t * = 1 μ t o υ f t o t *
In view of the above new quantities, the problem was modified as:
w ξ , t t = 1 Λ o Λ 1 1 + β 1 t 2 w ξ , t ξ 2 + Λ 2 Λ o G r   C o s δ   ϑ ξ , t 1 Λ o M   S i n θ + K e f f w ξ , t
Λ 3 P r ϑ ξ , t t = q ξ , t ξ       ;                   ξ , t > 0
q ξ , t = Λ 4 ϑ ξ , t ξ
with:
w ξ , 0 = 0 ,                 ϑ ξ , 0 = 0               ;             ξ > 0
w 0 , t b w ξ ξ = 0 = g t ,           ϑ 0 , t = t ,                     0 < t 1   1 ,                                   t > 1  
w ξ , t 0 ,                 ϑ ξ , t 0       ;             ξ ,   t > 0
where:
ρ n f = 1 φ ρ f + φ ρ s     ,               μ n f = μ f 1 φ 2.5 ρ β n f = 1 φ ρ β f + φ ρ β s     ,               ρ C p n f = 1 φ ρ C p f + φ ρ C p s k n f k f = k s + 2 k f 2 φ k f k s k s + 2 k f + 2 φ k f k s ,             σ n f σ f = 1 + 3 σ s σ f 1 φ σ s σ f + 2 σ s σ f 1 φ 1 Λ o = 1 φ + φ ρ s ρ f   ,               Λ 1 = 1 1 φ 2.5   ,                   Λ 2 = 1 φ + φ ρ β T s ρ β T f Λ 3 = 1 φ + φ ρ C p s ρ C p f   ,               Λ 4 = k n f k f ,             P r = μ C p f k f G r = g β ν f T w T U o 3 ,     M = υ f ν o 2 σ n f B o 2 ρ f υ f   ,         β 1 = α 1 υ f υ f ν o 2
The classical thermal features of water, kerosene oil, copper oxide, and silver particles are presented in Table 1.

3. Basics of Fractional Simulations

Definition 1.
The definition of Atangana–Baleanu (AB) for the g ξ , t   function is presented as:
D   A B t β g ξ , t = 1 1 β 0 t E β β t z β 1 β g ξ , t d t     ;     0 < β < 1
where the AB operator is D   A B t β . The Mittage-Leffer function, E β z , is explained as:
E β z = r = 0 z β Γ r β + 1     ;      0 < β < 1 ,     z
Definition 2.
For Equation (13), the definition of the Laplace transform is:
L D   A B t β g ξ , t = q β L g ξ , t q β 1 g ξ , 0 1 β q β + β
and
lim β 1 D   A B t β g ξ , t = g ξ , t t
Definition 3.
For any function  g ξ , t , the definition of the Caputo fractional (CF) is associated with the following definition:
D   C F t α g ξ , t = 1 1 α 0 t e x p α 1 τ 1 α g ξ , t d τ 0 < α < 1
where   D   C F t α is the Caputo–Fabrizio operator.
Definition 4.
The Laplace transformation for CF is:
L D   C F t α g ξ , t = s L g ξ , t g ξ , 0 1 α s + α
and with the Laplace-transformed variable s of time   t   for CF, we obtain:
lim α 1 D   C F t α g ξ , t = g ξ , t t

4. Solution by Atangana–Baleanu (AB) Operator

The implementation of the AB fractional model with operator D   A B t β is suggested for Equation (9) as:
q ξ , t = Λ 4 D   A B t β ϑ ξ , t ξ

4.1. AB Operator for Temperature Profile

Using the Laplace transform in Equations (8) and (17) yields:
Λ 3 Pr q ϑ ¯ ξ , q = q ¯ ξ , q ξ       ;                   ξ , t > 0
q ¯ ξ , q = Λ 4 q β 1 β q β + β ϑ ¯ ξ , q ξ
Introducing Equation (19) into Equation (18) yields:
2 ϑ ¯ ξ , q ξ 2 Λ 3 P r Λ 4 q 1 β 1 β q β + β ϑ ¯ ξ , q = 0      
where   ϑ ¯ ξ , q   is the Laplace transform of function   ϑ ξ , t with the solution:
ϑ ¯ 0 , q = 1 e q q 2                         a n d             ϑ ¯ ξ , q 0     a s         ξ
The solution via the AB approach for the temperature profile is:
ϑ ¯ ξ , q = 1 e q q 2 e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β
In a more extended form, Equation (21) can be written as:
ϑ ¯ ξ , q = 1 e q q 2 e ξ c 1 q + c 2 q 1 β
where:
ϰ 1 = Λ 3 P r   1 β Λ 4 ,                     ϰ 2 = β 1 β
Using the exponential function for determining the Laplace inverse of Equation (22), one obtains:
ϑ ¯ ξ , q = h ¯ 1 q * 1 q + a 1 = 1 a 2 = 0 ξ ϰ 1 a 1 a 1 !   a 2 ! ϰ 2 a 2 q 1 + a 2   β a 1 2 Γ a 1 2 + 1 Γ a 1 2 a 2 + 1
The solution for temperature change by following the Laplace inverse is:
ϑ ξ , t = h 1 t * 1 + a 1 = 1 a 2 = 0 ξ ϰ 1 a 1 a 1 !   a 2 ! ϰ 2 a 2 t a 2   β a 1 2 Γ 1 + a 2   β a 1 2 Γ a 1 2 + 1 Γ a 1 2 a 2 + 1 h 1 t = L 1 1 e q q = 1 H t 1
where   H .   is the Heaviside function.

4.2. Velocity Field via AB-Fractional Derivative

Using the Laplace transform on Equation (7), the equation becomes:
q   w ¯ ξ , q = 1 Λ o Λ 1 1 + β 1 q 2 w ¯ ξ , q ξ 2 + Λ 2 Λ o G r   C o s δ ϑ ¯ ξ , q 1 Λ o M   S i n θ + K e f f w ξ , t
with the conditions:
w ¯ 0 , q b w ¯ ξ , q ξ ξ = 0 = G q     and         w ¯ ξ , q 0         a s         ξ
The solution is:
w ¯ ξ , q = 1 1 + b q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q
( Λ 2 G r   C o s δ Λ o q 1 Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 ( 1 + b Λ 3 P r Λ 4 q 1 β 1 β q β + β ) + G q ) e ξ q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 Λ 2 G r   C o s δ Λ o q e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q

5. Modeling via CF Time-Fractional Derivative

In the above section, the simulations for the thermal and momentum profiles were found by utilizing CF-fractional derivatives. Now, the solution by a CF-fractional derivative can be formulated by using the CF-fractional derivative operative   D   C F t α in Fourier’s law of thermal conductivity as follows:
q ξ , t = Λ 4 D   C F t α ϑ ξ , t ξ

5.1. Energy Field via CF-Fractional Derivative

As the thermal profile is autonomous of the momentum profile, the simulations of the thermal field can originate from Equations (8) and (26) by utilizing the Laplace transformation technique as follows:
Λ 3 Pr s ϑ ¯ ξ , s = q ¯ ξ , s ξ       ;                   ξ , t > 0
q ¯ ξ , s = Λ 4 s 1 α s + α ϑ ¯ ξ , s ξ
where   ϑ ¯ ξ , s is the Laplace transformation of   ϑ ξ , t with the following conditions:
ϑ ¯ 0 , s = 1 e s s 2                 a n d             ϑ ¯ ξ , s 0     a s         ξ
Now, presenting Equation (28) into Equation (27), we obtain:
2 ϑ ¯ ξ , s ξ 2 Λ 3 P r Λ 4 1 α s + α ϑ ¯ ξ , s = 0    
By inputting the conditions, we yield:
ϑ ¯ ξ , s = 1 e s s 2 e ξ Λ 3 P r Λ 4 1 α s + α
Equation (29), in a more general form, can be formulated as:
ϑ ¯ ξ , s = 1 e s s e ξ d 1 s + d 2 s
where:
1 = Λ 3 P r   1 α Λ 4 ,                         2 = α 1 α
For the Laplace inverse of Equation (30), we use the summation form of the exponential function as follows:
ϑ ¯ ξ , s = h ¯ 1 s * 1 s + a 1 = 1 a 2 = 0 ξ 1 a 1 a 1 !   a 2 ! s a 2 1 2 a 2 a 1 2 Γ a 1 2 + 1 Γ a 1 2 a 2 + 1
By utilizing the Laplace inverse, we yield:
ϑ ξ , t = h 1 t * 1 + a 1 = 1 a 2 = 0 ξ 1 a 1 a 1 !   a 2 !   t a 2 2 a 2 a 1 2   Γ a 1 2 + 1 Γ a 1 2 a 2 + 1 Γ 1 a 2
which is the required solution of the temperature field by the CF-fractional derivative.

5.2. Velocity Field via CF-Fractional Derivative

Taking the Laplace transform of velocity field Equation (8) and its corresponding conditions, such that:
s w ¯ ξ , s = 1 Λ o Λ 1 1 + β 1 s 2 w ¯ ξ , s ξ 2 + Λ 2 Λ o G r   ϑ ¯ ξ , s C o s δ 1 Λ o M   S i n θ + K e f f w ξ , t
with the succeeding transformed conditions:
w ¯ 0 , s b w ¯ ξ , s ξ ξ = 0 = G s         and           w ¯ ξ , s 0         a s         ξ
the following solution of Equation (31) can be yielded by inputting the respective conditions in Equation (32):
w ¯ ξ , s = 1 1 + b s + K e f f Π 1 + Π 2 1 + s β 1 Π 3 s Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 1       + b Π 4 1 α s + α + G s e ξ s + K e f f Π 1 + Π 2 1 + s β 1       1 e s s 2 Π 3 e ξ Λ 3 P r Λ 4 1 α s + α Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1  
Π 1 =   Λ o Λ 1 ,                   Π 2 = Λ 1 M S i n θ                       Π 3 = Λ 2 G r   C o s δ Λ o ,                   Π 4 = Λ 3 P r Λ 4
As for the Laplace inverse, many authors have utilized different numerical schemes, so here, in this grinding stone, we also utilized the Stehfest and Tzous schemes for the solution of momentum fields attained by both the AB and CF-fractional definitions. The mathematical formulas of both algorithms can be summarized as follows [29,30,31,32]:
w ξ , t = l n 2 t n = 1 M v n w ¯ ξ , n l n 2 t
where M is a positive integer, and
v n = 1 n + M 2 p = q + 1 2 m i n q , M 2 p M 2 2 p ! M 2 p ! p !   p 1 !   q p !   2 p q !
and
w ξ , t = e 4.7 t 1 2 w ¯ r , 4.7 t + R e   j = 1 N 1 k w ¯ r , 4.7 + k π i t

6. Special Cases

Because the dynamics velocity profile solutions via the AB and CF definitions in Equations (25) and (33) are more comprehensive, in order to provide some further physical understanding of the problem, we address certain exceptional situations for the momentum profile function   g t , whose physical connotation is well-documented in the literature.
Case 1: g t = t
If we assume   g t = t , the equations for the velocity field via AB and the CF derivative using Equations (25) and (33) are as follows.
w ¯ ξ , q = 1 1 + b q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q
( Λ 2 G r   C o s δ Λ o q 1 Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 ( 1 + b Λ 3 P r Λ 4 q 1 β 1 β q β + β ) + 1 q 2 ) e ξ q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 Λ 2 G r   C o s δ Λ o q e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q  
and
w ¯ ξ , s = 1 1 + b s + K e f f Π 1 + Π 2 1 + s β 1 ( Π 3 s Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 ( 1 + b Π 4 1 α s + α ) + 1 s 2 ) e ξ s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 Π 3 e ξ Λ 3 P r Λ 4 1 α s + α Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1  
Case 2: g t = S i n ω t
Assuming   g t = S i n ω t   ( ω   is the shear stress strength), the velocity solution is:
w ¯ ξ , q = 1 1 + b q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q
Λ 2 G r   C o s δ Λ o q 1 Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 1            + b Λ 3 P r Λ 4 q 1 β 1 β q β + β + ω ω 2 + q 2   e ξ q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q            1 e q q 2 Λ 2 G r   C o s δ Λ o q e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q  
and
w ¯ ξ , s = 1 1 + b s + K e f f Π 1 + Π 2 1 + s β 1 Π 3 s Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 1      + b Π 4 1 α s + α + ω ω 2 + s 2   e ξ s + K e f f Π 1 + Π 2 1 + s β 1      1 e s s 2 Π 3 e ξ Λ 3 P r Λ 4 1 α s + α Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1  
Case 3: g t = t   C o s t
Replacing the   g t = t   C o s   t   with the Laplace transform   G q = q 2 1 q 2 + 1 2 , the velocity determination in view of Equations (25) and (33) is:
w ¯ ξ , q = 1 1 + b q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q
( Λ 2 G r   C o s δ Λ o q 1 Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 ( 1 + b Λ 3 P r Λ 4 q 1 β 1 β q β + β ) + q 2 1 1 + q 2 2 )   e ξ q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 Λ 2 G r   C o s δ Λ o q e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q  
and
w ¯ ξ , s = 1 1 + b s + K e f f Π 1 + Π 2 1 + s β 1 ( Π 3 s Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 ( 1 + b Π 4 1 α s + α ) + s 2 1 1 + s 2 2 )   e ξ s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 Π 3 e ξ Λ 3 P r Λ 4 1 α s + α Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1  
Case 4: g t = t   e t
Finally, for g t = t e t with Laplace G q = 1 q 1 2 , the velocity solution for the AB and CF models in view of Equations (25) and (33) leads to:
w ¯ ξ , q = 1 1 + b q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q
Λ 2 G r   C o s δ Λ o q 1 Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q 1 e q q 2 1          + b Λ 3 P r Λ 4 q 1 β 1 β q β + β + 1 q 1 2   e ξ q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q          1 e q q 2 Λ 2 G r   C o s δ Λ o q e ξ Λ 3 P r Λ 4 q 1 β 1 β q β + β Λ 3 P r Λ 4 q 1 β 1 β q β + β q + K e f f Λ 0 Λ 1 + Λ 1 M S i n θ 1 + β 1 q  
and
w ¯ ξ , s = 1 1 + b s + K e f f Π 1 + Π 2 1 + s β 1 Π 3 s Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1 1 e s s 2 1     + b Π 4 1 α s + α + 1 s 1 2 e ξ s + K e f f Π 1 + Π 2 1 + s β 1     1 e s s 2 Π 3 e ξ Λ 3 P r Λ 4 1 α s + α Π 4 1 α s + α s + K e f f Π 1 + Π 2 1 + s β 1  

7. Validation of Results

The fractional computations are verified in Figure 2 with the work of Basit et al. [30]. It is clearly noted that the obtained results show fine accuracy with these results.

8. Physical Analysis of Results

After performing the fractional simulations, the physical reflection of the parameters is presented in this section. For assessing the physical dynamic, numerical values were assigned to physical parameters, such as P r = 6.2 ,   φ = 0.02 ,   t = 1.0 ,   α = 0.5 , β = 0.5 ,   P r = 2.4 ,   G r = 2.4 ,   θ = π 4 ,   w = 0.9 ,   b = 0.4 ,   and   δ = π 4 . Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 illustrate the impacts of various parameters on the velocity and temperature field dynamics. It was remarked that, for nonlocal derivatives, we could analyze only one profile at α ,   β = 1 , but with the virtue of fractional-order derivative models, more than one profile of the subsequent problems could be analyzed. Moreover, the non-integer-order derivatives models could explain different layers for flowing fluid. Furthermore, due to this advantage of fractional models, the comparisons of different fluid layers could be examined individually. Figure 3 presents the comparative analysis of the velocity profile against the water and kerosene-oil base fluid with the suspension of copper oxide and silver nanoparticles. The water-based nanofluid showed more enhancement in the velocity profile compared to the kerosene-oil-based nanoparticle. Furthermore, the nanofluid mixed with copper oxide   C u O   nanoparticles had a more significant impact on the velocity compared to the silver   A g   nanoparticles. Figure 4 reports the changes in velocity for α and β . An increasing influence of velocity due to both factors was observed. The influence of the volume fraction parameter,   φ ,   on velocity is highlighted Figure 5. The declining velocity change was noted for increasing the volume fraction parameter. The effects of parameter   M   on the velocity field are illustrated in Figure 6. As expected, the magnetic parameter slowed the fluid moment due to the influence of drag force, namely the Lorentz force. The Lorentz force could produce some resistance to fluid velocity due to resistive behavior.
The comparison of different nanofluids and nanoparticles for the thermal profile is plotted in Figure 7a. It is clear from the figure that the kerosene-oil-based nanofluid had a greater tendency to enhance the thermal conductivity than the water-based nanofluid. Figure 7b illustrates the impact of the Prandtl number,   P r ,   on the temperature profile. A falling behavior of the temperature profile can be seen against a larger Prandtl constant. Physically, with the increment in   Pr   values, the thermal forces were weakened. Figure 7c aims to present the thermal aspect of the fractional parameters α and β on the temperature profile. The temperature profile rose for fractional parameters α and β .  Table 2 aims to address the numerical fluctuations in the velocity and temperature profile by using the Stehfest [31] and Tzou’s [32] integration tools. The lower temperature rate due to a larger space parameter, ξ , was noted. However, the velocity increased for larger ξ values using both techniques. After careful inspection, it was noted that the numerical results obtained via Tzou’s [32] technique had better relative convergence compared to Stehfest’s [31] tools. Table 3 presents the numerical values of the Nusselt number and the wall shear force. With an increment in ξ , the Nusselt number decreased, while an increasing impact of skin friction was observed.

9. Concluding Remarks

The fractional model for the nanofluids with kerosene oil and water base materials is presented. The thermal endorsement of the nanofluid particles was additionally impacted by slip effects, mixed convection, and magnetic forces. The computations were simulated with the Atangana–Baleanu (AB) and Caputo–Fabrizio (CF) fractional tools. The main results from current investigation are:
The thermal characteristics of heat transfer could be magnificently improved in the presence of copper oxide and silver nanoparticles.
A relatively improved profile of velocity was observed for the copper–water-based suspension compared to the kerosene-oil suspension.
A lower velocity rate was observed for the volume fraction and increasing values of fractional parameters.
The increasing thermal outcomes were predicted for the suspension of kerosene oil and silver nanoparticles compared to copper oxide and water-based material.
Decaying thermal outcomes were observed when the fractional parameters and Prandtl number were varied.
The numerical values for the temperature profile were more impressive compared to the velocity profile using Stehfest [31] and Tzou’s [32] integration techniques.
These results can be further extended by incorporating various thermal features such as entropy generation, exponential heat sources, activation energy, joule heating, and bioconvection and by using different non-Newtonian models.

Author Contributions

Conceptualization, O.T.B., S.U.K.; methodology, A.R.; software, A.R.; validation, M.I.K., A.N.; formal analysis, N.B.K.; investigation, S.U.K.; resources, M.I.K.; data curation, O.T.B.; writing—original draft preparation, A.R., S.U.K.; writing—review and editing, M.I.K., E.S.M.T.-E.; visualization, E.S.M.T.-E.; supervision, M.I.K.; project administration, M.I.K., S.U.K.; funding acquisition, E.S.M.T.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are clearly available in the manuscript.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work (grant code: 22UQU4340531DSR04).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

w Velocity m / s
t Time s
U o Constant Velocity m / s
ρ n f Density Kg / m 3
T Temperature K
k n f Thermal conductivity W / m 2 k
P r Prandtl number
G r Grashof number
μ n f Dynamic viscosity Kg / ms
ρ s Density of nanoparticles Kg / m 3
M Magnetic parameter
n f Nanofluid
β 1 slip parameter
β AB-fractional derivative operator
α CF-fractional derivative operator
N u Nusselt number
C f Skin Friction
CFCaputo–Fabrizio fractional derivative
q Laplace variable by AB
s Laplace variable by CF
g Acceleration due to gravity m / s 2
ABAtangana–Baleanu fractional derivative
ρ f Density of base fluid Kg / m 3

References

  1. Choi SUS 1995 Enhancing thermal conductivity of fluids with nanoparticles Proc Int Mech Eng Congress San Francisco USA ASME FED 231/MD66 99–105.
  2. Mahdavi, M.; Sharifpur, M.; Ahmadi, M.H.; Meyer, J.P. Nanofluid flow and shear layers between two parallel plates: A simulation approach. Eng. Appl. Comput. Fluid Mech. 2020, 14, 1536–1545. [Google Scholar] [CrossRef]
  3. Dawar, A.; Wakif, A.; Thumma, T.; Shah, N.A. Towards a new MHD non-homogeneous convective nanofluid flow model for simulating a rotating inclined thin layer of sodium alginate-based Iron oxide exposed to incident solar energy. Int. Commun. Heat Mass Transf. 2022, 130, 105800. [Google Scholar] [CrossRef]
  4. Uddin, M.J.; Rasel, S.K.; Adewole, J.K.; al Kalbani, K.S. Finite element simulation on the convective double diffusive water-based copper oxide nanofluid flow in a square cavity having vertical wavy surfaces in presence of hydro-magnetic field. Results Eng. 2022, 13, 100364. [Google Scholar] [CrossRef]
  5. Ghasemi, S.E.; Mohsenian, S.; Gouran, S.; Zolfagharian, A. A novel spectral relaxation approach for nanofluid flow past a stretching surface in presence of magnetic field and nonlinear radiation. Results Phys. 2022, 32, 105141. [Google Scholar] [CrossRef]
  6. Li, Y.M.; Al-Khaled, K.; Gouadria, S.; El-Zahar, E.R.; Usman Khan, S.U.; Khan, M.I.; Malik, M.Y. Numerical simulations for three-dimensional rotating porous disk flow of viscoelastic nanomaterial with activation energy, heat generation and Nield boundary conditions. Waves Random Complex Media 2022. [Google Scholar] [CrossRef]
  7. Usman; Bhatti, M.M.; Ghaffari, A.; Doranehgard, M.H. The role of radiation and bioconvection as an external agent to control the temperature and motion of fluid over the radially spinning circular surface: A theoretical analysis via Chebyshev spectral approach. Math. Meth. Appl. Sci. 2022, 1–18. [Google Scholar] [CrossRef]
  8. Mohsenian, S.; Gouran, S.; Ghasemi, S.E. Evaluation of weighted residual methods for thermal radiation on nanofluid flow between two tubes in presence of magnetic field. Case Stud. Therm. Eng. 2022, 32, 101867. [Google Scholar] [CrossRef]
  9. Klazly, M.; Bognar, G. Heat transfer enhancement for nanofluid flows over a microscale backward-facing step. Alex. Eng. J. 2022, 61, 8161–8176. [Google Scholar] [CrossRef]
  10. Acharya, N.; Mabood, F.; Shahzad, S.A.; Badruddin, I.A. Hydrothermal variations of radiative nanofluid flow by the influence of nanoparticles diameter and nanolayer. Int. Commun. Heat Mass Transf. 2022, 130, 105781. [Google Scholar] [CrossRef]
  11. Muthtamilselvan, M.; Suganya, S.; Al-Mdallal, Q.M. Stagnation-Point Flow of the Williamson Nanofluid Containing Gyrotactic Micro-organisms. Proc. Natl. Acad. Sci. India Sect. A Phys. Sci. 2021, 91, 633–648. [Google Scholar] [CrossRef]
  12. Al-Mdallal, Q.M.; Renuka, A.; Muthtamilselvan, M.; Abdalla, B. Ree-Eyring fluid flow of Cu-water nanofluid between infinite spinning disks with an effect of thermal radiation. Ain Shams Eng. J. 2021, 12, 2947–2956. [Google Scholar] [CrossRef]
  13. Milici, C.; Draganescu, G.; Machado, J.T. Introduction to Fractional Differential Equations; Springer Nature Switzerland AG: Cham, Switzerland, 2019. [Google Scholar]
  14. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  15. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  16. Waqas, H.; Oreijah, M.; Guedri, K.; Khan, S.U.; Yang, S.; Yasmin, S.; Khan, M.I.; Bafakeeh, O.T.; Tag-ElDin, E.S.M.; Galal, A.M. Gyrotactic motile microorganisms impact on pseudoplastic nanofluid flow over a moving Riga surface with exponential heat flux. Crystals 2022, 12, 1308. [Google Scholar] [CrossRef]
  17. Shahid, M.; Javed, H.M.A.; Ahmad, M.I.; Qureshi, A.A.; Khan, M.I.; Alnuwaiser, M.A.; Ahmed, A.; Khan, M.A.; Tag-Eldin, E.; Shahid, A.; et al. A brief assessment on recent developments in efficient electrocatalytic Nitrogen reduction with 2D non-metallic nanomaterials. Nanomaterials 2022, 12, 3413. [Google Scholar] [CrossRef]
  18. Manzoor, N.; Qasim, I.; Khan, M.I.; Ahmed, M.W.; Guedri, K.; Bafakeeh, O.T.; Tag-Eldin, E.S.M.; Galal, A.M. Antibacterial applications of low pressure plasma on degradation of multidrug resistant V. cholera. Appl. Sci. 2022, 12, 9737. [Google Scholar] [CrossRef]
  19. Mamatha, S.U.; Devi, R.L.V.R.; Ahammad, N.A.; Shah, N.A.; Rao, B.M.; Raju, C.S.K.; Khan, M.I.; Guedri, K. Multi-linear regression of triple diffusive convectively heated boundary layer flow with suction and injection: Lie group transformations. Int. J. Mod. Phys. B 2022, in press. [CrossRef]
  20. Kiranakumar, H.V.; Thejas, R.; Naveen, C.S.; Khan, M.I.; Prasanna, G.D.; Reddy, S.; Oreijah, M.; Guedri, K.; Bafakeeh, O.T.; Jameel, M. A review on electrical and gas-sensing properties of reduced graphene oxide-metal oxide nanocomposites. Biomass Convers. Biorefinery 2022, in press. [CrossRef]
  21. Abbasi, A.; Farooq, W.; Tag-ElDin, E.S.M.; Khan, S.U.; Khan, M.I.; Guedri, K.; Elattar, S.; Waqas, M.; Galal, A.M. Heat transport exploration for hybrid nanoparticle (Cu, Fe3O4)-based blood flow via tapered complex wavy curved channel with slip features. Micromachines 2022, 13, 1415. [Google Scholar] [CrossRef]
  22. Wahid, N.S.; Arifina, N.M.; Khashi’ie, N.S.; Pop, I.; Bachok, N.; Hafidzuddin, M.E.H. MHD mixed convection flow of a hybrid nanofluid past a permeable vertical flat plate with thermal radiation effect. Alex. Eng. J. 2022, 61, 3323–3333. [Google Scholar] [CrossRef]
  23. Ashorynejad, H.R.; Shahriari, A. MHD natural convection of hybrid nanofluid in an open wavy cavity. Results Phys. 2018, 9, 440–455. [Google Scholar] [CrossRef]
  24. Puneetha, V.; Khan, M.I.; Jameel, M.; Geudri, K.; Galal, A.M. The convective heat transfer analysis of the casson nanofluid jet flow under the influence of the movement of gyrotactic microorganisms. J. Indian Chem. Soc. 2022, 9, 100612. [Google Scholar] [CrossRef]
  25. Oztop, H.F.; Abu-Nada, E. Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  26. Raza, A.; Ghaffari, A.; Khan, S.U.; Haq, A.U.; Khan, M.I.; Khan, M.R. Non-singular fractional computations for the radiative heat and mass transfer phenomenon subject to mixed convection and slip boundary effects. Chaos Solitons Fractals 2022, 155, 111708. [Google Scholar] [CrossRef]
  27. Raza, A.; Khan, S.U.; Al-Khaled, K.; Khan, M.I.; Haq, A.U.; Alotaibi, F.; Mousa, A.A.A.; Qayyum, S. A fractional model for the kerosene oil and water-based Casson nanofluid with inclined magnetic force. Chem. Phys. Lett. 2022, 787, 139277. [Google Scholar] [CrossRef]
  28. Raza, A.; Khan, S.U.; Khan, M.I.; Farid, S.; Muhammad, T.; Khan, M.I.; Galal, A.M. Fractional order simulations for the thermal determination of graphene oxide (GO) and molybdenum disulphide (MoS2) nanoparticles with slip effects. Case Stud. Therm. Eng. 2021, 28, 101453. [Google Scholar] [CrossRef]
  29. Guo, B.; Raza, A.; Al-Khaled, K.; Khan, S.U.; Farid, S.; Wang, Y.; IjazKhan, M.; Malik, M.Y.; Saleem, S. Fractional-order simulations for heat and mass transfer analysis confined by elliptic inclined plate with slip effects: A comparative fractional analysis. Case Stud. Therm. Eng. 2021, 28, 101359. [Google Scholar] [CrossRef]
  30. Basit, A.; Asjad, M.I.; Akgül, A. Convective flow of a fractional second grade fluid containing different nanoparticles with Prabhakar fractional derivative subject to non-uniform velocity at the boundary. Math. Methods Appl. Sci. 2021. [Google Scholar] [CrossRef]
  31. Stehfest, H. Algorithm 368: Numerical inversion of Laplace transforms [D5]. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
  32. Tzou, D.Y. Macro-to Microscale Heat Transfer: The Lagging Behavior; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
Figure 1. Flow configuration.
Figure 1. Flow configuration.
Applsci 12 10860 g001
Figure 2. Comparison of the results with the work of Basit et al. [30].
Figure 2. Comparison of the results with the work of Basit et al. [30].
Applsci 12 10860 g002
Figure 3. Comparative velocity profile for water- and kerosene-based nanofluids.
Figure 3. Comparative velocity profile for water- and kerosene-based nanofluids.
Applsci 12 10860 g003
Figure 4. Change in velocity for 𝛼 and 𝛽 .
Figure 4. Change in velocity for 𝛼 and 𝛽 .
Applsci 12 10860 g004
Figure 5. Change in velocity for φ .
Figure 5. Change in velocity for φ .
Applsci 12 10860 g005
Figure 6. Change in velocity for M .
Figure 6. Change in velocity for M .
Applsci 12 10860 g006
Figure 7. (ac): (a) Comparative temperature profile for water- and kerosene-based nanofluids, (b) change in temperature for P r , and (c) change in temperature for β .
Figure 7. (ac): (a) Comparative temperature profile for water- and kerosene-based nanofluids, (b) change in temperature for P r , and (c) change in temperature for β .
Applsci 12 10860 g007
Table 1. The numerical values for water, kerosene oil, copper oxide, and silver nanoparticles.
Table 1. The numerical values for water, kerosene oil, copper oxide, and silver nanoparticles.
Materials ρ   K g / m 3 C p   J / K g   K k   w / m .   K β T × 10 5   1 / K
H 2 O 997.141790.61321
kerosene oil88419100.11470
C u O 6320531.876.51.80
A g 10,5002354291.89
Table 2. The numerical outcomes for temperature and velocity profiles via Stehfest [31] and Tzou’s [32] integration tools.
Table 2. The numerical outcomes for temperature and velocity profiles via Stehfest [31] and Tzou’s [32] integration tools.
ξ Temperature Change by Stehfest [31]Temperature Change by Tzou [32]Velocity Change by Stehfest [31]Velocity Change by Tzou [32]
0.10.93140.93130.51590.5130
0.30.63800.63790.57330.5703
0.50.43680.43670.57590.5730
0.70.29890.29880.54710.5443
0.90.20440.20430.50190.4993
1.10.13970.13970.44960.4473
1.30.09540.09540.39620.3940
1.50.06510.06510.34480.3429
1.70.04450.04440.29430.2957
1.90.03030.03030.25460.2532
Table 3. Numerical analysis for Nusselt number and skin friction by both numerical algorithms.
Table 3. Numerical analysis for Nusselt number and skin friction by both numerical algorithms.
ξ Nusselt Number by Stehfest [31]Nusselt Number by Tzou [32]Skin Friction by Stehfest [31]Skin Friction by Tzou [32]
0.10.7683530.7621760.9646970.9650723
0.30.7066660.7016320.9850970.9855834
0.50.6581850.6541531.24114431.247146
0.70.6142010.6111551.4108351.417757
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bafakeeh, O.T.; Raza, A.; Khan, S.U.; Khan, M.I.; Nasr, A.; Khedher, N.B.; Tag-Eldin, E.S.M. Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives. Appl. Sci. 2022, 12, 10860. https://doi.org/10.3390/app122110860

AMA Style

Bafakeeh OT, Raza A, Khan SU, Khan MI, Nasr A, Khedher NB, Tag-Eldin ESM. Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives. Applied Sciences. 2022; 12(21):10860. https://doi.org/10.3390/app122110860

Chicago/Turabian Style

Bafakeeh, Omar T., Ali Raza, Sami Ullah Khan, Muhammad Ijaz Khan, Abdelaziz Nasr, Nidhal Ben Khedher, and El Sayed Mohamed Tag-Eldin. 2022. "Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives" Applied Sciences 12, no. 21: 10860. https://doi.org/10.3390/app122110860

APA Style

Bafakeeh, O. T., Raza, A., Khan, S. U., Khan, M. I., Nasr, A., Khedher, N. B., & Tag-Eldin, E. S. M. (2022). Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives. Applied Sciences, 12(21), 10860. https://doi.org/10.3390/app122110860

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

Article Metrics

Back to TopTop