Next Article in Journal
Orthonormal Ultraspherical Operational Matrix Algorithm for Fractal–Fractional Riccati Equation with Generalized Caputo Derivative
Next Article in Special Issue
Oscillation Criteria of Solutions of Fourth-Order Neutral Differential Equations
Previous Article in Journal
Quarter-Sweep Preconditioned Relaxation Method, Algorithm and Efficiency Analysis for Fractional Mathematical Equation
Previous Article in Special Issue
Qualitative Behavior of Unbounded Solutions of Neutral Differential Equations of Third-Order
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advancement of Non-Newtonian Fluid with Hybrid Nanoparticles in a Convective Channel and Prabhakar’s Fractional Derivative—Analytical Solution

by
Muhammad Imran Asjad
1,
Noman Sarwar
1,
Muhammad Bilal Hafeez
2,
Wojciech Sumelka
2,* and
Taseer Muhammad
3
1
Department of Mathematics, University of Management and Technology, Lahore 54770, Pakistan
2
Institute of Structural Analysis, Poznan University of Technology, Piotrowo 5 Street, 60-965 Poznan, Poland
3
Department of Mathematics, College of Sciences, King Khalid University, Abha 61413, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(3), 99; https://doi.org/10.3390/fractalfract5030099
Submission received: 2 July 2021 / Revised: 11 August 2021 / Accepted: 14 August 2021 / Published: 17 August 2021

Abstract

:
The present paper deals with the advancement of non-Newtonian fluid containing some nanoparticles between two parallel plates. A novel fractional operator is used to model memory effects, and analytical solutions are obtained for temperature and velocity fields by the method of Laplace transform. Moreover, a parametric study is elaborated to show the impact of flow parameters and presented in graphical form. As a result, dual solutions are predicted for increasing values of fractional parameters for short and long times. Furthermore, by increasing nanoparticle concentration, the temperature can be raised along with decreasing velocity. A fractional approach can provide new insight for the analytical solutions which makes the interpretation of the results easier and enable the way of testing possible approximate solutions.

1. Introduction

Channel flow plays an important part in an inclusive range of industrial applications, such as heat exchangers in power plants and chemical reactors in the development industry. Although many processes of practical interest can be considered as two-phase flows by Newtonian behavior in both phases, there is a large number of related applications where the continuous liquid phase exhibits non-Newtonian flow features. Many examples can be found in the biochemical and biomedical as well as the food processing industries [1]. Zheng et al. [2] investigated the impact of a vortex generator shape on liquids and the heat transformation of hybrid nanofluids in a channel. D’Ippolito et al. [3] studied the resistance of open channel flow because of vegetation at the research scale.
Nanofluids (term familiarized by Choi [4]) are considered identical in significant studies due to their thermal transfer development application. Some recent studies showed that the capability of heat transfer of nanofluids is more improved than of consistent fluids. Therefore, it is reasonable to exchange regular fluids with nanofluids. Because of their higher heat capability, the attention of numerous researchers is drawn toward the study of nanofluids. They play a dynamic role in industrial sectors such as electronics, transportation, foods, nuclear reactor, as well as biomedicine. The size of nanoparticles is very small, and these particles enhance the conductivity of base fluids when added to them. The construction of nanoparticles consists of metal oxide, carbon tubes, nitride, metal, and carbide [5]. Reddy et al. [6] addressed the performance of a micropolar thermal flux model magnesium oxide and molybdenum disulfide nanoparticle. The impact of a nanoparticle’s mixture ratio, as well as the temperature and concentration on the thermal conductivity of hybrid nanofluids, was established by Wole-Osho et al. [7]. Haque et al. [8] discussed the laminar nanoparticles having heat transfer, using convection force in the presence of noncircular ducts. Khashi’ie et al. [9] deliberated a stream of hybrid nanoparticles due to melting heat transfer in the presence of a moving surface. Goldanlou et al. [10] introduced a heat exchanger pipe prepared by the edge-form of turbulators of the heat transfer of hybrid nanoparticles in an explosive. Montazer et al. [11] described the development of the original relationships and the heat transfer of the stormy nanoparticles. Benkhedda et al. [12] defined the horizontal tube considering the forms of nanoparticles, resulting in convective heat transfer presentation of the hybrid nanoparticles. Bhattad et al. [13] conferred the plate in the heat exchanger of hydrothermal presentation of dissimilar alumina hybrid nanoparticles. Izadi et al. [14] exclusively debated the three-sided cavity and the properties of porous solid on fleeting normal convection of the heat transfer of nanofluids.
Concerning the hybrid nanofluid, which is a combination of two diverse nanosized particles with base liquids, which influences the convective heat transfer, one can follow the following progress in the literature. Waini et al. [15] examined steady fully developed mixed convection flow of a lengthwise an upright surface, fixed in a porous medium with hybrid nanoparticles was examined. Pandya et al.’s [16] study focuses on the numerical model, developed to maximize the thermal performance of axial lined heat tube, working on CeO2+ MWCNT/water-based hybrid nano-liquid. Babazadeh et al. [17] investigated the hybrid nano-liquid of a free convection transport inside a permeable average in the appearance of an outside applied magnetics force, modeling it through control volume, creating a fixed element technique. In Asadi et al.’s study [18], the discussion about the influence of hybrid nano-liquid n on the system performance was presented. Ikram et al. [19] studied clay nanoparticles for diverse-based liquids with a modern fractional scientific model that extracts nanofluids. Huminic et al.’s study [20] addressed dissimilar kinds of thermal structures aimed at dissimilar boundary situations and physical conditions, and contained a review of the contributions on entropy generation of the nano-liquid and hybrid nano-liquids. In their study, Nadeem et al. [21] addressed the base fluid ending in an exponentially stretching curved surface. The investigation evaluates the stimulus of including constraints of the hybrid nano-element.
Several fractional operators like the Grünwald–Letnikov derivative, Liouville and Riemann–Liouville derivatives, Caputo derivative, Marchaud derivative, Hadamard derivative, and Riesz derivative were applied as well [22,23]. Sene [24] studied the design of second-grade liquids, using the fractional differential equation through the Caputo Liouville fractional derivative. Mozafarifarda et al. [25] analyzed the fractional heat transfer equation for thin metal films utilizing the Caputo fractional derivative. Thabet et al. [26] studied the numerical findings of a novel disease of COVID-19 by using the mathematical method related to the existence and stability of ABC.
In the recent literature, researchers have paid attention to the vertical geometry for fractional models with different operators but they ignored the fluid flow problem between two plates and obtained the fractional model through generalized constitutive laws with the Prabhakar fractional model. The main advantage of a fractional operator is only the kernel that consists of three parameters, but, namely, the generalized Mittag-Leffler function. The Prabhakar fractional operator with a fractional coefficient can be a useful tool for determining an appropriate mathematical model that can result in a good agreement between theoretical and real data. For a detailed discussion about the Prabhakar integral and its properties, see the following references [27,28,29,30]. Giusti and Colombaro [31] proposed a generalized fractional model for non-Newtonian fluid by using the generalized constitutes laws. The Prabhakar type fractional derivative contains a generalized Mittage-Leffler function as a kernel that is responsible for the history of fluid properties and controls the thermal and momentum boundary layers. Elnaqeeb et al. [32] studied an infinite vertical heated plate via Prabhakar-like thermal carriage. The generalized fractional constitutive equation aimed at shear stress proportion and thermal flux density vector of natural convective flows of Prabhakar-like fractional viscoelastic fluids via an infinite vertical frenzied wall was studied by Shah et al. [33]. Finally, in Akgül et al.’s study [34], the magnetohydrodynamics effects on heat transfer phenomena was examined. Wang et al. [35] discussed the Casson nanofluid with a modified Mittag-Leffler kernel of a Prabhakar form of unsteady thermal transportation flow.
In the presented state of the art, no work was applying a fractional operator in an important class of fluids, like non-Newtonian fluid between two vertical plates. The Prabhakar fractional derivative was used to model the channel flow problem, containing hybrid nanoparticles. Therefore, the prime interest is to find the analytical solutions of energy and momentum with the Laplace transform method and parametric analysis to show the flexibility of the proposed model.

