The improved thermal efficiency of Prandtl–Eyring hybrid nanofluid via classical Keller box technique

Prandtl–Eyring hybrid nanofluid (P-EHNF) heat transfer and entropy generation were studied in this article. A slippery heated surface is used to test the flow and thermal transport properties of P-EHNF nanofluid. This investigation will also examine the effects of nano solid tubes morphologies, porosity materials, Cattaneo–Christov heat flow, and radiative flux. Predominant flow equations are written as partial differential equations (PDE). To find the solution, the PDEs were transformed into ordinary differential equations (ODEs), then the Keller box numerical approach was used to solve the ODEs. Single-walled carbon nanotubes (SWCNT) and multi-walled carbon nanotubes (MWCNT) using Engine Oil (EO) as a base fluid are studied in this work. The flow, temperature, drag force, Nusselt amount, and entropy measurement visually show significant findings for various variables. Notably, the comparison of P-EHNF's (MWCNT-SWCNT/EO) heat transfer rate with conventional nanofluid (SWCNT-EO) results in ever more significant upsurges. Spherical-shaped nano solid particles have the highest heat transport, whereas lamina-shaped nano solid particles exhibit the lowest heat transport. The model's entropy increases as the size of the nanoparticles get larger. A similar effect is seen when the radiative flow and the Prandtl–Eyring variable-II are improved.

