New Solutions of Fractional Jeffrey Fluid with Ternary Nanoparticles Approach

The existing work deals with the Jeffrey fluid having an unsteady flow, which is moving along a vertical plate. A fractional model with ternary, hybrid, and nanoparticles is obtained. Using suitable dimensionless parameters, the equations for energy, momentum, and Fourier’s law were converted into non-dimensional equations. In order to obtain a fractional model, a fractional operator known as the Prabhakar operator is used. To find a generalized solution for temperature as well as a velocity field, the Laplace transform is used. With the help of graphs, the impact of various parameters on velocity as well as temperature distribution is obtained. As a result, it is noted that ternary nanoparticles approach can be used to increase the temperature than the results obtained in the recent existing literature. The obtained solutions are also useful in the sense of choosing base fluids (water, kerosene and engine oil) for nanoparticles to achieved the desired results. Further, by finding the specific value of fractional parameters, the thermal and boundary layers can be controlled for different times. Such a fractional approach is very helpful in handling the experimental data by using theoretical information. Moreover, the rate of heat transfer for ternary nanoparticles is greater in comparison to hybrid and mono nanoparticles. For large values of fractional parameters, the rate of heat transfer decreases while skin friction increases. Finally, the present results are the improvement of the results that have already been published recently in the existing literature. Fractional calculus enables us to control the boundary layers as well as rate of heat transfer and skin friction for finding suitable values of fractional parameters. This approach can be very helpful in electronic devices and industrial heat management system.


