Applications of Prabhakar-like Fractional Derivative for the Solution of Viscous Type Fluid with Newtonian Heating Effect

: This article examines a natural convection viscous unsteady ﬂuid ﬂowing on an oscillating inﬁnite inclined plate. The Newtonian heating effect, slip effect on the boundary wall, and constant mass diffusion conditions are also considered. In order to account for extended memory effects, the semi-analytical solution of transformed governed partial differential equations is attained with the help of a recent and more efﬁcient fractional deﬁnition known as Prabhakar, like a thermal fractional derivative with Mittag-Lefﬂer function. Fourier and Fick’s laws are also considered in the thermal proﬁle and concentration ﬁeld solution. The essentials’ preliminaries, fractional model, and execution approach are expansively addressed. The physical impacts of different parameters on all governed equations are plotted and compared graphically. Additionally, the heat transfer rate, mass diffusion rate, and skin friction are examined with different numerical techniques. Consequently, it is noted that the variation in fractional parameters results in decaying behavior for both thermal and momentum proﬁles while increasing with the passage of time. Furthermore, in comparing both numerical schemes and existing literature, the overlapping of both curves validates the attained solution of all governed equations.


Introduction
Natural convection viscous fluid flows flowing on a vertical plate are widely considered in the literature because of their massive solicitations in different fields of engineering and ecological processes. They are also concentrated in manufacturing submissions such as nuclear reactors, filtration methods, spaceship design, fiber insulation, geothermic schemes, etc. Numerous investigators have examined unsteady natural convective flows past an oscillating and vertical plate with different thermal imposed conditions. Georgantopoulos [1] was the first who find a closed-form solution for natural convection impacts on the flow of a viscous fluid along with a vertical plate. The naturalconvection flows of a viscoelastic fluid with an accelerated plate were discussed by Raptis and Singh [2]. The impacts of the magnetic field were also reserved for contemplation. Free convection impacts considering an exponentially accelerated and infinite vertical plate were deliberated in [3]. Natural convection wavering flow past an infinite permeable plate and continuous suction were considered through Soundalgekar [4]. The communication of natural convection and the vibrating flow with vertical plate and thermal radiation can be seen in [5]. The impacts of heat contamination on the boundary layer with a horizontal plate were examined by Ishak [6,7]. Haq et al. [8] calculated the closed-form solutions of MHD natural convection flow past along a vertical and oscillating plate with thermal flux in a permeable medium with integral transforms.In [9], the authors studied the flow between two parallel plates with the non-isothermal slip effect of a non-linear model. Different parameters with slip coefficients in thermal profiles have also been considered. A non-uniformly heated viscous fluid flow model flowing through a bounded domain was studied by Baranovskii et al. [10].
Natural convection flows have fascinated extensive consideration because of their importance in many fields of science containing biomedical engineering and fluid dynamics. It is often momentous in refrigeration, atmospheric as well as oceanic circulation, concentration, spaceship structure, dehydration, filtration, ventilation of building plan, processing of permeable materials in textile workshops, conserving systems and electronic things for nuclear-powered plants, nuclear reactors, cooled or heated storing places, electric power strategies, and solar collectors [11]. The natural convection mechanism happens because of significant temperature alterations, which can impact the density of the fluid and then source relative buoyancy of the fluid. It has abundant applications in numerous scientific as well as industrial concerns. Natural convection in various systems, actions, and tools is incredibly momentous in exothermal reactions tasks. It is suitable for security attention under different circumstances, where the typical mode continues to flop and when the procedure trusts natural convection to dispose of the designed temperature. It has unusual significance in energy production and digital devices and systems in which such an arrangement is compulsory to avoid extreme heat [12]. In the last few years, extensive analytical and theoretical studies on natural convection have been conducted to understand several environmental features and phenomena in considered frequent scientific circumstances. Javaid et al. [13] deliberated the free convection flowed based on second-grade fluid and stated that with the Grashof number and Prandtl number, the velocity for viscous fluid was more than the velocity estimations of second-grade fluid. The electro-osmotic flow analysis for second-grade fluids and slip conditions was effectively discoursed through Wang et al. [14]. Closed-form solutions were found utilizing integral transform. Numerous particular cases prevailing in the literature were also improved. The natural convection flow of fractional incompressible and second-grade fluid along two upright plates was discussed in the radiation effect between both plates, and an analytical solution was attained by Nisa et al. [15]. The prior studies may perceive countless research topics on the natural convection flow models [16][17][18][19][20][21][22].
In the last few years, fractional calculus has been promising because of its massive implication in different fields of engineering and applied sciences that are not present in the fields of non-fractional calculus, which agrees with a random order of differentiation as well as integration. Ali et al. [23,24], Nehad et al. [25], and Zafar et al. [26] exploited different new definitions of fractional derivatives for the solutions of different Newtonian and non-Newtonian type fluids. Imran et al. [27] utilized the AB-fractional model to excavate out the semi-analytical solution of the Maxwell fractional model. Nadeem et al. [28] examined the AB and CF factional model and attained an analytical solution of Cassontype nanofluid. Various physical phenomena constructed on differential equations with fractional calculus are immensely used for the modeling of numerous everyday life physical problems, as fractional calculus takes memory impacts, such as problems in relaxation, oscillation, diffusion, fluid flow, retardation, dynamical processes, engineering processes, where classical models cannot predict the previous state of processes [29][30][31][32][33][34].
The construction of numerous physical phenomena consequences increases the complexity of partial differential equations. The calculation of solutions for such complicated equations is significant in investigating such physical phenomena. Fractional calculus is the critical division of mathematics that effectively reports diverse methods for the simulation of such types of problems. The complex mathematical expressions, including the local and non-local kernels, can be efficiently tackled through a fractional approach. The Caputo-Fabrizio (CF) and Atangana-Baleanu (AB) methods are the most well-known apparatuses in fractional calculus that different investigators have extensively applied in past years. Tialk Raj Prabhakar was an Indian mathematician who proposed a novel three-parameter fractional operator:the Mittag-Leffer function with three different fractional operators. This Mittag-Leffer operator successfully applies conventional kernels [35]. Sulaiman et al. [36] found the solution for Burger's equation through diverse fractional derivatives comprising MLF kernels. AB and CF fractional derivatives and Mittag-Leffer and exponential functions were compared [37]. Singh et al. [38] studied the Cattaneo-Christov derivatives and solved fractional diffusion equations by exploiting the Hilfer-Prabhakar operator. Samraiz et al. [39] deliberated the (k,s)-Hilfer-Prabhakr fractional operator for diverse types of mathematical fractional problems. In [40], Basit et al. employed the Prabhakar fractional technique to examine the solution of second-grade fractional nanofluid, taking different types of nanoparticles. Rehman et al. [41] studied the free convection Maxwell fluid flow under the Newtonian heating effect using the Prabhakar fractional derivative technique of the Mittag-Leffler function, in which they also considered Fourier and Fick's law to investigate the solution of thermal and concentration profiles. In [42], the authors inspected the free convection fluid flow with generalized thermal transport and carbon nanotubes as nanoparticles using a recent fractional derivative definition.
The main focus of current research is to study the natural convection, incompressible, viscous, and unsteady flow on an infinite oscillating inclined plate under the influence of an applied inclined magnetic field in the sense of the Prabhakar fractional derivative operator. Newtonian heating influences are engaged into focus in the fractional thermal flow model to check the thermal performance. The Prabhakar fractional derivative operator technique with Laplace transform is followed to accomplish the fractional and numerical simulations. To achieve the inverse of the Laplace method, two important diverse methods T'zous method and Stehfest, are applied. In the end, the thermal results are enumerated for different flow characters.