Velocity component in x , y direction (m s −1 ) U w The velocity of the stretching sheet V ς Vertical velocity x, y Dimensional space coordinates (m) Greek symbols The fluid temperature of the surface The Single-walled carbon nanotubes MWCNT Multi-walled carbon nanotubes Liquid mechanics' limits are defined by the thin fluid or liquid layer in contact with the pipe's or an aircraft wing's surface. In the boundary layer, shear forces can damage the liquid. Given that the fluid is in touch with the surface, a range of speeds exists between the maximum and zero boundary layer speeds. Limits on the trailing edge of an aeroplane wing, for example, are smaller and thicker. A thickening of the flow occurs at the front or upstream end of these boundaries. In 1904, Prandtl proposed the concept of boundary layers to describe the flow behaviour of viscous fluid near a solid barrier (see Aziz et al. 1 ). Using the Navier Stoke equations, Prandtl constructed and inferred boundary layer equations for large Reynolds number flows. As a necessary simplification of the original Navier-Stokes equation, the boundary layer theory equations were critical. Studying wall jets, free jets, fluid jets, flow over a stretched platform/surface, and inductive flow from a shrinking plate helps develop the equations for these phenomena. Boundary layer equations are often solved using a variety of boundary conditions that are specific to a given physical model. For a magnetohydrodynamics (MHD) fluid flow with gyrotactic microorganisms, Sankad et al. 2 found that the magnetic and Peclet numbers may be utilised to reduce the thermal boundary layer thickness. After that, Hussain et al. 3 discovered that the thickness of the thermal boundary layer increases as a Casson liquid flows towards the growing porous wedge due to convective heat transfer. The literature has several experiments with various physical parameter impacts on boundary layer flow [4][5][6] and multiple liquids 7,8 .
A hybrid nanofluid is now attracting the attention of many researchers. Hybrid nanofluids are cutting-edge nanofluids that combine two different types of nanoparticles in a single fluid. The thermal properties of the hybrid nanofluid are better than those of the primary liquid and nanofluids. In machining and manufacturing, hybrid nanofluids are commonly utilised in solar collectors, refrigeration, and coolants. According to Suresh et al. 9 , copper nanoparticles in the alumina matrix mixed at most modest and sufficient levels may preserve the hybrid nanofluid's strength, first introduced in 9 . Despite having a lower thermal conductivity than copper nanoparticles, alumina nanoparticles have excellent chemical inactivity and stability. Yildiz et al. 10 developed an equivalence between theoretical and experimental thermal conductivity models for heat transfer performance in hybrid-nanofluid. In comparison to a mono nanofluid, the hybridisation of nanoparticles improved heat transfer at a lower particle percentage (Al 2 O 3 ). Waini et al. 11 investigated a hybrid nanofluid's unsteady flow and heat transfer using a curved surface. As the surface curvature changed, the presence of dual solutions resulted in intensification in the volume percentage of copper nanoparticles. Many years later, Qureshi et al. 12  www.nature.com/scientificreports/ increasing the barrier's radius improves heat transfer by as much as 119%. In addition, the horizontal orientation of the cylinder only supports a heat transfer efficiency of 2.54%. Mabood and Akinshilo 13 investigated the influence of uniform magnetic and radiation on the heat transfer flow of Cu-Al 2 O 3 /H2O hybrid nanofluid flowing over the stretched surface. Discoveries such as those made at the science fair show how radiation speeds up heat transport while magnetic forces slow it down. Further hybrid nanofluid studies and experiments have been conducted by these researchers [14][15][16][17][18][19][20][21][22][23][24] . The Cattaneo-Christov heat flux model describes the heat transfer in viscoelastic flows caused by an exponentially expanding sheet. There may be a relationship between thermal relaxation time and the boundaries of this model. Dogonchi and Ganji 25 researched unstable squeezing MHD nanofluid flow across parallel plates using a Cattaneo-Cristov heat flux model some years ago. The thermal relaxation parameter, they found, slowed heat transfer. Additionally, Muhammad et al. 26 discovered that when thermal relaxation increased, the fluid temperature decreased. Other researchers have used the Cattaneo-Christov heat flux model to examine fluid flow and determine the physical features that thermal relaxation affects. Scholars like 27-32 may be found in the literature as examples of this group. Even the temperature of a nanofluid may be reduced by the thermal relaxation parameter, according to Ali et al. 33 . This finding is critical to the contemporary food, medicinal, and aerospace industries. Waqas et al. 34 introduced mathematical modelling using the Cattaneo-Christov model for hybrid nanofluid flow in a rocket engine. The finding exposed that the temperature is reduced when thermal relaxation and melting parameters vary, but the Biot number increases. Other types of hybrids nanofluid characteristics using the Cattaneo-Christov model have been discussed by Haneef et al. 35 . The vital discovery uncovered an escalation causes shrinkage in wall shear stress in momentum relaxation time. Different encounters were found by Reddy et al. 36 in the Cattaneo-Christov model problem for hybrid dusty nanofluid flow. It reveals that dusty hybrid nanofluid has a better heat transfer method than hybrid nanofluid.
Nevertheless, a few years back, a new type of fluid was found called hybrid nanofluid, and many researchers have been eager to search for the characteristics of this type of fluid since then. The research for finding the aspect of non-Newtonian hybrid nanofluid also needed to be done. Latterly, Yan et al. 37 have conducted an investigation towards the rheological behaviour of non-Newtonian hybrid nanofluid for a powered pump. They reported at the highest volume fraction hybrid nanofluid, the viscosity reduced at most 21%. Nabwey and Mahdy 38 are doing an inclusive exploration of micropolar dusty hybrid nanofluid. The finding indicates that the temperature fluctuation in both the micropolar hybrid nanofluid and dust phases is strengthened by increased thermal relaxation. Several investigations have been carried out for the different types of non-Newtonian hybrid nanofluid, such as aluminium alloy nanoparticles by Madhukesh et al. 39 , MWCNT-Al 2 O 3 /5W50 by Esfe et al. 40 and ZnO-Ag/ H 2 O by He et al. 41 in the literature. Despite that, only a few research available in the literature investigating the viscoelastic hybrid nanofluid behaviour. Several models can be used to examine the physical properties of the viscoelastic fluid, including the power-law model, the Prandtl fluid model, and the Prandtl-Eyring model. The power-law model predicts the non-linear relationship between deformation rate and shear stress. It has been hypothesised that shear stress is connected to the sine inverse function of deformation rate by the Prandtl model and that it is related to the hyperbolic sine function of deformation rate by the Prandtl-Eyring model. Hussain et al. 42 have investigated the physical aspect of MHD Prandtl-Eyring fluid flow and reported that at all positions in the flow domain, a substantial rise in momentum transportation had been seen against an increase in the fluid parameter. Rehman et al. 43 added in the findings that Prandtl-Eyring liquid particles are subjected to drag forces in a flow when their skin friction coefficients are high (or low). A similar discovery has been conveyed by Khan et al. 44 which the skin friction improves for the Prandtl-Eyring nanofluid. Later, Akram et al. 45 model a MHD Prandtl-Eyring nanofluid peristaltic pumping in an inclined channel. This study demonstrates that the wall tension and mass parameters have a rising influence on axial velocity, whereas the wall damping parameter has a decreasing impact. Li et al. 46 have explored the entropy of the Prandtl-Eyring fluid flow model over a rotating cone. The result demonstration the velocity and temperature have been shown to behave differently when the viscosity parameter increases in magnitude. Latest study for the Prandtl-Eyring hybrid nanofluid model being carried out by Jamshed et al. 47 . The outcome was mentioning the entropy upsurged with radiative flux and Prandtl-Eyring parameter.
The famous numerical technique for solving non-linear boundary layer equations in fluid mechanics is derived by Keller and Cebeci 48 called Keller Box Method (KBM). It is being popularised by Cebeci and Bradshaw 49 . The technique is known for highly accurate and time computation in solving non-linear problems. A lot of investigations of fluid dynamics have been solved using KBM in the literature. Bilal et al. 50 implemented the KBM for solving Williamson fluid flow towards a cylindrical surface and found the results are comparable with other published results. Similar numerical computation was reported by Swalmeh et al. 51 in solving the micropolar nanofluid over a solid sphere using KBM. The computed solution being reported as having a good agreement with the solution computed by bvp4c (MATLAB). The KBM is a universal solver since it is proven can solve another type of mathematical modelling, for instance, Carreau fluid model (Salahuddin 52 ), micropolar fluid (Singh et al. 53 ), viscous fluid model (Bhat and Katagi 54 ), Prandtl nanofluid (Habib et al. 55 ), MHD nanofluid (Zeeshan et al. 56 ) and third-grade nanofluid (Abbasi et al. 57 ).
Size and distribution descriptors should be chosen to offer the most significant discrimination for particulate quality concerning specific attributes or characterisation of a manufacturing process, depending on their use. If particle form affects these attributes, the shape and distribution of the particles should be studied in addition to their size. Qualitative terminology like fibres or flakes can be used, or quantitative terms like elongation, roundness, and angularity can also be used. Other quantitative terms include percentages of certain model forms and fractal dimensions. Despite the importance of the particle shape, only a few research can be found in the literature, such as [58][59][60] . The latest research has been done by Sahoo 61 , which claimed that the particle shapes heavily influence the thermo-hydraulic performance of a ternary hybrid nanofluid. Similar findings have been illustrated by Elnaqeeb et al. 62 in hybrid nanofluid flow with the impact of suction and stretching surface. www.nature.com/scientificreports/ Meanwhile, Rashid et al. 63 suggested that the temperature and Nusselt number profiles demonstrate the sphere shape nanoparticles have superior temperature disturbance and heat transmission on hybrid nanofluid flow with the influence of relevant factors. A few publications have examined the impact of the porosity material, viscid dissipative flow, Cattaneo-Christov heat flow and thermal radiative flow shape-factor along the elongated surface using nanofluid Tiwari-Das type on P-EHNF entropy generation. However, none of these papers has addressed these issues. In the Tiwari-Das (monotonic model), the fluid, speed, and temperature are all the same. As a result, the model is simpler and easier to solve when using the single-phase technique numerically. However, this technique has the drawback of resulting in numerical effects that differ from experimental results in some cases. Nanoparticle concentrations in this model volume range from 3 to 20%. Numerical results could only mimic the effects of SWCNT-EO, MWCNT-EO hybrids, and conventional nanofluids in this study. Thus, in order to bridge the gap, the current research focuses on the solid-fluid characteristics impacts and the level of chaos in the boundary layer using the Keller-box technique of P-EHNF.

