Computational Analysis of Heat Transfer Intensification of Fractional Viscoelastic Hybrid Nanofluids

In the current article, we have performed computational analysis on convection heat transfer of a hybrid nanofluid in occurrences where porous media and the effect of magnetic force are involved. In order to assess the time-fractional derivatives, Caputo’s notion is utilized while the Darcy–Forchheimer model is applied due to the involvement of the porous medium. Moreover, the boundary conditions are assumed to be nonuniform through the equilibrium between the surface tension and shear stress over a semi-infinite permeable flat surface. Keeping in view the complexity of the fractional derivative model and nonuniform boundary conditions, the problem in question is a complicated one. Accordingly, the coupled momentum and energy equation is linearized and the finite difference scheme is then applied and implemented in MATLAB Code R2020b. Furthermore, we have also offered a comprehensive analysis regarding error and convergence of the proposed numerical method. (e newly introduced numerical technique to determine the numerical solutions and some unique and interesting deductions are established. From the computational results, one can conclude that the fluid motion in both hybrid and single nanofluids slows down due to magnetic field, porosity, and inertia coefficient as the magnetic and electric fields are synchronized due to the formation of the Lorentz force and viscous interference. We believe that our proposed numerical technique regarding employment of the fractional model for heat transfer application to the hybrid nanofluid over a semi-infinite nonuniform permeable surface along with variable heat flux is not found in the literature so far. Furthermore, the obtained results will be a valuable addition to fractional calculus from an engineering point of view.


Introduction
It has always remained one of the foremost objectives of researchers to formulate innovative methods to enhance the heat transfer rate in thermal systems, which will resultantly improve the productivity of heating systems while cutting off the expense simultaneously. Owing to a wide range of applications in various industries, convective heat transfer analysis inside a porous medium is one of the most appealing subjects. e employment of the porous medium in thermal systems results in two main advantages. First, the heat exchange rate is improved due to the enlargement of the heat transfer area. Second, the flow mixing is considerably increased owing to the chaotic movement of the fluid through appropriate pores. Accordingly, porous media have attracted a significant number of researchers to employ the same in heat transfer equipment in the recent past [1][2][3][4][5].
us, the characteristics of improved heat transfer can be referred to as facts from the literature. A porous medium is defined as a medium comprising a solid matrix with interconnected vacuums like metal foams and rocks. Figure 1 can be referred to have an illustration of carbonbased porous materials. Alibakhsh et al. [6], reviewed the nanofluid flow and heat transfer through porous media. Although the tendency of physical quantities in the scale of pores is not methodical, these are assessed in engineering applications in considerably large areas. us, a good number of researchers have analyzed various surfaces partially filled with the porous medium, and they established that maximum heat transfer could be obtained without filling the channel with the porous medium [2]. Without doubt, convection heat transfer in porous media is a well-known and imperative area of research. Nevertheless, porous media most often have excellent thermal efficiency. However, the weak thermal conductivity of fluids for heat transport causes severe limitations in enhancing heat transfer in cases where fluid flow is involved. From [7][8][9], it is evident that the addition of solid nanoparticles to conventional fluids results in a visible increase in thermal conductivity.
us, instantaneous engagement of nanomaterials (i.e., suspensions of nanostructures) and permeable amalgamations may be referred as an effective technique to achieve substantial heat transfer. e scientific developments in the discipline of nanofluids lead to significant progress in various applications. Based on scientific research findings, nanofluids enhance the thermal conductivity and conductive heat transfer in relation to the base liquids.One cannot deny the fact that this modern kind of fluids received a remarkable acceptance from industry as these are deemed to be a good alternative of many traditional fluids, heat exchangers, and heat transport liquids [10] since nanofluids possess good temperature efficiency over liquids so that they can be utilized in numerous heating and ventilation applications (Figure 2(a)). e issue of water consumption in huge quantities and the production of excessive industrial waste can be addressed by increasing the heat transfer and simultaneously reducing water consumption.
Keeping in view the ultimate objective of enhancing convective heat transfer in different appliances, researchers have applied and tested multiple techniques like geometry modification, alteration of the boundary conditions, and enhancement of the thermophysical characteristics of the fluid in the recent past. e researchers and scientists have employed the technique of diffusing tiny solid particles in conventional fluids to maximize fluids' thermal characteristics. Previously, the particles with micrometer or millimeter measurements were used for this purpose. e use of these particles led to serious problems such as fouling, degradation, and loss of elevated pressure and, therefore, was not widely considered by researchers for future investigation. With the advancement in technology, nanotechnology was introduced, which further led to the discovery of nanoparticles which were of nanometer scale. Figure 3, can be referred for an illustration of nanoparticles. e modern type of dispersions named nanofluids is obtained by suspension of solid nanoparticles in ordinary liquids. ese suspensions were named as "nanofluid" for the very first time in [11]. Choi presented a new type of fluid for heat transport containing tiny amounts of nanostructures that were metallic or nonmetallic. e particles described above were uniformly dispersed in a consistent manner and at a regular phase.
e key advancements in nanofluid technology revealed the extraordinary capability of nanofluids in heat transfer applications, which provided a ground for research in the manufacturing area and academic research. e particles in nanocomposites vary from 1 to 100 nm on average. However, new research avenues, with simultaneous pros and cons, have been generated with the introduction of nanomaterials in fluid dynamics and heat transport. Investigators of nanomaterials need to have a thorough knowledge of the collective and rheological behaviors of nanofluids. is indicates that many researchers have made attempts to enhance the thermal performance utilizing the addition of various nanofluids to several heat exchangers (Figure 2(a)). Based on the base fluid and type of nanoparticles, nanofluids are divided into different categories. Quite a few researchers have also tried to investigate the thermal conductivity by adding metal nanomaterials to basic fluids like water (H 2 O), ethylene glycol (EG), and oil [12][13][14][15]. Nonmetals and nonmetal oxides are the nanoparticles that are also used to formulate different nanofluids [16,17]. Recently, the researchers have also paid attention to carbon nanotubes [18,19]. From [20][21][22][23], it is evident that  comprehensive studies addressing the experimental, numerical, and theoretical aspects are conducted by researchers to investigate the convective heat exchange of nanocomposites in various thermal environments. Figure 2(b) offers a pictorial view of nanofluid applications. One can find several research articles wherein it is established that the existing models are not capable of predicting the thermal conductivity of nanofluids, and by the addition of nanoparticles in ordinary liquids, it can be enhanced multiple times. Here, Table 1 Figure 2: Common heat exchangers used in various industries [49] and primary applications of some nanofluids [50].
dealing with multiple kinds of nanomaterials at a time, termed as "hybrid nanofluids." e hybrid nanofluids can be considered the most modern kind of nanofluids where the main or base fluid contains more than one type of nanoparticles. More than one form of nanoparticles is utilized in such a category of nanocomposites. e type of nanofluids mentioned earlier possesses better properties when compared with mono nanofluids, and a number of researchers have paid attention to these fluids [24][25][26]. Values of the thermal conductivity of base fluids and different nanofluids are displayed in Figures 4(a) and 4(b), respectively [19,[27][28][29]. Very few research articles offer a comprehensive analysis of free convection characteristics of hybrid nanofluids in porous media. In the recent past, the study in [30] presented the numerical analysis of MHD hybrid nanofluids inside a T-shaped cavity. ey studied the thermal enhancement of hybrid nanofluids by considering the mixture of MWCNT-Fe 3 O 4 . e viscosity of the nanofluid responds to the power of the magnetic force.
e Darcy-Brinkman-Forchheimer (DBF) integrative version was utilized to obtain the governing equations based on the dimensionless variables. Natural convection heat transfer owes a wide range of applications in several fields, including, but not limited to, renewable radiation, earth sciences, and development of thermal mass and nuclear power. In particular, free convection is substantial in instruments positioned at insufficient spaces, which is a common phenomenon in technological situations. In this case, heat transfer is caused by the fluid motion occurring due to buoyancy force and no external factor contributes towards heat transfer generation.
us, the mechanisms involving free convection do not face the challenges of mechanical failures, which often occur in forced convective regimes. Moreover, free convection involves relatively low cost and produces immaterial turbulence from the shortage of established synchronous generators. Keeping in view the significance of natural convection in porous structures, several researchers have evaluated the impact of deployment of nanofluids in these media [31][32][33]. In the literature, a number of studies can be referred to that deal with the natural convection of nanomaterials in porous media with varying magnetic fields. For instance, magnetohydrodynamic (MHD) effects have attracted a number of scholars due to its widespread refrigerating capabilities of power plants and nuclear power plants, management of industrial waste, regulation of the spread of contaminants, etc. [34,35]. e mathematical field that interacts with the study of defining real or complex order of the derivatives is called fractional calculus. In the recent past, significant developments have been witnessed in this area and this is being utilized in several engineering applications [40][41][42][43][44]. Researchers made several attempts to present a precise definition of the fractional derivatives. e definitions offered by Riemann-Liouville, Grunwald-Letnikov, and Caputo are the most widely accepted ones. e authors of [45] carried out a study to offer a comparative analysis of these definitions as follows: Grunwald-Letnikov definition is Riemann-Liouville definition is Caputo definition is  Figure 3: Comparison of different particle sizes [49]. In reference to the above, the literature offers attributes of Caputo and Reimann-Liouville derivatives, implementations of the time-fractional Navier-Stokes equations, and examination of many structural mathematical problems in the CFD (computational fluid dynamics) [46,47]. e above analysis leads us to believe that the hybrid nanofluid has appeared as a promising heat transfer fluid despite having some unpredictability in the established outcomes and inadequate knowledge of the heat transfer tool in nanofluids. Researchers have extended research studies in the recent past by considering hybrid nanofluids, which are formed through suspension of unlike nanoparticles either in mixture or compound form.
is provided us sufficient motivation to carry out the current study to present the fractional partial differential equations of the mixed convective flow over a nonuniform permeable semi-infinite surface filled with hybrid nanofluids subject to an applied magnetic field with variable heat flux. Here, the Darcy-Forchheimer model is utilized for the porous media, while Caputo's notion is used to approximate time-fractional derivatives. Moreover, two models, namely, fractional shear stress and Cattaneo heat flux, are presented in the constitutive relations with nonuniform boundary conditions. us, with mixed time-space derivatives in the convection terms, the nonlinear principal equations are formed. e finite difference method, together with L 1 algorithm is used to acquire the numerical solutions of the proposed problem. With the help of graphs and tables, it is shown how the parameters of interest affect velocity and temperature fields. Furthermore, the subject parameters are investigated in detail to affirm surface tension properties and heat transfer performance. In the current study, it is observed that the fractional-order derivative viscoelastic model gives more accurate and precise results than the integer-order viscoelastic model. Keeping in view the distinctive application like guided wave testing technique, one should have comprehensive knowledge of the subject model. Furthermore, it is pertinent to mention here that the uniform suction or injection is investigated widely by a number of researchers, but any researcher rarely studies nonuniform suction/injection over a semi-infinite flat surface filled with hybrid nanofluid. In industrial applications, the nonuniform suction/injection is more repeatedly observed than the uniform ones. Poor permeability is witnessed in a pipe-type filter system due to the accumulation of impurities at the inlet of the pipeline. In thermal protection applications where the wall has nonuniform slots, the skin friction is reduced by the adverse pressure gradient inside the boundary layer. e framework of the paper is given as follows: Section 2 deals with a comprehensive overview of the flow problem. Sections 3 and 4 refer to some significant notions related to engineering parameters and nondimensionless parameters. A finite difference method for the model equations has been developed in Section 5 and a precise numerical scheme is derived accordingly. A detailed discussion regarding the convergence of our proposed numerical scheme along with error analysis of the L 2 error norm is given in Section 6. Section 7 is dedicated to candid discussion regarding results Mathematical Problems in Engineering obtained through this study. In Section 8, a precise conclusion regarding the outcomes of the current study is given.

Description of the Mathematical Problem
Here in this article, we have taken into consideration the mixed convection flow of a generalized incompressible hybrid nanofluid over a dense porous leveled surface with flexible heat flux where the velocity of the even mainstream is taken as u ∞ . According to the boundary layer theory, the viscous effects have a leading role in the momentum equation, whereas pressure gradients have a negligible effect. e fluid is infused at the surface of the plate with uneven velocity v w (x) when the impurities in the fluids have partially blocked the holes. e fluid field outside the boundary layer zone is not altered as the magnitude of v w (x) is considerably small. Moreover, changes were observed in the power-law q h (x) � Ax n and heat flux is generated and sustained at all time t > 0. At the same time, free convection occurs due to temperature gradient and the uniform temperature of the atmosphere far from the surface is T ∞ . e initial condition is supposed to be fixed with ambient temperature T ∞ . A porous medium occupies the flow domain and the fluid employed is the Al 2 O 3 − ZrO 2 /BF hybrid nanofluid, for which Table 2 gives thermophysical characteristics. e Forchheimer model presents the porous matrix, and the local thermal equilibrium model between the mixture and the medium is maintained. A square velocity term representing the nonlinear drag impact that occurs because of a porous medium is introduced in Darcy law to obtain the Forchheimer model, the first-ever non-Darcy model. Furthermore, in a restricted temperature domain, the fluid characteristics are assumed to be constant.
is model presents the momentary conduct of the heat transfer course of a stationary heat source. e fractional time derivatives α and c contribute to controlling fluid flow and heat propagation. All other physical parameters are constant, and viscous dissipation in the energy equation is almost insignificant. We have utilized the fractional Caputo time derivative to develop a mathematical model [51]. Assuming all the above conditions, the principal boundary layer equations of momentum and energy for the fractional hybrid nanofluid are as follows [25,[52][53][54]: where u and v represent the flow velocity, T and ρ hnf denote the hybrid fluid density and fluid temperature, respectively, (c p ) hnf is the hybrid nanofluid heat capacity at constant pressure, g is the acceleration, and β hnf is the thermal expansion coefficient of hybrid nanofluid, and σ xy and q are the shear stress and heat flux. In Table 3, ϕ 1 and ϕ 2 , respectively, denote the volume fractions of Al 2 O 3 and ZrO 2 nanomaterials, where ϕ 1 � ϕ 2 � 0 is the normal fluid, μ is the dynamic viscosity, ρ is the standard liquid, μ reflects the elastic viscosity, ρ is the density, c p is the specific heat at constant strain, ρc p is the heat capacity, k is the thermal conductivity, and σ is the electrical conductivity in which hnf, nf, f, and s1 and s2 subscripts represent hybrid nanofluid, nanofluid, fluid, and solid nanoparticles for Al 2 O 3 and ZrO 2 , respectively.
In the momentum equation, the constitutive relation of the nanoliquid with the Caputo fractional model is incorporated as where λ 1 , μ hnf , and α are the relaxation time for momentum equation, hybrid nanofluid viscosity, and memory term for velocity, respectively. For the heat transfer enhancement, we have added the Darcy-Forchheimer porosity medium to the momentum equation. e strain of Darcy-Forchheimer model to the liquid that passes via the porous medium is described as [60] e permeability and porosity of the porous medium are k and ψ, respectively.

Nanofluid phase
Hybrid nanofluid Heat capacity Dynamic viscosity Electrical conductivity Mathematical Problems in Engineering utilizing the latest Taylor series of fractional order, which establishes thermoelasticity with time-delay and memorydependent derivative; the fractional Cattaneo heat conduction model is [61] q + λ where λ 2 and c are relaxation time for energy equation and the control parameter for the same equation; this model removes the dilemma of the infinite speed of energy distribution in thermoelectric liquid and increases the performance of the thermoelectric substance as per the specification of permissible fractional derivative order values c . Analogously, the fractional temperature equation can be derived by incorporating equation (10) into (6) as For the above-considered equations, the required initial and boundary conditions are stated as follows:

Significant Parameters of Practical Interest
e drag coefficient and heat transfer rate are the important variables for the above models, which can be considered the measurement of the surface shear stress and heat transport ability of the convection flow system simultaneously. From the constitutive fractional equations (7) and (10), the local skin friction coefficient and the Nusselt number can be interpreted as follows:

Nondimensional Caputo Fractional Flow Equations
To nondimensionalize the above equations, the corresponding dimension variables are added to where L is the characteristic length. Governing equations (4), (9), and (11) are (for simplicity, the dimensionless mark * is omitted): together with the following flow conditions: where Mathematical Problems in Engineering e local drag factor and heat transfer rate in steady state are achieved by solving (13) at each time step, as t k approaches the convergence time: us, the average drag coefficient and heat transfer rate satisfy, respectively,

Numerical Technique
is section deals with the numerical solution of fractional diffusion equations. In order to determine the numerical solution of nonlinear equations (15) and (16) mixed with space and time derivatives, the earlier formulated finite difference method together with the L 1 algorithm is used.
e L 1 scheme developed in [62] is utilized to solve Caputo fractional derivative and is considered one of the most valuable and effective methods for discretizing one-dimensional space-time derivative.

Discretization Method.
e unsteady coupled system of nonlinear fractional time and space derivative equations (15) and (16) are solved by a well-known finite difference process (FDM) along with L 1 in nonuniform boundary condition (16). We have x i � iδx, where i � 0, 1, 2, 3, . . . , M; y j � jδy, where j � 0, 1, 2, 3, . . . , N; and t k � kτ, where k � 0, 1, 2, . . . , R, in which δx � L x /M and δy � L y /N are space steps and τ is the time step. e numerical solution of velocity and temperature at the mesh point (x i , y j , t k ) is denoted as u k ij and T k ij . To discretize the time-fractional derivative, the L 1 approach is applied to (0 < α < 1) and (0 < c < 1) as [62] where α s � (s + 1) 1− α − (s) 1− α , in which s � 0, 1, 2, . . . , r. By the backward difference, the time derivative is discretized: By linearization processing, the implicit difference scheme is considered for the space derivatives: e finite difference approximations of equation (15) using Crank-Nicolson method is defined below : en, equations (22) and (23) are substituted into the fractional derivatives: Eventually, we obtain the linear iteration scheme and equation for velocity: e associated temperature terms are approached by a similar technique of difference. Ultimately, a system of linear iteration equations is acquired: where Mathematical Problems in Engineering , where c s � (s + 1) 1− c − s 1− c , in which s � 0, 1, 2, . . . , r. Likewise, the thermal condition is computed as follows at y � 0 [54]: For simplification, equation (28) is substituted by iteration equations.

Convergence Analysis
e numerical convergence of the proposed numerical implicit method will be demonstrated as follows: let u(x i , y j , t k ) and T(x i , y j , t k ) (i � 1, 2, . . . , Mx − 1; j � 1, 2, . . . , My − 1; k � 1, 2, . . . , n) be the exact solution of the problem (15)-(16) at mesh point (x i , y j , t k ). We considered the problem We now considered the boundary conditions for u and T and negligible effect for v since convergence analysis for the above model equations is as follows: e right-hand-side terms f 1 (x, y, t) and f 2 (x, y, t) are calculated and given in Appendix (A.1). e exact solution is u(x, y, t) � 100(1 + t 2 )y(y − L y )x(x − L x ) and T(x, y, t) � tx n (1 − y)((y − L y )/L y )(x − L x ). Comparison between the exact results and the numerical solutions of u and T are presented in Figure 5 with ZrO 2 , and t � 1.0. e good agreement of results with the experimental data supports the accuracy of our proposed numerical technique.
Moreover, the L 2 norm is used to calculate the L 2 error, and analysis of Table 4 shows that the error at n x is minimized by refining mesh size; n y � 20 corresponds to the optimal error and a L 2 -norm of order 10 − 4 and 10 − 5 . e error is considerably reduced with mesh refinement. e L 2 norm of the finite difference scheme has a greater order of accuracy for almost all the refinements; however, supreme accuracy is obtained at n x and n y � 200.

Discussion
e current section is dedicated to discussing the outcomes acquired when the finite difference method is applied to the discretized momentum and energy equations. e different values of parameters of interest as required in several industrial and biomedical applications are numerically calculated by formulating a MATLAB code. It is shown through graphs followed by a detailed discussion of how the velocity and temperature fields are affected by the fractional derivative, porosity, mass convective parameter, generalized Reynold number, permeable parameter, and power-law exponent. Moreover, a graphical and tabular presentation regarding enhanced thermal conductivity of the proposed hybrid nanofluids over a flat permeable surface is also offered. e impact of nonuniform permeation on the viscoelastic hybrid nanofluid, along with mass transfer over a permeable surface, is also studied in detail. In order to have simplicity, we have fixed the temperature and velocity relaxation times as 0.1. Furthermore, the Reynolds number is assumed to have a fixed value of 0.95, the Prandtl number is 6.2, and the time is 5.0. e quantities of interest like skin friction coefficient and Nusselt number are evaluated. e deviation in the velocity of hybrid and single nanofluids against various values of the fractional parameter α is shown in Figure 6(a). A decreasing trend in the velocity of hybrid and single nanofluids is witnessed against an increment in the fractional parameter α. At a maximum time t � 5.0, the behavior of α is analyzed. It is observed that the viscosity of fluids and boundary layer thickness increase against an increment in the value of the fractional parameter, which further results in decreasing the velocity of the fluid. Moreover, a comparative analysis between Al 2 O 3 − ZrO 2 /BF (2.0 + 2.0%) and the ZrO 2 − BF (4.0%) is offered, and it is observed that the velocity of the single nanofluid ZrO 2 /BF faces less viscosity effect than that of the hybrid nanofluid. Furthermore, it is shown that an increment in the fractional parameter α for a small time (t � 0.01) significantly contributes to enhancing the velocity of the hybrid nanofluid. e impact of the velocity relaxation parameter on velocity distribution is also attributed through Figure 6(b).
is figure also reveals that corresponding to an increase in the relaxation parameter λ 1 , the velocity and boundary layer thickness tend to increase and decrease, respectively. e enhanced viscosity of the fluid prevents the occurrence of the delaying characteristic of viscoelastic fluid. Also, in our case, the fluid mixture comprising Al 2 O 3 and ZrO 2 will take maximum time to resume the normal state when the value of λ 1 is increased from 0.0 to 1.2. In addition to this, the graphical presentation gives us a comparative view to witness that hybrid nanofluids assume higher velocity than single nanofluid. Furthermore, the effect of the porosity parameter λ on the velocity profile is presented through Figure 6(c). It is observed that corresponding to higher values of the porosity parameter λ, the velocity shows a decreasing trend because the increased value of porosity offers resistance to the flow. On the other hand, corresponding to higher values of the porosity parameter λ, the boundary layer thickness tends to decrease. e current study also indicates that the viscosity of the fluids increases against an increment in the values of 0 ≤ α ≤ 1 and λ 1 , which can also restrict the velocity of fluid in the pores of the porous medium. e velocity of the hybrid nanofluid Al 2 O 3 − ZrO 2 /BF (2.0 + 2.0%) is somewhat slower than that of the single nanofluid ZrO 2 − BF (4.0%) as it possess higher density. Along similar lines, corresponding to higher values of the inertia coefficient parameter (Forchheimer parameter) (F r ), the hybrid fluid velocity shows a decreasing behavior (see Figure 6(d)), whereas the thickness of the fluid is not altered by F r . Since we have assumed free stream velocity at y ⟶ ∞, the F r has no impact. Physically, it can be inferred as a state of the porous gap where the pore size is big enough such that the porous medium is stretched with fluid-solid interaction, which further increases the viscous intervention. As better flow resistance is produced against an increase in F r , the fluid velocity of both the hybrid nanofluid Al 2 O 3 − ZrO 2 /BF (2.0 + 2.0%) and single nanofluid ZrO 2 − BF is restricted. Figure 7(a) shows the impact of the magnetic parameter Ha on the velocity profile of both fluids. Here, the Lorentz force is produced, in the presence of a transverse magnetic field, due to the magnetic field's attraction throughout the fluid's motion over a flat permeable surface. e positive values of the magnetic field parameter cause the velocity to decrease because resistive force is increased. Accordingly, the velocity is decreased in the x direction with boundary layer thickness. Upon comparison of velocities of both the nanofluids, it is observed that the resistance forces have a significant impact on the hybrid nanofluid Al 2 O 3 − ZrO 2 /BF (2.0 + 2.0%) due to its high density than that on the single nanofluid ZrO 2 /BF. Several studies suggest different ranges of the parameter like 1.0 ≤ Re ≤ 2.0 and 0.8 ≤ α ≤ 1.0. It can be noticed from Figure 7(b) that corresponding to an increment in the generalized Reynolds number (Re), the velocity shows an increasing trend while a total opposite behavior is witnessed for boundary layer thickness which also agrees with the integer-order model. Hence, we can say that the velocity of the flow can also be characterized by the magnitude of the fractional Reynolds number. Furthermore, the effect of Λ 1 on velocity fields is presented in Figure 7(c). In convection problems, the heat transfer is governed by the thermal Grashof number. It can be physically interpreted    Figure 7(d). It can be seen that change in v w illustrates the improvement in mass transfer via the permeable surface at a low pace where v w is the permeable parameter. In addition to this, it is noted that the optimum permeability of the surface is attained when the length of the permeable plate is fixed. In order to demonstrate the influence of varying heat transports on the boundary layer through the permeable surface, we computed the velocity profile horizontally. We then compared it with that of the conventional layer. Here, we are not comparing vertical velocity as the order of magnitude is not the same as the horizontal velocity. e variation in u(x, y, t) against different permeable parameters (v w ) is reflected through Figure 7(d). It also depicts the change in the horizontal direction with a fixed vertical direction. Figure 7(d) also offers a pictorial view of u(x, y, t)-profiles against various permeable parameters (v w ). e function v w (x) increases along the surface in the x direction in the case of all the curves. erefore, the mass transfer is significantly higher via the permeable surface of both the nanofluids and by an increasing x. A higher amount of fluid is inoculated into the boundary layer. How the temperature profile is affected by c is shown through Figure 8(a). A substantial temperature gradient that endorses the convection flow is presented through the smooth profiles. It is also observed that c holds a significant part in the thermal boundary layer as the temperature distribution and the thickness of the thermal boundary layer drastically increase with c. e temperature of the hybrid nanofluid Al 2 O 3 − ZrO 2 /BF (2.0 + 2.0%) is higher to some extent when compared with that of the single nanofluid ZrO 2 /BF. e thermal flow is regulated with the help of the order of fractional derivative (c), whereas the velocity flow is controlled by α. Moreover, an increment in λ 2 results in a time gap for the transfer of heat to the flow area from the heat source, and resultantly, the temperature profile tends to decrease in that specific region (see Figure 8(b)). In addition, λ 2 and heat transfer rate are directly proportional as the increase in λ 2 results in an increment in the heat transfer rate. us, the temperature profile shows a decreasing behavior corresponding to an increase in λ 2 . e effect of power-law exponent n on temperature is shown through Figure 8(c). e temperature gradient decreases against an increment in the power-law exponent, and temperature distributions of both nanofluids tend to decline. e performance of heat transfer varies at different temperatures, which conveys an evident impact of temperature. In the case of the same nanoparticle volume concentration, the temperature profiles of both the nanofluids tend to decrease when the generalized Reynolds number (Re) is increased. Consequently, the extent to which the heat transfer can be enhanced is extended. e results obtained in the current study suggest that our proposed numerical approach is efficient irrespective of the permeation function. Hence, a comparison of the three different permeation functions with the same permeable parameter (v w ) is offered for both the hybrid and single nanofluids in Figure 9(a). Furthermore, we have assumed the three parameters (v w (x)) as follows: e three different scenarios are as follows. Case I: v w (x) � v w · x, a basic function that deals with permeability at the surface; Case II: v w (x) � v w · x 2 , a quadratic function that presents a precise nonlinear permeability at the surface; Case III: v w (x) � v w · sin(Π · x/2), a well-known permeability function that is considered in a number of numerical studies performed earlier. Figure 9(a) shows the impact of the nonuniform permeable sheet on the velocity boundary layer. In real-world problems, this phenomenon occurs due to an anisotropic membrane. Moreover, the velocity distribution in the boundary layer can be altered by changing the permeability on the surface. e said technique makes the current approach different from the other classic methods where inlet velocity, the Reynolds number (Re), or the physical properties of the fluids are altered to change the velocity. In the current study, we have assumed three     different permeation functions following the identical pattern where the permeability at the right side of the surface differs from that at the left side. We have assumed the same pattern for comparison purposes. Figure 9(b) suggests that against an increment in the volume concentration of the nanoparticles, the thermal conductivity of the hybrid nanofluid shows an increasing and linear trend. It is also noticed that the hybridization of the nanofluid considerably increases the thermal conductivity. e variance in viscosity between Al 2 O 3 − ZrO 2 /BF and ZrO 2 /BF is widened with concentration of nanoparticles. Moreover, Table 5 illustrates that the thermal conductivity is being underrated by all these models when compared to the acquired correlation and the variation becomes more substantial with enhanced nanoparticle volume concentration. e table, as mentioned earlier, also offers a numerical comparison between calculated values of thermal conductivity ratios in [63] and the values projected by many other models. It can be noticed that the increment in gradient responds to an increase in nanoparticle concentration. Figure 10 presents the local skin friction coefficient and Nusselt number distributions against different values of α and Re. It can be witnessed that both the surface tension and local Nusselt number tend to increase against an increment in fractional derivative and generalized Reynold number. Furthermore, an important observation is witnessed that the surface tension shows an increasing behavior. At the same time, the velocity exhibits a decreasing trend for the value of α due to the high viscosity fluid effect. Similarly, the heat transfer rate shows a considerable increase corresponding to an increase in the generalized Reynold number Re. Objectively, it can be said that the heat transfer is increased to 74% as compared to a single nanofluid where the heat transfer rate escalates but is slower than that of hybrid nanofluid.

Conclusion
Heat efficiency and heat storage always remain as appealing subjects for researchers when it comes to renewable energy. Various fluid classes have been offered, and several types of phase-change materials have been utilized for energy storage. Nanofluids, liquids having nonmetallic particles, for instance, water, have been gaining significant attention in the recent past. In the current article, the application of fractional calculus is elaborated and fractional partial differential equations governing the mixed convection flow of hybrid nanofluids over a semi-infinite permeable surface subject to the applied magnetic field and variable heat flux are offered. e hybrid nanofluid is prepared with dispersion of nanoparticles (Al 2 O 3 and ZrO 2 ) in a viscoelastic surfactant solution, aqueous solution of cetyltrimethylammonium chloride (CTAC) and sodium salicylate (NaSal). Furthermore, we have considered a flat surface filled with a porous medium, and the notion of Caputo fractional derivatives is used. Commercial aluminum and zirconium particles, having diameters of < 50 nm and < 20 nm, respectively, were distributed in water. e solution of the governing equations is determined with the help of the finite difference method and the L 1 algorithm. Furthermore, a detailed analysis of the effects of parameters of interest on velocity and temperature fields is offered. e following explicit and important inferences are obtained from the current numerical study: (i) Variation in the magnetic parameter together with the inertia coefficient F r and simultaneously changing generalized Reynold number Re have caused variation in the approximate solutions for both velocity and temperature equations, and these variations follow the same pattern as obtained from different numerical schemes, which endorses the accuracy of our proposed scheme. us, our proposed numerical scheme can be referred to as an alternate and efficient solver. (ii) Comparative analysis of the offered scheme suggests the accuracy of 10 − 04 for the solution u(x, y, t) when the mesh size is increased along x and y directions. In contrast, the accuracy for the solutions of T(x, y, t) is somewhat near 10 − 05 . (iii) e temperature profiles considerably increase against an increase in fractional derivative parameters (c), and accordingly, the viscoelastic effects of the fluids are enhanced. A remarkable rise in temperature is observed for Al 2 O 3 − ZrO 2 /BF as compared with ZrO 2 /BF. (iv) Furthermore, when the generalized Reynolds number (Re) is increased, the temperature profile of both nanofluids decreases and hence produces the optimum level of heat transfer over the permeable surface. (v) In investigating the effect of nanoparticles, hybrid nanoparticles have had a more efficient role in heat transfer enhancement than a single nanofluid over a nonuniform permeable surface. (vi) During the analysis of the impact of nanoparticles, it is revealed that hybrid nanoparticles play a more dominant role in enhancing heat transfer than single nanoparticles over a nonuniform permeable surface. It is noted that the thermal conductivity of nanofluid Al 2 O 3 − ZrO 2 /BF is increased by 6.24% (diff) when the volume concentration is increased. (vii) A visible increment in heat transfer rate is observed with an increase in Reynold number Re for hybrid and single nanofluids. e average skin friction coefficient follows the same trend when the value of α is increased.
Finally, it is established from the data analysis that hybrid nanofluids have significant implications in the applications of thermodynamics along with considerable scale control over the behavior of the synthetic fluid. Considering a precise application and particular engineering nanofluids for that application, the hybrid nanofluids can entirely change the field of heat transfer. us, we can say that in order to have a more efficient thermal system, it is suggested that hybrid nanofluids should be used instead of a single type of nanoparticle.
us, the challenge of manufacturing economic and stable products can be more efficiently addressed, leading to novel research avenues. Appendix f 1 (x, y, t) , 2tx n x − L x L y . (A.1)

Data Availability
e data used to support the findings of this study are included within the article.