Next Article in Journal
Artificial Intelligence Methodologies for Data Management
Next Article in Special Issue
Solvability Analysis of a Mixed Boundary Value Problem for Stationary Magnetohydrodynamic Equations of a Viscous Incompressible Fluid
Previous Article in Journal
Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law

by
Noman Sarwar
1,
Muhammad Imran Asjad
1,
Thanin Sitthiwirattham
2,*,
Nichaphat Patanarapeelert
3 and
Taseer Muhammad
4
1
Department of Mathematics, University of Management and Technology, Lahore 54770, Pakistan
2
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
3
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
4
Department of Mathematics, College of Sciences, King Khalid University, Abha 61413, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(11), 2039; https://doi.org/10.3390/sym13112039
Submission received: 7 September 2021 / Revised: 13 October 2021 / Accepted: 22 October 2021 / Published: 29 October 2021
(This article belongs to the Special Issue Applied Mathematics and Fluid Dynamics)

Abstract

:
In the present work, an unsteady convection flow of Casson fluid, together with an oscillating vertical plate, is examined. The governing PDEs corresponding to velocity and temperature profile are transformed into linear ODEs with the help of the Laplace transform method. The ordinary derivative model generalized to fractional model is based on a generalized Fourier law. The solutions for energy and velocity equations are obtained after making the equations dimensionless. To check the insight of the physical parameters, especially the symmetric behavior of fractional parameters, it is found that for small and large values of time, fluid properties show dual behavior. Since the fractional derivative exhibits the memory of the function at the chosen value of time, therefore the present fractional model is more suitable in exhibiting memory than the classical model. Such results can be useful in the fitting of real data where needed. In the limiting case when fractional parameters are taken β = γ = 0 and α = 1 for both velocity and temperature, we get the solutions obtained with ordinary derivatives from the existing literature.

1. Introduction

Casson proposed the Casson fluid model in 1959 to predict flow patterns of suspended pigment–oil mixtures. These fluids are known as viscoelastic fluids, and they have a broad range of uses in the drilling process. Honey jelly soup, concentrated fruit juices and artificial fibers are all forms of Casson fluid [1]. The Casson fluid model has many uses in biomechanics and processes. In comparison to previous viscoelastic models, the Casson fluid model is more suitable for studying materials such as blood, chocolate, and other rheological details. At zero shear rate, Casson fluid is a shear-thinning liquid with infinite viscosity. The fluid behaves like a solid when the yield stress is greater than the shear stress [2,3]. Based on generalized Fourier’s and Fick’s laws, as well as mass and heat transfer, new fractional derivative model Casson fluid has been developed by Sheikh et al. [4]. A new method for solving the fractional generalized Casson fluids model of the Caputo fractional operator has been developed by Ndolane [5]. Arif et al. [6] deliberated applications of motor oil using a fractional model of Casson fluids of ramped wall temperature. Ali et al. [7] described Wright’s function solutions for time fractional Casson fluid flow in free convection. Raza et al. [8] conducted thermal movement of conservative Casson nanoparticles through ramped temperature because of an infinite perpendicular plate by fractional derivative. Khan et al. [9] examined Casson fluid fractional derivative generalized unsteady flow because of the infinite plate. Saqib et al. [10] expressed nonlocal fractional model of BTF incorporating hybrid nanoparticles with symmetric MHD channel flow.
Oscillating flow is a well-known reason for higher mass and heat transfer. Oscillating movement was realized moreover by fluids vibration around a fixed article. While fluids vibration around a fixed object needs more energy. Oscillating flows are widely used in piston engines, pulsating burners, heat exchangers, cryogenic refrigeration, chemical reactors, and in numerous applications in the atmosphere military and industry fields [11,12]. Song et al. [13] discussed the significance of exponentially Darcy’s and heating law for second grade fluid flow with oscillating plate by utilizing the Atangana–Baleanu time fractional operator. Ali et al. [14] analyzed the flow of second grade fluids on an oscillating perpendicular plate for unsteady free convection closed arrangement solutions. Haider et al. [15] deliberated heat transfer of micro polar hybrid nanofluids through an oscillating vertical plate of Newtonian heating with thermal radiation. Shah et al. [16] observed thermal analysis of second grade fluids over oscillation perpendicular flat plate by consuming time fractional Caputo–Fabrizio derivatives. The impacts of heat generation and thermal diffusion on magnetohydrodynamics (MHD) nanofluids flow through an oscillating vertical plate with permeable medium, as addressed by Sheikholeslami et al. [17].
It is observed that the numerical model of integer-order derivative, counting the non-linear model, does not work satisfactorily in numerous cases. Fractional calculus has many applications in the field of viscoelasticity, signals processing, electromagnetics, fluid mechanics, and optics. It has been utilized to illustrate the model physical and designing forms that are found by fractional differential conditions. Fang et al. [18] discussed Caputo fractional derivatives along with applications of time-fractional distribution equations through a fast technique for variable-order. Dassios and Baleanu [19] introduced Caputo’s fractional-differential equation for singular linear structure with an optimal solution. Ali et al. [20] described hybrid nanofluids along with constant comparative Caputo fractional derivatives of power laws memory because of a pressure gradient. Baleanu et al. [21] studied mathematical displaying of the human liver along with Caputo–Fabrizio fractional derivatives. Alizadeh et al. [22] analyzed the transitory reaction of the equivalent RCL circuit by utilizing Caputo–Fabrizio fractional derivatives. A non-linear fractional system model of dengue illness under Caputo–Fabrizio derivatives was discussed by Shah et al. [23]. Abdeljawad et al. [24] deliberated ABC-fractional analysis of approximately generalized logistic models. Thabet et al. [25] investigated transference dynamics of COVID-19 numerical models under ABC-fractional derivatives. Atangana–Baleanu fractional derivatives are a new application in designing of ABC-fractional mask image processing conversed by Ghanbari and Atangana [26]. Benchohra et al. [27] analyzed the Hilfer–Katugampola fractional derivative terminal value problem for differential equations. Using the Caputo–Fabrizio derivative, the unsteady flow of fractional fluid between two parallel walls with arbitrary wall shear stress was discussed by Saqib et al. [28].
The Prabhakar work is named after Indian mathematician Tilak Raj Prabhakar, who proposed a generalization of the Mittag–Leffler function to three parameters in 1971. The utilization of Prabhakar administrators with specific fractional coefficients could be a helpful path for a decision of satisfactory numerical models that settled on a decent arrangement among hypothetical and trial results [29,30]. Elnaqeeb et al. [31] debated the carbon nanotube nanofluid with natural convection flow and Prabhakar-like thermal transportation. Shah et al. [32] deliberated generalized thermal transport and natural convection streams of Prabhakar-like fractional of Maxwell fluids. Dynamics and stability of Prabhakar fractional differential structures that are integrodifferential and neutral were regularized by Eshaghi et al. [33]. Tanveer et al. [34] designated the Prabhakar fractional model used to study the thermal behavior of viscous carbon nanotube nanofluid in free convection flows with generalized thermal transport. Alidousti [35] investigated the fractional differential classifications through the Prabhakar derivative with a stable field. Derakhshan [36] defined, in the context of Hilfer–Prabhakar derivatives, a new numerical procedure for solving variable order fractional integrodifferential equations. Prabhakar fractional derivatives of non-Newtonian fluid in convection channel with hybrid nanoparticles were discussed by Asjad et al. [37]. Basit et al. [38] introduced various nanoparticles of convection flow of fractional second grade fluid subjects to the non-uniform velocity at the edge. Unstable Casson nanofluid of thermal transport flow generalized with Prabhakar types Mittage–Leffler kernel was analyzed by Wang et al. [39]. Sene and Srivastava [40] discussed the generalized Mittag–Leffler input stability of the fractional differential equations. Modeling, analysis, and solutions of time-fractional models for COVID-19 were expressed by Iyiola et al. [41].
Previously, Khalid et al. [42] described their model for unsteady Casson fluid across the oscillating vertical plate with constant temperature and found their results for small-time values. The present investigation is focused on the study of the Prabhakar fractional model of Casson fluid through generalized Fourier law for oscillating surface that controls the momentum and thermal boundary surface with Laplace transform methods. In this paper, the governing equations for non-Newtonian fluid are generalized to Casson fluid, and the governing equations are found by constitutive laws after making them dimensionless. Furthermore, analytical solutions for the fluid properties temperature and velocity are found using the Laplace transform approach. For insight on the behavior of flow particles, some graphs are plotted in the graphical section.