Flow model formulations
The mathematical flow equations shows the moved horizontal plate with the irregular expanding velocity 64 : where Y = ∞ signify the temperature variation amount, heat of surface, and surrounds correspondingly. The plate is supposed to be slippery, and the surface is subjected to a temperature variation.
Primary addition SWCNT nano solid-particles synthesise the hybrid nanofluid in the EO-based liquid at an interaction volume fraction ( φ ST ) and it is fixed at 0.09 during the examination. MWCNT nano molecules have been extended in combination to obtain a hybrid nanofluid at the concentrated size ( φ MT ).  www.nature.com/scientificreports/ the appropriate connection conditions are as follows, which can be located in Aziz et al. 67 :

Prandtl-Eyring
We formulate the Y = as a fluid heat. Other vital parameters are surface permeability V ς , heat transfer coefficient h ς , porosity (k) and heat conductivity of firm k ς . Physical features identical, Convectional animated surface experienced its heat loss through conductive (Newtonian thermal) and flowing swiftness close to the sheet is comparative to the cut stress exerts in it (slippy form) are deliberate.
Heat-physical possessions of P-ENF. Nano solid particles dispersed in EO induce improved thermophysical characteristics. The next Table 1 equations summarize P-ENF substance variables 68,69 .
φ is the nano solid-particle size coefficient. µ f , ρ f , (C p ) f and κ f are dynamical viscidness, intensity, functioning thermal capacity, and thermal conductivity of the standard fluid, respectively. The additional characteristics ρ s , (C p ) s and κ s are the concentration, effective heat capacitance, and heat conductance of the nano molecules, correspondingly.
Thermo-physical properties of P-EHNF. The primary assumption of hybrid nanofluids is the suspension of two distinct forms of nano solid particles inside the basis fluid 70 . This assumption improves the capacity for heat transmission of common liquids and is a higher heat interpreter than nanofluids. P-EHNF variables content is summarised in Table 2 71,72 .
In Table 2, µ hnf , ρ hnf , ρ(C p ) hnf and κ hnf are mixture nanofluid functional viscidness, concentration, exact thermal capacitance, and thermal conductance. φ is the volume of solid nano molecules coefficient for mono (4) www.nature.com/scientificreports/ nanofluid and φ hnf = φ ST + φ MT is the nano solid particles magnitude measurement for the combination nanofluid. µ f , ρ f , (C p ) f , κ f and σ f are functional viscidness, density, exact thermal capacity, and heat conductivity of the base fluid. ρ p 1 , ρ p 2 , (C p ) p 1 , (C p ) p 2 , κ p 1 and κ p 2 are the density, specific heat capacity, and thermal conductivity of the nano-molecules.
Nano solid-particle shape-factor m. The scale of the multiple nano solid-particles is defined as the shaped-nanoparticles factor. Table 3 shows the importance of the experiential form factor for different particle forms (for instance, see Xu and Chen 73 ).