Problem Description
Suppose a natural convection, incompressible, viscous, and unsteady flow on an infinite oscillating inclined plate under the impact of the applied inclined magnetic field with strength B o . As the inclined plate is in the xy-plane, therefore, all the governed equations are functions of y and t. Initially, at t = 0, the inclined plate and fluid are both in a constant situation with a constant temperature T ∞ and concentration of C ∞ . With time t > 0 + , the fixed plate starts to oscillate with particular continual velocity f (t), with f (0) = 0, and temperature and concentration levels also rise with time, as shown in Figure 1. In light of the above conditions and neglecting the viscous dissipation and pressure gradient, the governed partial differential equations with their physical boundary conditions can be modeled by Boussinesp's approximation as follows [43].
With its dependable physical boundary conditions, as follows into the proceeding governed equations and forgetting the "*" notation. We attain the succeeding non-dimensional formulas as  ∂v (y,t) With its dependable physical boundary conditions, as follows v (y,0) = 0, T (y,0) = T ∞ , C (y,0) = C ∞ ; y ≥ 0 where δ (y,t) , J (y,t) identifies the thermal flux rate by Fourier's law and Fick's law, respectively. Now, the non-dimensional governing Equations (1)-(5) and the corresponding consistent conditions present the subsequent non-dimensional values: into the proceeding governed equations and forgetting the "*" notation. We attain the succeeding non-dimensional formulas as ) with the following dimensionless physical conditions where , Pr e f f = Pr Re , Sc e f f = Sc Re The Navier slip coefficient and Newtonian heating effect are also considered in the above-mentioned physical conditions of the flowing fluid and the fluid temperature, concentration, and velocity will be zero at y → ∞ . This article will discuss a generalized mathematical model in which the memory effect and shear stress are also under consideration in the sense of the Prabhakar fractional operator. The respective mathematical preliminaries will define the basic concept of the fractional scheme and the relation of the Prabhakar fractional technique with other fractional models.