Introduction
The fluid is a type of matter that is constantly distorted when an insignificant quantity of force is applied outwardly [1]. Non-Newtonian fluids in comparison to Newtonian fluids have become more suitable for scientific as well as technical disciplines. Non-Newtonian fluids have created a lot of interest in the field of engineering as well as in science because of their massive range of industrial applications, food processing, lubricant performance and many more. Fluids not obeying Newton's law are termed as non-Newtonian fluids. In order to define the viscous behavior of these kinds of fluids, numerous methods were suggested. A subclass of these fluids, specifically, Jeffrey fluid, attained much interest in recent years due to the fact that the constitutive equation for Jeffrey can be condensed to that of the Newtonian model. The effects of retardation and also relaxation are demonstrated with the help of the Jeffrey fluid [2][3][4]. Jeffrey is not only a basic theoretical idea, but is also used to explain various practical problems [5]. In peristalsis, Jeffrey fluid has been also nanofluids are improved heat transfer fluids in experimental research [30]. The changing aspects of ternary hybrid nanofluids fit the precise explanation of Newtonian fluid. Due to increased temperature, the density of ternary hybrid nanofluids decreases linearly [31]. For industrial use, ternary hybrid nanofluids are vital to select. Scientists appreciate how the thermophysical properties and rheological behavior are affected by various parameters like temperature, mixture ratio, and volume fraction of such nanofluids [32]. Saleem et al. inspected the changing aspects of ternary-hybrid nanofluids with an accent on the variation of heat, friction, and mass transfer rates crossways the domain. The growth of species that shape water-based ternary-hybrid nanofluids is declining because of the rising unpredictable motion and thermomigration of diverse nanoparticles [33]. In thermal systems, the use of ternary hybrid nanoparticles have is the least acknowledged and all the earlier numerical and experimental research mainly focused on binary or mono nanofluids [34]. In this paper, Yb-doped and pure and ZnO [Zn1-xYbxO (x = 0, 0.01, 0.03, 0.05)] gels were set via the sol-gel procedure and then they were screen-printed on a glass substrate; as a result, UV visible spectroscopy (300-800 nm) proves that ZnO film demonstrates practical transparency in the visible range, having a supreme transmission of 60 percent, and the band gap diminishes from 3.29 to 3.22 eV by using Yb-doping [35]. In the previous two years, academicians have focused on the nanofluid squeezing flow in between the parallel plates. The applications of nanofluids have seen the attention shifted to the hybrid nanofluid, but still, there is so much to be investigated. This examination is performed to explore the symmetry of MHD squeezing hybrid nanofluid (MoS 2 -SiO 2 /H 2 O-C 2 H 6 O 2 ) flow and nanofluid (MoS 2 /H 2 O) flow between the parallel plates. The consequence shows that at the bottom of the plate, the Nusselt number is greater for the (MoS 2 -SiO 2 /H 2 O-C 2 H 6 O 2 ) hybrid nanofluid as compared to the (MoS 2 /H 2 O) nanofluid flow [36].
The Falkner-Skan problem for stretching is comprehensive for nanoparticle aggregation effects. The model is made in the occurrence of injection/suction effects, thermal radiation, and magnetic field. It is seen that the effect of nanoparticle aggregation is durable even when the concentration of nanoparticles is less. The TiO 2 /EG nanofluid heat transmission rate is found to be higher with nanoparticle aggregation effects in contrast to when it is absent [37]. In this, magnetic ferrite nanoparticles were manufactured and glazed with silica ( f errite@SiO 2 NPs) by means of a sol-gel method. Afterward, silica propylmethylimidazolium chloride ionic liquid [Sipmim]Cl was equipped and related with the above-equipped ferrite@SiO 2 NPs to manufacture ferrite silica propylmethylimidazolium chloride [Fesipmim]Cl catalyst [38]. The foremost focus of this examination is to reveal the mathematical modeling of chemically reactive Maxwell nanofluid for axisymmetric flow case by seeing Cattaneo-Christo heat flux model and reviewed nanofluid model. The result shows that associated boundary layer thickness and temperature produces due to the enrichment of the heat source and temperature ratio parameter. Moreover, advanced approximation of radiation and magnetic parameter shows devaluation of wall heat transfer. Additionally, the concentration of nanoparticles is disparaged for raising the values of reaction rate parameters [39]. Azam et al. presented a mathematical demonstration and numerical solution for unsteady Bioconvection flow having chemically sensitive Sutterby nanofluid under the impact of nonlinear radiation and gyrotactic microorganism. The results obtained are very interesting, as it is noted that the microorganism field is elevated for advanced estimation of the microorganism alteration parameter and Peclet number. Furthermore, the greater bioconvection Rayleigh number disparages the fluid velocity [40]. Azam et al. also devise novel mathematical models for an unsteady bioconvection flow of chemically responsive Casson nanoliquid with gyrotactic microorganisms and non-linear for axisymmetric flow case. Additionally, viscous dissipation ramifications as well as activation energy are observed. It is observed that the density number for microorganisms is enounced for advanced approximation of the microorganism difference parameter and Peclet number [41]. Azam et al. proposed a model of transient bioconvection of Maxwell nanofluid under the inferences of variable thermal conductivity and viscosity for the case of axisymmetric flow and gyrotactic microorganisms. It is found that the microorganism density number is wired because of the rising bioconvection Lewis number and Peclet number, but the opposite behavior is dejected for a strengthened Maxwell parameter [42]. Azam et al. investigate numerical exploration concerning the bioconvection ramifications and activation energy in the growth of thermal extrusion of a chemically delicate Casson nanoliquid over a holey cylinder under the consequences of viscous dissipation and gyrotactic microorganism. The consequences of this study disclose that the density number of microorganisms is deplored for advanced approximation of the microorganism difference parameter and Peclet number. Additionally, wall heat flux is belittled due to the adulation of the Rayleigh number [43]. In light of the above recent literature, researchers have examined applications of different types of nanoparticles and found the solutions for different geometries and fluids of Newtonian and non-Newtonian nature. In the above mentioned literature, researchers have ignored the fractional approach of such models. The main advantage of the fractional approach can also help us to improve the fluid properties for various values of fractional parameters and times. Usually, the mentioned models in the introduction are represented in terms of traditional partial differential equations (PDEs), and these PDEs are unable to decode the complex behavior of physical flow parameters and memory effects. The classical calculus only measures the instant rate of change of the output, when the input level changes. Therefore, it is not able to include the previous state of the system called memory effect. On the other hand, in fractional calculus, the rate of change is affected by all the points of the considered interval, so it is able to incorporate the previous history/memory effects of any system. Fractional calculus enables us to control the boundary layers as well as the rate of heat transfer and skin friction for finding suitable value of fractional parameter. This approach can be very helpful in electronic devices and industrial heat management systems. The thermophysical properties of nanofluids are given in Table 1. Table 1. Thermodynamics properties of nanomaterials.

Research questions
How to do fractional modeling in heat transfer? How is the rate of heat transfer of ternary nanoparticles greater than hybrid and mono nanoparticles, respectively? What is the advantage of choosing different base fluids for three types of nanoparticles?

Mathematical Formulation
Let us consider the unsteady flow of Jeffrey fluid containing ternary nanoparticles over an infinite vertical plate situated at y = 0 along the x-direction as shown in the Figure 1. At the beginning, both the plate and the nanofluid are stationary with constant wall temperature T ∞ . When t > 0, the plate remains at rest while the temperature is raised to T w . Under the Boussinesq approximation, temperature and velocity are the functions of y and t only, and the governing equations for Jeffrey fluid containing ternary nanoparticles are given as [2].
where S is the extra stress tensor, λ 1 and λ are material parameters and -pI is indeterminate part of the Jeffrey fluid.