2. Preliminaries of Fractional Calculus

Definition 1
(The Prabhakar kernel). Function [43]. For t ∈ R,
e α , β γ ( a ; t ) = t β 1 E α , β γ ( a t α ) , α , β , γ , a C , Re ( α ) > 0 ,
is known as the Prabhakar kernel.
Definition 2
(The Prabhakar integral). Function [43]. For t ∈ (0, b) and a function h L 1 (0, b)
E α , β , a γ h ( t ) = h ( t ) e α , β γ ( a ; t ) = 0 t h ( τ ) ( t τ ) β 1 E α , β γ ( a ( t τ ) α ) d τ ,
is called the Prabhakar integral. Laplace transform of Equation (2) is given:
L { E α , β , a γ h ( t ) } ( s ) = L { h ( t ) } L { e α , β γ ( a ; t ) } = L { h ( t ) } s a γ β ( s α a ) γ .
Definition 3
(The Prabhakar derivative). Let g L 1 (0, b) and τ (0, b), then for m = [ ν ] . The Prabhakar derivative is defined by [43].
D α , ν , β , a + γ g = d m d t m E α , n ν , β , a + γ g ( τ ) ,
where α , β , γ , ν are the elements of complex number C.
Definition 4
(The regularized Prabhakar derivative). Let 0 < β n , where n belongs to Z and h∈ A C n (0, b). The regularized Prabhakar derivative is defined as [43,44]
c D α , β , a γ h ( t ) = h ( n ) ( t ) E α , n β , a γ = h n ( t ) e α , n β γ ( a ; t ) = 0 t h ( n ) ( τ ) ( t τ ) n β 1 E α , n β γ ( a ( t τ ) α ) d τ ,
where h ( n ) represents the nth differential of h(t) and A C n (0, b) represents the set of functions h(t) with real values. Laplace transform of the Equation (5) is give by
L ( c D α , β , a γ h ( t ) ) = L { h ( n ) ( t ) e α , n β γ } = L { h ( n ) ( t ) } L { e α , n β γ ( a ; t ) } = s β n ( 1 a s α ) γ L { h ( n ) ( t ) } .

3. Problem Formulation