Nano solid-particles and basefluid lineaments.
In this analysis, the material characteristics of the primary oil-based liquid of the engine are specified in Table 4 74,75 .

Rosseland approximation.
Radiative flow only passes a shortened distance because its non-Newtonian P-EHNF is thicker. Because of this, the approximation for radiative fluxing from Rosseland 76 is utilised in formula (4).
herein, σ * signifies the constant worth of Stefan-Boltzmann and k * symbols the rate.

Dimensionless formulations model
Given the similarity technology that transforms the governing PDEs into ODEs, the BVP formulas (2)-(6) are modified. Familiarising stream function ψ in the formula 75 The specified similarity quantities are into Eqs. (2)-(4). We get with (7)

Prandtl-Eyring parameter-I
Drag-force and Nusselt number. The drag-force C f combined with the Nusselt amount (Nu x ) are the interesting physical amounts that controlled the flowing and specified as 66 where τ w and q w determine as The dimensionless transmutations (9) are implemented to obtain where Nu x means Nusselt aggregate and C f states drag force constant. Re x = u w x ν f is local Re built in the extended swiftness u w (x).

Classical Keller box technique
Because of its rapid convergence, the Keller-box approach (KBM) 77 is used to find solutions for model formulas (Fig. 2). KBM is used to find the localised solve of (10) and (11) with constraints (12). The policy of KBM is specified as next: Stage 1: ODEs adaptation. In the early stage, all of the ODEs must be changed into 1st-order ODEs (10)-(12) where j is used for the spacing in h in a horizontal direction to show the position of the coordinates. The solution to the problem is to be found without any initial approximation. It is very crucial for finding velocity, temperatures, temperature variations, and entropy to make a preliminary assumption between = 0 and = ∞ . The frameworks from the result have been approximated solutions provided as they can happen the boundary conditions of the problem. It is imperative to remark that the results must be equalled with different preliminary estimations are chosen, but the replication computation and time are varied which have been taken for conducting the calculations (see Fig. 3): By implementing significant differences, difference equivalences are figured, and functions are used to replace the mean values. The 1st-order ODEs (19)- (23) have been modified to algebraic formulas which are non-linear. By replacing this in (25) to (29) and after overlooking the higher-elevated bounds of θ i j a linear tri-diagonal equation scheme has been resulting as follows: where 5 × 5 block-sized matrix is denoted by F that corresponds to the size of J × J . However, the vector of order J × 1 is represented by θ and p . An esteemed LU factorising method is used for solving θ later. The equation Fθ = p denotes that F with an array θ is used to yield a production array marked by p . Further, F is splinted into lower and upper trigonal matrices, i.e., F = LU can be written as LUθ = p . Let Uθ = y tends to Ly = p, which is used to provide the solution of y . Further, the values of y computed are replaced into the equation Uθ = y for solving θ . The technique of back-substitution has been implemented as this is the easy method to find a solution.