Rivlin-Ericksen tensor is
The velocity, stress, and temperature fields for one-dimensional and unidirectional flow is T =T(t,ȳ), i indicates unit vector in x-direction and the stress field satisfying the constraints S(ȳ, 0). From this, we obtain S xx = S xz = S zx = 0, S yy = S yz = 0, S zz = S zy = 0 and where S xy indicates non-trivial tangential stress. By means of Equations (1)-(5) and Boussinesq approximation equation of momentum, energy, and also Fourier's law are found in [13].
Generalized Fourier law [2].q The boundary constraints for temperature and velocity are as follows: where ρ mn f , µ mn f , (ρβ T ) mn f and (ρC p ) mn f and k mn f are density, viscosity, volumetric expansion, specific heat transfer and thermal conductivity for ternary nanoparticles. The properties of nanoparticles are defined in Table 2 [12]. Table 2. Properties of mono, hybrid, and ternary nanoparticles [12].

Mono nanoparticles
In order to obtain the problem free from the flow region, some dimensionless variables are introduced from Equation (12): into (6)- (11), we have dimensionless governing equations: where Pr = Dimensionless constraints are as follows: .

Solution of Temperature with Ternary Nanoparticles
Applying Laplace transform on Equations (14) and (15) with boundary conditions of (16) and using Prabhakar fractional derivative. The temperature field is Introducing Equation (20) into (19): The general solution of Equation (22), which satisfies the boundary constraints given in Equation (21), is To obtain inverse Laplace transform, Equation (23) can be written in the following form by using formula (1 − x) −p = ∑ ∞ q=0 x q ( Γ(p+q) q! (Γp)) [2].
The inverse Laplace of Equation (24) by using formula L −1 1

Solution of Temperature and Velocity with Hybrid Nanoparticles
By following the same procedure of Section 3 solution of temperature with hybrid nanoparticles is given as The inverse Laplace of Equation (31) is . Solution of velocity with hybrid nanoparticles is given as where v, c 0 , c 1 , c 2 and c 3 are defined as v(s) = Prsc 0 Inverse Laplace of (32) can be obtained numerically by using Tzou's and Stehfest's algorithms [44,45].

Solution of Temperature and Velocity with Mono Nanoparticles
By following same procedure of Sections 3 and 4 the temperature equation for mono nanoparticles is as follows: In order to get the inverse Laplace transform of Equation (33), it can be written in the following form by using formula (1 − x) −p = ∑ ∞ q=0 x q ( Γ(p+q) q! (Γp)) [2]. The inverse Laplace of Equation (33) is . Velocity field for mono nanoparticles is as follows: where w, d 0 , d 1 , d 2 and d 3 are defined as w(s) = Prsd 0 The inverse Laplace of Equation (34) can be obtained analytically by using Tzou's and Stehfest's algorithms [44,45].

Nusselt Number for Ternary Nanoparticles
The heat transfer rate in terms of the Nusselt number for ternary nanoparticles is given by the following relation (Table 3). where

Nusselt Number for Hybrid Nanoparticles
The heat transfer rate in terms of Nusselt number for hybrid nanoparticles is given by the following relation (Table 3). where

Nusselt Number for Mono Nanoparticles
The heat transfer rate in terms of the Nusselt number for mono nanoparticles is given in the following relation (Table 3). where

Skin Friction for Ternary Nanoparticles
The skin friction for ternary nanaoparticles is given in the following relation (Table 4).

Skin Friction for Hybrid Nanoparticles
The skin friction relation for hybrid nanaoparticles is given by ( Table 4).

Skin Friction for Mono Nanoparticles
The skin friction relation for mono nanoparticels is given by ( Table 4). where