2. Materials and Methods

Let us consider MHD convection flow through the microchannel of electrically conductive (Ag-TiO2-H2O) hybrid nanofluid as shown in the Figure 1 under the following constraints:
(a)
Microchannel length is infinite with width L;
(b)
The channel is along x-axis and normal to y-axis;
(c)
At t ≤ 0, the temperature of the system is T 0 ;
(d)
After t = 0 + , the temperature increases from T 0 to T w ;
(e)
Fluid accelerates in the x-direction;
(f)
The magnetic field of strength B 0 works transversely to the flow direction.
The flow of electrically conductive Ag-TiO2 -H2O hybrid nanofluid causes the electromotive force, which yields current. Simultaneously, the induced magnetic field is ignored because of the hypothesis of a very small Reynolds number. Moreover, the electromagnetic force changes the intensity of the electric flux [36].
Within the above discussion, the governing equations for the stated problem are [37].
  • The momentum equation
    ρ h n f u y , t t + β b * u y , t = μ h n f 2 u y , t y 2 σ h n f B 0 2 u y , t + g ρ β T h n f T T 0 ,
  • The energy balance equation
    ρ C p h n f T y , t t = q y , t y ,
  • The generalized Fourier’s law for thermal flux
    q y , t = k h n f C D α , β , α γ T y , t y .
We consider the following initial and boundary conditions:
u y , 0 = 0 ,              T y , 0 = T 0 ,              y 0 ,
u 0 , t = 0 ,              T 0 , t = T 0              t > 0 ,
u L , t = 0 ,              T L , t = T w ,              t > 0 ,
where the definition of the regularized Prabhakar derivative is [32,38]
C D α , β , α γ f ( t ) = E α , m β , α γ f ( m ) ( t ) = e α , m β γ ( a ; t ) f ( m ) ( t ) = 0 t ( t τ ) m β 1 E α , m β γ a ( t τ ) α f ( m ) ( τ ) d τ ,
where “*” represent the convolution product, f m denotes the mth derivative of f(t) ∈A C m (0, b) stands for the set of real-valued functions f(t) whose derivatives are continuous up to order (m − 1) on the interval (0, b) where f m 1 (t) is absolutely continuous function, and m = [β] represents the integer part of the parameter β.
For Equations (1)–(3), the thermophysical boundaries of the considered nanofluid and hybrid nanofluid are listed in Table 1.
Finally, for generality, we introduce the dimensionless form as
τ = v f L 2 t ,   V = L v f u , Y = y L , θ = T T 0 T w T 0 ,   q * = q q 0 , q 0 = k h n f T w T u 0 v h n f .
Therefore (neglecting the star notations), we have
A 0 * V Y , τ τ + β b * V Y , τ = A 1 * 2 V Y , τ Y 2 A 2 * M V Y , τ + A 3 * G r θ Y , τ ,
P r θ Y , τ τ = q Y , τ Y ,
q Y τ = C D α , β , α γ θ Y , τ Y
under the constraints
V Y , 0 = 0 ,              θ Y , 0 = 0 ,              Y 0 ,
V 0 , τ = 0 ,              θ 0 , τ = 0 ,              τ > 0 ,
V 1 , τ = 0 ,              θ 1 , τ = 1 ,              τ > 0 ,
where
β b * = β b * L 2 v f ,   M = L 2 σ f B 0 2 μ f ,   G r = g β T f T w T 0 L 3 v f 2 ,   P r = μ C p f k f ,   λ h n f = k h n f k f ,   A 1 * = 1 1 Φ A g + Φ T i O 2 2.5 , A 0 * = 1 Φ h n f + Φ A g ρ A g + Φ T i O 2 ρ T i O 2 ρ f ,   A 2 * = σ h n f σ f ,   A 3 * = 1 Φ h n f + Φ A g ρ β T A g + Φ T i O 2 ρ β T T i O 2 ρ β T f .
In the above equations, β b * denotes Brinkman parameter, M is the magnetic field, P r is the Prandtl number, and G r is the Grashof number, respectively. It is essential that to accomplish the heat memory impacts, we additionally postulate the summed-up Fourier’s law dependent on Prabhakar’s fractional derivative, namely.