Code verification
On the other, by measuring the heat transmission rate outcomes from the current technique against the recent results available in the literature 78,79 , the validity of the method was evaluated. Table 5 summarises the comparing of reliabilities current during the researches. Nevertheless, the outcomes of the current examination are exceedingly accurate.

Second law of thermodynamics
Porous media generally increase the entropy of the system. Jamshed et al. 80 and Jamshed 81 described the nanofluid entropy production by: The non-dimensional formulation of entropy analysis is as follows [82][83][84] , www.nature.com/scientificreports/ By formula (9), the non-dimensional entropy formula is: Here R e is the Reynolds number, B ς signifies Brinkmann amount and symbols the non-dimensional variation of the temperature.

Results and discussion
An adequate discussion is indicated by numerical results that reach the model described before. As a result of these potential parameters, the values for A * 1 , A * 2 , K ς , φ , � ς , S , N ς ,ε ς , E ς H ς , R e and B ς been illustrated. These parameters show the physical performance of the non-dimensional quantities in Figs. 4, 5, 6, 7, 8, 9, 10, 11, 12, Table 5.   www.nature.com/scientificreports/ 13,14,15,16,17,18,19,20,21,22,23,24,25 and 26, such as velocity, energy, and entropy production. The results are obtained for Cu-EO normal P-ENF and MWCNT-SWCNT/EO non-Newtonian P-EHNF. The coefficient of skin friction and temperature variations are shown in Table 6. For example, the default values were 1.0 for A * 1 and 0.4 for A * 2 , K ς was set to be equal to 0.1, and φ = 0.18, φ MT was set to 0.09, � ς was set to 0.3, S was set to 0.4, and N ς was set to 0.3, ε ς was set to 0.1, E ς was set to 0.3, H ς was set to 0.3, and R e and B ς was set to 5.    Fig. 4. As the value of A * 1 was elevated, so was the velocity profile for both fluids. The physical reason for this occurrence is that it causes the fluid's viscosity to decrease, reducing resistance while boosting fluid velocity. MWCNT-SWCNT nanofluid, on the other hand, has faster   www.nature.com/scientificreports/ acceleration than SWCNT nanofluid. It can be explained as the hybrid nanofluid have an enormous density impact rather than the nanofluid. The temperature curve for the Prandtl-Eyring parameter A * 1 is shown in Fig. 5. MWCNT-SWCNT hybrid nanofluid had a lower temperature profile since the value of A * 1 was raised, while the Cu nanofluid had a higher temperature profile. More heat can be conveyed faster when the caused in this lowered manner due to velocity improve and expand. Another important distinction is that the hybrid nanofluid   www.nature.com/scientificreports/ exhibits significantly reduced thermal conductivity when compared to pure nanofluid. Figure 6 depicted the Prandtl-Eyring hybrid nanofluid entropy fluctuation based on its parameter A * 1 . The quantity of entropy produced decreased as the amount of A * 1 enhanced. MWCNT-SWCNT fluid exhibited a lower entropy value than SWCNT hybrid nanofluid, even though their values were the same at one point in the graph. This phenomenon   www.nature.com/scientificreports/ occurs due to the low temperatures reducing hybrid nanofluid mobility, causing the system's entropy to proliferation.
Influence of Prandtl-Eyring parameter A * 2 . There was an influence of Prandtl-Eyring Parameter A * 2 on the Prandtl-Eyring hybrid nanofluid temperature, velocity, and entropy production profile (see Figs. 7, 8). Fig-Figure 18. Entropy change with � ς .  www.nature.com/scientificreports/ ure 7 depicts the varying A * 2 with velocity . The velocity profile narrows as A * 2 rises, with MWCNT-SWCNT/EO achieving a higher top speed than SWCNT-EO. Hybrid nanofluid particles have resistance due to the fact that they vary inversely with momentum diffusivity. As a result, the flow's velocity will be reduced with A * 2 . This phenomenon is because SWCNT-EO has a higher density and hence has a thicker flow than MWCNT-SWCNT/EO, making the fluid challenging to transport. Figure 8 shows the temperature change after A * 2 has had its impact. As   www.nature.com/scientificreports/ the value of A * 2 grew, so did the temperature, with SWCNT-EO quickly reaching the desired temperature. The occurrence happens because the flow velocity dropped, and as a result, the heat transmission from the surface was degraded. Figure 9 shows the change in entropy according to the Prandtl-Eyring parameter A * 2 . The entropy profile grew as the value of A * 2 grew, showing a clear connection between the two. It suggested that A * 2 amplifying the impediment in the system, resulting in the entropy of the developing system being elevated.    Figure 11 displays the expansion of the porous medium variable ( K ς ) results in better heat dispersion throughout the domain. When a hole is made in a porous medium, the flow slows down, allowing more time to collect heat from the surface. This phenomenon improves the thermal distribution around the area. Since particle motions across porous media are sluggish, the porosity aids in the irreversibility of energy transfer across the domain during entropy production (N G ) (Fig. 12).
Effect of nanomolecules size φ and φ hnf . The efficacy of the nanofluid and hybrid versions appears to be determined by the fractional nanoparticle size in the base fluid. The more excellent fractional range of nanoparticles reduces flowability because of the additional load it adds. For some reason, the fractional upgrade prefers the hybrid nanofluid over the single-nanofluid, which flows lower in Fig. 13. This incidence displayed the primary reason for utilising nano-and hybrid-based fluid mixtures because of their exceptional heat transmission properties. This degradation occurs as a result of excessive nanoparticle surface area and higher hybrid nanofluid density. As the fractional volume of both kinds of flow fluids improved, so did the resultant thermal distribution, as shown in Fig. 14. Because of the temperature difference, when the nano molecule size is reduced, the molecules will disperse in the far-field flow. The thermal boundary layer's thickness will rise as a result of this change. The minimal size of nano molecules can be utilised to create the lowest possible temperature profile, as determined through experimentation. Figure 15 exhibits the leading nanofluid varies in the middle and settles down to the hybrid nanofluid at the far end, with energy entropy fluctuations also intensifying for fractional volume. SWCNT-EO has a greater entropy than MWCNT-SWCNT/EO because the hybrid nanofluid has a far higher thermal conductivity than nanofluid. Table 6.
Values of C f Re 1/2 x and Nu x Re  Figure 16 illustrates the flow conditions in Prandtl-Eyring fluid mixtures are primarily centred upon the viscous behaviour. Due to this occurrence, slip conditions become incredibly critical in fluids as a whole. For a hybrid suspended Prandtl-Eyring nanofluid, the viscous nature and higher levels of flow slip generate more complex fluidity circumstances, with the result that the fluidity of the single nanofluid drops even more rapidly. Due to the flow hierarchy, the SWCNT-EO nanofluid maintains a higher temperature state than the MWCNT-SWCNT/EO hybrid nanofluid, which is depicted in Fig. 17. The improvement in boundary layer viscosity due to the decline in velocity will have a similar effect. As a result, it will have skyrocketed the flow's temperature. Because the hybrid nanofluid has less viscosity than the conventional nanofluid, it is predicted MWCNT-SWCNT/EO to have a lower temperature than SWCNT-EO. A descending trend in entropy formation can be seen for higher slip parameters because the slipped flow acts against the domain's entropy formation.
Thermal radiative variable N ς and relaxation time parameter ( ε ς ) influence. Figures 19 and 20 highlights the actual status of thermal diffusion and entropy generation under enhanced heat radiative flow limitation N ς . Thermally diffusing nanofluids have a propensity to rise in temperature past the interesting domain, boosting the heat transmission burden for radiation constrictions on the transient nanofluid. This temperature rise may be explained in a physical sense by supposing that thermal radiation is converted into electromagnetic energy. As a result, the distance from the surface from which radiation is emitted rises, ultimately superheating the boundary layer flow. As a result, the thermal radiative variable is critical in determining the system's temperature profile. A limit on radiative flow N ς via entropy generation is illustrated in Fig. 20 by the overfilled dispersions. For different N ς values, the entropic side-by-side leans toward developing more in MWCNT-SWCNT/ EO than in SWCNT-EO nanofluid. A reasonable explanation for this occurrence is the system's irreparable heat transfer mechanism is entirely irreversible. According to Fig. 21, greater values of the relaxation time parameter cause a rise in the temperature of the Tangent hyperbolic hybrid nanofluids, as seen in the graph. As the temperature drops, the thickness of the thermal boundary layer reduces. Table 5 shows that when the rate of heat output efficiency, the effectiveness of the thermal system improves as well. Figure 22 shows the impact of engine oil-based nanofluid entropy profiles. The velocity profile, on the other hand, shows no change, while the entropy of the system increases with varying values ε ς .
Effect of the diverse solid particle shape m. It is well-known that NPs have high thermal conductivity and transfer rates under a variety of physical conditions. In porous medium difficulties, such nano-level particles become an issue, modelled using the shape variable (m) in this study. From spherical (m = 3) to lamina (m = 16.176), the forms considered here ranged. To improve the thermal state, Fig. 23 indicates that nanoparticle shapes impact it. In comparison to SWCNT-EO mono nanofluid, the MWCNT-SWCNT/EO hybrid nanofluid has a more significant form impact. Hybrid nanofluid has a broader thermal layer boundary and a more excellent thermal distribution than nanofluid. Even in the MWCNT-SWCNT/EO hybrid nanofluid, the lamina (m = 16.176) shaped particles remain ahead of the others. The main physical reason for this phenomenon is the lamina shape particles has the most remarkable viscosity while the sphere has the minimum viscosity. It is also noted that at a higher temperature, the viscosity of the particles will be diminished. This phenomenon happens because of the temperature-dependent shear-thinning characteristic. The profiles in Fig. 24 indicate the form factors have a more substantial influence in MWCNT-SWCNT/EO NHF, which has a higher entropy rate than SWCNT-EO mono nanofluid, even though the morphologies of the particles have a much less impact.