Graphical Results and Discussion
The prevailing work deals with the Jeffrey fluid moving along a vertical plate. Laplace transform is applied to obtain an analytical solution. The graphical investigation is made by means of Mathcad software. Figure 2 shows the effect of fractional parameters on the temperature profile. It can be seen that by increasing the values of fractional parameters α, β, γ, the temperature is decreasing. Physically, the fractional solution behaves as dual for small and large times for different values of fractional parameters. It is evident that for large values of time, the fluid temperature decreases. Therefore, we have a choice to achieve the desired results of the fractional model.  Figure 3 is showing the influence of the concentration of nanoparticles φ 1 , φ 2 , φ 3 on the temperature field. It is observed that with increasing the concentration of nanoparticles that the temperature is also increasing. Physically, increasing the concentration of nanoparticles leads to increase thermal conductivity and ultimately the thickness of thermal conductivity and hence temperature increases.  Figure 4 is showing the behavior of Pr with respect to the temperature field. When Pr values increase, temperature distribution declines. Pr number is a ratio between momentum diffusivity to heat diffusivity. Fluids with inferior Pr number values have greater thermal conductivity, which causes heat diffuse far from the heated surface more swiftly and quicker in comparison to increased values of Pr. Consequently, Figure 4 represents the ternary nanoparticles by choosing different base fluids, for example water, kerosene, and engine oil, respectively. It is very clear that by taking engine oil as the base fluid, the temperature decreases due to high values of Prandtl number.  Figure 5 shows the impact of fractional parameters on the velocity profile. It can be seen that by increasing the values of fractional parameters, the velocity profile increases. With time, the boundary layer thickness increases and the velocity is extreme in that area of the plate. Physically, the fractional solution behaves dual for small and large times for different values of fractional parameters. It is evident that for large values of time, the fluid velocity increases. Therefore, we have a choice to achieve the desired results of the fractional model. Figure 5. u(y, t) with α, β, γ when t = 1, a = 0.1, λ = 0.5, λ 2 = 0.9, Pr = 6.2, β 1 = 0.9, φ 1 = φ 2 = φ 3 = 0.05. Figure 6 is intended to show the effect of the concentration of nanoparticles on the velocity field. It can be seen that by increasing the concentration of nanoparticles velocity profile is decreasing. Figure 6 depicted the impact of the concentration of ternary nanoparticles on velocity. It is found that larger the values of φ 1 , φ 2 and φ 3 stronger the viscous force makes the boundary layer denser and hence decreases. Figure 6. u(y, t) with φ 1 , φ 2 , φ 3 when t = 1, a = 0.01, λ = 0.5, λ 2 = 0.9, Pr = 6.2, β 1 = 0.9, Figure 7 shows the effect of Pr on the velocity field. By raising Pr values, the velocity distribution is decreasing. Due to the difference between the thermal conductivities of different base fluids, the thermal conductivity of water is comparatively higher than kerosene and engine oil. Velocity became higher, and a more distinctive peak is observed near the plate. Therefore, an upsurge in the thickness of the boundary layer condenses velocity distribution.     . Comparison graph of temperature between present work and Basit et al. [2] when t = 0.8, a = 0.9, λ = 0.5, Pr = 6.2, β 1 = 0.9. Figure 10 shows that temperature distribution declines in the presence of nanoparticles. It is found that the presence of ternary nanoparticles can be used to enhance the temperature.      Figure 12 is also the comparison of velocity with and without the presence of nanoparticles. It is found that velocity in the absence of ternary nanoparticles velocity is higher. This is due to the fact that fluid becomes less thick and the distance between boundary layers increases and boosts up the flow. To see the influence of fractional parameters and time on the rate of heat transfer and skin friction, we have constructed Tables 3 and 4. From Table 3, it can be seen that the heat transfer rate in terms of the Nusselt number is decreasing by enhancing fractional parameter values and increasing with respect to time. On further notice, it is observed that heat transfer for ternary nanoparticles is greater as compared to hybrid and mono nanoparticles. Therefore, this approach can be used in practical situations, especially in heat management systems and electronic appliances. Similarly, from Table 4, we can see that the effect of fractional parameter α on the skin friction is estimated numerically and indicates that the skin friction is an increasing function when values of α are increased. Additionally, it is observed that skin friction for ternary nanoparticles is greater in comparison to hybrid and mono nanoparticles, respectively.

Conclusions
In this paper, the unsteady flow of a Jeffrey fluid moving along a vertically placed plate is resolved analytically with the help of the Laplace transform. Effects of fractional parameters, the concentration of nanoparticles, Prandtl number Pr and thermal Grashof number Gr on temperature as well as velocity distribution are obtained graphically. The main results of this paper are as follows: 1.
The fractional and ternary nanoparticles approach is helpful in improving the exiting results.

2.
Choosing a specific value of fractional parameters, thermal and momentum boundary layers can be controlled for different time. Such a fractional approach is very helpful in handling the experimental data by using theoretical information. 3.
The rate of heat transfer for ternary nanoparticles is greater in comparison of hybrid and mono nanoparticles, respectively. 4.
The obtained results are also useful in the sense of choosing base fluids (water, kerosene and Engin oil) for three types of naoparticles to achieved the desired targets in any transport problem.

5.
The new approach with ternary nanoparticles is more reliable and effective, and was validated with the existing literature. Funding: This research received no external funding.
Data Availability Statement: All the data are given within the article.

Acknowledgments:
The authors extend their appreciation to the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University for supporting this work.

Conflicts of Interest:
The authors declare no conflict of interest.