3. Solution of the Problem

3.1. Temperature Field

Applying the Laplace transform to Equations (10) and (11) with requirements 13 2 and 14 2 and utilizing Prabhakar fractional derivative, we get for the temperature field
P r s θ ¯ Y , s = q ¯ Y , s y ,
q ¯ Y , s = s β 1 a s α γ θ ¯ Y , s Y .
In the following, using Equations (16) and (17), we have
2 θ ¯ Y , s Y 2   P r s θ ¯ Y , s s β 1 a s α γ = 0 ,
which satisfies the following boundary conditions
θ ¯ 0 , s = 0 ,                    θ ¯ 1 , s = 1 s ,                    τ > 0 .
The general solution of Equation (18) with conditions Equation (19) is
θ ¯ Y , s = 1 s s i n h Y P r s s β 1 a s α γ s i n h P r s s β 1 a s α γ .
It is important that Equation (20) can be written in the equivalent form
θ ¯ Y , s = 1 s n = 0 e 2 n + 1 Y P r s s β 1 a s α γ e 2 n + 1 + Y P r s s β 1 a s α γ .
Moreover, Equation (21) can also be expressed in series form so that we can find the inverse of Laplace transform analytically, namely,
θ ¯ Y , s = 1 s + [ n = 0 m = 1 k = 0 2 n 1 + Y m a k P r m 2 Γ γ m 2 + k m ! k !         Γ γ m 2   1 s α k + β m 2 m 2 + 1 ] [ n = 0 l = 0 p = 0 2 n 1 Y l a p P r l 2   Γ γ l 2 + p l ! p !         Γ γ l 2   1 s α p + β l 2 l 2 + 1 ]
with its inverse Laplace transform
θ Y , τ =   1 + n = 0 m = 1 k = 0 2 n 1 + Y m a k P r m 2   Γ γ m 2 + k m ! k !       Γ γ m 2   t α k + β m 2 m 2 Γ α k + β m 2 m 2 + 1 n = 0 l = 0 p = 0 2 n + 1 + Y l a p P r l 2   Γ γ l 2 + p   l ! p !       Γ γ l 2   t α p + β l 2 l 2 Γ α p + β l 2 l 2 + 1 .

3.2. Velocity Field

