Significance of Free Convection Flow over an Oscillating Inclined Plate Induced by Nanofluid with Porous Medium: The Case of the Prabhakar Fractional Approach

Given the importance and use of electrically conducted nanofluids, this work aims to examine an engine-oil-based nanofluid including various nanoparticles. In the current study, a fractional model for inspecting the thermal aspect of a Brinkman-type nanofluid, composed of (molybdenum disulfide (MOS2) and graphene oxide (GO) nanoparticles flows on an oscillating infinite inclined plate, which characterizes an asymmetrical fluid flow, heat, and mass transfer. Furthermore, the Newtonian heating effect, magnetic field, and slip boundary conditions were taken into account. The objectives for implementing the Prabhakar-like fractional model are justified because this fractional algorithm has contemporary definitions with no singularity restrictions. Furthermore, the guided fractional model was solved using the Laplace transform and several inverse methods. The obtained symmetrical solutions have been visually analyzed to investigate the physics of several relevant flow parameters on the governed equations. Some exceptional cases for the momentum field are compiled to see the physical analysis of the flowing fluid symmetry. The results show that the thermal enhancement can be progressively improved with the interaction of the molybdenum disulfide-engine oil-based nanofluid suspension, rather than with the graphene oxide mixed nanoparticle fluid. Furthermore, the temperature and momentum profiles enhance due to the factional parameters for molybdenum disulfide and the graphene oxide-engine oil-based nanofluid suspension. This study’s graphical and numerical comparison with the existing literature has shown a very close resemblance with the present work, which provides confidence that the unavailable results are accurate. The results show that an increase improved the heat transmission in the solid nanoparticle volume fractions. In addition, the increment in the mass and heat transfer was analyzed in the numerical evaluation, while the shear stress was enhanced with the enhancement in the Prabhakar fractional parameter α.


Introduction
Due to the numerous uses for thermal energy across all industries, such as the biomedical and engineering fields, and the current exceptionally high demand for thermal energy, researchers have been motivated to accelerate the rate of thermal perturbation from the pending sources by the extraordinary vision of nano-science. Enhancement in heat effectiveness can have an advantage over solicitations incorporating the analysis of plasma, microelectronic equipment, such as nuclear reactors, space refrigeration, microchips in computers, power generation, and many more [1]. Non-Newtonian fluids possess broad applications in science and technology, such as in cosmetics, biochemical, and medical applications. Because of their extensive applications in different fields, investigators are constantly fascinated by these fluids [2,3]. Therefore, different models in the past for non-Newtonian fluids are studied. One of them is known as the Brinkman-type fluid, which Darcy introduced for those fluids that pass through small permeable surfaces [4]. More precisely, this fluid model is appropriate for the flow on a body with a negligible permeability. However, Darcy's law does not apply to the flow on a body or surfaces with a very high porosity. Using a fluid model across the porous surfaces was an idea put forth by Brinkman [5,6]. A viscous, incompressible fluid with a significant amount of porosity, is known as a Brinkman-type fluid. Minimal research in this regard, has been achieved, through the Brinkman-type fluid. In particular, convection flows with Brinkman-type fluid and heat and mass transmission are much less studied, while these investigations have abundant applications in many industries [7,8].
The sharp rise in energy claims worldwide has directed increased efforts to accomplish energy-effective arrangements. The use of innovative schemes to expand the thermal concert of heat exchangers is one of the furthermost serious and modern subjects in the modern world. However, most of these approaches are accompanied by a rise in the pressure drop, called defect. Specifically, the enhancement procedure is hydraulic and thermally supportive [9]. Recently, nano-structure carbon mixtures, such as graphite, graphene, and some metals, such as molybdenum and magnesium, have been utilized to stock hydrogen gas. Nanofluids are simply latent heat-transmitting flowing liquids with enhanced thermophysical properties, and a heat transfer execution that may be used in various devices for more efficient operations. They are used to improve the thermal conductivity of base fluids, such as oil, ethylene glycol, water, propylene glycol, and so on, because they are poor heat-transmitting fluids. Temperature reduction, cancer treatment, and manufacturing use various biomedical and scientific engineering fields. Thermic conduction of a schematic heat-transporting fluid is enriched by the interruptions of significant molecules, which increases the heat transfer coefficient. Because solid metallic components have more thermic conduction than base fluids, suspended constituents are used to promote thermic conduction and heat transmission. The significant advantages of nanofluids are that they have a sufficient viscosity and are auxiliary stable with enhanced wetting, flowing, and diffusion characteristics along the solid surfaces, even for tiny immersion nanoparticles [10]. Srivastava et al. [11] investigated a biological population fractional model. The blended homotopy methods relating to the Sumudu transform are operated to find the solutions and show that the prey's population abruptly declines with time. Kumar et al. [12] derived a new numerical technique to obtain the numerical solution of fractional PDEs, including the Caputo-Fabrizio derivatives, compared their obtained numerical results with the analytical results, and determined that their anticipated numerical technique attains accurate results. Mahanthesh [13] numerically studied the nanofluid flow and heat transfer using the modified Buongiorno model. The author proved that the interruption of the nanoparticle's growth increases the thermal conduction and consequently expands the temperature. Rana et al. [14] investigated the flow of ethylene glycol-based titania nano liquid using the finite difference method, and noted that the accumulation of nanoparticles increases the temperature, and the velocity field is reduced. The finite difference scheme was used to examine the radiative flow of a polar nonliquid, along with a stretchable sheet in [15]. It was observed that the heat generation parameter has a decreasing impact on the Nusselt number. The impact of the nanoparticle combination upon the 3D flow of the titanium nonliquid, because of an exponentially extended surface, was studied in [16], by the shooting method and showed that the impact of the magnetism is much more in the ordinary fluid than in the nanofluid. Abro et al. [17] studied the influences of the magnetized nanoparticles, analytically, on the natural convection flow of the nanofluid using the Laplace method. They proved that the improvement in volume fraction increases the thickness of the thermal boundary layer. For this article's sake, we incorporate new studies, based on the exact and numerical fractional approaches [18][19][20][21][22][23][24][25]. Shafee et al. [26] used computational methods to forecast the turbulent nanofluid (CuO/H 2 O) flow and heat transfer, and showed that the inclusion of helical tabulator measurement enhances the core nanofluid movement. Farshad and Sheikholeslami [27] discussed the turbulent nanomaterial flow using the finite volume method and showed that by enhancing friction because of the increasing inlet velocity, the quantity of the Bejan number decreases. Hussanan et al. discussed and evaluated the MHD's inherently unstable flow of nanofluids passing through an accelerating upright plate, set in a poriferous medium [28]. Sheikholeslami et al. [29] investigated the effect of the persuaded magnetic flux on heat conduction using the KKL (Koo Kleinstreuer-Li) relationship on the nanofluid flow. Akyürek et al. [30] examined the experimentally turbulent forced convective thermal transmission, as well as the pressure drop properties of the nanofluids through a concentric tube. They proved that nanofluids with small particle concentrations did not demonstrate a significant impact on the pressure drop variation.
Fractional calculus is an appropriate tool to explain various anomalous dynamics and processes in complex media, in physics, rheology, and electrical engineering [31][32][33][34][35][36]. Fractional calculus techniques are also applied to evaluate the anomalous relaxation phenomena [37]. In current years, it has been determined that some fractional operators and non-singular kernels are anticipated. For instance, the Atangana-Baleanu derivative and the Mittag-Leffler kernel [38] have attracted a lot of attention. A fractional derivative with a Prabhakar-like kernel is recommended by Giusti and Colombaro [39], which is the most extensive form of the Atangana-Baleanu derivative. Mahanthesh et al. [40] studied an unsteady nonlinear Casson fluid and an oscillating plate, analytically, with the Laplace technique. They proved that the nonlinear convective flow has a substantial impact on heat, as well as mass transfer properties. Abro et al. [41] studied a fractional model of the unsteady natural convection radiating flow, along with the Fourier sine transform, and showed that enhancing the fractional parameter velocity has no symmetrical behavior. Saqib et al. [42] analyzed the Brinkman-type nanofluid with CNTs and ramped the temperature conditions. Khan et al. [43] discussed an advection-reaction diffusion fractional model numerically using the Laplace, forward Euler, and Lagrange interpolation methods. They delivered a base for their considered model's uniqueness, existence, and the HU-stability analysis. Abro et al. [44] investigated the thermal influences of the MHD micropolar fluid which was secreted by the Caputo-Fabrizio derivative. They found the exact solutions by applying the Fourier sine and the Laplace approach. They concluded that the velocity decreases while growing the magnetic parameters. Siddiqui et al. [45] studied the waterbased ferrofluids, numerically, by employing a unique finite difference technique. They noted that the magnitude of the secondary vortices declines if the magnetic strength is boosted. Pandey et al. [46] discussed a space-time advection-reaction-diffusion fractional equation with the homotopy perturbation method. They presented that solute concentration when the considered system tended to move from the fractional-order to integer-order, in the existence of the sink and source terms. They claimed that their proposed method is more suitable for dealing with a fractional model, taking Mittag-Leffler's fractional derivative and estimating the motion of the diverse models ascending in various science and engineering disciplines. Ahokposi et al. [47] recently investigated the groundwater fractal flow with the fractional differentiation and the Mittag-Leffler rule. They employed three distinct numerical strategies: implicit, explicit, and the Crank-Nicholson approach. Khan et al. [48] focused on the flow of generalized Casson fluids with the fractional derivatives, and discovered the closed-form solutions in the Wright function. Khan et al. [49] used the Caputo-Fabrizio derivatives in the heat transfer study of a Maxwell fluid in, another paper. This is an appealing research subject because of the major uses of nanofluids in cooling loads, solar thermal energy, electronic cooling, and solar energy engineering.
According to a review of the literature, as mentioned above, no studies have been carried out on the MHD Brinkman-type nanofluids (such as graphene oxide and molybdenum disulfide nanoparticles with a viscoelastic base fluid), using the Newtonian heating effect, slip boundary conditions, and our suggested setting, in combination with the Prabhakar time-fractional approach to obtain the better energy transfer phenomenon. To address this research gap, a Prabhakar time-fractional Brinkman-type fluid model for the flow of (MoS 2 ) and (GO) based nanoparticles under the influence of an externally inclined magnetic field, was investigated. According to the available literature, the integer order derivatives are local in nature, but the time-fractional derivatives are non-local and preserve the memory property. With such motivations, the current research attempts to offer Prabhakar with a fractional model for a nanofluid problem, involving multiple types of nanoparticles. A close thermal inspection of the modified nanofluid model is offered by applying molybdenum disulfide and graphene oxide nanoparticles with a viscoelastic base fluid. The comparative thermal data for these various nanoparticles is provided to further understand the energy thermal transfer mechanisms. The primary phrases and concepts of the Prabhakar fractional model are presented first, followed by the application of this novel technique in the simulated situation. In two respects, this study is innovative. It first demonstrates the application of the Prabhakar fractional model to various coupled differential systems. Second, this model confirms the thermal applications of different nanoparticles to improve the heating transportation. Furthermore, temperature, concentration, and velocity differences are depicted graphically to demonstrate the effect of the controlling factors.

Mathematical Formulation
Consider a free convection Brinkman-type nanofluid, mixed with (molybdenum disulfide, graphene oxide) nanoparticles and base engine oil fluid, moving over a poured infinite plate has been shown in Figure 1. The x-axis is drawn vertically up the plate, whereas the y-axis is orthogonal to it. The pressure gradient is supposed to be ignored in the presented model, and the moving fluid is electrically conducted. An angled magnetic field, sliding, and the Newtonian heating effects are also considered. Because the magnetic Reynolds number is assumed to be very small, the induced magnetic field is ignored, the base is assumed to be fluid, and the suspension nanoparticles are in the thermal equilibrium. At t = 0, both the plate and the combined fluid are at rest, with the ambient temperatures and concentrations T ∞ and C ∞ . Due to the increase in the temperature variation and the oscillations of the pored oscillating plate, the free convection develops after a brief interval of time t > 0 + , and the constant nanofluid begins to move on the inclined plate. The constant plate vibrates with the constant velocity f (t), having some of its Laplace. We made the following assumptions.

•
Except for the impact of the body action term, all fundamental fluid parameters are supposed to be fix.

•
An applied magnetic field with a strength of B 2 o is inclined with θ 1 as the inclination of the magnetic field.

•
Because the fluid's conductivity is considered low, the magnetic Reynolds number is less than one, and the induced field is small, compared to the transverse magnetic field.

•
It is also assumed that the temperature, concentration, and velocity depend on y and t.

•
It is also assumed that there is still no applied voltage, since the electric field is nonexistent [50]. • Based on the above presumptions and Boussinesq's approximations [51], the governed equations can be stated as follows [42,50]. field. • It is also assumed that the temperature, concentration, and velocity depend on y and t.

•
It is also assumed that there is still no applied voltage, since the electric field is nonexistent [50]. • Based on the above presumptions and Boussinesq's approximations [51], the governed equations can be stated as follows [42,50].
Thermal equation: Fourier law of thermal flux: Diffusion balance equation: Momentum equation: Thermal equation: Fourier law of thermal flux: Diffusion balance equation: Fick's law: Furthermore, the identical physical conditions are specified as [52,53]: where w (y,0) is the fluid velocity, ρ n f symbolizes the effective density of the nanofluid, µ n f is the effective dynamic viscosity of the nanofluid, B indicates the magnetic field, β Tn f is the effective thermal volumetric coefficient, β Cn f is the effective solutal volumetric coefficient, k n f is the effective thermal conductivity, and D indicates the diffusion coefficient. The mathematical relations for ρ f , µ f , σ f , k f regular nanofluid parameters with different thermal characteristics are presented in Table 1. Table 1. Model for the thermophysical properties of the nanofluid quantities [54].

Thermal Features Nanofluid
Density The subscripts n f , f , s in the above table denote nanofluid, base fluid, and solid, all regarded as separate nanoparticles. We implement the following non-dimensional quantities to non-dimensionalize the given unsteady problem: The non-dimensional restrictions indicated above, to reduce the number of variables in the controlled Equations (1)-(5) and ignoring the " * " notation, one obtains Sc where the equivalent conditions are incorporated as Equation (9) represents the constant fluid motion, thermal profile, and concentration field in the above-transformed conditions at t = 0, as supposed in the problem formulation section. In Equation (10), the rate of the fluid flow with the constant fluid velocity f (t) with slipping, the boundary effect is implemented. Furthermore, the thermal field rate with a constant concentration in the Newtonian heating effect is considered. Moreover, in the last boundary equation (11), the zero-fluid movement, temperature, and concentration at y → ∞ are considered. Finally, the thermal features of the considered base fluid (engine oil) and the nanoparticles (molybdenum disulfide, graphene oxide) are examined in Table 2. Table 2. Thermal properties of the base liquids and nanoparticles [55].

Prabhakar Fractional Derivative Scheme
A practical and contemporary mathematical fractional approach has been used in this work, from which the thermal memory effect may be examined. As a result, the Prabhakar fractional derivative is presented here, which is based chiefly on the modified Fourier's and Fick's laws of thermal conductivity [56] where C D γ α,β,α is the Prabhakar fractional derivative, which may be derived numerically [56,57]  where h (n) represents the nth derivative of h(t) ∈ AC n (0, b), AC n (0, b) and means a realvalued function that has continuous derivatives of the (n − 1) order with interval (0, b).
and by taking β = γ = 0 the classical Fourier law can also be obtained.

Solution of the Temperature Profile
Because the criterion for is assumed to be β ∈ [0, 1 ), use m = 0 in the preceding formula of Equation (14). Furthermore, for the temperature field solution, the Laplace transformation scheme on Equations (7) and (12) and its accompanying conditions are used By plugging Equation (16) into (15) and using the above-transformed conditions, the thermal profile may be solved as follows Tables 3 and 4 show that the Laplace inverse of the aforementioned thermal profile solution will be numerically investigated using the Stehfest and Tzou techniques.

Solution of the Concentration Profile
For the solution of the concentration of the boundary layers, using the LT on the non-dimensionalized and converted with the leading Equations (8) and (13), and using the result of Equation (11), the solution of C (y,s) overcomes that.  Inserting Equation (20) into (19) and applying the above-transformed conditions, the concentration field will be Again the Laplace inverse of the above-attained solution will be performed, numerically, in Tables 3 and 4.

Solution for the Velocity Profile
This section examines the semi-analytical solution of the momentum profile with its physically transformed conditions. Bearing in mind the result of Equation (13), employing the LT on Equation (6), the non-homogeneous ordinary differential equation from the momentum field will be as follows w (0,s) − ∂w (y,s) ∂y y=0 = F(s); w (y,s) → 0 as y → ∞ By interpolating the above-considered conditions and simplifying Equation (20), the momentum profile will be w (y,s) = 1 1+h where We can't find the inverse transform analytically in the complex transform spaces for several real-world applications. To produce the inverse Laplace, authors have also employed specific numerical techniques. Many scholars have successfully applied Stehfest's technique and Tzou's [58] for the numerical Laplace method, to solve the fractional differential equations, efficiently and effectively. The formula for the Stehfest and Tzou algorithms is as follows, respectively: where N is a positive integer.

Limiting Case:
If, in the attained solution of the momentum field Equation (21), the physical quantities Gm = Sc = K = 0, θ 1 = π 2 become zero, then The above-attained solution of the velocity profile is quite similar to the attained solution by Saqib et al. [45]. Furthermore, the Prabhakar fractional constraints β, γ can be taken as zero for the generalized results of the governed equations.

Special Cases:
The following are essential exceptional cases, with a technical significance that is widely recognized and well defined in the existing literature to understand the problem better.

Case 1: f (t) = Sin(ωt)
In the first particular case, consider the function f (t) = Sin(ωt), in which ω signifies the amplitude of the oscillating plate, then the momentum profile for this case will be w (y,s) = 1 1+h

Case 2: f (t) = Cos(ωt)
In the second particular case, consider the function f (t) = Cos(ωt), then the momentum profile for this case will be

Case 3: f (t) = te t
In the third particular case, consider the function f (t) = te t having its Laplace F(s) = 1 (s−1) 2 , then the momentum field for this case will be w (y,s) = 1 1+h

Validity of the Fractional Model
Both numerical approaches, Stehfest's and Tzou's, were compared by drawing Figure 2a,b for the concentration and temperature profiles. There is a slight overlap of the findings between the two curves. Figure 3a,b compares both numerical methods and solutions for the velocity field, using the Prabhakar fractional methodology with Saqib et al. [42]. The simulations obtained by employing the Prabhakar fractional model have a good accuracy, compared to Saqib et al.'s study [42].

Results with Discussion
The mixed convection Brinkman-type nanofluid model's thermal study, due to an oscillating inclined plate in an applied magnetic field, is examined using fractional simulations. Different nanoparticles, i.e., molybdenum disulfide ( ) and graphene-oxide ( ) are utilized for the nanofluid with a base engine oil fluid. A close thermal inspection of the modified nanofluid model is provided using molybdenum disulfide and graphene oxide nanoparticles with a viscoelastic base fluid. In the preceding part, the Prabhakar fractional derivative framework was successfully executed for the solution of the governed fractional equations. Some exceptional cases for the momentum profile are also described, improving the flowing fluid's physical significance. The physical impact of the different constraints on the thermal, concentration, and momentum profiles are also analyzed graphically in Figures 4-11, with different nanoparticles and different ranges, such as 0.1 < ', !, • < 0.9, 4 < g < 16, 0.01 < F < 0.04, 0.5 <^< 2.0, 2.0 < • < 10, 0.5 < b < 2.0, 0.5 < < 2.0 and 0.1 < j < 1.5. This section outlined the thermal dynamics of the flow model when the flow parameters were varied.

Results with Discussion
The mixed convection Brinkman-type nanofluid model's thermal study, due to an oscillating inclined plate in an applied magnetic field, is examined using fractional simulations. Different nanoparticles, i.e., molybdenum disulfide (MOS 2 ) and graphene-oxide (GO) are utilized for the nanofluid with a base engine oil fluid. A close thermal inspection of the modified nanofluid model is provided using molybdenum disulfide and graphene oxide nanoparticles with a viscoelastic base fluid. In the preceding part, the Prabhakar fractional derivative framework was successfully executed for the solution of the governed fractional equations. Some exceptional cases for the momentum profile are also described, improving the flowing fluid's physical significance. The physical impact of the different constraints on the thermal, concentration, and momentum profiles are also analyzed graphically in Figures 4-11, with different nanoparticles and different ranges, such as 0.1 < α, β, γ < 0.9, 4 < Gr < 16, 0.01 < ϕ < 0.04, 0.5 < ω < 2.0, 2.0 < Gm < 10, 0.5 < K < 2.0, 0.5 < M < 2.0 and 0.1 < Sc < 1.5. This section outlined the thermal dynamics of the flow model when the flow parameters were varied. nanofluid with the increasing M. Similarly, the inclination of a magnetic field reduces the influence of M on the Lorentz force. The Lorentz force has the most significant impact at = Ï/2 (normal magnetic field), as seen in the figures. While the slip parameter ℎ improved, the fluid motion and momentum field increased with the enhancement in the slip factor, as predicted in the figures. Figure 11a,b predicts the comparative change in the flow rate with the interaction of all types of nanoparticles º)ˆ with the engine oil as based liquid and the effects of time on the momentum field. It is interesting to note that the interaction of the nanoparticles are more practical in improving the flowing fluid rate, as compared to the nanoparticles. The addition of the different nanoparticles means an enhancement in the fluid density, increasing the flowing nanofluid's thickness. Therefore, the ( ) nanoparticle-based suspension of the nanofluid has a lower flowing rate, as compared to the ( ) nanoparticles mixed with the nanofluid. Furthermore, as illustrated in Figure 11b, the thickness of the boundary layer grows thicker than the thickness of the momentum boundary layer with time, thus enhancing the fluid motion.
Finally, the numerical analysis of the temperature, concentration, and velocity profiles by Stehfest's and Tzou's numerical schemes, is analyzed in Table 3. To validate our numerical inverse Laplace, transform solutions using Stehfest's and Tzou's methods, we performed numerical calculations for the temperature and velocity fields, which were found to be in good agreement, as shown in Table 3. The temperature profile, concentration, and momentum fields declined as they moved away from the inclined plate with the increasing 7 and became zero for 7 → ∞, which signifies the boundary conditions. Furthermore, in the numerical analysis of both Stehfest's and Tzou's algorithms, very closed values of all governed parameters by the different numerical techniques, also validate this study. At instant rises, the decreasing change in the heat transmission and mass transfer is progressive. According to Table 4, the local Nusselt number and the wall shear force decreased as it increased, indicating a significant agreement for the temperature profile. Table 5 also shows the percentage difference in the momentum field analysis.                       By changing (α, β, γ), the temperature profile for both MOS 2 − EO and (GO − EO) nanomaterial suspensions showed a declining tendency. The thermal profile nearer the plate is maximal for the lower values of the fractional parameters, while the most significant fall in the profile is seen for the bigger values. The higher values of (α, β, γ) related to the thickness of the thermal and momentum boundary layers represent the physical point. The higher values of reducing (α, β, γ) the thickness of the thermal boundary layer, result in a declining trend in the temperature and velocity profiles. Moreover, the improvement in thermal rate, due to MOS 2 − EO is more progressive, as compared to GO − EO, due to the physical characteristics of the selected nanoparticles. Figure 4b conveys the effects of Pr, showing a declining change in the temperature field. Mathematically, as Pr is the ratio of the kinematic viscosity to the thermal diffusivity and, physically, varying the value of the Prandtl number, the thermal conductivity reduces, and the fluid becomes thicker, due to which the temperature reduces, as depicted in Figure 4b. Figure 5a,b indicates the impact of the fractional constraints α, β, γ, and the Schmidt number Sc on the boundary layer concentration for both types of nanofluids (MOS 2 − EO) and (GO − EO). Again, similar to the thermal profile, the concentration field also decays by increasing the fractional constraints and Schmidt number values. Physically, the molecular diffusivity will decrease with the Schmidt number values' enhancement. Furthermore, with the thermal field, the concentration profile also declares more values in the graphical representation for the molybdenum disulfide (MOS 2 ) based nanofluid, as compared to the graphene oxide (GO) mixed suspension. Figure 6a,b highlights the physical impact of the fractional constraints α, β, γ and the volume fraction on the momentum field for both the molybdenum disulfide-engine oil (MOS 2 − EO) and graphene oxide-engine oil (GO − EO) type nanofluids. Increasing the values of (α, β, γ), lowers the mobility of the EO-based nanofluid, due to a decrease in the thickness of the momentum barrier layer. Assigning the variation in the fractional parameters, the fluid motion slows down. Moreover, with the variation in the volume fraction parameters, the velocity profile is improved, due to the effective density of the considered nanoparticles. The fluid velocity is also improved for the (MOS 2 ) based nanoparticles, as compared to the (GO) suspension, due to the heat conduction. As the volume fraction increases, the fluid becomes more viscous, indicating that the GO nanoparticles' fluid velocity will decrease, owing to the heat conduction. These novel observations may present many applications in improving the thermal capacitance of various engineering and industrial processes. Figure 7a,b predicts the influence of the Prandtl number Pr and the Schmidt number Sc on the momentum profile of the Brinkman-type nanofluid. The flowing rate is observed as the decreasing value for both parameters Pr and Sc. As the enhancement in Pr creates some more hinderances and increases the flowing fluid viscosity, it results in a decrement in the fluid velocity, and the fluid flows more slowly. Similarly, similar to the impact of the Prandtl number, the Schmidt number also slows down the fluid motion, due to a decrement in the diffusivity rate of the flowing fluid. It is concluded that the velocity profile of the nanoparticles (GO) decreases while approaching the plate, but the velocity of the nanoparticles (MoS 2 ) rises owing to the nanoparticle features.
The graphical impact of the porosity parameter K and the effects of the oscillations of the inclined plate, is represented in Figure 8a,b. For the general function f (t), we have considered the function Sin(ωt) as a particular function, where ω represents the frequency of the oscillations. Moreover, this can be seen that the fluid flow moves fast with the increment in the frequency of the oscillations. Figure 9a,b depicts the graphical behavior of the heat Grashof number Gr and the mass Grashof number Gm on a velocity profile. The velocity profile has been enhanced with large values of Gr and Gm. Physically, the heat Grashof number Gr creates natural convection, owing to the buoyant force, and an increase in the Gm improves the buoyancy forces, increasing the fluid velocity. This is the weighted average of the buoyancy and viscous forces. A higher Gr generates an increase in the buoyant forces, which causes the generated flows to increase. Figure 10a,b depicts the impact of the slip parameter h and the magnetic number M on the fluid velocity. The velocity field was reduced by raising the magnetic parameter value. Physically, an increase in the magnetic field results in an increase in the Lorentz force, which decreases the velocity of the fluid. Furthermore, the maximum effect of the Lorentz forces is when the angle of inclination of the applied magnetics is at the right angle to the oscillating plate. M is a dimensionless number associated with the Lorentz force, which opposes the nanofluid velocity. The more excellent the M, the greater the Lorentz force, which resists motion. As a result, the velocity was slowed in both cases of the engine oil-based nanofluid with the increasing M. Similarly, the inclination of a magnetic field reduces the influence of M on the Lorentz force. The Lorentz force has the most significant impact at θ 1 = π/2 (normal magnetic field), as seen in the figures. While the slip parameter h improved, the fluid motion and momentum field increased with the enhancement in the slip factor, as predicted in the figures. Figure 11a,b predicts the comparative change in the flow rate with the interaction of all types of nanoparticles MOS 2 and Go with the engine oil as based liquid and the effects of time t on the momentum field. It is interesting to note that the interaction of the MOS 2 nanoparticles are more practical in improving the flowing fluid rate, as compared to the GO nanoparticles. The addition of the different nanoparticles means an enhancement in the fluid density, increasing the flowing nanofluid's thickness. Therefore, the (GO) nanoparticle-based suspension of the nanofluid has a lower flowing rate, as compared to the (MoS 2 ) nanoparticles mixed with the nanofluid. Furthermore, as illustrated in Figure 11b, the thickness of the boundary layer grows thicker than the thickness of the momentum boundary layer with time, thus enhancing the fluid motion.
Finally, the numerical analysis of the temperature, concentration, and velocity profiles by Stehfest's and Tzou's numerical schemes, is analyzed in Table 3. To validate our numerical inverse Laplace, transform solutions using Stehfest's and Tzou's methods, we performed numerical calculations for the temperature and velocity fields, which were found to be in good agreement, as shown in Table 3. The temperature profile, concentration, and momentum fields declined as they moved away from the inclined plate with the increasing y and became zero for y → ∞ , which signifies the boundary conditions. Furthermore, in the numerical analysis of both Stehfest's and Tzou's algorithms, very closed values of all governed parameters by the different numerical techniques, also validate this study. At instant rises, the decreasing change in the heat transmission and mass transfer is progressive. According to Table 4, the local Nusselt number and the wall shear force decreased as it increased, indicating a significant agreement for the temperature profile. Table 5 also shows the percentage difference in the momentum field analysis.