Basic preliminaries:
The one-parametric MittagLeffler function with mathematical form as was studied by Mittag-Leffler [44]. Then after some time, Wiman [45] explored the more generalization form of the one-parametric function, known as two-parametric Mittag-Leffler functions with the mathematical form as follows: In [46], authors introduced the three-parametric Mittag-Leffler function, which is commonly known as the Prabhakar fractional derivative ; α, β, γ, z ∈ C, Re(α) > 0 with the basic properties is identified as the Prabhakar kernel.
(Prabhakar Integral) The Prabhakar integral can be defined as [47,48] Some applicable fractional constraint cases can be encapsulated as when β > 0, α = 0, then the property three repeats as L e γ α,β (0) = q −β (The regularized Prabhakar derivative) In [47,48], the regularized Prabhakar derivative is distinct as where C D γ α,β,α signifies the Prabhakar derivative operator and g m represents the mth derivative of g(t). The Laplace transformation of generalized Prabhakar and its kernel can be derived as The primary relations and connections with different fractional operators can be encapsulated as follows.

•
When β ≥ 0, γ = 0, the Prabhakar derivative will transform into the Caputo derivative When α = β = 1, γ = −1, then the relation between Prabhakar and CF derivative will , then the connection between the AB derivative and Prabhakar derivative will develop as C D −1 We may obtain the traditional Fourier's law by taking β = γ = 0. Moreover, because Fourier's law of thermal conductivity primarily determines the Prabhakar fractional derivative, Fourier and Fick's laws in the context of the Prabhakar derivative are as follows.

Solution of the Energy Profile
Appling the LT on transformed Equations (10) and (21) and utilizing its matching conditions for the solution of the energy profile By inserting these transmuted conditions, the solution of the thermal profile is as follows: We employ numerical approaches, namely, Stehfest and Tzou's schemes in Tables 1 and 2, to analyze the inversion of Equation (25).  Classical Solution of the Energy Field (β = γ = 0) For the classical solution of thermal profile take β = γ = 0, so δ(t) denotes the Dirac's Delta distribution. Due to this, the generalized Fourier's law is changed into classical Fourier's law. Additionally, 1

Solution of the Concentration Profile
With the same methodology as used for the energy field, utilizing the LT on Equations (12) and (22), Using the above conditions and after simplification, the solution of Equation (28) becomes Again, the inversion of Laplace of the above equation is examined numerically in Tables 1 and 2.
Classical Solution of the Concentration Profile (β = γ = 0) For the classical solution of the concentration profile, again, take β = γ = 0. The generalized Fick's law is changed into classical Fick's law as a result of this.
√ KSc e f f y Erfc