By applying the Laplace transform to Equation (9) with constraints , 12 1 , 13 1 and 14 1 , we obtain
2 Y 2 B 1 s B 2 V ¯ Y , s = B 3 θ ¯ Y , s ,
which satisfies the following conditions
V ¯ 0 , s = 0 ,                    V ¯ 1 , s = 0 ,                    τ > 0 .
The solution of Equation (24) subjected to constraint (25) is
V ¯ Y , s = B 4 s n = 0 e 2 n P r s s β 1 a s α γ e 2 n + 2 P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ × m = 0 e 2 m + 1 Y B 2 + B 1 s m = 0 e 2 m + 1 + Y B 2 + B 1 s +                                                       B 4 s n = 0 e 2 n + 1 Y P r s s β 1 a s α γ e     2 n + 1 + Y P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ ,
where
B 1 = A 0 * A 1 * ,        B 2 = A 2 * M + β b * A 1 * ,        B 3 = A 3 * G r A 1 * ,        B 4 = B 3 B 2 .
Equation (26) can be written as
V ¯ Y , s = V 1 ¯ Y , s + V 2 ¯ Y , s + V 3 ¯ Y , s + V 4 ¯ Y , s + V 5 ¯ Y , s + V 6 ¯ Y , s ,
where
V 1 ¯ Y , s = B 4 s n = 0 e 2 n P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ × m = 0 e 2 m + 1 Y B 2 + B 1 s ,
V 2 ¯ Y , s = B 4 s n = 0 e 2 n P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ × m = 0 e 2 m + 1 + Y B 2 + B 1 s ,
V 3 ¯ Y , s = B 4 s n = 0 e 2 n + 2 P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ × m = 0 e 2 m + 1 Y B 2 + B 1 s ,
V 4 ¯ Y , s = B 4 s n = 0 e 2 n + 2 P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ × m = 0 e 2 m + 1 + Y B 2 + B 1 s ,
V 5 ¯ Y , s = B 4 s n = 0 e 2 n + 1 Y P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ ,
V 6 ¯ Y , s = B 4 s n = 0 e 2 n + 1 + Y P r s s β 1 a s α γ 1 + B 1 s B 2 P r s B 2 s β 1 a s α γ .
It is challenging to find the inverse of Laplace transform of Equations (28)–(33), so we can rewrite it in a suitable series form.
V 1 ¯ ( Y , s ) = B 4 n = 0 m = 0 δ 1 = 0 δ 2 = 0 ω 1 = 0 ω 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0   ( Y 2 m 1 ) ω 1 ( B 1 ) z 1 ( P r ) δ 1 2 δ 1 ! δ 2 ! ω 1 ! ω 2 ! z 2 ! z 3 ! ( B 1 ) z 2 ω 2   ( B 2 ) ω 2 + z 1 ω 1 2 × ( 2 n ) δ 1 ( P r ) z 2 ( a ) δ 2 + z 3 s α δ 2 + β δ 1 2 δ 1 2 ω 2 + β z 2 + α z 3 z 1 + 1     × Γ ( γ δ 1 2 + δ 2 ) Γ ( ω 1 2 + 1 ) Γ ( z 1 + 1 ) Γ ( γ z 2 + z 3 ) Γ ( γ δ 1 2 ) Γ ( ω 1 2 + 1 ω 2 ) Γ ( z 1 + 1 z 2 ) Γ ( γ z 2 ) ,
V 2 ¯ ( Y , s ) = B 4 n = 0 m = 0 δ 1 = 0 δ 2 = 0 q 1 = 0 q 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0   ( ( Y + 2 m + 1 ) ) q 1 ( B 1 ) z 1 ( P r ) δ 1 2 δ 1 ! δ 2 ! q 1 ! q 2 ! z 2 ! z 3 ! ( B 1 ) z 2 q 2   ( B 2 ) q 2 + z 1 q 1 2 × ( 2 n ) δ 1 ( P r ) z 2 ( a ) δ 2 + z 3 s α δ 2 + β δ 1 2 δ 1 2 q 2 + β z 2 + α z 3 z 1 + 1     × Γ ( γ δ 1 2 + δ 2 ) Γ ( q 1 2 + 1 ) Γ ( z 1 + 1 ) Γ ( γ z 2 + z 3 ) Γ ( γ δ 1 2 ) Γ ( q 1 2 + 1 q 2 ) Γ ( z 1 + 1 z 2 ) Γ ( γ z 2 ) ,
V 3 ¯ Y , s = B 4 n = 0 m = 0 j 1 = 0 j 2 = 0 ω 1 = 0 ω 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y 2 m 1 ω 1 B 1 z 1 P r j 1 2 j 1 ! j 2 ! ω 1 ! ω 2 ! z 2 ! z 3 ! B 1 z 2 ω 2   B 2 ω 2 + z 1 ω 1 2 ×     2 n 2 j 1 P r z 2 a j 2 + z 3 s α j 2 + β j 1 2 j 1 2 ω 2 + β z 2 + α z 3 z 1 + 1 × Γ γ j 1 2 + j 2 Γ ω 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ j 1 2 Γ ω 1 2 + 1 ω 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 4 ¯ Y , s = B 4 n = 0 m = 0 j 1 = 0 j 2 = 0 q 1 = 0 q 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y + 2 m + 1 q 1 B 1 z 1 P r j 1 2 j 1 ! j 2 ! q 1 ! q 2 ! z 2 ! z 3 ! B 1 z 2 q 2   B 2 q 2 + z 1 q 1 2 ×   2 n 2 j 1 P r z 2 a j 2 + z 3 s α j 2 + β j 1 2 j 1 2 q 2 + β z 2 + α z 3 z 1 + 1 × Γ γ j 1 2 + j 2 Γ q 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ j 1 2 Γ q 1 2 + 1 q 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 5 ¯ Y , s =   B 4 n = 0 t 1 = 0 t 2 = 0 z 1 = 0 z 2 = 0 z 3 Y 2 n 1 t 1 P r t 1 2 B 1 z 1 t 1 ! t 2 ! z 2 ! z 3 !     B 2 z 1 B 1 z 2 × P r z 2 a t 2 + z 3 s α t 2 + β t 1 2 t 1 2 + β z 2 + α z 3 z 1 + 1 × Γ γ t 1 2 + t 2 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ t 1 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 6 ¯ Y , s = B 4 n = 0 v 1 = 0 v 2 = 0 z 1 = 0 z 2 = 0 z 3 Y + 2 n + 1 v 1 P r v 1 2 B 1 z 1 v 1 ! v 2 ! z 2 ! z 3 !     B 2 z 1 B 1 z 2 × P r z 2 a v 2 + z 3 s α v 2 + β v 1 2 v 1 2 + β z 2 + α z 3 z 1 + 1   ×                                             Γ γ v 1 2 + v 2 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ v 1 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
Next, taking the inverse of Equations (34)–(39) Laplace transform, we have
V Y , τ = V 1 Y , τ + V 2 Y , τ + V 3 Y , τ + V 4 Y , τ + V 5 Y , τ + V 6 Y , τ ,
V 1 Y , τ = B 4 n = 0 m = 0 δ 1 = 0 δ 2 = 0 ω 1 = 0 ω 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y 2 m 1 ω 1 B 1 z 1 P r δ 1 2 δ 1 ! δ 2 ! ω 1 ! ω 2 ! z 2 ! z 3 ! B 1 z 2 ω 2   B 2 ω 2 + z 1 ω 1 2 ×                                       2 n δ 1 P r z 2 a δ 2 + z 3   τ α δ 2 + β δ 1 2 δ 1 2 ω 2 + β z 2 + α z 3 z 1 Γ α δ 2 + β δ 1 2 δ 1 2 ω 2 + β z 2 + α z 3 z 1 + 1 × Γ γ δ 1 2 + δ 2 Γ q 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ δ 1 2 Γ q 1 2 + 1 q 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 2 Y , τ = B 4 n = 0 m = 0 δ 1 = 0 δ 2 = 0 q 1 = 0 q 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y + 2 m + 1 q 1 B 1 z 1 P r δ 1 2 δ 1 ! δ 2 ! q 1 ! q 2 ! z 2 ! z 3 ! B 1 z 2 q 2   B 2 q 2 + z 1 q 1 2 ×                                       2 n δ 1 P r z 2 a δ 2 + z 3 τ α δ 2 + β δ 1 2 δ 1 2 q 2 + β z 2 + α z 3 z 1 Γ α δ 2 + β δ 1 2 δ 1 2 q 2 + β z 2 + α z 3 z 1 + 1 × Γ γ δ 1 2 + δ 2 Γ q 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ δ 1 2 Γ q 1 2 + 1 q 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 3 Y , τ = B 4 n = 0 m = 0 j 1 = 0 j 2 = 0 ω 1 = 0 ω 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y 2 m 1 ω 1 B 1 z 1 P r j 1 2 j 1 ! j 2 ! ω 1 ! ω 2 ! z 2 ! z 3 ! B 1 z 2 ω 2   B 2 ω 2 + z 1 ω 1 2 × 2 n 2 j 1 P r z 2 a j 2 + z 3 τ α j 2 + β j 1 2 j 1 2 q 2 + β z 2 + α z 3 z 1 Γ α j 2 + β j 1 2 j 1 2 q 2 + β z 2 + α z 3 z 1 + 1 × Γ γ δ 1 2 + δ 2 Γ q 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ δ 1 2 Γ q 1 2 + 1 q 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
V 4 Y , τ = B 4 n = 0 m = 0 j 1 = 0 j 2 = 0 q 1 = 0 q 2 = 0 z 1 = 0 z 2 = 0 z 3 = 0 Y + 2 m + 1 q 1 B 1 z 1 P r j 1 2 j 1 ! j ! q 1 ! q 2 ! z 2 ! z 3 ! B 1 z 2 q 2   B 2 q 2 + z 1 q 1 2 ×                                     2 n 2 j 1 P r z 2 a j 2 + z 3 τ α j 2 + β j 1 2 j 1 2 q 2 + β z 2 + α z 3 z 1 Γ α j 2 + β j 1 2 j 1 2 q 2 + β z 2 + α z 3 z 1 + 1 × Γ γ j 1 2 + j 2 Γ q 1 2 + 1 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ j 1 2 Γ q 1 2 + 1 q 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
    V 6 Y , τ = B 4 n = 0 v 1 = 0 v 2 = 0 z 1 = 0 z 2 = 0 z 3 Y + 2 n + 1 v 1 P r v 1 2 B 1 z 1 v 1 ! v 2 ! z 2 ! z 3 !     B 2 z 1 B 1 z 2 × P r z 2 a v 2 + z 3 τ α v 2 + β v 1 2 v 1 2 + β z 2 + α z 3 z 1 Γ α v 2 + β v 1 2 v 1 2 + β z 2 + α z 3 z 1 + 1 ×       Γ γ v 1 2 + v 2 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ v 1 2 Γ z 1 + 1 z 2 Γ γ z 2 ,
    V 6 Y , τ = B 4 n = 0 v 1 = 0 v 2 = 0 z 1 = 0 z 2 = 0 z 3 Y + 2 n + 1 v 1 P r v 1 2 B 1 z 1 v 1 ! v 2 ! z 2 ! z 3 !     B 2 z 1 B 1 z 2 × P r z 2 a v 2 + z 3 τ α v 2 + β v 1 2 v 1 2 + β z 2 + α z 3 z 1 Γ α v 2 + β v 1 2 v 1 2 + β z 2 + α z 3 z 1 + 1 ×       Γ γ v 1 2 + v 2 Γ z 1 + 1 Γ γ z 2 + z 3 Γ γ v 1 2 Γ z 1 + 1 z 2 Γ γ z 2 .