Consider the heat transform of Casson fluid flowing past, an infinite vertical smooth plate y > 0 , where y is the coordinate measured in the surface’s normal direction. Both the plate and the fluid are assumed to be at rest with a constant temperature T at the start of the simulation t = 0 . In its plane, the plate start to oscillate y = 0 at time t = 0 + , according to
V ˜ = U ˜ s i n s ( w t ) i ^ ; o r V ˜ = U ˜ H ( t ) c o s ( w t ) i ^ ; t > 0 ,
where the constant U ˜ is the amplitude of the plate oscillations, H ( t ) is the unit step function, i ^ is the unit vector in the vertical flow direction, and ω is the frequency of oscillation of the plate. At the same time, the plate temperature is raised to T w . Which is then maintained at the same level (Figure 1). We assume that the rheological equation of state for an isotropic flow of Casson fluid can be written as [42]
τ ˜ = μ ˜ γ * + τ ˜ 0 ,
or
τ ˜ i j = { 2 μ ˜ B + p ˜ y 2 π ˜ e ˜ i j π ˜ > π ˜ c 2 μ ˜ B + p ˜ y 2 π ˜ c e ˜ i j π ˜ < π ˜ c ,
where π ˜ = e ˜ i j e ˜ i j and e ˜ i j is ( i , j ) i h component of deformations rate. π ˜ is a product of component deformation rate, π ˜ c is the critical value of this product based on the non-Newtonian model, μ ˜ B is the plastic dynamic viscosity of the non-Newtonian fluid, and p ˜ y is yield stress of the fluid. We derive the following set of partial differential equations under these conditions, and the energy equations of the viscous dissipation factor is assumed to be ignored:
Momentum equation:
ρ ( u ˜ t ( y , t ) ) = μ ˜ β ( 1 + 1 γ 1 ) u ˜ y y ( y , t ) + g ( ρ β ) ( T ˜ ( y , t ) T ˜ ) .
Energy equation:
( ρ C p ) T ˜ t ( y , t ) = q ˜ y ( y , t ) .
Generalized Fourier’s law for thermal flux:
q ˜ ( y , t ) = k ˜ C D α , β , a γ T ˜ y ( y , t ) ,
where C D α , β , a γ is the regularized Prabhakar derivative defined [30,31].
Constraints:
0 > t : u ˜ = 0 , T ˜ T ˜ = 0 y > 0 ,
0 t : u ˜ = U ˜ s i n ( ω t ) o r u ˜ = U ˜ H ˜ ( t ) c o s ( ω t ) , T ˜ T ˜ w = 0 a t y = 0 ,
u ˜ 0 , T ˜ T ˜ a s y .
introducing dimensionless variables.
τ ˜ * = τ ˜ ρ u 2 , Y = U ˜ ν y , V ˜ = u ˜ U ˜ , t * = U ˜ 2 ν t , ω * = ω ν U ˜ 2 , θ ˜ = T ˜ T ˜ T ˜ w T ˜ , q ˜ * = q ˜ q ˜ 0 , q ˜ 0 = k ˜ ( T ˜ w T ˜ ) u ˜ 0 ν ,
into Equations (10)–(15) from Equation (16) and we obtain (* symbols have been removed for clarity).
Dimensionless form of momentum:
V ˜ t ( Y , t ) = ( 1 + 1 γ 1 ) V ˜ Y Y ( Y , t ) + Gr θ ˜ ( Y , t ) ,
Dimensionless form of energy:
Pr θ ˜ t ( Y , t ) = q ˜ Y ( Y , t ) ,
Generalized Fourier law for thermal flux in dimensionless form:
q ˜ ( Y , t ) = C D α , β , a γ θ ˜ Y ( Y , t ) ,
with associated constraints are
t < 0 : V ˜ = 0 , θ ˜ = 0 Y > 0 ,
V ˜ = s i n ( ω t ) o r t 0 : H ˜ ( t ) c o s ( ω t ) , θ ˜ = 1 a s Y = 0 ,
V ˜ 0 , θ ˜ 0 a t Y ,
where Pr = μ C p k ˜ dimensionless Prandtl number, Gr = ν g β ( T ˜ w T ˜ ) U ˜ 3 dimensionless Grashof number and γ 1 = μ ˜ β 2 π ˜ c p ˜ y Casson parameter.

4. Results of the Problem

In this section, the Prabhakar derivative is used to construct a fractional model of the physical system.

4.1. Temperature Field’s Outcome

By using the Laplace transform of Equations (18) and (19), we get
Pr s θ ˜ ¯ ( Y , s ) = q ˜ ¯ Y ( Y , s ) ,
and
q ˜ ¯ ( Y , s ) = s β ( 1 a s α ) γ θ ˜ ¯ Y .
By introducing Equation (24) into Equation (23), we get
θ ˜ ¯ Y Y ( Y , s ) θ ˜ ¯ ( Y , s ) Pr s s β ( 1 a s α ) γ = 0 ,
satisfy the following constraints:
θ ˜ ¯ ( , s ) = 0 , θ ˜ ¯ ( 0 , s ) = 1 s .
With the general solution of Equation (25) subject to the constraints (26), we get
θ ˜ ¯ ( Y , s ) = 1 s e Y P r s s β ( 1 a s α ) γ .
We cannot find the inverse Laplace transform directly from the table of Laplace transform. Therefore, Equation (27) is expressed in series expression as follows:
θ ˜ ¯ ( Y , s ) = 1 s + = 1 p = 0 ( Y ) ( Pr ) 2 ( a ) p ! p ! s β 2 2 + α p + 1 Γ ( γ 2 + p ) Γ ( γ 2 ) .
Taking the Laplace inverse transform of Equation (28), we have
θ ˜ ( Y , t ) = 1 + = 1 p = 0 ( Y ) ( Pr ) 2 ( a ) p t β 2 2 + α p ! p ! Γ ( β 2 2 + α p + 1 ) Γ ( γ 2 + p ) Γ ( γ 2 ) ,
which satisfy the boundary condition.

4.2. Velocity Field’s Outcome

By applying the Laplace transform of Equation (17), we have
D 2 s a 1 V ˜ ¯ ( Y , s ) = Gr a 1 θ ˜ ¯ ( Y , s ) ,
satisfy following constraints:
V ˜ ¯ ( , s ) = 0 , V ˜ ¯ ( 0 , s ) = ω s 2 + ω 2 .
From Equations (30) and (31), we acquire the solution of velocity:
V ˜ ¯ ( Y , s ) = ω s 2 + ω 2 e Y s a 1 Gr s 2 e Y s a 1 ( 1 a 1 Pr s β ( 1 a s α ) γ ) + Gr s 2 e Y Pr s s β ( 1 a s α ) γ ( 1 a 1 Pr s β ( 1 a s α ) γ ) ,
where a 1 = ( 1 + 1 γ 1 ) .
Equation (32) can be expressed as,
V ˜ ¯ ( Y , s ) = V ˜ ¯ 1 ( Y , s ) + V ˜ ¯ 2 ( Y , s ) + V ˜ ¯ 3 ( Y , s ) + V ˜ ¯ 4 ( Y , s ) ,
where V ˜ ¯ 1 ,   V ˜ ¯ 2 ,   V ˜ ¯ 3   and   V ˜ ¯ 4 are expressed in a suitable form.
V ˜ ¯ 1 ( Y , s ) = 1 2 ω k = 0 k 1 = 0 ( Y ) k k ! ( i ω ) k 1 ( a 1 ) k 2 s k 2 k 1 ,
V ˜ ¯ 2 ( Y , s ) = 1 2 ω k = 0 k 2 = 0 ( Y ) k k ! ( i ω ) k 2 ( a 1 ) k 2 s k 2 k 2 ,
V ˜ ¯ 3 ( Y , s ) = Gr k = 0 ρ 1 = 0 ρ 2 = 0 ( Y ) k ( a 1 ) ρ 1 k 2 ( Pr ) ρ 1 ( a ) ρ 2 k ! ρ 2 ! s β ρ 1 + α ρ 2 k 2 + 2 Γ ( γ ρ 1 + ρ 2 ) Γ ( γ ρ 1 ) ,
V ˜ ¯ 4 ( Y , s ) = Gr λ 1 = 0 λ 2 = 0 ρ 1 = 0 ρ 2 = 0 ( Y ) λ 1 ( a 1 ) ρ 1 ( Pr ) ρ 1 + λ 1 2 ( a ) λ 2 + ρ 2 λ 1 ! λ 2 ! ρ 2 ! s β λ 1 2 λ 1 2 + α λ 2 + β ρ 1 + α ρ 2 + 2 Γ ( γ λ 1 2 + λ 2 ) Γ ( γ λ 1 2 ) Γ ( γ ρ 1 + ρ 2 ) Γ ( γ ρ 1 ) .
Next, taking the inverse of Equations (33) and (37),
V ˜ ( Y , t ) = V ˜ 1 ( Y , t ) + V ˜ 2 ( Y , t ) + V ˜ 3 ( Y , t ) + V ˜ 4 ( Y , t ) ,
V ˜ 1 ( Y , t ) = 1 2 ω k = 0 k 1 = 0 ( Y ) k k ! ( i ω ) k 1 ( a 1 ) k 2 t k 2 k 1 1 Γ ( k 2 k 1 ) ,
V ˜ 2 ( Y , t ) = 1 2 ω k = 0 k 2 = 0 ( Y ) k k ! ( i ω ) k 2 ( a 1 ) k 2 t k 2 k 2 1 Γ ( k 2 k 2 ) ,
V ˜ 3 ( Y , t ) = Gr k = 0 ρ 1 = 0 ρ 2 = 0 ( Y ) k ( a 1 ) ρ 1 k 2 ( Pr ) ρ 1 ( a ) ρ 2 t β ρ 1 + α ρ 2 k 2 + 1 k ! ρ 2 ! Γ ( β ρ 1 + α ρ 2 k 2 + 2 ) Γ ( γ ρ 1 + ρ 2 ) Γ ( γ ρ 1 ) ,
V ˜ 4 ( Y , t ) = Gr λ 1 = 0 λ 2 = 0 ρ 1 = 0 ρ 2 = 0 ( Y ) λ 1 ( a 1 ) ρ 1 ( Pr ) ρ 1 + λ 1 2 ( a ) λ 2 + ρ 2 t β λ 1 2 λ 1 2 + α λ 2 + β ρ 1 + α ρ 2 + 1 λ 1 ! λ 2 ! ρ 2 ! Γ ( β λ 1 2 λ 1 2 + α λ 2 + β ρ 1 + α ρ 2 + 2 ) Γ ( γ λ 1 2 + λ 2 ) Γ ( γ λ 1 2 ) Γ ( γ ρ 1 + ρ 2 ) Γ ( γ ρ 1 ) .

5. Graphically Results and Discussion

The Prabhakar fractional derivative was used to analyze the behavior of Casson fluid using the Laplace transform method. To get a physical understanding of flow parameters, particularly fractional parameters, and to look at some additional graphs from [2,3,4,5,6,7,8,9,10,11,12,13].
For accuracy and validation of temperature and velocity, Figure 2 and Figure 3 are presented. These figures portrayed the comparison between classical solutions obtained by Khalid et al. [42] and the present results when fractional parameters β = γ = 0 and α = 1 . It is found that our general results in the limiting form are validated and correct. Figure 4 and Figure 5 depict the comparison between classical solutions obtained by Khalid et al. [42] and the present results for different values of fractional parameters. The temperature of the fluid for different values of fractional parameters is obtained by considering small and large values of the time. It is evident from the figures that for a small time, the fluid properties can be increased, and for a large time, fluid properties can be decreased.Moreover, it is concluded that the fractional parameters control the boundary layer thickness and show the dual behavior for different values of the time. In Figure 6, we have another comparison that shows maximum velocity can be achieved in the mainstream region away from the plate with the larger values of fractional parameters.
Figure 7 and Figure 8 are plotted for the influence of fractional parameters on temperature and velocity, respectively. The fractional approach is responsible for exhibiting the memory of the fluid properties. It is clear from the figures that for smaller fractional parameters that maximum memory can be achieved at the present time, and for larger values of fractional parameters, the boundary layers approach the classical model. However, it is concluded that the fractional solutions show behavior for a different time. Therefore, we can say that the fractional approach has an advantage in controlling the boundary layers, and this fact can be used in the fitting of real data. As we increase the value of Gr in Figure 9, the velocity increases. This is because when Gr increases, the buoyancy forces become much greater, resulting in more convection and ultimately velocity increases.
When the value of Pr is increased in Figure 10, the temperature drops. Physically, increasing the value of Pr makes the fluid thicker, resulting in the lowest thickness of the thermal boundary layer. Figure 11 depicts the effect of the velocity profile on Pr . Velocity decreases as Pr increases. The temperature and velocity profiles for different larger values of Prandtl number Pr are shown in Figure 12 and Figure 13 when the other parameters are fixed. It is observed that the temperature of the fluid decreases with increasing Prandtl number and velocity of fluid increases with increases of the values of Prandtl number Pr . The thickness of the thermal boundary layer is found to be a minimum near the plate. Fluids with a high Pr number have a higher viscosity and low thermal conductivity limits because of thicker fluid. It is concluded that the Prabhakar fractional model of temperature and velocity field holds for smaller and larger Prandtl numbers. Figure 14 depicts the graphical results for the fractional parameter α and phase angle ω t . The fluid appears to be oscillating between 0 and 1. These variations near the plate are extreme and rise for further values of independent variable y. We have concentrated more on such results to illustrate them on the values of α = 0.2, 0.4, 0.6, 0.8 and ω t = 0, π 4 , π 2 and π .

6. Conclusions

The analytical solutions of fractional partial differential equations that emerge in Casson fluid phenomena are the subject of this research. The Prabhakar fractional model of Casson fluid through generalized Fourier law for oscillating surface with Laplace transform methods. Graphical representations have been used to discuss some of the physical effects of flow characteristics:
  • For fractional parameters, the energy and momentum boundary layer can be controlled with different values of the time.
  • The Prabhakar fractional model holds for small and larger values of the Prandtl number.
  • For validation, recent results are recovered and presented graphically in the graphic section and show good agreement.
  • The fractional derivative approach can be useful in fitting the real data where needed instead of classical approach.
  • In the limiting case, when fractional parameters are taken β = γ = 0 and α = 1 for both velocity and temperature, we find the solutions obtained with ordinary derivatives from the existing literature.
  • This model is a linear fractional model with oscillation at the boundary displaying the laminar flow. It is possible when the fractional model is non-linear with turbulence and not handled with integral transform. Then, in order to handle the turbulence, some numerical methods are required, such as finite element method, RK-4, etc.

Author Contributions

N.S.: writing-original draf, Methodology; M.I.A.: Data curation, Formal analysis, Conceptualization, Supervision; T.S.: Resources, Funding acquisition, Investigation; N.P.: Project administration, Writing—review & editing; T.M.: Visualization, Validation, Software. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by King Mongkut’s University of Technology North Bangkok. Contract No. KMUTNB-64-KNOW-43.

Data Availability Statement

All the data contains with in the manuscript.

Acknowledgments

The authors are grateful to the University of Management and Technology Lahore, Pakistan for facilitating and affirming this research study.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolExplaination
Pr Prandtl number
Gr Groshof number
γ 1 Casson Parameter
p ˜ y Yield stress of fluid
π ˜ c Critical value of non-Newtonian
ρ Density
H ˜ ( t ) Unit step function
μ ˜ β Plastic dynamic viscosity
ω Frequency of oscillation
C p Specific heat
k ˜ Thermal conductivity
iUnit vector
u ˜ Velocity field
β Fractional parameter
τ ˜ Heat capacity ratio
tTime
U ˜ Amplitude of plate oscillation
γ Fractional parameter
T ˜ Temperature
T ˜ w Fluid temperature plate at raised
T ˜ 0 Fluid temperature of plate at rest
α Fractional parameter
θ ˜ Dimensionless temperature

References

  1. Casson, N. A flow equation for pigment-oil suspensions of the printing ink type. In Rheol. Disperse Syst.; Pergamon Press: Oxford, UK, 1959. [Google Scholar]
  2. Pramanik, S. Casson fluid flow and heat transfer past an exponentially porous stretching surface in presence of thermal radiation. Ain Shams Eng. J. 2014, 5, 205–212. [Google Scholar] [CrossRef] [Green Version]
  3. Dash, R.K.; Mehta, K.N.; Jayaraman, G. Casson fluid flow in a pipe filled with a homogeneous porous medium. Int. J. Eng. Sci. 1996, 34, 1145–1156. [Google Scholar] [CrossRef]
  4. Sheikh, N.A.; Ching, D.L.C.; Khan, I.; Kumar, D.; Nisar, K.S. A new model of fractional Casson fluid based on generalized Fick’s and Fourier’s laws together with heat and mass transfer. Alex. Eng. J. 2020, 59, 2865–2876. [Google Scholar] [CrossRef]
  5. Ndolane, S.E.N.E. A new approach for the solutions of the fractional generalized Casson fluid model described by Caputo fractional operator. Adv. Theory Nonlinear Anal. Its Appl. 2020, 4, 373–384. [Google Scholar]
  6. Arif, M.; Kumam, P.; Kumam, W.; Khan, I.; Ramzan, M. A Fractional Model of Casson Fluid with Ramped Wall Temperature: Engineering Applications of Engine Oil. Comput. Math. Methods 2021, e1162. [Google Scholar] [CrossRef]
  7. Ali, F.; Sheikh, N.A.; Khan, I.; Saqib, M. Solutions with Wright function for time fractional free convection flow of Casson fluid. Arab. J. Sci. Eng. 2017, 42, 2565–2572. [Google Scholar] [CrossRef]
  8. Raza, A.; Khan, S.U.; Farid, S.; Khan, M.I.; Sun, T.C.; Abbasi, A.; Khan, M.I.; Malik, M.Y. Thermal activity of conventional Casson nanoparticles with ramped temperature due to an infinite vertical plate via fractional derivative approach. Case Stud. Therm. Eng. 2021, 27, 101191. [Google Scholar] [CrossRef]
  9. Khan, I.; Shah, N.A.; Vieru, D. Unsteady flow of generalized Casson fluid with fractional derivative due to an infinite plate. Eur. Phys. J. Plus 2016, 131, 1–12. [Google Scholar] [CrossRef]
  10. 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]
  11. Koffi, M.; Andreopoulos, Y.; Jiji, L. Heat transfer enhancement by induced vortices in the vicinity of a rotationally oscillating heated plate. Int. J. Heat Mass Transf. 2017, 112, 862–875. [Google Scholar] [CrossRef]
  12. Akcay, S.; Akdag, U.; Palancioglu, H. Experimental investigation of mixed convection on an oscillating vertical flat plate. Int. Commun. Heat Mass Transf. 2020, 113, 104528. [Google Scholar] [CrossRef]
  13. Song, Y.Q.; Raza, A.; Al-Khaled, K.; Farid, S.; Khan, M.I.; Khan, S.U.; Shi, Q.H.; Malik, M.Y.; Khan, M.I. Significances of exponential heating and Darcy’s law for second grade fluid flow over oscillating plate by using Atangana-Baleanu fractional derivatives. Case Stud. Therm. Eng. 2021, 27, 101266. [Google Scholar] [CrossRef]
  14. Ali, F.; Khan, I.; Shafie, S. Closed form solutions for unsteady free convection flow of a second grade fluid over an oscillating vertical plate. PLoS ONE 2014, 9, e85099. [Google Scholar] [CrossRef] [Green Version]
  15. Haider, M.I.; Asjad, M.I.; Ali, R.; Ghaemi, F.; Ahmadian, A. Heat transfer analysis of micropolar hybrid nanofluid over an oscillating vertical plate and Newtonian heating. J. Therm. Anal. Calorim. 2021, 144. [Google Scholar] [CrossRef]
  16. Shah, N.A.; Khan, I. Heat transfer analysis in a second grade fluid over and oscillating vertical plate using fractional Caputo–Fabrizio derivatives. Eur. Phys. J. C 2016, 76, 1–11. [Google Scholar] [CrossRef] [Green Version]
  17. Sheikholeslami, M.; Kataria, H.R.; Mittal, A.S. Effect of thermal diffusion and heat-generation on MHD nanofluid flow past an oscillating vertical plate through porous medium. J. Mol. Liq. 2018, 257, 12–25. [Google Scholar] [CrossRef]
  18. Fang, Z.W.; Sun, H.W.; Wang, H. A fast method for variable-order Caputo fractional derivative with applications to time-fractional diffusion equations. Comput. Math. Appl. 2020, 80, 1443–1458. [Google Scholar] [CrossRef]
  19. Dassios, I.; Baleanu, D. Optimal solutions for singular linear systems of Caputo fractional differential equations. Math. Methods Appl. Sci. 2021, 44, 884–7896. [Google Scholar] [CrossRef] [Green Version]
  20. Ali, R.; Akgül, A.; Asjad, M.I. Power law memory of natural convection flow of hybrid nanofluids with constant proportional Caputo fractional derivative due to pressure gradient. Pramana 2020, 94, 1–11. [Google Scholar] [CrossRef]
  21. Baleanu, D.; Jajarmi, A.; Mohammadi, H.; Rezapour, S. A new study on the mathematical modelling of human liver with Caputo–Fabrizio fractional derivative. Chaos Solitons Fractals 2020, 134, 109705. [Google Scholar] [CrossRef]
  22. Alizadeh, S.; Baleanu, D.; Rezapour, S. Analyzing transient response of the parallel RCL circuit by using the Caputo–Fabrizio fractional derivative. Adv. Differ. Equ. 2020, 2020, 55. [Google Scholar] [CrossRef]
  23. Shah, K.; Jarad, F.; Abdeljawad, T. On a nonlinear fractional order model of dengue fever disease under Caputo-Fabrizio derivative. Alex. Eng. J. 2020, 59, 2305–2313. [Google Scholar] [CrossRef]
  24. Abdeljawad, T.; Hajji, M.A.; Al-Mdallal, Q.M.; Jarad, F. Analysis of some generalized ABC–fractional logistic models. Alex. Eng. J. 2020, 59, 2141–2148. [Google Scholar] [CrossRef]
  25. 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] [PubMed]
  26. Ghanbari, B.; Atangana, A. A new application of fractional Atangana–Baleanu derivatives: Designing ABC-fractional masks in image processing. Phys. A Stat. Mech. Appl. 2020, 542, 123516. [Google Scholar] [CrossRef]
  27. 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]
  28. Benchohra, M.; Bouriah, S.; Nieto, J.J. Terminal value problem for differential equations with Hilfer–Katugampola fractional derivative. Symmetry 2019, 11, 672. [Google Scholar] [CrossRef] [Green Version]
  29. Asif, M.; Ul Haq, S.; Islam, S.; Abdullah Alkanhal, T.; Khan, Z.A.; Khan, I.; Nisar, K.S. Unsteady flow of fractional fluid between two parallel walls with arbitrary wall shear stress using Caputo–Fabrizio derivative. Symmetry 2019, 11, 449. [Google Scholar] [CrossRef] [Green Version]
  30. Polito, F.; Tomovski, Z. Some properties of Prabhakar-type fractional calculus operators. arXiv 2015, arXiv:1508.03224. [Google Scholar] [CrossRef]
  31. 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]
  32. 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]
  33. Eshaghi, S.; Ghaziani, R.K.; Ansari, A. Stability and dynamics of neutral and integro-differential regularized Prabhakar fractional differential systems. Comput. Appl. Math. 2020, 39, 1–21. [Google Scholar] [CrossRef]
  34. Tanveer, M.; Ullah, S.; Shah, N.A. Thermal analysis of free convection flows of viscous carbon nanotubes nanofluids with generalized thermal transport: A Prabhakar fractional model. J. Therm. Anal. Calorim. 2021, 144, 2327–2336. [Google Scholar] [CrossRef]
  35. Alidousti, J. Stability region of fractional differential systems with Prabhakar derivative. J. Appl. Math. Comput. 2020, 62, 135–155. [Google Scholar] [CrossRef]
  36. Derakhshan, M. New Numerical Algorithm to Solve Variable-Order Fractional Integrodifferential Equations in the Sense of Hilfer-Prabhakar Derivative. Abstr. Appl. Anal. 2021, 2021, 8817794. [Google Scholar]
  37. 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. [Google Scholar] [CrossRef]
  38. 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]
  39. 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, 14, 1292–1300. [Google Scholar] [CrossRef]
  40. Sene, N.; Srivastava, G. Generalized Mittag–Leffler input stability of the fractional differential equations. Symmetry 2019, 11, 608. [Google Scholar] [CrossRef] [Green Version]
  41. Iyiola, O.; Oduro, B.; Zabilowicz, T.; Iyiola, B.; Kenes, D. System of time fractional models for COVID-19: Modeling, analysis and solutions. Symmetry 2021, 13, 787. [Google Scholar] [CrossRef]
  42. Khalid, A.; Khan, I.; Shafie, S. Exact solutions for unsteady free convection flow of Casson fluid over an oscillating vertical plate with constant wall temperature. In Abstract and Applied Analysis; Hindawi: London, UK, 2015; Volume 2015. [Google Scholar]
  43. Asjad, M.I.; Basit, A.; Iqbal, A.; Shah, N.A. Advances in transport phenomena with nanoparticles and generalized thermal process for vertical plate. Phys. Scr. 2021, 96, 114001. [Google Scholar] [CrossRef]
  44. Garra, R.; Garrappa, R. The Prabhakar or three parameter Mittag–Leffler function: Theory and application. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 314–329. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Physical model.