Solution of Momentum Field
In this part, the velocity equation solution is determined using the same approach as the energy equation solution. Using the LT technique on Equation (9) and its accompanying conditions, we obtain Employing these conditions, the solution of the momentum equation turns out to be 1+h Pr e f f q

Req
1+h Sc e f f (K+q) Nusselt number, Sherwood number, and skin friction are as follows: Various writers have utilized various numerical inverse methods to calculate the Laplace inverse. Recently, Ali et al. [33], and Tiwanaet al. [49] used the Stehfest and Tzous algorithms for the inversion of Laplace of governed equations solutions. Similarly, Aleem et al. [50], Chu et al. [51], and Asjadet al. [52] also used different numerical techniques for the solution of different hybrid nanofluid, Brinkman type nanofluid, and natural convection flowing fluids. As a result, we additionally employed Stehfestand Tzou's methods to quantitatively examine the temperature, concentration, and velocity profile solution. Gaver-Stehfest and Tzou's [53] algorithm can mathematically be inscribed as where N is a positive integer, and

Results and Discussion
This paper investigates the flow of a viscous, incompressible, and unstable fluid under the influence of an angled magnetic field. In order to account for extended memory effects, a fractional model was built using the recently proposed Prabhakar fractional function and a Mittag-Leffler kernel in the governing equations. Then, the solution of all governed equations containing thermal, concentration, and momentum profiles in the sense of Fourier and Fick's laws is derived using the integral transform scheme, namely Laplace transformation, and various inverse numerical algorithms, i.e., Stehfest and Tzou's algorithms are used to obtain the Laplace inverse of guided equations. Finally, graphs 2-15 are drawn to investigate the physical impact of different constraints such as α, β, γ, Pr e f f , Sc e f f , Gr, Gm, Re, M and quantitatively investigate the heat transfer rate, mass diffusion rate, and skin friction, which are examined with different numerical techniques in Tables 1 and 2. Figure 2 highlights the impact of fractional constraints of the Prabhakar fractional derivatives on the thermal profile for different time values. It can be seen that by increasing the values of fractional parameters, the thermal profile represents the decaying behavior. Using fractional derivatives, it is significant to check the influence of time on governed equations. With time, the fluid temperature increases asymptotically near the plate and has a more significant effect for large values of time. Figure  Like the thermal profile, again, the variation in fractional parameters and Schmidt number results in a decline in the concentration profile. For large values of the Schmidt number Sc, viscous forces overwhelm the diffusional effects. The Schmidt number represents mass transfer via diffusion and the relative effectiveness of momentum in concentration and speed in free convection flow regimes. As a result, momentum diffusion will effectively counterbalance large values of Sc, since molecular diffusivity will be decelerated, and the viscosity effects will be enhanced. As a result, with high Sc values, the concentration profile of the boundary layer is reduced.  The effects of fractional parameters on the momentum profile of the flowing fluid are highlighted in Figure 6. It can be seen that the fluid velocity decreases by increasing the value of fractional constraints and increases near the plate with increasing the value oftime and reaches its maximum point; it slowly decreases along the y-axis and reaches zero for higher values of y. The fluid velocity increases near the plate, and the boundary thickness increases with the enhancement in the values of time t. Physically, the increment in the values of factional constraints relates to the thickness of the momentum and thermal profiles. The variation in the values of fractional parameters reduces the thickness of the boundary layer of flowing fluids, due to which the thermal and momentum field shows a decrement trend for this parameter. Figure 7 explains the influence of adequate Prandtl numbers on the momentum profile of the flowing fluid. Again, the momentum field represents decaying behavior with the increment in the value of Prandtl numbers. Physically, the enhancement in the Prandtl numbers means the increment in the viscosity of the fluid and decrement in the thermal conductivity of the flowing fluid, which results in decay in the momentum field of the fluid flow. The variation in velocity profile due to the heat Grashof number and mass Grashof number can be seen in Figures 8 and 9, respectively. Due to the increment in the Grashof number, the buoyancy effect also increases, which increases the thickness of the boundary layer of the fluid, so by varying the values of both parameters, the fluid velocity also increases and shows a unique maximum point near the plate, varying the value of time. In addition, the velocity increases with higher time values and the highest finding near the plate. The fluid motion slows down asymptotically as y increases, satisfying the boundary conditions. Figures 10 and 11 are plotted to see the influence of Sc e f f and Re for the velocity profile for different values of the time. An increment in the Schmidt number means a decrease in the molecular diffusion of the boundary layer, which decreases the momentum profile.       The effects of fractional parameters on the momentum profile of the flowing fluid are highlighted in Figure 6. It can be seen that the fluid velocity decreases by increasing the value of fractional constraints and increases near the plate with increasing the value oftime and reaches its maximum point;it slowly decreases along the y-axis and reaches zero for higher values of .The fluid velocity increases near the plate, and the boundary thickness increases with the enhancement in the values of time t. Physically, the incre-               The effects of the applied magnetic field and the angle of inclination of the applied magnetic field are examined in Figures 12 and 13. It can be seen that with the variation of both parameters, the momentum profile decays. Physically, the variation in the magnetic field enhances the Lorentz force, resulting in a decrease in the momentum field, and the maximum strength of the applied magnetic field is at the right angle, as shown in the figures. The numerical comparison of concentration, temperature, and momentum profiles with different numerical inversion techniques can be examined in Table 1. They are found to be in good agreement due to the very close results of both numerical techniques. Moreover, the numerical analysis of the heat transfer rate, mass transfer rate, and skin friction at different times is examined in Table 2. It can be predicted that both the Nusselt number (heat transfer rate) and skin friction decrease continuously while increasing with time.
The effects of the applied magnetic field and the angle of inclination of the applied magnetic field are examined in Figures 12 and 13. It can be seen that with the variation of both parameters, the momentum profile decays. Physically, the variation in the magnetic field enhances the Lorentz force, resulting in a decrease in the momentum field, and the maximum strength of the applied magnetic field is at the right angle, as shown in the figures. The numerical comparison of concentration, temperature, and momentum profiles with different numerical inversion techniques can be examined in Table 1. They are found to be in good agreement due to the very close results of both numerical techniques. Moreover, the numerical analysis of the heat transfer rate, mass transfer rate, and skin friction at different times is examined in Table 2. It can be predicted that both the Nusselt number (heat transfer rate) and skin friction decrease continuously while increasing with time. The effects of the applied magnetic field and the angle of inclination of the applied magnetic field are examined in Figures 12 and 13. It can be seen that with the variation of both parameters, the momentum profile decays. Physically, the variation in the magnetic field enhances the Lorentz force, resulting in a decrease in the momentum field, and the maximum strength of the applied magnetic field is at the right angle, as shown in the figures. The numerical comparison of concentration, temperature, and momentum profiles with different numerical inversion techniques can be examined in Table 1. They are found to be in good agreement due to the very close results of both numerical techniques. Moreover, the numerical analysis of the heat transfer rate, mass transfer rate, and skin friction at different times is examined in Table 2. It can be predicted that both the Nusselt number (heat transfer rate) and skin friction decrease continuously while increasing with time.