4. Discussion

Analytical solutions of heat transfer flow of a Brinkman fluid with a Prabhakar fractional derivative are obtained with the Laplace transform method. In order to see the physical insight of flow parameters, and especially fractional parameters, some graphs are presented in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, and Figure 13. Moreover, we have applied a time-fractional operator in a heat transfer problem. We plotted Figure 2 and Figure 3 for the temperature and velocity fluids for different values of fractional parameters by considering small and large values of the time. It is evident from the Figures that for a large time the fluid properties can be enhanced by considering greater values of fractional parameters. However, in Figure 4 and Figure 5, an opposite trend was observed for small values of the time. It is concluded that the fractional solutions show dual behavior for a different time. The fractional parameters can be used to control the thermal and momentum boundary layer, respectively.
The influence of ϕ h n f volume fraction of hybrid nanoparticles is depicted in Figure 6 and Figure 7. It is found that the temperature can be enhanced for greater values of ϕ h n f and velocity depicted declines, respectively. For example, from ϕ h n f = 0 , the base fluid is water and has greater velocity, but on increasing ϕ h n f = 0.01 ,   0.02 ,   0.04 , decay was observed in velocity. Physically, for large values of ϕ h n f , the viscous forces dominate and finally velocity decreases. As a result, water-based hybrid nanoparticles are denser than pure water.
Temperature decreases when increasing the value of Pr in Figure 8. Physically, the thermal conductivity increasing the values of Pr, making fluid thicker, causes the minimum thickness of the thermal boundary layer. The effect of the velocity profile on Pr is shown in Figure 9. Velocity reduces for bigger values of Pr. Actually, the fluids have bigger viscosity for a large Pr number, and minimum thermal conductivity reduces the velocity due to thicker fluid. In the following, the influence of time to temperature and velocity can be seen in Figure 10 and Figure 11. Time advance causes the increase of the temperature as well as velocity.
Velocity increases as we increase the values of Gr in Figure 12. This is because when Gr is increased, the buoyancy forces become stronger, which results in more convection taking place. Therefore, after that, the velocity increases. Figure 13 shows that the velocity decreases as we rise the values of M. Physically, it can occur that it answers to the drag force, which has an influence on the velocity profile that faces the fluid motion. As a result, velocity decreases.

5. Conclusions

This paper deals with the analytical solutions of fractional partial differential equations appearing in heat transport phenomena. The analysis includes base fluid modeled as water together with different nanoparticles. The solutions are obtained via the Laplace transform method, and the following conclusions can be drawn:
  • Dual solutions for both temperature and velocity fields are predicted for different values of fractional parameters describing generalized. Fourier’s law (based on Prabhakar’s fractional derivative);
  • Fluid property, i.e., temperature, can be enhanced by increasing the concentration level of nanoparticles;
  • Finally, the fractional parameters (which accomplish the heat memory impacts) can control the thermal and momentum boundary layer thickness;
  • The obtained solutions can be beneficial for proper analysis of real data and provide a tool for testing possible approximate solutions where needed.

Author Contributions

M.I.A. Formal analysis, Conceptualization, N.S. Visualization, Validation. M.B.H. original draft: Software, Investigation, W.S. Funding acquisition, Supervision, Project administration, T.M. Resources, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by National Science Center grant number 2017/27/B/ST8/0351.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available within the article.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, Abha, Saudi Arabia, for funding this work through the general research groups program under grant number GRP/342/42.

Conflicts of Interest

There is no conflict of interest.