Entropy variations for Reynolds number ( R e ) and Brinkman number ( B ς ). Figure 25 depicted as
Reynolds number proliferations (R e ) , the entropy rate (N G ) improves as well. An aggregate Reynolds number supports nanoparticle mobility in porous media because of the dominance of inertial over viscous forces in the system. Consequently, entropy can be generated over the domain. MWCNT-SWCNT/EO HNF generated a higher entropy rate than MWCNT-SWCNT/EO nanofluid because of the combined efficiency of the particles. The Brinkman number B ς was used to describe the heat created by viscous properties because it enhances the generated heat above and beyond other thermal inputs. The heightened heat-inducing ability of such viscosity enhancement promotes entropy production in the system as a whole (N G ) . Figure 26 illustrated the elevated entropy layers, which improved the Brinkman number B ς values. The primary feature of viscous dislocation heat produces a decrease in the escalating Brinkman numbers, which theoretically leads to a higher rate of entropy development.

Final results and future guidance
The entropy production and heat transmission by a Prandtl-Eyring hybrid nanofluid (P-EHNF) over a stretched sheet is studied. By utilising a single model phase, a computational model may be developed. Several physical characteristics are used to derive the results. These include changes in velocity, energy, and entropy. Cattaneo-Christov heat flux ε ς is also discussed in this study, as are the effects of Prandtl-Eyring parameters A * 1 and A * 2 as well as K ς for porous medium, nanomolecular size φ and φ hnf , � ς for velocity slip, thermal radiative variable N ς and Biot number H ς as well as various solid particle shapes m , R e and B ς . The following are the study's significant findings: 1. In comparison to traditional Prandtl-Eyring nanofluid (SWCNT/-EO), hybrid Prandtl-Eyring nanofluid (MWCNT-SWCNT/EO) is shown to be a superior thermal conductor.