Validation of Attained Results
The comparison of the two numerical systems, Stehfest and Tzou's, is examined by drawing Figure 14a-c for all temperature, concentration, and momentum profiles. A slight overlap of findings between the two curves validates the attained results. In Figure  15a,b, a comparison of the solutions for temperature and velocity fields using the Prabhakar fractional methodology is shown with the work of Imran et al. [54]. The simulations generated by employing the Prabhakar fractional model have good accuracy when compared to the study of Imran et al. [54].

Validation of Attained Results
The comparison of the two numerical systems, Stehfest and Tzou's, is examined by drawing Figure 14a-c for all temperature, concentration, and momentum profiles. A slight overlap of findings between the two curves validates the attained results. In Figure  15a,b, a comparison of the solutions for temperature and velocity fields using the Prabhakar fractional methodology is shown with the work of Imran et al. [54]. The simulations generated by employing the Prabhakar fractional model have good accuracy when compared to the study of Imran et al. [54].

Validation of Attained Results
The comparison of the two numerical systems, Stehfest and Tzou's, is examined by drawing Figure 14a-c for all temperature, concentration, and momentum profiles. A slight overlap of findings between the two curves validates the attained results. In Figure 15a,b, a comparison of the solutions for temperature and velocity fields using the Prabhakar fractional methodology is shown with the work of Imran et al. [54]. The simulations generated by employing the Prabhakar fractional model have good accuracy when compared to the study of Imran et al. [54]. Fractal Fract.2022, 6,