References

  1. Bräuer, F.; Trautner, E.; Hasslberger, J.; Cifani, P.; Klein, M. Turbulent Bubble-Laden Channel Flow of Power-Law Fluids: A Direct Numerical Simulation Study. Fluids 2021, 6, 40. [Google Scholar] [CrossRef]
  2. Zheng, Y.; Yang, H.; Mazaheri, H.; Aghaei, A.; Mokhtari, N.; Afrand, M. An investigation on the influence of the shape of the vortex generator on fluid flow and turbulent heat transfer of hybrid nanofluid in a channel. J. Therm. Anal. Calorim. 2021, 143, 1425–1438. [Google Scholar] [CrossRef]
  3. D’Ippolito, A.; Calomino, F.; Alfonsi, G.; Lauria, A. Flow Resistance in Open Channel Due to Vegetation at Reach Scale: A Review. Water 2021, 13, 116. [Google Scholar] [CrossRef]
  4. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles (No. ANL/MSD/CP-84938.CONF-951135-29); Argonne National Lab.: Argonne, IL, USA, 1995. [Google Scholar]
  5. Nadeem, S.; Ahmad, S.; Khan, M.N. Mixed convection flow of hybrid nanoparticle along a Riga surface with Thomson and Troian slip condition. J. Therm. Anal. Calorim. 2021, 143, 2099–2109. [Google Scholar] [CrossRef]
  6. Reddy, M.G.; Shehzad, S.A. Molybdenum disulfide and magnesium oxide nanoparticle performance on micropolar Cattaneo-Christov heat flux model. Appl. Math. Mech. 2021, 42, 541–552. [Google Scholar] [CrossRef]
  7. Wole-Osho, I.; Okonkwo, E.C.; Adun, H.; Kavaz, D.; Abbasoglu, S. An intelligent approach to predicting the effect of nanoparticle mixture ratio, concentration and temperature on thermal conductivity of hybrid nanofluids. J. Therm. Anal. Calorim. 2021, 144, 671–688. [Google Scholar] [CrossRef] [Green Version]
  8. Haque, M.E.; Hossain, M.S.; Ali, H.M. Laminar forced convection heat transfer of nanofluids inside noncircular ducts: A review. Powder Technol. 2020. [Google Scholar] [CrossRef]
  9. Khashi’ie, N.S.; Arifin, N.M.; Pop, I.; Nazar, R. Melting heat transfer in hybrid nanofluid flow along a moving surface. J. Therm. Anal. Calorim. 2020. [Google Scholar] [CrossRef]
  10. Goldanlou, A.S.; Sepehrirad, M.; Papi, M.; Hussein, A.K.; Afrand, M.; Rostami, S. Heat transfer of hybrid nanofluid in a shell and tube heat exchanger equipped with blade-shape turbulators. J. Therm. Anal. Calorim. 2021, 143, 1689–1700. [Google Scholar] [CrossRef]
  11. Montazer, E.; Shafii, M.B.; Salami, E.; Muhamad, M.R.; Yarmand, H.; Gharehkhani, S.; Badarudin, A. Heat transfer in turbulent nanofluids: Separation flow studies and development of novel correlations. Adv. Powder Technol. 2020, 31, 3120–3133. [Google Scholar] [CrossRef]
  12. Benkhedda, M.; Boufendi, T.; Tayebi, T.; Chamkha, A.J. Convective heat transfer performance of hybrid nanofluid in a horizontal pipe considering nanoparticles shapes effect. J. Therm. Anal. Calorim. 2020, 140, 411–425. [Google Scholar] [CrossRef]
  13. Bhattad, A.; Sarkar, J.; Ghosh, P. Hydrothermal performance of different alumina hybrid nanofluid types in plate heat exchanger. J. Therm. Anal. Calorim. 2020, 139, 3777–3787. [Google Scholar] [CrossRef]
  14. Izadi, M. Effects of porous material on transient natural convection heat transfer of nanofluids inside a triangular chamber. Chin. J. Chem. Eng. 2020, 28, 1203–1213. [Google Scholar] [CrossRef]
  15. Waini, I.; Ishak, A.; Groşan, T.; Pop, I. Mixed convection of a hybrid nanofluid flow along a vertical surface embedded in a porous medium. Int. Commun. Heat Mass Transf. 2020, 114, 104565. [Google Scholar] [CrossRef]
  16. Pandya, N.S.; Desai, A.N.; Tiwari, A.K.; Said, Z. Influence of the geometrical parameters and particle concentration levels of hybrid nanofluid on the thermal performance of axial grooved heat pipe. Therm. Sci. Eng. Prog. 2021, 21, 100762. [Google Scholar] [CrossRef]
  17. Babazadeh, H.; Shah, Z.; Ullah, I.; Kumam, P.; Shafee, A. Analysis of hybrid nanofluid behavior within a porous cavity including Lorentz forces and radiation impacts. J. Therm. Anal. Calorim. 2021, 143, 1129–1137. [Google Scholar] [CrossRef]
  18. Asadi, A.; Alarifi, I.M.; Foong, L.K. An experimental study on characterization, stability and dynamic viscosity of CuO-TiO2/water hybrid nanofluid. J. Mol. Liq. 2020, 307, 112987. [Google Scholar] [CrossRef]
  19. Ikram, M.D.; Asjad, M.; Ahmadian, A.; Ferrara, M. A new fractional mathematical model of extraction nanofluids using clay nanoparticles for different based fluids. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  20. Huminic, G.; Huminic, A. Entropy generation of nanofluid and hybrid nanofluid flow in thermal systems: A review. J. Mol. Liq. 2020, 302, 112533. [Google Scholar] [CrossRef]
  21. Nadeem, S.; Abbas, N.; Malik, M.Y. Inspection of hybrid based nanofluid flow over a curved surface. Comput. Methods Programs Biomed. 2020, 189, 105193. [Google Scholar] [CrossRef]
  22. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods; World Scientific: London, UK, 2012; Volume 3. [Google Scholar]
  23. Ortigueira, M.D.; Machado, J.T. Fractional derivatives: The perspective of system theory. Mathematics 2019, 7, 150. [Google Scholar] [CrossRef] [Green Version]
  24. Sene, N. Second-grade fluid model with Caputo–Liouville generalised fractional derivative. Chaos Solitons Fractals 2020, 133, 109631. [Google Scholar] [CrossRef]
  25. Mozafarifard, M.; Toghraie, D. Time-fractional subdiffusion model for thin metal films under femtosecond laser pulses based on Caputo fractional derivative to examine anomalous diffusion process. Int. J. Heat Mass Transf. 2020, 153, 119592. [Google Scholar] [CrossRef]
  26. Thabet, S.T.; Abdo, M.S.; Shah, K.; Abdeljawad, T. Study of transmission dynamics of COVID-19 mathematical model under ABC fractional order derivative. Results Phys. 2020, 19, 103507. [Google Scholar] [CrossRef]
  27. Prabhakar, T.R. A Singular Integral Equation with a Generalized Mittag Leffler Function in the Kernel. 1971. Available online: https://www.semanticscholar.org/paper/A-SINGULAR-INTEGRAL-EQUATION-WITH-A-GENERALIZED-IN-Prabhakar/cf26590a6a758228228045275591231430db614f (accessed on 2 July 2021).
  28. Yang, X.J. General Fractional Derivatives: Theory, Methods and Applications; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  29. Garra, R.; Gorenflo, R.; Polito, F.; Tomovski, Ž. Hilfer–Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef] [Green Version]
  30. Fernandez, A.; Baleanu, D. Classes of operators in fractional calculus: A case study. Math. Methods Appl. Sci. 2021. [Google Scholar] [CrossRef]
  31. Giusti, A.; Colombaro, I. Prabhakar-like fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 138–143. [Google Scholar] [CrossRef] [Green Version]
  32. Elnaqeeb, T.; Shah, N.A.; Mirza, I.A. Natural convection flows of carbon nanotubes nanofluids with Prabhakar-like thermal transport. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  33. Shah, N.A.; Fetecau, C.; Vieru, D. Natural convection flows of Prabhakar-like fractional Maxwell fluids with generalized thermal transport. J. Therm. Anal. Calorim. 2021, 143, 2245–2258. [Google Scholar] [CrossRef]
  34. Akgül, A.; Siddique, I. Novel applications of the magnetohydrodynamics couple stress fluid flows between two plates with fractal-fractional derivatives. Numer. Methods Part. Differ. Equ. 2021, 37, 2178–2189. [Google Scholar] [CrossRef]
  35. Wang, F.; Asjad, M.I.; Zahid, M.; Iqbal, A.; Ahmad, H.; Alsulami, M.D. Unsteady thermal transport flow of Casson nanofluids with generalized Mittag-Leffler kernel of Prabhakar’s type. J. Mater. Res. Technol. 2021. [Google Scholar] [CrossRef]
  36. Tashtoush, B.; Magableh, A. Magnetic field effect on heat transfer and fluid flow characteristics of blood flow in multi-stenosis arteries. Heat Mass Transf. 2008, 44, 297–304. [Google Scholar] [CrossRef]
  37. Saqib, M.; Shafie, S.; Khan, I.; Chu, Y.M.; Nisar, K.S. Symmetric MHD channel flow of nonlocal fractional model of BTF containing hybrid nanoparticles. Symmetry 2020, 12, 663. [Google Scholar] [CrossRef]
  38. Polito, F.; Tomovski, Z. Some properties of Prabhakar-type fractional calculus operators. Fract. Differ. Calc. 2016, 6, 73–94. [Google Scholar] [CrossRef]