Figure 1. Physical model.
Symmetry 13 02039 g001
Figure 2. Comparison of the present results with Khalid et al. [42] on temperature field. When β = γ = 0 and t = 2, α = 1 , Pr = 20, a = 0.01.
Figure 2. Comparison of the present results with Khalid et al. [42] on temperature field. When β = γ = 0 and t = 2, α = 1 , Pr = 20, a = 0.01.
Symmetry 13 02039 g002
Figure 3. The current outcomes are compared with Khalid et al. [42] on velocity field. When β = γ = 0 and t = 0.2, Pr = 15, α = 1 , a = 0.01, ω = π 4 ,   γ 1 = 45 ,   Gr = 0.2 .
Figure 3. The current outcomes are compared with Khalid et al. [42] on velocity field. When β = γ = 0 and t = 0.2, Pr = 15, α = 1 , a = 0.01, ω = π 4 ,   γ 1 = 45 ,   Gr = 0.2 .
Symmetry 13 02039 g003
Figure 4. Comparison of the present results with Khalid et al. [42] on temperature field. When t = 0.2 and a = 0.01, Pr = 20.
Figure 4. Comparison of the present results with Khalid et al. [42] on temperature field. When t = 0.2 and a = 0.01, Pr = 20.
Symmetry 13 02039 g004
Figure 5. Comparison of the present results with Khalid et al. [42] on temperature field. When t = 2 and a = 0.01, Pr = 20.
Figure 5. Comparison of the present results with Khalid et al. [42] on temperature field. When t = 2 and a = 0.01, Pr = 20.
Symmetry 13 02039 g005
Figure 6. The current outcomes are compared with Khalid et al. [42] on velocity field. When t = 0.2 and a = 0.8, Pr = 20, ω = π 4 ,   γ 1 = 45 ,   Gr = 60 .
Figure 6. The current outcomes are compared with Khalid et al. [42] on velocity field. When t = 0.2 and a = 0.8, Pr = 20, ω = π 4 ,   γ 1 = 45 ,   Gr = 60 .
Symmetry 13 02039 g006
Figure 7. The influence of fractional parameters ( α , β , γ ) on temperature field. When t = 0.4, a = 0.002, Pr = 20.
Figure 7. The influence of fractional parameters ( α , β , γ ) on temperature field. When t = 0.4, a = 0.002, Pr = 20.
Symmetry 13 02039 g007
Figure 8. The influence of fractional parameters ( α , β , γ ) on velocity field. When t = 0.7, Pr = 20, ω = π 4 , γ 1 = 0.8 , a = 0.1, Gr = 15.
Figure 8. The influence of fractional parameters ( α , β , γ ) on velocity field. When t = 0.7, Pr = 20, ω = π 4 , γ 1 = 0.8 , a = 0.1, Gr = 15.
Symmetry 13 02039 g008
Figure 9. The impacts of Grashof number Gr on velocity field. When t = 0.3, a = 0.002, Pr = 20, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 .
Figure 9. The impacts of Grashof number Gr on velocity field. When t = 0.3, a = 0.002, Pr = 20, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 .
Symmetry 13 02039 g009
Figure 10. The influence of smaller Prandtl number Pr on temperature field. When t = 1, a = 0.3, α = 0.7 ,   β = 0.7 ,   γ = 0.7 .
Figure 10. The influence of smaller Prandtl number Pr on temperature field. When t = 1, a = 0.3, α = 0.7 ,   β = 0.7 ,   γ = 0.7 .
Symmetry 13 02039 g010
Figure 11. The influence of smaller Prandtl number Pr on velocity field. When t = 0.3, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 , a = 0.1, γ 1 = 0.8 , Gr = 15.
Figure 11. The influence of smaller Prandtl number Pr on velocity field. When t = 0.3, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 , a = 0.1, γ 1 = 0.8 , Gr = 15.
Symmetry 13 02039 g011
Figure 12. The influence of larger Prandtl number Pr on temperature field. When t = 1, a = 0.3, α = 0.7 ,   β = 0.7 ,   γ = 0.7 .
Figure 12. The influence of larger Prandtl number Pr on temperature field. When t = 1, a = 0.3, α = 0.7 ,   β = 0.7 ,   γ = 0.7 .
Symmetry 13 02039 g012
Figure 13. The influence of larger Prandtl number Pr on velocity field. When t = 0.7, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 , a = 0.1 γ 1 = 3.8 ,   Gr = 15.
Figure 13. The influence of larger Prandtl number Pr on velocity field. When t = 0.7, α = 0.5 ,   β = 0.5 ,   γ = 0.5 , ω = π 4 , a = 0.1 γ 1 = 3.8 ,   Gr = 15.
Symmetry 13 02039 g013
Figure 14. The influence of α and ω t on velocity field when t = 0.4, β = 0.5 ,   γ = 0.5 , a = 0.1, γ 1 = 0.5 , Pr = 12, Gr = 2.5.
Figure 14. The influence of α and ω t on velocity field when t = 0.4, β = 0.5 ,   γ = 0.5 , a = 0.1, γ 1 = 0.5 , Pr = 12, Gr = 2.5.
Symmetry 13 02039 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sarwar, N.; Asjad, M.I.; Sitthiwirattham, T.; Patanarapeelert, N.; Muhammad, T. A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law. Symmetry 2021, 13, 2039. https://doi.org/10.3390/sym13112039

AMA Style

Sarwar N, Asjad MI, Sitthiwirattham T, Patanarapeelert N, Muhammad T. A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law. Symmetry. 2021; 13(11):2039. https://doi.org/10.3390/sym13112039

Chicago/Turabian Style

Sarwar, Noman, Muhammad Imran Asjad, Thanin Sitthiwirattham, Nichaphat Patanarapeelert, and Taseer Muhammad. 2021. "A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law" Symmetry 13, no. 11: 2039. https://doi.org/10.3390/sym13112039

APA Style

Sarwar, N., Asjad, M. I., Sitthiwirattham, T., Patanarapeelert, N., & Muhammad, T. (2021). A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law. Symmetry, 13(11), 2039. https://doi.org/10.3390/sym13112039

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