Conclusions
The Prabhakar-like thermal fractional technique is used in this study to examine the problem of a viscous, incompressible, and unsteady fluid flowing across an oscillating inclined plate. In order to account for extended memory effects, a recent and more efficient fractional definition, namely the Prabhakar fractional derivative, is utilized with the Mittag-Leffler kernel. Non-dimensional fractional governed equations are solved using the LT approach; an integral transform method, namely the Laplace inverse of governed equations, is computed using a variety of numerical approaches. The impact of various restrictions on leading equations is visually and numerically investigated. The primary outputs derived from graphical and numerical representation can be bulleted as follows: • The impact of larger values of fractional parameter and an adequate Prandtl number declines the profiles of temperature distributions.

•
The boundary layer concentration also decays with the enhancement in fractional parameter and Schmidt number.

Conclusions
The Prabhakar-like thermal fractional technique is used in this study to examine the problem of a viscous, incompressible, and unsteady fluid flowing across an oscillating inclined plate. In order to account for extended memory effects, a recent and more efficient fractional definition, namely the Prabhakar fractional derivative, is utilized with the Mittag-Leffler kernel. Non-dimensional fractional governed equations are solved using the LT approach; an integral transform method, namely the Laplace inverse of governed equations, is computed using a variety of numerical approaches. The impact of various restrictions on leading equations is visually and numerically investigated. The primary outputs derived from graphical and numerical representation can be bulleted as follows: • The impact of larger values of fractional parameter and an adequate Prandtl number declines the profiles of temperature distributions.

•
The boundary layer concentration also decays with the enhancement in fractional parameter and Schmidt number.

Conclusions
The Prabhakar-like thermal fractional technique is used in this study to examine the problem of a viscous, incompressible, and unsteady fluid flowing across an oscillating inclined plate. In order to account for extended memory effects, a recent and more efficient fractional definition, namely the Prabhakar fractional derivative, is utilized with the Mittag-Leffler kernel. Non-dimensional fractional governed equations are solved using the LT approach; an integral transform method, namely the Laplace inverse of governed equations, is computed using a variety of numerical approaches. The impact of various restrictions on leading equations is visually and numerically investigated. The primary outputs derived from graphical and numerical representation can be bulleted as follows:

•
The impact of larger values of fractional parameter and an adequate Prandtl number declines the profiles of temperature distributions.

•
The boundary layer concentration also decays with the enhancement in fractional parameter and Schmidt number.

•
The momentum profile is an increasing function for Re, Gr, Gm ,while it decreases with the variation in Pr e f f , Sc e f f , M, θ 1 and Prabhakar fractional parameters.
• Thermal profile, concentration, and momentum profiles asymptotically increase with time.

•
The overlapping of both numerical schemes validates the attained solution of all governed equations.

•
The momentum profile is maximal near the plate. It approaches its distinctive peak values in the stream region and then decreases away along the y-axis.

•
The rate of heat transfer, mass transfer, and skin friction varies with the increment in time values.

•
In the comparison of numerical techniques and with the attained results of Imran et al. [54], the overlapping of both curves validates the attained results of this study.

•
As the Prabhakar fractional derivative is the more recent definition of the fractional derivatives technique, it has more efficient and accurate results as compared to other fractional operators as depicted in the comparison of Imran et al. [54].