Figure 1. Flow geometry.
Figure 1. Flow geometry.
Fractalfract 05 00099 g001
Figure 2. Temperature distribution against y due to fractional parameters α , β , γ for large time, when t = 3, Pr = 6.2,   ϕ h n f = 0.02.
Figure 2. Temperature distribution against y due to fractional parameters α , β , γ for large time, when t = 3, Pr = 6.2,   ϕ h n f = 0.02.
Fractalfract 05 00099 g002
Figure 3. Velocity distribution against y due to for fractional parameters α , β , γ for large time, when t = 3, Pr = 6.2 ϕ h n f = 0.04, Gr = 15, M = 0.1.
Figure 3. Velocity distribution against y due to for fractional parameters α , β , γ for large time, when t = 3, Pr = 6.2 ϕ h n f = 0.04, Gr = 15, M = 0.1.
Fractalfract 05 00099 g003
Figure 4. Temperature distribution against y due to fractional parameters α , β , γ for small time, when t = 0.1, Pr = 6.2, ϕ h n f = 0.02.
Figure 4. Temperature distribution against y due to fractional parameters α , β , γ for small time, when t = 0.1, Pr = 6.2, ϕ h n f = 0.02.
Fractalfract 05 00099 g004
Figure 5. Velocity distribution against y due to fractional parameters α , β , γ for small time, when t = 0.3, Pr = 6.2 ϕ h n f = 0.04, Gr = 25, M = 0.1.
Figure 5. Velocity distribution against y due to fractional parameters α , β , γ for small time, when t = 0.3, Pr = 6.2 ϕ h n f = 0.04, Gr = 25, M = 0.1.
Fractalfract 05 00099 g005
Figure 6. Temperature distribution against y due to ϕ h n f , when t = 3, Pr = 6.2 α = β = γ = 0.5 .
Figure 6. Temperature distribution against y due to ϕ h n f , when t = 3, Pr = 6.2 α = β = γ = 0.5 .
Fractalfract 05 00099 g006
Figure 7. Velocity distribution against y due to ϕ h n f , when t = 40, Pr = 6.2, Gr = 15, M = 2, α = β = γ = 0.5 .
Figure 7. Velocity distribution against y due to ϕ h n f , when t = 40, Pr = 6.2, Gr = 15, M = 2, α = β = γ = 0.5 .
Fractalfract 05 00099 g007
Figure 8. Temperature distribution against y due to Pr, when t = 3 ϕ h n f = 0.02, α = β = γ = 0.5 .
Figure 8. Temperature distribution against y due to Pr, when t = 3 ϕ h n f = 0.02, α = β = γ = 0.5 .
Fractalfract 05 00099 g008
Figure 9. Velocity distribution against y due to Pr, when t = 3 ϕ h n f = 0.04, Gr = 15, M = 0.1, α = β = γ = 0.5 .
Figure 9. Velocity distribution against y due to Pr, when t = 3 ϕ h n f = 0.04, Gr = 15, M = 0.1, α = β = γ = 0.5 .
Fractalfract 05 00099 g009
Figure 10. Temperature distribution against y due to t, when Pr = 6.2, ϕ h n f = 0.02, α = β = γ = 0.5 .
Figure 10. Temperature distribution against y due to t, when Pr = 6.2, ϕ h n f = 0.02, α = β = γ = 0.5 .
Fractalfract 05 00099 g010
Figure 11. Velocity distribution against y due to t, when Pr = 6.2, ϕ h n f = 0.04, Gr = 15, M = 0.1, α = β = γ = 0.5 .
Figure 11. Velocity distribution against y due to t, when Pr = 6.2, ϕ h n f = 0.04, Gr = 15, M = 0.1, α = β = γ = 0.5 .
Fractalfract 05 00099 g011
Figure 12. Velocity distribution against y due to Gr, when Pr = 6.2, ϕ h n f = 0.04, t = 3, M = 0.1, α = β = γ = 0.5 .
Figure 12. Velocity distribution against y due to Gr, when Pr = 6.2, ϕ h n f = 0.04, t = 3, M = 0.1, α = β = γ = 0.5 .
Fractalfract 05 00099 g012
Figure 13. Velocity distribution against y due to M, when Pr = 6.2, ϕ h n f = 0.04, t = 3, Gr = 15, α = β = γ = 0.5 .
Figure 13. Velocity distribution against y due to M, when Pr = 6.2, ϕ h n f = 0.04, t = 3, Gr = 15, α = β = γ = 0.5 .
Fractalfract 05 00099 g013
Table 1. The thermophysical boundaries of the considered nanofluids.
Table 1. The thermophysical boundaries of the considered nanofluids.
NanofluidHybrid Nanofluid
ρ n f = 1 ϕ ρ f + ϕ ρ s ρ h n f = ( 1 ϕ h n f ) ρ f + ϕ A g ρ A g + ϕ T i O 2 ρ T i O 2
μ n f = μ f 1 ϕ 2.5 μ h n f = μ f 1 ϕ A g + ϕ T i O 2 2.5
ρ C p n f = 1 ϕ ρ C p f + ϕ ρ C p s ( ρ C p ) h n f = ( 1 ϕ h n f ) ( ρ C p ) f + ϕ A g ( ρ c p ) A g + ϕ T i O 2 ρ C p T i O 2
ρ B T n f = 1 ϕ ρ B T f + ϕ ρ B T s ( ρ β T ) h n f = ( 1 ϕ h n f ) ( ρ β T ) f + ϕ A g ( ρ β T ) A g + ϕ T i O 2 ( ρ β T ) T i O 2
σ n f σ f = 1 + 3 σ s σ f 1 ϕ σ s σ f + 2 σ s σ f 1 ϕ σ h n f σ f = 1 + 3 ϕ A g σ A g + ϕ T i O 2 σ T i O 2 σ f ϕ h n f ϕ A g σ A g + ϕ T i O 2 σ T i O 2 ϕ h n f σ f + 2 ϕ A g σ A g + ϕ T i O 2 σ T i O 2 σ f ϕ h n f
K n f K f = k s + 2 k f 2 ϕ k s k f k s + 2 k f + ϕ k s k f K h n f K f = ϕ A g k A g + ϕ T i O 2 k T i O 2 ϕ h n f + 2 k f + 2 ϕ A g k A g + ϕ T i O 2 k T i O 2 2 ϕ h n f K f ϕ A g k A g + ϕ T i O 2 k T i O 2 ϕ h n f + 2 k f + ϕ A g k A g + ϕ T i O 2 k T i O 2 ϕ h n f K f
where ρ h n f , u(y, t), β b * , μ h n f , σ h n f , B 0 , β h n f , C p h n f , T(y, t), and k h n f are the density, velocity, Brinkmann parameter, dynamic viscosity, electrical conductivity, magnetic field, thermal expansion, specific heat, temperature, and thermal conductivity, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Asjad, M.I.; Sarwar, N.; Hafeez, M.B.; Sumelka, W.; Muhammad, T. Advancement of Non-Newtonian Fluid with Hybrid Nanoparticles in a Convective Channel and Prabhakar’s Fractional Derivative—Analytical Solution. Fractal Fract. 2021, 5, 99. https://doi.org/10.3390/fractalfract5030099

AMA Style

Asjad MI, Sarwar N, Hafeez MB, Sumelka W, Muhammad T. Advancement of Non-Newtonian Fluid with Hybrid Nanoparticles in a Convective Channel and Prabhakar’s Fractional Derivative—Analytical Solution. Fractal and Fractional. 2021; 5(3):99. https://doi.org/10.3390/fractalfract5030099

Chicago/Turabian Style

Asjad, Muhammad Imran, Noman Sarwar, Muhammad Bilal Hafeez, Wojciech Sumelka, and Taseer Muhammad. 2021. "Advancement of Non-Newtonian Fluid with Hybrid Nanoparticles in a Convective Channel and Prabhakar’s Fractional Derivative—Analytical Solution" Fractal and Fractional 5, no. 3: 99. https://doi.org/10.3390/fractalfract5030099

APA Style

Asjad, M. I., Sarwar, N., Hafeez, M. B., Sumelka, W., & Muhammad, T. (2021). Advancement of Non-Newtonian Fluid with Hybrid Nanoparticles in a Convective Channel and Prabhakar’s Fractional Derivative—Analytical Solution. Fractal and Fractional, 5(3), 99. https://doi.org/10.3390/fractalfract5030099

Article Metrics

Back to TopTop