Conclusions
The fractional-supported Brinkman-type nanofluid was examined using different nanoparticles and an oscillating inclined poured plate. Molybdenum disulfide and graphene oxide nanoparticles with a viscoelastic base liquid are investigated to increase the thermal performance. The motivations for implementing the Prabhakar-like fractional model are justified, as this fractional algorithm contains modern definitions without any restriction of singularities. The correctness of the fractional scheme is checked by setting up a comparison task with another fractional approach, and several unique and limiting circumstances are given to obtain the physical insight into the examined fractional model. Finally, the thermal determination of the used molybdenum-disulfide (MoS 2 ) and graphene-oxide (GO) nanoparticles are depicted in comparison. The following are the most critical findings from the present model:

•
The fractional parameters α, β, γ declined the thermal, concentration, and momentum profiles for the MoS 2 − EO and GO − EO nanoparticle suspension and the rate decrement are thicker as the time instants increase.

•
The heat transfer rate for both nanoparticles susceptible to the engine oil base material may be regulated for the specific Prandtl number values.

•
Because of the thermal conductivity factor, the thermal performances of the molybdenum disulfide MoS 2 nanoparticles with engine oil base fluid are more progressive than the graphene oxide GO nanoparticles.

•
The changing velocity is caused by the variations in the heat and mass Grashof constants for the molybdenum disulfide engine oil and the graphene oxide engine oil suspensions. • Using the Prabhakar operators with the fractional coefficient's parameter settings, might help to select an appropriate computational formula that produces a good consistency between the experimental and theoretical values.

•
In the graphical comparison of the attained solution of the momentum profile with Saqib et al. [42], the overlapping of both curves validates the attained results of this study. • Implementing the Prabhakar fractional derivative technique is also a powerful tool to simulate the analytical expressions for any coupled and nonlinear partial differential system.

•
The patterns and characteristics of all physical flow metrics coincide perfectly with the published studies.

•
The overlaying of both curves in assessing both numerical techniques confirms the obtained solutions of the governed equations.
Our computational effort has effectively clarified the parametric implications of the flow of two different nanoparticle phases. The results for the two different nanoparticles show significant monotonic differences. This study might be expanded to three or more hybrid phases to identify which hybrid phase is the most effective.